Papers for Friday, Jun 30 2023

The Laser Interferometer Space Antenna (LISA), due for launch in the mid 2030s, is expected to observe gravitational waves (GW)s from merging massive black hole binaries (MBHB)s. These signals can last from days to months, depending on the masses of the black holes, and are expected to be observed with high signal to noise ratios (SNR)s out to high redshifts. We have adapted the PyCBC software package to enable a template bank search and inference of GWs from MBHBs. The pipeline is tested on the LISA data challenge (LDC)'s Challenge 2a ("Sangria"), which contains MBHBs and thousands of galactic binaries (GBs) in simulated instrumental LISA noise. Our search identifies all 6 MBHB signals with more than $98\%$ of the optimal signal to noise ratio. The subsequent parameter inference step recovers the masses and spins within their $90\%$ confidence interval. Sky position parameters have 8 high likelihood modes which are recovered but often our posteriors favour the incorrect sky mode. We observe that the addition of GBs biases the parameter recovery of masses and spins away from the injected values, reinforcing the need for a global fit pipeline which will simultaneously fit the parameters of the GB signals before estimating the parameters of MBHBs.

Classical pulsating stars such as Cepheid and RR Lyrae variables exhibit well-defined Period-Luminosity relations at near-infrared wavelengths. Despite their extensive use as stellar standard candles, the effects of metallicity on Period-Luminosity relations for these pulsating variables, and in turn, on possible biases in distance determinations, are not well understood. We present ongoing efforts in determining accurate and precise metallicity coefficients of Period-Luminosity-Metallicity relations for classical pulsators at near-infrared wavelengths. For Cepheids, it is crucial to obtain a homogeneous sample of photometric light curves and high-resolution spectra for a wide range of metallicities to empirically determine metallicity coefficient and reconcile differences with the predictions of the theoretical models. For RR Lyrae variables, using their host globular clusters covering a wide range of metallicities, we determined the most precise metallicity coefficient at near-infrared wavelengths, which is in excellent agreement with the predictions of the horizontal branch evolution and stellar pulsation models.

The Lilly-Madau plot is commonly interpreted as the history of the cosmic star formation of the Universe by showing the co-moving star formation rate density (SFRD) over cosmic time. Therefore, the Lilly-Madau plot is not only sensitive to the star formation history (SFH) but also to the number density of galaxies. Assessing the Catalogue of Neighbouring Galaxies, we reconstruct the SFHs of galaxies located in the Local Volume (LV) based on delayed-$\tau$ and power-law SFH models. Galaxies with stellar masses of $M_{*} \gtrsim 10^{10}\,\rm{M_{\odot}}$ typically evolve according to the delayed-$\tau$ model by having first increasing followed by exponentially declining SFRs, while the majority of less massive star-forming galaxies has an almost constant or increasing SFH. Deducing the cosmic SFRD evolution of the LV reveals that the SFHs of local galaxies are inconsistent with the Lilly-Madau plot. The SFRDs of the LV are significantly lower at redshifts of $z \lesssim 3$ underestimating the Lilly-Madau peak at $z = 1.86$ by a factor of $2.16\pm0.32$ (delayed-$\tau$) and $5.90\pm0.88$ (power-law model). Assuming the delayed-$\tau$ model for galaxies with $M_{*} \geq 10^{10}\,\rm{M_{\odot}}$ and a power-law model for less massive galaxies, the SFRD is $2.22\pm0.33$ lower than measured at $z = 1.86$. This inconsistency between the evolution of the local and global SFRD has cosmological implications. Since the Lilly-Madau plot also constrains the cosmological matter field, the near-constancy of SFHs of LV galaxies could imply that the peak of the Lilly-Madau plot at $z = 1.86$ is the imprint of a $\approx~5$ Gpc-scale inhomogeneity.

Mergers of galaxies are a ubiquitous phenomenon in the Universe and represent a natural consequence of the ``bottom-up'' mass accumulation and galaxy evolution cosmological paradigm. It is generally accepted that the peak of AGN accretion activity occurs at nuclear separations of $\lesssim10$ kpc for major mergers. Here we present new NuSTAR and XMM-Newton observations for a subsample of mid-IR preselected dual AGN candidates in an effort to better constrain the column densities along the line-of-sight for each system. Only one dual AGN candidate, J0841+0101, is detected as a single, unresolved source in the XMM-Newton and NuSTAR imaging, while the remaining three dual AGN candidates, J0122+0100, J1221+1137, and J1306+0735, are not detected with NuSTAR; if these non-detections are due to obscuration alone, these systems are consistent with being absorbed by column densities of log($N_{\rm{H}}/\rm{cm}^{-2}$) $\geq$ 24.9, 24.8, and 24.6, which are roughly consistent with previously inferred column densities in these merging systems. In the case of J0841+0101, the analysis of the 0.3-30 keV spectra reveal a line-of-sight column density of $N_{\rm{H}}\gtrsim10^{24}$ cm$^{-2}$, significantly larger than the column densities previously reported for this system and demonstrating the importance of the higher signal-to-noise XMM-Newton spectra and access to the $>10$ keV energies via NuSTAR. Though it is unclear if J0841+0101 truly hosts a dual AGN, these results are in agreement with the high obscuring columns expected in AGNs in late-stage mergers.

Context: A non-linear relation between quasar monochromatic luminosities at 2500A and 2 keV holds at all observed redshifts and luminosities, and it has been used to derive quasar distances and to build a Hubble Diagram of quasars. The choice of the X-ray and UV indicators has so far been somewhat arbitrary, and has typically relied on photometric data. Aims: We want to determine the X-ray and UV proxies that provide the smallest dispersion of the relation, in order to obtain more precise distance estimates, and to confirm the reliability of the X-ray to UV relation as a distance indicator. Methods: We performed a complete UV spectroscopic analysis of a sample of $\sim$1800 quasars with SDSS optical spectra and XMM- Newton X-ray serendipitous observations. In the X-rays, we analysed the spectra of all the sample objects at redshift z $>$1.9, while we relied on photometric measurements at lower redshifts. As done in previous studies, we analysed the relation in small redshift bins, using fluxes instead of luminosities. Results: We show that the monochromatic fluxes at 1 keV and 2500A are, respectively, the best X-ray and UV continuum indicators among those that are typically available. We also find a tight relation between soft X-ray and Mg ii2800A line fluxes, and a marginal dependence of the X-ray to UV relation on the width of the Mg ii line. Conclusions: Our analysis suggests that the physical quantities that are more tightly linked to one another are the soft X-ray flux at $\sim$1 keV and the ionizing UV flux blueward of the Lyman limit. However, the "usual" monochromatic fluxes at 2 keV and 2500A estimated from photometric data provide an almost as-tight X-ray to UV relation, and can be used to derive quasar distances. The Hubble diagram obtained using spectroscopic indicators is fully consistent with the one presented in previous papers, based on photometric data.

Surface densities of gas, dust and stars provide a window into the physics of star-formation that, until the advent of high-resolution far-infrared/sub-millimeter observations, has been historically difficult to assess amongst dusty galaxies. To study the link between infrared (IR) surface densities and dust properties, we leverage the Atacama Large Millimetre/Submillimetre Array (ALMA) archive to measure the extent of cold dust emission in 15 $z\sim2$ IR selected galaxies selected on the basis of having available mid-IR spectroscopy from Spitzer. We use the mid-IR spectra to constrain the relative balance between dust heating from star-formation and active galactic nuclei (AGN), and to measure emission from Polycylic Aromatic Hydrocarbons (PAHs) -- small dust grains that play a key role in the photoelectric heating of gas. In general, we find that dust-obscured star-formation at high IR surface densities exhibits similar properties at low- and high-redshift, namely: local luminous IR galaxies have comparable PAH luminosity to total dust mass ratios as high-$z$ galaxies, and star-formation at $z\sim0-2$ is more efficient at high IR surface densities despite the fact that our sample of high$-z$ galaxies are closer to the main-sequence than local luminous IR galaxies. High star-formation efficiencies are coincident with a decline in the PAH/IR luminosity ratio reminiscent of the deficit observed in far-infrared fine-structure lines. Changes in the gas and dust conditions arising from high star-formation surface densities might help drive the star-formation efficiency up. This could help explain high efficiencies needed to reconcile star-formation and gas volume densities in dusty galaxies at cosmic noon.

Artificial neural network emulators have been demonstrated to be a very computationally efficient method to rapidly generate galaxy spectral energy distributions (SEDs), for parameter inference or otherwise. Using a highly flexible and fast mathematical structure, they can learn the nontrivial relationship between input galaxy parameters and output observables. However, they do so imperfectly, and small errors in flux prediction can yield large differences in recovered parameters. In this work, we investigate the relationship between an emulator's execution time, uncertainties, correlated errors, and ability to recover accurate posteriors. We show that emulators can recover consistent results to traditional fits, with precision of $25\!-\!40\%$ in posterior medians for stellar mass, stellar metallicity, star formation rate, and stellar age. We find that emulation uncertainties scale with an emulator's width $N$ as $\propto N^{-1}$ while execution time scales as $\propto N^2$, resulting in an inherent tradeoff between execution time and emulation uncertainties. We also find that emulators with uncertainties smaller than observational uncertaities are able to recover accurate posteriors for most parameters without a significant increase in catastrophic outliers. Furthermore, we demonstrate that small architectures can produce flux residuals that have significant correlations, which can create dangerous systematic errors in colors. Finally, we show that the distributions chosen for generating training sets can have a large effect on emulators' ability to accurately fit rare objects. Selecting the optimal architecture and training set for an emulator will minimize the computational requirements for fitting near-future large-scale galaxy surveys.

Gravitational potentials of the Milky Way and extragalactic structures can influence the propagation of the cosmic neutrino background (CNB). Of particular interest to future CNB observatories, such as PTOLEMY, is the CNB number density on Earth. In this study, we have developed a simulation framework that maps the trajectories of relic neutrinos as they move through the local gravitational environment. The potentials are based on the dark matter halos found in state-of-the-art cosmological N-body simulations, resulting in a more nuanced and realistic input than the previously employed analytical models. We find that the complex dark matter distributions, along with their dynamic evolution, influence the abundance and anisotropies of the CNB in ways unaccounted for by earlier analytical methods. Importantly, these cosmological simulations contain multiple instances of Milky Way-like halos that we employ to model a variety of gravitational landscapes. Consequently, we notice a variation in the CNB number densities that can be primarily attributed to the differences in the masses of these individual halos. For neutrino masses between $0.01$ and $0.3$ eV, we note clustering factors within the range of $1+\mathcal{O}(10^{-3})$ to $1+\mathcal{O}(1)$. Furthermore, the asymmetric nature of the underlying dark matter distributions within the halos results in not only overdense, but intriguingly, underdense regions within the full-sky anisotropy maps. Gravitational clustering appears to have a significant impact on the angular power spectra of these maps, leading to orders of magnitude more power on smaller scales beyond multipoles of $\ell = 3$ when juxtaposed against predictions by primordial fluctuations. We discuss how our results reshape our understanding of relic neutrino clustering and how this might affect observability of future CNB observatories such as PTOLEMY.

