Papers for Monday, Mar 30 2026

The majority of terrestrial direct detection experiments for Dark Matter (DM) rely on the Standard Halo Model (SHM), which assumes the local DM velocity distribution follows a Maxwell-Boltzmann distribution. However, galaxy mergers can deposit DM that remains kinematically clustered today, inducing deviations from the smooth SHM prediction. Previous studies have suggested that the local stellar velocity distribution may serve as a tracer for DM populations originating from the same progenitor systems. In this work, we systematically investigate how merger mass and accretion time affect the correlation between local stellar and DM velocity distributions in Milky Way-like galaxies from the FIRE-2 simulations. We find a strong correlation between traceable DM components and their stellar counterparts, with the tightest correspondence arising from lower-mass mergers accreted at earlier cosmic times. For the remaining DM that lacks an identifiable stellar counterpart, which dominate the full DM fraction, we find that its velocity distribution is well described by a component-wise generalized Gaussian. Combining these two ingredients, we reconstruct the full local DM velocity distribution. This framework captures merger-induced features-such as co-rotation of accreted material with the galactic disk-that are entirely absent in the SHM. Finally, we propagate uncertainties through the reconstruction and show that they are dominated by the stellar mass-halo mass relation, which is unlikely to improve substantially in the near term. We therefore argue that this framework approaches the current limit of our ability to characterize the local DM velocity distribution.

Small-scale physics in the intergalactic medium (IGM) plays a crucial role in shaping the progress of cosmic reionization and several high-redshift observables that probe this period. Several recent studies have characterized the complex, dynamical response of the IGM to reionization down to kilo-parsec scales, including its effect on observables such as the Ly$\alpha$ forest. However, there has been no concentrated attempt to simulate and characterize these effects across the full parameter space of realistic large-scale IGM environments during reionization. To meet this need, we introduce the SAGUARO simulation suite, sub-titled ``Simulating IGM Evolution and Environments At High Resolution''. SAGUARO is a suite of over two hundred high-resolution, coupled radiative-hydrodynamics simulations of IGM gas dynamics during and after reionization. The suite spans a grid of photoionization rates, redshifts of reionization, and box-scale densities. We also simulate other physical effects, such as X-ray pre-heating, recombination radiation, baryon-dark matter free-streaming, and alternative dark matter cosmologies. Our suite includes box sizes of $2$ and $0.25$ $h^{-1}$Mpc, extending to volumes large enough to begin capturing halos above the atomic cooling limit and resolutions high enough to fully resolve the IGM Jeans scale in the cold, neutral universe. We present a detailed description of the setup and first results from SAGUARO, descriptions of the IGM gas dynamics and thermal structure, opacity, self-shielding properties, the effect of the IGM on the reionization photon budget, and the halo mass function, and Ly$\alpha$ transmission properties. SAGUARO will help facilitate detailed studies of small-scale IGM structure and its effects that will help inform the next generation of reionization simulations and data interpretation.

We present map-making methodologies and preliminary spectral data cubes for SPHEREx, a NASA Explorer mission that launched in March 2025 and has been performing an all-sky near-infrared spectral survey. The SPHEREx instrument observes from 0.75 to 5.0 microns with a spectral resolution ranging from 35 to 130 and a pixel size of 6.15". We define a nominal set of 102 wavelength channels, each of which maps the entire sky approximately twice per year. Among the main mission goals is an investigation of the cosmic history of galaxy formation through intensity mapping of the extragalactic background light (EBL), which is a primary motivation for the map maker described in this work. The SPHEREx dataset contains a wealth of additional mapping targets, e.g., resolved galaxies and nebulae and diffuse clouds of Galactic dust and gas, which display strong spectral features such as hydrogen recombination lines, molecular-hydrogen lines and emission from polycyclic aromatic hydrocarbons (PAHs). We describe how our map maker handles these various cases, how to mitigate foregrounds such as zodiacal light and upper-atmospheric emission and how to monitor and mitigate systematics and signal loss. Our maps are produced both in tangent-plane projection and in full-sky HEALPix format. Specialized maps will be released to accompany future publications from the SPHEREx Science Team, and a public mosaic tool will be made available by the NASA/IPAC Infrared Science Archive (IRSA).

In this paper we present the public release of the Yale Rotating Evolution Code (YREC). YREC is a stellar evolution code that covers brown dwarfs and stars across a wide range of masses, and evolutionary states from the pre-MS through helium burning. We summarize the key ingredients of the code, document the code performance, and discuss its strengths and limitations. We present libraries of input files, documentation, sample use cases, and scripts. In addition to usage as a research tool, we highlight the utility of the code for educational purposes.

The LIGO-Virgo-KAGRA (LVK) candidate event S251112cm suggests the presence of at least one compact object with sub-solar masses. Since such objects cannot be produced through standard stellar evolution, this observation provides a potential indication of non-standard formation channels. Primordial black holes (PBHs), formed from the collapse of primordial density fluctuations in the early Universe, are a well-motivated candidate. We investigate the interpretation of S251112cm as the merger of two PBHs with masses in the range 0.1-1$\,M_\odot$. Combining analytic estimates of the PBH merger rate with current observational constraints on their abundance and the sensitivity of LVK searches, we compute the probability of observing such an event. Within a relaxed constraint scenario, the probability reaches unity in the range $M_{\rm PBH} \sim 0.5$-$1\,M_\odot$, while it remains sizable, $\sim \mathcal{O}(0.5)$, in more conservative scenarios and at lower masses. Our results show that a PBH interpretation of S251112cm is viable within current bounds. Owing to the dependence of our results on astrophysical uncertainties, such as those affecting the constraints on the abundance of PBHs, they cannot be regarded as conclusive with respect to the nature of the detected event. At the same time, our analysis highlights the potential of sub-solar gravitational wave events as a probe of PBHs and their contribution to dark matter.

The dynamics and observable signatures of neutron star mergers are governed by physics under the most extreme conditions. They are particularly impacted by the high-density equation of state, which for the most sophisticated models is usually available in the form of tables. Numerical relativity codes usually evolve particularly well-behaved numerical ("conservative") variables, but at the price that the physically interesting ("primitive") variables need to be found at every computational element and at every integration sub-step by means of expensive (and not always successful) root-finding algorithms. We have recently developed the Lagrangian numerical relativity code SPHINCS_BSSN which evolves the spacetime on an adaptive mesh with well tested methods, but the fluid is evolved by means of freely moving particles. Since our evolution equations differ from those of conventional numerical relativity, we need to develop new conservative-to-primitive algorithms if we want to use tabulated equations of state. We present here three such algorithms: a 3D and a 2D Newton-Raphson method and a 1D root-finding algorithm based on Ridders' method. We find the 3D method to be very fast and robust with an average failure fraction in a full-blown neutron star merger simulation (with the DD2 equation of state) well below 1%. While we do not find obvious advantages for the 2D method, the 1D Ridders' method is slow, but essentially fail-safe. Therefore, we choose the 3D Newton-Raphson as default and fall back to the 1D Ridders' method as a safe "parachute".

The sensitivity of radio telescopes is becoming increasingly limited by the presence of radio frequency interference (RFI), which will worsen as the radio spectrum becomes more crowded. One context where this poses a challenge is the field of fast radio burst (FRB) science, where there is increasing scientific interest in capturing as large of a population of bursts as possible and accurately measuring their celestial coordinates using interferometry. With several modern radio facilities actively collecting data for large FRB surveys that will be transformative to the field, properly mitigating unwanted interference is essential for the science goals of these surveys to be met. In this work, we present variations of a spatial filter based on the Karhunen-Loeve (KL) Transform to enhance the sensitivity of radio interferometers and demonstrate its applicability to FRB detection and localization. We derive a particular variation of the filter for the case of point-like radio pulses, which we show reduces to the maximum-signal-to-noise beamformer. We apply this filter to CHIME/FRB baseband data and demonstrate its capability to enhance the sensitivity and overall localization rate of CHIME/FRB Outriggers. We compare the cross-correlation signal-to-noise obtained using the spatial filter with that obtained using a spectral-kurtosis RFI flagger for a sample of 100 FRBs recorded by CHIME and its Outriggers, and show that this filter will double the total number of FRBs successfully localized with the CHIME/FRB Outrigger telescopes. While demonstrated here in the context of CHIME/FRB Outriggers, the spatial filter presented in this work--which we have made publicly available--is broadly applicable to other interferometric radio facilities engaged in FRB science and transient detection, including next-generation telescopes such as CHORD, DSA-2000, BURSTT, and CHARTS.

Magnetized and rotating core-collapse supernovae (CCSNe) are promising candidates for producing long gamma-ray bursts and hypernovae. In this project, we present 34 two-dimensional magnetized core-collapse supernova simulations with self-consistent neutrino transport, systematically exploring the parameter space of initial magnetic field strengths ($B_0 = 0$--$3.5 \times 10^{12}$~G) and rotation rates ($\Omega_0 = 0$--$0.5$~rad~s$^{-1}$) for a 40~$M_\odot$ progenitor. Our simulations reveal four distinct explosion morphologies: failed explosions leading to black hole formation, monopolar jet explosions, bipolar jet explosions, and neutrino-driven explosions. We find that the 40 $M_\odot$ progenitor model failed to explode without magnetic fields in two dimensions, even with rapid rotation. The non-rotating models require strong seed magnetic fields ($B_0 \gtrsim 1.5 \times 10^{12}$~G) to launch magnetically driven explosions, while the introduction of rotation substantially lowers this threshold. The explosion timescale decreases systematically with both increasing magnetic field strength and rotation rate, ranging from $>500$~ms in marginally successful models to $<150$~ms in strongly magnetized, rapidly rotating systems. Diagnostic explosion energies in the most extreme models approach $\sim 10^{51}$~erg within 250~ms and continue growing in time, making them potential hypernovae and long gamma-ray burst progenitors. Finally, we analyze the gravitational wave signatures associated with each morphology and find that the gravitational wave frequencies mainly depend on the rotation rates but are less sensitive to the magnetic field strengths and explosion morphologies. However, the gravitational wave amplitudes strongly depend on the explosion morphologies and magnetic fields, making searches for gravitational waves from magnetorotational core-collapse supernovae more challenging.

