Papers for Friday, Dec 16 2022

We investigate the main tensions within the current standard model of cosmology from the perspective of the void size function in BOSS DR12 data. For this purpose, we present the first cosmological constraints on the parameters $S_8\equiv \sigma_8\sqrt{\Omega_{\rm m}/0.3}$ and $H_0$ obtained from voids as a stand-alone probe. We rely on an extension of the popular volume-conserving model for the void size function, tailored to the application on data, including geometric and dynamic distortions. We calibrate the two nuisance parameters of this model with the official BOSS collaboration mock catalogs and propagate their uncertainty through the statistical analysis of the BOSS void number counts. We focus our analysis on the $\Omega_{\rm m}$--$\sigma_8$ and $\Omega_{\rm m}$--$H_0$ parameter planes and derive the marginalized constraints $S_8 = 0.78^{+0.16}_{-0.14}$ and $H_0=65.2^{+4.5}_{-3.6}$ $\mathrm{km} \ \mathrm{s}^{-1} \ \mathrm{Mpc}^{-1}$. Our estimate of $S_8$ is fully compatible with constraints from the literature, while our $H_0$ value slightly disagrees by more than $1\sigma$ with recent local distance ladder measurements of type Ia supernovae. Our results open up a new viewing angle on the rising cosmological tensions and are expected to improve notably in precision when jointly analyzed with independent probes.

Halo inhabitants are individual stars, stellar streams, star and globular clusters, and dwarf galaxies. Here we compare the two last categories that include objects of similar stellar mass, which are often studied as self-dynamical equilibrium systems. We discover that the half-light radius of globular clusters depends on their orbital pericenter and total energy, and that Milky Way (MW) tides may explain the observed correlation. We also suggest that the accretion epoch of stellar systems in the MW halo can be calibrated by the total orbital energy, and that such a relation is due to both the mass growth of the MW and dynamical friction affecting mostly satellites with numerous orbits. This calibration starts from the bulge, to Kraken, Gaia Sausage Enceladus, Sagittarius stellar systems, and finally to the new coming dwarfs, either or not linked to the vast-polar structure. The most eccentric globular clusters and dwarfs have their half-light radius scaling as the inverse of their binding energy, and this over more than two decades. This means that earlier arriving satellites are smaller due to the tidal effects of the MW. Therefore, most halo inhabitants appear to have their structural parameters shaped by MW tides and also by ram-pressure for the most recent arrivals, the dwarf galaxies. The correlations found in this study can be used as tools to further investigate the origin of globular clusters and dwarfs, as well as the assembly history of our Galaxy.

We conduct hydrodynamical MOND simulations of isolated disc galaxies over the stellar mass range $M_{\star}/M_\odot = 10^7 - 10^{11}$ using the adaptive mesh refinement code \textsc{phantom of ramses} (\textsc{por}), an adaptation of the \textsc{ramses} code with a Milgromian gravity solver. The scale lengths and gas fractions are based on observed galaxies, and the simulations are run for 5~Gyr. The main aim is to see whether existing sub-grid physics prescriptions for star formation and stellar feedback reproduce the observed main sequence and reasonably match the Kennicutt-Schmidt relation that captures how the local and global star formation rates relate to other properties. Star formation in the models starts soon after initialisation and continues as the models evolve. The initialized galaxies indeed evolve to a state which is on the observed main sequence, and reasonably matches the Kennicutt-Schmidt relation. The available formulation of sub-grid physics is therefore adequate and leads to galaxies that largely behave like observed galaxies, grow in radius, and have flat rotation curves $-$ provided we use Milgromian gravitation. Furthermore, the strength of the bars tends to be inversely correlated with the stellar mass of the galaxy, whereas the bar length strongly correlates with the stellar mass. Irrespective of the mass, the bar pattern speed stays constant with time, indicating that dynamical friction does not affect the bar dynamics. The models demonstrate Renzo's rule and form structures at large radii, much as in real galaxies. In this framework, baryonic physics is thus sufficiently understood to not pose major uncertainties in our modelling of global galaxy properties.

FG Sge has evolved from the hot central star of the young planetary nebula Hen 1-5 to a G-K supergiant in the last 100 years. It is one of the three born-again objects that has been identified as of yet, and they are considered to have undergone a thermal pulse in the post-asymptotic giant branch evolution. FG Sge was observed with MIDI at the Very Large Telescope Interferometer at baselines of 43 and 46 m between 8 and 13 micron. The circumstellar dust environment of FG Sge was spatially resolved, and the Gaussian fit to the observed visibilities results in a full width at half maximum of ~10.5 mas. The observed mid-infrared visibilities and the spectral energy distribution can be fairly reproduced by optically thick (tauV ~ 8) spherical dust shell models consisting of amorphous carbon with an inner radius Rin of ~30 Rstar (corresponding to a dust temperature of 1100 +/- 100 K). The dust shell is characterized with a steep density profile proportional to r^{-3.5+/-0.5} from the inner radius Rin to (5-10) x Rin, beyond which it changes to r^{-2}. The dust mass is estimated to be ~ 7 x 10^{-7} solar mass, which translates into an average total mass-loss rate of ~9 x 10^{-6} solar mass/yr as of 2008 with a gas-to-dust ratio of 200 being adopted. In addition, the 8--13 micron spectrum obtained with MIDI with a field of view of 200 mas does not show a signatureof the polycyclic aromatic hydrocarbon (PAH) emission, which is in marked contrast to the spectra taken with the Spitzer Space Telescope six and 20 months before the MIDI observations with wide slit widths of 3.6--10 arcsec. This implies that the PAH emission originates from an extended region of the optically thick dust envelope, i.e., from the material ejected before the central star became H-deficient.

Direct photometric measurements of the cosmic optical background (COB) provide an important point of comparison to both other measurement methodologies and models of cosmic structure formation, and permit a cosmic consistency test with the potential to reveal additional diffuse sources of emission. The COB has been challenging to measure from Earth due to the difficulty of isolating it from the diffuse light scattered from interplanetary dust in our solar system. We present a measurement of the COB using data taken by the Long-Range Reconnaissance Imager (LORRI) on NASA's New Horizons mission, considering all data acquired to 47 AU. We employ a blind methodology where our analysis choices are developed against a subset of the full data set, which is then unblinded. Dark current and other instrumental systematics are accounted for, including a number of sources of scattered light. We fully characterize and remove structured and diffuse astrophysical foregrounds including bright stars, the integrated starlight from faint unresolved sources, and diffuse galactic light. For the full data set, we find the surface brightness of the COB to be $\lambda I_{\lambda}^{\mathrm{COB}}$ $=$ 21.98 $\pm$ 1.23 (stat.) $\pm$ 1.36 (cal.) nW m$^{-2}$ sr$^{-1}$. This result supports recent determinations that find a factor of 2 ${-}$ 3 $\times$ more light than expected from the integrated light from galaxies and motivate new diffuse intensity measurements with more capable instruments that can support spectral measurements over the optical and near-IR.

The planets of the TOI-216 system have been previously observed to exhibit large transit timing variations (TTVs), which enabled precise mass characterization of both transiting planets. In the first year of TESS observations, TOI-216 b exhibited grazing transits, precluding a measurement of its radius. In new observations, we demonstrate the orbit of the planet has precessed and it is now fully transiting, so we can accurately measure its radius. TOI-216 b is a puffy Neptune-mass planet, with a much larger radius that is now well constrained to $7.28^{+0.14}_{-0.13}$ $R_{\oplus}$ and a density of $0.253^{+0.017}_{-0.016}$ g cm$^{-3}$. We numerically integrate the system across the TESS observations to update and refine the masses and orbits of both planets, finding the uncertainty in the masses are now dominated by uncertainties in the stellar parameters. TOI-216 b represents a growing class of super-puff planets in orbital resonances and with a companion in a nearly circular orbit, suggesting these planets early evolution is driven by smooth disk migration.

