2026-09-16 | | Total: 99
Galactic foregrounds represent a major contamination to the measurement of primordial $B$-modes from observations of the Cosmic Microwave Background polarisation. Even after the application of component separation algorithms, foreground residuals may potentially still bias the estimate of the tensor-to-scalar ratio $r$, causing a false detection. In this work, we present the methodology of a robustness test for the validation of an eventual detection of primordial $B$-modes, as obtained by a future, LiteBIRD-like satellite experiment. The goal of the test is to identify the foreground residuals contamination by looking for non-Gaussian properties in the CMB $B$-modes map, recovered through blind component separation algorithms. We adopt scattering transforms (ST) as our summary statistics sensitive to the non-Gaussian features of foreground residuals and to their correlation with foregrounds tracer maps. We characterise and validate the methodology on realistic sky simulations with different levels of foregrounds complexity. The proposed test is able to identify a bias on the tensor-to-scalar ratio of $\gtrsim 10^{-3}$ in $\sim 90\%$ of our simulations, with this bias value being of the same order of the accuracy targeted by LiteBIRD. Additionally, for our particular experimental configuration, the test is passed when the bias is lower than the sensitivity on the $r$ parameter, and no warning is raised. These results provide an important step forward in the development of statistical tools for validating future measurement of cosmological parameters, against foregrounds contamination.
Traditionally, the disc is believed to evolve under the influence of turbulent viscosity, although the importance of this mechanism has recently been questioned. Alternatively, magnetic disc wind or gravitational instability could be considered, but which of these three mechanisms contributes the most at a given evolutionary stage remains an open question. We present global numerical hydrodynamic simulations of protoplanetary disc evolution in the thin-disc limit on disc lifetime timescales, including all three mechanisms of mass and angular momentum transport. We calculate the gravitational, viscous, and magnetic torques to assess the contribution of each mechanism. We found that these transport mechanisms occupy distinct zones of influence in the disc depending on the stage of evolution: the dominance of disc self-gravity during the early phase is replaced by the prevalence of magnetic disc wind and viscosity, with the latter being especially relevant at the disc periphery. We show that a noticeable decrease in disc mass occurs only after the termination of the embedded phase. The suppression of disc viscous spreading is achieved only for the most intense wind; otherwise, the disc continues to grow in size till the end of simulations. The spiral structure contributes to the temporary retention of dust in the disc. However, this is insufficient to prevent the depletion of dust on long timescales. The results emphasize the importance of a comprehensive approach to long-term simulations of protoplanetary discs with no single mechanism of mass and angular momentum transport regarded as exclusive.
Recent observations have identified a new class of long-duration gamma-ray bursts (GRBs) likely produced by compact binary mergers. Using radio observations and lightcurve modelling, we investigate whether a magnetar remnant can power the long-lasting prompt emission observed in this merger-driven GRB sample. We derive constraints on the magnetars' rotational energy and magnetic-field strength, and assess whether the magnetar central-engine model is viable for the observed events. We conducted 2.1 GHz radio observations of seven nearby merger-driven GRBs ($z<0.25$) using the Australia Telescope Compact Array (ATCA), obtaining sensitive upper limits on late-time radio emission 520-6900 days post-burst. We modelled the expected radio light curves from magnetar-energized ejecta interacting with the circumburst medium and compared them to our observations. We separately tested the magnetar hypothesis against the extended gamma-ray emission using a fallback-accreting magnetar model. No radio counterpart is detected in any of the observed GRBs, with 3$σ$ flux density upper limits of 33-72 $μ$Jy at 2.1 GHz. Modelling the expected synchrotron emission from magnetar-energized ejecta interacting with the circumburst medium, we find that the radio non-detections remain compatible with energetic outflows in low-density environments. The analysis of the prompt gamma-ray emission provides complementary constraints on the spin period and magnetic field of the magnetar engine. Once fallback accretion and jet baryon-loading are considered, the parameter space is narrowly confined.
Measuring M dwarf ages is a central challenge in stellar astrophysics. These stars evolve so slowly that isochrone fitting gives little constraint, requiring alternative methods based on rotation, activity, and kinematics. Reliable ages help constrain the formation and evolution of planetary systems around M dwarfs, including the temperate sub-Neptunes now being characterised with the James Webb Space Telescope (JWST). We present an age analysis of six nearby M dwarf planet hosts (TOI-732, TOI-736, TOI-270, TOI-1468, TOI-1231, and K2-18) using lithium absorption at 6708 Å, rotation periods, and Galactic kinematics. The Li I 6708 Åfeature is absent above 3$σ$ in all targets, implying substantial depletion and ages exceeding 200 Myr. Rotation periods of 39-145 days yield ages of 2.8-8.6 Gyr, calibrated against open cluster M dwarfs and including the intrinsic dispersion of the rotation-age relation. Kinematic ages span 0.8-13 Gyr with uncertainties of 4-6 Gyr. These are less precise per star but help identify outliers such as TOI-1231 (~13 Gyr) and K2-18 ( ~1-2 Gyr) and place the stars in Galactic context. Rotation yields age constraints 2-6 times tighter than kinematics where periods are measured, and comparable constraints where periods must be inferred from chromospheric activity. Combining multiple indicators enables more secure classification of these stars' evolutionary states, providing reference points for future comparative studies.
Braking index measurements of hundreds of pulsars, despite observational caveats, show a clear trend: $n<3$ at early ages, $n>3$ at middle ages, and increasingly large, fluctuating values at larger characteristic ages. Crustal magnetic fields in isolated neutron stars evolve through Ohmic dissipation and non-linear Hall redistribution across scales, with additional early-time contributions from the chiral magnetic effect and post-burial re-emergence. Initial conditions are usually assumed dipole-dominated, whereas more complex and realistic fields, with energy spread more evenly over a broad range of scales, have not been explored in this context. We assess how the internal dynamics of initially non-trivial fields change the dipolar field strength and obliquity, and how their combined contribution to the braking index $n$, together with that of the alignment, compares with the observed values. To this end we perform long-term 3D magneto-thermal simulations, for configurations ranging from the dipole-dominated case to tangled, small-scale-dominated fields. Small-scale dominated configurations with weak dipolar components show rich dynamics: the dipolar mode is continuously fed by the small scales and may grow or decay depending on the initial spectrum, driving $n<3$ or $n>3$, while the obliquity varies on timescales as short as ${\cal O}$(kyr). These spontaneous changes often dominate the alignment torque and can account for the large values of $n-3$ of either sign observed at characteristic ages $τ_c\gtrsim10^5$ yr, unlike the dipole-dominated models. However, short-term magnetospheric fluctuations or superfluid-driven variations of the torque may still be needed to explain the measurements and the structured timing residuals quantitatively. Future work comparing these models with the observed timing properties could fine-tune the most promising initial configurations.
Recent measurements of the distribution of galaxies by the Dark Energy Spectroscopic Instrument (DESI) are in mild tension with a standard $Λ$-Cold Dark Matter cosmological model inferred from the anisotropies of the Cosmic Microwave Background (CMB) when extrapolated to low redshifts. Several solutions to explain this discrepancy in terms of new physics have been put forward, notably from new dynamics in a complex Dark Sector. Here we report on a resolution to this discrepancy that postulates that the inertial and gravitational masses of dark matter particles are not the same. Our best fit to CMB data combined either with Baryon Acoustic Oscillation (BAO) and Full Shape (FS) data from DESI DR1 (Data Release 1), or with BAO DR2 data alone, suggests the presence of a new long-range force in the Dark Sector with strength of $β\sim 0.5\% $ relative to gravity, at the $\sim 2.5σ$ level. A preliminary analysis combining FS DR1 and BAO DR2 data confirms the statistical power of DR2 and raises the evidence for dark forces to $\sim 3σ$, albeit ignoring the not yet available cross-correlation between the two datasets.
We discuss the interpretation of an observed running Hubble constant with the redshift, in terms of the possible underlying physical scenario. Assuming, like in the $Λ$CDM model, that the matter and dark energy components do not directly interact, we arrive at the conclusion that the physical content of the cosmological dynamics can always be represented with an evolutionary dark energy paradigm. We discuss both the case of the power-law running Hubble constant, and that of an optimized fitting model, both studied with the binned Pantheon Sample of the Type Ia Supernovae data. These two cases turn out to be associated with a phantom-like evolutionary dark energy.
We present a component-resolved analysis of the IR--X-ray connection in a sample of 104 spectroscopically confirmed Type~1 AGN. By combining \xmm, \swift/UVOT, and new high--angular resolution \jwst/NIRCam imaging for a subset of the sources, we constructed detailed SEDs in which the IR emission was explicitly decomposed into accretion-disk, torus, and polar-dust components. This approach enables us to isolate the intrinsic AGN-IR luminosities with minimal host-galaxy contamination, providing a clean assessment of the physical relationship between the accretion flow, circumnuclear dust, and the X-ray corona. By incorporating JWST/NIRCam observations into our SED decomposition, we significantly reduce host-galaxy contamination and isolate the nuclear infrared emission with greater fidelity, particularly at near-infrared (NIR) wavelengths. This improvement is most apparent at 1 $μ\mathrm{m}$, where sources with JWST coverage exhibit a markedly steeper relation between $ν$L$_ν$ (2 keV) and the NIR luminosity than those without JWST data. This behavior is naturally explained by starlight dilution, whereby host-galaxy emission dominates lower-resolution measurements.
The terrestrial ionosphere is a major limitation for high angular resolution imaging at decametre wavelengths. Its rapidly varying and direction-dependent phase errors can defocus sources, reduce recovered flux densities, and produce strong field-to-field variations in image quality. We investigate these effects using three overlapping pointings from the LOFAR Decametre Sky Survey (LoDeSS), observed under different ionospheric conditions. Using phase-calibration solutions on a bright source, we derive the spatial diffractive scale, the temporal decorrelation time, and the angular scale over which direction-dependent calibration solutions can be transferred before residual phase errors become significant. These metrics are reflected in the final image quality. The field with the weakest ionospheric phase structure has the sharpest compact sources and the most reliable flux-density scale, while the field with the strongest phase structure shows substantial source broadening and flux-density loss. The comparison indicates that spatial ionospheric structure is the dominant limitation in these observations. Calibration-derived ionospheric metrics, therefore, provide a practical way to diagnose and predict image-quality variations in decametre survey data.
It is important to locate water-ice in the lunar regolith to build futuristic, sustainable bases there. But communication delays from Earth, limited bandwidth, and power constraints demand an onboard, autonomous, time-efficient decision-making framework for locating volatiles such as water-ice in the lunar regolith. In this work, we present a framework, DewTwin-Coin, that can locate potential water-ice sites in real time without drilling. Its decision-making principle is based on a find-S learning algorithm that uses surface temperature and the Hydrogen-to-Oxygen intensity ratio at a location. The core idea is a physics-based rule $-$ water-ice is stable only if the temperature is below 110 K and the Hydrogen-to-Oxygen intensity ratio is between 1.7 and 2.3. We validate DewTwin-Coin on Chandrayaan-3's 3165 laser-induced breakdown spectroscopy (LIBS) elemental data with 387 Chandra's surface thermophysical experiment (ChaSTE) data. On a 16 GB RAM system, it classified all 3165 locations within 615 seconds, yielding no strong signatures of water-ice, which matches the Chandrayaan-3 in-situ analysis. This framework is very important for future lunar and other planetary missions in which an onboard rover can autonomously determine feasible drilling locations for water-ice without issuing false drilling commands in warm terrains. Thus, the mission's onboard setup can stretch its technical limits in both power and memory capacity.
Future near-infrared (NIR) astronomy missions dedicated to photon-starved science cases, in particular an all-sky astrometric survey extending the legacy of Gaia into the NIR, as identified within the ESA Voyage~2050 programme, require large-format detectors combining sub-electron effective read noise with very low dark current. Linear-mode avalanche photodiodes (LmAPDs) based on HgCdTe meet this need by multiplying the photo-generated charge by an avalanche gain before the readout stage, thereby reducing the relative contribution of the read noise. We report the first electro-optical characterization of an IBEX detector, a $2048\times2048$, \SI{15}{\micro\meter}-pitch HgCdTe LmAPD array developed by Leonardo with the European Space Agency and operated at \SI{80}{\kelvin} on a dedicated bench at CEA-IRFU. We first discuss the central difficulty of characterizing an APD array, the degeneracy between avalanche gain, quantum efficiency (QE), and conversion gain in the measured response, and we then separate the measurements into two categories. Among the quantities that are directly measurable without assumptions, we report a signal-to-noise ratio for a CDS measurement that exceeds that of a Euclid-type H2RG above \SI{8}{\volt} pixel bias under identical low-flux conditions, a photo-response non-uniformity stable at the $\sim$10\% level, and a quantum efficiency-to-excess-noise ratio (QEFR) of \num{0.40} at \SI{10.5}{\volt}. Among the quantities derived under an explicit set of assumptions, we obtain a conversion gain corresponding to a sense-node capacitance of \SI{27}{\femto\farad}, a quantum efficiency of $46\pm12\%$ at \SI{2.5}{\volt}, and an excess noise factor $F=1.15\pm0.14$ at \SI{10.5}{\volt}. These results establish IBEX as a promising European large-format detector for future ultra-low-flux NIR instruments.
Magnetars are young, isolated neutron stars that possess an exceptionally strong magnetic field, with surface dipolar strengths on the order of $10^{15}$ G. One of the plausible scenarios for generating such a strong field is an exponential amplification by a turbulent convective dynamo during the protoneutron star phase. However, the short expected duration of the convection ($\sim 10$ s) imposes a stringent constraint on the dynamo growth rate. We perform an extensive set of 82 three-dimensional convective dynamo simulations in the anelastic approximation and investigate the kinematic phase to quantify the dynamo growth rate $γ$. We find that $γ$ increases with both the magnetic Prandtl number ${\rm Pm}$ and the Rayleigh number ${\rm Ra}$, with the most unstable mode becoming highly non-axisymmetric and multipolar. We further observe a gradual transition from large-scale to small-scale dynamo as the magnetic Reynolds number ${\rm Rm}$ increases, resulting in a magnetic field that is predominantly concentrated at small scales. The trend remains unchanged when the outer magnetic boundary condition is varied. Since resolving the increasingly small scales becomes numerically impractical, we employ the theoretical small-scale Kazantsev dynamo model to explore the large ${\rm Pm}$ regime characteristic of protoneutron stars. The model qualitatively captures the growth rate behaviour observed in simulations and, upon extrapolation to the large ${\rm Pm}$ limit, indicates that $γ$ is only weakly dependent on the resistivity in a PNS. Under conditions relevant to the PNS, this model predicts a magnetic energy growth rate of the order of $\sim 1$ ms$^{-1}$.
GG Tau $A$ is a young triple system with the close pair, $Ab_1$--$Ab_2$, and $Aa$, surrounded by a massive circumtriple disk with a large inner cavity difficult to explain in a binary framework. We aim to determine the orbital architecture and individual stellar masses of GG Tau $A$ and assess how the available astrometric and disk-based constraints restrict the range of admissible solutions. We performed a joint fit using Oracle, developed specifically for hierarchical stellar systems. All astrometric measurements were placed in a common reference frame, since the historical wide-orbit astrometry is given relative to the unresolved photocenter of the $Ab$ subsystem. The fit included one new wide-orbit astrometric epoch and a prior on the total stellar mass derived from disk kinematics. We then applied, in post-processing, an additional geometrical constraint based on the observed center of the circumtriple disk. The fit yields orbital solutions compatible with the available astrometric and disk-based constraints and provides estimates of the individual stellar masses. The additional wide-orbit epoch only marginally reduces the range of admissible solutions. By contrast, the disk-center constraint leaves the favored orbital architectures largely unchanged but significantly tightens the stellar-mass partition. This yields posterior masses of $0.521^{+0.069}_{-0.051}$, $0.106^{+0.017}_{-0.013}$, and $0.79^{+0.10}_{-0.10}$ $M_\odot$ for $Ab_1$, $Ab_2$, and $Aa$, respectively. The reported values are posterior medians with 16th--84th percentile intervals. A joint treatment of the two orbital levels is required to recover a physically meaningful architecture and constrain the individual stellar masses. The resulting solutions provide a basis for future dynamical modeling of the circumtriple disk and for testing whether the observed cavity can further constrain the system architecture.
We use the proton spectrum resolved by LHAASO across the knee to test superheavy dark matter in the nucleon channel at PeV energies. Both the hardening at $0.34$ PeV and the knee at $3.3$ PeV are reproduced by nucleons from a relic of mass $5\cdot 10^{7}$ GeV decaying in the Galactic halo, with the mass as the only shape parameter of the spectrum and the background index profiled under an independent prior on the proton slope below $0.1$ PeV. It works equally well in all three hadronic-interaction-model reconstructions. On the other hand, the photons of the same cascade rule out this interpretation by about two orders of magnitude against the photon fraction of the cosmic radiation, and by about one order of magnitude against the LHAASO diffuse Galactic emission.
Fast radio bursts (FRBs) are bright, millisecond-duration transient radio signals of extragalactic origin with unknown source and emission mechanisms. While several thousand FRBs have been discovered at frequencies above 400 MHz, detections at lower frequencies have been exceptionally rare. Contributing to their rarity are the magnified dispersive delay, in the tens of seconds, and signal scattering that make low-frequency searches with established techniques computationally prohibitive. In particular, the only confirmed detection below 300 MHz came from a Low-Frequency Array observing campaign targeted at FRB 20180916B, a source already known to emit at higher frequencies. Here we present the first non-targeted detection of an FRB below 300 MHz, achieved in a wide-field, imaging-based survey with the Murchison Widefield Array. The implied FRB rate is consistent with that observed at higher frequencies, demonstrating that non-targeted wide-field surveys can uncover an unbiased, large low-frequency FRB population. This result challenges models invoking severe scattering or absorption to explain previous non-detections, and strengthens the case of low-frequency observations as promising avenue to study the local environment of FRBs.
PyExoCross is a Python-based spectroscopic post-processing framework for converting large atomic and molecular line lists into scientifically useful quantities, including partition function, specific heats, cooling functions, lifetimes, oscillator strengths, line intensities, stick spectra, and absorption and emission cross sections. It is intended for applications in astrophysics, planetary atmospheres, laboratory spectroscopy, and other high-temperature environments where large modern databases require efficient and reproducible computational treatment. The new release provides both a configuration-file-driven command-line interface (CLI) and a Python application programming interface (API) within a unified computational framework. The original input file workflow is retained, while the Python API enables direct use in scripts, notebooks, and automated pipelines. The version PyExoCross 2.0 also adds explicit support for non-local thermodynamic equilibrium (non-LTE) calculations, including two-temperature models and user-defined density and population treatments for both absorption and emission spectra. In addition, database compatibility has been expanded beyond ExoMol, HITRAN, and HITEMP line lists to include high-resolution molecular database ExoMolHR and atomic database ExoAtom. GPU acceleration is also introduced for computationally intensive intensity and cross-section calculations and is available for all currently supported database formats. These developments improve workflow flexibility, reproducibility, and integration with modern data-analysis environments, while extending the physical modelling capabilities of the code. PyExoCross therefore provides a more general and extensible platform for large-scale spectroscopic simulations based on modern atomic and molecular databases.
Numerous attempts to detect glycine, a simple amino acid and a key building block of proteins in all life on Earth, in the interstellar medium have so far been unsuccessful in spite of the fact that it is present in meteorites, comets and asteroids. Glycine is conformationally diverse, existing in eight different forms, and the aim is to show that interconversions, in the dark, are governed by quantum mechanical tunnelling at temperatures from 300 K down to 10 K and to compute the rate at which these changes occur resulting in recommendations to maximise a successful outcome.
Much of the ordinary matter in the present-day Universe still escapes direct detection. These missing baryons are thought to reside in a warm-hot intergalactic medium (WHIM) threading the filaments of the cosmic web, whose turbulent and thermal state remains poorly constrained yet underpins any attempt to observe it. Beyond its observational stakes, how this turbulence emerges, as gravitational collapse converts ordered inflow into a disordered cascade, is a question of structure formation in its own right. We characterise the emergence of turbulence in the diffuse gas of a cosmic filament, from its assembly to the present day, in the first simulation of the Web Evolution in Filament Targeted zoom simulations (WEFT) project. The cosmological zoom-in technique is applied to a single filament, evolved from z=63 to z=0 with the moving-mesh code AREPO, reaching a median gas cell size of ~8 kpc. We trace the filament assembly, the thermodynamic state of its diffuse gas, the generation of vorticity at its accretion shocks, and the growth of turbulent motions, then quantify the intermittency of the cascade through high-order velocity structure functions and their relative scaling exponents from extended self-similarity. The filament assembles by the hierarchical merging of several proto-filaments rather than by laminar accretion. Vorticity is seeded baroclinically where the rotating strands meet their accretion shocks obliquely. The flow evolves from a supersonic, shock-dominated state, statistically close to the bifractal Burgers limit, into a developed, mildly supersonic, intermittent cascade whose longitudinal exponents lie closest to the sheet-like She-Lévêque model at z=0. Turbulence is thus an intrinsic, quantifiable property of the diffuse gas of filaments. This pathfinder run provides the turbulent input for forecasts of WHIM observability.
X-ray emission in solar flares is produced by both multi-thermal plasma and accelerated electrons. Classical imaging approaches reconstruct X-ray intensity maps which contain contributions from multiple spectral components (e.g., hot, superhot, and non-thermal), and do not allow retrieving the morphology of the different components separately. We introduce a novel imaging technique, called "Count-based spectral Component Imaging (CCI)", to jointly reconstruct spatially resolved emission measure maps of the thermal components, and the electron flux distribution of the non-thermal component from data provided by the Spectrometer/Telescope for Imaging X-rays (STIX) aboard Solar Orbiter. We formulate a linear model linking the Differential Emission Measure (DEM), approximated by two thermal components, and the non-thermal electron flux to the observed counts. The resulting inverse problem is solved with the Richardson-Lucy algorithm. We apply CCI to STIX observations of SOL2024-10-01T22 and compare it with the previously developed Spectral Component Imaging (SCI) method, as well as classical imaging approaches. The reconstructed thermal and non-thermal components show good agreement with those obtained using SCI. This proof-of-concept study shows that CCI obtains results consistent with SCI but with fewer inputs. In contrast with SCI, CCI can also be applied to hard X-ray focusing optics imaging.
The chemical compositions of planets and their hosts stars are intrinsically linked, having formed from the same protostellar material. Characterising this relationship provides key constraints on the processes governing planetary formation and evolution. Focusing on host stars with near-solar chemical abundances limits our understanding of how planetary composition varies across a broader chemical parameter space. Expanding the sample to include stars with compositions markedly different from the Sun is therefore essential for building a complete picture of star--planet compositional connections. Here, we focus specifically on iron-poor hosts since they are more likely to be alpha-enhanced and represent some of the most chemically distinct stars relative to the Sun (i.e., thick disc stars). We present a sample of 45 stars hosting 64 planets across the super-Earth and sub-Neptune regimes, for which we homogeneously obtain new stellar parameters including abundances, re-derive planetary mass and radius, and model the bulk interior compositions of their planetary companions. No statistically significant evidence for a correlation between stellar and planetary composition was found in our sample. We suggest that this null result is primarily driven by the large uncertainties inherent to both compositional proxies, which may obscure an underlying relationship. Quantifying these uncertainties is therefore a critical step that previous studies have not fully addressed. Furthermore, even with improved precision, uncovering such a relation may require a higher-dimensional treatment that accounts for additional parameters such as planetary equilibrium temperature.
Paleo-detectors provide a unique avenue to reconstruct the multi-million-year history of cosmic-ray (CR) flux, preserving signatures of transient high-energy events such as nearby supernovae. This technique aims to use natural minerals as particle detectors, looking at the persistent damage tracks created by CR-induced nuclear recoils, accumulated over the minerals' geological lifespan, offering a geological archive of past particle fluxes. Building on our study of Messinian Salinity Crisis evaporites, which demonstrated that minerals with specific geological histories may enable the detection of primary CR flux variations, we have now expanded to diverse terrestrial records. This contribution presents our recent publication proposing olivine xenoliths from Auvergne, France, where eruption chronosequences could allow for the differentiation of CR flux scenarios over the last 50 kyr. This phenomenological work is supported by the INFN-funded PRImuS experiment, which aims to prove the efficacy of high-throughput optical microscopy and plasma etching to analyze these mineral targets. By validating theoretical track-length spectra and refining the estimates of experimental effects, PRImuS's goal is to establish paleo-detectors as a powerful tool for very-long-range time-domain astrophysics.
Laboratory simulations of meteor ablation provide a critical quantitative link between the chemical composition of meteoroids and their observed spectral features. In this work, we analyzed high-resolution Echelle spectra (wavelength range 380-780 nm) of 22 diverse meteorites from the dataset presented in our previous work (Matlovič et al., 2024), representing the largest collection of laboratory meteor analogs to date. Using a radiative transfer model assuming local thermodynamic equilibrium (LTE) and accounting for self-absorption in optically thick plasma, we derived plasma parameters and elemental abundances for both major (Fe, Mg, Cr, Mn, Si, Na, Ni, Li, and K) and minor (Co, Cu, and V) species. Comparison with known bulk meteorite compositions allowed us to validate the modeling approach and assess chemical biases resulting from laboratory-induced ablation. Our analysis suggested plasma temperatures between 5220 and 5810 K and revealed systematic discrepancies in the elemental abundances compared to the original chemical composition. Specifically, we observed a significant enhancement of volatile species (Na, K) relative to Fe, accompanied by a depletion of the moderately volatile element Mg, while refractory elements (Al, Ca, Ti) remained undetected in the plasma radiation. These trends are consistent with the equilibrium vaporization model and demonstrate that under the simulated entry conditions ($\sim$12 km/s at $\sim$80 km altitude), the ablation process is dominated by incomplete and fractional vaporization. We conclude that while laboratory spectra of plasma from ablated meteorites do not fully reflect the original bulk composition, radiative transfer modeling effectively characterizes the state of the radiating plasma, offering a more robust approach for interpreting compositional properties from meteor observations.
Disk fractions in star-forming clusters decline with stellar age, and are routinely read as the lifetime distribution of primordial protoplanetary disks. We argue that this reading is not warranted, because late infall can replenish or reform disks after collapse. To formalize this, we model the disk fraction as a superposition of primordial and environmentally replenished components, with region-to-region dispersion in the supply rate. Disk fractions cannot distinguish a long-lived disk from a short-lived one that is repeatedly replenished. In Sco--Cen both reproduce the data, the replenished model marginally better. Stellar age therefore need not equal disk age, and stars can host successive disk generations with distinct mass, orientation, and chemistry. A key test is whether regions of the same stellar age show a genuine spread in disk fractions, as environmental variation in the replenishment supply predicts.
Binding energies of molecules on ice mantles are important to understand the evolution of molecular complexity in molecular clouds. Binding energies are often computed using density functional theory (DFT) calculations, typically on either small amorphous ice clusters or crystalline slabs. Since these calculations require an accurate description of the electronic structure, hybrid functionals with dispersion corrections, basis set superposition error corrections and zero point energy corrections are typically employed. This, however, comes at a high computational cost, so most often small ice clusters are considered, frequently containing no more than twenty water molecules or so. While several recent studies have explored binding energy distributions, few have addressed the importance of substrate size. We here perform DFT calculations, using six different functionals, on ice clusters containing 10 to 100 H2O molecules to quantify the ice cluster size effect and to separate the electronic from the geometrical contribution. As probe molecules, we use CO, CO2 and NH3. These calculations demonstrate that, irrespective of the molecule and functional used, interaction energies only start to converge from thirty to forty water molecules onwards. The dispersion energy flattens out earlier, whereas induction and polarization effects require larger clusters to stabilize, particularly at structurally confined (cavity) sites. We conclude that ice cluster sizes of at least 30-40 water molecules are needed to obtain reliable binding energies, and is as important as the choice of the functional.
Empirical mass-radius relations have long been used in exoplanetology to study planet demographics and to predict masses and radii for objects where only one of these is known. Inferred transitions between planet classes can also provide empirical evidence in support of theory. However, mass-radius relations fitted in two-dimensions can be complicated when additional factors influence the planetary parameters, as is the case with the radius inflation of giant planets. In this work, we present an empirical mass-radius relation derived by fitting a broken power-law to a cleaned PlanetS sample from which inflated giants have been removed. We find that when inflated giants are excluded, the preferred model contains three breakpoints and four segments, showing the emergence of a distinct Saturnian regime between $54\pm3~{\rm M_\oplus} <M< 258\pm11~{\rm M_\oplus}$. In this regime, planets are still growing but at a lower rate than the Neptunian regime, consistent with the onset of gravitational self-compression. We also update the irradiation threshold for giant planet inflation, $S_{\rm thr}=99~\rm S_\oplus$ ($1.3~\rm ergs\,s^{-1}\,cm^{-2}$); below $S_{\rm thr}$, 90% of giant planets have a radius excess $Δ\log{R}<2σ$.