The Hubble Missing Globular Cluster Survey (MGCS) has taken one of the last opportunities to complete the census of Galactic globular clusters (GCs) started by past Hubble Space Telescope (HST) programs, securing high-resolution data for 34 GCs never observed before by HST. The previous papers of the series have highlighted the astrometric and photometric potential of the project by analyzing a subsample of targets. We present, and release to the community, the official astro-photometric catalogs of the MGCS for all GCs imaged by this project. We describe the data reduction using state-of-the-art techniques designed for HST. We discuss the photometric calibration and show, for the first time, the synergy with the Gaia catalog to ensure homogeneous photometry across our data set. We compute artificial-star tests that can be used to assess systematics and the completeness level of our data. We combined HST and Gaia data to refine the absolute PMs of our GCs, reaching a precision $\sim$3 times better than that of Gaia alone. We used these new PMs to update (and to determine for the first time for six systems) the associations between GCs and their putative galaxy progenitors. This work continues decades-long efforts of large Treasury programs in sharing precise and accurate atlases to the community for studying GCs across a wide range of scientific endeavors.

We investigate the redshift evolution of the matter density parameter, $\Omega_m(z)$, using galaxy cluster gas mass fraction measurements combined with cosmic chronometer $H(z)$ data and type Ia supernova luminosity distances. Our approach employs Gaussian Process Regression to reconstruct $\Omega_m(z)$ in a non-parametric way, remaining only weakly dependent on a specific background cosmology. The reconstructed evolution is consistent with the standard $\rho_m \propto (1+z)^3$ scaling predicted by the $\Lambda$CDM model. We obtain $\Omega_{m0}=0.296 \pm 0.044$ from the 44-cluster sample, and $\Omega_{m0}=0.271 \pm 0.016$, $0.253 \pm 0.017$, and $0.210 \pm 0.013$ for the 103-cluster compilation, depending on the assumed mass calibration. While $\Omega_m(z)$ follows the expected redshift behaviour, the inferred value of $\Omega_{m0}$ shows a strong dependence on the cluster mass calibration. Within this framework, mass bias emerges as the dominant source of uncertainty, exceeding statistical errors.

We propose a possible quantum signature of the early Universe that could lead to observational imprints of the quantum nature of the inflationary period. Graviton production from the presence of a classical, coherent state of the inflaton scalar field results in entangled states in the gravitons' polarizations. At horizon crossing, interactions between the gravitons and (lower scale) inflatons, together with the gathering of ``which-path information'' from the cosmological horizon, perform the required Bell experiments leading to a definitive measure, which can be imprinted in the scalar correlation four-point function. This is because of a non-trivial effect due to the derivatives on two scalar fluctuations, and it provides a fingerprint that depends on the polarization of the graviton that Alice and/or Bob measured in their patch. We hint how this signature could be measured in the high-order correlation function of galaxies, in particular on the halo bias and the intrinsic alignment.

We propose a method that provides an observational signature of the quantum origin of primordial fluctuations generated during inflation. The method gives a prescription for testing a Bell inequality constructed exclusively from the standard scalar and tensor perturbations of minimal single-field inflation. We consider an inflationary spacetime populated by pairs of gravitons entangled in their polarization states. Third-order interactions between two scalars and one graviton transfer polarization information to the scalar sector through the product of spatial derivatives of scalars with the tensor polarization factors. Rather than performing the full multidimensional momentum integrations, we isolate and compute the tensor polarization structure of the primordial scalar eight-point correlation function. This eight-point correlation function factorizes into the product of four scalar two-point functions associated with opposite (mirrored) momentum configurations in Fourier space. This factorization falls from the fact that the two gravitons are spatially well-separated within the cosmological horizon of inflation, replicating the setup of standard Bell experiments. Through these interactions, we track how non-local correlations between both gravitons from polarization entanglement are imprinted on the scalar sector. We show that, for specific configurations of the scalar momenta after the end of inflation (detailed in the text), this observable can be used to construct a Bell-violating quantity in a way that matches the well-known Clauser-Horne-Shimony-Holt inequality definition. In principle, this offers a route to probe the quantum nature of primordial fluctuations through observables accessible today.

Owing to their proximity to the Milky Way, the Large and Small Magellanic Clouds (L/SMC) uniquely probe the evolution of low-mass galaxies undergoing mutual interactions. In this work, we investigate the connection between the star formation histories (SFHs) of the L/SMC measured from HST imaging in the Scylla survey and APOGEE chemical abundances. We model the chemical evolution of the L/SMC in the [Mg/Fe]-[Fe/H] plane within a robust statistical framework to predict chemical abundance signatures resulting directly from starbursts in Scylla SFHs. Both the L/SMC rapidly enrich to high metallicity ([Fe/H] $\gtrsim$ $-1$) within 3 Gyr, followed by slower chemical evolution regulated by sequential starbursts, where the SMC may require higher Fe yields from Type Ia supernovae than the LMC. We also model the [Mg/Fe]-[Fe/H] plane to infer starburst properties across distinct spatial regions in the L/SMC. We identify dominant starbursts in the L/SMC with broadly similar timing, though the SMC may host an earlier burst, and larger burst strength in the LMC. The global starburst properties are nearly uniform across the LMC disk, whereas the dominant SMC population experiences a stronger and later-onset burst in its eastern wing compared to the main body. We also find evidence for a chemically distinct secondary population in the SMC that preferentially traces the foreground and may originate from the LMC. We discuss the implications of these results for the evolutionary history of the L/SMC and for starbursts in interacting low-mass galaxy pairs.

In this paper we present a newly detected companion to the Be star, HD 52244 (B2IVnpe), using the Fine Guidance Sensors (FGSs) on the Hubble Space Telescope (HST). In fall 2021, HST became momentarily unavailable to support nominal operations, and we used the operational FGS to carry out a multiplicity survey of 6 Be stars. We were able to resolve a companion to HD 52244, with a separation of 42.7 +/- 1.1 mas (74 AU) and a position angle of 144.2 +/- 0.3 with a differential magnitude in the F583W filter of 1.91 +/- 0.02 mag. This study presents the results to the newly detected companion of HD 52244 and lays the groundwork for future studies looking for wide or third companions to Be stars.

We discuss the results of using large language models (LLMs) to conduct original scientific research in an unfamiliar subject area during the Fall 2025 semester. Students in a graduate astronomy and astrophysics course were asked to test whether LLMs could help them complete research tasks faster and at a level of detail and accuracy required for scientific publication. Most students employed LLMs for a total of 5-10 hours. While all students completed a draft paper on an unsolved problem related to galaxies by semester's end, their impressions of the models' value varied. About half thought that the models saved them time. Many noted that LLMs failed to provide appropriately detailed insights or steps to addressing open, niche questions over a several-month timeframe. The LLMs also frequently (about 20% of the time) returned false citations, links, or summaries of papers. The models struggled with generating complex functional code, accessing online packages or Application Programming Interfaces (APIs), and retrieving astronomical datasets from existing archives. In writing code and in chats, the LLMs made implicit, overly simplifying assumptions and often doubled down even after being corrected. Given the rapid pace of LLM development, new models may soon address at least some of these issues and thus significantly enhance research productivity. Yet students expressed concerns about how LLM use might dampen creativity and reflection during the research process. To improve learning experiences in future semesters, the class will first discuss LLM best practices and limitations. Students will be encouraged to explore free online resources for tips for generative model applications and will decide for themselves whether to use LLMs for their research project. This white paper was not written using LLMs.

Blue straggler stars (BSSs) are believed to form through mass transfer in binary systems or stellar collisions. The reported presence of double BSS sequences in some globular clusters (GCs) has been interpreted as evidence that these two formation channels produce distinct sequences in color-magnitude diagram (CMD). We reassess this claim using HST UV Globular Cluster Survey (HUGS) photometry of 56 Galactic GCs. We used the Hartigan Dip Test to test bimodality, and Akaike model comparison to test whether BSS distance distributions are better described by a mixture of two unskewed Gaussians or a skewed unimodal Gaussian model. We find no strong statistical evidence for bimodality; no cluster yields a dip test p-value below 0.15, and Akaike model comparison favors the skewed unimodal model in 94 out of 112 cases. We re-examine NGC 7099 (M30), the prototypical case of a double BSS sequence, using three reductions of HST data. We find bimodality is detected at p = 4 x 10-3 , vs. the originally reported p ~ 10-5 , in the original photometry. The observed uncertainties derived from the subgiant branch widths are comparable to the suggested separation between the proposed BSS sequences, making the detection of statistically significant bimodality challenging. Our results suggest that the dip between two BSS sequences in M30 photometry is a coincidence, and that later bifurcation claims can be explained as skew in the BSS color distribution, rather than two separate distributions.