The spin-distribution of accreting neutron stars in low-mass X-ray binary (LMXB) systems shows a concentration of pulsars well below the Keplarian break-up limit. It has been suggested that their spin frequencies may be limited by the emission of gravitational waves, due to the presence of large-scale asymmetries in the internal temperature profile of the star. These temperature asymmetries have been demonstrated to lead to a non-axisymmetric mass-distribution, or `mountain', that generates gravitational waves at twice the spin frequency. The presence of a toroidal magnetic field in the interior of accreting neutron stars has been shown to introduce such anisotropies in the star's thermal conductivity, by restricting the flow of heat orthogonal to the magnetic field and establishing a non-axisymmetric temperature distribution within the star. We revisit this mechanism, extending the computational domain from (only) the crust to the entire star, incorporating more realistic microphysics, and exploring different choices of outer boundary condition. By allowing a magnetic field to permeate the core of the neutron star, we find that the likely level of temperature asymmetry in the inner crust ($\rho \sim 10^{13}$ g cm$^{-3}$) can be up to 3 orders of magnitude greater than the previous estimate, improving prospects for one day detecting continuous gravitational radiation. We also show that temperature asymmetries sufficiently large to be interesting for gravitational wave emission can be generated in strongly accreting neutron stars if crustal magnetic fields can reach $\sim 10^{12}$ G.

We present the Extragalactic Serendipitous Swift Survey (ExSeSS), providing a new well-defined sample constructed from the observations performed using the Swift X-ray Telescope. The ExSeSS sample consists of 79,342 sources detected in the medium (1-2 keV), hard (2-10 keV) or total (0.3-10 keV) energy bands, covering 2086.6 deg$^{2}$ of sky across a flux range of $f_\mathrm{0.3-10keV}\sim10^{-15}-10^{-10}$ erg s$^{-1}$ cm$^{-2}$. Using the new ExSeSS sample we present measurements of the differential number counts of X-ray sources as a function of 2-10 keV flux that trace the population of Active Galactic Nuclei (AGN) in a previously unexplored regime. We find that taking the line-of-sight absorption column density into account has an effect on the differential number count measurements and is vital to obtain agreement with previous results. In the hard band, we obtain a good agreement between the ExSeSS measurements and previous, higher energy data from NuSTAR and Swift/BAT when taking into account the varying column density of the ExSeSS sample as well as the X-ray spectral parameters of each of the samples we are comparing to. We also find discrepancies between the ExSeSS measurements and AGN population synthesis models, indicating a change in the properties of the AGN population over this flux range that is not fully described by current models at these energies, hinting at a larger, moderately obscured population at low redshifts ($z\lesssim0.2$) that the models are not currently taking into account.

Hydrodynamical instabilities, such as the standing accretion-shock instability (SASI), play an essential role in the dynamics of core-collapse supernovae, with observable imprints in the neutrino and gravitational wave signals. Yet, the impact of stellar rotation on the development of SASI is poorly explored. We investigate the conditions favoring the growth of SASI in the presence of rotation through a perturbative analysis. The properties of SASI are compared in two stationary configurations, cylindrical and spherical equatorial, which mainly differ by their advection timescales from the shock to the proto-neutron star surface. Without rotation, the mode $m=1$, corresponding to a one-armed spiral SASI deformation, can be significantly more unstable in the spherical equatorial configuration. In fact, the shorter advection time in the spherical equatorial geometry allows for a larger contribution of the entropic-acoustic coupling from the region of adiabatic compression near the surface of the proto-neutron star. The angular momentum of the collapsing core favors the growth of prograde spiral modes $m=1$ and $m = 2$ in both geometries. Although the growth rate of the spiral instability is systematically faster in spherical geometry, its oscillation frequency is remarkably insensitive to the geometry. Such a contrast with non-rotating flows calls for a deeper understanding of the role of advection in the mechanism of spiral SASI. Our findings suggest that the resonant coupling of acoustic waves with their corotation radius may play a major role in the instability mechanism of collapsing cores with rotation. Elucidating this physical mechanism is essential to interpret the signal from future multi-messenger supernova observations.

In this work, we present a photospheric solar silicon abundance derived using CO5BOLD model atmospheres and the LINFOR3D spectral synthesis code. Previous works have differed in their choice of a spectral line sample and model atmosphere as well as their treatment of observational material, and the solar silicon abundance has undergone a downward revision in recent years. We additionally show the effects of the chosen line sample, broadening due to velocity fields, collisional broadening, model spatial resolution, and magnetic fields. CO5BOLD model atmospheres for the Sun were used in conjunction with the LINFOR3D spectral synthesis code to generate model spectra, which were then fit to observations in the Hamburg solar atlas. We present a sample of 11 carefully selected lines (from an initial choice of 39 lines) in the optical and infrared, made possible with newly determined oscillator strengths for the majority of these lines. Our final sample includes seven optical Si I lines, three infrared Si I lines, and one optical Si II line. We derived a photospheric solar silicon abundance of $\log \epsilon_\mathrm{Si} = 7.57 \pm 0.04$, including a $-0.01$ dex correction from Non-Local Thermodynamic Equilibrium (NLTE) effects. Combining this with meteoritic abundances and previously determined photospheric abundances results in a metal mass fraction Z/X = $0.0220 \pm 0.0020$. We found a tendency of obtaining overly broad synthetic lines. We mitigated the impact of this by devising a de-broadening procedure. The over-broadening of synthetic lines does not substantially affect the abundance determined in the end. It is primarily the line selection that affects the final fitted abundance.

Neutron stars (NSs) and black holes (BHs) are born when the final collapse of the stellar core terminates the lives of stars more massive than about 9 Msun. This can trigger the powerful ejection of a large fraction of the star's material in a core-collapse supernova (CCSN), whose extreme luminosity is energized by the decay of radioactive isotopes such as 56Ni and 56Co. When evolving in close binary systems, the compact relics of such infernal catastrophes spiral towards each other on orbits gradually decaying by gravitational-wave emission. Ultimately, the violent collision of the two components forms a more massive, rapidly spinning remnant, again accompanied by the ejection of considerable amounts of matter. These merger events can be observed by high-energy bursts of gamma rays with afterglows and electromagnetic transients called kilonovae, which radiate the energy released in radioactive decays of freshly assembled rapid neutron-capture elements. By means of their mass ejection and the nuclear and neutrino reactions taking place in the ejecta, both CCSNe and compact object mergers (COMs) are prominent sites of heavy-element nucleosynthesis and play a central role in the cosmic cycle of matter and the chemical enrichment history of galaxies. The nuclear equation of state (EoS) of NS matter, from neutron-rich to proton-dominated conditions and with temperatures ranging from about zero to ~100 MeV, is a crucial ingredient in these astrophysical phenomena. It determines their dynamical processes, their remnant properties even at the level of deciding between NS or BH, and the properties of the associated emission of neutrinos, whose interactions govern the thermodynamic conditions and the neutron-to-proton ratio for nucleosynthesis reactions in the innermost ejecta. This chapter discusses corresponding EoS dependent effects of relevance in CCSNe as well as COMs. (slightly abridged)

Gravitational-wave astronomy allows us to study objects and events invisible to electromagnetic waves. So far, only signals triggered by coalescing binaries have been detected. However, as the interferometers' sensitivities improve over time, we expect to observe weaker signals in the future, e.g. emission of continuous gravitational waves from spinning, isolated neutron stars. Parallax is a well-known method, widely used in electromagnetic astronomical observations, to estimate the distance to a source. In this work, we consider the application of the parallax method to gravitational-wave searches and explore possible distance estimation errors. We show that detection of parallax in the signal from a spinning down source can constrain the neutron star moment of inertia. For instance, we found that the relative error of the moment of inertia estimation is smaller than $10\%$ for all sources closer than 300 pc, for the assumed birth frequency of 700 Hz and for two years of observations by the Einstein Telescope.

Low-mass galaxies are highly susceptible to environmental effects that can efficiently quench star formation. We explore the role of ram pressure in quenching low-mass galaxies ($M_{*}\sim10^{5-9}\,\rm{M}_{\odot}$) within 2 Mpc of Milky Way (MW) hosts using the FIRE-2 simulations. Ram pressure is highly variable across different environments, within individual MW haloes, and for individual low-mass galaxies over time. The impulsiveness of ram pressure -- the maximum ram pressure scaled to the integrated ram pressure prior to quenching -- correlates with whether a galaxy is quiescent or star-forming. The time-scale between maximum ram pressure and quenching is anticorrelated with impulsiveness, such that high impulsiveness corresponds to quenching time-scales $<1$ Gyr. Galaxies in low-mass groups ($M_\mathrm{*,host}\sim10^{7-9}\,\rm{M}_{\odot}$) outside of MW haloes experience typical ram pressure only slightly lower than ram pressure on MW satellites, helping to explain effective quenching via group pre-processing. Ram pressure on MW satellites rises sharply with decreasing distance to the host, and, at a fixed physical distance, more recent pericentre passages are typically associated with higher ram pressure because of greater gas density in the inner host halo at late times. Furthermore, the inner gas density of Local Group-like paired host haloes is larger at small angles ($\lesssim30^\circ$) off the host galaxy disc, compared to isolated hosts. The ram pressure and quiescent fraction of satellites within these low latitude regions are correspondingly elevated, signaling anisotropic quenching via ram pressure around paired hosts.

We study the young star cluster populations in 23 dwarf and irregular galaxies observed by the HST Legacy ExtraGalactic Ultraviolet Survey (LEGUS), and examine relationships between the ensemble properties of the cluster populations and those of their host galaxies: star formation rate (SFR) density ($\Sigma_{\rm SFR}$). A strength of this analysis is the availability of SFRs measured from temporally resolved star formation histories which provide the means to match cluster and host-galaxy properties on several timescales (1-10, 1-100, and 10-100~Myr). Nevertheless, studies of this kind are challenging for dwarf galaxies due to the small numbers of clusters in each system. We mitigate these issues by combining the clusters across different galaxies with similar $\Sigma_{\rm SFR}$ properties. We find good agreement with a well-established relationship ($M_{V}^{brightest}$-SFR), but find no significant correlations between $\Sigma_{\rm SFR}$ and the slopes of the cluster luminosity function, mass function, nor the age distribution. We also find no significant trend between the the fraction of stars in bound clusters at different age ranges ($\Gamma_{1-10}$, $\Gamma_{10-100}$, and $\Gamma_{1-100}$) and $\Sigma_{\rm SFR}$ of the host galaxy. Our data show a decrease in $\Gamma$ over time (from 1-10 to 10-100~Myr) suggesting early cluster dissolution, though the presence of unbound clusters in the youngest time bin makes it difficult to quantify the degree of dissolution. While our data do not exhibit strong correlations between $\Sigma_{\rm SFR}$ and ensemble cluster properties, we cannot rule out that a weak trend might exist given the relatively large uncertainties due to low number statistics and the limited $\Sigma_{\rm SFR}$ range probed.

The adiabatic growth of a central massive black hole could compress the surrounding dark matter halo, leading to a steeper profile of the dark matter halo. This phenomenon is called adiabatic compression. We investigate the adiabatic compression of wave dark matter - a light bosonic dark matter candidate with its mass smaller than a few eV. Using the adiabatic theorem, we show that the adiabatic compression leads to a much denser wave dark matter halo similar to the particle dark matter halo in the semiclassical limit. The compressed wave halo differs from that of the particle halo near the center where the semiclassical approximation breaks down, and the central profile depends on dark matter and the central black hole mass as they determine whether the soliton and low angular momentum modes can survive over the astrophysical time scale without being absorbed by the black hole. Such a compressed profile has several astrophysical implications. As one example, we study the gravitational waves from the inspiral between a central intermediate-mass black hole and a compact solar-mass object within the wave dark matter halo. Due to the enhanced mass density, the compressed wave dark matter halo exerts stronger dynamical friction on the orbiting object, leading to the dephasing of the gravitational waves. The pattern of dephasing is distinctive from that of inspirals in the particle dark matter halo because of the difference in density profile and because of the relatively suppressed dynamical friction force, originating from the wave nature of dark matter. We investigate the prospects of future gravitational wave detectors, such as Laser Interferometer Space Antenna, and identify physics scenarios where the wave dark matter halo can be reconstructed from gravitational wave observations.

Using CANUCS imaging we found an apparent major merger of two $z\sim5$ ultra-low-mass galaxies ($M_\star\sim10^{7}M_\odot$ each) that are doubly imaged and magnified $\sim$12-15$\times$ by the lensing cluster MACS 0417. Both galaxies are experiencing young ($\sim$100 Myr), synchronised bursts of star formation with $\log({\rm sSFR/Gyr^{-1}} )\sim$1.3-1.4, yet SFRs of just $\sim 0.2 M_\odot\ {\rm yr}^{-1}$. They have sub-solar ($Z\sim0.2Z_\odot$) gas-phase metallicities and are connected by an even more metal-poor star-forming bridge. The galaxy that forms from the merger will have a mass of at least $M_\star\sim 2\times10^7 M_\odot$, at least half of it formed during the interaction-induced starburst. More than half of the ionizing photons produced by the system (before and during the merger) will have been produced during the burst. This system provides the first detailed look at a merger involving two high-$z$ ultra-low-mass galaxies of the type believed to be responsible for reionizing the Universe. It suggests that such galaxies can grow via a combination of mass obtained through major mergers, merger-triggered starbursts, and long-term in-situ star formation. If such high-$z$ mergers are common, then merger-triggered starbursts could be significant contributors to the ionizing photon budget of the Universe.

New optical photometric, spectrocopic and imaging polarimetry data are combined with publicly available data to study some of the physical properties of the two H-poor superluminous supernovae (SLSN) SN 2021bnw and SN 2021fpl. For each SLSN, the best-fit parameters obtained from the magnetar model with \texttt{MOSFiT} do not depart from the range of parameter obtained on other SLSNe discussed in the literature. A spectral analysis with \texttt{SYN++} shows that SN 2021bnw is a W Type, Fast evolver, while SN 2021fpl is a 15bn Type, Slow evolver. The analysis of the polarimetry data obtained on SN 2021fpl at four epochs (+1.8, +20.6, +34.1 and +43.0 days, rest-frame) shows $> 3\sigma$ polarization detections in the range 0.8--1 $\%$. A comparison of the spectroscopy data suggests that SN 2021fpl underwent a spectral transition a bit earlier than SN 2015bn, during which, similarly, it could have underwent a polarization transition. The analysis of the polarimetry data obtained on SN 2021bnw do not show any departure from symmetry of the photosphere at an empirical diffusion timescale of $\approx$ 2 (+81.1 days rest-frame). This result is consistent with those on the sample of W Type SLSN observed at empirical diffusion timescale $\le$ 1 with that technique, even though it is not clear the effect of limited spectral windows varying from one object to the other. Measurements at higher empirical diffusion timescale may be needed to see any departure from symmetry as it is discussed in the literature for SN 2017egm.

HCN, HNC, and their isotopologues are ubiquitous molecules that can serve as chemical thermometers and evolutionary tracers to characterize star-forming regions. Despite their importance in carrying information that is vital to studies of the chemistry and evolution of star-forming regions, the collision rates of some of these molecules have not been available for rigorous studies in the past. We perform an up-to-date gas and dust chemical characterization of two different star-forming regions, TMC 1-C and NGC 1333-C7, using new collisional rates of HCN, HNC, and their isotopologues. We investigated the possible effects of the environment and stellar feedback in their chemistry and their evolution. With millimeter observations, we derived their column densities, the C and N isotopic fractions, the isomeric ratios, and the deuterium fractionation. The continuum data at 3 mm and 850 $\mu$m allowed us to compute the emissivity spectral index and look for grain growth as an evolutionary tracer. The H$^{13}$CN/HN$^{13}$C ratio is anticorrelated with the deuterium fraction of HCN, thus it can readily serve as a proxy for the temperature. The spectral index $(\beta\sim 1.34-2.09)$ shows a tentative anticorrelation with the H$^{13}$CN/HN$^{13}$C ratio, suggesting grain growth in the evolved, hotter, and less deuterated sources. Unlike TMC 1-C, the south-to-north gradient in dust temperature and spectral index observed in NGC 1333-C7 suggests feedback from the main NGC 1333 cloud. With this up-to-date characterization of two star-forming regions, we found that the chemistry and the physical properties are tightly related. The dust temperature, deuterium fraction, and the spectral index are complementary evolutionary tracers. The large-scale environmental factors may dominate the chemistry and evolution in clustered star-forming regions.

The 21 cm intensity mapping (IM) technique provides us with an efficient way to observe the cosmic large-scale structure (LSS). From the LSS data, one can use the baryon acoustic oscillation and redshift space distortion to trace the expansion and growth history of the universe, and thus measure the dark energy parameters. In this paper, we make a forecast for cosmological parameter estimation with the synergy of three 21 cm IM experiments. Specifically, we adopt a novel joint survey strategy, FAST ($0<z<0.35$)+SKA1-MID ($0.35<z<0.8$)+HIRAX ($0.8<z<2.5$), to measure dark energy. We simulate the 21 cm IM observations under the assumption of excellent foreground removal. We find that the synergy of three experiments could place quite tight constraints on cosmological parameters. For example, it provides $\sigma(\Omega_{\rm m})=0.0039$ and $\sigma(H_0)=0.27\ \rm km\ s^{-1}\ Mpc^{-1}$ in the $\Lambda$CDM model. Notably, the synergy could break the cosmological parameter degeneracies when constraining the dynamical dark energy models. Concretely, the joint observation offers $\sigma(w)=0.019$ in the $w$CDM model, and $\sigma(w_0)=0.085$ and $\sigma(w_a)=0.32$ in the $w_0w_a$CDM model. These results are better than or equal to those given by CMB+BAO+SN. In addition, when the foreground removal efficiency is relatively low, the strategy still performs well. Therefore, the 21 cm IM joint survey strategy is promising and worth pursuing.

Horava gravity has been proposed as a renormalizable, higher-derivative, Lorentz-violating quantum gravity model without ghost problems. A Horava gravity based dark energy (HDE) model for dynamical dark energy has been also proposed earlier by identifying all the extra (gravitational) contributions from the Lorentz-violating terms as an effective energy-momentum tensor in Einstein equation. We consider a complete CMB, BAO, and SNe Ia data test of the HDE model by considering general perturbations over the background perfect HDE fluid. Except from BAO, we obtain the preference of non-flat universes for all other data-set combinations. We obtain a positive result on the cosmic tensions between the Hubble constant H0 and the cosmic shear S8, because we have a shift of H0 towards a higher value, though not enough for resolving the H0 tension, but the value of S8 is unaltered. This is in contrast to a rather decreasing H0 but increasing S8 in a non-flat LCDM. For all other parameters, like Omega_m and Omega_Lambda, we obtain quite comparable results with those of LCDM for all data sets, especially with BAO, so that our results are close to a cosmic concordance between the datasets, contrary to the standard non-flat LCDM. We also obtain some undesirable features but we propose several promising ways for improvements by generalizing our analysis.

We present a new technique to perform passive bistatic subsurface radar probes on airless planetary bodies. This technique uses the naturally occurring radio impulses generated when high-energy cosmic rays impact the body's surface. As in traditional radar sounding, the downward-beamed radio emission from each individual cosmic ray impact will reflect off subsurface dielectric contrasts and propagate back up to the surface to be detected. We refer to this technique as Askaryan radar after the fundamental physics process, the Askaryan effect, that produces this radio emission. This technique can be performed from an orbiting satellite, or from a surface lander, but since the radio emission is generated beneath the surface, an Askaryan radar can completely bypass the effects of surface clutter and backscatter typically associated with surface-penetrating radar. We present the background theory of Askaryan subsurface radar and show results from both finite-difference time-domain (FDTD) and Monte Carlo simulations that confirm that this technique is a promising planetary radar sounding method, producing detectable signals for realistic planetary science applications.

Type II solar radio bursts are caused by magnetohydrodynamics (MHD) shocks driven by solar eruptive events such as Coronal Mass Ejections (CMEs). Often both fundamental and harmonic bands of type II bursts are split into sub-bands, generally believed to be coming from upstream and downstream regions of the shock; however this explanation remains unconfirmed. Here we present combined results from imaging analysis of type II radio burst band-splitting and other fine structures, observed by the Murchison Widefield Array (MWA) and extreme ultraviolet observations from Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA) on 2014-Sep-28. The MWA provides imaging-spectroscopy in the range of 80-300 MHz with a time resolution of 0.5 s and frequency resolution of 40 kHz. Our analysis shows that the burst was caused by a piston-driven shock with a driver speed of $\sim$112 km s$^{-1}$ and shock speed of $\sim$580 km s$^{-1}$. We provide rare evidence that band-splitting is caused by emission from multiple parts of the shock (as opposed to the upstream/downstream hypothesis). We also examine the small-scale motion of type II fine structure radio sources in MWA images. We suggest that this small-scale motion may arise due to radio propagation effects from coronal turbulence, and not because of the physical motion of the shock location. We present a novel technique that uses imaging spectroscopy to directly determine the effective length scale of turbulent density perturbations, which is found to be 1 - 2 Mm. The study of the systematic and small-scale motion of fine structures may therefore provide a measure of turbulence in different regions of the shock and corona.

In multiple stellar systems interactions among the companion stars and their discs affect planet formation. In the circumstellar case tidal truncation makes protoplanetary discs smaller, fainter and less long-lived than those evolving in isolation, thereby reducing the amount of material (gas and dust) available to assemble planetary embryos. On the contrary, in the circumbinary case the reduced accretion can increase the disc lifetime, with beneficial effects on planet formation. In this chapter we review the main observational results on discs in multiple stellar systems and discuss their possible explanations, focusing on recent numerical simulations, mainly dealing with dust dynamics and disc evolution. Finally, some open issues and future research directions are examined.

Although many exoplanets have been indirectly detected over the last years, direct imaging of them with ground-based telescopes remains challenging. In the presence of atmospheric fluctuations, it is ambitious to resolve the high brightness contrasts at the small angular separation between the star and its potential partners. Post-processing of telescope images has become an essential tool to improve the resolvable contrast ratios. This paper contributes a post-processing algorithm for fast-cadence imaging, which deconvolves sequences of telescopes images. The algorithm infers a Bayesian estimate of the astronomical object as well as the atmospheric optical path length, including its spatial and temporal structures. For this, we utilize physics-inspired models for the object, the atmosphere, and the telescope. The algorithm is computationally expensive but allows to resolve high contrast ratios despite short observation times and no field rotation. We test the performance of the algorithm with point-like companions synthetically injected into a real data set acquired with the SHARK-VIS pathfinder instrument at the LBT telescope. Sources with brightness ratios down to $6\cdot10^{-4}$ to the star are detected at $185$ mas separation with a short observation time of $0.6\,\text{s}$.

The Cosmological Principle, that the Universe is homogeneous and isotropic on sufficiently large scales, underpins the standard model of cosmology. However, a recent analysis of 1.36 million infrared-selected quasars has identified a significant tension in the amplitude of the number-count dipole compared to that derived from the CMB, thus challenging the Cosmological Principle. Here we present a Bayesian analysis of the same quasar sample, testing various hypotheses using the Bayesian evidence. We find unambiguous evidence for the presence of a dipole in the distribution of quasars with a direction that is consistent with the dipole identified in the CMB. However, the amplitude of the dipole is found to be 2.7 times larger than that expected from the conventional kinematic explanation of the CMB dipole, with a statistical significance of $5.7\sigma$. To compare these results with theoretical expectations, we sharpen the $\Lambda$CDM predictions for the probability distribution of the amplitude, taking into account a number of observational and theoretical systematics. In particular, we show that the presence of the galactic plane mask causes a considerable loss of dipole signal due to a leakage of power into higher multipoles, exacerbating the discrepancy in the amplitude. By contrast, we estimate using probabilistic arguments that the source evolution of quasars improves the discrepancy, but only mildly so. These results support the original findings of an anomalously large quasar dipole, independent of the statistical methodology used.

Early afterglow observations of gamma-ray bursts (GRBs) are valuable for exploring the properties of their jets and ambient medium. We report our photometric and spectroscopic observations of GRB 210104A and discuss its jet properties with multiwavelength data. Our spectroscopic observation reveals several absorption features and a tentative redshift of 0.46 is identified. A bright optical flare that has a peak brightness of $R=13$ mag at $112\pm 7$~s was observed in the $R$ band during $67\sim 165$ seconds post the GRB trigger. The flux of the $R$-band afterglow decays with a slope of $\alpha_{\rm O}={-0.91\pm 0.03}$ at $t>650$~s. The early X-ray afterglow lightcurve is a smooth bump, and it decays with a slope of $\alpha_{\rm X}=-1.18\pm 0.01$ at late epoch. Our joint spectral fit to the optical-X-ray afterglows during $(1.1-1.3)\times 10^{4}$~s yields a photon index $\Gamma_{\rm O,X}=-1.82\pm 0.04$. The derived host galaxy extinction is $A_{R}=0.87$. Attributing the early optical flare to the reverse-shock (RS) emission and the late optical-X-ray emission to the forward shock emission; the optical and X-ray lightcurves at $t<3\times 10^4$~s can be well fit adopting an Markov Chain Monte Carlo algorithm. Comparing the properties of GRB 210104A with other GRBs that have detection of bright RS emission, we show that its jet is mildly magnetized ($R_{\rm B}=28$), with high radiation efficiency ($77\%$), is sub-energetic ($E_{\rm k, iso}=4.5\times 10^{51}$ erg), and moderately relativistic ($\Gamma_0\sim 35$) in a density medium ($n_{0}\sim 417\;{\rm cm}^{-3}$). It follows the tight $L_{\gamma,\rm iso}-E_{\rm p,z}-\Gamma_{0}$ relation as with typical GRBs.

We carry out sensitive searches for the CO J=1-0 and J=2-1 lines in the giant extragalactic HI ring in Leo to investigate the star formation process within environments where gas metallicities are close to solar but physical conditions are different than those typical of bright galaxy disks. Our aim is to check the range of validity of known scaling relations. We use the IRAM-30m telescope to observe eleven regions close to HI gas peaks or where sparse young massive stars have been found. For all pointed observations we reached a spectral noise between 1 and 5~mK for at least one observed frequencies at 2~km/s spectral resolution. We marginally detect two CO J=1-0 lines in the star forming region Clump~1 of the Leo ring, whose radial velocities are consistent with those of Halpha lines but line widths are much smaller than observed for virialized molecular clouds of similar mass in galaxies. The low signal-to-noise ratio, the small line widths and the extremely low number densities suggest that a more standard population of molecular clouds, still undetected, might be in place. Using upper limits to the CO lines, the most sensitive pointed observations show that the molecular gas mass surface density is lower than expected from the extrapolation of the molecular Kennicutt-Schmidt relation established in the disk of galaxies. The sparse stellar population in the ring, possibly forming ultra diffuse dwarf galaxies, might then be the result of a short molecular gas depletion time in this extreme environment.}

The cosmic web consists of a complex configuration of voids, walls, filaments, and clusters, which formed under the gravitational collapse of Gaussian fluctuations. Understanding under what conditions these different structures emerge from simple initial conditions, and how different cosmological models influence their evolution, is central to the study of the large-scale structure. Here, we present a general formalism for setting up initial random density and velocity fields satisfying non-linear constraints for specialized N-body simulations. These allow us to link the non-linear conditions on the eigenvalue and eigenvector fields of the deformation tensor, as specified by caustic skeleton theory, to the current-day cosmic web. By extending constrained Gaussian random field theory, and the corresponding Hoffman-Ribak algorithm, to non-linear constraints, we probe the statistical properties of the progenitors of the walls, filaments, and clusters of the cosmic web. Applied to cosmological N-body simulations, the proposed techniques pave the way towards a systematic investigation of the evolution of the progenitors of the present-day walls, filaments, and clusters, and the embedded galaxies, putting flesh on the bones of the caustic skeleton. The developed nonlinear constrained random field theory is valid for generic cosmological conditions. For ease of visualization, the case study presented here probes the two-dimensional caustic skeleton.

In 2019 while launching a multidisciplinary research project aimed at developing the Puna de Atacama region as a natural laboratory, investigators within the University of Atacama (Chile) conducted a bibliographic search identifying previously studied geographical points of the region and of potential interest for planetary science and astrobiology research. This preliminary work highlighted a significant absence in foreign publications consideration of local institutional involvement. In light of this, a follow-up study was carried out to confirm or refute these first impressions, by comparing the search in two bibliographic databases: Web of Science and Scopus. The results show that almost 60% of the publications based directly on data from the Puna, the Altiplano or the Atacama Desert with objectives related to planetary science or astrobiology do not include any local institutional partner (Argentina, Bolivia, Chile and Peru). Indeed, and beyond the ethical questioning of international collaborations, Latin-American planetary science deserve a strategic structuring, networking, as well as a road map at a national and continental scale, not only to enhance research, development and innovation but also to protect an exceptional natural heritage sampling extreme environmental niches on Earth. Examples of successful international collaborations such as the field of meteorites, terrestrial analogues and space exploration in Chile or astrobiology in Mexico are given as illustrations and possible directions to follow in order to develop planetary sciences in South America.

We present the measurement of the energy dependence of the boron flux in cosmic rays and its ratio to the carbon flux \textcolor{black}{in an energy interval from 8.4 GeV$/n$ to 3.8 TeV$/n$} based on the data collected by the CALorimetric Electron Telescope (CALET) during $\sim 6.4$ years of operation on the International Space Station. An update of the energy spectrum of carbon is also presented with an increase in statistics over our previous measurement. The observed boron flux shows a spectral hardening at the same transition energy $E_0 \sim 200$ GeV$/n$ of the C spectrum, though B and C fluxes have different energy dependences. The spectral index of the B spectrum is found to be $\gamma = -3.047\pm0.024$ in the interval $25 < E < 200$ GeV$/n$. The B spectrum hardens by $\Delta \gamma_B=0.25\pm0.12$, while the best fit value for the spectral variation of C is $\Delta \gamma_C=0.19\pm0.03$. The B/C flux ratio is compatible with a hardening of $0.09\pm0.05$, though a single power-law energy dependence cannot be ruled out given the current statistical uncertainties. A break in the B/C ratio energy dependence would support the recent AMS-02 observations that secondary cosmic rays exhibit a stronger hardening than primary ones. We also perform a fit to the B/C ratio with a leaky-box model of the cosmic-ray propagation in the Galaxy in order to probe a possible residual value $\lambda_0$ of the mean escape path length $\lambda$ at high energy. We find that our B/C data are compatible with a non-zero value of $\lambda_0$, which can be interpreted as the column density of matter that cosmic rays cross within the acceleration region.

Detection of delayed sub-TeV photons from Gamma-Ray Bursts (GRBs) by MAGIC and HESS has proved the promising future of GRB afterglow studies with the Cherenkov Telescope Array (CTA), the next-generation gamma-ray observatory. With the unprecedented sensitivity of CTA, afterglow detection rates are expected to increase dramatically. In this paper, we explore the multi-dimensional afterglow parameter space to see the detectability of sub-TeV photons by CTA. We use a one-zone electron synchrotron and synchrotron self-Compton model to obtain the spectral energy distribution. We consider bursts going off in a medium of homogenous density. The blast wave is assumed to be radiatively inefficient and evolving adiabatically. Considering that the electron acceleration is not efficient if the acceleration timescale exceeds the radiative cooling timescale, we find that the Sub-TeV emission is always due to the self-Compton process. We find that jets with high kinetic energy or large bulk Lorentz factor decelerating into a dense ambient medium offer better detection prospects for CTA. For relatively lower values of the downstream magnetic field, electrons are slow cooling and the emitted radiation is positively correlated with magnetic field. For larger magnetic fields, electron population enters the fast cooling phase where the radiated flux is inversely proportional to magnetic field. We apply our results in the context of bright TeV afterglows detected in recent years. Our results indicate that cosmological short GRBs have only moderate prospects of detection by CTA while local Neutron Star merger counterparts can be detected if the jet is launched towards the observer.

Analyzing the clustering of galaxies at the field level in principle promises access to all the cosmological information available. Given this incentive, in this paper we investigate the performance of field-based forward modeling approach to galaxy clustering using the effective field theory (EFT) framework of large-scale structure (LSS). We do so by applying this formalism to a set of consistency and convergence tests on synthetic datasets. We explore the high-dimensional joint posterior of LSS initial conditions by combining Hamiltonian Monte Carlo sampling for the field of initial conditions, and slice sampling for cosmology and model parameters. We adopt the Lagrangian perturbation theory forward model from [1], up to second order, for the forward model of biased tracers. We specifically include model mis-specifications in our synthetic datasets within the EFT framework. We achieve this by generating synthetic data at a higher cutoff scale $\Lambda_0$, which controls which Fourier modes enter the EFT likelihood evaluation, than the cutoff $\Lambda$ used in the inference. In the presence of model mis-specifications, we find that the EFT framework still allows for robust, unbiased joint inference of a) cosmological parameters - specifically, the scaling amplitude of the initial conditions - b) the initial conditions themselves, and c) the bias and noise parameters. In addition, we show that in the purely linear case, where the posterior is analytically tractable, our samplers fully explore the posterior surface. We also demonstrate convergence in the cases of nonlinear forward models. Our findings serve as a confirmation of the EFT field-based forward model framework developed in [2-7], and as another step towards field-level cosmological analyses of real galaxy surveys.

Active galactic nuclei (AGN) are supermassive black holes with luminous accretion disks found in some galaxies, and are thought to play an important role in galaxy evolution. However, traditional optical spectroscopy for identifying AGN requires time-intensive observations. We train a convolutional neural network (CNN) to distinguish AGN host galaxies from non-active galaxies using a sample of 210,000 Sloan Digital Sky Survey galaxies. We evaluate the CNN on 33,000 galaxies that are spectrally classified as composites, and find correlations between galaxy appearances and their CNN classifications, which hint at evolutionary processes that affect both galaxy morphology and AGN activity. With the advent of the Vera C. Rubin Observatory, Nancy Grace Roman Space Telescope, and other wide-field imaging telescopes, deep learning methods will be instrumental for quickly and reliably shortlisting AGN samples for future analyses.

Massively parallel multi-object spectrographs are on the leading edge of cosmology instrumentation. The highly successful Dark Energy Spectroscopic Instrument (DESI) which begun survey operations in May 2021, for example, has 5,000 robotically-actuated multimode fibers, which deliver light from thousands of individual galaxies and quasars simultaneously to an array of high-resolution spectrographs off-telescope. The redshifts are individually measured, thus providing 3D maps of the Universe in unprecedented detail, and enabling precise measurement of dark energy expansion and other key cosmological parameters. Here we present new work in the design and prototyping of the next generation of fiber-positioning robots. At 6.2 mm center-to-center pitch, with 1-2 um positioning precision, and in a scalable form factor, these devices will enable the next generation of cosmology instruments, scaling up to instruments with 10,000 to 25,000 fiber robots.

We investigate the statistics of the available Pantheon+ dataset. Noticing that the $\chi^2$ value for the best-fit $\Lambda$CDM model to the real data is small, we quantify how significant its smallness is by calculating the distribution of $\chi^2$ values for the best-fit $\Lambda$CDM model fit to mock Pantheon+-like datasets, using the provided covariance matrix. We further investigate the distribution of the residuals of the Pantheon+ dataset, with respect to the best-fit $\Lambda$CDM model, and notice they scatter smaller than would be expected from the covariance matrix but find no significant amount of kurtosis. These results point to the conclusion that the Pantheon+ covariance matrix is over-estimated. One simple interpretation of these results is a $\sim$5\% overestimation of errors on SN distances in Pantheon+ data. When the covariance matrix is reduced by subtracting an intrinsic scatter term from the diagonal terms of the covariance matrix, the best-fit $\chi^2$ for the $\Lambda$CDM model achieves a normal value of 1580 and no deviation from $\Lambda$CDM is detected. We further quantify how consistent the $\Lambda$CDM model is with respect to the modified data with the subtracted covariance matrix using model independent reconstruction techniques such as the iterative smoothing method and we find that the standard model is consistent with the data.

It has been conjectured that a specific set of Cauchy data may lead to an unexpected growth of the gauge fields in conformally flat cosmological backgrounds. After introducing a class of models where the (hyper)magnetic fields are produced quantum mechanically from the early variation of the gauge coupling, we show that the standard evolution for wavelengths larger than the Hubble radius fully accounts for the anomalous scaling which is however unable to increase the magnetic power spectra. This dynamical evolution has nothing to do with a peculiar choice of the initial conditions but rather with the rate of variation of the gauge coupling during inflation. While the correct scaling follows from the properly normalized solutions of the mode functions and from the late dominance of the conductivity, a further amplification of the magnetic power spectra can only occur when the post-inflationary expansion rate is slower than radiation and it is again unrelated to the initial data possibly set during the inflationary stage.

Building upon existing signal processing techniques and open-source software, this paper presents a baseline design for an RF System-on-Chip Frequency Division Multiplexed readout for a spatio-spectral focal plane instrument based on low temperature detectors. A trade-off analysis of different FPGA carrier boards is presented in an attempt to find an optimum next-generation solution for reading out larger arrays of Microwave Kinetic Inductance Detectors (MKIDs). The ZCU111 RF SoC FPGA board from Xilinx was selected, and it is shown how this integrated system promises to increase the number of pixels that can be read out (per board) which enables a reduction in the readout cost per pixel, the mass and volume, and power consumption, all of which are important in making MKID instruments more feasible for both ground-based and space-based astrophysics. The on-chip logic capacity is shown to form a primary constraint on the number of MKIDs which can be read, channelised, and processed with this new system. As such, novel signal processing techniques are analysed, including Digitally Down Converted (DDC)-corrected sub-maximally decimated sampling, in an effort to reduce logic requirements without compromising signal to noise ratio. It is also shown how combining the ZCU111 board with a secondary FPGA board will allow all 8 ADCs and 8 DACs to be utilised, providing enough bandwidth to read up to 8,000 MKIDs per board-set, an eight-fold improvement over the state-of-the-art, and important in pursuing 100,000 pixel arrays. Finally, the feasibility of extending the operational frequency range of MKIDs to the 5 - 10 GHz regime (or possibly beyond) is investigated, and some benefits and consequences of doing so are presented.

We present the STScI WFC3 project webpage, Transiting Exoplanets List of Space Telescope Spectroscopy, TrExoLiSTS. It tabulates existing observations of transiting exoplanet atmospheres, available in the MAST archive made with HST WFC3 using the stare or spatial scan mode. A parallel page is available for all instruments aboard JWST using the spectral Time Series Observation (TSO) mode. The web pages include observations obtained during primary transits, secondary eclipses and phase curves. TREXOLISTS facilitates proposal preparation for programs that are highly-complementary to existing programs in terms of targets, wavelength coverage, as well as reduces duplication and redundant effort. Reference for the quality of the HST WFC3 visits taken more than 1.5 years ago are made available via including diagrams of the direct image, white light curve and drift of the spectral time series across the detector. Future improvements to the webpage will include: Expanding program query to other HST instruments and reference for the quality of JWST visits.

The accretion of matter onto primordial black holes (PBHs) during the dark ages and the subsequent energy injection in the medium should have left imprints on the cosmic microwave background (CMB) anisotropies. Recent works have claimed stringent CMB limits on the PBH abundance, hardly compatible with a PBH interpretation of the gravitational-wave observations of binary BH mergers. By using a more realistic accretion model based on hydrodynamical simulations and conservative assumptions for the emission efficiency, we show that CMB limits on the PBH abundance are up to two orders of magnitude less stringent than previously estimated between $10$ and $10^4$ M$_\odot$. This reopens the possibility that PBHs might explain at the same time (at least a fraction of) the dark matter, some of the LIGO-Virgo-KAGRA binary BH mergers and the existence of super-massive BHs. More generally, we emphasize that PBH accretion can be a rather complex physical process with velocity dependences that are hard to assess, which introduces large uncertainties in accretion-based limits on the PBH abundance.

Gamma-Ray Bursts (GRBs) are panchromatic, highly energetic transients whose energy emission mechanism is still debated. One of the possible explanations is the standard fireball model, which can be tested with the closure relations (CRs), or relations between the temporal and spectral indices of a GRB. To test these, we compile an extensive sample of radio afterglow light curves (LCs) that span from 1997 to 2020, the most comprehensive analysis of GRBs with radio observations to date. We fit 202 LCs from 82 distinct GRBs with a broken power law, obtaining a sample of 26 that display a clear break and a subsample of 14 GRBs that present a radio plateau. We test these samples against CRs corresponding to a constant-density interstellar medium (ISM) or a stellar wind medium in both fast- and slow-cooling regimes, as well as three additional density profiles, $k = 1, 1.5, 2.5,$ following $n \propto r^{-k}$ , and consider sets of CRs both with and without energy injection. We find that 12 of the 26 GRBs (46%), of which 7/12 present a radio plateau, fulfill at least one CR in the sets tested, suggesting our data is largely incompatible with the standard fireball model. Of the fulfilled CRs, the most preferred environment is the ISM, SC, $\nu_m < \nu < \nu_c$ without energy injection. Our results are consistent with previous studies that test the standard fireball model via the CRs in radio.

In the past years, we have undertaken an extensive investigation of LMC and SMC star clusters based on HST data. We present photometry and astrometry of stars in 101 fields observed with the WFC/ACS, UVIS/WFC3 and NIR/WFC3 cameras. These fields comprise 113 star clusters. We provide differential-reddening maps and illustrate various scientific outcomes that arise from the early inspection of the photometric catalogs. In particular, we provide new insights on the extended main-sequence turn-off (eMSTO) phenomenon: i) We detected eMSTOs in two clusters, KMHK361 and NGC265, which had no previous evidence of multiple populations. This finding corroborates the conclusion that the eMSTO is a widespread phenomenon among clusters younger than ~2 Gyr. ii) The homogeneous color-magnitude diagrams (CMDs) of 19 LMC clusters reveal that the distribution of stars along the eMSTO depends on cluster age. iii) We discovered a new feature along the eMSTO of NGC1783, which consists of a distinct group of stars going on the red side of the eMSTO in CMDs composed of ultraviolet filters. Furthermore, we derived the proper motions of stars in the fields of view of clusters with multi-epoch images. Proper motions allowed us to separate the bulk of bright field stars from cluster members and investigate the internal kinematics of stellar populations in various LMC and SMC fields. As an example, we analyze the field around NGC346 to disentangle the motions of its stellar populations, including NGC364 and BS90, young and pre-MS stars in the star-forming region associated with NGC346, and young and old field stellar populations of the SMC. Based on these results and the fields around five additional clusters, we find that young SMC stars exhibit elongated proper-motion distributions that point toward the LMC, thus bringing new evidence for a kinematic connection between the LMC and SMC.

Joint analysis of CMB and large-scale structure at high redshifts provide new and unique windows into unexplored epochs of early structure formation. Here, we demonstrate how cosmic infrared background and high-redshift galaxies can be jointly analysed with CMB to probe the epoch of helium reionization ($2<z<4$) on the light cone using kinetic Sunyaev Zel'dovich tomography. Characterising this epoch has great potential significance for understanding astrophysics of galaxy formation, quasar activity and formation of the super-massive black holes. We find a detection at $8-10\sigma$ can be expected from combinations of data from CCAT-prime, Vera Rubin Observatory and CMB-S4 in the upcoming years.

We study a possibility of constraining the isotropic cosmic birefringence with help of cosmic microwave background polarisation data without relying on any assumptions about the Galactic foreground angular power spectra and in particular on their EB correlation. We propose a new analysis framework aiming at measuring isotropic cosmic birefringence, based on a generalised parametric component separation approach, which accounts simultaneously on the presence of galactic foregrounds, relevant instrumental effects and external priors. We find that in the context of an upcoming multi-frequency CMB instrument, assuming in-lab calibration priors we are able to constrain instrumental polarisation angle for each frequency band and correct the observed sky component maps accordingly. We then produce an instrumental-effect-corrected and foreground-cleaned CMB map, which we use to estimate the isotropic birefringence angle and the tensor-to-scalar ratio, accounting on statistical and systematic uncertainties incurred during the entire procedure. In particular, in the case of a Simons Observatory-like, three Small Aperture Telescopes, we derive a uncertainty on the birefringence angle of $\sigma(\beta_{b}) = 0.07^\circ$, assuming calibration priors for all frequency channels with the precision of $\sigma_{\alpha_i}= 0.1^\circ$ and the standard cosmology. This implies that, using our method and given the calibration precision expected for current, the near future ground-based multi-frequency experiments could confirm or disprove the recently detected value of $\beta_b=0.35^\circ$ with a significance of up to $5 \sigma$. [abridged version]

(Abridged) Studies have shown that a Jovian mass planet embedded in a viscous protoplanetary disc (PPD) can accrete gas efficiently through the gap and doubles its mass in $\sim 0.1$ Myr. The planet also migrates inwards on a timescale of $\sim 0.1$ Myr. These timescales are short compared to PPD lifetimes, and raise questions about the origins of cold giant exoplanets. However, PPDs are unlikely to be globally turbulent, and instead they may launch magnetised winds such that accretion towards the star occurs in laminar accretion flows located in narrow layers near the surfaces of the disc. The aim of this study is to examine the rate at which gas accretes onto Jovian mass planets that are embedded in layered PPDs. We use 3D hydrodynamical simulations of planets embedded in PPDs, in which a constant radial mass flux towards the star of ${\dot m} = 10^{-8}$ M$_{\odot}$ yr$^{-1}$ is sustained. We consider a classical viscous alpha model, and also models in which an external torque is applied in narrow surface layers to mimic the effects of a magnetised wind. The accreting layers are parameterised by their column densities $\Sigma_{\rm A}$, and we consider values in the range 0.1 to 10 g cm$^{-2}$. The viscous model gives results in agreement with previous studies. We find the accretion rate onto the planet in the layered models crucially depends on the planet's ability to block the wind-induced mass flow. For $\Sigma_{\rm A}=10$ g cm$^{-2}$, the planet torque can block the mass flow through the disc, accretion onto the planet is slow, and a mass doubling time of 10 Myr is obtained. For $\Sigma_{\rm A}=0.1$ g cm$^{-2}$, accretion is fast and the mass doubling time is 0.2 Myr. Although the radial mass flow through the layered disc models is always $10^{-8}$ M$_{\odot}$ yr$^{-1}$, adopting different values of $\Sigma_{\rm A}$ leads to very different gas accretion rates onto gas giant planets.

We present a spectro-photometric analysis of 2880 cool white dwarfs within 100 pc of the Sun and cooler than Teff = 10,000 K, with grizy Pan-STARRS photometry and Gaia trigonometric parallaxes available. We also supplement our data sets with near-infrared JHK photometry, when available, which is shown to be essential for interpreting the coolest white dwarfs in our sample. We perform a detailed analysis of each individual object using state-of-the-art model atmospheres appropriate for each spectral type including DA, DC, DQ, DZ, He-rich DA, and the so-called IR-faint white dwarfs. We discuss the temperature and mass distributions of each subsample, as well as revisit the spectral evolution of cool white dwarfs. We find little evidence in our sample for the transformation of a significant fraction of DA stars into He-atmosphere white dwarfs through the process of convective mixing between Teff = 10,000 K and 6500 K, although the situation changes drastically in the range Teff = 6500 - 5500 K where the fraction of He-atmosphere white dwarfs reaches 45%. However, we also provide strong evidence that at even cooler temperatures (Teff < 5200 K), most DC white dwarfs have H atmospheres. We discuss a possible mechanism to account for this sudden transformation from He- to H-atmosphere white dwarfs involving the onset of crystallization and the occurrence of magnetism. Finally, we also argue that DQ, DZ, and DC white dwarfs may form a more homogeneous population than previously believed.

We present fast (~200s sampling) ugriz photometry of the low mass AGN NGC 4395 with the Liverpool Telescope, followed by very fast (3s sampling) us, gs, rs, is and zs simultaneous monitoring with HiPERCAM on the 10.4m GTC. These observations provide the fastest ever AGN multiband photometry and very precise lag measurements. Unlike in all other AGN, gs lags us by a large amount, consistent with disc reprocessing but not with reprocessing in the Broad Line Region (BLR). There is very little increase in lag with wavelength at long wavelengths, indicating an outer edge (Rout) to the reprocessor. We have compared truncated disc reprocessing models to the combined HiPERCAM and previous X-ray/UV lags. For the normally accepted mass of 3.6E5 solar, we obtain reasonable agreement with zero spin, Rout ~1700 Rg, and the DONE physically-motivated temperature-dependent disc colour correction factor (fcol). A smaller mass of 4E4 solar can only be accomodated if fcol=2.4, which is probably unrealistically high. Disc self gravity is probably unimportant in this low mass AGN but an obscuring wind may provide an edge. For the small mass the dust sublimation radius is similar to Rout, so the wind could be dusty. However for the more likely large mass the sublimation radius is further out so the optically-thick base of a line-driven gaseous wind is more likely. The inner edge of the BLR is close to Rout in both cases. These observations provide the first good evidence for a truncated AGN disc and caution that truncation should be included in reverberation lag modelling.

The Red-Giant Branch Bump (RGBB) is one of the most noteworthy features in the red-giant luminosity function of stellar clusters. It is caused by the passage of the hydrogen-burning shell through the composition discontinuity left at the point of the deepest penetration by the convective envelope. When crossing the discontinuity the usual trend in increasing luminosity reverses for a short time before it increases again, causing a zig- zag in the evolutionary track. In spite of its apparent simplicity the actual physical reason behind the decrease in luminosity is not well understood and several different explanations have been offered. Here we use a recently proposed simple toy-model for the structure of low-mass red giants, together with previous results, to show beyond reasonable doubt that the change in luminosity at the RGBB can be traced to the change in the mean molecular weight of the layers on top of the burning shell. And that these changes happen on a nuclear timescale. The change in the effective mean molecular weight, as the burning shell approaches the discontinuity, causes a drop in the temperature of the burning shell which is attenuated by the consequent feedback contraction of the layers immediately below the burning shell. Our work shows that, when applied correctly, including the feedback on the structure of the core together with of the increase in the mass of the core, shell-source homology relations do a great quantitative job in explaining the properties of full evolutionary models at the RGBB.

We present JWST-MIRI MRS spectra of the protoplanetary disk around the low-mass T Tauri star GW Lup from the MIRI mid-INfrared Disk Survey (MINDS) GTO program. Emission from $^{12}$CO$_{2}$, $^{13}$CO$_{2}$, H$_{2}$O, HCN, and C$_{2}$H$_{2}$ is identified in this disk with $^{13}$CO$_{2}$ being detected for the first time in a protoplanetary disk. We characterize the chemical and physical conditions in the inner few au of the GW Lup disk using these molecules as probes. The spectral resolution of JWST-MIRI MRS paired with high signal-to-noise data is essential to identify these species and determine their column densities and temperatures. The $Q$-branches of these molecules, including those of hot-bands, are particularly sensitive to temperature and column density. We find that the $^{12}$CO$_{2}$ emission in the GW Lup disk is very strong and is coming from optically thick emission at a temperature of $\sim$500 K. $^{13}$CO$_{2}$ is optically thinner, and tracing a larger emitting area than $^{12}$CO$_{2}$ and deeper into the disk based on a lower temperature of $\sim$200 K. The derived $N_{\rm{CO_{2}}}$/$N_{\rm{H_{2}O}}$ ratio of 0.25 is relatively high compared to other targets determined from \textit{Spitzer}-IRS data. This abundance ratio may be due to an inner cavity in the gas and dust with a radius in between the H$_{2}$O and CO$_{2}$ snowlines and/or a dust trap at this location keeping water-rich icy grains from passing the water snowline. This paper demonstrates the unique ability of JWST to probe inner disk structures and chemistry through weak, previously unseen molecular features.

Line-intensity mapping (LIM) is an emerging technique to probe the large-scale structure of the Universe. By targeting the integrated intensity of specific spectral lines, it captures the emission from all sources and is sensitive to the astrophysical processes that drive galaxy evolution. Relating these processes to the underlying distribution of matter introduces observational and theoretical challenges, such as observational contamination and highly non-Gaussian fields, which motivate the use of simulations to better characterize the signal. In this work we present SkyLine, a computational framework to generate realistic mock LIM observations that include observational features and foreground contamination, as well as a variety of self-consistent tracer catalogs. We apply our framework to generate realizations of LIM maps from the MultiDark Planck 2 simulations coupled to the UniverseMachine galaxy formation model. We showcase the potential of our scheme by exploring the voxel intensity distribution and the power spectrum of emission lines such as 21 cm, CO, CII, and Lyman-$\alpha$, their mutual cross-correlations, and cross-correlations with galaxy clustering. We additionally present cross-correlations between LIM and sub-millimeter extragalactic tracers of large-scale structure such as the cosmic infrared background and the thermal Sunyaev-Zel'dovich effect, as well as quantify the impact of galactic foregrounds, line interlopers and instrument noise on LIM observations. These simulated products will be crucial in quantifying the true information content of LIM surveys and their cross-correlations in the coming decade, and to develop strategies to overcome the impact of contaminants and maximize the scientific return from LIM experiments.

We present a numerical implementation of the guiding center approximation to describe the relativistic motion of charged test particles in the PLUTO code for astrophysical plasma dynamics. The guiding center approximation (GCA) removes the time step constraint due to particle gyration around magnetic field lines by following the particle center of motion rather than its full trajectory. The gyration can be detached from the guiding center motion if electromagnetic fields vary sufficiently slow compared to the particle gyration radius and period. Our implementation employs a variable step-size linear multistep method, more efficient when compared to traditional one-step Runge Kutta schemes. A number of numerical benchmarks is presented in order to assess the validity of our implementation.

Dark matter (DM) in the Milky Way halo may annihilate or decay to photons, producing monochromatic gamma rays. We search for DM-induced spectral lines using 14 years of data from the Large Area Telescope onboard the Fermi Gamma-ray Space Telescope ($\textit{Fermi}$-LAT) between $10\,\mathrm{GeV}$ and $2\,\mathrm{TeV}$ in the inner Milky Way leveraging both the spatial and spectral morphology of an expected signal. We present new constraints as strong as $\langle \sigma v \rangle \lesssim 6\times 10^{-30}\, \mathrm{cm}^3/\mathrm{s}$ for the two-to-two annihilations and $\tau \gtrsim 10^{30}\,\mathrm{s}$ for one-to-two decays, representing leading sensitivity between $10\,\mathrm{GeV}$ and $\sim$$500\,\mathrm{GeV}$. We consider the implications of our line-constraints on the Galactic Center Excess (GCE), which is a previously-observed excess of continuum $\sim$GeV gamma-rays that may be explained by DM annihilation. The Higgs portal and neutralino-like DM scenarios, which have been extensively discussed as possible origins of the GCE, are constrained by our work because of the lack of observed one-loop decays to two photons. More generally, we interpret our null results in a variety of annihilating and decaying DM models, such as neutralinos, gravitinos, and glueballs, showing that in many cases the line search is more powerful than the continuum, despite the continuum annihilation being at tree level.

We investigate the contribution of higher spin particles in the signal of direct detection experiments searching for dark matter. We consider a bosonic or fermionic higher spin dark matter (HSDM) candidate which interacts with the Standard Model via a dark U(1) mediator. For a particular subclass of interactions, spin-polarized targets may be used for spin determination: The angular dependence of scatterings can distinguish integer (spin-$s$) vs. half-integer (spin-$s + 1/2$), while the recoil energy dependence of the signal determines $s$. We consider also the signal of a supersymmetric higher spin dark sector, which suggests a characteristic signal (''SUSY Rilles'') for directional direct detection.

A driven evolution of a magnetic field has three aspects: field line topology, magnetic energy, and magnetic helicity. An ideal evolution with minimal energy input can produce magnetic field line chaos, which makes the preservation of field-line topology exponentially sensitive to non-ideal effects on an evolution timescale, but has no direct effect on energy or helicity dissipation. Resistive dissipation of the power input requires highly localized current densities $j\approx vB/\eta$, where $v$ is the velocity given by the evolution drive and $\eta$ is the plasma resistivity. Dissipation of the magnetic helicity input cannot be balanced when the magnetic Reynolds number is large compared to unity. Current densities $j\approx vB/\eta$ are consistent with those required to produce the solar corona with the observed height of the transition region by the Dreicer electron runaway effect.

The interaction of defects can lead to a phenomenon of erasure. During this process, a lower-dimensional object gets absorbed and dissolved by a higher-dimensional one. The phenomenon is very general and has a wide range of implications, both cosmological and fundamental. In particular, all types of strings, such as cosmic strings, QCD flux tubes, or fundamental strings, get erased when encountering a defect, either solitonic or a $D$-brane that deconfines their fluxes. This leads to a novel mechanism of cosmic string break-up, accompanied by gravitational and electromagnetic radiations. The arguments based on loss of coherence and the entropy count suggest that the erasure probability is very close to one, and strings never make it through the deconfining layer. We confirm this by a numerical simulation of the system, which effectively captures the essence of the phenomenon: a $2+1$-dimensional problem of interaction between a Nielsen-Olesen vortex of a $U(1)$ Higgs model and a domain wall inside which the $U(1)$ gauge group is unHiggsed and the magnetic flux is deconfined. In accordance with the entropy argument, in our simulation, the vortex never makes it across the wall.

We investigate rarely explored details of supercooled cosmological first-order phase transitions at the electroweak scale, which may lead to strong gravitational wave signals or explain the cosmic baryon asymmetry. The nucleation temperature is often used in phase transition analyses, and is defined through the nucleation condition: on average one bubble has nucleated per Hubble volume. We argue that the nucleation temperature is neither a fundamental nor essential quantity in phase transition analysis. We illustrate scenarios where a transition can complete without satisfying the nucleation condition, and conversely where the nucleation condition is satisfied but the transition does not complete. We also find that simple nucleation heuristics, which are defined to approximate the nucleation temperature, break down for strong supercooling. Thus, studies that rely on the nucleation temperature $\unicode{x2014}$ approximated or otherwise $\unicode{x2014}$ may misclassify the completion of a transition. Further, we find that the nucleation temperature decouples from the progress of the transition for strong supercooling. We advocate use of the percolation temperature as a reference temperature for gravitational wave production, because the percolation temperature is directly connected to transition progress and the collision of bubbles. Finally, we provide model-independent bounds on the bubble wall velocity that allow one to predict whether a transition completes based only on knowledge of the bounce action curve. We apply our methods to find empirical bounds on the bubble wall velocity for which the physical volume of the false vacuum decreases during the transition. We verify the accuracy of our predictions using benchmarks from a high temperature expansion of the Standard Model and from the real scalar singlet model.

In this paper, we find new scalarized black holes by coupling a scalar field with the Gauss-Bonnet invariant in Teleparallel gravity. The Teleparallel formulation of this theory uses torsion instead of curvature to describe the gravitational interaction and it turns out that, in this language, the usual Gauss-Bonnet term in four dimensions, decays in two distinct boundary terms, the Teleparallel Gauss-Bonnet invariants. Both can be coupled individually, or in any combination, to a scalar field, to obtain a Teleparallel Gauss-Bonnet extension of the Teleparallel equivalent of general relativity. The theory we study contains the familiar Riemannian Einstein-Gauss-Bonnet gravity theory as a particular limit and offers a natural extension, in which scalarization is triggered by torsion and with new interesting phenomenology. We demonstrate numerically the existence of asymptotically flat scalarized black hole solutions and show that, depending on the choice of coupling of the boundary terms, they can have a distinct behaviour compared to the ones known from the usual Einstein-Gauss-Bonnet case. More specifically, non-monotonicity of the metric functions and the scalar field can be present, a feature that was not observed until now for static scalarized black hole solutions.

We study self-gravitating bound states of a complex vector field, known as Proca stars, with a new type of quartic-order self-interaction which does not exist in the case of either a complex scalar field or a real vector field. Depending on the sign of the coupling constant, this quartic self-interaction can yield a distinct feature of Proca stars from the previously investigated self-interaction of the vector field. We find that self-gravitating solutions can be so compact that the photon sphere could form. However, we also show that the self-interaction gives rise to a ghost instability for the stars whose compactness is close that for the formation of a photon sphere, which might invalidate the formation of the photon sphere.

This work is devoted to the capabilities analysis of constellation and small spacecraft developed using CubeSat technology to solve promising problems of the Earth remote sensing in the area of greenhouse gases emissions. This paper presents the scientific needs for such tasks, followed by descriptions and discussions of the micro-technology application both in the small satellite platform design and in the payload design. The overview of analogical spacecraft is carried out. The design of a new spacecraft for determination the oxygen and carbon dioxide concentration in the air column along the line of sight of the spacecraft when it illuminated by reflected sunlight is introduced. A mock-up of the device was made for greenhouse gases remote sensing a Fourier Transform Infrared (FTIR) spectroradiometer is placed in the small spacecraft design. The results of long-term measurements of greenhouse gas concentrations using the developed Fourier spectrometer mock-up is presented.

The age of exascale computing has arrived and the risks associated with neutron and other atmospheric radiation are becoming more critical as the computing power increases, hence, the expected Mean Time Between Failures will be reduced because of this radiation. In this work, a new and detailed calculation of the neutron flux for energies above 50 MeV is presented. This has been done by using state-of-the-art Monte Carlo astroparticle techniques and including real atmospheric profiles at each one of the next 23 exascale supercomputing facilities. Atmospheric impact in the flux and seasonal variations were observed and characterised, and the barometric coefficient for high-energy neutrons at each site was obtained. With these coefficients, potential risks of errors associated with the increase in the flux of energetic neutrons, such as the occurrence of single event upsets or transients, and the corresponding failure-in-time rates, can be anticipated just by using the atmospheric pressure before the assignation of resources to critical tasks at each exascale facility. For more clarity, examples about how the rate of failures is affected by the cosmic rays are included, so administrators will better anticipate which more or less restrictive actions could take for overcoming errors.

The aim of the research is to look into a new solution for isotropic compact stars in the context of the $f(R,\,T)$ theory of gravity. We used the Buchdahl [H.A. Buchdahl, Phys. Rev. {\bf 116} (1959) 1027] metric potentials as input to deal with the field equations in the $f(R,\,T)$ framework. For different values of the coupling parameter $\chi$, graphical representation of the model parameters have been shown to canvass the analytical results more clearly. Interestingly, we have proven that for $\chi=0$, the standard General Relativity (GR) results can be recovered. A comparison of our obtained solutions with the GR results is also discussed. To study the effect of the coupling parameter $\chi$, the numerical values of the different physical variables have been tabulated for the values of the coupling parameter $\chi=0,\,0.25,\,0.5,\,0.75,1,\,1.25$. We used the compact stars candidate LMC X-4 with mass$=(1.04 \pm 0.09)M_{\odot}$; Radius $= 8.301_{-0.2}^{+0.2}$ km. respectively, for graphical analysis. To determine the physical acceptability of the model, we looked into the necessary physical properties such as energy conditions, causality, hydrostatic equilibrium, and pressure-density ratio etc. and found that our system satisfies all of these criteria, indicating that the model is physically reasonable.

It is known for some time that the cosmological bound on the invisible axion scale can be avoided by an early phase of strong QCD. While most approaches rely on theories where the strong coupling constant is determined through the expectation value of some scalar field, we show that this mechanism can also be implemented into the benchmark KSVZ and DFSZ models when the early phase of strong QCD emerges by the modification of the running coupling during inflation. For both models the physics that are responsible for making QCD strong do not displace the axions minimum by too much, so that the efficiency of the relaxation is controlled by parameters of the theory and the number of inflationary e-folds. In particular, we consider the case of very efficient relaxation where the axion abundance is dominated by inflationary quantum and post-inflationary thermal fluctuations. Within this situation we identify the parameter space compatible with all cosmological constrains and derive conditions on the reheating temperature and the QCD scale during inflation that result in the axion making up all the dark matter. Due to duality below the Peccei-Quinn scale and a minor influence of the KSVZ and DFSZ fields on the running, our findings also apply to a minimal 2-form implementation of the axion.

Nuclei that are unstable with respect to double beta decay are investigated in this work for a novel Dark Matter (DM) direct detection approach. In particular, the diagram responsible for the neutrinoless double beta decay will be considered for the possible detection technique of a Majorana DM fermion inelastically scattering on a double beta unstable nucleus, stimulating its decay. The exothermic nature of the stimulated double beta decay would allow the direct detection also of a light DM fermion, a class of DM candidates that are difficult or impossible to investigate with the traditional elastic scattering techniques. The expected signal distribution for different DM masses and the upper limits on the nucleus scattering cross sections, are shown and compared with the existing data for the case of $^{136}$Xe nucleus.

The high-energy atmospheric neutrino flux is dominated by neutrinos from the decays of charmed hadrons produced in the forward direction by cosmic ray interactions with air nuclei. We evaluate the charm contributions to the prompt atmospheric neutrino flux as a function of the center-of-mass energy $\sqrt{s}$ of the hadronic collision and of the center-of-mass rapidity $y$ of the produced charm hadron. Uncertainties associated with parton distribution functions are also evaluated as a function of $y$. We find that the $y$ coverage of LHCb for forward heavy-flavour production, complemented by the angular coverage of present and future forward neutrino experiments at the LHC, bracket the most interesting $y$ regions for the prompt atmospheric neutrino flux. At $\sqrt{s}=14$ TeV foreseeen for the HL-LHC phase, nucleon collisions in air contribute to the prompt neutrino flux prominently below $E_\nu\sim 10^7$ GeV. Measurements of forward charm and/or forward neutrinos produced in hadron collisions up to $\sqrt{s}=100$ TeV, which might become possible at the FCC, are relevant for the prompt atmospheric neutrino flux up to $E_\nu=10^8$ GeV and beyond.

We search for stochastic gravitational wave background generated by domain wall networks in the Data Release-2 of Parkes Pulsar Timing Array and find that the observed strong common power-law process can be explained by domain wall networks for the wall tension $\sigma_{\rm{DW}}\sim (29-414~\rm{TeV})^3$ and the wall-decay temperature $T_d\sim 26-363~\rm{MeV}$. Interestingly, the same parameter region can largely alleviate the Hubble tension, if the free particles generated from domain wall networks further decay into dark radiation. In addition, the preferred parameter space corresponds to the axion mass range $m_a \sim 10^{-13}-10^{-8}\ {\rm eV}$ for QCD axion. On the other hand, assuming that the common power-law process is not due to domain wall networks, we can put stringent constraints on the wall tension and decay temperature around the energy scale of QCD phase transition.

We propose a mass-varying dark matter (MVDM) model consisting of a scalar field and a fermionic field interacting via a simple Yukawa coupling, and containing an exponential self-interaction potential for the scalar field. Analyzing the evolution of this coupled scalar-fermion system in an expanding Universe, we find that it initially behaves like radiation but then undergoes a phase transition after which it behaves like pressureless dark matter. The two free parameters of this model are the temperature at which the phase transition occurs and the current relic density of dark matter; the mass of the dark matter particle, given by the mass of the fermion, is derived from this. For a phase transition temperature between 10 MeV and $10^7$ GeV, the current dark matter relic density is achieved for a fermion mass in the range of 1 GeV to $10^9$ GeV. In this dark matter model, the scalar becomes a sub-dominant unclustered component of dark matter that can lower the amplitude of structure formation by up to a few percent. Another feature is that the mass-varying fermion component can lead to discrepant measurements of the current dark matter density of about ten percent inferred from early and late-time probes assuming LCDM.

The upcoming Hyper-Kamiokande (HyperK) experiment is expected to detect the Diffuse Supernova Neutrino Background (DSNB). This requires to ponder all possible sources of background. Sub-GeV dark matter (DM) which annihilates into neutrinos is a potential background that has not been considered so far. We simulate DSNB and DM signals, as well as backgrounds in the HyperK detector. We find that DM-induced neutrinos could indeed alter the extraction of the correct values of the parameters of interest for DSNB physics. Since the DSNB is an isotropic signal, and DM originates primarily from the Galactic centre, we show that this effect could be alleviated with an on-off analysis.

We develop a framework to derive consistency constraints on gravitational Regge amplitudes based on the finite energy sum rules (FESRs), which directly connect gravitational Regge amplitudes at a finite ultraviolet scale with infrared physics without suffering from super-Planckian physics. For illustration, we consider four-point scattering of an identical massless scalar coupled to gravity. First, we derive multiple FESRs without relying on the $s\text{-}t\text{-}u$ permutation invariance. We then make use of FESRs, crossing symmetry, and other principles such as unitarity, to derive bounds on the Regge parameters. The bounds result in infrared finite gravitational positivity bounds in four spacetime dimensions.

Horndeski gravity is the most general scalar-tensor theory with one scalar field leading to second-order Euler-Lagrange field equations for the metric and scalar field, and it is based on Riemannian geometry. In this paper, we formulate an analogue version of Horndeski gravity in a symmetric teleparallel geometry which assumes that both the curvature (general) and torsion are vanishing and gravity is only related to nonmetricity. Our setup requires that the Euler-Lagrange equations for not only metric and scalar field but also connection should be at most second order. We find that the theory can be always recast as a sum of the Riemannian Horndeski theory and new terms that are purely teleparallel. Due to the nature of nonmetricity, there are many more possible ways of constructing second-order theories of gravity. In this regard, up to some assumptions, we find the most general $k$-essence extension of Symmetric Teleparallel Horndeski. We also formulate a novel theory containing higher-order derivatives acting on nonmetricity while still respecting the second-order conditions, which can be recast as an extension of Kinetic Gravity Braiding. We finish our study by presenting the FLRW cosmological equations for our model.

We present a generalization of the curative initial data construction derived for equal-mass compact binaries in Refs.[1,2] to arbitrary mass ratios. We demonstrate how these improved initial data avoid substantial spurious artifacts in the collision dynamics of unequal-mass boson-star binaries in the same way as has previously been achieved with the simpler method restricted to the equal-mass case. We employ the improved initial data to explore in detail the impact of phase offsets in the coalescence of equal- and unequal-mass boson star binaries.