Locally authored papers of the past 5 days

The OSIRIS-REx sample return capsule hypersonic re-entry into the atmosphere is a rare opportunity to test a variety of sonic boom source models since the projectile dimensions are well characterized. While the as-flown flight path is unknown, the predicted flight path enables a rough approximation of the source Mach number and location. Six infrasound microphones deployed in the boom carpet along the predicted flight path recorded impulsive signals from the OSIRIS-REx re-entry. Using a suite of atmosphere profiles and the geometric acoustics approximation, we estimate locations with uncertainty estimates along the flight path from which the signals were emitted. Acoustic overpressure and signal duration predictions from Whitham's far field theory, Carlson's simplified sonic boom prediction method, and a drag-dominated hypersonic model are analyzed with uncertainty estimates from the location estimate. While the Carlson simplified sonic boom prediction method could be accurate, our preference is for the drag-dominated source model. Using this source model with an inviscid Burgers' equation solver for propagation, we obtained an excellent match to the recorded data. These results will help better inform future sample return capsule re-entry observation campaigns as well as contribute to a better understanding of high altitude infrasonic sources.

Astrophysical searches for dark matter in the Milky Way require a reliable model for its density distribution, which in turn depends on the influence of baryonic feedback on the Galaxy. In this work, we utilize a new suite of Milky Way-mass halos from the DREAMS Project, simulated with Cold Dark Matter (CDM),to quantify the influence of baryon feedback and intrinsic halo-to-halo variance on dark matter density profiles. Our suite of 1024 halos varies over supernova and black hole feedback parameters from the IllustrisTNG model, as well as variations in two cosmological parameters. We find that Milky Way-mass dark matter density profiles in the IllustrisTNG model are largely insensitive to astrophysics and cosmology variations, with the dominant source of scatter instead arising from halo-to-halo variance. However, most of the (comparatively minor) feedback-driven variations come from the changes to supernova prescriptions. By comparing to dark matter-only simulations, we find that the strongest supernova wind energies are so effective at preventing galaxy formation that the halos are nearly entirely collisionless dark matter. Finally, regardless of physics variation, all the DREAMS halos are roughly consistent with a halo contracting adiabatically from the presence of baryons, unlike models that have bursty stellar feedback. This work represents a step toward assessing the robustness of Milky Way dark matter profiles, with direct implications for dark matter searches where systematic uncertainty in the density profile remains a major challenge.

The correlation between galaxy stellar mass and gas-phase metallicity, known as the mass-metallicity relation (MZR), gives key insights into the processes that govern galaxy evolution. However, unquantified observational and selection biases can result in systematic errors in attempts to recover the intrinsic MZR, particularly at higher redshifts. We characterize the MZR at z~3-6 within a fully Bayesian framework using JWST NIRSpec spectra of 193 galaxies from the RUBIES survey. We forward model the observed mass-metallicity surface using prospector-generated spectra to account for two selection biases: the survey selection function and success in observing high signal-to-noise emission lines. We demonstrate that the RUBIES selection function, based on F444W magnitude and F150W-F444W color, has a negligible effect on our measured MZR. A correct treatment of the non-Gaussian metallicity uncertainties from strong-line calibrations lowers the derived MZR normalization by 0.2 dex and flattens the slope by ~20%; forward-modeling the effect of emission line observability steepens the slope by ~15%. Both of these biases must be taken into account in order to properly measure the intrinsic MZR. This novel forward modeling process motivates careful consideration of selection functions in future surveys, and paves the way for robust, high-redshift chemical enrichment studies that trace the evolution of the mass-metallicity relation across cosmic time.

Radio-loud high-redshift quasars (RHRQs) provide crucial insights into the evolution of relativistic jets and their connection to the growth of supermassive black holes. Beyond the extensively studied population at $z \ge 5$, the cosmic morning epoch ($3 \lesssim z \lesssim 5$) marks the peak of active galactic nucleus (AGN) activity and black hole accretion, yet remains relatively unexplored. In this work, we compiled the radio high-redshift quasar catalog (RHzQCat) by cross-matching the SDSS DR16Q catalog with four major radio surveys -- FIRST,NVSS, RACS, and GLEAM. Our tier-based cross-matching framework and visual validation ensured reliable source identification across surveys with diverse beam sizes. The catalog included 1629 reliable and 315 candidate RHRQs, with radio luminosities uniformly spanning $10^{25.5}$ -- $10^{29.3}$ W Hz$^{-1}$. About 95\% of the confirmed sources exhibited compact morphologies, consistent with Doppler-boosted or young AGN populations at high redshifts. Our catalog increases the number of known RHRQs at $z\ge3$ by an order of magnitude, representing the largest and most homogeneous catalog of radio quasars at cosmic morning, filling the observational gap between the early ($z>6$) and local Universe. It provides a robust reference for future statistical studies of jet evolution, AGN feedback, and cosmic magnetism with next-generation facilities such as the Square Kilometer Array (SKA).

Fast Radio Bursts (FRBs) are bright, millisecond-duration extragalactic radio transients that probe extreme astrophysical environments. Many FRBs exhibit multi-component structures, which encode information about their emission mechanisms or progenitor systems and thus provide important clues to their origins. In this work, we systematically analyze the burst morphology of FRB 20190520B and compare component distributions across four active FRBs observed with FAST: FRB 20121102A, FRB 20190520B, FRB 20201124A, and FRB 20240114A. We find that multi-component burst-clusters show spectral properties similar to single-peak bursts, and no periodicity is detected in their temporal behavior. The component-count distributions follow a power law, revealing scale-free behavior consistent with self-organized criticality (SOC) processes. Multi-component clusters account for 12-30% of all detected bursts, regardless of source activity, providing new insights into burst-to-burst variability and the physical processes driving FRB emission.