An evolving turbulent flow such as the solar wind can be meaningfully characterized by its "turbulence age" -- an estimate of the number of nonlinear times that have elapsed during a plasma parcel's propagation from the Sun to a given point in space. Recent observations of the near-Sun solar wind by the \textit{Parker Solar Probe} (\textit{PSP}) indicate high correlation between velocity and magnetic fluctuations (i.e., cross helicity, $\sigma_c$), which is known to impede development of magnetohydrodynamic (MHD) turbulence. Here we propose a new formulation of the turbulence age ($A_\text{t}$) of the solar wind that explicitly accounts for the Alfvénic nature of the fluctuations in the inner heliosphere. $A_\text{t}$ is then evaluated for slow and fast wind streams using a variety of data sources -- observations from the \textit{PSP, Advanced Composition Explorer}, and \textit{Voyager} missions, and a global solar wind simulation that includes turbulence transport. Compared to the formulation employed in previous work that neglected Alfvénicity, the present approach yields smaller values of $A_\text{t}$ in medium-to-high $\sigma_c$ solar wind; similar turbulence ages are then obtained for slow and fast wind in the ecliptic. The radial evolution of $A_\text{t}$ between heliocentric distances of $r\sim 0.2$ to 40 AU is examined. The rate of increase of $A_\text{t}$ is found to decrease until $\sim 5$ AU, indicating a gradual slowing of the \textit{in situ} development of turbulence in the inner heliosphere. Beyond $\sim 5$ AU this rate begins to increase, likely due to turbulence driving by pick-up ions. This paper highlights the important role of cross helicity in modulating MHD turbulence, and the results will aid in further interpretations of observations of the radial evolution of various turbulence parameters in the solar wind.

The big bang quantum cosmology model introduces the trace $J$ of the Schouten tensor as a form of dynamic dark energy. Together with cold dark matter, these components form the so-called $J$CDM cosmology model, proposed by M.H.P.M. van Putten (J. High Energy Astrophys., 45, 2025, 194), which offers a potential resolution to the Hubble tension. We derive the constraints on the $J$CDM cosmology model, utilizing early- and late-time cosmological data including cosmic microwave background (CMB), baryon acoustic oscillations (BAO) released by the Dark Energy Spectroscopic Instrument (DESI), cosmic chronometers (CC), and type Ia supernovae (SNIa). For a flat universe, the $J$CDM model yields \( H_0 = 66.95 \pm 0.51 \, \rm{km~s^{-1}~Mpc^{-1}} \) and \( \Omega_m = 0.3419 \pm 0.0065 \), results that are consistent with early-universe observations but exhibit a higher \( \Omega_m \) compared to the $\Lambda$CDM model. In the case of a non-flat universe, $J$CDM favors a slightly curved geometry with \( \Omega_k = 0.0154 \pm 0.0027 \), leading to \( H_0 = 69.13 \pm 0.56 \, \rm {km~s^{-1}~Mpc^{-1}} \) and \( \Omega_m = 0.3477 \pm 0.0074 \). The increase in \( H_0 \) in the non-flat scenario suggests a geometric degeneracy between spatial curvature and \( H_0 \). We also investigate the internal inconsistencies present in DESI data and evaluate their impacts on cosmological parameter constraints. Our analysis shows that while the $J$CDM model, which is constructed from first principles without free parameters beyond those of $\Lambda$CDM, agrees excellently with late-time cosmology, it struggles to simultaneously match early-universe observations in a fully self-consistent manner.

Although stellar bars are prevalent in local galaxies, unbarred galaxies constitute a significant fraction, particularly at high redshifts. While some galaxies are unbarred by nature due to stability against the bar instability, several mechanisms capable of transforming barred galaxies into unbarred systems have also been proposed, such as central mass concentration, specific dark matter halo morphologies and tidal interactions. Regarding galactic interactions, mergers can undoubtedly disrupt bars while potentially destroying the entire disk. However, the effects of pure collisions (non-merging scenarios) on bars remain unclear, with limited existing studies yielding contradictory conclusions. Here we aim to systematically investigate the disruptive effects of collisions on bars hosted by Milky Way-like galaxies using N-body/SPH simulations. We model collisions between the barred galaxy and a spherical intruder, conducting multiple simulations by varying interaction parameters, with mass ratios set at 1:3, 1:5, and 1:15. We find that bars are remarkably robust, with most interactions failing to significantly reduce their strength or pattern speed. Only off-center high-inclination retrograde collisions can effectively destroy bars, while central high-inclination collisions can substantially decrease the pattern speed. Such destruction and deceleration primarily arise from gravitational forces rather than gas-related processes. Notably, compared to collisions occurring along the bar's major axis, those along the minor axis cause greater weakening but can slow the bar's natural deceleration. Furthermore, changes in mass resolution do not significantly affect the results when the resolution is better than ~10^5 Solar mass.

Solar radio bursts (SRBs) are important signatures of dynamic processes in the solar corona, including particle acceleration and shock propagation associated with solar flares and coronal mass ejections. Among the missions that report solar radio bursts within 24 hours, the e-CALLISTO archive is the largest and most successful, with more than 150 stations worldwide. The archive generates large volumes of FITS data that are often affected by radio-frequency interference and background noise. Irregular frequency setups in different stations are also a limitation of statistical analysis of SRBs. Each CALLISTO observation is a 15-minute frame, which often causes a single burst to split over multiple frames, making event-level analysis difficult. This work presents the e-CALLISTO FITS Analyzer, a unified, interactive, cross-platform application for processing and analyzing e-CALLISTO dynamic spectra on Windows, macOS, and Linux. The application supports time and frequency merging to produce a continuous spectrum, applies mean background subtraction with user-controlled threshold clipping, and isolates burst regions through an interactive polygon mask in the time--frequency plane. It also extracts the maximum-intensity backbone, allows interactive outlier removal, and performs power-law fitting to estimate drift rates and derive shock height and speed using the Newkirk model, including $n$-fold scaling. For a Type II burst observed by Arecibo Observatory on 2 March 2022, the analyzer yielded an average drift rate of $-0.0400 \pm 0.0003\, MHz/s$ and an average shock speed of $449 \pm 1 \, km/s$ at a height of $1.715 \pm 0.002\, R_{\odot}$. The e-CALLISTO FITS Analyzer supports more reproducible, event-focused SRB analysis and improves access to physically meaningful measurements from e-CALLISTO FITS data.

Evolutionary tracks for 2-8M$_\odot$ models, covering a [Fe/H]=$-$1.0 ($Z=0.0014$) to [Fe/H]=+0.2 ($Z=0.02$) metallicity range are computed with Modules for Experiments in Stellar Astrophysics, MESA, to investigate evolutionary and pulsation properties of classical, fundamental mode Cepheids. We examine in detail the effects of convective overshooting from the Main Sequence core, as well as from the convective envelope on the Red Giant Branch. Mass loss is also included in a few model sets. Linear pulsation properties are derived consistently with a module of MESA, Radial Stellar Pulsation, RSP. We provide edges of the classical Instability Strip, as well as ages, crossing times through the Instability Strip and period change rates. Period-Luminosity, Mass-Luminosity, Period-Radius and Period-Age relations are provided, both in analytical and tabular form. Their dependence on metallicity, crossing number and overshooting parameters are investigated. Qualitative comparisons with classical Cepheids in the Milky Way and Magellanic Clouds as well as other theoretical relations are presented. We find satisfactory agreement for most of the observables and good match with other theoretical work, however reproducing short-period Cepheids in the Small Magellanic Cloud as well as Cepheid mass discrepancy pose a challenge for the presented models. Considering metallicity effect of the Period-Luminosity relation, we find $\gamma\approx -0.20$ mag dex$^{-1}$, nearly independent on photometric pass band and in good agreement with recent observational studies. The magnitude of this effect depends on the underlying mass-luminosity relation, being stronger for relations that predict higher luminosities at a given mass.

The electric fields inside thunderstorms can significantly modify the intensity of secondary cosmic ray muons at the ground level, producing measurable variations in their intensity ($\Delta$I$_{\mu}$). By utilizing the decade-long observations of thunderstorms (April 2011-December 2020) by the GRAPES-3 muon telescope (G3MT), a directional asymmetry in $\Delta$I$_{\mu}$ is observed, with nearly six times more events being detected from the east than the west directions. Using detailed CORSIKA Monte Carlo simulations, it is shown that this asymmetry is caused by the variations of the muon charge ratio R$_\mu$ (N$_{\mu^+}$/N$_{\mu^-}$). The anisotropic R$_\mu$ in turn, is caused by the systematic changes in geomagnetic cutoff rigidities, and subsequent selective filtering of predominantly positively charged primary cosmic rays. As a consequence, the R$_\mu$ increases systematically from west to east across the G3MT field of view, enhancing the sensitivity of east directions to positively charged thunderstorm top. Monte Carlo simulations with constant R$_\mu$ show that the directional asymmetry disappears, demonstrating the muon charge imbalance to be the dominant driver of the observed asymmetry. The dependence of R$_\mu$ on the hadronic interaction is also studied by comparing seven combinations high-, and low-energy hadronic interaction generators, which show a $\lesssim$7% spread in R$_\mu$, and $\lesssim$14% variation in the derived thunderstorm potentials. These results provide the first quantitative link between the muon charge asymmetry caused by the geomagnetic field, and the directional distribution of thunderstorms, reinforcing the role of muon observations as a probe of gigavolt potentials in atmospheric electrical structures.

Following long periods of quiescence, low-mass X-ray binaries can exhibit intense X-ray outbursts triggered by instabilities within the accretion disk. These outbursts can sometimes be detected in optical wavelengths before being detected in X-rays, acting as an early onset warning and enabling a deep study of accretion disk properties informed by the lag between optical and X-ray rise. We explore the potential of Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) to detect these outbursts early through optical observations. We evaluate the capabilities of LSST based on currently planned survey cadence, filter-specific depth, and other observational factors that affect early detection. We develop and apply an extended metric to assess outburst detectability and recovery fraction. We find that despite inherent potential for early detection of XRB outbursts, the currently planned survey strategy makes it challenging to detect early onset of XRBs. Lastly, we demonstrate how this estimate can be used to infer the wider LMXB population in the Galaxy as the LSST progresses.

