Papers for Tuesday, Jun 07 2022

The stellar initial mass function (IMF) is predicted to depend upon the temperature of gas in star-forming molecular clouds. The introduction of an additional parameter, $T_{IMF}$ , into photometric template fitting, suggest most galaxies obey an IMF top-heavier than the Galactic IMF. The implications of these revised fits on mass functions, quiescence and turnoff are discussed. At all redshifts the highest mass galaxies become quiescent first with the turnoff mass decreasing towards the present. The synchronous turnoff mass across galaxies suggests quiescence is driven by universal mechanisms rather than by stochastic or environmental processes.

(Abridged) The question of how most stars in the Universe form remains open. While star formation predominantly occurs in young massive clusters, the current framework focuses on isolated star formation. One way to access the bulk of protostellar activity within star-forming clusters is to trace signposts of active star formation with emission from molecular outflows. These outflows are bright in water emission, providing a direct observational link between nearby and distant galaxies. We propose to utilize the knowledge of local star formation as seen with molecular tracers to explore the nature of star formation in the Universe. We present a large-scale statistical galactic model of emission from galactic active star-forming regions. Our model is built on observations of well-resolved nearby clusters. By simulating emission from molecular outflows, which is known to scale with mass, we create a proxy that can be used to predict the emission from clustered star formation at galactic scales. We evaluated the impact of the most important global-star formation parameters (i.e., initial stellar mass function (IMF), molecular cloud mass distribution, star formation efficiency (SFE), and free-fall time efficiency) on simulation results. We observe that for emission from the para-H2O 202 - 111 line, the IMF and molecular cloud mass distribution have a negligible impact on the emission, both locally and globally, whereas the opposite holds for the SFE and free-fall time efficiency. Moreover, this water transition proves to be a low-contrast tracer of star formation. The fine-tuning of the model and adaptation to morphologies of distant galaxies should result in realistic predictions of observed molecular emission and make the galaxy-in-a-box model a tool to analyze and better understand star formation throughout cosmological times.

Well-characterized binary systems will provide valuable opportunities to study the conditions that are necessary for the onset of both auroral and non-auroral magnetospheric radio emission in the ultracool dwarf regime. We present new detections of non-auroral "quiescent" radio emission at 4-8 GHz of the three ultracool dwarf binary systems GJ 564 BC, LP 415-20, and 2MASS J21402931+1625183. We also tentatively detect a highly circularly polarized pulse at 4-6 GHz that may indicate aurorae from GJ 564 BC. Finally, we show that the brightest binary ultracool dwarf systems may be more luminous than predictions from single-object systems.

We report on a new capability added to our general relativistic radiation-magnetohydrodynamics code, Cosmos++: an implicit Monte Carlo (IMC) treatment for radiation transport. The method is based on a Fleck-type implicit discretization of the radiation-hydrodynamics equations, but generalized for both Newtonian and relativistic regimes. A multiple reference frame approach is used to geodesically transport photon packets (and solve the hydrodynamics equations) in the coordinate frame, while radiation-matter interactions are handled either in the fluid or electron frames then communicated via Lorentz boosts and orthonormal tetrad bases attached to the fluid. We describe a method for constructing estimators of radiation moments using path-weighting that generalizes to arbitrary coordinate systems in flat or curved spacetime. Absorption, emission, scattering, and relativistic Comptonization are among the matter interactions considered in this report. We discuss our formulations and numerical methods, and validate our models against a suite of radiation and coupled radiation-hydrodynamics test problems in both flat and curved spacetimes.

We present a comprehensive optical and near-infrared census of the fields of 90 short gamma-ray bursts (GRBs) discovered in 2005-2021, constituting all short GRBs for which host galaxy associations are feasible ($\approx$ 60% of the total Swift short GRB population). We contribute 245 new multi-band imaging observations across 49 distinct GRBs and 25 spectra of their host galaxies. Supplemented by literature and archival survey data, the catalog contains 335 photometric and 40 spectroscopic data sets. The photometric catalog reaches $3\sigma$ depths of $\gtrsim 24-27$ mag and $\gtrsim 23-26$ mag for the optical and near-infrared bands, respectively. We identify host galaxies for 84 bursts, in which the most robust associations make up 54% (49/90) of events, while only a small fraction, 6.7%, have inconclusive host associations. Based on new spectroscopy, we determine 17 host spectroscopic redshifts with a range of $z\approx 0.15-1.6$ and find that $\approx$ 25-44% of Swift short GRBs originate from $z>1$. We also present the galactocentric offset catalog for 83 short GRBs. Taking into account the large range of individual measurement uncertainties, we find a median of projected offset of $\approx 7.9$ kpc, for which the bursts with the most robust associations have a smaller median of $\approx 4.9$ kpc. Our catalog captures more high-redshift and low-luminosity hosts, and more highly-offset bursts than previously found, thereby diversifying the population of known short GRB hosts and properties. In terms of locations and host luminosities, the populations of short GRBs with and without detectable extended emission are statistically indistinguishable. This suggests that they arise from the same progenitors, or from multiple progenitors which form and evolve in similar environments. All of the data products are available on the BRIGHT website.

We present the stellar population properties of 68 short gamma-ray burst (GRB) host galaxies, representing the largest uniformly-modeled sample to date. Using the Prospector stellar population inference code, we jointly fit the photometry and/or spectroscopy of each host galaxy. We find a population median redshift of $z=0.60^{+0.90}_{-0.34}$ (68% confidence), including 10 new or revised photometric redshifts at $z\gtrsim1$. We further find a median mass-weighted age of $t_m=0.85^{+2.49}_{-0.57}$ Gyr, a stellar mass of $\log(M_*/M_\odot)=9.69^{+0.95}_{-0.66}$, a star formation rate of SFR=$1.46^{+10.55}_{-1.38}M_\odot$/yr, a stellar metallicity of $\log(Z_*/Z_\odot)=-0.4^{+0.46}_{-0.42}$, and a dust attenuation of $A_V=0.51^{+0.94}_{-0.43}$ mag (68% confidence). Overall, the majority of short GRB hosts are star-forming ($\approx85$%), with small fractions that are either transitioning ($\approx6$%) or quiescent ($\approx$9%); however, we observe a much larger fraction ($\approx60$%) of quiescent and transitioning hosts at $z\lesssim0.25$, commensurate with galaxy evolution. We find that short GRB hosts populate the star-forming main sequence of normal field galaxies, but do not include as many high-mass galaxies, implying that their binary neutron star (BNS) merger progenitors are dependent on a combination of host star formation and stellar mass. The distribution of ages and redshifts implies a broad delay-time distribution, with a fast-merging channel at $z>1$ and a decreased BNS formation efficiency at lower redshifts. Moreover, we find that short GRBs originate in environments with a wide range of stellar metallicities. If short GRB hosts are representative of BNS merger hosts within the horizon of current gravitational wave detectors, these results can inform future searches for electromagnetic counterparts. All of the data and modeling products are available on the BRIGHT website.

