Astrophysics
Cosmology, galaxies, stellar physics, and astronomical instrumentation
6,198 papers
Unveiling the 3D structure of the central molecular zone from stellar kinematics and photometry: The 50 and 20 km/s clouds
The central molecular zone (CMZ), surrounding the Galactic centre, is the largest reservoir of dense molecular gas in the Galaxy. Despite its relative proximity, the 3D structure of the CMZ remains poorly constrained, primarily due to projection effects. We aim to constrain the line-of-sight location of two molecular clouds in the CMZ -- the 50 and 20 km/s clouds -- and to investigate their possible physical connection using stellar kinematics and photometry. This study serves as a pilot for future applications across the full CMZ. We estimated the line-of-sight position of the clouds by analysing stellar kinematics, stellar densities, and stellar populations towards the cloud regions and a control field. We find an absence of westward moving stars in the cloud regions, which indicates that they lie on the near side of the CMZ. This interpretation is supported by the stellar density distributions. The similar behaviour observed in the two clouds, as well as in the region between them (the ridge), suggests that they are located at comparable distances and are physically linked. We also identified an intermediate-age stellar population (2-7 Gyr) in both regions, consistent with that observed on the near side of the CMZ. We estimated the line-of-sight distances at which the clouds and the ridge become kinematically detectable (i.e. where the proper motion component parallel to the Galactic plane differs from that of the control field at the 3 sigma level) by converting their measured proper motions parallel to the Galactic plane using a theoretical model of the stellar distribution. We find that the 50 and 20 km/s clouds are located at $43\pm8$ pc and $56\pm11$ pc from Sgr A*, respectively, and that the ridge lies at $56\pm11$ pc; this supports the idea that the clouds are physically connected through the ridge.
2601.05252Jan 2026
View
Mimicking Phantom Dark Energy with Evolving Dark Matter Mass
We present a general method to reproduce a given cosmological background through energy exchange between dark energy (DE) and dark matter (DM). This can be simply realized with a standard quintessence scalar field that controls the DM mass. In particular a background with phantom crossing can be effectively realized without introducing ghosts or other pathologies. For example one can reproduce exactly the background that gives the best fit to the recent DESI+CMB+DESY5 data, within the Chevallier-Polarski-Linder (CPL) parametrization of DE. Although the background evolution is identical, the perturbations differ, leading to modified growth of structures. If the DM mass varies at late times, early-time observables are not modified and can reproduce the main predictions of the target model, but late-time observables are affected. We discuss in particular the effects on the matter power spectrum, CMB lensing and ISW effect. When reproducing the best fit CPL background model, this scenario generically predicts $\mathcal{O}(10\%)$ deviations in such observables. However, for suitable choices of parameters, effects on the matter power spectrum can be smaller, motivating a detailed study. In general, energy exchange between DE and DM generates a mismatch between the matter power spectrum and the gravitational potential amplitudes compared to the decoupled case, that can lead to deviations observable in future experiments.
2601.05235Jan 2026
View
Three-dimensional scene reconstruction using Roman slitless spectra
The Nancy Grace Roman Space Telescope will carry out a wide-field imaging and slitless spectroscopic survey of Type Ia Supernovae to improve our understanding of dark energy. Crucial to this endeavor is obtaining supernova spectra uncontaminated by light from their host galaxies. However, obtaining such spectra is made more difficult by the inherent problem in wide-field slitless spectroscopic surveys: the blending of spectra of close objects. The spectrum of a supernova will blend with the host galaxy, even from regions distant from the supernova on the sky. If not properly removed, this contamination will introduce systematic bias when the supernova spectra are later used to determine intrinsic supernova parameters and to infer the parameters of dark energy. To address this problem we developed an algorithm that makes use of the spectroscopic observations of the host galaxy at all available observatory roll angles to reconstruct a three-dimensional (3d; 2d spatial, 1d spectral) representation of the underlying host galaxy that accurately matches the 2d slitless spectrum of the host galaxy when projected to an arbitrary rotation angle. We call this ``scene reconstruction''. The projection of the reconstructed scene can be subtracted from an observation of a supernova to remove the contamination from the underlying host. Using simulated Roman data, we show that our method has extremely small systematic errors and significantly less random noise than if we subtracted a single perfectly aligned spectrum of the host obtained before or after the supernova was visible.
2601.05233Jan 2026
View
Mitigating Simulator Dependence in AI Parameter Inference for the Epoch of Reionization: The Importance of Simulation Diversity
The 21cm signal of neutral hydrogen contains a wealth of information about the poorly constrained era of cosmological history, the Epoch of Reionization (EoR). Recently, AI models trained on EoR simulations have gained significant attention as a powerful and flexible option for inferring parameters from 21cm observations. However, previous works show that AI models trained on data from one simulator fail to generalize to data from another, raising doubts about AI models' ability to accurately infer parameters from observation. We develop a new strategy for training AI models on cosmological simulations based on the principle that increasing the diversity of the training dataset improves model robustness by averaging out spurious and contradictory information. We train AI models on data from different combinations of four simulators, then compare the models' performance when predicting on data from held-out simulators acting as proxies for the real universe. We find that models trained on data from multiple simulators perform better on data from a held-out simulator than models trained on data from a single simulator, indicating that increasing the diversity of the training dataset improves a model's ability to generalize. This result suggests that future EoR parameter inference methods can mitigate simulator-specific bias by incorporating multiple simulation approaches into their analyses.
2601.05229Jan 2026
View
Addressing Known Challenges in Solar Flare Forecasting I: Limb-Flare Prediction with a 4-pi Full-Heliosphere Framework
A demonstrated failure mode for operational solar flare forecasting is the inability to forecast flares that occur near, or just beyond, the solar limb. To address this shortcoming, we develop a "4pi" full-heliosphere event forecasting framework and evaluate its statistical classification ability against this specific challenge. A magnetic surface flux transport model is used to generate full-sun maps of the photospheric radial magnetic field from which active regions (ARs) are identified and tracked using a new labeling scheme that is observer-location agnostic and allows for post-facto modifications. Flare-relevant magnetic parameters couple to a "visibility" index that specifies AR location relative to the visible solar limb and expected flare detection. Flare labels are assigned according to peak Soft X-ray flux, and a statistical classification is performed using nonparametric discriminant analysis. A version where new or emerging ARs on the far ("invisible" side of the Sun are incorporated into the model by way of far-side helioseismology, is also tested. We evaluate the new framework by its performance specifically including the limb areas using Brier Skill Score and ROC Skill Score, finding improvement at the 2-sigma level or less. However, we do find that the number of False Negatives, or "missed" forecasts decreases, and find strong evidence that the additional information provided by the far-side helioseismology can help predict near- and just-beyond-limb flares, particularly for East-limb events. While individual components of this framework could be improved, we demonstrate that a known failure mode for solar flare forecasting can be mitigated with available resources.
2601.05209Jan 2026
View
Symbolically regressing dark matter halo profiles using weak lensing
The structure of dark matter haloes is often described by radial density profiles motivated by cosmological simulations. These are typically assumed to have a fixed functional form (e.g. NFW), with some free parameters that can be constrained with observations. However, relying on simulations has the disadvantage that the resulting profiles depend on the dark matter model and the baryonic physics implementation, which are highly uncertain. Instead, we present a method to constrain halo density profiles directly from observations. This is done using a symbolic regression algorithm called Exhaustive Symbolic Regression (ESR). ESR searches for the optimal analytic expression to fit data, combining both accuracy and simplicity. We apply ESR to a sample of 149 galaxy clusters from the HSC-XXL survey to identify which functional forms perform best across the entire sample of clusters. We identify density profiles that statistically outperform NFW under a minimum-description-length criterion. Within the radial range probed by the weak-lensing data ($R \sim 0.3 - 3$ h$^{-1}$ Mpc), the highest-ranked ESR profiles exhibit shallow inner behaviour and a maximum in the density profile. As a practical application, we show how the best-fitting ESR models can be used to obtain enclosed mass estimates. We find masses that are, on average, higher than those derived using NFW, highlighting a source of potential bias when assuming the wrong density profile. These results have important knock-on effects for analyses that utilise clusters, for example cosmological constraints on $σ_8$ and $Ω_m$ from cluster abundance and clustering. Beyond the HSC dataset, the method is readily applicable to any data constraining the dark matter distribution in galaxies and galaxy clusters, such as other weak lensing surveys, galactic rotation curves, or complementary probes.
2601.05203Jan 2026
View
How deep can a cosmic void be? Voids-informed theoretical bounds in Galileon gravity
We establish a general, void-based consistency test for Galileon scalar-tensor theories. We show that the previously reported unphysical breakdown of the predicted Newtonian force in certain Galileon models is controlled by a single condition linking non-linear void dynamics to the cosmic expansion history. This connection yields a redshift-dependent upper bound on the allowed depth of voids and promotes this requirement to a new viability condition, complementary to standard stability criteria. As an example, we apply this void-based criterion to a linear parameterization in the scale factor constrained by theoretical and observational bounds; we find that $\sim 60\%$ of the parameter space is excluded, with most problematic models failing by $z\lesssim 10$. These results position cosmic voids as sharp, broadly applicable, theory-informed filters for viable modified gravity, enabling more informed priors and parameter-space choices in future cosmological inference.
2601.05145Jan 2026
View
Detection of time delay between UV and X-ray variability in Mrk 1044 using AstroSat observations
Active galactic nuclei are known to exhibit flux variations across the entire electromagnetic spectrum. Among these, correlations between UV/optical and X-ray flux variations serve as a key diagnostics for understanding the physical connection between the accretion disk and the corona. In this work, we present the results of analysis of ultraviolet (UV) and X-ray flux variations in the narrow line Seyfert 1 galaxy Mrk 1044. Simultaneous observations in the far-UV band (FUV: 1300$-$1800 Å) and the X-ray band (0.5$-$7 keV) obtained during 31 August $-$ 8 September 2018 with the Ultraviolet Imaging Telescope and the Soft X-ray Telescope onboard \textit{AstroSat} were used for this study. Significant flux variability was detected in both FUV and X-ray bands. The fractional root mean square variability amplitude ($F_{\rm var}$) was found to be 0.036 $\pm$ 0.001 in the FUV band and 0.384 $\pm$ 0.004 in the X-ray band. To explore potential time lag between the two bands, cross-correlation analysis was performed using both the interpolated cross-correlation function (ICCF) and just another vehicle for estimating lags in nuclei (JAVELIN) methods. Results from both approaches are consistent within 2$σ$ uncertainty, indicating that X-ray variations lead the FUV variations, with measured lags of 2.25$\pm$0.05 days (ICCF) and $2.35_{-0.01}^{+0.02}$ days (JAVELIN). This is the first detection of a time delay between UV and X-ray variations in Mrk 1044. The observed UV lag supports the disk reprocessing scenario, wherein X-ray emission from the corona irradiates the accretion disk, driving the observed UV variability.
2601.05135Jan 2026
View
Multigroup Radiation Diffusion on a Moving Mesh: Implementation in RICH and Application to Tidal Disruption Events
Radiation-hydrodynamics (RHD) determines the bulk evolution and observable emission in a wide variety of high-energy astrophysical phenomena. Due to their complexity, RHD problems must usually be studied through numerical simulation. We have extended the publicly available RICH code, which previously solved the equations of RHD in the limit of grey flux-limited diffusion (FLD), to operate with a multigroup FLD solver. RICH is a semi-Lagrangian code that solves the equations of RHD on an unstructured moving mesh, and is the first multigroup RHD moving mesh code, making it uniquely applicable to problems with extreme dynamic range and dynamically important radiation forces. We validate our multigroup module against multiple analytic benchmarks, including a novel test of the RHD Doppler term. The computational efficiency of the code is aided by a novel scheme to accelerate convergence in optically thick cells. Finally, we apply multigroup RICH in a pilot study of a stellar tidal disruption event (TDE), using a $10^4 M_\odot$ intermediate-mass black hole. Our simulations self-consistently produce a bright early-time X-ray flash prior to peak optical/UV light, in qualitative agreement with post-processing of (grey) RICH simulations of supermassive black hole TDEs, as well as X-ray observations of the TDE AT 2022dsb.
2601.05120Jan 2026
View
Surveying exogenous species in Saturn with ALMA I. Detecting and Mapping CO
The origin of carbon monoxide (CO) in Saturn's stratosphere remains uncertain, with proposed sources including internal thermochemical production, cometary impacts, and exogenic material from the rings and icy moons (i.e. Enceladus). We aim to constrain the vertical and meridional distribution of stratospheric CO and assess the relative contributions of these potential sources. Here, we analysed high-spectral-resolution ALMA observations of the CO (J=3-2) line obtained on 25 May 2018, sampling Saturn's limb from 20°S to 69°N. CO vertical profiles were retrieved using a line-by-line radiative transfer model combined with spectral inversion techniques, testing multiple prior scenarios representative of different source hypotheses. CO is confined to a narrow layer between 0.1 and 1 mbar, with a robust negative vertical gradient and mean abundances of (3.7+/- 0.8) x 10$^{-8}$ at 0.1 mbar and (7.2 +/- 0.9) x 10$^{-8}$ at 1 mbar. The meridional distribution is statistically homogeneous, with a marginal enhancement near 60° N plausibly related to Enceladus. No significant equatorial enhancement is detected. The absence of a strong equatorial enhancement rules out a long-lived steady source associated with ring infall. The observations are most consistent with a relatively recent ($\approx$200-year-old or younger) cometary impact whose material has since been horizontally mixed, while any Cassini Grand Finale ring influx was either too recent or inefficient to affect CO abundances at the probed pressure levels.
2601.05090Jan 2026
View
Constraining the Primoridal Black Hole Abundance with Space-Based Detectors
Overdense regions can collapse into primordial black holes (PBHs) in the early universe, which are a compelling candidate for dark matter. Current constraints leave the asteroid-mass window the only possible one for PBH to account for all the dark matter, which can only be probed indirectly by the scalar-induced gravitational waves (GWs) sourced by the curvature perturbation which forms PBH. In this work, we explore the capabilities of future space-based gravitational wave detectors, including LISA, Taiji, and TianQin, to constrain such induced GWs as well as the PBH abundance. We systematically account for the width of the primordial curvature power spectrum, and find that the asteroid-mass window can be fully probed by all three space-based interferometers. If PBHs constitute the majority of dark matter, the induced GW leaves a strong signal in the mHz band with a signal-to-noise ratio of $10^3\sim10^4$.
2601.05069Jan 2026
View
Energetic particles accelerated via turbulent magnetic reconnection in protoplanetary discs -- I: ionisation rates
Ionisation controls the chemistry, thermal balance, and magnetic coupling in protoplanetary discs. However, standard ionisation vectors such as stellar UV, X-rays, Galactic Cosmic Rays (GCRs) might not be efficient enough, as UV/X-rays are attenuated rapidly with depth, while GCRs are modulated. Turbulence-induced magnetic reconnection in disc atmospheric layers offers a physically motivated, in-situ source of energetic particles (EPs) that has never been considered. We quantify the ionisation and heating produced by EPs accelerated by turbulent reconnection, identify where they dominate over X-rays and GCRs, and determine energetic thresholds for their relevance. We provide scalable diagnostics tied to the local energy budget. We adopt a Fermi-like acceleration model with parameters linked to a turbulent reconnection geometry trigger by the magneto-rotational instability, yielding a steady-state energy distribution of the EP forming a power-law of index $p=2.5$. We propagate electrons and protons through the disc and compute primary and secondary ionisation and associated heating on a fiducial T Tauri disc model background. The non-thermal normalisation is set by the fraction of local viscous accretion energy dissipation channelled to EPs, parametrised by $κ$. For $κ\gtrsim 0.4\%$, EPs ionisation overpass standard sources like X-rays and GCRs in the disc atmosphere and intermediate/deep layers out to radii of a few tens of AU. Even at $κ\sim 0.025\%$, EPs contribute at the few-percent level, thus are chemically and dynamically relevant. EP-induced heating complements UV/X-ray heating in the atmosphere and persists deeper. These results identify EPs accelerated by turbulence-induced magnetic reconnection as a rather robust, disc-internal ionisation channel that should be included in thermo-chemical and dynamical models of protoplanetary discs.
2601.05060Jan 2026
View
Jet launching from the Kerr black hole magnetosphere: An electrogeodesic approach
The launch of relativistic jets of plasma on astrophysical to cosmological scales are observed in a variety of astrophysical sources, from active galactic nuclei to X-rays binaries. While these jets can be reproduced by the general relativistic magneto-hydrodynamic (GRMHD) and particle-in-cells (GRPIC) simulations of the dynamical Kerr magnetosphere, the development of analytic models to describe the physics of the jets has remained limited. A key challenge is to analytically describe the individual trajectories of accelerated charged particles which ultimately build up the jet and emit radiation. In this work, we provide a first simple but fully analytical model of jet launching from the Kerr magnetosphere based on the motion of charged particles. To that end, we use the integrability of electrogeodesic motion in the Kerr monopole magnetosphere to study the ejection of charged particles near the poles. This enables us to derive (i) a criterion for the rotation axis to constitute a stable latitunal equilibrium position, thereby representing an idealized jet, (ii) the expression for the magnetic frame-dragging effect, and (iii) the condition for an asymptotic observer to measure blueshifted particles emanating from the black hole surroundings. Our study reveals that particles can be accelerated only in a specific region whose maximal radius depends on the spin and magnetization of the black hole. Alongside these results, we provide a detailed review of the construction of test magnetospheres from (explicit and hidden) symmetries of the Kerr geometry and the condition for the separability of the electrogeodesic motion in a test magnetosphere, which serves as a basis for the model we study in this work.
2601.05048Jan 2026
View
2601.05006Mutual Orbit Alignment in Resolved Triple Systems
A sample of 278 triple systems with outer separations under 300 au and resolved inner pairs is studied, focusing on the mutual alignment between inner and outer orbits. The degree of alignment increases with (i) decreasing outer separation, (ii) decreasing ratio of outer and inner separations, (iii) decreasing mass of the inner primary component, and (iv) increasing inner mass ratio. There is no dependence on the outer mass ratio. The average mutual inclination is ~40deg for the full sample and ~10deg for 38 triples with primary components less massive than 1 solar and outer separations below 50 au. Inner eccentricities in aligned triples are smaller compared to misaligned ones. In another sample of 371 hierarchies with known outer orbits and inner eclipsing subsystems, only 22% show mutual alignment within 20deg, while the rest are aligned randomly. These findings match qualitatively current understanding of the formation of hierarchical systems, where the N-body dynamics dominates at large scales, while the accretion and migration shape systems closer than $\sim$100 au. Fragmentation of isolated cores apparently produces approximately aligned low-mass hierarchies.
2601.05006Jan 2026
View
Probing millisecond magnetar formation in binary neutron star mergers through X-ray follow-up of gravitational wave alerts
The nature of the remnant of a binary neutron star (BNS) merger is uncertain. Though certainly a black hole (BH) in the cases of the most massive BNSs, X-ray lightcurves from gamma-ray burst (GRB) afterglows suggest a neutron star (NS) as a viable candidate for both the merger remnant as well as the central engine of these transients. When jointly observed with gravitational waves (GWs), X-ray lightcurves from BNS merger events could provide critical constraints on the remnant's nature. We aim to assess the current and future capabilities to detect a NS remnant through X-ray observations following GW detections. To this end, we simulate GW signals from BNS mergers and the subsequent X-ray emission from newborn millisecond magnetars. The GW detectability is modeled for both current and next-generation interferometers, while the X-ray emission is reproduced using a dedicated numerical code that models magnetar spin-down and ejecta dynamics informed by numerical-relativity simulations. In our simulations, 2% - 16% of BNS mergers form millisecond magnetars. Among these, up to 70% could be detectable, amounting to up to 1 millisecond magnetar detection per year with SVOM/MXT-like instruments during the LIGO Virgo KAGRA LIGO India (LVKI) O5 run, with optimal detectability occurring about 2 hours post-merger. For next-generation GW interferometers, this rate could increase by up to three orders of magnitude, with peak detectability 3 to 4 hours post-merger. We also explore how the magnetar's magnetic field strength and observer viewing angle affect detectability and discuss optimized observational strategies. Although more likely with upcoming GW interferometers, detecting the spin-down emission of a millisecond magnetar may already be within reach, warranting sustained theoretical and observational efforts given the profound implications for mergers, GRBs, and NS physics of a single detection.
2601.04990Jan 2026
View
Through Thick and Thin: The Cosmic Evolution of Disk Scale Height
To investigate the formation and evolution of vertical structures in disk galaxies, we measure global $\operatorname{sech}^2$ scale heights, averaging thin and thick components when present, for 2631 edge-on disk galaxies with $M_*>10^{10} M_\odot$ at $0< z < 3.5$ from the JWST COSMOS-Web survey. We show that dust extinction systematically overestimates scale heights at shorter rest-frame wavelengths, and therefore adopt a fixed rest-frame wavelength of 1 $μ$m. After further correcting for projection-induced bias using a new accurate method, we find that the median disk scale height increases from $0.56\pm0.03$ kpc at $z=3.25$ to $0.84\pm0.04$ kpc at $z=1.25$, and subsequently decreases to $0.67\pm0.06$ kpc at $z=0.25$. The disk length-to-height ratio remains constant at $2.7\pm0.2$ for $z>1.5$, but rises to $4.0\pm0.4$ at $z=0.25$. These results imply that the high-redshift progenitors of present-day thick disks were of intermediate thickness, neither thin nor thick, yet dynamically hot and dense. The observed radial variation of scale height is consistent with the artificial flaring expected from observational effects, disfavoring minor mergers as the primary mechanism of disk thickening. Instead, we suggest that the high-redshift intermediate-thickness disks were single-component systems that increased their vertical scale height through decreasing surface mass density and/or violent gravitational instabilities, eventually producing thick disks. Thin-disk growth begins at $z\approx2$ and dominates at $z\lesssim1$, yielding a vertically more compact system with decreasing scale heights from $z\approx1$ to $0$. The inferred thin-disk mass fraction increases from $0.1\pm0.03$ at $z=1$ to $0.6\pm0.1$ at $z=0$. Together, these findings reveal a continuous evolutionary link between high-redshift single-component disks and present-day thick thin disk systems.
2601.04988Jan 2026
View
Rapid emergence of overmassive black holes in the early Universe
The origin of supermassive black holes (SMBHs) remains a long-standing problem in astrophysics. Recent JWST observations reveal an unexpectedly abundant population of overmassive black holes at z>4-6, where the BH masses lie far above local scaling relations and not reproduced by current cosmological models. How such overmassive black holes form and rapidly grow within young galaxies has remained unclear. Here we present fully cosmological radiation-hydrodynamic simulations that, for the first time, self-consistently follow the birth, early growth, and emergent observable signatures of SMBHs in proto-cluster environments. We find that heavy seeds of order $10^6 M_\text{sun}$ naturally form, exceeding typical theoretical expectations by an order of magnitude. These seeds rapidly develop dense, optically thick disks whose strong electron scattering produces broad H$α$ emission comparable to that seen in little red dots (LRDs). Sustained super-Eddington accretion then drives fast growth to $\sim 3 \times 10^7 ~M_\text{sun}$ by $z \sim 8$. These results provide a unified physical scenario in which LRDs correspond to a short-lived, enshrouded phase of heavy-seed formation, naturally evolving into the overmassive quasars detected by JWST and ultimately the progenitors of today's SMBHs.
2601.04955Jan 2026
View2601.04953Triaxial Magnetars as Sources of Fast Radio Bursts
I suggest that some of the mysterious temporal properties of Fast Radio Bursts (FRB) may be explained if they are produced by dynamically triaxial magnetars. If the bursts are narrowly collimated along open field lines, then observed repeating FRB may be those in which the moment of inertia tensor is only slightly triaxial and the rotation axis, open field lines and radiation point nearly to the observer. Apparently non-repeating FRB may be triaxial with the direction of open field lines and radiation wandering across the sky, reducing their duty factors by several orders of magnitude. A slightly triaxial moment tensor in repeaters moves the line of sight into or out of the radiation pattern or within it, explaining periods of greater or lesser (or absent) activity, and making the probability of detecting a burst vary aperiodically. The dynamics of triaxial bodies may also thwart the coherent integration of gravitational signals from fast-rotating accreting neutron stars.
2601.04953Jan 2026
View2601.04939Multi-Messenger Studies with High-Energy Neutrinos and Gamma Rays: The WST Opportunity
The search for the sources of ultra-high-energy cosmic rays (UHECRs) using high-energy neutrinos represents a frontier in high-energy astrophysics. However, a critical bottleneck remains: the ability to rapidly survey the sizable sky areas defined by the localization uncertainties of neutrino detectors and to provide rapid spectroscopic classification of the multitude of optical transients found within them. By deploying a large field-of-view with high-multiplex Multi-Object Spectroscopy (MOS) on a large aperture telescope, one can instantaneously cover neutrino error circles, thus providing crucial spectroscopic classifications of potential counterparts discovered, for example, by the Vera C. Rubin Observatory (LSST) with unprecedented efficiency. Furthermore, simultaneous operation of a giant panoramic central Integral Field Spectrograph (IFS) would allow for detailed kinematic and environmental characterization of primary candidates. This facility would unlock deep synergies between next-generation neutrino telescopes (IceCube-Gen2, KM3NeT) and gamma-ray observatories (CTAO), transforming unique multi-messenger alerts into a comprehensive physical understanding.
2601.04939Jan 2026
View2601.04931Three-equation turbulent convection models in classical variables
Context. Turbulent convection models in nonlinear radial stellar pulsation models rely on an extra equation for turbulent kinetic energy and fail to adequately explain mode-selection problems. Since multidimensional calculations are computationally expensive, it is reasonable to search for generalizations of physically grounded 1D models that approximate multidimensional results with sufficient accuracy, at least in a given parameter range. A natural way of progressing from one-equation models is to use additional nonlocal equations. While these types of models also exist in the literature, they have not been adopted for this type of object. Aims. We aim to adapt the three-equation turbulent convection model from Kuhfuss to radial stellar pulsation modeling. Methods. We use a Reynolds-stress one-point closure approach to derive our extensions alongside the model, while using additional models from the literature to close the anisotropy and dissipation terms. Results. We provide five extensions to the original model. These include an enhanced dissipation correction to the mixing length, a local anisotropy model replacing eddy viscosity, a second-order correction for turbulent ion transport in the atmosphere (alongside opacity effects), and turbulent damping of entropy fluctuations and convective flux.
2601.04931Jan 2026
ViewPage 1 of 310