Angular momentum removal is a fundamental requirement for star and planet formation, yet the mechanisms driving this process remain debated. Magnetohydrodynamic disk winds, launched along magnetic field lines from extended disk regions, offer a promising solution, particularly in regions where magnetorotational turbulence is weak. Here we present high-resolution Atacama Large Millimeter/submillimeter Array observations of the Class 0 protostar HOPS 358, revealing a rotating, nested outflow structure traced by H2CO, SO, and CH3OH emission. The outflow preserves the disk's rotational sense and is aligned with the disk axis, providing direct observational evidence for a magnetically launched disk wind. From the measured kinematics, we derive a dimensionless magnetic lever arm of approximately 2.3 and constrain the wind-launching region to radii of 10-18 astronomical units within the planet-forming zone. These results demonstrate that magnetohydrodynamic disk winds operate during the deeply embedded phase, efficiently extracting angular momentum while shaping disk evolution and establishing initial conditions for planet formation.

Stellar-mass black holes (BHs) surrounded by neutrino-dominated accretion flows (NDAFs) are a leading central engine of gamma-ray bursts (GRBs). In this work, we investigate the electron fraction distribution in NDAFs with or without disk outflows for different accretion rates, BH spins, and outflow rates. As the results, for the cases of the massive disks at relatively low accretion rates, the outer boundary of the disks are predominantly advection-cooled, yielding electron fractions of \(Y_{\rm e} \sim 0.5\), as expected for massive collapsar progenitors. By contrast, in the cases of lower-mass disk at high accretion rates, neutrino cooling becomes highly efficient and mildly electron-degenerate disks emerge, characterized by \(Y_{\rm e} \lesssim 0.38\) at the outer boundary of the disk, even for the strong outflows, which is consistent with materials from compact object merger scenarios. Moreover, we find that these trends remain robust across different BH spins. Consequently, the self-consistent agreement between the electron fraction properties at the outer boundaries of NDAFs and those expected from GRB progenitors provides effectively support for NDAFs serving as the GRB central engines.

Using two observations obtained with the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we present a detailed single-pulse analysis of the high-nulling pulsar PSR J1820-0509. We measure an exceptionally high nulling fraction of approximately 81.78%, significantly exceeding previous estimates from Parkes observations. The single-pulse energy distribution exhibits a clear bimodal structure, consistent with classical nulling behavior. However, stacking the identified null pulses reveals a statistically significant residual profile above the noise level, indicating that the nulls correspond to a very weak emission state rather than a complete cessation of radio emission. The pulsar shows clustered burst activities spanning several hundred rotation periods, with prominent quasi-periodicities at 1191 +/- 81 and 590 +/- 15 pulse periods in the two observations. Based on temporal clustering and integrated profile morphology, we identify three distinct emission modes (A, B, and C) and a pseudo-null state (D). These modes exhibit systematic differences in pulse morphology, polarization, and energy statistics. The pulse width-energy relations reveal clear transitions between low- and high-energy regimes. The energy distributions of Modes A and C are well described by lognormal functions, while Mode B follows a composite Gaussian-lognormal distribution. These results suggest that the radio emission of PSR J1820-0509 is governed by multiple quasi-stable magnetospheric states. The presence of weak emission during pseudo-nulls, together with systematic mode-dependent variations, supports the interpretation that pulsar nulling reflects transitions between different magnetospheric activity levels rather than a complete shutdown of emission.

Internal rotation has emerged as a fundamental feature of globular clusters (GCs), yet its origin and long-term evolution remain poorly understood. We explore the evolution of rotating GCs over a Hubble time under the combined influence of two-body relaxation, tidal field, and stellar evolution. We introduce the ROLLIN' simulations, a suite of 25 N-body models characterized by a realistic number of stars from 250k to 1.5M, ran with the direct N-body code NBODY6++GPU and evolved for 14 Gyr. With present-day masses of 5 x 10^4 - 5x10^5 M_sun, the models cover the parameter space of low-density MW GCs. Our analysis reveals that rapidly rotating GCs experience earlier and more pronounced core collapse, efficiently segregating massive objects and remnants in their centers within the first few 100 Myr. In the long-term, internal rotation declines and a correlation emerges between rotation and GC mass, in agreement with observations. The primary driver of this evolution is mass loss, capturing both internal (stellar evolution, evaporation) and external processes (tidal stripping). The velocity anisotropy also evolves in response to mass loss: GCs initially near isotropy develop radial anisotropy, peaking around 40% mass loss, before progressing toward isotropy or tangentiality. The GC orbital history also plays a role, as retrograde rotators retain rotation more effectively than prograde rotators. Finally, we quantify the long-term changes of GCs after 12 Gyr: (1) The surface density decreases by up to 2 orders of magnitude. (2) The half-mass radius increases by a factor of 3-5. (3) The rotation decreases by a factor >5 for GCs that have lost >50% of their mass. The ROLLIN' simulations demonstrate that angular momentum is crucial to understand the origin, evolution, and survival of GCs. These models provide a benchmark for interpreting GC observations in the local and high-z Universe.

Recently, binary pulsar timing has been proposed as a viable probe of dark matter subhalos with masses of $\sim 10^7\,M_{\odot}$ in the solar neighborhood. We present a comprehensive analytical framework that incorporates the subhalo mass function, projection effects of line-of-sight acceleration, and the spatiotemporal geometric requirements for joint detection by binary systems, enabling a quantitative evaluation of the detectability of nearby subhalos. Applying this framework to the current binary pulsar sample, we find a probability $\leq 1.6 \times 10^{-4}$ of detecting at least one subhalo within the effective volume. An independent timing residual analysis shows no statistically significant excess in line-of-sight accelerations beyond predictions from data-driven Galactic gravitational potential models. These results place stringent constraints on detecting $<10^8~M_{\odot}$ dark matter subhalos with existing pulsar timing data, aligning with the theoretical expectation that such subhalos have a low survival probability in the solar neighborhood. A low detection prospect still holds even for future Square Kilometre Array observations.

We analyse the sizes and specific angular momentum content of galaxies in the Colibre cosmological hydrodynamical simulations spanning two orders of magnitude in mass resolution. We compare the predicted size-mass and angular momentum-mass relations to a broad range of observational measurements spanning redshifts $z=0$ to $4$. At $z=0$, Colibre reproduces observed size-mass relations over the sampled mass range $10^8 \lesssim M_\star/{\rm M_\odot}\lesssim 10^{11.5}$, and for multiple size definitions, including two- and three-dimensional stellar half-mass radii, half-light radii across several wavelengths, as well as alternative measures such as baryonic half-mass radii and characteristic radii defined by stellar surface density thresholds. The simulations also recover the observed segregation of galaxies in the size-mass plane by morphological type and star formation rate, and reproduce the distinct, approximately parallel sequences followed by star-forming discs and quenched spheroids in the stellar specific angular momentum-mass plane. The angular momentum content of star-forming Colibre galaxies match that of observed systems out to $z\approx 1.5$. At higher redshifts, massive galaxies ($ 10^{9.5}\lesssim M_\star/{\rm M_\odot}\lesssim 10^{11}$) in the simulations are somewhat smaller than observed, and the separation between star-forming and passive populations in the size-mass plane is reduced relative to observations, while at lower masses the agreement remains good. This apparent discrepancy may reflect the effects of dust attenuation, which is neglected in our analysis and may preferentially obscure the central regions of observed systems. Overall, our findings highlight the close connection between galaxy size, angular momentum, and morphology over cosmic time.

Close-in rocky planets are the most common type of exoplanets around late M dwarfs, ranging from more temperate worlds to highly irradiated lava planets with molten surfaces, and many theoretical studies have attempted to explain their formation. However, the origin of rocky planets with orbital periods shorter than one day, known as ultra-short-period (USP) planets, remains uncertain. We aim to investigate whether the formation and survival of USP planets is connected to the location of the inner edge of the protoplanetary disk, considering different disk edge prescriptions. We use N-body simulations that include planet-disk interactions, star-planet tidal interactions, and relativistic corrections, applied to a sample of lunar-mass planetary seeds growing via pebble accretion in a low-viscosity disk ($\alpha_t = 10^{-4}$). The inner edge of the disk is modeled in three ways: as a fixed close-in edge, as an outward-evolving edge set by the magnetospheric truncation radius, and as an inward-evolving edge defined by the corotation radius. USP planet formation appears to be tightly controlled by the location of the disk's inner edge. Our simulations show that only the close-in-fixed-edge Scenario and the inward-evolving-edge Scenario are capable of producing USP planets, as planets tend to follow the movement of the disk's inner edge. This suggests that USP planet formation is favored when the inner edge remains close to the corotation radius of a rapidly rotating star.