Extremely large telescopes (ELTs) provide an opportunity to observe surface inhomogeneities for ultracool objects including M dwarfs, brown dwarfs (BDs), and gas giant planets via Doppler imaging and spectro-photometry techniques. These inhomogeneities can be caused by star spots, clouds, and vortices. Star spots and associated stellar flares play a significant role in habitability, either stifling life or catalyzing abiogenesis depending on the emission frequency, magnitude, and orientation. Clouds and vortices may be the source of spectral and photometric variability observed at the L/T transition of BDs and are expected in gas giant exoplanets. We develop a versatile analytical framework to model and infer surface inhomogeneities which can be applied to both spectroscopic and photometric data. This model is validated against a slew of numerical simulations. Using archival spectroscopic and photometric data, we infer star spot parameters (location, size, and contrast) and generate global surface maps for Luhman 16B (an early T dwarf and one of our solar system's nearest neighbors at a distance of approximately 2 pc). We confirm previous findings that Luhman 16B's atmosphere is inhomogeneous with time-varying features. In addition, we provide tentative evidence of longer timescale atmospheric structures such as dark equatorial and bright mid-latitude to polar spots. These findings are discussed in the context of atmospheric circulation and dynamics for ultracool dwarfs. Our analytical model will be valuable in assessing the feasibility of using ELTs to study surface inhomogeneities of gas giant exoplanets and other ultracool objects.

We present a unique implementation of Python coding in an asynchronous object-oriented programming (OOP) framework to fully automate the process of collecting data with the George Mason University (GMU) Observatory's 0.8-meter telescope. The goal of this project is to perform automated follow-up observations for the Transiting Exoplanet Survey Satellite (TESS) mission, while still allowing for human control, monitoring, and adjustments. Prior to our implementation, the facility was computer-controlled by a human observer through a combination of webcams, TheSkyX, ASCOM Dome, MaxIm DL, and a weather station. We have automated slews and dome movements, CCD exposures, saving FITS images and metadata, initial focusing, guiding on the target, using the ambient temperature to adjust the focus as the telescope cools through the rest of the night, taking calibration images (darks and flats), and monitoring local weather data. The automated weather monitor periodically checks various weather data from multiple sources to automate the decision to close the observatory during adverse conditions. We have organized the OOP code structure in such a way that each hardware device or important higher-level process is categorized as its own object class or "module" with associated attributes and methods, with inherited common methods across modules for code reusability. To allow actions to be performed simultaneously across different modules, we implemented a multithreaded approach where each module is given its own CPU thread on which to operate concurrently with all other threads. After the initial few modules (camera, telescope, dome, data I/O) were developed, further development of the code was carried out in tandem with testing on sky on clear nights. The code, in its current state, has been tested and used for observations on 171 nights, with more planned usage and feature additions.

The NASA InSight mission has helped to measure the deep interior of Mars using observations of seismic waves excited by marsquakes. Currently, installation of seismometers on the moon is foreseen. We review the case for seismic experiments on all major planetary bodies of the solar system. We discuss scientific goals in accordance with the Decadal survey for planetary science and astrobiology and the ESA Voyage 2050 program as well as technical challenges and potential mission concepts, to answer the question: Where could we do seismology on other planets and why should we do it?

Magnetic fields in the turbulent interstellar medium (ISM) are a key element in understanding Galactic dynamics, but there are many observational challenges. One useful probe for studying the magnetic field component parallel to the line of sight (LoS) is Faraday rotation of linearly polarized radio synchrotron emission, combined with H$\alpha$ observations. HII regions are the perfect laboratories to probe such magnetic fields as they are localized in space, and are well-defined sources often with known distances and measurable electron densities. We chose the HII region Sharpless 2-27 (Sh 2-27). By using a map of the magnetic field strength along the LoS ($B_{\parallel}$) for the first time, we investigate the basic statistical properties of the turbulent magnetic field inside Sh 2-27. We study the scaling of the magnetic field fluctuations, compare it to the Kolmogorov scaling, and attempt to find an outer scale of the turbulent magnetic field fluctuations. We estimate the median value of $n_e$ as $7.3\pm0.1$ cm$^{-3}$, and the median value of $B_{\parallel}$ as $-4.5\pm0.1$ $\mu$G, which is comparable to the magnetic field strength in diffuse ISM. The slope of the structure function of the estimated $B_{\parallel}$-map is found to be slightly steeper than Kolmogorov, consistent with our Gaussian-random-field $B_{\parallel}$ simulations revealing that an input Kolmogorov slope in the magnetic field results in a somewhat steeper slope in $B_{\parallel}$. These results suggest that the lower limit to the outer scale of turbulence is 10 pc in the HII region, which is comparable to the size of the computation domain. This may indicate that the turbulence probed here could actually be cascading from the larger scales in the ambient medium, associated with the interstellar turbulence in the general ISM, which is illuminated by the presence of Sh 2-27.

The dynamics of the inner Galaxy contain crucial clues for untangling the evolutionary history of the Milky Way. However, the inner Galaxy's gravitational potential is poorly constrained, partly because the length of the Galactic bar is currently under debate with length estimates ranging from 3.5-5 kpc. We present a novel method for constraining the length and pattern speed of the Galactic bar using 6D phase space information to directly integrate orbits. We verify our method with N-body simulations and find that the maximal extent of orbits in the bar is not always consistent with that of the potential used to calculate the orbits. It is only consistent when the length of the bar in said potential is similar to the N-body model from which the initial positions and velocities of the stars are sampled. When we apply the orbit integration method to $\approx$210,000 stars in APOGEE DR17 and $Gaia$ eDR3 data, we find a self-consistent result only for potential models with a dynamical bar length of $\approx$3.5 kpc and pattern speed of 39 km/s/kpc. We find the Milky Way's trapped bar orbits extend out to only $\approx$3.5 kpc, but there is also an overdensity of stars at the end of the bar out to 4.8 kpc which could be related to an attached spiral arm. We also find that the measured orbital structure of the bar is strongly dependent on the properties of the assumed potential.

The conservation of wave action in moving plasmas has been well-known for over half a century. However, wave action is not conserved when multiple wave modes propagate and coexist close to degeneration condition (Sound speed equals Alfv\'en speed, i.e. plasma $\beta \sim 1$). Here we show that the violation of conservation is due to wave mode conversion, and that the total wave action summed over interacting modes is still conserved. Though the result is general, we focus on MHD waves and identify three distinctive mode conversion mechanisms, i.e. degeneracy, linear mode conversion, and resonance, and provide an intuitive physical picture for the mode conversion processes. We use 1D MHD simulations with the Expanding Box Model to simulate the nonlinear evolution of monochromatic MHD waves in the expanding solar wind. Simulation results validate the theory; total wave action therefore remains an interesting diagnostic for studies of waves and turbulence in the solar wind.

In November 2019 the nearby single, isolated DQ-type white dwarf LAWD 37 (LP 145-141) aligned closely with a distant background source and caused an astrometric microlensing event. Leveraging astrometry from Gaia and followup data from the Hubble Space Telescope we measure the astrometric deflection of the background source and obtain a gravitational mass for LAWD 37. The main challenge of this analysis is in extracting the lensing signal of the faint background source whilst it is buried in the wings of LAWD 37's point spread function. Removal of LAWD 37's point spread function induces a significant amount of correlated noise which we find can mimic the astrometric lensing signal. We find a deflection model including correlated noise caused by the removal of LAWD 37's point spread function best explains the data and yields a mass for LAWD 37 of $0.56\pm0.08 M_{\odot}$. This mass is in agreement with the theoretical mass-radius relationship and cooling tracks expected for CO core white dwarfs. Furthermore, the mass is consistent with no or trace amounts of hydrogen that is expected for objects with helium-rich atmospheres like LAWD 37. We conclude that further astrometric followup data on the source is likely to improve the inference on LAWD 37's mass at the $\approx3$ percent level and definitively rule out purely correlated noise explanations of the data. This work provides the first semi-empirical test of the white dwarf mass-radius relationship using a single, isolated white dwarf and supports current model atmospheres of DQ white dwarfs and white dwarf evolutionary theory.

We present an analysis of lifetimes and resonances of Earth Trojan Asteroids (ETAs) in the MEGASIM data set (Yeager & Golovich 2022). Trojan asteroids co-orbit the Sun with a planet but remain bound to the Lagrange points, L4 (60{\deg} leading the planet) or L5 (60{\deg} trailing). In the circular three-body approximation, the stability of a Trojan asteroid depends on the ratio of the host planet mass and the central mass. For the inner planets, the range of stability becomes increasingly small, so perturbations from the planets have made primordial Trojans rare. To date there have been just two ETAs (2010 TK7 and 2020 XL5), several Mars Trojans, and a Venus Trojan discovered. The estimated lifetimes of the known inner system Trojans are less than a million years, suggesting they are interlopers rather than members of a stable and long-lasting population. With the largest ETA n-body simulation to date, we are able to track their survival across a wide initialized parameter space. We find the remaining fraction of ETAs over time is well fit with a stretched exponential function that when extrapolated beyond our simulation run time predicts zero ETAs by 2.33 Gyr. We also show correlations between ETA ejections and the periods of the Milankovitch cycles. Though Earth's orbital dynamics dominate the instabilities of ETAs, we provide evidence that ETA ejections are linked to resonances found in the variation of the orbital elements of many, if not all of the planets.

We describe a simple extension to existing models for the tidal heating of dark matter subhalos which takes into account second order terms in the impulse approximation for tidal heating. We show that this revised model can accurately match the tidal tracks along which subhalos evolve as measured in high-resolution N-body simulations. We further demonstrate that, when a constant density core is introduced into a subhalo, this model is able to quantitatively reproduce the evolution and artificial disruption of N-body subhalos arising from finite resolution effects. Combining these results we confirm prior work indicating that artificial disruption in N-body simulations can result in a factor two underestimate of the subhalo mass function in the inner regions of host halos, and a 10-20% reduction over the entire virial volume.

The aim of this work is to evaluate the performance of photometric metallicity [Fe/H], determined based on V-band light-curves photometrically transformed from the gr-band light-curves. We tested this by using a set of homogeneous sample of fundamental mode RR Lyrae located in the Kepler field. It was found that the color-term is necessary in such photometric transformation. We demonstrated that including the color-term the determined photometric [Fe/H] are in good agreement to the spectroscopic [Fe/H], either based on the calibrated or the transformed V-band light-curves. We also tested the impact of Blazhko RR Lyrae in determining the photometric Fe/H], and found that Blazhko RR Lyrae can give consistent photometric [Fe/H]. Finally, we derived independent gVr-band [Fe/H]-$\phi_{31}$-P relations (where $\phi_{31}$ and P are Fourier parameter and pulsation period, respectively) using our light-curves. The V-band relation is in good agreement with the most recent determination given in the literature.

We use the {\it Gaia} EDR3 to explore the Galactic supernova remnant SNR G272.2-3.2, produced by the explosion of a Type Ia supernova (SNIa), about 7,500 years ago, to search for a surviving companion. From the abundances in the SNR ejecta, G272.2-3.2 is a normal SNIa. The {\it Gaia} parallaxes allow to select the stars located within the estimated distance range of the SNR, and the {\it Gaia} proper motions to study their kinematics. From the {\it Gaia} EDR3 photometry, we construct the HR diagram of the selected sample, which we compare with the theoretical predictions for the evolution of possible star companions of SNIa. We can discard several proposed types of companions by combining kinematics and photometry. We focus our study on the kinematically most peculiar star, {\it Gaia} EDR3 5323900215411075328 (hereafter MV-G272), is a very clear outlier in the motion along the Galactic plane, both as compared with the sample and with the Besan\c{c}on model of the Galaxy. Its trajectory on the sky locates it at the center of the SNR, 6,000--8,000 years ago, which constitutes a unique characteristic among the 3,082 stars in the sample. Spectra have been obtained allowing a stellar parameters determination and a chemical abundance analysis. The star now appears as an M1-M2 dwarf. Its chemistry is consistent with that of a M dwarf companion, where a deep convective envelope could have diluted the material captured from the SN ejecta. We have also explored a wide area around the SNR, looking for a possible hypervelocity star ejected by the SN explosion. None has been found. In conclusion, we have a candidate to be the surviving companion of the SNIa that resulted in SNR G272.2-3.2. It is supported by all its kinematical characteristics and by its trajectory within the SNR. This opens the possibility for SNe Ia originating in binary systems with a M-dwarf mass donor.

We present a method that derives dust temperatures and infrared (IR) luminosities of high-redshift galaxies assuming the radiation equilibrium in a simple dust and stellar distribution geometry. Using public data in the Atacama Large Millimeter/submillimeter Array (ALMA) archive, we studied dust temperatures in a clumpy interstellar medium (ISM) model for high-redshift galaxies, and tested consistency of the results with what derived from other methods. We find the dust distribution model assuming clumpiness of ${\rm log}\,\xi_{\rm clp}=-1.02\pm0.41$ can represent the ISM of high-redshift star-forming galaxies. By assuming the value of $\xi_{\rm{clp}}$, our method enables to derive dust temperatures and IR luminosities of high-redshift galaxies from only a single-band ALMA observation that gives two ALMA measurements required for the method: the dust continuum flux and the dust continuum emission size. Using this method, we demonstrate to find the dust temperature ($T_{\rm d}=95^{+13}_{-17}\,\rm{K}$) of a $z\sim8.3$ star-forming galaxy, MACS0416-Y1. The method only requires a single-band dust observation to derive a dust temperature. It is more easily accessible than multi-band observations or emission line searches at high-redshift. Thus can be applicable to large samples of galaxies in future studies using high-resolution interferometers such as ALMA.

We present the average [CII] $158\,\rm{\mu m}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{\mu m}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective radius $r_e$ of $2.2\pm0.2\,\rm{kpc}$. It is $\gtrsim2\times$ larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for $z\lesssim6$ galaxies. Additionally, we compared the average [CII] size with $4<z<6$ galaxies observed by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). By analysing [CII] sizes of $4<z<6$ galaxies in two redshift bins, we find an average [CII] size of $r_{\rm e}=2.2\pm0.2\,\rm{kpc}$ and $r_{\rm e}=2.5\pm0.2\,\rm{kpc}$ for $z\sim5.5$ and $z\sim4.5$ galaxies, respectively. These measurements show that star-forming galaxies, on average, show no evolution in the size of the [CII] $158\,{\rm \mu m}$ emitting regions at redshift between $z\sim7$ and $z\sim4$. This finding suggest that the star-forming galaxies could be morphologically dominated by gas over a wide redshift range.