We present a generalized analytical Bayesian framework for calculating the occurrence rate of steady emission (or absorption) in astrophysical objects. As a proof-of-concept, we apply this framework to non-flaring quiescent radio emission in ultracool ($\leq$ M7) dwarfs. Using simulations, we show that our framework recovers the simulated radio occurrence rate to within 1-5% for sample sizes of 10-100 objects when averaged over an ensemble of trials and simulated occurrence rates for our assumed luminosity distribution models. In contrast, existing detection rate studies may under-predict the simulated rate by 51-66% because of sensitivity limits. Using all available literature results for samples of 82 ultracool M dwarfs, 74 L dwarfs, and 23 T/Y dwarfs, we find that the maximum-likelihood quiescent radio occurrence rate is between $15^{+4}_{-4}$ - $20^{+6}_{-5}$%, depending on the luminosity prior that we assume. Comparing each spectral type, we find occurrence rates of $17^{+9}_{-7}$ - $25^{+13}_{-10}$% for M dwarfs, $10^{+5}_{-4}$ - $13^{+7}_{-5}$% for L dwarfs, and $23^{+11}_{-9}$ - $29^{+13}_{-11}$% for T/Y dwarfs. We rule out potential selection effects and speculate that age and/or rotation may account for tentative evidence that the quiescent radio occurrence rate of L dwarfs may be suppressed compared to M and T/Y dwarfs and phenomenon. Finally, we discuss how we can harness our occurrence rate framework to carefully assess the possible physics that may be contributing to observed occurrence rate trends.

Millimeter sized dust grains experience radial velocities exceeding the gas velocities by orders of magnitude. The viscous evolution of the accretion disk adds disk material onto the central star's convective envelope, influencing its elemental abundances, [X/H]. At the same time, the envelope mass shrinks over time, amplifying the rate of abundance change. Therefore, the elemental abundances of the star are sensitive to disk processes that alter the composition and timing of disk accretion. We perform numerical 1D log-radial simulations integrating the disk advection-diffusion equation, accounting for phase transitions of chemical species at the evaporation fronts. They reveal a peak of refractory abundance within the first 2 Myr of $\Delta\mathrm{[X/H]}\sim 5\times 10^{-2}$ if grain growth is significant, but subsequent accretion diminishes previous refractory abundance increases for long-lived disks. Planet formation can reduce the abundance of dust species whose evaporation fronts lie within the planet's orbit by opening a gap and consequently blocking inward drifting pebbles. We expect the accretion of the Solar protoplanetary disk with Jupiter present to have changed the Sun's elemental abundances by ${\sim}10^{-2}$ throughout its lifetime. These considerations are also applied to the HD106515 wide binary system. We find that measurements of $\Delta\mathrm{[X/H]}$ are in reasonable agreement with results from simulations where the observed giant planet around HD106515 A is included and if HD106515B's disk formed planetesimals more efficiently. Simulations where the planet formed inside the water ice line are more favorable to agree with observations. Even though the general changes in the stellar abundances due to disk accretion are small, they are detectable at current sensitivities, indicating that the here presented methods can be used to constrain the planet formation pathway.