Context: A statistically significant sampling of H$_2$ rotational excitation in the diffuse interstellar medium (ISM) is essential to identifying its excitation mechanisms and assessing the importance of H$_2$ in the cooling of the gas and the regulation of thermal pressure. Aims: To complement the statistics provided by ancillary telescopes, we conducted a search for pure rotational H$_2$ emission lines in all publicly available background observations obtained with the Medium Resolution Spectrometer (MRS) aboard the JWST. Methods: The sample consists of 276 background observations acquired over the past three years. Departing from the standard pipeline, each uncalibrated MRS background file was reprocessed, enabling the analysis of H$_2$ pure rotational emission. Lines of sight likely associated with star-forming complexes were excluded to focus on emission from the diffuse ISM. The results were compared with FUSE absorption data and were analyzed in relation to the column densities of H and H$_2$ and to dust emission derived from HI4PI, Planck, and WISE data. Results: This analysis reveals widespread H$_2$ emission throughout the Galaxy. We report the first detections of the pure rotational S(4), S(5), and S(7) lines in the diffuse ISM. The S(1) line is detected along 84 lines of sight, corresponding to a detection rate of 41%. Its integrated intensity decreases steeply with Galactic latitude, spanning nearly two orders of magnitude, in remarkable agreement with absorption measurements. The $T_{34}$ and $T_{35}$ excitation temperatures vary between 200 and $\sim$1000 K, are correlated with each other, and are anticorrelated with the column density of H$_2$ , as expected from ancillary data. All lines of sight in the sample have undergone the H-H$_2$ transition, at $N_{\rm{H}} \gtrsim 10^{20} \ \rm{cm}^{-2}$, and are partly molecular, with $f_{\rm H_2} \gtrsim 0.1$. Under these conditions, the cooling rate associated with the S(1) line, expressed per hydrogen atom, is found to be remarkably constant, with a characteristic value of $\sim 4\times10^{-27}$ erg s$^{-1}$ H$^{-1}$. Conclusions: This study demonstrates that the high sensitivity of the JWST enables measurements that both strengthen and complement those from absorption studies. Observations collected over just a fraction of JWST's lifetime have already yielded detections along dozens of lines of sight, significantly expanding the statistical sample of H$_2$ rotational excitation in the diffuse ISM.

Relativistic jets are among the most fascinating objects in the Universe, and recent high-resolution Very Long Baseline Interferometric (VLBI) observations, including the Global mm-VLBI Array and the Event Horizon Telescope (EHT), are able to resolve their structure close to their launching site. These observations reveal strongly limb-brightened jet structures for Centaurus A (Cen A), M 87 and 3C 84. Thus, the question arises which physical mechanism can generate the limb-brightened structure, and if this structure is common for jets from low-luminosity active galactic nuclei (LLAGN) seen under large viewing angles. Therefore, as a pilot study, we aim to model the EHT observations of Cen A. We performed a 3D two-temperature general-relativistic magnetohydrodynamic (GRMHD) simulation of accreting supermassive black holes (SMBHs) and jet launching to study the plasma dynamics and computed the connected emission via general relativistic radiative transfer (GRRT) calculations considering possible anisotropies in the distribution of the radiating particles. In order to adjust our simulations to the EHT observations of Cen A, we carried out a Bayesian fitting in the Fourier plane. We find that GRMHD simulations of magnetically arrested disks (MADs) combined with anisotropically emitting particle distributions along the direction of the magnetic field, parametrized by a value {\eta}=0.07, are able to mimic the recent EHT observations of Cen A. In addition, we extracted a black hole mass of {M_\mathrm{BH} = 6\times10^7 M_\odot} and a viewing angle of {\vartheta}=72°. Our obtained model can reproduce key features of the EHT and Atacama Large Millimeter/submillimeter Array (ALMA) observations in total and polarized emission. Finally, we predict that the black hole shadow in Cen A will be observable at a frequency of $\sim$ 3 THz.

The impulsive phase of a solar flare is known to generate strong turbulence and to transfer magnetic energy into accelerated electrons. Recognizing the importance of angular diffusion on the dynamics of the accelerated electrons, we extend previous treatments by deriving analytic solutions for the electron flux and associated energy deposition in two regimes: scattering dominated by inelastic Coulomb collisions and scattering dominated by elastic interactions with turbulent scattering centers. We show that the turbulence-dominated scattering term strongly reshapes the spatial distribution of the plasma heating: compared to the traditional collisional thick-target approach, turbulent scattering could lead to an order-of-magnitude increase in coronal heating and an even greater suppression of chromospheric heating. Scattering also acts to reduce the anisotropy of the electron distribution and so reduces the net current associated with the nonthermal electrons. The return-current Ohmic heating is accordingly reduced to a level that renders it negligible compared to direct collisional heating. The results have significant implications for models of atmospheric response to impulsive phase energy release, in particular chromospheric evaporation, flare-driven coronal heating, the formation of loop-top hard X-ray sources, and the longstanding discrepancy between modeled and observed soft X-ray line profiles.

Recently, the GeV--sub-PeV spectrum of supernova remnant (SNR) Cassiopeia A (Cas A), one of the youngest and most well-studied SNRs in our Galaxy, has been updated by observations of Fermi-LAT and LHAASO. We revisit Cas A with our previous shell-plus-jet asymmetric model and investigate its particle acceleration ability. The broadband fitting results suggest that the double-peaked gamma-ray spectrum can be well attributed to proton-proton (PP) collisions and inverse Compton scattering within the SNR shell, while the synchrotron emission from a jet component with velocity of $\sim0.1c$ can account for the hard X-ray emission up to 220 keV. Furthermore, the PP collisions in the jet can produce a sub-PeV emission, but constrained by the LHAASO-KM2A limit to a flux below $\sim 1\times10^{-14}\rm erg/(cm^2s)$ at 100 TeV. The energy of accelerated protons in the jet of Cas A could be up to $5\times10^{47}$ erg, which, assuming that the PeV cosmic ray distribution is clumpy in the Galaxy with the clump size comparable to the thickness of the Galactic plane, derives a proton flux consistent with the observed one at 1 PeV, implying that the Cas A-like SNRs can still be PeVatrons in the Galaxy. It is encouraging for LHAASO and future telescopes to detect or constrain Cas A spectrum above 100-TeV more precisely.

In this paper, part of a series on the effects of X-ray sources in promoting Population III (Pop III) star formation, we investigate the ionisation and heating of the intergalactic medium (IGM) and the consequent enhancement of molecular hydrogen (H$_{2}$) and Pop III formation using cosmological zoom-in simulations. We adopt a minimal X-ray feedback model in which X-rays originate solely from Pop III supernovae, and compute the global X-ray and Lyman-Werner (LW) radiation backgrounds on-the-fly during the simulation of a mean-density region of the Universe. This approach self-consistently captures the feedback loop between Pop III stars and the radiation backgrounds they produce. Pop III supernovae generate a weak X-ray background (J$_{\mathrm{X,21}} \sim 10^{-5}$) and a moderate LW background (J$_{\mathrm{LW,21}} \sim 10^{-1}$); the latter intensifies below $z \approx 12$ (J$_{\mathrm{LW,21}} \sim 10^{1}-10^{2}$) with the onset of Pop II star formation. Applying these backgrounds to regions of varying mean density produces a net positive X-ray feedback that increases the Pop III number density, with stronger enhancement in underdense regions. The positive feedback is more pronounced when the X-ray background is computed on-the-fly rather than by post-processing, demonstrating the importance of the feedback loop. The X-ray background also raises the Thomson scattering optical depth at high redshift, while the total optical depth remains consistent with Planck 2018 constraints. Because our model includes only Pop III supernovae as X-ray sources, it represents the most conservative scenario; stronger X-ray feedback is expected when additional sources are included, as will be explored in future work.

An X-ray background in the early Universe enhances molecular hydrogen formation, the main coolant of primordial gas, thereby lowering the threshold for Pop III star formation. Continuing our series on X-ray impacts on Pop III star formation, we investigate how a soft X-ray background promotes Pop III star formation using cosmological zoom-in simulations of ten cosmic volumes spanning a range of halo number densities. Each volume is irradiated by the Lyman-Warner (LW) H$_{2}$ dissociating background and a weak (J$_{21} \sim 10^{-5}$), soft ($E \sim 0.2-2.0$ keV) X-ray background produced by pair-instability SNe (PISNe) from Pop III stars and calculated self-consistently as described in a companion paper. We also compare the same simulations with and without X-rays to isolate the X-ray effect. The background promotes Pop III star formation in two ways: (1) by reducing the mean host halo mass by a factor of $\sim 2-3$, and (2) by enabling Pop III star formation in haloes that would otherwise remain sterile, thereby increasing the halo occupation fraction. The resulting gain in Pop III number density is largest in underdense regions (a factor of $\approx 3$ on average, reaching up to 7). In the most extreme case, Pop II stars form only in the presence of X-rays and the gas-phase metallicity rises by an order of magnitude, suggesting that dwarf galaxies in underdense regions may be significantly influenced by an early X-ray background. We also provide fitting functions for the halo occupation probability of Pop III stars as a function of redshift for both X-ray and LW-only simulations, which can serve as inputs for semi-analytic models.

We analyse the gravitational-wave emission from 60 two-dimensional core-collapse supernova simulations. The models cover a range of progenitors and equations of state. We focus on the narrow frequency interval in the gravitational-wave spectrum where the emitted power is strongly suppressed (the power gap) and how its central frequency relates to the physical properties of the simulations. We find that the power-gap frequency exhibits strong and systematic correlations with the properties of the inner core of the forming neutron star, for example the sound speed, suggesting that the gap encodes information about the behaviour of matter at extreme densities. We further examine how well several mechanisms proposed in the literature account for the presence and evolution of the gap in our simulations. Finally, we explore a scenario in which the gap arises from destructive interference between a narrow oscillation mode and a broadband background signal, demonstrating that such an interaction can produce a sharp minimum in the emitted gravitational-wave power.

