This is the list of the papers for the past 5 days that include local authors affiliated with Princeton University's Astrophysical Sciences department.
The swift assembly of the earliest galaxies poses a significant challenge to our understanding of galaxy formation. Ultra-massive quiescent galaxies at intermediate redshifts ($3 < z < 5$) currently present one of the most pressing problems for theoretical modeling, since very few mechanisms can be invoked to explain how such galaxies formed so early in the history of the Universe. Here, we exploit the fact that these galaxies all reside within significant overdensities to explain their masses. To this end, we construct and release a modified version of the Extreme Value Statistics (EVS) code which takes into account galaxy environment by incorporating clustering in the calculation. With this new version of EVS, we find that ultra-massive quiescent galaxies at $3<z<5$ do not present as serious a tension with simple models of galaxy formation when the analysis of a given galaxy is conditioned on its environment.
A fraction of active galactic nuclei (AGN) have double-peaked H$\alpha$, H$\beta$ and Mg II broad lines attributed to emission from rotating gas in the accretion disk. Using optical spectroscopy of a flux-limited sample of AGN selected via ultrahard X-rays from the BAT AGN Spectroscopic Survey (BASS), we systematically identify 71 double-peaked emitters amongst 343 broad-line AGN with redshifts $0.004<z<0.297$ and 2-10 KeV X-ray luminosities of log 40-45.7 (erg/s), and provide their best-fit accretion disk geometry parameters. We find that ~21% of X-ray selected broad-line AGN are double-peaked emitters (DPEs), consistent with rates previously reported for $z<0.2$ broad-line AGN selected for strong optical variability in ZTF. 11 of 71 DPEs (15%) exhibited a single-peaked Gaussian component to the broad line profile in addition to the double-peaked disk profile. In this sample, DPEs have intrinsically higher masses by ~0.4 dex and lower Eddington ratios by ~0.3 dex than other broad-line AGN, and have a preference for elliptical host galaxies, higher X-ray luminosities, and higher [OI] $\lambda$6302 to narrow H$\alpha$ flux ratios than other broad-line AGN. We find that DPEs are not segregated from other broad-line AGN in the $L_{\rm bol}$ vs $M_{\rm BH}$ relation or their X-ray to radio luminosity ratios, and do not show a preference for intermediate Seyfert types over Seyfert 1s. We do not find differences in a wide range of multi-wavelength properties when comparing DPEs to other broad-line AGN, including optical and mid-IR variability levels, and the rate of changing-look events. We discuss the two populations in the context of multi-component disk-wind models of the AGN broad line region and consider how unrecognized contributions of disk emission to the broad lines introduce biases in virial SMBH mass estimates.
Cosmological simulations provide a wealth of data in the form of point clouds and directed trees. A crucial goal is to extract insights from this data that shed light on the nature and composition of the Universe. In this paper we introduce CosmoBench, a benchmark dataset curated from state-of-the-art cosmological simulations whose runs required more than 41 million core-hours and generated over two petabytes of data. CosmoBench is the largest dataset of its kind: it contains 34 thousand point clouds from simulations of dark matter halos and galaxies at three different length scales, as well as 25 thousand directed trees that record the formation history of halos on two different time scales. The data in CosmoBench can be used for multiple tasks -- to predict cosmological parameters from point clouds and merger trees, to predict the velocities of individual halos and galaxies from their collective positions, and to reconstruct merger trees on finer time scales from those on coarser time scales. We provide several baselines on these tasks, some based on established approaches from cosmological modeling and others rooted in machine learning. For the latter, we study different approaches -- from simple linear models that are minimally constrained by symmetries to much larger and more computationally-demanding models in deep learning, such as graph neural networks. We find that least-squares fits with a handful of invariant features sometimes outperform deep architectures with many more parameters and far longer training times. Still there remains tremendous potential to improve these baselines by combining machine learning and cosmology to fully exploit the data. CosmoBench sets the stage for bridging cosmology and geometric deep learning at scale. We invite the community to push the frontier of scientific discovery by engaging with this dataset, available at this https URL
Understanding how the dynamical state of the interstellar medium (ISM) changes across spatial scales can provide important insights into how the gas is organized and ultimately collapses to form stars. To this end, we present ALMA $^{12}\mathrm{CO}(2-1)$ observations at $7$ pc ($0.4''$) spatial resolution across a $1.4~\mathrm{kpc}\times5.6~\mathrm{kpc}$ ($1'.3\times1'.3$) region located in the disk of the nearby ($D = 3.5$ Mpc), massive, star-forming galaxy NGC 253. We decompose this emission with a hierarchical, multiscale dendrogram algorithm to identify 2463 structures with deconvolved sizes ranging from $\sim3$ to $300$ pc, complete to a limiting mass of $10^4~\mathrm{M_\odot}$. By comparing the virial parameter of these structures against physical properties including size, mass, surface density, velocity dispersion, and hierarchical position, we carry out a comprehensive search for a preferred scale at which gravitationally bound structures emerge. Ultimately, we do not identify any emergent scale for bound objects in our data, nor do we find a correlation between the virial parameter and structure sizes. These findings suggest that gravitational binding cannot be used to define molecular clouds and emphasize the need for multiscale approaches to characterize the ISM.
The China Space Station Telescope (CSST) is a next-generation Stage-IV sky survey telescope, distinguished by its large field of view (FoV), high image quality, and multi-band observation capabilities. It can simultaneously conduct precise measurements of the Universe by performing multi-color photometric imaging and slitless spectroscopic surveys. The CSST is equipped with five scientific instruments, i.e. Multi-band Imaging and Slitless Spectroscopy Survey Camera (SC), Multi-Channel Imager (MCI), Integral Field Spectrograph (IFS), Cool Planet Imaging Coronagraph (CPI-C), and THz Spectrometer (TS). Using these instruments, the CSST is expected to make significant contributions and discoveries across various astronomical fields, including cosmology, galaxy and active galactic nuclei (AGN), the Milky Way and nearby galaxies, stars, exoplanets, Solar System objects, astrometry, and transients and variable sources. This review aims to provide a comprehensive overview of the CSST instruments, observational capabilities, data products, and scientific potential.
Juno and Cassini have shown that Jupiter and Saturn likely contain extended gradients of heavy elements. Yet, how these gradients can survive over billions of years remains an open question. Classical convection theories predict rapid mixing and homogenization, which would erase such gradients on timescales far shorter than the planets' ages. To address this, we estimate the energy required to erode both dense and fuzzy cores, and compare it to what the planet can realistically supply. If the entire cooling budget is available to drive mixing, then even a compact core can, in principle, be destroyed. But if mixing is limited to the thermal energy near the core, which is another plausible scenario, the energy falls short. In that case, Jupiter can erode a fuzzy core by up to approximately $10~\mearth$, but a compact one remains intact. Saturn's core is more robust. Even in the fuzzy case, only about $1~\mearth$ is lost, and if the core is compact, erosion is negligible. The outcome depends sensitively on the assumed initial temperature and entropy profiles. Hotter and more superadiabatic interiors are more prone to mixing. We suggest that 3D simulations of convection driven from above, with realistic stratification and enough depth (i.e., many density scale heights) would be of great interest to further constrain the energy budget for core erosion.
Hot Jupiters (HJs) are commonly thought to host the strongest dynamo-generated magnetic fields among exoplanets, up to one order of magnitude larger than Jupiter. Thus, they have often been regarded as the most promising exoplanets to display magnetic star-planet interaction signals and magnetically-driven coherent radio emission, which unfortunately remains elusive, despite many diversified observational campaigns. In this work, we investigate the evolution of the internal convection and dynamo properties of HJs via one-dimensional models. We explore the dependency on orbital distance, planetary and stellar masses, and types of heat injection. We employ one-dimensional evolutionary models to obtain internal convective structures. Specifically, we obtain the Rossby number $\mathrm{Ro}$ as a function of planetary depth and orbital period, after showing that tidal synchronization is likely valid for all HJs. When the heat is applied uniformly, the convective layers of almost all HJs remain in the fast rotator regime, $\mathrm{Ro} \lesssim 0.1$, except possibly the most massive planets with large orbital distances (but still tidally locked). We recover magnetic field strengths for inflated HJs by applying well-known scaling laws for fast rotators. When strong heat sources are applied mostly in the outer envelope and outside the dynamo region, as realistic Ohmic models predict, convection in the dynamo region often breaks down. Consequently, the heat flux and the derived surface magnetic fields can be greatly reduced to or below Jovian values, contrary to what is commonly assumed, thus negatively affecting estimates for coherent radio emission, and possibly explaining the failure in detecting it so far.
In this study, we report multi-year polarization measurements of four repeating FRBs initially discovered by CHIME: FRBs~20190117A, 20190208A, 20190303A, and 20190417A. We observed the four repeating FRBs with FAST, detecting a total of 66 bursts. Two bursts from FRB~20190417A exhibit a circular polarization signal-to-noise ratio greater than 7, with the highest circular polarization fraction recorded at 35.7%. While the bursts from FRBs 20190208A and 20190303A are highly linearly polarized, those from FRBs~20190117A and 20190417A show depolarization due to multi-path propagation, with \sigma_{\mathrm{RM}} = 2.78 \pm 0.05 rad m$^{-2}$ and 5.19 \pm 0.09 rad m$^{-2}$, respectively. The linear polarization distributions among five repeating FRB--FRBs~20190208A, 20190303A, 20201124A, 20220912A, and 20240114A--are nearly identical but show distinct differences from those of non-repeating FRBs. FRBs~20190117A, 20190303A, and 20190417A exhibit substantial rotation measure (RM) variations between bursts, joining other repeating FRBs in this behavior. Combining these findings with published results, 64% of repeating FRBs show RM variations greater than 50 rad m$^{-2}$, and 21\% exhibit RM reversals. A significant proportion of repeating FRBs reside in a dynamic magneto-ionic environment. The structure function of RM variations shows a power-law index of $\gamma \sim (0-0.8)$, corresponding to a shallow power spectrum $\alpha = -(\gamma + 2) \sim -(2.0-2.8)$ of turbulence, if the RM variations are attributed to turbulence. This suggests that the variations are dominated by small-scale RM density fluctuations. We perform K-S tests comparing the RMs of repeating and non-repeating FRBs, which reveal a marginal dichotomy in the distribution of their this http URL caution that the observed dichotomy may be due to the small sample size and selection biases.
In this study, we report multi-year polarization measurements of four repeating FRBs initially discovered by CHIME: FRBs~20190117A, 20190208A, 20190303A, and 20190417A. We observed the four repeating FRBs with FAST, detecting a total of 66 bursts. Two bursts from FRB~20190417A exhibit a circular polarization signal-to-noise ratio greater than 7, with the highest circular polarization fraction recorded at 35.7%. While the bursts from FRBs 20190208A and 20190303A are highly linearly polarized, those from FRBs~20190117A and 20190417A show depolarization due to multi-path propagation, with \sigma_{\mathrm{RM}} = 2.78 \pm 0.05 rad m$^{-2}$ and 5.19 \pm 0.09 rad m$^{-2}$, respectively. The linear polarization distributions among five repeating FRB--FRBs~20190208A, 20190303A, 20201124A, 20220912A, and 20240114A--are nearly identical but show distinct differences from those of non-repeating FRBs. FRBs~20190117A, 20190303A, and 20190417A exhibit substantial rotation measure (RM) variations between bursts, joining other repeating FRBs in this behavior. Combining these findings with published results, 64% of repeating FRBs show RM variations greater than 50 rad m$^{-2}$, and 21\% exhibit RM reversals. A significant proportion of repeating FRBs reside in a dynamic magneto-ionic environment. The structure function of RM variations shows a power-law index of $\gamma \sim (0-0.8)$, corresponding to a shallow power spectrum $\alpha = -(\gamma + 2) \sim -(2.0-2.8)$ of turbulence, if the RM variations are attributed to turbulence. This suggests that the variations are dominated by small-scale RM density fluctuations. We perform K-S tests comparing the RMs of repeating and non-repeating FRBs, which reveal a marginal dichotomy in the distribution of their this http URL caution that the observed dichotomy may be due to the small sample size and selection biases.
Intermediate mass ratio inspirals (IMRIs) of binary black holes with mass ratios $10^{-4}\lesssim q \lesssim 0.1$ are astrophysically interesting sources of gravitational waves. Mergers of intermediate-mass black holes (IMBHs) with stellar-mass black holes would be IMRIs, so their detection can help us probe the formation mechanisms of IMBHs. They can also help us perform precise tests of general relativity due to the presence of strong higher-order mode emission. We perform a search for aligned-spin IMRIs within the data of the two LIGO detectors in the third observing run (O3) of the LIGO-Virgo-KAGRA (LVK) collaboration, including higher modes in the template banks for the first time. We use the IAS-HM pipeline for our search and construct template banks in the range $1/100 < q<1/18$ using the SEOBNRv5HM waveform model. Our banks retain a similar level of effectualness for IMRPhenomXHM and BHPTNRSur2dq1e3 waveforms, making our search results relatively robust against waveform systematics. We show that the sensitivity volume of the search increases by up to $\sim 500\%$ upon inclusion of higher modes. We do not find any significant candidates with inverse false alarm rate (IFAR) $> 1$ year in the O3 data. This gives us upper limits on the IMRI merger rate in the local Universe, ranging from $\sim 30$ to $10^3$ Gpc$^{-3}$ yr$^{-1}$ depending on the masses of the black holes in the binary. These constraints are consistent with rate predictions in the literature. Our projections indicate that we would be able to detect IMRIs or constrain some of their proposed formation channels in the fourth (O4) and fifth (O5) observing runs.
this https URL . A Frozen Version of the Codes with Produced Data are Also Available on Zenodo: this https URL
One striking feature of binary black hole (BBH) mergers observed in the first decade of gravitational-wave astronomy is an excess of events with component masses around $35\,\mathrm{M}_{\odot}$. Multiple formation channels have been proposed to explain this excess. To distinguish among these channels, it is essential to examine their predicted population-level distributions across additional parameters. In this work, we focus on BBH mergers near the $35\,\mathrm{M}_{\odot}$ peak and infer the population distributions of primary mass ($m_1$), mass ratio ($q$), effective spin ($\chi_{\rm eff}$), and redshift ($z$). We observe a gradual increase in the merger rate with $m_1$, rising by a factor of $3$ from $20\,\mathrm{M}_{\odot}$ to a peak around $34\,\mathrm{M}_{\odot}$, followed by a sharp, order-of-magnitude decline by $50\,\mathrm{M}_{\odot}$. This population also shows a weak preference for equal-mass mergers and has a $\chi_{\rm eff}$ distribution skewed toward positive values, with a median of zero excluded at approximately $90\%$ confidence. We find no significant $q-\chi_{\rm eff}$ correlation in the $35\, \mathrm{M}_{\odot}$ peak population, suggesting that lower-mass systems ($m_1<20\,\mathrm{M}_{\odot}$) likely drive the $q-\chi_{\rm eff}$ anti-correlation observed in the full BBH merger catalog. The redshift evolution of the merger rate is consistent with the cosmic star formation rate. We compare our findings with predictions from a wide range of formation channels. We find that common variants of the pair-instability supernova scenario, as well as hierarchical merger models, are incompatible with the observed features of the $35\,\mathrm{M}_{\odot}$ population. Ultimately, none of the formation channels we consider can explain all or even most of the features observed in this population. The ''mid-thirties'' of black hole mergers are in crisis.
We explore the effects of the neutrino collisional flavor instability (CFI) based on 1D and 2D core-collapse supernova (CCSN) simulations done using the sophisticated radiation-hydrodynamic code F{\sc ornax}. We compare the growth rates of homogeneous CFI {(hCFI)} modes calculated by numerically solving the multi-group dispersion relation to those calculated using the monochromatic approximation. We find that the widely-used monochromatic approximation {leads to incorrect growth rates} when applied in multi-group scenarios. {As opposed to the $\sim10^5$ s$^{-1}$ values given by the monochromatic approximation,} the actual growth rates of non-resonance multi-group {hCFI} are at most $\sim$200 s$^{-1}$ in all our models and they are too slow to affect CCSN outcomes. We adopt a BGK flavor conversion scheme in the simulations to include the effects of resonance-like {hCFI}. We find that the CCSN dynamics and neutrino emission properties are only weakly influenced, and the intrinsic stochasticity due to convection and neutrino-driven turbulence can naturally lead to comparable effects. Hence, our analysis of the non-resonance and resonance-like {hCFI} into CCSN simulations suggests that the effects of neutrino flavor conversion triggered by {hCFI} modes are in general small.
Changing-Look Active Galactic Nuclei (CLAGN) are a rare subset of AGN that show significant changes to the flux of broad Balmer emission lines. Recent studies of CLAGN, such as 1ES 1927+654 and Mrk 590, have revealed that changes in the optically observed accretion rate are accompanied by changes in radio activity. We present a time-domain population study of 474 spectroscopically confirmed CLAGN at radio wavelengths using the Australia SKA Pathfinder Variable and Slow Transients Survey and the Very Large Array Sky Survey. We compare the radio properties of this CLAGN sample to a control sample of AGN that have not had recent changing-look events, and to AGN that were found to have transitioned from radio-quiet to radio-loud over 10-year timescales in VLASS. For 20 newly studied CLAGN detected in ASKAP VAST, we do not detect Mrk 590 or 1ES 1927+654-like fading of the radio flux in the 10 years following changing-look events. For 6 CLAGN with a sufficiently low redshift and high enough mass, we rule out a Mrk 590-like flare. We find that at the population level, CLAGN have higher VAST/VLASS detection rates, lower fractions of radio loudness, and higher variability rates in the 1 GHz frequency compared to the control AGN. Through VLA observations of radio SEDs and Magellan spectroscopic observations, we do not find evidence of a link between CLAGN and AGN that transitioned from radio-loud to radio-quiet in VLASS. We discuss the implications of this study for the physical mechanisms that drive enhanced accretion episodes.
Light reprocessed by dust grains emitting in the infrared allows the study of the physics at play in dusty, embedded regions, where ultraviolet and optical wavelengths are attenuated. Infrared telescopes such as JWST have made it possible to study the earliest feedback phases, when stars are shielded by cocoons of gas and dust. This phase is crucial for unravelling the effects of feedback from young stars, leading to their emergence and the dispersal of their host molecular clouds. Here we show that the transition from the embedded to the exposed phase of star formation is short (< 4 Myr) and sometimes almost absent (< 1 Myr), across a sample of 37 nearby star-forming galaxies, covering a wide range of morphologies from massive barred spirals to irregular dwarfs. The short duration of the dust-clearing timescales suggests a predominant role of pre-supernova feedback mechanisms in revealing newborn stars, confirming previous results on smaller samples and allowing, for the first time, a statistical analysis of their dependencies. We find that the timescales associated with mid-infrared emission at 21 {\mu}m, tracing a dust-embedded feedback phase, are controlled by a complex interplay between giant molecular cloud properties (masses and velocity dispersions) and galaxy morphology. We report relatively longer durations of the embedded phase of star formation in barred spiral galaxies, while this phase is significantly reduced in low-mass irregular dwarf galaxies. We discuss tentative trends with gas-phase metallicity, which may favor faster cloud dispersal at low metallicities.
We present the discovery of 30 transiting giant planets that were initially detected using data from NASA's Transiting Exoplanet Survey Satellite (TESS) mission. These new planets orbit relatively bright ($G \leq 12.5$) FGK host stars with orbital periods between 1.6 and 8.2 days, and have radii between 0.9 and 1.7 Jupiter radii. We performed follow-up ground-based photometry, high angular-resolution imaging, high-resolution spectroscopy and radial velocity monitoring for each of these objects to confirm that they are planets and determine their masses and other system parameters. The planets' masses span more than an order of magnitude ($0.17\,M_J < M_p < 3.3\,M_J$). For two planets, TOI-3593 b and TOI-4961 b, we measured significant non-zero eccentricities of $0.11^{+0.05}_{-0.03}$ and $0.18^{+0.04}_{-0.05}$ respectively, while for the other planets, the data typically provide a 1-$\sigma$ upper bound of 0.15 on the eccentricity. These discoveries represent a major step toward assembling a complete, magnitude-limited sample of transiting hot Jupiters around FGK stars.
We present results on extending the strong lens discovery space down to much smaller Einstein radii ($\theta_E\lesssim0.03''$) and much lower halo mass ($M_\mathrm{halo}<10^{11}M_\odot$) through the combination of JWST observations and machine learning (ML) techniques. First, we forecast detectable strong lenses with JWST using CosmoDC2 as the lens catalog, and a source catalog down to 29th magnitude. By further incorporating the VELA hydrodynamical simulations of high-redshift galaxies, we simulate strong lenses. We train a ResNet on these images, achieving near-100\% completeness and purity for ``conventional" strong lenses ($\theta_E\gtrsim 0.5''$), applicable to JWST, HST, the Roman Space Telescope and Euclid VIS. For the first time, we also search for very low halo mass strong lenses ($M_{halo}<10^{11}M_\odot$) in simulations, with $\theta_E\ll 0.5''$, down to the best resolution ($0.03''$) and depth (10,000~sec) limits of JWST using ResNet. A U-Net model is employed to pinpoint these small lenses in images, which are otherwise virtually impossible for human detection. Our results indicate that JWST can find $\sim 17$/deg$^2$ such low-halo-mass lenses, with the locations of $\sim 1.1$/deg$^2$ of these detectable by the U-Net at $\sim100$\% precision (and $\sim 7.0$/deg$^2$ at a 99.0\% precision). To validate our model for finding ``conventional" strong lenses, we apply it to HST images, discovering two new strong lens candidates previously missed by human classifiers in a crowdsourcing project (Garvin et al. 2022). This study demonstrates the (potentially ``superhuman") advantages of ML combined with current and future space telescopes for detecting conventional, and especially, low-halo-mass strong lenses, which are critical for testing CDM models.