Context. Young, solar analogue stars provide key insights into the early stages of stellar evolution, particularly in terms of magnetic activity and rotation. Their rapid rotation, high flaring rate, and enhanced surface activity make them ideal laboratories for testing stellar models or even the solar dynamo. Aims. Using long-term photometric data, we investigated the cyclic behaviour of EK Dra over the last century. We analyze its short- term activity based on 13 Transiting Exoplanet Survey Satellite (TESS) sectors. Applying Doppler imaging on high-resolution spectral data we investigate short and long-term spot evolution and surface differential rotation. Methods. We use Short-term Fourier-transform on a 120 years long archival photometric data in order to search for activity cycles. The short-term space photometry data is fitted with an analytic three-spot model, and we hand-select flares from it to analyze their phase and frequency distribution. Spectral synthesis is used to determine the astrophysical parameters of EK Dra. Using the iMap multi-line Doppler imaging code, we reconstruct 13 Doppler images. Differential rotation is derived by cross-correlating consecutive Doppler maps. Results. Long-term photometric data reveal a 10.7-12.1 year cycle that was persistently present for 120 years. In the more recent half of the light curve a 7.3-8.2 years-long signal is also visible. The distribution of the 142 flares in the TESS data shows no correlation with the rotational phase or with the spotted longitudes. The reconstructed Doppler images show a surface that varies from rotation to rotation, putting the lower limit of the spot lifetime between 10-15 days. Based on the cross-correlation of the Doppler maps, EK Dra has a solar-type differential rotation with a surface shear parameter of $\alpha_{DR} = 0.030 \pm 0.008$.

Spectropolarimetry provides a unique probe of ejecta asphericities, offering direct insights into the underlying explosion physics of Type Ia supernovae (SNe Ia). We analyze the statistical properties of pre-maximum spectropolarimetric data for 24 SNe Ia observed with VLT/FORS, focusing on the Si II $\lambda$6355 Åline. Previous studies have revealed a correlation between the peak Si II polarization degree and the expansion velocity. Here, we combine these observations with multi-dimensional non-LTE radiative transfer simulations. We consider two asphericity classes: (i) lopsided abundance distributions produced by off-center delayed-detonation transitions in near-$M_{Ch}$ white dwarfs or, for example, WD collisions (Class I), and (ii) global, axisymmetric density asphericities such as those arising from explosions of rapidly rotating WDs or mergers (Class II). Our model grid spans normal to subluminous SNe Ia and successfully reproduces the observed Si II velocity-polarization trend, with higher velocities associated with stronger asphericities. Consistent with observations, transitional SNe Ia and the faint end of the normal SNe Ia population show the highest Si II polarization and are best explained by Class I scenarios. In contrast, subluminous SNe Ia are dominated by Class II asphericities, characterized by lower Si II polarization but significant continuum polarization. The observed distribution of Si II polarization depends on both the observer's viewing angle $\theta$ and the intrinsic asphericity. Statistical analysis of these spectropolarimetric snapshots enables the separation of Class I and Class II contributions and highlights the intrinsic diversity among SNe Ia. Our results imply viewing-angle-dependent luminosities in our local sample, which may have implications when using high-redshift SNe Ia as evidence for the need of non-standard cosmology.

The mid-infrared spectrum of star-forming, high metallicity galaxies is dominated by emission features from aromatic and aliphatic bonds in small carbonaceous dust grains, often referred to as polycyclic aromatic hydrocarbons (PAHs). In metal-poor galaxies, the abundance of PAHs relative to the total dust sharply declines, but the origin of this deficit is unknown. We present JWST observations that detect and resolve emission from PAHs in the 7% Solar metallicity galaxy Sextans A, representing the lowest metallicity detection of PAH emission to date. In contrast to higher metallicity galaxies, the clumps of PAH emission are compact (0.5-1.5'' or 3-10 pc), which explains why PAH emission evaded detection by lower resolution instruments like Spitzer. Ratios between the 3.3, 7.7, and 11.3 $\mu$m PAH features indicate that the PAH grains in Sextans A are small and neutral, with no evidence of significant processing from the hard radiation fields within the galaxy. These results favor inhibited grain growth over enhanced destruction as the origin of the low PAH abundance in Sextans A. The compact clumps of PAH emission are likely active sites of in-situ PAH growth within a dense, well-shielded phase of the interstellar medium. Our results show that PAHs can form and survive in extremely metal-poor environments common early in the evolution of the Universe.

We analyze the properties of satellite galaxies around 1,024 Milky Way-mass hosts from the DREAMS Project, simulated within a $\Lambda$CDM cosmology. Utilizing the TNG galaxy-formation model, the DREAMS simulations incorporate both baryonic physics and cosmological uncertainties for a large sample of galaxies with diverse environments and formation histories. We investigate the relative impact of the physical uncertainty from the galaxy-formation model on predicted satellite properties using four metrics: the satellite stellar mass function, radial distribution, inner slope of dark matter density profile, and stellar half-light radius. We compare these predictions to observations from the SAGA Survey and the DREAMS N-body simulations and find that uncertainties from baryonic physics modeling are subdominant to the scatter arising from halo-to-halo variance. Where baryonic modeling does affect satellites, the supernova wind energy has the largest effect on the satellite properties that we investigate. Specifically, increased supernova wind energy suppresses the stellar mass of satellites and results in more extended stellar half-light radii. The adopted wind speed has only a minor impact, and other astrophysical and cosmological parameters show no measurable effect. Our findings highlight the robustness of satellite properties against uncertainties in baryonic physics modeling.