Detecting faint objects in cislunar space using ground-based optical telescopes is difficult because of their low brightness, strong lunar background, and complex, nonlinear apparent motion. Traditional shift-and-stack techniques based on linear motion assumption suffer signal trailing loss due to significant nonlinear motion during long integrations, thus producing a degraded signal-to-noise ratio (SNR). In this paper, we first derive a theoretical criterion based on the point spread function to determine the maximum applicable integration time for linear-motion stacking. We then propose a quadratic shift-and-stack (QSS) method to correct for the first-order nonlinear motion, namely the angular acceleration of cislunar targets. Simulations of typical cislunar orbits verify this theoretical criterion and show that the QSS method significantly improves SNR from stacking and can enhance the detection limit by up to 1 stellar magnitude compared with the linear-motion stacking method. Furthermore, tests using observational data of the cislunar object Tiandu-1 confirm that while linear stacking degrades after a 29-minute integration due to trajectory curvature, the QSS method achieves continuous SNR improvement over a 46-minute integration, outperforming the peak SNR of the linear method by 31%.

In this paper, a new tool for Doppler tomography, Tomo-V (this https URL) that is developed based on the algebraic reconstruction technique (ART) has been presented. Previously, the ART method has not been widely used in tomography, as its direct implementation was computationally complex. The author has developed a fast version of this algorithm, which allowed it to be implemented within a web application that runs at acceptable speed in a browser on a personal computer. This method can be used to obtain sharp tomographic images from blurred profiles. Furthermore, the method has demonstrated excellent results in reconstructing images from noisy data, from a small number of profiles, and from profiles contaminated by absorption lines and emission from the expanding envelope. Tomo-V also includes tools for analyzing the resulting tomograms, allowing the position of accretion disks and Roche lobes to be displayed on the tomogram, as well as back-projecting the tomographic image onto flow elements in spatial coordinates. The paper is partially based on a report presented at the Modern Stellar Astronomy 2025 conference.

Black holes (BHs) and wormholes (WHs) are characterized by distinct spacetime geometries, whose differences become pronounced close to the central objects. A useful way to probe such differences is via the dynamics of stellar tidal disruption events in the regime of strong gravity. Here, using a general relativistic smoothed particle hydrodynamics code inspired from an algorithm developed by Liptai and Price, we perform a suite of numerical simulations of solar mass polytropic stars in the background of supermassive Schwarzschild BHs and similar mass exponential WHs. Important differences between the two geometries near the BH event horizon or the WH throat is provided by the distinct outcomes of such events. For a given impact parameter, BH backgrounds lead to greater tidal stripping compared to WHs ones and further, the critical impact parameter, beyond which the star undergoes full tidal disruption is higher for WH backgrounds compared to BHs. We further study the differences in observable peak fallback rates in the two backgrounds. We also provide a quantitative explanation for the tendency of stars in partial tidal disruptions to retain larger cores around more massive centers, by computing tidal stresses in a Fermi normal coordinate system and introducing an appropriate measure of stellar compactness. Finally, we suggest a way to observationally distinguish BH and WH backgrounds, based on the properties of different observables.

We present the first spectral analysis of Scorpius X-1 (Sco X-1) using intentional stray light (SL) observations taken by NuSTAR. Unlike focused observations that have high telemetry load when observing bright sources, intentional SL observations can help reduce the telemetry and reduce the effect of dead time, thereby maximizing the on-source exposure time; all of which are critical for extremely bright sources that exhibit short timescale variability like Sco X-1. The intentional SL observation of Sco X-1, taken in 2023, captured the source primarily in the flaring branch (FB) of the Z track. We performed spectral modeling of the continuum and reprocessed emission. A combination of thermal and Comptonization components (modeled with thcomp) provided a robust fit to the continuum. We test both scenarios for Comptonized emitting regions arising from the accretion disk and close to the neutron star, which provides comparable fit statistics. Reflection was modeled with the relxillNS model, enabling measurements of disk inclination consistent with prior radio and IXPE studies and comparison of inner disk radius to the emission radii of the thermal components. Overall, the results from the intentional SL data provide comparable results to literature on the focused NuSTAR data of Sco X-1 in the FB or taken contemporaneously. The success of this observation demonstrates the capability of SL data to yield high-quality spectral constraints comparable to focused observations, offering a powerful avenue for studying bright X-ray binaries with NuSTAR.

The atmospheric characterisation of temperate exoplanets is becoming accessible with JWST, providing a critical connection between Solar System planets and the more commonly observed hot-Jupiters. K2-18 b, a temperate sub-Neptune orbiting an M dwarf, has emerged as a benchmark case following extensive JWST observations and ongoing debate regarding its atmospheric composition. We investigate K2-18 b using a self-consistent forward model to constrain its metallicity, composition, and thermal structure, with particular emphasis on the role of disequilibrium chemistry, photochemical hazes and clouds. For the first time in this context, we also assess the impact of photoelectrons on the atmospheric chemistry of an exoplanet. We explore a wide range of metallicities and intrinsic temperatures, evaluate haze and cloud formation, and compare the resulting transmission spectra with available JWST observations from multiple independent pipelines. We find that a high metallicity (200-400xsolar) H2-rich atmosphere consistently reproduces the observed transit spectra, independent of the data reduction pipeline used. The atmospheric composition is strongly shaped by disequilibrium chemistry, with CH4 dominating the spectrum alongside contributions from CO2 and OCS, and a potential contribution from C2H4 at mid-infrared wavelengths. Photoelectrons enhance the production of several disequilibrium species, particularly nitrogen-bearing molecules. Photochemical hazes play a key role in shaping the thermal structure, producing a temperature minimum near the 10-100 mbar range that enables efficient H2O condensation, suppressing its gaseous abundance. Under sufficiently strong haze cooling, NH4SH condensation provides a natural explanation for the apparent absence of NH3 in the observed spectra. No additional molecular species beyond those considered here are required to reproduce the observed spectra.

For frontier experiments operating at the edge of detectability, instrument design directly determines the probability of discovery. We introduce Conditional Neural Bayes Ratio Estimation (cNBRE), which extends neural Bayes ratio estimation by conditioning on design parameters, enabling a single trained network to estimate Bayes factors across a continuous design space. Applied to 21-cm radio cosmology with simulations representative of the REACH experiment, the amortised nature of cNBRE enables systematic design space exploration that would be intractable with traditional point-wise methods, while recovering established physical relationships. The analysis demonstrates a ~20 percentage point variation in detection probability with antenna orientation for a single night of observation, a design decision that would be trivial to implement if determined prior to antenna construction. This framework enables efficient, globally-informed experimental design optimisation for a wide range of scientific applications.

Lightcone observations are the natural data format of galaxy surveys, but their evolving geometry breaks the translational symmetry assumed by standard convolutional neural networks (CNNs). In particular, applying CNNs to 3D gridded lightcone data implicitly treats the line-of-sight direction as translationally invariant, despite encoding cosmic time evolution. We propose a simple alternative (CNN+2D) that divides the lightcone into redshift slices, projects each onto a HEALPix sphere, and analyzes them with a 2D CNN. Using \texttt{AbacusSummit} halo lightcone mocks ($0.3<z<0.8$, $40^\circ\times40^\circ$), we compare this approach with fully connected networks (FC) applied to different summary statistics, including spherical harmonic coefficients ($a_{\ell m}$), wavelet scattering transform (WST) coefficients, and the angular two-point correlation function (2PCF), along with standard 2PCF likelihood and Fisher forecasts. We find that multiple statistics beyond CNNs can achieve competitive performance: FC networks combined with $a_{\ell m}$ and especially WST significantly outperform 2PCF-based methods, with FC+WST yielding the best overall parameter constraints across cosmologies. Meanwhile, for a fiducial cosmology with multiple realizations, the CNN+2D approach achieves the smallest statistical uncertainties. These results demonstrate that both learned features and carefully constructed summary statistics can effectively extract cosmological information from lightcone data, providing flexible and robust analysis strategies for upcoming surveys such as DESI, Euclid, and CSST.

We present the first radiation-hydrodynamical simulations of the formation of a born-again planetary nebula (PN) triggered by a late thermal pulse (LTP). The 2D radiation-hydrodynamic simulations, performed with the {\sc pluto} code, have been consistently coupled to stellar evolution calculations using the Modules for Experiments in Stellar Astrophysics ({\sc mesa}) code. Very particularly the stellar evolution model uses (i) updated opacity tables for H-deficient, C-rich mixtures during the LTP, and (ii) a mass-loss prescription tailored for H-deficient [Wolf-Rayet]([WR])-type winds during the post-LTP phase. Our stellar model reproduces the nearly complete depletion of H expected after an LTP event, while matching the observed abundances and spectral types of iconic [WR]-type central stars of PNe. The simulations show for the first time that the H-deficient LTP ejecta forms a transient double-shell structure which, after $\sim$1000 yr, becomes fully mixed with the H-rich PN. The ejecta mass ($\sim3.4\times10^{-4}$~M$_\odot$) is too small to leave a lasting imprint on the nebular abundances, predicting H-rich PNe around [WR] central stars. The injection of LTP material into the hot bubble drives turbulence, clump formation, and enhanced mixing, providing an explanation to the larger expansion velocities and larger turbulent nebular structures of PNe with [WR] central stars compared to those with H-rich central stars. These results provide robust support for the born-again scenario as the origin of H-deficient [WR] central stars within H-rich PNe.

Rotation and magnetic fields in the cores of evolved massive stars in their final phase are thought to play an important role in the subsequent supernova explosion and the formation of a compact object, especially in hyperenergetic explosions. However, the interplay between rotation, magnetic fields, and convection up to the final collapse is a nonlinear, multidimensional effect that is difficult to capture with standard one-dimensional (1D) stellar evolution models. We quantify the magnetic angular momentum (AM) transport in the convective oxygen burning shell in a three-dimensional (3D) rotating core-collapse progenitor model. We find that the radial direction of magnetic AM transport is directly related to the Rossby number of the convective oxygen shell. We also analyze the magnetic energy, which sets the amplitude of the magnetic AM flux. The magnetic energy is determined both by rotation and the nuclear energy generation rate analogously to low-mass stars like the Sun. Based on these results, we construct a 1D model of magnetic AM transport in the convection zone for the first time in terms of properties of a given stellar evolution model. This model successfully reproduces the AM transport in the 3D model when the convective dynamo is in a quasi-steady state. Notably, our model for radial AM transport is the first to account for inward AM flux. This may result in interesting differences compared to the conventional treatment of magnetic AM transport in stellar evolution models, which assume AM is transported outward by a purely diffusive process.