We present an analysis of oscillation mode variability in the hot B subdwarf star EPIC~220422705, a new pulsator discovered from $\sim78$~days of {\em K}2 photometry. The high-quality light curves provide a detection of 66 significant independent frequencies, from which we identified 9 incomplete potential triplets and 3 quintuplets. Those {\sl g-} and {\sl p-}multiplets give rotation periods of $\sim$ 36 and 29 days in the core and at the surface, respectively, potentially suggesting a slightly differential rotation. We derived a period spacing of 268.5\,s and 159.4\,s for the sequence of dipole and quadruple modes, respectively. We characterized the precise patterns of amplitude and frequency modulations (AM and FM) of 22 frequencies with high enough amplitude for our science. Many of them exhibit intrinsic and periodic patterns of AM and FM, with periods on a timescale of months as derived by the best fitting and \texttt{MCMC} test. The nonlinear resonant mode interactions could be a natural interpretation for such AMs and FMs after other mechanisms are ruled out. Our results are the first step to build a bridge between mode variability from {\em K}2 photometry and nonlinear perturbation theory of stellar oscillation.

The very careful Event Horizon Telescope estimate of the mass of the supermassive black hole at the center of the Giant CD galaxy M87, allied with recent high quality photometric and spectroscopic measurements, yields a proper dark/luminous mass decomposition from the galaxy center to its virial radius. That provides us with decisive information on crucial cosmological and astrophysical issues. The dark and the standard matter distributions in a wide first time detected galaxy region under the supermassive black hole gravitational control. The well known supermassive black hole mass vs stellar dispersion velocity relationship at the highest galaxy masses implies an exotic growth of the former. This may be the first case in which one can argue that the supermassive black hole mass growth was also contributed by the Dark Matter component. A huge dark matter halo core in a galaxy with inefficient baryonic feedback is present and consequently constrains the nature of the dark halo particles. The unexplained entanglement between dark/luminous structural properties, already emerged in disk systems, also appears.

Changing-look active galactic nuclei (CLAGNs) show the disappearance and reappearance of broad emission lines in a few years, which challenges the orientation-based AGN unification model. We reduce the X-ray data for five well-studied CLAGNs that show a strong change in broad emission lines in the past several decades. We find that the X-ray photon index, $\Gamma$, and the Eddington-scaled X-ray luminosity, $L_{\rm 2-10 keV}/L_{\rm Edd}$, normally follow negative and positive correlations when the Eddington ratio is lower and higher than a critical value of $\sim 10^{-3}$. We find that the CLAGNs observed with broad H$\beta$ emission lines stay in the positive part of the $\Gamma-L_{\rm 2-10 keV}/L_{\rm Edd}$ correlation, while the broad H$\beta$ lines become weak or disappear in the anticorrelation part of the $\Gamma-L_{\rm 2-10 keV}/L_{\rm Edd}$ correlation, which suggests that the evolution of the broad lines should be correlated with the evolution of the underlying accretion process. We further find that the CLAGNs are consistent with the other different types of AGNs in the $L_{\rm bol}-L_{\rm bol}/L_{\rm Edd}$ correlation. These results support that the CLAGNs are belong to a special stage of AGNs with a bolometric Eddington ratio $\sim$1\%, where the broad emission lines are easily affected by the strong variation in ionization luminosity that is caused by the transition of accretion modes.

Galaxy clusters have the potential to accelerate cosmic rays (CRs) to ultra-high energies via accretion shocks or embedded CR acceleration sites. CRs with energies below the Hillas condition will be confined within the cluster and will eventually interact with the intracluster medium (ICM) gas to produce secondary neutrinos and $\gamma$ rays. Using 9.5 years of muon-neutrino track events from the IceCube Neutrino Observatory, we report the results of a stacking analysis of 1094 galaxy clusters, with masses $\gtrsim 10^{14}$ \(\textup{M}_\odot\) and redshifts between 0.01 and $\sim$1, detected by the {\it Planck} mission via the Sunyaev-Zeldovich (SZ) effect. We find no evidence for significant neutrino emission and report upper limits on the cumulative unresolved neutrino flux from galaxy clusters after accounting for the completeness of the catalog, assuming three different weighting scenarios for the stacking and three different power-law spectra. Weighting the sources according to mass and distance, we set upper limits at $90\%$ confidence level that constrain the flux of neutrinos from galaxy clusters to be no more than $4.6\%$ of the diffuse IceCube observations at 100 TeV, assuming an unbroken $E^{-2.5}$ power-law spectrum.

In its long-duration observation phase, the PLATO satellite will observe two non-overlapping fields for a total of 4 yr. The exact duration of each pointing will be determined 2 yr before launch. Previous estimates of PLATO's yield of Earth-sized planets in the habitable zones (HZs) around solar-type stars ranged between 6 and 280. We use the PLATO Solar-like Light curve Simulator (PSLS) to simulate light curves with transiting planets around bright (m_V > 11) Sun-like stars at a cadence of 25 s, roughly representative of the >15,000 targets in PLATO's high-priority P1 sample (mostly F5-K7 dwarfs and sub-dwarfs). Our study includes light curves generated from synchronous observations of 6, 12, 18, and 24 of PLATO's 12 cm aperture cameras over both 2 yr and 3 yr of continuous observations. Automated detrending is done with the Wotan software and post-detrending transit detection is performed with the Transit Least Squares (TLS) algorithm. We scale the true positive rates (TPRs) with the expected number of stars in the P1 sample and with modern estimates of the exoplanet occurrence rates and predict the detection of planets with 0.5 R_E <= R_p <= 1.5 R_E in the HZs around F5-K7 dwarf stars. For the (2 yr + 2 yr) long-duration observation phase strategy we predict 11-34 detections, for the (3 yr + 1 yr) strategy we predict 8-25 discoveries. Our study of the effects of stellar variability on shallow transits of Earth-like planets illustrates that our estimates of PLATO's planet yield, which we derive using a photometrically quiet star like the Sun, must be seen as upper limits. In conclusion, PLATO's detection of about a dozen Earth-sized planets in the HZs around solar-type stars will mean a major contribution to this yet poorly sampled part of the exoplanet parameter space with Earth-like planets.