We report the discovery of a doubly-imaged Little Red Dot (LRD) candidate behind the galaxy cluster Abell 383, which we dub A383-LRD1. Initially classified as a dropout galaxy in HST imaging with several ground-based emission line detections placing it at $z_{\mathrm{spec}}=6.027$, new JWST/NIRCam observations taken as part of the cycle 4 VENUS survey now reveal that the source consists of two underlying components: A red point-source with a V-shaped SED consistent with LRD selection criteria, and a nearby ($\sim 380$ pc) compact blue companion which was the main contributor to the previous rest-frame UV detections. Based on lensing symmetry and its SED, the LRD appears to lie at a similar redshift as well. The magnification of the two images of A383-LRD1 is $\mu_{\mathrm{A}}=16.2\pm1.2$ and $\mu_\mathrm{B}=9.0\pm0.6$, respectively, and the predicted time delay between them is $\Delta t_{\mathrm{grav}}=5.20\pm0.14$ yr ($\sim0.7$ yr in the rest-frame). After correcting for the lensing magnification, we derive an absolute magnitude of $M_{\mathrm{UV,LRD}}=-16.8\pm 0.3$ for the LRD, and $M_{\mathrm{UV,BC}}=-18.2\pm 0.2$ for the blue companion. We perform SED fits to both components, revealing the LRD to be best fitted with a black hole star (BH*) model and a substantial host galaxy, and the blue companion with an extremely young, emission-line dominated star-forming nebula. A383-LRD1 represents the second known multiply-imaged LRD detected to date, following A2744-QSO1, and to our knowledge, the first LRD system with a confirmed detection of [C $_{II}$]$\lambda158 \ \mu$m emission from ALMA observations. Thanks to lensing magnification, this system opens a unique door to study the relation between a LRD, its host galaxy, and its environment, and represents a prime candidate for deep JWST spectroscopy and high-resolution ALMA follow-up observations.

We present the quasar catalog from Data Releases 10 to 12 of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Quasar Survey, comprising quasars observed between September 2021 and June 2024. We robustly identified $11,346$ quasars, of which $5,386$ are newly discovered objects not present in the Million Quasars catalog. This release brings the total number of quasars identified by the 12-year LAMOST survey to $67,521$, of which $29,513$ are newly discovered. While the absolute flux calibration for LAMOST quasar spectra from Data Releases 6 to 9 was previously performed using the SDSS/PanSTARRS1 multi-band photometric data, the inherent variability of quasars can affect the flux accuracy. To address this limitation, we recalibrated the LAMOST spectra using (quasi-)simultaneous photometric data from Zwicky Transient Facility (ZTF), which has conducted high-cadence sky monitoring since March 2018. Based on the recalibrated single-epoch spectra, we estimated the emission line fluxes, continuum fluxes, and virial black hole masses. These improved spectra facilitate direct comparison with the spectra of common quasars from the Sloan Digital Sky Survey (SDSS), enabling searches for rare quasars, such as changing-look quasars exhibiting the appearance or disappearance of broad emission lines and broad absorption line quasars. The combined dataset of photometry and multi-epoch spectra will enhance the detections of AGN-related transients, such as Bowen fluorescence flares and extreme variability quasars, thereby improving our understanding of quasar variability.

Baryon feedback redistributes gas relative to the underlying dark matter distribution and suppresses the matter power spectrum on small scales, but the amplitude and scale dependence of this effect are uncertain. We constrain the impact of baryon feedback on the matter power spectrum by jointly analysing X-ray gas mass fractions from the eROSITA and HSC-XXL samples and SDSS/DESI+ACT kinetic Sunyaev-Zel'dovich (kSZ) effect profiles; the samples are characterised with galaxy-galaxy lensing and together span group and cluster masses at $0<z<1$. Using the baryonification framework, our joint eROSITA and kSZ model gives precise constraints on the suppression of the matter power spectrum: $10 \pm 2\%$ at $k=1~h~\mathrm{Mpc}^{-1}$. The inferred gas profiles are more extended and the power suppression is stronger than predicted by the fiducial models of recent hydrodynamical simulation suites, including FLAMINGO and BAHAMAS. The HSC-XXL gas mass fractions, which the fiducial simulations were calibrated to reproduce, prefer more moderate power suppression than the kSZ and eROSITA data: $5 \pm 4\%$ at $k=1~h~\mathrm{Mpc}^{-1}$. With a simulated LSST Year 1 weak lensing analysis, we demonstrate a framework for next-generation surveys: calibrating feedback models with multi-wavelength gas observables to recover the small-scale statistical power of cosmic shear.

The $y$-type distortion of the blackbody spectrum of the cosmic microwave background radiation probes the pressure of the gas trapped in galaxy groups and clusters. We reanalyze archival data of the FIRAS instrument with an improved astrophysical foreground cleaning technique, and measure a mean $y$-distortion of $\langle y\rangle = (1.2\pm 2.0) \times 10^{-6}$ ($\langle y\rangle\lesssim 5.2\times 10^{-6}$ at 95\% C.L.), a factor of $\sim 3$ tighter than the original FIRAS results. This measurement directly rules out many models of baryonic feedback as implemented in cosmological hydrodynamical simulations, mostly using information in objects with mass $M\lesssim 10^{14} {\rm M}_{\odot}$. We discuss its implications for the analysis of cosmic shear and kinetic Sunyaev-Zel'dovich effect data, and future spectral distortion experiments.

We introduce a new suite of 1,024 cosmological and hydrodynamical zoom-in simulations of Milky Way-mass halos, run with Cold Dark Matter, as part of the DREAMS Project. Each simulation in the suite has a unique set of initial conditions and combination of cosmological and astrophysical parameters. The suite is designed to quantify theoretical uncertainties from halo-to-halo variance, as well as stellar and black hole feedback. We develop a novel weighting scheme that prioritizes regions of the input parameter space, yielding galaxies consistent with the observed present-day stellar mass--halo mass relation. The resulting galaxy population exhibits a wide diversity in structural properties that encompasses those of the actual Milky Way, providing a powerful statistical sample for galactic archaeology. To demonstrate the suite's scientific utility, we investigate the connection between a galaxy's merger history, focusing on Gaia-Sausage-Enceladus~(GSE) analogs, and its present-day properties. We find that galaxies with a GSE analog have lower star formation rates, more compact disks, and more spherical stellar halos. Crucially, significant halo-to-halo scatter remains, demonstrating that matching more than the most significant events in the Milky Way's past is necessary to recover its present-day properties. Our results highlight the necessity for large statistical samples to disentangle the stochastic nature of galaxy formation and robustly model the Milky Way's unique history.