Cosmological models exhibiting crossing of the phantom divide improve the fit to current data, suggesting late-time dark energy (DE) dynamics at $\sim3\sigma$ CL. However, they favor low values of $H_0$, in tension with SH0ES. This may point to the presence of new physics prior to the decoupling era. In this work, we reconstruct the background DE functions using the Weighted Function Regression (WFR) method, introducing three main improvements compared to our previous JCAP 12 (2025) 049. First, we adopt the Frequentist-Bayesian approach for the weights. Second, we combine CMB and BAO with the DES-Dovekie SNIa sample and compare our findings with those derived from Pantheon+, still assuming standard recombination. Third, we study in a model-independent manner the viability of early-time ``solutions'' to the Hubble tension and how they affect the evidence for dynamical DE at late times, under the influence of the SH0ES and the more conservative CCHP calibration of the cosmic ladders, separately. We find that, if the physics prior to decoupling is unmodified, the probability of phantom crossing is $\sim 96.7\text{--}98.5\%$, with $\Lambda$CDM excluded at $\sim 2.5\sigma$ and $\sim 3\sigma$ CL. New physics before recombination can alleviate the Hubble tension, but requires extremely large values of the reduced matter density parameter when the SH0ES calibration is employed, in strong tension with those inferred from full CMB analyses. This raises serious concerns about the actual viability of these models to explain the SH0ES measurement. We find that phantom crossing, while not excluded, is no longer required, with only a very mild preference for quintessence. Nevertheless, given the aforesaid tension in $\omega_m$, it would be rash to draw firm conclusions about how the dynamical DE signal is affected in these scenarios. [abridged]

We investigate static, spherically symmetric halo configurations within Unified Dark Matter (UDM) scalar-field models, developing a systematic mapping between standard cold dark matter (CDM) density profiles and their UDM counterparts. Exploiting the equivalence-class structure of UDM models, we show that, in principle, different Lagrangian realisations can share the same weak-field rotation curve while exhibiting distinct field properties. We reconstruct the effective energy density, radial and tangential pressures from a phenomenological circular velocity profile, ensuring the absence of ghosts and instabilities and the preservation of the Null Energy Condition (NEC). Applying our procedure to several commonly used CDM halo profiles -- including Persic, Salucci \& Stel, NFW, and Burkert models -- we demonstrate that their phenomenological success can be retained within a relativistic UDM framework, reproducing the observed flatness of rotation curves without introducing separate dark matter and dark energy components.

We report the discovery of 31 new open clusters (OCs) identified in \textit{Gaia}~DR3 data through a systematic search over 220 adjacent $1^\circ\times1^\circ$ fields towards the Galactic anticentre, in the direction of the Perseus arm gap. Eight of them display low-density structures, possibly indicating open cluster remnants properties. The objects were identified and characterized through a combined analysis of photometric, kinematic, and spatial distributions, a methodology successfully applied in our previous works. Their structural properties, mean proper motions, ages, distances and reddening were derived and their centres cross-matched with the available catalogues. The clusters are low-concentrated systems and are mostly located within $3<d<5$ kpc, exhibiting reddening up to $E(B-V)\approx1.5$, and ages from $\sim$20 Myr to 1 Gyr. The new OCs represent a significant increase in the anticentre cluster census: $31\%$ for $3<d<4$ kpc and $12\%$ for $d>4$ kpc. They do not belong to the Perseus arm, but may be associated with the Outer Norma arm. The Gulf of Camelopardalis region appears as an interruption in the Perseus arm, possibly reflecting low star-formation activity, dust obscuration, or that the Milky Way is a flocculent, rather than a grand-design spiral galaxy.

Asteroseismology has become a powerful tool in stellar astrophysics, offering unprecedented insights into the internal structures and dynamics of stars. It enables precise characterization of stellar interiors across a wide range of stellar masses and of evolutionary phases, from the main sequence to the white dwarf phase. At the same time, the number of detected close stellar and planetary companions throughout the entire stellar evolutionary phases has increased significantly, prompting key questions about the interplay between stellar evolution and binarity. We investigate the competition between gravity waves (IGW) excited by internal convection and those excited by tides in stars that host a companion. By modelling the energy and angular momentum luminosities transported by IGWs stochastically excited by convection and by tides, we seek to quantify their relative contributions and identify the key parameters that govern their efficiency. We compute the energy and angular momentum luminosities transported by both types of waves for a range of stellar masses and evolutionary stages, with a particular focus on understanding how the presence of a companion influences the angular momentum transport of the radiative layers of the host star. The competition between the two excitation mechanisms is sensitive to the mass and orbital properties of the companion, as well as the internal structure of the host star. We find that for a Jupiter-mass companion, the stochastic excitation dominates over tidal excitation during all evolutionary phases. Only for close-in stellar companions around late-type stars does the tidal excitation become more efficient. The presence of a companion is unlikely to significantly alter the internal angular momentum transport in the radiative layers of the host star, simplifying the modelling of IGW-driven angular momentum transport in stars that host a companion.

We analyze a sample of low surface brightness dwarf galaxies (mu_e,g > 24.2 mag arcsec^-2), detected using interpretable machine learning tools from the DES survey. We use the Tanoglidis et al. (2021) sample, identified with machine learning, supplemented by Thuruthipilly et al. (2024). We focus on the Fornax-Eridanus Supercluster, where our group determined its 3D filamentary spine using massive galaxies. We study the effect of the large-scale environment on dwarfs in the Fornax-Eridanus Complex. To do this, we compare the properties of dwarfs in clusters, groups, and the field, and examine how these properties vary with distance to the spine of the Fornax Wall. We check if dwarfs trace the Fornax Wall spine, defined by massive galaxies. We identify Fornax Wall members from our photometric dwarf catalog, dividing them into i) within one virial radius of a galaxy group or cluster and ii) outside this radius (field galaxies). We assume dwarfs near the Fornax Wall are at the same distance as the massive galaxies. We then study their distribution within the complex. We probe the morphology-density relation and examine galaxy properties versus distance from the Fornax Wall spine. Red dwarfs are mostly in or near groups close to the Fornax Wall, dominating the population, while blue dwarfs dominate the field. Larger-sized red dwarfs tend to reside in group environments, with significantly larger effective radii than those in the field. Red dwarfs are more concentrated towards the Fornax Wall than blue dwarfs. This suggests that the group environment plays a significant role in the evolution of dwarfs. Mass density distribution in field and group/cluster is similar, indicating the group/cluster population could be an aged version of the field. The group/cluster objects with excess sizes must have been made through interactions in the groups/clusters.

We investigate the evolution of dark matter halo shapes in cosmological N-body simulations both in scale free Einstein-De Sitter (EdS) and $\Lambda$CDM cosmologies. We compute the axis ratios ($q=b/a,s=c/a$) of well resolved central halos using the shape tensor. These halos are identified using two different halo finding algorithms, SUBFIND and ROCKSTAR. We find that at fixed mass, halos become more spherical with decreasing redshift. The distribution $P(q,s)$ along with their median values ($q$ and $s$) shows self-similar behaviour as a function of mass scaled by the non-linear mass, $(M/M_{nl})$ across power-law spectral indices for scale free EdS models. However the median $q$ and $s$ show a tighter self-similar evolution as a function of peak height $\nu=\delta_c/\sigma(M,z)$. We find that the median $q(\nu)$ and $s(\nu)$ are consistent with an evolution along a universal curve described by $y=\alpha-\delta \tanh \left[ \omega \left(\log_{10}(\nu) - \mu\right)\right]$ across the spectral indices ranging from $n=-1.0$ to $n=-2.2$. Our results hold for both SUBFIND and ROCKSTAR, although there are some differences between them. The universality of the evolution of median $q(\nu)$ and $s(\nu)$ also holds for the $\Lambda$CDM runs, although with a different behaviour at small $\nu$ compared to the scale free models. The width of the distributions of $P(q)$ and $P(s)$ in both, scale-free and $\Lambda$CDM, classes of simulations can be reduced further by classifying halos as oblate, triaxial and prolate, each of which also follows a universal behaviour. Although oblate halos are relatively rare at all redshifts, their fraction increases over time at the expense of the other two populations.

When blazars flare, their optical position moves. We show this by combining Gaia DR3 proper motions with epoch photometry for blazars with strong optical jet emission. In 60 of 74 sources with significant proper motion, rising flux drives the centroid upstream while fading flux drives it downstream - a near-universal pattern captured by a simple two-component model of constant extended emission and a flaring region. Using this connection, we geometrically localize the optical flares to within <1 mas of the VLBI position - a few parsecs at typical blazar distances - placing them in the innermost jet or accretion disk. This purely geometric method requires no multi-wavelength correlations or model-dependent assumptions, and provides an independent spatial anchor for localizing higher-energy flares. Per-epoch astrometry from Gaia DR4 is set to tighten our constraints even further.