The important role played by magnetic reconnection in the early acceleration of coronal mass ejections (CMEs) has been widely discussed. However, as CMEs may have expansion speeds comparable to their propagation speeds in the corona, it is not clear whether and how reconnection contributes to the true acceleration and expansion separately. To address this question, we analyze the dynamics of a moderately fast CME on 2013 February 27, associated with a continuous acceleration of its front into the high corona, even though its speed had reached $\sim$700~km~s$^{-1}$ and larger than the solar wind speed. The apparent CME acceleration is found to be due to the CME expansion in the radial direction. The CME true acceleration, i.e., the acceleration of its center, is then estimated by taking into account the expected deceleration caused by the solar wind drag force acting on a fast CME. It is found that the true acceleration and the radial expansion have similar magnitudes. We find that magnetic reconnection occurs after the CME eruption and continues during the CME propagation in the high corona, which contributes to the CME dynamic evolution. Comparison between the apparent acceleration related to the expansion and the true acceleration that compensates the drag shows that, for this case, magnetic reconnection contributes almost equally to the CME expansion and to the CME acceleration. The consequences of these measurements for the evolution of CMEs as they transit from the corona to the heliosphere are discussed.

Highly variable Markarian 421 is a bright high synchrotron energy peaked blazar showing wide featureless non-thermal spectrum making it a good candidate for our study of intraday flux and spectral variations over time. We analyse its X-ray observations of over 17 years taken with the EPIC-PN instrument to probe into the intraday variability properties. The photon energy band of 0.3 - 10.0 keV, and its sub-bands, soft 0.3-2.0 keV and hard 2.0-10.0 keV. To examine flux variability, fractional variability amplitude and the minimum variability timescale have been calculated. We also probed into the spectral variability by studying hardness ratio for each observation and the correlation between the two energy bands using discrete correlation function and inspecting the normalized light curves. The parameters obtained from these methods have been studied for any correlation or non-random trends. From this work, we speculate on the constraints for possible particle acceleration and emission processes in the jet, for better understanding of the processes involving a turbulent behaviour except of shocks. A positive discrete correlation function between the two sub-bands indicates the role of the same electron population in the emission of photons in the two bands. The correlation between the parameters of flux variability and parameters of spectral variation and lags in sub-energy bands provide the constraints to be considered for any modelling of emission processes.

The orbital period of V Sge is decreasing at a rate which increased from dP/dt=-(4.11\pm 0.33)\times 10^{-10} in 1962 to -(5.44\pm 0.61)\times 10^{-10} in 2022. This implies that the mass trahsfer from the secondary component is accelerating. From the evidence based on the orbital period variations, combined with estimates of the mass loss from the system based on radio observations, it follows that (1) the mass transfer rate from the secondary component is larger than dM2/dt=-5\times 10^{-6}M{\odot }/yr, possibly as large as dM2/dt=-2.5\times 10^{-5}M{\odot }/yr, and (2) the mass loss rate from the primary component is dM1/dt=-4\times 10^{-7}M{\odot }/yr or larger. Close similarity of V Sge to binary Wolf-Rayet stars supports the model with primary component being a hot, evolved star loosing its mass. Several arguments are presented which exclude the alternative model with primary component being a white dwarf with an accretion disk.

We present an observational study of the plasma dynamics at the base of a solar coronal jet, using high resolution extreme ultraviolet imaging data taken by the Extreme Ultraviolet Imager on board Solar Orbiter, and by the Atmospheric Imaging Assembly on board Solar Dynamics Observatory. We observed multiple plasma ejection events over a period of $\sim$1 hour from a dome-like base that is ca.~4 Mm wide and is embedded in a polar coronal hole. Within the dome below the jet spire, multiple plasma blobs with sizes around 1--2 Mm propagate upwards to dome apex with speeds of the order of the sound speed (ca.~120~km~s$^{-1}$ ). Upon reaching the apex, some of these blobs initiate flows with similar speeds towards the other footpoint of the dome. At the same time, high speed super-sonic outflows ($\sim$230~km~s$^{-1}$) are detected along the jet spire. These outflows as well as the intensity near the dome apex appear to be repetitive. Furthermore, during its evolution, the jet undergoes many complex morphological changes including transitions between the standard and blowout type eruption. These new observational results highlight the underlying complexity of the reconnection process that powers these jets and also provide insights into the plasma response when subjected to rapid energy injection.

Kilonovae, one source of electromagnetic emission associated with neutron star mergers, are powered by the decay of radioactive isotopes in the neutron-rich merger ejecta. Models for kilonova emission consistent with available modeling and the electromagnetic counterpart to GW170817 also predict characteristic abundance patterns, determined by the relative balance of different types of material in the outflow. Assuming the observed source is prototypical, this inferred abundance pattern in turn must match \emph{r}-process abundances deduced by other means, such as what is observed in the solar system. We report on analysis comparing the input mass-weighted elemental compositions adopted in our radiative transfer simulations to the mass fractions of elements in the Sun. We characterise the extent to which our parameter inference results depend on our assumed composition for the dynamical and wind ejecta and examine how the new results compare to previous work. We find that a mass ratio of $M_w/M_d$ = 2.81 reproduces the observed AT2017gfo kilonova light curves while also producing the abundance of neutron-capture elements in the solar system.

We investigated the decades' long-term X-ray variations in bright low-mass X-ray binaries containing a neutron star (NS-LMXB). The light curves of MAXI/GSC and RXTE/ASM covers $\sim$ 26 yr, and high-quality X-ray light curves are obtained from 33 NS-LMXBs. Among them, together with Ginga/ASM, two sources (GX 3$+$1 and GX 9$+$1) showed an apparent sinusoidal variation with the period of $\sim 5$ yr and $\sim 10$ yr in the 34 yr light curve. Their X-ray luminosities were $(1-4)\times10^{37}$ erg s$^{-1}$ in the middle of the luminosity distribution of the NS-LMXB. Other seven sources (Ser X-1, 4U 1735--444, GX 9$+$9, 4U 1746$-$37, 4U 1708$-$40, 4U 1822$-$000, and 1A 1246$-$588) have also similar sinusoidal variation, although the profiles (amplitude, period, and phase) are variable. Compering the 21 sources with known orbital periods, a possible cause of the long-term sinusoidal variation might be the mass transfer cycles induced by the irradiation to the donor star.

We performed wave-optics-based numerical simulations at mid-infrared wavelengths to investigate how the presence or absence of entrance slits and optical aberrations affect the spectral resolving power $R$ of a compact, high-spectral-resolving-power spectrometer containing an immersion-echelle grating. We tested three cases of telescope aberration (aberration-free, astigmatism and spherical aberration), assuming the aberration budget of the Space Infrared Telescope for Cosmology and Astrophysics (SPICA), which has a 20-$\mathrm{\mu m}$-wavelength diffraction limit. In cases with a slit, we found that the value of $R$ at around 10--20 $\mathrm{\mu m}$ is approximately independent of the assumed aberrations, which is significantly different from the prediction of geometrical optics. Our results also indicate that diffraction from the slit improves $R$ by enlarging the effective illuminated area on the grating window and that this improvement decreases at short wavelengths. For the slit-less cases, we found that the impact of aberrations on $R$ can be roughly estimated using the Strehl ratio.

We present an analysis of the subpulse drift in PSR J1750-3503, which is characterized by abrupt transitions of drift direction. As the pulsar does not exhibit other mode changes or clear nulling, it is an ideal candidate system for studying the phenomenon of drift direction change. For $\sim 80\%$ of the time the subpulses are characterized by positive drift - from early to later longitudes - while the drift direction is negative in the other $\sim 20\%$. The subpulse separation for single pulses with positive drift, $P_2=(18.8\pm 0.1)^{\circ}$, is higher then for single pulses with negative drift, $P_2=(17.5\pm 0.2)^{\circ}$. When the drift is stable, the measured repetition time of the drift pattern is $P_3^{\rm obs}=(43.5 \pm 0.4) P$, where $P$ is pulsar period. We show that the observed data can be reproduced by a carousel models with subpulse rotation around the magnetic axis using purely dipolar configuration of surface magnetic field. The observed drift characteristics can be modeled assuming that the actual repetition time $P_3<2P$, such that we observe its aliased value. A small variation in $P_3$, of the order of $6\%$ (or less assuming higher alias orders), is enough to reproduce the characteristic drift direction changes we observe.