The orders of magnitude variation in lithium abundances of evolved stars have long been a puzzle. Diluted signals, ambiguous evolutionary states and unknown masses have made it challenging to both map the expected lithium signals and explain the anomalously lithium-rich stars. We show here using a set of asteroseismically characterized evolved stars that the base lithium abundance in red giant stars is mass dependent, with higher mass stars having higher `normal' lithium abundances, while highly lithium enhanced stars may cluster around 0.8 or 1.8 M$_\odot$. We confirm previous studies that have shown that lithium enhancement and rapid rotation are often coincident, but find that the actual correlation between lithium abundance and the rotation rate, whether surface rotation, internal rotation, or radial differential rotation, is weak. Our data support previous assertions that most lithium rich giants are in the core-helium burning phase. We also note a tentative correlation between the highest lithium abundances and unusual carbon to nitrogen ratios, which is suggestive of binary interactions, though we find no simple correlation between lithium richness and indicators of binarity.

We have performed a 5-colour surface photometric study of the high-galactic-latitude area of dark nebula LDN 1642. Scattered light properties are presented of diffuse, translucent and opaque dust over the range of 3500 -- 5500 A. Far infrared absolute photometry at 200 um improves the precision of and provides a zero point to the extinction. The intensity of the scattered light depends on dust column density in a characteristic way: for optically thin dust the intensity first increases linearly, then turns to a saturation value; at still larger extinctions the intensity turns down to a slow decrease. The $A_V$ value of the saturated intensity maximum shifts in a systematic way, from $A_V\approx$ 1.5 mag at 3500 A, to $\sim 3$ mag at 5500 A. The intensity curves offer a straight-forward explanation for the behaviour of the scattered-light colours. At the intensity peak the colour agrees with the integrated starlight colour, while it is bluer at the low- and redder at the high-column-density side of the peak, respectively. These colour changes are a direct consequence of the wavelength dependence of the extinction. We have compared the colours of the LDN 1642 area with other relevant observational studies: high-latitude diffuse/translucent clouds, wide-field cirrus dust; and externally illuminated AGB-star envelopes. For extragalactic low-surface-brightness sources cirrus is an unwanted foreground contaminant. Our results for cirrus colours can help to distinguish cases where a diffuse plume or stream, apparently associated with a galaxy or a group or cluster, is more likely a local cirrus structure. Keywords: ISM: dust, extinction -- ISM: clouds, individual LDN 1642 -- Galaxy: solar neighbourhood -- Astronomical instruments, methods and techniques: methods -- Physical data and processes: scattering

Characterizing rocky exoplanet atmospheres is a key goal of exoplanet science, but interpreting such observations will require understanding the stellar UV irradiation incident on the planet from its host star. Stellar UV mediates atmospheric escape, photochemistry, and planetary habitability, and observations of rocky exoplanets can only be understood in the context of the UV SED of their host stars. Particularly important are SEDs from observationally favorable but poorly understood low-mass M-dwarf stars, which are the only plausible targets for rocky planet atmospheric characterization for the next 1-2 decades. In this work, we explore the utility of AstroSat UVIT for the characterization of the UV SEDs of low-mass stars. We present observations of the nearby M0 star HIP 23309 in the FUV and NUV gratings of UVIT. Our FUV spectra are consistent with contemporaneous HST data and our NUV spectra are stable between orbits, suggesting UVIT is a viable tool for the characterization of the SEDs of low-mass stars. We apply our measured spectra to simulations of photochemistry and habitability for a hypothetical rocky planet orbiting HIP 23309 and elucidate the utility and limitations of UVIT in deriving UV SEDs of M-dwarf exoplanet hosts. Our work validates UVIT as a tool to complement HST in the characterization of exoplanet host stars and carries implications for its successor missions like INSIST.

We present bolometric luminosities, black hole masses and Eddington ratios for 42 luminous quasars at z>6 using high signal-to-noise ratio VLT/X-Shooter spectra, acquired in the enlarged ESO Large Programme XQR-30. In particular, we derive bolometric luminosities from the rest-frame 3000 A, luminosities using a bolometric correction from the literature, and the black hole masses by modelling the spectral regions around the CIV 1549A and the MgII 2798A emission lines, with scaling relations calibrated in the local universe. We find that the black hole masses derived from both emission lines are in the same range, and the scatter of the measurements agrees with expectations from the scaling relations. The MgII-derived masses are between ~(0.8-12) x 10^9 Msun, and the derived Eddington ratios are within ~0.13-1.73, with a mean (median) of 0.84 (0.72). By comparing the total sample of quasars at z>5.8, from this work and from the literature, to a bolometric luminosity distribution-matched sample at z~1.5, we find that quasars at high redshift host slightly less massive black holes which accrete slightly more rapidly than at lower-z, with a difference in the mean Eddington ratios of the two samples of ~0.27, in agreement with recent literature work.

The third data release (DR3) of Gaia has provided a five-fold increase in the number of radial velocity measurements of stars, as well as a stark improvement in parallax and proper motion measurements. To help with studies that seek to test models and interpret Gaia DR3, we present nine Gaia synthetic surveys, based on three solar positions in three Milky-Way-mass galaxies of the Latte suite of the Fire-2 cosmological simulations. These synthetic surveys match the selection function, radial velocity measurements, and photometry of Gaia DR3, adapting the code base Ananke, previously used to match the Gaia DR2 release in Sanderson et al. 2020. The synthetic surveys are publicly available and can be found at this http URL Similarly to the previous release of Ananke, these surveys are based on cosmological simulations and thus able to model non-equilibrium dynamical effects, making them a useful tool in testing and interpreting Gaia DR3.

Compact object binaries (a black hole or a neutron star orbiting a non-degenerate stellar companion) are key to our understanding of late massive star evolution, in addition to being some of the best probes of extreme gravity and accretion physics. Gaia has opened the door to astrometric studies of these systems, enabling geometric distance measurements, kinematic estimation, and the ability to find new previously unknown systems through measurement of binary orbital elements. Particularly puzzling are newly found massive black holes in wide orbits (~AU or more) whose evolutionary history is difficult to explain. Astrometric identification of such binaries is challenging for Gaia, with only two such examples currently known. Roman's enormous grasp, superb sensitivity, sharp PSF and controlled survey strategy can prove to be a game-changer in this field, extending astrometric studies of compact object binaries several mag deeper than Gaia. We propose to use the microlensing Galactic Bulge Time Domain Survey to identify new wide-orbit black hole compact object binaries, determine their prevalence and their spatial distribution, thus opening up new parameter space in binary population studies.

We propose a novel statistical method to extend Fermi-LAT catalogues of high-latitude $\gamma$-ray sources below their nominal threshold. To do so, we rely on a recent determination of the differential source-count distribution of sub-threshold sources via the application of deep learning methods to the $\gamma$-ray sky. By simulating ensembles of synthetic skies, we assess quantitatively the likelihood for pixels in the sky with relatively low-test statistics to be due to sources. Besides being useful to orient efforts towards multi-messenger and multi-wavelength identification of new $\gamma$-ray sources, we expect the results to be especially advantageous for statistical applications such as cross-correlation analyses.

A relatively new concern for the forthcoming massive spectroscopic sky surveys is the impact of contamination from low earth orbit satellites. Several hundred thousand of these satellites are licensed for launch in the next few years and it has been estimated that, in some cases, up to a few percent of spectra could be contaminated when using wide field, multi-fiber spectrographs. In this paper, a multi-staged approach is used to assess the practicality and limitations of identifying and minimizing the impact of satellite contamination in a WEAVE-like stellar spectral survey. We develop a series of convolutional-network based architectures to attempt identification, stellar parameter and chemical abundances recovery, and source separation of stellar spectra that we artificially contaminate with satellite (i.e. solar-like) spectra. Our results show that we are able to flag 67% of all contaminated sources at a precision level of 80% for low-resolution spectra and 96% for high-resolution spectra. Additionally, we are able to remove the contamination from the spectra and recover the clean spectra with a $<$1% reconstruction error. The errors in stellar parameter predictions reduce by up to a factor of 2-3 when either including contamination as an augmentation to a training set or by removing the contamination from the spectra, with overall better performance in the former case. The presented methods illustrate several machine learning mitigation strategies that can be implemented to improve stellar parameters for contaminated spectra in the WEAVE stellar spectroscopic survey and others like it.

We present MAXI and NuSTAR observations of the Be X-ray binary, MAXI J0655-013, in outburst. NuSTAR observed the source once early in the outburst, when spectral analysis yields a bolometric (0.1--100 keV), unabsorbed source luminosity of $L_{\mathrm{bol}}=5.6\times10^{36}\mathrm{erg\,s^{-1}}$, and a second time 54 days later, by which time the luminosity dropped to $L_{\mathrm{bol}}=4\times10^{34}\,\mathrm{erg\,s^{-1}}$ after first undergoing a dramatic increase. Timing analysis of the NuSTAR data reveals a neutron star spin period of $1129.09\pm0.04$ s during the first observation, which decreased to $1085\pm1$ s by the time of the second observation, indicating spin-up due to accretion throughout the outburst. Furthermore, during the first NuSTAR observation, we observed quasiperiodic oscillations with centroid frequency $\nu_0=89\pm1$ mHz, which exhibited a second harmonic feature. By combining the MAXI and NuSTAR data with pulse period measurements reported by Fermi/GBM, we are able to show that apparent flaring behavior in the MAXI light-curve is an artifact introduced by uneven sampling of the pulse profile, which has a large pulsed fraction. Finally, we estimate the magnetic field strength at the neutron star surface via three independent methods, invoking a tentative cyclotron resonance scattering feature at $44$ keV, QPO production at the inner edge of the accretion disk, and spin-up via interaction of the neutron star magnetic field with accreting material. Each of these result in a significantly different value. We discuss the strengths and weaknesses of each method and infer that MAXI J0655-013 is likely to have a high surface magnetic field strength, $B_{s}>10^{13}$ G.

The Local Group (LG) consists of two giant spiral galaxies, the Milky Way (MW) and Andromeda (M31), and several smaller galaxies. The MW and M31 are approaching each other at a radial velocity of about $-109\,$km\,s$^{-1}$. Observational evidence suggests that there is an overall infalling motion of gas and galaxies in the LG, dominated by the dynamics of its two main members. From our perspective, this flow imprints a velocity dipole pattern in the sky when Galactic rotation is removed. We investigate the kinematic properties of gas and galaxies in the LG using a suite of high-resolution simulations performed by the {\sc Hestia} (High-resolution Environmental Simulations of The Immediate Area) collaboration. Our simulations include the correct cosmography surrounding LG-like regions. We build sky maps from the local, Galactic and LG standard of rest reference frames. Our findings show that the establishment of a radial velocity dipole near the preferred barycentre direction is a natural outcome of simulation kinematics for material \textit{outside} the MW virial radius after removing galaxy rotation when the relative radial velocity of MW and M31 is similar to the observed value. These results favour a scenario where gas and galaxies stream towards the LG barycentre, producing the observed velocity dipole.

The detection and accurate astrometry of fast-moving near-Earth objects (NEOs) has been a challenge for the follow-up community. Their fast apparent motion results in streaks in sidereal images, thus affecting the telescope's limiting magnitude and astrometric accuracy. A widely adopted technique to mitigate trailing losses is non-sidereal tracking, which transfers the streaking to background reference stars. However, no existing publicly available astrometry software is configured to detect such elongated stars. We present Astreaks, a streaking source detection algorithm, to obtain accurate astrometry of NEOs in non-sidereal data. We validate the astrometric accuracy of Astreaks on 371 non-sidereally tracked images for 115 NEOs with two instrument set-ups of the GROWTH-India Telescope. The observed NEOs had V-band magnitude in the range [15, 22] with proper motion up to 140$^{\prime\prime}$/min, thus resulting in stellar streaks as high as 6.5$^\prime$ (582 pixels) in our data. Our method obtained astrometric solutions for all images with 100% success rate. The standard deviation in Observed-minus-Computed (O-C) residuals is 0.52$^{\prime\prime}$, with O-C residuals <2$^{\prime\prime}$(<1$^{\prime\prime}$) for 98.4% (84.4%) of our measurements. These are appreciable, given the pixel scale of $\sim$0.3$^{\prime\prime}$ and $\sim$0.7$^{\prime\prime}$ of our two instrument set-ups. This demonstrates that our modular and fully-automated algorithm helps improve the telescope system's limiting magnitude without compromising astrometric accuracy by enabling non-sidereal tracking on the target. This will help the NEO follow-up community cope with the accelerated discovery rates and improved sensitivity of the next-generation NEO surveys. Astreaks has been made available to the community under an open-source license.

An unidentified 3.55 keV X-ray line in stacked spectra of galaxies and clusters raises the interesting possibility that it originates from the decay of sterile neutrino dark matter. In this work, we explore mixed sterile neutrino dark matter models that combine cold dark matter and warmer sterile neutrino dark matter produced through lepton number-driven active-to-sterile neutrino transformation. We analyze the sensitivity of the sterile neutrino spectra on active-sterile mixing and on initial neutrino lepton numbers. Furthermore, we assess the viability of these models with estimates of the number of subhalos formed as the host sites of satellite galaxies.

The B[e] phenomenon is manifested by a heterogeneous group of stars surrounded by gaseous and dusty circumstellar envelopes with similar physical conditions. Among these stars, the FS CMa-type objects are suspected to be binary systems, which could be experiencing or have undergone a mass-transfer process that could explain the large amount of material surrounding them. We aim to contribute to the knowledge of a recently confirmed binary, MWC 645, which could be undergoing an active mass-transfer process. We present near-infrared and optical spectra, identify atomic and molecular spectral features, and derive different quantitative properties of line profiles. Based on publicly available photometric data, we search for periodicity in the light curve and model the spectral energy distribution. We have detected molecular bands of CO in absorption at 1.62 $\mu$m and 2.3 $\mu$m for the first time. We derive an upper limit for the effective temperature of the cool binary component. We found a correlation between the enhancement of the H$\alpha$ emission and the decrease in optical brightness that could be associated with mass-ejection events or an increase in mass loss. We outline the global properties of the envelope, possibly responsible for brightness variations due to a variable extinction, and briefly speculate on different possible scenarios.

We carried out 3D smoothed particle hydrodynamics simulations of the common envelope binary interaction using the approximation of Bowen to calculate the dust opacity in order to investigate the resulting dust-driven accelerations. We have simulated two types of binary star: a 1.7 and a 3.7 $M_{\odot}$ thermally-pulsating, asymptotic giant branch stars with a 0.6 $M_{\odot}$ companion. We carried out simulations using both an ideal gas and a tabulated equations of state, with the latter considering the recombination energy of the envelope. We found that the dust-driven wind leads to a relatively small increase in the unbound gas, with the effect being smaller for the tabulated equation of state simulations and for the more massive primary. Dust acceleration does contribute to envelope expansion with only a slightly elongated morphology, if we believe the results from the tabulated equation of state as more reliable. The Bowen opacities in the outer envelopes of the two models, at late times, are large enough that the photosphere of the post-inspiral object is about ten times larger compared to the same without accounting for the dust opacities. As such, the prediction of the appearance of the transient would change substantially if dust is included

In addition to the light curve and energy spectrum, polarization is also important for the study of Gamma-ray burst (GRB) prompt emission. Rotation of the polarization angle (PA) with time will cause depolarization of the time-integrated polarization degree. However, it is rarely studied before. Here, we use the magnetic reconnection model with a large-scale ordered aligned magnetic field in the emitting region to study the influence of the key parameters on the PA rotations. We find that half-opening angle of the jet $\theta_{j}$, the observational angle $\theta_{V}$, and the bulk Lorentz factor $\Gamma$ all have significant impacts on the PA rotations. For a fixed $\theta_{j}\Gamma_{0}$ value ($\Gamma_{0}$ is the normalization factor of $\Gamma$), regardless of concrete $\theta_{j}$ and $\Gamma_{0}$ values, PA rotation within $T_{90}$ ($\triangle$PA) remains roughly unchanged for a $q\equiv\theta_{V}/\theta_{j}$ value. As $\theta_{j}\Gamma_{0}$ value increases, the $q$ range for $\triangle$PA$>10^{\circ}$ becomes smaller. The most significant PA rotation with $\triangle$PA$\thicksim90^{\circ}$ will happen when $\theta_{j}\Gamma_{0}\thicksim100$ and $1.1\leq q\leq1.2$. For the top-hat jet, observations of the PA rotation within $T_{90}$ will imply a slightly off-axis observation.

The rotations of the polarization angle (PA) with time (energy) can lead to the depolarization of the time-integrated (energy-integrated) polarization. However, we don't know how and when it will rotate. Here, we consider the magnetic reconnection model to investigate the polarizations, especially the PA rotations of GRB prompt emission. For the large-scale ordered aligned magnetic field configuration, we find that PAs will evolve with time (energy) for off-axis observations. Our studies show that the rotations of the PAs are due to the changes of the ``observed shape'' of the emitting region (before averaged). We apply our models to the single pulse burst of GRB 170101A and GRB 170114A with time-resolved PA observations. We find it can interpret the violent PA variation of GRB 170101A. The model could not predict the twice $90^{\circ}$ PA changes in GRB 170114A. Detailed model should be considered.

The rapid variability of X-ray binaries produces a wide range of X-ray states that are linked to activity across the electromagnetic spectrum. It is particularly challenging to study a sample of sources large enough to include all types in their various states, and to cover the full range of frequencies that show flux density variations. Simultaneous observations with many telescopes are necessary. In this project we monitor 48 X-ray binaries with seven telescopes across the electromagnetic spectrum from 5 x 10^9 Hz to 10^19 Hz, including ground-based radio, IR, and optical observatories and five instruments on two spacecraft over a one-week period. We construct spectral energy distributions and matching X-ray color-intensity diagrams for 20 sources that have the most extensive detections. Our observations are consistent with several models of expected behavior proposed for the different classes: we detect no significant radio emission from pulsars or atoll sources, but we do detect radio emission from Z sources in the normal or horizontal branch, and from black holes in the high/soft, low/hard and quiescent states. The survey data provide useful constraints for more detailed models predicting behavior from the different classes of sources.

The local galaxy peculiar velocity field can be reconstructed from the surrounding distribution of large-scale structure and plays an important role in calibrating cosmic growth and expansion measurements. In this paper, we investigate the effect of the stochasticity of these velocity reconstructions on the statistical and systematic errors in cosmological inferences. By introducing a simple statistical model between the measured and theoretical velocities, whose terms we calibrate from linear theory, we derive the bias in the model velocity. We then use lognormal realisations to explore the potential impact of this bias when using a cosmic flow model to measure the growth rate of structure, and to sharpen expansion rate measurements from host galaxies for gravitational wave standard sirens with electromagnetic counterparts. Although our illustrative study does not contain fully realistic observational effects, we demonstrate that in some scenarios these corrections are significant and result in a measurable improvement in determinations of the Hubble constant compared to standard forecasts.

The correlation between close-in super Earths and distant cold Jupiters in planetary systems has important implications for their formation and evolution. In contrary to some earlier findings, a recent study conducted by Bonomo et al.\ suggests that the occurrence of cold Jupiter companions is not excessive in super Earth systems. Here we show that this discrepancy can be seen as a Simpson's paradox and is resolved once the metallicity dependence of the super Earth--cold Jupiter relation is taken into account. A common feature is noticed that almost all the cold Jupiter detections with inner super Earth companions are found around metal-rich stars. Focusing on the Sun-like hosts with super-solar metallicities, we show that the frequency of cold Jupiters conditioned on the presence of inner super Earths is $39_{-11}^{+12}\%$, whereas the frequency of cold Jupiters in the same metallicity range is no more than $20\%$. Therefore, the occurrences of close-in super Earths and distant cold Jupiters appear correlated around metal-rich hosts. The relation between the two types of planets remains unclear for stars with metal-poor hosts due to the limited sample size and the much lower occurrence rate of cold Jupiters, but a correlation between the two cannot be ruled out.

Gamma~Leo is a long-period visual binary system consisting of K0III (A) and G7III (B) giants, in which particular interest is attracted by the brighter A since the discovery of a planet around it. While detailed spectroscopic comparative study of both components would be worthwhile (e.g., for probing any impact of planet formation on chemical abundances), such a research seems to have been barely attempted as most available studies tend to be biased toward A. Given this situation, the physical properties of A and B along with their differences were investigated based on high-dispersion spectra in order to establish their stellar parameters, evolutionary status, and surface chemical compositions. The following results were obtained. (1) The masses were derived as ~1.7Msun and ~1.6Msun for A and B, respectively, both of which are likely to be in the stage of red clump giants after He-ignition. The mass of the planet around A has also been revised as m*sin(i) = 10.7M_Jupiter (increased by ~20%). (2) These are normal giants of subsolar metallicity ([Fe/H]~-0.4) belonging to the thin-disk population. (3) A as well as B show moderate C deficiency and N enrichment, which are in compatible with the prediction from the standard stellar evolution theory. (4) The chemical abundances of 26 elements are practically the same within <~0.1dex for both components, which implies that the surface chemistry is not appreciably affected by the existence of a planet in A.

Type II radio bursts are often associated with coronal shocks that are typically driven by coronal mass ejections (CMEs) from the Sun. Here, we conduct a case study of a type II radio burst that is associated with a C4.5 class flare and a blowout jet, but without the presence of a CME. The blowout jet is observed near the solar disk center in the extreme-ultraviolet (EUV) passbands with different characteristic temperatures. Its evolution involves an initial phase and an ejection phase with a velocity of 560 km/s. Ahead of the jet front, an EUV wave propagates at a projected velocity of 403 km/s in the initial stage. The moving velocity of the source region of the type II radio burst is estimated to be 641 km/s, which corresponds to the shock velocity against the coronal density gradient. The EUV wave and the type II radio burst are closely related to the ejection of the blowout jet, suggesting that both are likely the manifestation of a coronal shock driven by the ejection of the blowout jet. The type II radio burst likely starts lower than those associated with CMEs. The combination of the velocities of the radio burst and the EUV wave yields a modified shock velocity at 757 km/s. The Alfven Mach number is in the range of 1.09-1.18, implying that the shock velocity is 10%-20% larger than the local Alfven velocity.

We present a chemical and dynamical analysis of the leading arm (LA) and trailing arm (TA) of the Sagittarius (Sgr) stream, as well as for the Sgr dwarf galaxy core (SC), using red giant branch, main sequence, and RR Lyrae stars from large spectroscopic survey data. The different chemical properties among the LA, TA, and SC generally agree with recent studies, and can be understood by radial metallicity gradient established in the progenitor of the Sgr dwarf, followed by preferential stellar stripping from the outer part of the Sgr progenitor. One striking finding is a relatively larger fraction of low-eccentricity stars (e < 0.4) in the LA than in the TA and SC. The TA and SC exhibit very similar distributions. Considering that a tidal tail stripped off from a dwarf galaxy maintains the orbital properties of its progenitor, we expect that the e-distribution of the LA should be similar to that of the TA and SC. Thus, the disparate behavior of the e-distribution of the LA is of particular interest. Following the analysis of Vasiliev et al., we attempt to explain the different e-distribution by introducing a time-dependent perturbation of the Milky Way by the Large Magellanic Cloud (LMC)'s gravitational pull, resulting in substantial evolution of the angular momentum of the LA stars to produce the low-e stars. In addition, we confirm from RR Lyrae stars with high eccentricity (e > 0.6) that the TA stars farther away from the SC are also affected by disturbances from the LMC.

Focal Ratio Degradation (FRD) in fibres is a crucial factor to control in astronomical instruments in order to minimize light loss. As astronomical instrumentation has advanced, the integration of large populations of fibres has become common. However, determining FRD in multiplexed fibre systems has become a challenging and time-consuming task. The Integral Field Unit for the Fiber Arrayed Solar Optical Telescope (FASOT-IFU) represents the most densely arranged fibre-based IFU in a single unit. Due to the close packing of fibres in the V-groove of the slit end, measuring FRD is particularly challenging as the output spots are prone to overlapping with adjacent fibres. In this paper, a novel method based on the quasi-near field model is proposed to enable rapid FRD measurement in highly multiplexed fibre systems like IFUs and multi-object observation systems. The principle and uncertainties associated with the method are investigated. The method's validity is demonstrated by applying it to determine the FRD in FASOT-IFU, with the achieved FRD performance meeting the acceptable requirements of FASOT-IFU, where the output focal ratio primarily falls within the range of 5.0-7.0. The results indicate that the proposed method offers several advantages, including the simultaneous and rapid measurement of FRD in multiple fibres with high accuracy (error smaller than 0.35 in F-ratio). Furthermore, besides FRD, the method exhibits potential for extensive measurements of throughput, scrambling, and spectral analysis.

Here we report an ensemble study of 214 A- and F-type stars observed by \textit{Kepler}, exhibiting the so-called \textit{hump and spike} periodic signal, explained by Rossby modes (r~modes) -- the \textit{hump} -- and magnetic stellar spots or overstable convective (OsC) modes -- the \textit{spike} -- respectively. We determine the power confined in the non-resolved hump features and find additional gravity~modes (g~modes) humps always occurring at higher frequencies than the spike. Furthermore, we derive projected rotational velocities from FIES, SONG and HERMES spectra for 28 stars and the stellar inclination angle for 89 stars. We find a strong correlation between the spike amplitude and the power in the r and g~modes, which suggests that both types of oscillations are mechanically excited by either stellar spots or OsC modes. Our analysis suggests that stars with a higher power in $m=1$ r~modes humps are more likely to also exhibit humps at higher azimuthal orders ($m$ = 2, 3, or 4). Interestingly, all stars that show g~modes humps are hotter and more luminous than the observed red edge of the $\delta$ Scuti instability strip, suggesting that either magnetic fields or convection in the outer layers could play an important role.

Recent independent announcements by several collaborations have shown strong evidence of a Stochastic Gravitational-Wave Background (SGWB) detected through Pulsar Timing Arrays (PTAs). In this study, we investigate the implications of a first-order phase transition occurring within the early universe's dark quantum chromodynamics (dQCD) epoch, specifically within the framework of the mirror twin Higgs dark sector model. Our analysis indicates a distinguishable SGWB signal originating from this phase transition, which can explain the measurements obtained by PTAs. Remarkably, a significant portion of the parameter space within the mirror twin Higgs model that accounts for the SGWB signal also effectively resolves the existing tensions in both the $H_0$ and $S_8$ measurements in Cosmology. This intriguing correlation suggests a possible common origin for these three phenomena. Furthermore, the parameter region, $0.2 < \Delta N_{\rm eff} < 0.5$, where the mirror dark matter component constitutes less than $30\%$ of the total dark matter abundance, can accommodate all current cosmological observations and PTA measurements.

The proposed Daksha mission comprises of a pair of highly sensitive space telescopes for detecting and characterising high-energy transients such as electromagnetic counterparts of gravitational wave events and gamma-ray bursts (GRBs). Along with spectral and timing analysis, Daksha can also undertake polarisation studies of these transients, providing data crucial for understanding the source geometry and physical processes governing high-energy emission. Each Daksha satellite will have 340 pixelated Cadmium Zinc Telluride (CZT) detectors arranged in a quasi-hemispherical configuration without any field-of-view collimation (open detectors). These CZT detectors are good polarimeters in the energy range 100 -- 400 keV, and their ability to measure polarisation has been successfully demonstrated by the Cadmium Zinc Telluride Imager (CZTI) onboard AstroSat. Here we demonstrate the hard X-ray polarisation measurement capabilities of Daksha and estimate the polarisation measurement sensitivity (in terms of the Minimum Detectable Polarisation: MDP) using extensive simulations. We find that Daksha will have MDP of~$30\%$ for a fluence threshold of $10^{-4}$ erg cm$^2$ (in 10 -- 1000 keV). We estimate that with this sensitivity, if GRBs are highly polarised, Daksha can measure the polarisation of about five GRBs per year.

The nearest GRB 170817A provided an opportunity to probe the angular structure of the jet of this short gamma-ray burst (SGRB), by using its off-axis observed afterglow emission. It is investigated that whether the afterglow-constrained jet structures can be consistent with the luminosity of the prompt emission of GRB 170817A. Furthermore, by assuming that all SGRBs including GRB 170817A have the same explosive mechanism and jet structure, we apply the different jet structures into the calculation of the flux and redshfit distributions of the SGRB population, in comparison with the observational distributions of the Swift and Fermi sources. As a result, it is found that the single-Gaussian structure can be basically ruled out, whereas the power-law and two-Gaussian models can in principle survive.

Radio blazars have been linked both to individual high-energy neutrino events and to excesses in likelihood sky maps constructed from lower-energy neutrino data. However, the exact mechanism by which neutrinos are produced in these sources is still unknown. Here, we demonstrate that IceCube neutrinos with energies over 200 TeV, which were previously associated with bright radio blazars, are significantly more likely to be accompanied by flares of lower-energy events, compared to those lacking blazar counterparts. The parsec-scale core radio flux of blazars positioned within the error regions of energetic events is strongly correlated with the likelihood of a coincident day-scale lower-energy neutrino flare reported by IceCube. The probability of a chance correlation is 3.6*10^{-4}. This confirms the neutrino-blazar connection in a new and independent way, and provides valuable clues to understanding the origin of astrophysical neutrinos.

Galaxies in the Universe are distributed in a web-like structure characterised by different large-scale environments: dense clusters, elongated filaments, sheetlike walls, and under-dense regions, called voids. The low density in voids is expected to affect the properties of their galaxies. Indeed, previous studies have shown that galaxies in voids are on average bluer and less massive, and have later morphologies and higher current star formation rates than galaxies in denser large-scale environments. However, it has never been observationally proved that the star formation histories (SFHs) in void galaxies are substantially different from those in filaments, walls, and clusters. Here we show that void galaxies have had, on average, slower SFHs than galaxies in denser large-scale environments. We also find two main SFH types present in all the environments: 'short-timescale' galaxies are not affected by their large-scale environment at early times but only later in their lives; 'long-timescale' galaxies have been continuously affected by their environment and stellar mass. Both types have evolved slower in voids than in filaments, walls, and clusters.

The non-linear relationship between the monochromatic X-ray and UV luminosities in quasars offers the possibility of using high-z quasars as standard candles for cosmological testing. In this paper, we use a high-quality catalog of 1598 quasars extending to redshift 6, to compare the flat and uniformly expanding cosmological model, $R_h$ = ct and $\Lambda$CDM cosmological models which are the most debated. The quasar samples are mainly from the XMM-Newton and the Sloan Digital Sky Survey (SDSS). The final result is that the Akaike Information Criterion favors $\Lambda$CDM over $R_h$=ct with a relative probability of 86.30% versus 13.70%.

Artificial intelligence technology has been widely used in astronomy, and new artificial intelligence technologies and application scenarios are constantly emerging. There have been a large number of papers reviewing the application of artificial intelligence technology in astronomy. However, relevant articles seldom mention telescope intelligence separately, and it is difficult to understand the current development status and research hotspots of telescope intelligence from these papers. This paper combines the development history of artificial intelligence technology and the difficulties of critical technologies of telescopes, comprehensively introduces the development and research hotspots of telescope intelligence, then conducts statistical analysis on various research directions of telescope intelligence and defines the research directions' merits. All kinds of research directions are evaluated, and the research trend of each telescope's intelligence is pointed out. Finally, according to the advantages of artificial intelligence technology and the development trend of telescopes, future research hotspots of telescope intelligence are given.

The stability of astrophysical jets in the linear regime is investigated by presenting the methodology to find the growth rates of the various instabilities. We perturb a cylindrical axisymmetric steady jet, linearize the relativistic ideal magnetohydrodynamic (MHD) equations, and analyze the evolution of the eigenmodes of the perturbation by deriving the differential equations that need to be integrated subject to the appropriate boundary conditions, in order to find the dispersion relation. We also apply the WKBJ approximation and additionally give analytical solutions in some subcases corresponding to unperturbed jets with constant bulk velocity along the symmetry axis.

Missing cascades from TeV blazar beams indicate that collective plasma effects may play a significant role in their energy loss. It is possible to mimic the evolution of such highly energetic pair beams in laboratory experiments using modern accelerators. The fate of the beam is governed by two different processes, energy loss through the unstable mode and energetic broadening of the pair beam through diffusion in momentum space. We chalk out this evolution using a Fokker-Planck approach in which the drift and the diffusion terms respectively describe these phenomena in a compact form. We present particle-in-cell simulations to trace the complete evolution of the unstable beam-plasma system for a generic narrow Gaussian pair beam for which the growth rate is reactive. We show that the instability leads to an energetic broadening of the pair beam, slowing down the instability growth in the linear phase, in line with the analytical and numerical solutions of the Fokker-Planck equation. Whereas in a laboratory experiment the change in the momentum distribution is an easily measured observable as a feedback of the instability, the consequence of diffusive broadening in an astrophysical scenario can be translated to an increase in the opening angle of the pair beam.

We investigate planetary migration in the dead zone of a protoplanetary disk where there are a set of spiral waves propagating inward due to the turbulence in the active zone and the Rossby wave instability (RWI), which occurs at the transition between the dead and active zones. We perform global 3D unstratified magnetohydrodynamical (MHD) simulations of a gaseous disk with the FARGO3D code, using weak gradients in the static resistivity profiles that trigger the formation of a vortex at the outer edge of the dead zone. We find that once the Rossby vortex develops, spiral waves in the dead zone emerge and interact with embedded migrating planets by wave interference, which notably changes their migration. The inward migration becomes faster depending on the mass of the planet, due mostly to the constructive (destructive) interference between the outer (inner) spiral arm of the planet and, the destruction of the dynamics of the horseshoe region by means of the set of background spiral waves propagating inward. The constructive wave interference produces a more negative Lindblad differential torque which inevitably leads to an inward migration. Lastly, for massive planets embedded in the dead zone, we find that the spiral waves can create an asymmetric wider and depeer gap than in the case of $\alpha$-disks, and can prevent the formation of vortices at the outer edge of the gap. The latter could generate a faster or slower migration compared to the standard type-II migration.

Embedded planets are potentially the cause of substructures like gaps and cavities observed in several protoplanetary disks. Thus, the substructures observed in the continuum and in line emission encode information about the presence of planets in the system and how they interact with the natal disk. The pre-transitional disk around the star PDS 70 is the first case of two young planets imaged within a dust depleted gap that was likely carved by themselves. We aim to determine the spatial distribution of the gas and dust components in the PDS 70 disk. The axisymmetric substructures observed in the resulting profiles are interpreted in the context of planet-disk interactions. We develop a thermo-chemical forward model for an axisymmetric disk to explain a subset of the Atacama Large Millimeter/Submillimeter Array (ALMA) band 6 observations of three CO isotopologues plus the continuum towards PDS 70. Combining the inferred gas and dust distributions, the model results in a variable gas-to-dust ratio profile throughout the disk that spans two orders of magnitude within the first $130$ au and shows a step gradient towards the outer disk, which is consistent with the presence of a pressure maxima driven by planet-disk interactions. We find a gas density drop factor of ${\sim} 19$ at the location of the planet PDS 70 c with respect to the peak gas density at $75$ au. Combining this value with literature results on the hydrodynamics of planet-disk interactions, we find this gas gap depth to be consistent with independent planet mass estimates from infrared observations. Our findings point towards gas stirring processes taking place in the common gap due to the gravitational perturbation of both planets.

We use Bayesian inference and nested sampling to develop a non-parametric method to reconstruct the primordial power spectrum $P_{\mathcal{R}}(k)$ from Large Scale Structure (LSS) data. The performance of the method is studied by applying it to simulations of the clustering of two different object catalogues, low-$z$ (ELGs) and high-$z$ (QSOs), and considering two different photometric errors. These object clusterings are derived from different templates of the primordial power spectrum motivated by models of inflation: the Standard Model power law characterized by the two parameters $A_s$ and $n_s$; a local feature template; and a global oscillatory template. Our reconstruction method involves sampling $N$ knots in the log $\{k,P_{\mathcal{R}}(k)\}$ plane. We use two statistical tests to examine the reconstructions for signs of primordial features: a global test comparing the evidences and a novel local test quantifying the power of the hypothesis test between the power law model and the marginalized probability over $N$ model. The method shows good performance in all scenarios considered. In particular, the tests show no feature detection for the SM. The method is able to detect power spectrum deviations at a level of $\approx 2\%$ for all considered features, combining either the low-$z$ or the high-$z$ redshift bins. Other scenarios with different redshift bins, photometric errors, feature amplitudes and detection levels are also discussed. In addition, we include a first application to real data from the Sloan Digital Sky Survey Luminous Red Galaxy Data Release 4 (SDSS LRG 04), finding no preference for deviations from the primordial power law. The method is flexible, model independent, and suitable for its application to existing and future LSS catalogues.

The recent Gaia Data Release 3 has unveiled a catalog of over eight hundred thousand binary systems, providing orbital solutions for half of them. Since most of them are unresolved astrometric binaries, several astrophysical parameters that can be only derived from their relative orbits together with spectroscopic additional data, such as the individual stellar masses, remain unknown, and only the mass of the primary, m1, and a wide interval, [m2_lower, m2_upper], for the secondary companion of main-sequence binaries have been derived. To obtain the correct values for both components, in this paper, we describe an independent analytic algorithm to estimate the two most probable relative orbits and magnitude differences of a certain main-sequence or subgiant astrometric binary using all available Gaia data. Subsequently, both possible solutions are constrained to the one that is consistent with m1, m2_lower and m2_upper. Moreover, we deduce not only the correct values of the individual masses of each binary but also the size of the telescope necessary to resolve their components. The workflow of our algorithm as well as the ESMORGA (Ephemeris, Stellar Masses, and relative ORbits from GAia) catalog with more than one hundred thousand individual masses, spectral types, and effective temperatures derivated from its application are also presented.

Theoretical and observational challenges to standard cosmology such as the cosmological constant problem and tensions between cosmological model parameters inferred from different observations motivate the development and search of new physics. A less radical approach to venturing beyond the standard model is the simple mathematical reformulation of our theoretical frameworks underlying it. While leaving physical measurements unaffected, this can offer a reinterpretation and even solutions of these problems. In this spirit, metric transformations are performed here that cast our Universe into different geometries. Of particular interest thereby is the formulation of cosmology in Minkowski space. Rather than an expansion of space, spatial curvature, and small-scale inhomogeneities and anisotropies, this frame exhibits a variation of mass, length and time scales across spacetime. Alternatively, this may be interpreted as an evolution of fundamental constants. As applications of this reframed cosmological picture, the naturalness of the cosmological constant is reinspected and promising candidates of geometric origin are explored for dark matter, dark energy, inflation and baryogenesis. An immediate observation thereby is the apparent absence of the cosmological constant problem in the Minkowski frame. The formalism is also applied to identify new observable signatures of conformal inhomogeneities, which have been proposed as simultaneous solution of the observational tensions in the Hubble constant, the amplitude of matter fluctuations, and the gravitational lensing amplitude of cosmic microwave background anisotropies. These are found to enhance redshifts to distant galaxy clusters and introduce a mass bias with cluster masses inferred from gravitational lensing exceeding those inferred kinematically or dynamically.

We studied the hub filament system G323.46-0.08 based on archival molecular line data from the SEDIGISM 13CO survey and infrared data from the GLIMPSE, MIPS, and Hi-GAL surveys. G323.46-0.08 consists of three filaments, F-north, F-west, and F-south, that converge toward the central high_mass clump AGAL 323.459-0.079. F-west and Part 1 of the F-south show clear large-scale velocity gradients 0.28 and 0.44 km s-1 pc-1, respectively. They seem to be channeling materials into AGAL 323.459-0.079. The minimum accretion rate was estimated to be 1216 M Myr-1. A characteristic V-shape appears around AGAL 323.459-0.079 in the PV diagram, which traces the accelerated gas motions under gravitational collapse. This has also been supported by model fitting results. All three filaments are supercritical and they have fragmented into many dense clumps. The seesaw patterns near most dense clumps in the PV diagram suggests that mass accretion also occurs along the filament toward the clumps. Our results show that filamentary accretion flows appear to be an important mechanism for supplying the materials necessary to form the central high-mass clump AGAL 323.459-0.079 and to propel the star forming activity taking place therein.

We study the long-term evolution of the Milky Way (MW) over cosmic time by modeling the star formation, cosmic rays, metallicity, stellar dynamics, outflows and inflows of the galactic system to obtain various insights into the galactic evolution. The mass accretion is modeled by the results of cosmological $N$-body simulations for the cold dark matter. We find that the star formation rate is about half the mass accretion rate of the disk, given the consistency between observed Galactic Diffuse X-ray Emissions (GDXEs) and possible conditions driving the Galactic wind. Our model simultaneously reproduces the quantities of star formation rate, cosmic rays, metals, and the rotation curve of the current MW. The most important predictions of the model are that there is an unidentified accretion flow with a possible number density of $\sim10^{-2}$ cm$^{-3}$ and the part of the GDXEs originates from a hot, diffuse plasma which is formed by consuming about 10 % of supernova explosion energy. The latter is the science case for future X-ray missions; XRISM, Athena, and so on. We also discuss further implications of our results for the planet formation and observations of externalgalaxies in terms of the multimessenger astronomy.

We study the production of stochastic gravitational wave background from early dark energy (EDE) model. It is caused by resonant amplification of scalar field fluctuations, which easily takes place for typical EDE potential based on the string axion or $\alpha$-attractor model. The resultant spectrum of gravitational wave background is computed by performing 3D lattice simulations. We show that, specifically in some class of generalized $\alpha$-attractor EDE model, a significant amount of gravitational waves can be produced via tachyonic instability with a peak around femto-Hz frequency range. Models predicting such gravitational waves can be constrained by the cosmic microwave background observations.

Various pulsar timing array (PTA) experiments (NANOGrav, EPTA, PPTA, CPTA, including data from InPTA) very recently reported evidence for excess red common-spectrum signals in their latest datasets, with inter-pulsar correlations following the Hellings-Downs pattern, pointing to a stochastic gravitational wave background (SGWB) origin. Focusing for concreteness on the NANOGrav signal (given that all signals are in good agreement between each other), I inspect whether it supports an inflationary SGWB explanation, finding that such an interpretation calls for an extremely blue tensor spectrum, with spectral index $n_T \simeq 1.8 \pm 0.3$, while Big Bang Nucleosynthesis limits require a very low reheating scale, $T_{\rm rh} \lesssim 10\,{\rm GeV}$. While not impossible, an inflationary origin for the PTA signal is barely tenable: within well-motivated inflationary models it is hard to achieve such a blue tilt, whereas models who do tend to predict sizeable non-Gaussianities, excluded by observations. Intriguingly, ekpyrotic models naturally predict a SGWB with spectral index $n_T=2$, although with an amplitude too suppressed to be able to explain the signal detected by PTA experiments. Finally, I provide explicit expressions for a bivariate Gaussian approximation to the joint posterior distribution for the intrinsic-noise amplitude and spectral index of the NANOGrav signal, which can facilitate extending similar analyses to different theoretical signals.

We introduce a novel approach to boost the efficiency of the importance nested sampling (INS) technique for Bayesian posterior and evidence estimation using deep learning. Unlike rejection-based sampling methods such as vanilla nested sampling (NS) or Markov chain Monte Carlo (MCMC) algorithms, importance sampling techniques can use all likelihood evaluations for posterior and evidence estimation. However, for efficient importance sampling, one needs proposal distributions that closely mimic the posterior distributions. We show how to combine INS with deep learning via neural network regression to accomplish this task. We also introduce NAUTILUS, a reference open-source Python implementation of this technique for Bayesian posterior and evidence estimation. We compare NAUTILUS against popular NS and MCMC packages, including EMCEE, DYNESTY, ULTRANEST and POCOMC, on a variety of challenging synthetic problems and real-world applications in exoplanet detection, galaxy SED fitting and cosmology. In all applications, the sampling efficiency of NAUTILUS is substantially higher than that of all other samplers, often by more than an order of magnitude. Simultaneously, NAUTILUS delivers highly accurate results and needs fewer likelihood evaluations than all other samplers tested. We also show that NAUTILUS has good scaling with the dimensionality of the likelihood and is easily parallelizable to many CPUs.

We present the spectra of Complex Organic Molecules (COMs) detected in HOPS 373SW with the Atacama Large Millimeter/submillimeter Array (ALMA). HOPS 373SW, which is a component of a protostellar binary with a separation of 1500 au, has been discovered as a variable protostar by the JCMT Transient monitoring survey with a modest ~30% brightness increase at submillimeter wavelengths. Our ALMA Target of Opportunity (ToO) observation at ~345 GHz for HOPS 373SW revealed extremely young chemical characteristics with strong deuteration of methanol. The dust continuum opacity is very high toward the source center, obscuring line emission from within 0.03 arcsec. The other binary component, HOPS 373NE, was detected only in C17O in our observation, implying a cold and quiescent environment. We compare the COMs abundances relative to CH3OH in HOPS 373SW with those of V883 Ori, which is an eruptive disk object, as well as other hot corinos, to demonstrate the chemical evolution from envelope to disk. High abundances of singly, doubly, and triply deuterated methanol (CH2DOH, CHD2OH, and CD3OH) and a low CH3CN abundance in HOPS 373SW compared to other hot corinos suggest a very early evolutionary stage of HOPS 373SW in the hot corino phase. Since the COMs detected in HOPS 373SW would have been sublimated very recently from grain surfaces, HOPS 373SW is a promising place to study the surface chemistry of COMs in the cold prestellar phase, before sublimation.

The pattern of neutrino flavor oscillations could be altered by the influence of noisy perturbations such as those arising from a gravitational wave background (GWB). A stochastic process that is consistent with a GWB has been recently reported by the independent analyses of pulsar timing array (PTA) data sets collected over a decadal timescale by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), the European Pulsar Timing Array (EPTA), the Parkes Pulsar Timing Array (PPTA), and the Chinese Pulsar Timing Array (CPTA) collaborations. We investigate the modifications in the neutrino flavor oscillations under the influence of the GWB reported by the PTA collaborations and we discuss how such effects could be potentially revealed in near-future neutrino detectors, possibly helping the discrimination of different models for the GWB below the nHz frequency range.

Two arguments have suggested the presence of subsurface water ice at latitudes lower than 30\textdegree~on Mars. First, the absence of CO2 frost on pole-facing slopes was explained by the presence of subsurface ice. Second, models suggested that subsurface ice could be stable underneath these slopes. We revisit these arguments with a new slope microclimate model. Our model shows that below 30{\deg} latitude, slopes are warmer than previously estimated as the air above is heated by warm surrounding plains. This additional heat prevents the formation of CO2 and subsurface water ice for most slopes. Higher than 30{\deg}S, our model suggests the presence of subsurface water ice. In sparse cases (steep dusty slopes), subsurface ice may exist down to 25{\deg}S. While hypothetical unstable ice deposits cannot be excluded by our model, our results suggest that water ice is rarer than previously thought in the +- 30{\deg} latitude range considered for human exploration.

The NANOGrav Collaboration has recently announced evidence for nHz gravitational waves (GWs), in the form of a Hellings-Downs angular correlation in the common-spectrum process that had been observed previously by them and other Pulsar Timing Arrays (PTAs). We analyze the possibility that these GWs originate from binary supermassive black holes (SMBHs) with total masses $\gtrsim 10^9\, M_{\odot}$. The spectral index of the GW signal differs at 95 % CL from that predicted for binary evolution by GW emission alone, and we find $> 3 \sigma$ evidence that environmental effects such as dynamical friction with gas, stars, and dark matter may be affecting the binary evolution. We estimate the required magnitude and spectrum of such environmental effects and comment on their possible implications for measurements of GWs at higher frequencies.

Here we consider the results of direct measurements of muons in extensive air showers with zenith angles $\theta \le 45^{\circ}$ and energy above $10^{17}$ eV, obtained at the Pierre Auger Observatory and Yakutsk array. In both experiments muons were registered with underground scintillation detectors with $\approx 1.0 \times \sec\theta$ GeV energy threshold. Measured density values were compared to theoretical predictions calculated within the framework of the QGSJet-II.04 hadron interaction model. They differ by factor $1.53 \pm 0.13$(stat). We demonstrate that this difference is due to overestimation of muon densities by 1.22 times and underestimation of primary energy by 1.25 times in the Auger experiment.

M-dwarfs are common stellar hosts of habitable-zone exoplanets. NUV radiation can severely impact the atmospheric and surface conditions of such planets, making characterization of NUV flaring activity a key aspect in determining habitability. We use archival data from the GALEX and XMM-Newton telescopes to study the flaring activity of M-dwarfs in the NUV. The GALEX observations form the most extensive dataset of M-dwarfs in the NUV to date, with exploitation of this data possible due to the new gPhoton2 pipeline. We run a dedicated algorithm to detect flares in the pipeline produced lightcurves and find some of the most energetic flares observed to date within the NUV bandpass, with energies of $\sim 10^{34}$ ergs. Using GALEX data, we constrain flare frequency distributions for stars from M0 to M6 in the NUV up to $10^5$ s in equivalent duration and $10^{34}$ ergs in energy, orders of magnitude above any previous study in the UV. We estimate the combined effect of NUV luminosities and flare rates of stars later than M2 to be sufficient for abiogenesis on habitable zone exoplanets orbiting them. As a counterpoint, we speculate the high frequencies of energetic UV flares and associated coronal mass ejections would inhibit the formation of an ozone layer, possibly preventing genesis of complex Earth-like lifeforms due to sterilizing levels of surface UV radiation. We also provide a framework for future observations of M-dwarfs with ULTRASAT, a wide FoV NUV telescope to be launched in 2026.

JWST has shown that CO2 and CO are common on the surfaces of objects in the Kuiper belt and have apparent surface coverages even higher than that of water ice, though water ice is expected to be significantly more abundant in the bulk composition. Using full Mie scattering theory, we show that the high abundance and the unusual spectral behaviour around the 4.26 micron v1 band of CO2 can be explained by a surface covered in a few micron thick layer of ~ 1-2 micron CO2 particles. CO is unstable at the temperatures in the Kuiper belt, so the CO must be trapped in some more stable species. While hydrate clathrates or amorphous water ice are often invoked as a trapping mechanism for outer solar system ices, the expected spectral shift of the absorption line for a CO hydrate clathrates or trapping in amorphous ice is not seen, nor does the H2O abundance appear to be high enough to explain the depth of the CO absorption line. Instead, we suggest that the CO is created via irradiation of CO2 and trapped in the CO2 grains during this process. The presence of a thin surface layer of CO2 with embedded CO suggests volatile differentiation driving CO2 from the interior as a major process driving the surface appearance of these mid-sized Kuiper belt objects, but the mechanisms that control the small grain size and depth of the surface layer remain unclear.

Axions coupled to nucleons might be copiously emitted from core-collapse supernovae (SNe). If the axion-nucleon coupling is strong enough, axions would be emitted from the SN as a burst and, reaching Earth, may excite the oxygen nuclei in water Cherenkov detectors (${}^{16}{\rm O} + a \to {}^{16}{\rm O}^{*}$). This process will be followed by radiative decay(s) of the excited state (${}^{16}{\rm O}^* \rightarrow {}^{16}{\rm O}+\gamma $) providing a strategy for a direct detection of axions from a Galactic SN in large underground neutrino Cherenkov detectors. Motivated by this possibility, we present an updated calculation of axion-oxygen cross section obtained by using self-consistent continuum Random Phase Approximation. We calculate the branching ratio of the oxygen nucleus de-excitation into gamma-rays, neutrons, protons and $\alpha$-particles. These results are used to revisit the detectability of axions from SN 1987A in Kamiokande-II.

Swift J0549.7-6812 is an Be/X-ray binary system (BeXRB) in the Large Magellanic Cloud (LMC) exhibiting a 6s pulse period. Like many such systems the variable X-ray emission is believed to be driven by the underlying behaviour of the mass donor Be star. In this paper we report on X-ray observations of the brightest known outburst from this system which reached a luminosity of 8 x 10^37 erg/s. These observations are supported by contemporaneous optical photometric observations, the first reported optical spectrum, as well as several years of historical data from OGLE and GAIA. The latter strongly suggest a binary period of 46.1d. All the observational data indicate that Swift J0549.7-6812 is a system that spends the vast majority of its time in X-ray quiescence, or even switched off completely. This suggests that occasional observations may easily miss it, and many similar systems, and thereby underestimate the massive star evolution numbers for the LMC.

If the Universe has non-trivial spatial topology, observables depend on both the parameters of the spatial manifold and the position and orientation of the observer. In infinite Euclidean space, most cosmological observables arise from the amplitudes of Fourier modes of primordial scalar curvature perturbations. Topological boundary conditions replace the full set of Fourier modes with specific linear combinations of selected Fourier modes as the eigenmodes of the scalar Laplacian. We present formulas for eigenmodes in orientable Euclidean manifolds with the topologies $E_1$ - $E_6$, $E_{11}$, $E_{12}$, $E_{16}$, and $E_{18}$ that encompass the full range of manifold parameters and observer positions, generalizing previous treatments. Under the assumption that the amplitudes of primordial scalar curvature eigenmodes are independent random variables, for each topology we obtain the correlation matrices of Fourier-mode amplitudes (of scalar fields linearly related to the scalar curvature) and the correlation matrices of spherical-harmonic coefficients of such fields sampled on a sphere, such as the temperature of the cosmic microwave background (CMB). We evaluate the detectability of these correlations given the cosmic variance of the observed CMB sky. We find that topologies where the distance to our nearest clone is less than about 1.2 times the diameter of the last scattering surface of the CMB give a correlation signal that is larger than cosmic variance noise in the CMB. This implies that if cosmic topology is the explanation of large-angle anomalies in the CMB, then the distance to our nearest clone is not much larger than the diameter of the last scattering surface. We argue that the topological information is likely to be better preserved in three-dimensional data, such as will eventually be available from large-scale structure surveys.

The origin of interstellar and intergalactic magnetic fields is largely unknown, and the primordial magnetic fields (PMFs) produced by, e.g., phase transitions of the early Universe are expected to provide seeds for those magnetic fields. The PMFs affect the evolution of the Universe at an early time, resulting in a series of phenomena. In this work, we show that the PMF-induced turbulence can give rise to nanohertz (nHz) gravitational waves reported by several pulsar timing arrays, including NANOGrav, PPTA, EPTA, and CPTA. Using the nHz gravitational wave data, we obtain the constraints on the characteristic magnetic field strength ($B_{\rm ch}^* \sim \mathcal{O}(1)~\rm{\mu G}$) and coherent length scale ($\ell_{\rm ch}^* \sim \mathcal{O}(1)~\rm{pc}$) of PMFs, assuming a generation temperature of approximately the QCD temperature ($\sim 100$ MeV). In addition, the PMFs which evolve to the recombination era can induce baryon density inhomogeneities, and then alter the ionization process. This naturally results in an alleviation of the tension of the Hubble parameter $H_0$ and the matter clumpiness parameter $S_8$ between early and late-time measurements. Assuming an evolution form of $B_{\rm ch}\sim \ell_{\rm ch}^{-\alpha}$ from the epoch of the production of PMFs to the epoch of recombination, we find $0.91<\alpha<1.08$ (95\% credible region).

Chemically peculiar Ap and Bp stars host strong large-scale magnetic fields in the range of $200$~G up to $30$~kG, which are often considered to be the origin of fossil magnetic fields. We assess the evolution of such fossil fields during the star formation process and the pre-main sequence evolution of intermediate stars, considering fully convective models, models including a transition to a radiative protostar and models with a radiative core. We also examine the implications of the interaction between the fossil field and the core dynamo. We employ analytic and semi-analytic calculations combined with current observational constraints. For fully convective models, we show that magnetic field decay via convection can be expected to be very efficient for realistic parameters of turbulent resistivities. Based on the observed magnetic field strength - density relation, as well as the expected amount of flux loss due to ambipolar diffusion, it appears unlikely that convection could be suppressed via strong enough magnetic fields. On the other hand, a transition from a convective to a radiative core could very naturally explain the survival of a significant amount of flux, along with the presence of a critical mass. We show that in some cases, the interaction of a fossil field with a core dynamo may further lead to changes in the surface magnetic field structure. In the future, it will be important to understand in more detail how the accretion rate evolves as a function of time during the formation of intermediate-mass protostars, including its impact on the protostellar structure. The latter may even allow to derive quantitative predictions concerning the expected population of large scale magnetic fields in radiative stars.

Pulsar Timing Arrays (PTAs) have reported evidence for a stochastic gravitational wave (GW) background at nHz frequencies, possibly originating in the early Universe. We show that the spectral shape of the low-frequency (causality) tail of GW signals sourced at temperatures around $T\gtrsim 1$ GeV is distinctively affected by confinement of strong interactions (QCD), due to the corresponding sharp decrease in the number of relativistic species. A Bayesian analysis in the latest International PTA dataset reveals a significant improvement in the fit with respect to cubic power law spectra, previously employed for the causality tail. Comparison with the results of NANOGrav 15 years and European PTA Data Release 2 suggests that our inclusion of Standard Model effects on GWs can have a potentially decisive impact on model selection.

This research paper focuses on the implementation of radial Basis Function (RBF) Support Vector Machines (SVM) for classifying asteroid orbits. Asteroids are important astronomical objects, and their orbits play a crucial role in understanding the dynamics of the solar system. The International Astronomical Union maintains data archives that provide a playground to experiment with various machine-learning techniques. In this study, we explore the application of RBF SVM algorithm to classify asteroids. The results show that the RBF SVM algorithm provides a good efficiency and accuracy to the dataset. We also analyze the impact of various parameters on the performance of the RBF SVM algorithm and present the optimal parameter settings. Our study highlights the importance of using machine learning techniques for classifying asteroid orbits and the effectiveness of the RBF SVM algorithm in this regard.

Recently, the NANOGrav, PPTA, EPTA and CPTA collaborations reported compelling evidence of the existence of the Stochastic Gravitational-Wave Background (SGWB). The amplitude and spectrum of this inferred gravitational-wave background align closely with the astrophysical predictions for a signal originating from the population of supermassive black-hole binaries. In light of these findings, we explore the possibility to detect dark matter spikes surrounding massive black holes, which could potentially impact the gravitational-wave waveform and modulate the SGWB. We demonstrate that the SMBH binary evolution induced by the combined effects of GW radiation and the dynamical friction of the dark matter spike exhibits detectable manifestations within the nHz frequency range of the SGWB.

We analyze cosmic superstring models in light of NANOGrav 15-year pulsar timing data. A good fit is found for a string tension $G \mu \sim 10^{-12} - 10^{-11}$ and a string intercommutation probability $p \sim 10^{-3} - 10^{-1}$. Extrapolation to higher frequencies assuming standard Big Bang cosmology is compatible at the 68\% CL with the current LIGO/Virgo/KAGRA (LVK) upper limit on a stochastic gravitational wave background (SGWB) in the 10 to 100 Hz range. Most of the superstring parameter space would be accessible to LVK with design parameters, but could be rendered inaccessible by a period of matter-dominated cosmological expansion. However, even in this case a SGWB due to cosmic superstrings would be detectable by ET, AION-km, AEDGE, LISA, the Nancy Roman telescope, GAIA and SKA. A period of inflation could also suppress the superstring SGWB above PTA frequencies, but it would again be detectable by these detectors. We conclude that the superstring interpretation of the NANOGrav data would be robustly testable in these modified cosmological scenarios.

The recent data releases by multiple pulsar timing array experiments (NANOGrav, EPTA, PPTA and CPTA) show evidence for Hellings-Downs angular correlations indicating that the observed stochastic common spectrum can be interpreted as a stochastic gravitational wave background. In this letter, we study whether the signal may originate from gravitational waves induced by high-amplitude primordial curvature perturbations. Such large perturbations may be accompanied by the generation of a sizeable primordial black hole (PBH) abundance. We improve existing analyses of the PBH abundance by including non-Gaussianities typical of several scenarios such as curvaton and inflection-point models. We show that Gaussian scenarios for scalar-induced gravitational waves are disfavoured by more than 2{\sigma} as the sole explanation of the most constraining NANOGrav 15-year data by the overproduction of PBHs. This excludes most explanations relying on single-field inflation by more than 3{\sigma}. This tension, however, can be alleviated in models in which non-Gaussianites suppress the PBH abundance, for instance, in curvaton models with a large rdec or models with a negative fNL. On the flip side, the current NANOGrav data does not constrain the abundance of PBHs in the stellar mass range.

Very recently, the Pulsar Timing Array (PTA) experiments reported strong evidence for the presence of the nanohertz stochastic gravitational wave background (SGWB). In this work we show that the cosmic string loops can account for the nanohertz SGWB data with a $G\mu \sim 2\times 10^{-12}$ and the loops number density $N \sim 10^{4}$. Though the presence of cosmic string loops can also effectively enhance the number density of massive galaxies at high redshifts, we do not find a reasonable parameter space to self-consistently interpret both the SGWB data and the JWST observations. This implies either an extension of the model adopted in this work or the different physical origins of these two phenomena.

Supermassive black hole binaries source gravitational waves measured by Pulsar Timing Arrays. The frequency spectrum of this stochastic background is predicted more precisely than its amplitude. We argue that Dark Matter friction can suppress the spectrum around nHz frequencies, where it is measured, allowing to derive robust and significant bounds on the Dark Matter density, which, in turn, controls indirect detection signals from galactic centers. A precise spectrum of gravitational waves would translate in a tomography of the DM density profile, potentially probing DM particle-physics effects that induce a characteristic DM density profile, such as DM annihilations or de Broglie wavelength.

We report the identification of Lyman Break Galaxy (LBG) candidates around the most luminous Hot Dust-Obscured Galaxy (Hot DOG) known, WISE J224607.56$-$052634.9 (W2246$-$0526) at $z=4.601$, using deep \textit{r}-, \textit{i}-, and \textit{z}-band imaging from the Gemini Multi-Object Spectrograph South (GMOS-S). We use the surface density of LBGs to probe the Mpc-scale environment of W2246$-$0526 to characterize its richness and evolutionary state. We identify LBG candidates in the vicinity of W2246$-$0526 using the selection criteria developed by \cite{2004VOuchi} and \cite{2006Yoshida} in the Subaru Deep Field and in the Subaru XMM-Newton Deep Field, slightly modified to account for the difference between the filters used, and we find 37 and 55 LBG candidates, respectively. Matching to the $z$-band depths of those studies, this corresponds to $\delta = 5.8^{+2.4}_{-1.9}$ times the surface density of LBGs expected in the field. Interestingly, the Hot DOG itself, as well as a confirmed neighbor, do not satisfy either LBG selection criteria, suggesting we may be missing a large number of companion galaxies. Our analysis shows that we are most likely only finding those with higher-than-average IGM optical depth or moderately high dust obscuration. The number density of LBG candidates is not concentrated around W2246$-$0526, suggesting either an early evolutionary stage for the proto-cluster or that the Hot DOG may not be the most massive galaxy, or that the Hot DOG may be affecting the IGM transparency in its vicinity. The overdensity around W2246$-$0526 is comparable to overdensities found around other Hot DOGs and is somewhat higher than typically found for radio galaxies and luminous quasars at a similar redshift.

In this work we study the large-scale structure around a sample of non-fossil systems and compare the results with earlier findings for a sample of genuine fossil systems selected using their magnitude gap. We compute the distance from each system to the closest filament and intersection as obtained from a catalogue of galaxies in the redshift range $0.05 \le z \le 0.7$. We then estimate the average distances and distributions of cumulative distances to filaments and intersections for different bins of magnitude gap. We find that the average distance to filaments is $(3.0\pm 0.8)$ $R_{200}$ for fossil systems, whereas it is $(1.1\pm 0.1)\,R_{200}$ for non-fossil systems. Similarly, the average distance to intersections is larger in fossil than in non-fossil systems, with values of $(16.3\pm 3.2)$ and $(8.9\pm 1.1) \,R_{200}$, respectively. Moreover, the cumulative distributions of distances to intersections are statistically different between fossil and non-fossil systems. Fossil systems selected using the magnitude gap appear to be, on average, more isolated from the cosmic web than non-fossil systems. No dependence is found on the magnitude gap (i.e. non-fossil systems behave in a similar manner independently of their magnitude gap and only fossils are found at larger average distances from the cosmic web). This result supports a formation scenario for fossil systems in which the lack of infalling galaxies from the cosmic web, due to their peculiar position, favours the building of the magnitude gap via the merging of all the massive satellites with the central galaxy. Comparison with numerical simulations suggests that fossil systems selected using the magnitude gap are not old fossils of the ancient Universe, but systems located in regions of the cosmic web not influenced by the presence of intersections.

In a broad class of theories, the accumulation of ultralight dark matter (ULDM) with particles of mass $10^{-22}~\textrm{eV} < m_{\phi} < 1~\textrm{eV}$ leads the to formation of long-lived bound states known as boson stars. When the ULDM exhibits self-interactions, prodigious bursts of energy carried by relativistic bosons are released from collapsing boson stars in bosenova explosions. We extensively explore the potential reach of terrestrial and space-based experiments for detecting transient signatures of emitted relativistic bursts of scalar particles, including ULDM coupled to photons, electrons, and gluons, capturing a wide range of motivated theories. For the scenario of relaxion ULDM, we demonstrate that upcoming experiments and technology such as nuclear clocks as well as space-based interferometers will be able to sensitively probe orders of magnitude in the ULDM coupling-mass parameter space, challenging to study otherwise, by detecting signatures of transient bosenova events. Our analysis can be readily extended to different scenarios of relativistic scalar particle emission.

Atom gradiometers have emerged as compelling broadband probes of scalar ultralight dark matter (ULDM) candidates that oscillate with frequencies between approximately $10^{-2}$ Hz and $10^3$ Hz. ULDM signals with frequencies greater than $\sim 1$ Hz exceed the expected Nyquist frequency of atom gradiometers, and so are affected by aliasing and related phenomena, including signal folding and spectral distortion. To facilitate the discovery of super-Nyquist ULDM signals, in this work we investigate the impact of these effects on parameter reconstruction using a robust likelihood-based framework. We demonstrate that accurate reconstruction of ULDM parameters can be achieved as long as the experimental frequency resolution is larger than the ULDM signal linewidth. Notably, as ULDM candidates whose frequencies differ by integer multiples of the sampling frequency are identified at the same aliased frequency, our discovery analysis recovers discrete islands in parameter space. Our study represents the first comprehensive exploration of aliasing in the context of dark matter direct detection and paves the way for enhanced ULDM detection strategies with atom gradiometers.

The relativistic extension of the classic stellar structure equations is investigated. It is pointed out that the Tolman-Oppenheimer-Volkov (TOV) equation with the gradient equation for gravitational mass can be made complete as a closed set of differential equations by adding that for the Tolman temperature, and the set is proposed as the relativistic hydrostatic structure equations. The exact forms of the relativistic Poisson equation and the steady-state heat conduction equation in the curved spacetime are derived. The application to an ideal gas of particles with the conserved particle number current leads to a strong prediction that the heat capacity ratio almost becomes one in any Newtonian convection zone such as the solar surface. The steady-state heat conduction equation is solved exactly in the system and thermodynamic observables exhibit the power law behavior, which implies the possibility for the system to be a new model of stellar corona and a flaw in the earlier one obtained by using the non-relativistic stellar structure equations. The mixture with another ideal gas yields multilayer structure to a stellar model, in which the classic stellar structure equations are reproduced and analytic multilayer structure of luminous stars is revealed in a suitable approximation.

The paper has explored analogue of gravitational synchrotron massive particle and Penrose process in MOdified Gravity (MOG) known as Scalar-Tensor-Vector-Gravity (STVG). Investigation of the gravitational field around Kerr-MOG black hole showed that it has strong gravitational field with large horizon and can rotate faster than Kerr black hole due to the effect of STVG. We have studied influence of STVG in circular motion of massive particle around Kerr-MOG black hole and discussed the Innermost Stable Circular Orbit (ISCO) of massive test particle. It is shown that STVG plays a crucial role in energy extraction from a rotating black hole, with an energy efficiency of more than $100\%$ according to the Penrose process. Furthermore, we have explored the gravitational synchrotron radiation analogue produced by a massive particle orbiting around a Kerr-MOG black hole. It has been shown that the intensity of gravitational radiation from binary systems of stellar black holes (SBH) and supermassive black holes (SMBH).

The self-interacting dark matter (SIDM) paradigm provides a potential solution to the challenge faced by the cold dark matter model in explaining small-scale structure problems. This paradigm incorporates self-interactions among DM particles, typically mediated by a particle with a mass around MeV. The recent evidences of nano-Hertz gravitational waves from NANOGrav, EPTA, PPTA, and CPTA collaborations indicate a first-order phase transition (FOPT) occurring at a temperature of the MeV scale. Considering the close proximity between these two scales, we postulate that the mediator mass in the SIDM model originates from the spontaneous breaking of a $U(1)'$ symmetry, which is driven by the FOPT indicated by pulsar time array data. Consequently, the alignment of these two scales is believed to be deeply connected by the same underlying physics. Through a comprehensive survey of the parameter space, we identify the viable region favored by SIDM and simultaneously provide an explanation for the pulsar timing array data.

Very recently Pulsar Timing Array (PTA) collaborations have independently reported the evidence for a stochastic gravitational-wave background (SGWB), which can unveil the formation of primordial seeds of inhomogeneities in the early universe. With the SGWB parameters inferred from PTAs data, we can make a prediction of the Primordial Black Hole (PBH) clusters from the domain walls of axion-like particles (ALPs). These primordial seeds can naturally provide a solution to the early Active Galactic Nuclei (AGN) formation indicated by James Webb Space Telescope (JWST). Besides, the mass of ALP is also constrained, $m_a \sim 10^{-15}-10^{-14}$ eV, within the reach of upcoming cavity experiments.

For the typical forbidden dark matter (DM), the correct relic density is determined exclusively by kinetically forbidden DM annihilations which vanish at zero temperature. We present a model that contains the DM and a heavier but unstable scalar mediator in the hidden sector. When the temperature drops below $\sim m_{\rm DM}$, this hidden sector, thermally decoupled from the visible sector, enters a cannibal phase (with zero chemical potential), during which the DM density is depleted with the out-of-equilibrium decay of the scalar mediator. As such, the freeze-out process, described by forbidden DM annihilations to mediators, evolves with a temperature different from the SM bath. The DM candidate of having a mass in the range of tens of MeV can result in the correct relic density and sizable 2-to-2 self-interactions, which fit small structure problems. The future sensitivity of the NA62 beam dump experiment can probe the parameter space of the related scalar mediator.

We investigate the potential of the warped-extradimension framework as an explanation for the recently observed stochastic gravitational background at nHz frequencies in pulsar timing arrays (PTA). Our analysis reveals that the PTA data can be effectively accommodated by a first-order phase transition triggered by a radion at the MeV-GeV scale feebly coupled to the Standard Model. Remarkably, this outcome remains robust irrespective of the specific details of the warped extradimension embedding, providing a foundation for future investigations aiming to develop concrete extradimension descriptions of Nature. We also demonstrate that many existing embeddings are not viable, as their radion and graviton phenomenology clash with a MeV-GeV scale radion. As a possible way-out, we sketch a promising solution involving multiple branes, wherein the light radion, graviton, and ensuing light resonances remain consistent with collider bounds and gravity tests.

We explore the possibility that a confining first-order phase transition of a nearly-conformal dark sector generates the reported NANOGrav signal of a stochastic gravitational wave background. The visible Standard Model (SM) sector and the dark sector are initially thermally decoupled so that their temperatures are different. The nearly conformal phase transition is described by the shallow potential of a dilaton (or a radion in the 5D holographic perspective) generated by a new dark Yang-Mills field coupled to the conformal sector. For a dark sector only gravitationally connected with the visible sector, the NANOGrav signal is explained by the phase transition without contradicting the $\Delta N_{\rm eff}$ constraint, together with a contribution from supermassive black hole binaries. While the dilaton and dark glueballs can be produced after the phase transition, they immediately decay into dark radiation, which can help ameliorate the Hubble tension and be tested by the future CMB-S4 experiment. Alternatively, for a dark conformal sector decaying into the visible sector after the phase transition, the $\Delta N_{\rm eff}$ constraint is not applied and the phase transition can solely explain the NANOGrav signal.

Shock waves are common in astrophysical environments. On many occasions, they are collisionless, which means they occur in settings where the mean free path is much larger than the dimensions of the system. For this very reason, magnetohydrodynamic (MHD) is not equipped to deal with such shocks, be it because it assumes binary collisions, hence temperature isotropy, when such isotropy is not guaranteed in the absence of collisions. Here we solve a model capable of dealing with perpendicular shocks with anisotropic upstream pressure. The system of MHD conservation equations is closed assuming the temperature normal to the flow is conserved at the crossing of the shock front. In the strong shock sonic limit, the behavior of a perpendicular shock with isotropic upstream is retrieved, regardless of the upstream anisotropy. Generally speaking, a rich variety of behaviors is found, inaccessible to MHD, depending on the upstream parameters. The present work can be viewed as the companion paper of MNRAS 520, 6083-6090 (2023), where the case of a parallel shock was treated. Differences and similarities with the present case are discussed.

Recently, the Hellings Downs correlation has been observed by different pulsar timing array (PTA) collaborations, such as NANOGrav, European PTA, Parkes PTA, and Chinese PTA. These experimental studies through PTA of the most precise pulsars within the Milky Way show the first robust evidence for the stochastic gravitational wave background of our Universe. We study the ultralight axion interpretation of the new discovery by investigating the gravitational wave from the energy-level transition of the gravitational atoms, which is composed of cosmic populated Kerr black holes and their surrounding axion clouds from the superradiance process. We demonstrate that this new observation admits an axion interpretation for the ultralight axion mass in the range $10^{-21}\sim 10^{-20}$~eV.

Transient gravitational waves (aka gravitational wave bursts) within the nanohertz frequency band could be generated by a variety of astrophysical phenomena such as the encounter of supermassive black holes, the kinks or cusps in cosmic strings, or other as-yet-unknown physical processes. Radio-pulses emitted from millisecond pulsars could be perturbed by passing gravitational waves, hence the correlation of the perturbations in a pulsar timing array can be used to detect and characterize burst signals with a duration of $\mathcal{O}(1\text{-}10)$ years. We propose a fully Bayesian framework for the analysis of the pulsar timing array data, where the burst waveform is generically modeled by piecewise straight lines, and the waveform parameters in the likelihood can be integrated out analytically. As a result, with merely three parameters (in addition to those describing the pulsars' intrinsic and background noise), one is able to efficiently search for the existence and the sky location of {a burst signal}. If a signal is present, the posterior of the waveform can be found without further Bayesian inference. We demonstrate this model by analyzing simulated data sets containing a stochastic gravitational wave background {and a burst signal generated by the parabolic encounter of two supermassive black holes.

We show that the recently reported NANOGrav, EPTA, PPTA, and CPTA data suggesting the existence of stochastic gravitational waves in the nanohertz region can be explained by axion domain walls coupled to QCD. In this scenario, the non-perturbative effects of QCD generate a temperature-dependent bias for the domain wall around the QCD phase transition, leading to an immediate collapse of the domain walls. We perform dedicated lattice simulations of the axion domain walls, taking into account the temperature dependence of the bias, to estimate the gravitational waves emitted during the domain wall annihilation process. We also discuss the future prospects for accelerator-based searches for the axion and the potential for the formation and detection of primordial black holes.

For a discrete symmetry that is anomalous under QCD, the domain walls produced in the early universe from its spontaneous breaking can naturally annihilate due to QCD instanton effects. The gravitational waves generated from wall annihilation have their amplitude and frequency determined by both the discrete symmetry breaking scale and the QCD scale. The evidence of stochastic gravitational waves at nanohertz observed by pulsar timing array experiments suggests that the discrete-symmetry-breaking scale is around 100 TeV, assuming the domain-wall explanation. The annihilation temperature is about 100 MeV, which could naturally be below the QCD phase transition temperature. We point out that the QCD phase transition within some domains with an effective large QCD $\theta$ angle could be a first-order one. To derive the phase diagram in $\theta$ and temperature, we adopt a phenomenological linear sigma model with three quark flavors. The domain-wall explanation for the NANOGrav, EPTA, PPTA and CPTA results hints at a first-order QCD phase transition, which predicts additional gravitational waves at higher frequencies. If the initial formation of domain walls is also a first-order process, this class of domain-wall models predicts an interesting gravitational wave spectroscopy with frequencies spanning more than ten orders of magnitude, from nanohertz to 100 Hz.