The quadratic theory of gravity is the unique renormalizable theory of quantum gravity in 4 dimensions, as proved by K. S. Stelle in 1977. Over the decades, the theory has been understood to contain a massive tensor ghost, and several attempts have been made to evade its adverse effects by proposing new quantization prescriptions and interpretations. In this paper, we show that the additional spin--2 of quadratic gravity can be turned into a healthy inverted harmonic oscillator (IHO)-like instability, which can be quantized consistently with direct-sum quantum field theory (DQFT), which incorporates geometric superselection sectors. Such modes possess a well-defined quantum description yet do not admit a particle interpretation and are not part of the asymptotic spectrum, being characterized by hyperbolic evolution and spacelike momentum support. We argue that, as a consequence, the extra spin--2 degree of freedom remains off-shell and effectively decoupled from ordinary matter fields, avoiding unitarity violations in observable processes. We argue that this IHO instability is a prevalent feature of fundamental physics, whether it concerns quantum fields on curved spacetimes or the Higgs $\mathbb{Z}_2$ symmetry breaking in the Standard Model of particle physics. Thus, we demonstrate that our new understanding of quadratic gravity offers a fundamental pathway to quantum gravity and a safe beginning for the Universe. Furthermore, we derive key observational predictions of this construction in the view of primordial gravitational waves with new bounds on the tensor-to-scalar ratio and the parity asymmetric features on the large angular scales.

We derive expressions for the first and second derivatives of the quintessence potential $V(\phi)$, in terms of $\lambda = -V^{\prime}/V$ and $\Gamma = (V^{\prime \prime}/V)/(V^\prime/V)^2$, as functions of the quintessence density fraction $\Omega_\phi$ and the cosmographic parameters $q$, $j$, and $s$. Our mapping is not explicitly a function of the equation of state parameter $w$. We use these results, along with recent observational data, to derive expansions of $V(\phi)$ about the present-day value of the scalar field, $\phi_0$.

In 2025, the French Society for Astronomy \& Astrophysics (SF2A), gave the environmental transition commission the opportunity to share their considerations during a plenary session at the annual SF2A conference. This year, the presentation focused on some of the main results obtained from the survey entitled 'French astronomy and astrophysics research activities in the face of the environmental crisis, from 2019 to 2024'. The survey was initiated in 2019 by the group 'Environnement-Transition' (coordinated by P. Martin) at IRAP, whose results were presented during the SF2A annual conference 2019 in Nice. The survey was updated in 2024 by the CNRS INSU-AA prospective working group 'Climate and ecological challenge' (coordinated by S. Bontemp). The SF2A environmental transition commission took on the survey to the French institutes, sorted the answers and extracted the preliminary results. The full results will be published at the end of 2025 in the final CNRS INSU-AA 2024 prospective document. This publication presents a selection of pieces from the full survey, along with a few of the main discussions it triggers.

We present results from directed searches for continuous gravitational waves from a sample of 15 nearby supernova remnants, likely hosting young neutron star candidates, using data from the first eight months of the fourth observing run (O4) of the LIGO-Virgo-KAGRA Collaboration. The analysis employs five pipelines: four semi-coherent methods -- the Band-Sampled-Data directed pipeline, Weave and two Viterbi pipelines (single- and dual-harmonic) -- and PyStoch, a cross-correlation-based pipeline. These searches cover wide frequency bands and do not assume prior knowledge of the targets' ephemerides. No evidence of a signal is found from any of the 15 sources. We set 95\% confidence-level upper limits on the intrinsic strain amplitude, with the most stringent constraints reaching $\sim 4 \times 10^{-26}$ near 300 Hz for the nearby source G266.2$-$1.2 (Vela Jr.). We also derive limits on neutron star ellipticity and $r$-mode amplitudes for the same source, with the best constraints reaching $\lesssim 10^{-7}$ and $\lesssim 10^{-5}$, respectively, at frequencies above 400 Hz. These results represent the most sensitive wide-band directed searches for continuous gravitational waves from supernova remnants to date.

Radiation reaction cooling plays an important role in describing the extreme plasma conditions found in the magnetospheres of astrophysical compact objects. Strong electromagnetic fields, characteristic of these environments, can trigger the development of anisotropic ring-shaped plasma distributions with inverted Landau populations in momentum space. In this work, we present the first systematic investigation of this mechanism in realistic astrophysical configurations, by accounting for how non-uniform electromagnetic field geometries and general-relativistic effects modify the phase-space dynamics of radiatively cooled plasmas. We demonstrate analytically that drift velocities favour the formation of spiral-shaped momentum distributions that still display inverted Landau populations, and estimate the minimum and maximum plasma injection distances required for inverted momentum distributions to be able to power the emission of coherent radiation through kinetic instabilities. From numerical simulations, we conclude that curved spacetime increases the gradient of the distribution function responsible for the development of kinetic instabilities, and prolongs the persistence of the inverted momentum structure relative to flat spacetime, confirming that realistic astrophysical conditions preserve and enhance the conditions necessary for synchrotron-powered emission of coherent radiation to occur.

We consider circular particle motion under the action of an unspecified force in a static spherically symmetric spacetime. We derive the machinery that allows one to find the force acting on a circular particle and deduce whether its position is stable or not. This also allows one to extend the definition of ISCO to the case of a non-zero external force. By conducting the near-horizon expansion, we obtain that for any non-extremal black holes, the acceleration for extremal ones is finite, and for ultraextremal (multiple) horizons it tends to zero. Applying the derived machinery to the case of the Schwarzschild metric assuming that a force is constant, we scrutiny how the number of orbits for a given force depends on its value. In particular, if a force is big enough, an additional branch of solutions appears that was absent in the case of geodesic motion. Then, for various circular orbits, we numerically investigate their stability. A similar problem is solved for the Reissner-Nordstrom (RN) metric and uncharged particles. It appears that for the near-extremal and extremal RN black holes, there exist near-horizon circle trajectories (in contrast to the nonextremal case). For the ISCO, the dependence of the orbit radius on $\kappa$ (the surface gravity) is similar to that in the case of neutral particles moving in the background of rotating black holes. In addition, two scenarios of high-energy particle collisions near such orbits are considered, and it is found that dependence on $\kappa$ is also similar to that for rotating black holes.

Quantum electrodynamics (QED) becomes nonlinear when the magnetic field strength surpasses the critical Schwinger limit $B_Q \approx 4.41\cdot 10^{13}$ G. This limit is surpassed, for example, in the magnetospheres of a specific class of neutron stars known as magnetars, which has important consequences for magnetospheric plasma dynamics due to modifications in scattering cross sections. Using a formalism previously applied to the study of magnetic catalysis, I calculate the cross sections of all tree-level 1-to-2, 2-to-1, and 2-to-2 particle QED scattering processes that do not include a photon propagator. The calculations are done in a strong background magnetic field and the results are implemented into an open-source Python package. This article focuses on presenting the formalism and computational techniques required for the calculations, while the impact of the results on, e.g., magnetospheric plasma dynamics is discussed in a companion letter (Kiuru et al. 2026).

Large-eddy simulations (LES) with an appropriate subgrid-scale (SGS) model provide a powerful tool for investigating real-world turbulence. The Smagorinsky model, one of the simplest and most used SGS models, often shows an over-dissipative behavior even when using dynamic procedures to adjust the model coefficient. By incorporating the structural or geometrical information of turbulence provided by helicity (velocity-vorticity correlations), the helicity SGS model is expected to alleviate these issues in the standard Smagorinsky framework, in which only information of turbulence intensity is considered through the turbulent energy. The validity of helicity SGS models is investigated here with the aid of direct numerical simulations (DNSs). Using configurations with and without net rotation, and with large-scale helicity gradients sustained by a mechanical forcing, we show that to better model SGS turbulence, SGS helicity effects should be incorporated into the model together with the Smagorinsky-like eddy viscosity.

The optimization of large experiments in fundamental science, such as detectors for subnuclear physics at particle colliders, shares with the optimization of complex systems for industrial or societal applications the common issue of addressing the inter-relation between parameters describing the hardware used in data production and parameters used to analyse those data. While in many cases this coupling can be ignored -- when the problem can be successfully factored into simpler sub-tasks and the latter addressed serially -- there are situations in which that approach fails to converge to the absolute maximum of expected performance, as it results in a mis-alignment of the optimized hardware and software solutions. In this work we consider a few use cases of interest in fundamental science collected primarily from particle physics and related areas, and a pot-pourri of industrial and societal applications where the matter is similarly of relevance. We discuss the emergence of strong hardware-software coupling in some of those systems, as well as co-design procedures that may be deployed to identify the global maximum of their relevant utility functions. We observe how numerous opportunities exist to advance methods and tools for hardware-software co-design optimization, bridging fundamental science and industry through application- and challenge-driven projects, and shaping the future of scientific experiments and industrial systems.

We present Automatic Laplace Collapsed Sampling (ALCS), a general framework for marginalising latent parameters in Bayesian models using automatic differentiation, which we combine with nested sampling to explore the hyperparameter space in a robust and efficient manner. At each nested sampling likelihood evaluation, ALCS collapses the high-dimensional latent variables $z$ to a scalar contribution via maximum a posteriori (MAP) optimisation and a Laplace approximation, both computed using autodiff. This reduces the effective dimension from $d_\theta + d_z$ to just $d_\theta$, making Bayesian evidence computation tractable for high-dimensional settings without hand-derived gradients or Hessians, and with minimal model-specific engineering. The MAP optimisation and Hessian evaluation are parallelised across live points on GPU-hardware, making the method practical at scale. We also show that automatic differentiation enables local approximations beyond Laplace to parametric families such as the Student-$t$, which improves evidence estimates for heavy-tailed latents. We validate ALCS on a suite of benchmarks spanning hierarchical, time-series, and discrete-likelihood models and establish where the Gaussian approximation holds. This enables a post-hoc ESS diagnostic that localises failures across hyperparameter space without expensive joint sampling.