Recent discoveries of streamer-like structures around protostellar sources challenge the traditional picture of isolated, axisymmetric star formation. Here, we present new ALMA observations of [BHB2007]1, a flat-spectrum source connected to at least three such elongated structures. Two of these features are symmetrically located to the north and south of the disk, with velocities aligned with the disk on their respective sides. However, their unbound kinematics and curved morphology make it difficult to determine their origin. Possible explanations include outflows, interactions with the nearby BHB2 system, and hyperbolic infall, but none fully account for all observed properties. In contrast, a newly identified collimated structure to the west shows clear evidence of gravitationally bound infall. Estimates of its mass, mass infall rate, and angular momentum suggest that this infalling streamer would roughly double the mass budget available to form planets and tilt the disk by a few tens of degrees. Furthermore, its misalignment with the midplane of the disk and the lack of diffuse envelope emission indicate that the streamer may have formed due to gravitational capture of cloud material unrelated to the source's natal core. Together, these findings support a more dynamic picture of star formation, one where environmental interactions continue to shape conditions for building planetary systems.

Type II orbital migration is a key process to regulate the mass and semimajor axis distribution of exoplanetary giant planets. The conventional formula of type II migration generally predicts too rapid inward migration to reconcile with the observed pile-up of gas giant beyond 1 au. Analyzing the recent high-resolution hydrodynamical simulations by Li et al. (2024) and Pan et al. (2025) that show robust outward migration of a gas accreting planet, we here clarify the condition for the outward migration to occur and derive a general semi-analytical formula that can be applied for broad range of planet mass and disk conditions. The striking outward migration is caused by azimuthal asymmetry in corotation torque exerted from cicumplanetary disk regions (connecting to horseshoe flow) that is produced by the planetary gas accretion, while the conventional inward migration model is based on radial asymmetry in the torques from the circumstellar protoplanetry disk. We found that the azimuthal asymmetry dominates and the migration is outward, when the gap depth defined by the surface density reduction factor of $1/(1+K')$ is in the range of $0.03 \lesssim K' \lesssim 50$. Using simple models with the new formula, we demonstrate that the outward migration plays an important role in shaping the mass and semimajor axis distribution of gas giants. The concurrent dependence of planets' accretion rate and migration direction on their masses and disk properties potentially reproduces the observed pile-up of exoplanetary gas giants beyond 1 au, although more detailed planet population synthesis calculations are needed in the future.

Type Ia supernovae (SNe Ia) are among the most precise cosmological distance indicators used to study the expansion history of the Universe. The vast increase of SN Ia data due to large-scale astrophysical surveys has led to the discovery of a wide variety of SN Ia sub-classes, such as transitional and fast-declining SNe Ia. However, their distinct photometric and spectroscopic properties differentiate them from the population of normal SNe Ia such that their use as cosmological tools remains challenged. Here, we present a high-cadenced photometric and spectroscopic dataset of two SNe Ia, SNe 2020ue and 2020nlb, which were discovered in the nearby Virgo cluster of galaxies. Our study shows that SN 2020nlb is a normal SN Ia whose unusually red color is intrinsic, arising from a lower photospheric temperature rather than interstellar reddening, providing clear evidence that color diversity among normal SNe Ia can have a physical origin. In contrast, SN 2020ue has photometric properties, such as color evolution and light-curve decay rate, similar to those of transitional SNe, spectroscopically it is more aligned with normal SNe Ia. This is evident from spectroscopic indicators such as the pseudo-equivalent width of \ion{Si}{II} lines. Thus, such SNe Ia that are photometrically at the edge of the standard normal SNe Ia range may be missed in cosmological SNe Ia samples. Our results highlight that spectroscopic analysis of SNe Ia around peak brightness is crucial for identifying intrinsic color variations and constructing a more complete and physically homogeneous SN Ia sample for precision cosmology.

Sgr A*, the supermassive black hole at the center of the Milky Way, exhibits frequent short-duration flares with luminosity greater than 1e34 erg/s across multiple wavelengths. The origin of the flares is still unknown. We revisited the role of small planetary bodies, originally from the stellar disk, and their tidally disrupted fragments as a source of flaring activity in Sgr A*. We refined previous models by incorporating material strength constraints on the tidal disruption limit and by evaluating the evaporation dynamics of the resulting fragments. We analyzed the tidal fragmentation and gas-induced fragmentation of small planetary bodies with rubble-pile and monolithic structures. Using constraints from recent space missions (e.g., NASA OSIRIS-REx and JAXA Hayabusa2), we estimated the survivability of fragments under aerodynamic heating and computed their expected luminosity from ablation, modeled as fireball flares analogous to meteor events. We find that planetary fragments can approach as close as 8 gravitational radii, consistent with observed flare locations. The fireball model yields luminosities from 1e34 to 1e36 erg/s for fragments whose parent bodies are a few kilometers in size. The derived flare frequency vs. luminosity distribution follows a power law with index 1.83, in agreement with observed values (1.65 - 1.9), while the flare duration scales as L^(-1/3), consistent with observations. We consider the young stars around Sgr A* as the planetary reservoir. Given a small-body population analogous in mass to the primordial Kuiper belt and the common existence of close-in super-Earths and long-period Neptunes, we show that this planetary reservoir can supply the observed flares.

The Euclid system performance is defined in terms of image quality metrics tuned to the weak gravitational lensing (WL) cosmological probe. WL induces stringent requirements on the shape and stability of the VIS instrument system point spread function (PSF). The PSF is affected by error contributions from the telescope, the focal plane and image motion, and is controlled by a global error budget with error allocations to each contributor. Aims. During spacecraft development, we verified through a structural-thermal-optical performance (STOP) analysis that the built and verified telescope with its spacecraft interface meets the in-orbit steady-state and transient image quality requirements. Methods. For the purposes of the STOP analysis, a detailed finite-element mathematical model was set up and a standard set of test cases, both steady-state and transient, was defined, comprising combinations of worst-case boundary conditions. Results. The STOP analysis addressed the interaction of all spacecraft components in transmitting temperature-induced loads that lead to optical train deformation. The results of the prelaunch analysis demonstrated that temperature-induced optical perturbations will be well below the allowable limits for all permitted observing conditions. During the first year in orbit, we used the STOP analysis predictions to help interpret the measured performance as a function of environmental variables. Unpredicted disturbances were discovered and unexpected sensitivities were revealed. In-orbit temperature variations are small (<300 mK) and so are their effects on the telescope structure, but they are detected in the time histories of the image quality metrics and are a non-negligible factor in the PSF stability budget demanded by the WL science. Taking everything into account, our analysis confirms the excellent overall performance of the telescope.

We present 29 successfully recovered CIV time lags in Active Galactic Nuclei from the complete Dark Energy Survey Reverberation Mapping campaign. The AGN in this sample span a redshift range of 1.9<z<3.5. We successfully measure the velocity dispersion from the CIV spectral linewidth for 25 of these 29 sources, and use these to calculate new high-redshift black hole mass estimates, finding masses between 0.8 and 1.3 billion solar masses. We also identify a selection effect due to the duration of the survey that can impact the radius-luminosity relation derived from this and other (high-redshift) data. This paper represents the culmination of the OzDES CIV campaign.

Over the last decade, the Australian Dark Energy (OzDES) collaboration has used Reverberation Mapping to measure the masses of high redshift supermassive black holes. Here we present the final review and analysis of this OzDES reverberation mapping campaign. These observations use 6-7 years of photometric and spectroscopic observations of 735 Active Galactic Nuclei (AGN) in the redshift range $z\in [0.13, 3.85]$ and bolometric luminosity range $\log_{10}(L_{\mathrm{bol}})\in [44.3, 47.5] \mathrm{erg/s}$. Both photometry and spectra are observed in visible wavelengths, allowing for the physical scale of the AGN broad line region to be estimated from reverberations of the Hbeta, MgII and CIV emission lines. We successfully use reverberation mapping to constrain the masses of 62 super-massive black holes, and combine with existing data to fit a power law to the lag-luminosity relation for the Hbeta and MgII lines with a scatter of $\sim0.25$ dex, the tightest and most robust fit yet identified. We fit a similarly constrained relation for CIV, resolving a tension with the low luminosity literature AGN by accounting for selection effects. We also examine the impact of emission line width and luminosity (related to accretion rate) in reducing the scatter of these scaling relationships and find no significant improvement over the lag-only approach for any of the three lines. Using these relations, we further estimate the masses and accretion rates of 246 AGN. We also use these relations to estimate the relative sizes of the Hbeta, MgII and CIV emitting regions, and find evidence that the MgII emission may occur further out than Hbeta. In short, we provide a comprehensive benchmark of high redshift AGN reverberation mapping at the close of this most recent generation of surveys, including light curves, time-delays, and the most reliable radius-luminosity relations to date.

Polarized radiation serves as a vital diagnostic tool in astrophysics, providing unique insights into magnetic field geometries, scattering processes, and three-dimensional structures in diverse astrophysical scenarios. To address these applications, we present Kratos-polrad, a novel GPU-accelerated Monte Carlo Radiative Transfer code built upon the heterogeneous computing framework of Kratos, designed for self-consistent and efficient polarization calculations. It utlizes comprehensive treatment of Stokes parameters throughout photon propagation, featuring transforms the grain-lab frame transforms using quaternion algebra and consistent non-linear polarization extinction in cells, which are useful in modeling radiative transfer processes with scatterings by aligned dust grains. The code implements two-step polarimetry imaging that decouples Monte Carlo sampling of scattering physics from imaging geometry, enabling efficient synthesis maximizing the utilization of photon packets. Extensive validation against analytical solutions and established codes demonstrates accurate treatment of diverse polarization phenomena, including self-scattering polarization, dichroic extinction in aligned dust grains, and complex polarization patterns in twisted magnetic field configurations. By leveraging massive GPU parallelism, optimized memory access patterns, and analytical approaches for optically thick cells, Kratos-polrad achieves performance improvements of $\sim 10^{2}$ times compared to CPU-based methods, enabling previously prohibitive studies in polarimetric astrophysics.

Our comprehension of the history of star formation at $z>3$ relies on rest-frame UV observations, yet this selection misses the most dusty and massive sources, yielding an incomplete census at early times. Infrared facilities such as Spitzer and the James Webb Space Telescope have revealed a hidden population at $z=3$-$6$ with extreme red colours, named HIEROs (HST-to-IRAC extremely red objects), identified by the criterion $H_{\mathrm{E}}-\mathrm{ch2}>2.25$. Recently, Euclid Early Release Observations (ERO) have made it possible to further study such objects by comparing Euclid data with ancillary Spitzer/IRAC imaging. We investigate a $232$ arcmin$^2$ area in the Perseus field using VIS and NISP photometry, complemented by the four Spitzer channels and ground-based MegaCam bands ($u$, $g$, $r$, ${\rm H}\alpha$, $i$, $z$). Applying the colour cut yields $121$ HIEROs; after removing globular clusters, brown dwarfs, and unreliable cases through visual inspection of multiband cutouts, we obtain a final sample of $42$ robust HIEROs. Photometric redshifts and physical properties are estimated with the SED-fitting code Bagpipes. From the resulting $z_{\mathrm{phot}}$ and $M_*$ values, we compute the galaxy stellar mass function at $3.5<z<5.5$. Even after excluding possible AGN hosts or systems where the stellar mass may be overestimated, the high-mass end remains comparable to previous determinations, suggesting the true abundance could be higher. These results highlight the importance of further study of this obscured population to assess its role in the cosmic star-formation rate density and its consistency with galaxy-formation models, demonstrating Euclid's capability to advance our understanding of dust-hidden star formation across early epochs.

Recent observations of radio-quiet active galactic nuclei (RQAGN) have shown the presence of millimeter emission, whose origin remains unknown, from within parsec scales of the central black hole. We argue that the mm emission comes from a spatially extended region that is magnetically connected to the compact X-ray corona, in analogy to the solar wind and corona. We present an analytic model scaled to corona values in which non-equipartition electrons from multiple heights along an extended conical outflow shape the mm emission. In this model, the 100 GHz emission originates from within $\lesssim10^4$ gravitational radii ($r_g$) of the central black hole, though the projected distance from the black hole can be as low as $50r_g$ depending on the line-of-sight. Our model predicts a flat emission spectrum $F_{\nu}\sim{\rm const}$ and a mm-to-X-ray luminosity ratio $L_{\rm mm}/L_X\sim10^{-4}$, consistent with observations. These quantities depend weakly on the underlying electron power-law distribution function and black hole mass. We demonstrate this model's plausibility using a general relativistic magneto-hydrodynamic (GRMHD) simulation of a thin accretion disc as a case study. Our model highlights the need to study continual dissipation along the outflow to connect the X-ray- and mm-emitting regions.

We use the DAWN JWST Archive to construct and characterise a sample of 116 little red dots (LRDs) across 2.3<z<9.3, selecting all sources with v-shaped UV-optical continua from NIRSpec/PRISM spectra and compact morphologies in NIRCam/F444W imaging. We show that LRD continuum spectra are ubiquitously well described by modified blackbodies across ~$0.4-1.0\mu$m, with typical T~5000K or $\lambda_{peak}$~$0.65\mu$m across 2 dex in luminosity, and a tail toward T~2000K. LRDs therefore trace a locus in the Hertzsprung-Russell diagram that is directly analogous to stars on the Hayashi track, strongly supporting the picture that LRDs are AGN embedded in thermalised dense gas envelopes in approximate hydrostatic equilibrium. Hotter LRDs with $\lambda_{peak}<0.65\mu$m typically have strong Balmer breaks, redder UV slopes and high optical luminosities; other LRDs show weak or no Balmer breaks, and wide variety in $\beta_{UV}$ and $L_{5100}$. Crucially, we demonstrate that the UV-optical continuum shapes and luminosities are strongly linked to the $H\alpha,\ H\beta$, [OIII] and OI line properties. There is a tight linear relation between the H$\alpha$ and optical continuum luminosities, as well as H$\alpha$ and OI$_{8446}$, indicating that Balmer, OI and optical emission must primarily be powered by the same source. The Balmer decrement increases strongly toward higher $L_{H\alpha}$, $L_{5100}$ and Balmer break strength, providing key evidence for luminosity-dependent effects of collisional (de-)excitation and resonant scattering in the gaseous envelopes. In contrast, we show that [OIII] emission likely originates from star-forming host galaxies, and that its strong correlation with Balmer break strength arises naturally from variation in the AGN-to-host ratio. Our work presents an empirical description of the nature and structure of LRDs, defining a new benchmark for ongoing LRD model developments.

Large-scale $B$-mode polarization of the cosmic microwave background (CMB) is a prime target for current and future experiments in search of primordial gravitational waves (PGW). With increasingly sensitive instruments being deployed, secondary $B$-modes induced by the weak gravitational lensing of CMB photons are becoming an important limitation and need to be removed, a process known as delensing. In this work, we combine internally reconstructed CMB lensing maps from the Atacama Cosmology Telescope (ACT) DR6 with galaxy samples from unWISE and a map of the cosmic infrared background (CIB) fluctuations from Planck to produce a well-correlated tracer of the CMB lensing field. Our co-added tracer, shown to be 55-85% correlated with the true lensing convergence at multipoles $L \leq 2000$, is then convolved with ACT DR6 $E$-mode polarization to yield a template of the lensing $B$-modes. We assess its performance on a wide range of scales by using it to delens ACT DR6 and Planck $B$-modes over 23% of the sky, removing around 39% of the lensing power at $100\leq l \leq 1500$ and 47% at $30 \leq l \leq 300$, respectively. Our template achieves the highest delensing efficiency to date and will be useful for the analysis of early polarization maps from the Simons Observatory (SO). We finally outline prospects for further improvements by including additional large-scale structure tracers from upcoming cosmological surveys.

The WISE and NEOWISE missions have provided the only mid-infrared all-sky time-domain data, opening a unique observational window for variability studies. Yet, a comprehensive and systematic catalog of mid-infrared variable sources has remained unavailable. In this work, we construct the first large-scale mid-infrared variability catalog based on the unTimely coadded photometry, covering tens of millions of sources. By employing a Bayesian Gaussian mixture model with a Dirichlet process, we identified 8,256,042 variable sources in the W1 band and 7,147,661 in the W2 band, significantly expanding the landscape of known mid-infrared variables. In addition to robust variability metrics, our analysis highlights rare and extreme outliers through dedicated outlier-detection algorithms, enabling the discovery of unusual classes of objects such as eruptive young stellar objects, highly variable active galactic nuclei, and other rare transients. This unprecedented dataset provides a new foundation for time-domain astronomy in the mid-infrared, offering complementary insights to optical and near-infrared surveys, and opening the door to systematic investigations of stellar evolution, accretion processes, and dust-enshrouded astrophysical environments on a Galactic and extragalactic scale.

We analyzed the spatially resolved and global star formation histories (SFHs) for a sample of 25 TNG50-SKIRT Atlas galaxies to assess the feasibility of reconstructing accurate SFHs from Euclid-like data. This study provides a proof of concept for extracting the spatially resolved SFHs of local galaxies with Euclid, highlighting the strengths and limitations of SFH modeling in the context of next-generation galaxy surveys. We used the spectral energy distribution (SED) fitting code Prospector to model both spatially resolved and global SFHs using parametric and nonparametric configurations. The input consisted of mock ultraviolet--near-infrared photometry derived from the TNG50 cosmological simulation and processed with the radiative transfer code SKIRT. We show that nonparametric SFHs provide a more effective approach to mitigating the outshining effect by recent star formation, offering improved accuracy in the determination of galaxy stellar properties. Also, we find that the nonparametric SFH model at resolved scales closely recovers the stellar mass formation times (within 0.1~dex) and the ground truth values from TNG50, with an absolute average bias of $0.03$~dex in stellar mass and $0.01$~dex in both specific star formation rate and mass-weighted age. In contrast, larger offsets are estimated for all stellar properties and formation times when using a simple $\tau$-model SFH, at both resolved and global scales, highlighting its limitations. These results emphasize the critical role of nonparametric SFHs in both global and spatially resolved analyses, as they better capture the complex evolutionary pathways of galaxies and avoid the biases inherent in simple parametric models.

The Solar Close Observations and Proximity Experiments (SCOPE) mission will send a spacecraft into the solar atmosphere at a low altitude of just 5 R_sun from the solar center. It aims to elucidate the mechanisms behind solar eruptions and coronal heating, and to directly measure the coronal magnetic field. The mission will perform in situ measurements of the current sheet between coronal mass ejections and their associated solar flares, and energetic particles produced by either reconnection or fast-mode shocks driven by coronal mass ejections. This will help to resolve the nature of reconnections in current sheets, and energetic particle acceleration regions. To investigate coronal heating, the mission will observe nano-flares on scales smaller than 70 km in the solar corona and regions smaller than 40 km in the photosphere, where magnetohydrodynamic waves originate. To study solar wind acceleration mechanisms, the mission will also track the process of ion charge-state freezing in the solar wind. A key achievement will be the observation of the coronal magnetic field at unprecedented proximity to the solar photosphere. The polar regions will also be observed at close range, and the inner edge of the solar system dust disk may be identified for the first time. This work presents the detailed background, science, and mission concept of SCOPE and discusses how we aim to address the questions mentioned above.

Pulsar timing arrays (PTAs) have recently entered the detection era, quickly moving beyond the goal of simply improving sensitivity at the lowest frequencies for the sake of observing the stochastic gravitational wave background (GWB), and focusing on its accurate spectral characterization. While all PTA collaborations around the world use Fourier-domain Gaussian processes to model the GWB and intrinsic long time-correlated (red) noise, techniques to model the time-correlated radio frequency-dependent (chromatic) processes have varied from collaboration to collaboration. Here we test a new class of models for PTA data, Gaussian processes based on time-domain kernels that model the statistics of the chromatic processes starting from the covariance matrix. As we will show, these models can be effectively equivalent to Fourier-domain models in mitigating chromatic noise. This work presents a method for Bayesian model selection across the various choices of kernel as well as deterministic chromatic models for non-stationary chromatic events and the solar wind. As PTAs turn towards high frequency (>1/yr) sensitivity, the size of the basis used to model these processes will need to increase, and these time-domain models present some computational efficiencies compared to Fourier-domain models.

Synergies between large-scale radio-continuum and optical/near-infrared galaxy surveys are a powerful tool for cosmology. Cross-correlating these surveys can constrain the redshift distribution of radio sources, mitigate systematic effects, and place constraints on cosmological models. We perform the first measurement of the clustering cross-spectrum between radio-continuum sources in the Evolutionary Map of the Universe (EMU) survey and galaxies from the ESA Euclid satellite mission's Q1 release. Our goal is to detect and characterise the cross-correlation signal, test its robustness against systematic effects, and compare our measurements with theoretical predictions. We use data from the Australian SKA Pathfinder's EMU Main Survey, which overlaps with the Euclid Deep Field South. We generate two radio-source catalogues using different source finders to create galaxy maps. We measure the harmonic-space cross-correlation signal using a pseudo-spectrum estimator. The measured signal is compared to theoretical predictions based on a {\Lambda}CDM cosmology, using several models for the EMU source redshift distribution and bias. We report detection above 8{\sigma} of the cross-correlation signal consistent across all tested models and data sets. The measured cross-spectra from the two radio catalogues are in excellent agreement, demonstrating that the cross-correlation is robust against the choice of source-finding algorithm. The measured signal also agrees with theoretical models developed from previous cross-correlation studies and simulations. This pathfinder study establishes a statistically significant cross-correlation between EMU and Euclid. The robustness of the signal is a crucial validation of the methodology, paving the way for future large-scale analyses leveraging the full power of this synergy to constrain cosmological parameters and our understanding of galaxy evolution.

Feedback from active galactic nuclei (AGN) is believed to play a significant role in suppressing cooling flows in cool-core (CC) clusters. Turbulence in the intracluster medium (ICM), which may be induced by AGN activity or pre-existing motions, has been proposed as a potential heating mechanism based on analysis of Chandra X-ray surface brightness fluctuations. However, subsequent simulation results have found the subdominant role of turbulence in heating the ICM. To investigate this discrepancy, we perform three-dimensional hydrodynamic simulations of a Perseus-like cluster including both AGN feedback and pre-existing turbulence, which is stirred to the observationally constrained level in the Perseus cluster. Our results indicate that, although the velocity field is dominated by the pre-existing turbulence, AGN heating through bubbles and shocks remains significant. More importantly, analysis of the velocity structure function and the energy power spectrum shows that the turbulent heating rate is smaller than the radiative cooling rate, especially in the cluster core. Our results offer insights relevant for recent XRISM observations and indicate that turbulent heating alone cannot offset radiative cooling in CC clusters.

In this work, we explore in a consistent fashion the effects of fast flavor conversion (FFC) in 1D and 2D core-collapse supernova (CCSN) simulations. In addition, we investigate the impact of various angular reconstruction methods and compare the ``3-species'' and ``4-species'' neutrino transport schemes. We find that the FFC effects are insensitive to the different methods tested and that the FFC alters supernova hydrodynamics is only minor ways. We also present a ``quasi-equipartition'' approximation which can be used to estimate the FFC-altered neutrino properties by post-processing the neutrino signals extracted from no-oscillation CCSN simulations. The relative errors in neutrino number and energy luminosities of this phenomenological method are less than 2\% for 1D models, and less than 10\% for 2D models. This method provides a simple way to include the effects of FFC on neutrino signals without implementing a complex and expensive FFC scheme or redoing simulations.

Dust production is a fundamental aspect of the baryonic cycle of star formation. It is known that dust is injected into the interstellar medium during early star formation by supernovae and later on by evolved stars. From individual objects, these mechanisms are well understood, but the overall dust production in star clusters at different evolutionary stages is still challenging to quantify. We present 22 massive (> 105M$_{\odot}$) extra galactic star clusters with ages between 3 and 100 Myr exhibiting a compact dust morphology seen with JWST-MIRI. We only find PAH features associated with one star cluster and nineteen have already cleared themselves from their natal dust. Their main characteristic is a significant enhancement at 10${\mu}$m, which is likely due to silicate emission and cannot be explained by ionized gas. We discuss several possible explanations including dust production from evolved stars such as red super giants, more exotic star types like yellow hypergiants and luminous blue variable stars. Stochastic dust injection from supernovae or a single supernova in dense gas can also create significant silicate emission. However, for this scenario secondary tracers such as a X-ray signal are expected which we only observe in three star clusters. We find the most luminous 10${\mu}$m emitter to be the three most massive star clusters (> 106M$_{\odot}$) which is at least a magnitude stronger than any known stellar sources indicating a rare mechanism that only appears at extreme masses and a short lifetime.

We present the strategy to identify and mitigate potential sources of angular systematics in the Euclid spectroscopic galaxy survey, and we quantify their impact on galaxy clustering measurements and cosmological parameter estimation. We first survey the Euclid processing pipeline to identify all evident, potential sources of systematics, and classify them into two broad classes: angular systematics, which modulate the galaxy number density across the sky, and catastrophic redshift errors, which lead to interlopers in the galaxy sample. We then use simulated spectroscopic surveys to test our ability to mitigate angular systematics by constructing a random catalogue that represents the visibility mask of the survey; this is a dense set of intrinsically unclustered objects, subject to the same selection effects as the data catalogue. The construction of this random catalogue relies on a detection model, which gives the probability of reliably measuring the galaxy redshift as a function of the signal-to-noise ratio (S/N) of its emission lines. We demonstrate that, in the ideal case of a perfect knowledge of the visibility mask, the galaxy power spectrum in the presence of systematics is recovered, to within sub-percent accuracy, by convolving a theory power spectrum with a window function obtained from the random catalogue itself. In the case of only approximate knowledge of the visibility mask, we test the stability of power spectrum measurements and cosmological parameter posteriors by using perturbed versions of the random catalogue. We find that significant effects are limited to very large scales, and parameter estimation remains robust, with the most impacting effects being connected to the calibration of the detection model.

We present the observational evidence of the existence of a double-decker filament channel (FC) by using observations in extreme ultraviolet and H{\alpha} wavelengths. For both FCs, the east foot-point roots in the active region (AR), while the west one roots in the remote quiet region. The bottom FC (FC1) appears as intermittent filaments. Within the AR, the FC1 appears as an S-shaped filament (F1), which consisted of two J-shaped filaments (F1S/F1N for the south/north one). For the upper one (FC2), only the east part is filled with dark plasma and visible as a small filament (F2). Its east foot-point roots around the junction of F1S and F1N. Initially, due to the recurrent reconnections, F1N and F1S link to each other and form a new filament (F3) thread by thread. Meanwhile, the heated plasma, which appears as brightening features, flows from the east foot-point of F2 to the west, and becomes invisible about 1.1$\times$10^{5} km away. The failed eruption of F1S is triggered by the reconnection, which appears as the brightening threads changing their configuration from crossed to quasiparallel in between the F1S and F3, and is confined by the upper magnetic field. Associated with the eruption, the distant invisible plasma becomes visible as a brightening feature. It continuously flows to the remote foot-point, and becomes invisible before reaching it. The brightening plasma flow outlines the skeleton of FC2 gradually. The observations show the existence of a double-decker FC, as a magnetic structure, before they appear as a brightening/dark feature when fully filled with hot/cool plasma.

Cosmic rays (CRs) streaming in weakly magnetized plasmas can drive large-amplitude magnetic fluctuations via nonresonant streaming instability (NRSI), or Bell instability. Using one-dimensional kinetic simulations, we investigate how mono-energetic and power-law CR momentum distributions influence the growth and saturation of NRSI. The linear growth is governed solely by the CR current and is largely insensitive to the CR distribution. However, the saturation depends strongly on the CR distribution and is achieved through CR isotropization, which quenches the driving current. Mono-energetic CRs effectively amplify the magnetic field and isotropize. For power-law distributions, the lowest-energy CRs dominate current relaxation and magnetic growth, while the highest-energy CRs remain weakly scattered, limiting their contribution to saturation. In the absence of low-energy CRs, high-energy particles amplify magnetic fields effectively and isotropize. We provide a modified saturation prescription accounting for these effects and propose a layered CR-confinement scenario upstream of astrophysical shocks, relevant to particle acceleration to high energies.

Migration typically occurs during the formation of planets and is closely linked to the planetary formation process. In classical theories of non-accreting planetary migration, both type I and type II migration typically result in inward migration, which is hard to align with the architecture of the planetary this http URL this work, we conduct systematic, high-resolution 3D/2D numerical hydrodynamic simulations to investigate the migration of an accreting planet. Under different disk conditions, we compared the dynamical evolution of planets with different planet-to-star mass ratios. We find that accretion of planets can significantly diminish the inward migration tendency of planets, or even change the migration direction. The migration of low-/high-mass planets is classified as Type I/II inward migration, respectively, while intermediate-mass planets, which have the strongest accretion, show an outward migration trend. We confirm that the outward migration is mainly attributed to the positive torque from the azimuthal asymmetric structures around the accreting planet, similar to Li et al. (2024). The termination of planetary mass growth is thus synonymous with the transition from outward to inward migration. For the high viscosity $\alpha=0.04$ and disk aspect ratio height $h_0=0.05$ cases, the mass ratio range for planetary outward migration is $1\times10^{-4}\lesssim q\lesssim4\times10^{-3}$. For the low viscosity case with $\alpha=0.001$, and/or the low disk aspect ratio cases $h_0=0.03$, the mass ratio range for the outward migration will shift toward the lower end. Our parameter survey reveals that a simple gap opening parameter determines the outward migration condition; details of the analytical interpretation are presented in Ida et al. (2025).