Weak lensing studies typically require excellent seeing conditions for the purpose of maximizing the number density of well-resolved galaxy images. It is interesting to ask to what extent the seeing size limits the usefulness of the astronomical images in weak lensing. In this work, we study this issue with the data of the DECam Legacy Survey (DECaLS), which is a part of the target selection program for the Dark Energy Spectroscopic Instrument (DESI). Using the Fourier Quad shear measurement pipeline, we demonstrate that images with relatively poor seeing conditions (around 1.5 arcsec) can still yield accurate shear estimators. We do not find any correlation between systematic shear error and the image resolution.

The unrivalled astrometric and photometric capabilities of the Gaia mission have given new impetus to the study of young stars: both from an environmental perspective, as members of comoving star-forming regions, and from an individual perspective, as targets amenable to planet-hunting direct-imaging observations. In view of the large availability of theoretical evolutionary models, both fields would benefit from a unified framework that allows a straightforward comparison of physical parameters obtained by different stellar and substellar models. Methods. To this aim, we developed MADYS, a flexible Python tool for age and mass determination of young stellar and substellar objects. In this first release, MADYS automatically retrieves and cross-matches photometry from several catalogs, estimates interstellar extinction, and derives age and mass estimates for individual objects through isochronal fitting. Harmonizing the heterogeneity of publicly-available isochrone grids, the tool allows to choose amongst 16 models, many of which with customizable astrophysical parameters, for a total of $\sim$ 110 isochrone grids. Several dedicated plotting function are provided to allow an intuitive visual perception of the numerical output. After extensive testing, we have made the tool publicly available (https://github.com/vsquicciarini/madys). We demonstrate here the capabilities of madys, summarising already published results as well providing several new examples.

Aims. We focus on how to homogenize the comparison between galaxy samples having different characteristics. We check the projections of the fundamental metallicity relation (FMR) and the evolution of these projections between a sample selected at $z\sim0$ (SDSS) and $z\sim0.7$ (VIPERS). We check, in particular, whether and to what extent selection criteria can affect the results. Methods. We checked the influence of different biases introduced either by physical constraints (evolution of the luminosity function and differences in the fraction of blue galaxies) or data selection (the signal-to-noise ratio and quality of the spectra) on the FMR and its projections. To separate the differences occurring due to the physical evolution of galaxies with redshift from the false evolution mimed by these biases, we first analyzed the effects of these biases individually on the SDSS sample, and next, starting from the SDSS data, we built a VIPERS-equivalent $z \sim 0$ sample, replicating the main characteristics of VIPERS sample at $z\sim0.7$ for a fair comparison. Results. We found that the FMR projections are all sensitive to biases introduced by the selection on S/N and the quality flags of the emission line measurements in the spectra, especially the $\left[\text{O{\,\sc{iii}}}\right]\lambda 4959$ line. The exception is the metallicity versus the sSFR plane which is insensitive to these biases. The evolution of the luminosity function introduces a bias only in the plane metallicity versus the star formation rate (SFR) while the fraction of blue galaxies has no impact on results.

We present results from quasi-simultaneous multiwavelength observations of the Galactic black hole X-ray transient MAXI J1820$+$070 during the decay of the 2018 outburst and its entire subsequent mini-outburst in March 2019. We fit the X-ray spectra with phenomenological and Comptonizaton models and discuss the X-ray spectral evolution comparing with the multiwavelength behaviour of the system. The system showed a rebrightening in UV/Optical/NIR bands 7-days after the soft-to-hard transition during the main outburst decay while it was fading in X-rays and radio. In contrast, the mini-outburst occurred 165-days after the hard state transition of the initial outburst decay and was detected in all wavelengths. For both events, the measured timescales are consistent with those observed in other black hole systems. Contemporaneous hard X-ray/soft $\gamma$-ray observations indicate a non-thermal electron energy distribution at the beginning of the UV/Optical/NIR rebrightening, whereas a thermal distribution can fit the data during the hard mini-outburst activity. The broadband spectral energy distributions until the rebrightening are consistent with the irradiated outer accretion disc model. However, both the SEDs produced for the peak of rebrightening and close to the peak of mini-outburst provided good fits only with an additional power-law component in the UV/Optical/NIR frequency ranges which is often interpreted with a jet origin.

Solar jets are ubiquitous transient collimated mass outflows in the solar atmosphere over a wide range of sizes from small scale nanojets to a few solar radii, embedded in the solar chromosphere to solar corona. Jets are frequently accompanied by solar flares and these flares provide the force to propagate the plasma material upward and could be accompanied by coronal mass ejections. These jets could act as a source for transporting a significant mass and energy from the lower solar atmosphere to the upper coronal heights and consequently heating the solar corona and accelerating the solar wind. Magnetic reconnection is believed to be the triggering reason behind these jet activity. The thesis entitled: Study of Solar Jets and Related Flares, includes various case studies with different mechanisms to set off the jet initiation, associated large scale eruptions and mounts strong observational evidences to validate the numerical experiments for the magnetic flux emergence models. Such studies on solar jets along with their magnetic origin contribute to resolve the scandalous coronal heating problem and provide the evidences for the existing theoretical models and open a new window for the interplanetary science.

We present the results of a new stage of the long-term photometric study of FG Sge which is a quickly evolving central star of the planetary nebula Hen 1-5. Our new observations carried out on the SAI MSU telescopes in the optical ($BVR_CI_C$) and infrared (IR) ($JHKLM$) regions in 2008-2021 and 2013-2021, respectively, allowed us to trace the evolution of the star's brightness in recent years. The most significant observations were performed in 2019 when the star suffered a short clearing of the dust shell and became visible in $BVR_C$. Based on the spectral energy distribution of FG Sge in the 0.4-5 $\mu$m wavelength range we derived the dust shell parameters: the size of dust grains $a=0.01\mu$m, the inner radius temperature $T_{\text{dust}}=900$ K, optical depth $\tau (K)=0.5$ ($\tau (V)=4.5$), the total mass of dust $M_{\text{dust}}=7\cdot10^{-5} M_{\odot}$. After the short-term clearing of the dust shell in 2019, another dust structure was ejected that resulted in the star fading in all the observed bands. Based on the IR brightness and color curves, we estimated the dust depth growth in 2019-2020.

The Gaia spacecraft presents an unprecedented opportunity to reveal the population of long period ($a>1\,au$) exoplanets orbiting stars across the H-R diagram, including white dwarfs. White dwarf planetary systems have played an important role in the study of planetary compositions, from their unique ability to provide bulk elemental abundances of planetary material in their atmospheres. Yet, very little is known about the population of planets around white dwarfs. This paper predicts the population of planets that Gaia will detect around white dwarfs, evolved from known planets orbiting main-sequence stars. We predict that Gaia will detect $6\pm1$ planets around white dwarfs: $8\pm\,3\%$ will lie inside $3\,au$ and $30\pm10\,\%$ will be less massive than Jupiter. As surviving planets likely become dynamically detached from their outer systems, those white dwarfs with Gaia detected planets may not have planetary material in their atmospheres. Comparison between the predicted planet population and that found by Gaia will reveal the importance of dynamical instabilities and scattering of planets after the main-sequence, as well as whether photoevaporation removes the envelopes of gas giants during their giant branch evolution

The fact that the spatial velocity of pulsars is generally higher than that of their progenitor stars has bothered astronomers for nearly 50 years. It has been extensively argued that the high pulsar velocity should be acquired during a natal kick process on a timescale of 100ms - 10s in the supernova explosion, in which some asymmetrical dynamical mechanism plays a key role. However, a satisfactory picture generally is still lacking. In this study, it is argued that the neutrino rocket model can well account for the high speed as well as the long-term evolution behaviors of pulsars. The neutrinos are emitted from superfluid vortex neutrons through the neutrino cyclotron radiation mechanism. The unique characters of left-handed neutrinos and right-handed antineutrinos resulting from the nonconservation of parity in weak interactions play a major role in the spatial asymmetry. The continuous acceleration of pulsars can be naturally explained by this model, which yields a maximum velocity surpassing 1000 km s$^{-1}$. The alignment between the spinning axis and the direction of motion observed for the Crab pulsar (PSR 0531) and the Vela pulsar (PSR 0833) can be well accounted for. The observed correlation between the spin-down rate and the period of long-period pulsars with $P \gtrsim 0.5$s can also be satisfactorily explained.

We propose a novel model of the modified (Starobinsky-like) old-minimal-type supergravity coupled to a chiral matter superfield, that can {\it simultaneously} describe multi-field inflation, primordial black hole (PBH) formation, dark matter (DM), and spontaneous supersymmetry (SUSY) breaking after inflation in a Minkowski vacuum. The PBH masses in our supergravity model of double slow-roll inflation, with a short phase of "ultra-slow-roll" between two slow-roll phases, are close to $10^{18}$ g. We find that a significant PBH fraction in the allowed mass window requires the very high SUSY breaking scale with the gravitino mass close to the scalaron (inflaton) mass $M$ of the order $10^{13}$ GeV. Our supergravity model favors the {\it composite} nature of DM as a mixture of PBH and heavy gravitinos as the lightest SUSY particles. The composite DM significantly relaxes fine-tuning needed for the whole PBH-DM. The PBH-DM fraction is derived, and the second-order gravitational wave background induced by the enhanced scalar perturbations is calculated. Those gravitational waves may be accessible by the future space-based gravitational interferometers.

[Abridged] We map the stellar age distribution ($\lesssim 1$ Gyr) across a 6kpc$\,\times\,$6kpc area of the Galactic disc to constrain our Galaxy's recent star-formation history. Our modelling draws on the sample of Zari et al. (2021) that encompasses all presumed disc OBA stars ($\sim 500,000$ sources) with $G<16$. To be less sensitive to reddening, we do not forward model the detailed CMD distribution of these stars, but instead the K-band absolute magnitude distribution, $n(M_K)$, among stars with $M_K<0$ and $T_{\mathrm{eff}} > 7000$K at a certain positions $\vec{x}$ in the disc as a step function with five age bins, $b(\tau~|~\vec{x}, \vec{\alpha})$, logarithmically-spaced in age from $\tau = 5$ Myr to $\tau = 1$ Gyr. Given a set of isochrones and a Kroupa (2001) initial mass function, we sample $b(\tau\,|\,\vec{x}, \vec{\alpha})$ to maximize the likelihood of the data, accounting for the selection function. This results in a set of mono-age stellar density maps across a sizeable portion of the Galactic disc. These maps show that some, but not all, spiral arms are reflected in overdensities of stars younger than 500 Myr. The maps of the youngest stars ($<10$ Myr) trace major star forming regions. The maps at all ages exhibit an outward density gradient and distinct spiral-like spatial structure, which is qualitatively similar on large scales among all age intervals. When summing over the maps' area and extrapolating to the whole disc, we find an effective star-formation rate over the last 10 Myr of $\approx 3.3 \mathrm{M_{\odot}/yr}$, higher than previously published estimates. Remarkably, our stellar age distribution implies that the star-formation rate has been three times lower throughout most of the last Gyr, having risen distinctly only in the very recent past.

Molecular hydrogen is the most abundant molecule in the Galaxy and plays important roles for planets, their circumstellar environments, and many of their host stars. We have confirmed the presence of molecular hydrogen in the AU Mic system using high-resolution FUV spectra from HST-STIS during both quiescence and a flare. AU Mic is a $\sim$23 Myr M dwarf which hosts a debris disk and at least two planets. We estimate the temperature of the gas at 1000 to 2000 K, consistent with previous detections. Based on the radial velocities and widths of the H$_2$ line profiles and the response of the H$_2$ lines to a stellar flare, the H$_2$ line emission is likely produced in the star, rather than in the disk or the planet. However, the temperature of this gas is significantly below the temperature of the photosphere ($\sim$3650 K) and the predicted temperature of its star spots ($\gtrsim$2650 K). We discuss the possibility of colder star spots or a cold layer in the photosphere of a pre-main sequence M dwarf.

Cosmic rays represent one of the most important energy transformation processes of the universe. They bring information about the surrounding universe, our galaxy, and very probably also the extragalactic space, at least at the highest observed energies. More than one century after their discovery, we have no definitive models yet about the origin, acceleration and propagation processes of the radiation. The main reason is that there are still significant discrepancies among the results obtained by different experiments, probably due to some still unknown systematic uncertainties affecting the measurements. In this paper, we will focus on the detection of galactic cosmic rays in the 10$^{15}$ eV energy range, where the so-called \emph{`knee'} in the all-particle energy spectrum is observed. The measurement of the (p+He) energy spectrum is presented and discussed.

Deuterium fractionation can constrain the physical and chemical conditions at the early stage of brown dwarf formation. We present IRAM 30m observations over a wide frequency range of 213-279 GHz of singly and doubly deuterated species of formaldehyde (HDCO and D$_{2}$CO) towards Class 0/I proto-brown dwarfs (proto-BDs). Multiple low-excitation HDCO and D$_{2}$CO transition lines with upper energy level $\leq$ 40 K are detected. The D$_{2}$CO/HDCO, HDCO/H$_{2}$CO, D$_{2}$CO/H$_{2}$CO abundance ratios range between 0.01 and 2.5 for the proto-BDs, similar to the range seen in low-mass protostars. The highest ratios of D$_{2}$CO/HDCO $\sim$1.3-2.5 are measured for two Stage 0 proto-BDs. These objects could possess a warm corino, similar to the few hot corino cases reported among Class 0 protostars. The mean D$_{2}$CO/HDCO, D$_{2}$CO/H$_{2}$CO, and HDCO/H$_{2}$CO ratios for the proto-BDs are comparatively higher than the range predicted by the current gas-grain chemical models, indicating that HDCO and D$_{2}$CO are formed via grain surface reactions in the dense and cold interiors of the proto-BDs at an early formation stage.

Primordial black holes are under intense scrutiny since the detection of gravitational waves from mergers of solar-mass black holes in 2015. More recently, the development of numerical tools and the precision observational data have rekindled the effort to constrain the black hole abundance in the lower mass range, that is $M < 10^{23}$g. In particular, primordial black holes of asteroid mass $M \sim 10^{17}-10^{23}\,$g may represent 100\% of dark matter. While the microlensing and stellar disruption constraints on their abundance have been relieved, Hawking radiation of these black holes seems to be the only detection (and constraining) mean. Hawking radiation constraints on primordial black holes date back to the first papers by Hawking. Black holes evaporating in the early universe may have generated the baryon asymmetry, modified big bang nucleosynthesis, distorted the cosmic microwave background, or produced cosmological backgrounds of stable particles such as photons and neutrinos. At the end of their lifetime, exploding primordial black holes would produce high energy cosmic rays that would provide invaluable access to the physics at energies up to the Planck scale. In this review, we describe the main principles of Hawking radiation, which lie at the border of general relativity, quantum mechanics and statistical physics. We then present an up-to-date status of the different constraints on primordial black holes that rely on the evaporation phenomenon, and give, where relevant, prospects for future work. In particular, non-standard black holes and emission of beyond the Standard Model degrees of freedom is currently a hot subject.

An understanding of the dynamical evolution of binary star systems, and their effects on stellar and planetary evolution, requires well-characterized binary populations across stellar ages. However, the observational resources required to find and characterize binaries are expensive. With the release of high-precision Gaia astrometry, the re-normalized unit weight error (RUWE) statistic has been shown to reveal the presence of binary systems, with RUWE values greater than 1.2 indicating the presence of a stellar companion within $ \sim 1$". Our goal is to assess whether this new diagnostic, which was developed for field-age systems ($>$1 Gyr), applies to young systems; specifically, those that host circumstellar disks. With a control sample of single-star systems, compiled from high-contrast imagining surveys of the Taurus and Upper Scorpius star-forming regions, we compare the RUWE values for systems with and without circumstellar disks. We show that the presence of a protoplanetary disk alone can result in inflated RUWE values. Based on the distribution of the RUWE for disk-bearing single stars, we suggest a more conservative single-star -- binary threshold is warranted in the presence of disk material. We place this cutoff at the distribution's 95th percentile, with RUWE $= 2.5$.

Results of the mathematical treatment of the radiation by the superluminally moving current sheet in the magnetosphere of a neutron star, which was presented in Ardavan (2021, MNRAS, 507, 4530), are explained here in more transparent physical terms with the aid of illustrations. Not only do these results provide an all-encompassing explanation for the salient features of the radiation received from pulsars (its brightness temperature, polarization, spectrum and profile with microstructure and with a phase lag between the radio and gamma-ray peaks), but they also shed light on the putative energetic requirements of gamma-ray pulsars and magnetars and the sources of fast radio bursts and gamma-ray bursts.

Cosmic-ray accelerators capable of reaching ultra-high energies are expected to also produce very-high energy neutrinos via hadronic interactions within the source or its surrounding environment. Many of the candidate astrophysical source classes are either transient in nature or exhibit flaring activity. Using the Earth as a neutrino converter, suborbital and space-based optical Cherenkov detectors, such as POEMMA and EUSO-SPB2, will be able to detect upward-moving extensive air showers induced by decaying tau-leptons generated from cosmic tau neutrinos with energies $\sim 10$ PeV and above. Both EUSO-SPB2 and POEMMA will be able to quickly repoint, enabling rapid response to astrophysical transient events. We calculate the transient sensitivity and sky coverage for both EUSO-SPB2 and POEMMA, accounting for constraints imposed by the Sun and the Moon on the observation time. We also calculate both detectors' neutrino horizons for a variety of modeled astrophysical neutrino fluences. We find that both EUSO-SPB2 and POEMMA will achieve transient sensitivities at the level of modeled neutrino fluences for nearby sources. We conclude with a discussion of the prospects of each mission detecting at least one transient event for various modeled astrophysical neutrino sources.

One of the main pillars of the {\Lambda}CDM model is the Cosmological Principle, which states that our Universe is statistically isotropic and homogeneous on large scales. Here we test this hypothesis using the Astrophysical Gravitational Wave Background (AGWB) expected to be measured by the Einstein Telescope-Cosmic Explorer network; in particular we perform a numerical computation of the AGWB dipole, evaluating the intrinsic contribution due to clustering and the kinematic effect induced by the observer motion. We apply a component separation technique in the GW context to disentangle the kinematic dipole, the intrinsic dipole and the shot noise (SN), based on the observation of the AGWB at different frequencies. We show how this technique can also be implemented in matched-filtering to minimize the covariance which accounts for both instrumental noise and SN. Since GW detectors are essentially full-sky, we expect that this powerful tool can help in testing the isotropy of our Universe in the next future.

We carry out a classification of the glitch amplitudes of radio pulsars using Extreme Deconvolution based Gaussian Mixture Model, where the observed uncertainties in the glitch amplitude $\left(\frac{\Delta \nu}{\nu}\right)$ are taken into account. We then use information theory criteria such as AIC and BIC to determine the optimum number of glitch classes. We find that both AIC and BIC show that the pulsar glitch amplitudes can be optimally described using a bimodal distribution. The mean values of the fractional glitch amplitude for the two components are equal to $5.9 \times 10^{-9}$ and $1.3 \times 10^{-6}$, respectively. The unified dataset and analysis codes used in this work have been made publicly available.

We investigate the warm ionized gas in the blue compact galaxy (BCG) Haro 14 by means of integral field spectroscopic observations taken with the MUSE/VLT. The large FoV of MUSE and its unprecedented sensitivity enable observations of the galaxy nebular emission up to large galactocentric distances. This allowed us to trace the ionized gas morphology and ionization structure up to kiloparsec scales and, for the first time, to accurately investigate the excitation mechanism operating in the outskirts of a typical BCG. The intensity and diagnostic maps reveal at least two highly distinct components of ionized gas: the bright central regions, mostly made of individual clumps, and a faint component which extends up to kiloparsec scales and consists of widespread diffuse emission, well-delineated filamentary structures, and faint knots. Noteworthy are the two curvilinear filaments extending up to 2 and 2.3 kpc southwest, which likely trace the edges of supergiant expanding bubbles driven by galactic outflows. We find that while the central clumps in Haro 14 are HII-region complexes, the morphology and line ratios of the whole low-surface-brightness component are not compatible with star formation photoionization. In the spatially resolved emission-line-ratio diagnostic diagrams, spaxels above the maximum starburst line form the majority. Moreover, our findings suggest that more than one alternative mechanism is ionizing the outer galaxy regions. The properties of the diffuse component are consistent with ionization by diluted radiation and the large filaments and shells are most probably shocked areas at the edge of bubbles. The mechanism responsible for the ionization of the faint individual clumps observed in the galaxy periphery is more difficult to assess. These clumps could be the shocked debris of fragmented shells or regions where star formation is proceeding under extreme conditions.

Deviations from the blackbody spectral energy distribution of the CMB are a precise probe of physical processes active both in the early universe (such as those connected to particle decays and inflation) and at later times (e.g. reionization and astrophysical emissions). Limited progress has been made in the characterization of these spectral distortions after the pioneering measurements of the FIRAS instrument on the COBE satellite in the early 1990s, which mainly targeted the measurement of their average amplitude across the sky. Since at present no follow-up mission is scheduled to update the FIRAS measurement, in this work we re-analyze the FIRAS data and produce a map of $\mu$-type spectral distortion across the sky. We provide an updated constraint on the $\mu$ distortion monopole $|\langle\mu\rangle|<47\times 10^{-6}$ at 95\% confidence level that sharpens the previous FIRAS estimate by a factor of $\sim 2$. We also constrain primordial non-Gaussianities of curvature perturbations on scales $10\lesssim k\lesssim 5\times 10^4$ through the cross-correlation of $\mu$ distortion anisotropies with CMB temperature and, for the first time, the full set of polarization anisotropies from the Planck satellite. We obtain upper limits on $f_{\rm NL}\lesssim 3.6 \times 10^6$ and on its running $n_{\rm NL}\lesssim 1.4$ that are limited by the FIRAS sensitivity but robust against galactic and extragalactic foreground contaminations. We revisit previous similar analyses based on data of the Planck satellite and show that, despite their significantly lower noise, they yield similar or worse results to ours once all the instrumental and astrophysical uncertainties are properly accounted for. Our work is the first to self-consistently analyze data from a spectrometer and demonstrate the power of such instrument to carry out this kind of science case with reduced systematic uncertainties.

We have conducted a search for strong gravitational lensing systems in the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys Data Release 9. This is the third paper in a series (following Huang et al. 2020; Huang et al. 2021, Paper I & II, respectively). These surveys together cover $\sim$ 19,000 deg$^2$ visible from the northern hemisphere, reaching a z-band AB magnitude of $\sim$ 22.5. We use a deep residual neural network, trained on a compilation of known lensing systems and candidates as well as non-lenses in the same footprint. After applying our trained neural networks to the survey data, we visually inspect and rank images with probabilities above a threshold. We have found 1895 lens candidates. Out of these, 1512 are identified for the first time. Combining the discoveries from this work, Paper I (335) and II (1210), the total number of strong lens candidates from the Legacy Surveys that we have discovered is 3057.

Chemical abundances of stars in the Milky Way disk are empirical tracers of its enrichment history. However, they capture joint-information that is valuable to disentangle. In this work, we seek to quantify how individual abundances evolve across the present-day radius of the disk, at fixed supernovae contribution ([Fe/H], [Mg/Fe]). We use 18,135 APOGEE DR17 red clump stars and 7,943 GALAH DR3 main sequence stars to compare the abundance distributions conditioned on ([Fe/H], [Mg/Fe]) across $3-13$ kpc and $6.5-9.5$ kpc, respectively. In total we examine 15 elements: C, N, Al, K (light), O, Si, S, Ca, ($\alpha$), Mn, Ni, Cr, Cu, (iron-peak) Ce, Ba (s-process) and Eu (r-process). We find that the conditional neutron capture and light elements most significantly trace variations in the disk's enrichment history, with absolute conditional radial gradients $\leq 0.03$ dex/kpc. The other elements studied have absolute conditional gradients $\lesssim 0.01$ dex/kpc. We uncover structured conditional abundance variations as a function of [Fe/H] for the low-$\alpha$, but not the high-$\alpha$ sequence. The average scatter between the mean conditional abundances at different radii is $\sigma_\text{intrinsic} \approx$ 0.02 dex (with Ce, Eu, Ba $\sigma_\text{intrinsic} >$ 0.05 dex). These results serve as a measure of the magnitude via which different elements trace Galactic radial enrichment history once fiducial supernovae correlations are accounted for. Furthermore, we uncover subtle systematic variations in all moments of the conditional abundance distributions that will presumably constrain chemical evolution models of the Galaxy.

The Event Horizon Telescope (EHT) collaboration recently released horizon-scale images of the supermassive black hole M87*. These images are consistently described by an optically thin, lensed accretion flow in the Kerr spacetime. General relativity (GR) predicts that higher-resolution images of such a flow would present thin, ring-shaped features produced by photons on extremely bent orbits. Recent theoretical work suggests that these "photon rings" produce clear interferometric signatures whose observation could provide a stringent consistency test of the Kerr hypothesis, with scant dependence on the astrophysical configuration. Gralla, Lupsasca and Marrone (GLM) argued that the shape of high-order photon rings follows a specific functional form that is insensitive to the details of the astrophysical source, and proposed an experimental method for measuring this GR-predicted shape via space-based interferometry. We wish to assess the robustness of their prediction by checking that it holds for a variety of astrophysical profiles, black hole spins and observer inclinations. We repeat their analysis for hundreds of models and identify the width of the photon ring and its angular variation as a main obstacle to their method's success. We qualitatively describe how this width varies with the emission profile, black hole spin and observer inclination. At low inclinations, an improved method is robust enough to confirm the shape prediction for a variety of emission profiles; however, the choice of baseline is critical to the method's success. At high inclinations, we encounter qualitatively new effects that are caused by the ring's non-uniform width and require further refinements to the method. We also explore how the photon ring shape could constrain black hole spin and inclination.

Cosmological evolution of axion field in the early universe might be significantly affected by a thermal friction induced by the axion coupling to thermalized hidden sector. We examine the effects of such a thermal friction on axion dark matter density and its perturbation when the thermal friction dominates over the Hubble friction until when the axion field begins to oscillate around the potential minimum. We show that in the presence of sizable thermal friction there can be an exponential decay phase of the axion field before the oscillation phase, during which the axion energy density is efficiently dissipated into hidden thermal bath. Consequently, the previously excluded parameter region due to overclosing relic axion density becomes cosmologically viable with thermal friction. In particular, a QCD axion much lighter than $\mu$eV is viable without tuning the initial misalignment angle. We also find that thermal friction can affect the density perturbation of axion dark matter in various ways. For instance, it can alleviate the large-scale isocurvature bound on axion dark matter in the pre-inflationary PQ breaking scenario, which would make the pre-inflationary axion dark matter compatible with high scale inflation over a wide range of model parameters. In the post-inflationary PQ breaking scenario, thermal friction can also significantly change the scaling behavior of axionic strings, and therefore the typical size of the resultant axion miniclusters.

Auroras can be regarded as the most fascinating manifestation of space weather and they are continuously observed by ground-based and, nowadays more and more, also by space-based measurements. Investigations of auroras and geospace comprise the main research goals of the Suomi 100 nanosatellite, the first Finnish space research satellite, which has been measuring Earth's ionosphere since its launch on Dec. 3, 2018. In this work, we present a case study where the satellite's camera observations of an aurora in Northern Europe are combined with ground-based observations of the same event. The analyzed image is, to authors' best knowledge, the first auroral image ever taken by a cubesat. Our data analysis shows that a satellite vantage point provides complementary, novel information of an aurora. The 3D auroral location reconstruction of the analyzed auroral event demonstrates how information from a 2D image can be used to provide location information of aurora under study. The location simulation also suggests that the Earth's limb direction, which was the case in the analyzed image, is an ideal direction to observe faint aurora. Overall, the data analysis and reconstruction simulations demonstrate how even a small 1-Unit (size: 10 cm x 10 cm x 10 cm) CubeSat and its camera, build with cheap commercial of-the-shelf components, can open new possibilities for auroral research, especially, when its measurements are combined with ground-based observations.

Radio waves provide a useful diagnostic tool to investigate the properties of the ionosphere because the ionosphere affects the transmission and properties of High Frequency (HF) electromagnetic waves. We have conducted a transionospheric HF-propagation research campaign with a nanosatellite on a low-Earth polar orbit and the EISCAT HF transmitter facility in Troms{\o}, Norway, in December 2020. In the active measurement, the EISCAT HF facility transmitted sinusoidal 7.953 MHz signal which was received with the HEARER radio spectrometer onboard 1 Unit (size: 10 cm x 10 cm x 10 cm) Suomi 100 space weather nanosatellite. Data analysis showed that the EISCAT HF signal was detected with the satellite's radio spectrometer when the satellite was the closest to the heater along its orbit. Part of the observed variations seen in the signal was identified to be related to the heater's antenna pattern and to the transmitted pulse shapes. Other observed variations can be related to the spatial and temporal variations of the ionosphere and its different responses to the used transmission frequencies and to the transmitted O- and X-wave modes. Some trends in the observed signal may also be associated to changes in the properties of ionospheric plasma resulting from the heater's electromagnetic wave energy. This paper is, to authors' best knowledge, the first observation of this kind of "self-absorption" measured from the transionospheric signal path from a powerful radio source on the ground to the satellite-borne receiver.

Neutrino physics in one of the most active fields of research with important implications for particle physics, cosmology and astrophysics. On the other hand, motivated by some theories including string theory, formulation of physical theories in more than four space-time dimensions has been the subject of increasing attention in recent years. Interaction of neutrinos with gravitational fields is one of the interesting phenomena which can lead to transition between different helicity states (spin oscillations). We study neutrino spin oscillations in Schwarzschild and RN backgrounds in higher dimensional gravitational fields. We calculate the transition probability as a function of time and also study the dependence of the oscillation frequency on the orbital radius. The results help us to better understand the behavior of gravity and neutrinos in higher dimensions.

Black holes are conjectured to be the fastest quantum scramblers in nature, with the stretched horizon being the scrambling boundary. Under this assumption, we show that any infalling body must couple to virtually the entire black hole Hilbert space even prior to the Page time in order for there to be any hope of preserving the often-cited claim of the equivalence principle that such bodies should experience `no drama' as they pass a black hole's horizon. Further, under the scrambling assumption, we recover the usual firewall result at the black hole's Page time for an initially pure-state black hole without the need for any complexity or computational assumptions. For a black hole that is initially impure, we find that the onset of the firewall is advanced to times prior to the standard Page time. Finally, if black holes really do efficiently scramble quantum information, this suggests that, in order to preserve this claim of the equivalence principle even prior to the onset of a full-blown firewall, the quantum state of a black hole interior must be a Bose-Einstein condensate.

Neutron star (NS) is regarded as the natural laboratory for nuclear physics. The equation of state (EoS) extracted in flat spacetime is used chronically as an input to the Tolman-Oppenheimer-Volkoff (TOV) equation to constrain the structure of NS. However, using such EoS to characterize the NS with obvious gravitational effect seems controversial. In our work, we demonstrate the EoS of the same nuclear matter, either on earth or inside NS, ought to be in the same form due to the relativity principle. Gravity only enhances the temperature and the chemical potential, known as Tolman's law and Klein's law. We also clarify the self-consistency of the TOV equation, i.e., the equilibrium thermodynamics and gravity are included uniformly. The reason for conclusions in JCAP 02, 026 (2021) and Phys. Rev. D 104, 123005 (2021) is that the equilibrium thermodynamic relations protected by the equivalence principle in local spacetime are not taken into account.

$R^2$-corrected dark energy (DE) models in $F(R)$ gravity have been widely investigated in recent years, which not only removes the weak singularity potentially present in DE models but also provide us with a unified picture of the cosmic history, including the inflationary and DE epochs. Towards the unified interpretation of dynamical DE all over the cosmic history in the class of $R^2$-corrected DE models, we explore the universal features of the scalaron dynamics in the radiation-dominated epoch, along with the chameleon mechanism, by keeping our eyes on the inflationary and DE epochs. We show that the scalaron evolution does not follow a {\it surfing solution} and is mostly adiabatic before big bang nucleosynthesis (BBN), even properly including the {\it kick} by the nonperturbative QCD phase transition, hence a catastrophic consequence claimed in the literature is not applied to this class of DE models. This is due to the presence of the gigantic scale hierarchy between $R^2$ correction and DE, so is the universal feature for the class of $R^2$-corrected DE models. The prospects for the post- or onset-inflationary epoch would be pretty different from what the standard $R^2$ inflationary scenario undergoes due to the presence of the chameleon mechanism.

This thesis focuses on a variety of active research topics, such as nuclear matter, neutron stars, and phase transition within the framework of the RMF model. We use the previously successful effective field theory-driven Relativistic Mean Field (RMF) and density-dependent RMF (DD-RMF)formalisms for analyzing hadron matter to examine the infinite nuclear matter and neutron stars. The presence of exotic phases such as quarks has been investigated using the MIT Bag model and its variants, such as the vBag model, at various bag constants. The other exotic phases, such as hyperons, have also been studied under the influence of a strong magnetic field.

This paper describes the use of the Murchison Widefield Array, a low-frequency radio telescope at a radio-quiet Western Australian site, as a radar receiver forming part of a continent-spanning multistatic radar network for the surveillance of space. This paper details the system geometry employed, the orbit-specific radar signal processing, and the orbit determination algorithms necessary to ensure resident space objects are detected, tracked, and propagated. Finally, the paper includes the results processed after a short collection campaign utilising several FM radio transmitters across the country, up to a maximum baseline distance of over 2500 km. The results demonstrate the Murchison Widefield Array is able to provide widefield and persistent coverage of objects in low Earth orbit.

MOG as a modified gravity theory is designed to be replaced with dark matter. In this theory, in addition to the metric tensor, a massive vector is a gravity field where each particle has a charge proportional to the inertial mass and couples to the vector field through the four-velocity of a particle. In this work, we present the Hamiltonian formalism for the dynamics of particles in this theory. The advantage of Hamiltonian formalism is a better understanding and analyzing the dynamics of massive and massless particles. The massive particles deviate from the geodesics of space-time and photons follow the geodesics. We also study the dynamics of particles in the Newtonian and post-Newtonian regimes for observational purposes. An important result of Hamiltonian formalism is that while lensing on large scales is compatible with the observations, however the deflection angle from stellar size lensing is larger than General Relativity. This result can rule out this theory unless we introduce a screening mechanism to change the effective gravitational constant near compact objects like stars.

The Event Horizon Telescope (EHT), recently released the image of supermassive black hole Sgr A* showing an angular shadow diameter $d_{sh}= 48.7 \pm 7\,\mu$as, with an inferred black hole mass $M = 4.0^{+1.1}_{-0.6} \times 10^6 M_\odot $ and Schwarzschild shadow deviation $\delta = -0.08^{+0.09}_{-0.09}~\text{(VLTI)},-0.04^{+0.09}_{-0.10}~\text{(Keck)}$. The EHT image of Sgr A* is consistent with a Kerr black hole's expected appearance and the results directly prove a supermassive black hole in the center of the Milky Way. The Kerr hypothesis, a strong-field prediction of general relativity (GR), may not hold in the theories of gravity that admit Kerr-like black holes having an additional deviation parameter arising from the underlying theory. Here, we use the EHT observational results of Sgr A* to investigate the constraints on the deviation parameter whereby, such a rotating Kerr-like black hole can be an astrophysical black hole candidate, paying attention to three leading models. Modelling Kerr-like black holes as supermassive black hole Sgr A*, we observe that for it to be a viable astrophysical black hole candidate, the EHT results of Sgr A* put more stringent constraints on the parameter space than those put by the EHT results of M87*. However, a systematic bias analysis shows Kerr-like black hole shadows may capture Kerr black hole shadows over a good part of the constrained parameter space, making Kerr-like and Kerr black holes indistinguishable and one can't rule out a possibility of potential modifications of the Kerr metric or GR.