Peer Revenue Impact factor: volatile

Papers

Fresh from arXiv, typeset properly. Highlighted words trade on the market.

Physics Astrophysics

astro-ph Watch Astrophysics (sign in to use)

Thursday 14 May

Astrophysics of Galaxies

Observation of spontaneous N-bearing PAH formation using ion trap: a new formation pathway in the interstellar medium

Nitrogen-bearing polycyclic aromatic hydrocarbons (N-PAHs) are key precursors to complex organic molecules in both the interstellar medium and the nitrogen-rich planetary atmospheres. Despite the recent detections of nitrogen-functionalized astromolecules, their formation pathways remain an open question. The discrepancies between their predicted and observed abundances point to unknown mechanism that govern their evolution in the astrophysical environments. Employing an ion trap technique in conjuction with electronic structure calculations, we unravel multiple barrier-less reactions between gas-phase pyrimidine cations (C$_4$H$_4$N$_2^+$) and acetylene (C$_2$H$_2$) which form a hitherto unreported endocyclic- N-PAHs (C$_8$H$_7$N$_2^+$). The present measurements on reactions involving a double-nitrogen subsituted aromatic heterocycle have implications to the astrochemistry of both the Titan's atmosphere and the interstellar medium.
Read the full abstract

All authors Siddhartha S. Payra, Pratikkumar Thakkar, Shiv Gupta, Ruth Ann Mathews, Yash Lenka, Saurav Dutta, Nihar Ranjan Behera, Krishna R. Nandipati, G. Aravind

High Energy Astrophysical Phenomena · Cross-listed from gr-qc

Dyonic black holes supporting nearly-black self-gravitating thin shells

It has recently been revealed that dyonic black-hole spacetimes of a quasitopological non-linear electrodynamic field theory may be characterized by discrete radial regions with the property $dg_{tt}(r)/dr=0$ in which spherically symmetric massive test shells (Dyson shells with negligible self-gravity) can be supported in static equilibrium states. In the present paper we prove that the dyonic spacetimes of the non-linear electrodynamic field theory may also be characterized by the presence of radial regions with the dimensionless property $d[r\cdot g_{tt}(r)]/dr\to0^+$ in which massive self-gravitating thin shells that are on the verge of becoming black holes can be supported in static equilibrium states. Intriguingly, it is proved that the discrete radii of these self-gravitating nearly-black Dyson shells are universal in the sense that they are independent of the masses of the central supporting dyonic compact objects.
Read the full abstract

Instrumentation and Methods for Astrophysics

A $Z_1^2$ framework for rotational-parameter estimation and uncertainty quantification in high-energy pulsars

We present a $Z_1^2$-based framework for estimating the spin frequency and frequency derivative of high-energy pulsars from Poisson-limited photon event lists. The key point is that the width of a coherent detection peak is not, by itself, the statistical uncertainty on the recovered rotational parameters. We develop and compare three computationally efficient estimators: segmented frequency regression, a coherent derivative scan, and a localized two-dimensional coherent fit. For sinusoidal signals, we derive the local form of the Z-squared response as a function of frequency and frequency derivative, and show that expressing the frequency at the midpoint of the observation removes the leading-order covariance between the two parameters. This gives simple uncertainty estimates in terms of the fitted peak amplitude and local widths, without requiring an exhaustive Monte Carlo simulation for each observation. We test these estimates with Monte Carlo simulations over a range of observing spans, signal strengths, grid resolutions, and good-time-interval structures, and show that the predicted uncertainties reproduce the run-to-run scatter of the recovered parameters in the tested regimes. We then apply the framework to AstroSat/LAXPC event lists for the Crab pulsar, Swift J0243.6+6124, and SAX J1808.4-3658. The results provide a practical and statistically motivated route to rotational-parameter estimation for targeted high-energy pulsar searches.
Read the full abstract

All authors Akshat Singhal, Rohit Nair, Devendra Sahu, Gayathri Raman, Suman Bala

Astrophysics of Galaxies

An Updated Characterization of Luminous Ly{\alpha} emitters at the End of Reionization

We present a multi-wavelength physical characterization of 14 luminous Ly$\alpha$ emitters (LAEs) at $z\approx6$, integrating deep ground-based Magellan/M2FS spectroscopy with heterogeneous JWST/NIRCam broad- and medium-band imaging. Identified via strong Ly$\alpha$ lines with extreme Ly$\alpha$ luminosities of ${>}10^{42.6}$ erg s$^{-1}$, the sample exhibits very large rest-frame equivalent widths (${\gtrsim}100$ Å) and steeply blue UV continua ($\beta_{\rm median}\simeq-2.2$, $-18.2>M_{\rm 1500}>-20.2$ mag). Crucially, the integration of NIRCam medium-band photometry (F410M) breaks the degeneracy between strong rest-optical nebular emission and Balmer breaks, resolving prior mass overestimations. The tightly constrained spectral energy distribution modeling demonstrates that these luminous LAEs tend to be unequivocally low-mass, ultra-young dwarf starbursts; half the sample is characterized by stellar masses of $M_* < 10^8 M_{\odot}$, ages $\lesssim10$ Myr, and negligible dust attenuation. We also map the production efficiency of ionizing photons and Ly$\alpha$ escape fractions ($f_{\rm esc}^{\rm Ly\alpha}$). The $f_{\rm esc}^{\rm Ly\alpha}$ values are exceptionally high, with a median of ${\gtrsim}40$%, increasing for the bluer UV continua. Finally, analyzing spatial offsets between the Ly$\alpha$ centroid and the stellar counterpart, we demonstrate empirically that internal dust content, rather than neutral hydrogen gas, dominate the suppression of Ly$\alpha$ radiative transfer. Our study reveals that strong Ly$\alpha$ emission of the luminous LAEs are generally attributed to both the vigorous starburst activities and the high $f_{\rm esc}^{\rm Ly\alpha}$. Resembling Lyman continuum leakers, these extreme dwarf systems function as highly efficient ionizing engines at the conclusion of the Epoch of Reionization.
Read the full abstract

All authors Yuanhang Ning, Zheng Cai, Linhua Jiang, Qiong Li, Si-Yue Yu, Xiaodi Yu, Zhen-Ya Zheng

Astrophysics of Galaxies

Reconstructing the Stripping History of the Sagittarius Stream with Neural Networks

The Sagittarius (Sgr) Stream is produced by the ongoing disruption of the Sgr dwarf spheroidal (dSph) galaxy and is thought to contain multiple wraps that were stripped during different pericentric passages. In this study, we introduce a neural-network–based method trained on $N$-body simulations to infer the stripping time of Sgr Stream stars directly from their phase-space coordinates. We combine spectroscopic data from SEGUE, APOGEE DR17, and LAMOST DR7 LRS with Gaia EDR3 astrometry and distance estimates from the latest StarHorse catalog to identify high-quality Sgr Stream members. Applying our method to these stars, we measure a clear metallicity gradient with stripping time, well described by a linear relation with slope $\sim 0.3~\mathrm{dex~Gyr^{-1}}$. We further predict the stripping times of globular clusters previously suggested to originate from the Sgr dSph. M 54, Terzan 7, Terzan 8, and Arp 2 exhibit stripping times consistent with being currently bound to the Sgr remnant. Pal 12, Whiting 1, and NGC 2419 are inferred to have been stripped $0.9 \pm 0.1$, $1.1 \pm 0.2$, and $2.1 \pm 0.2$ Gyr ago, respectively. For NGC 4147 and NGC 5634, whose membership in the Sgr system remains uncertain, our analysis suggests stripping times of $1.1 \pm 0.4$ and $1.1 \pm 0.1$ Gyr, respectively, if they are ultimately confirmed as genuine Sgr members. These results demonstrate that data-driven models of dynamical stripping histories offer a promising approach for reconstructing the formation and chemical evolution of the Sgr Stream.
Read the full abstract

All authors Jian Zhang, Cuihua Du, Mingji Deng, Zhongcheng Li, Haoyang Liu

Instrumentation and Methods for Astrophysics

Opportunities for Gravitational Wave Physics at the South Pole

Atom interferometers represent a promising approach for gravitational wave detection in the decihertz frequency band, complementary to existing light-based detectors. The South Pole offers unique advantages for such experiments: exceptionally low seismic noise, established infrastructure for large scientific projects, and a location that strengthens gravitational wave source localization through global triangulation. Here we discuss the scientific case and practical considerations for deploying a long-baseline atom interferometer at the South Pole, which has the potential to expand the global network of gravitational wave detectors while enabling precision tests of fundamental physics.
Read the full abstract

All authors C. A. Argüelles, M. DuVernois, P. W. Graham, T. Kovachy, J. Mitchell

High Energy Astrophysical Phenomena

The compact neutron star in 4U 1746-37 revisited: Reassessing the mass and radius

A recent analysis of photospheric radius expansion X-ray bursts from the low-mass X-ray binary 4U 1746-37 reported unusually small mass and radius estimates for the neutron star, suggesting it could be a quark star or quark-cluster star. Here, we propose an alternative interpretation: the star's mass and radius could be underestimated from significant blocking of the X-ray flux. Significant blocking factors ($\mathcal{B} \gtrsim 0.8$, reducing the observed flux to $\sim17\%$ of the intrinsic emission) permit neutron star parameters consistent with the canonical values: $M = 1.59 \pm 0.69 M_{\odot}$, $R = 13.0 \pm 5.45\,\mathrm{km}$, or $M = 2.12 \pm 1.08 M_{\odot}$, $R = 9.80 \pm 4.13\,\mathrm{km}$. The blocking factor, which varies with the photospheric radius, provides a natural explanation for the anomalously large peak-to-touchdown flux ratio ($\sim2.0$) and highlights the importance of accounting for geometric system configuration in neutron star mass–radius estimates.
Read the full abstract

All authors Kwang Hyun Sung, Myungkuk Kim, Young-Min Kim, Kyujin Kwak, Chang-Hwan Lee

Earth and Planetary Astrophysics

Dynamical Evolution of V-Shaped Collision Debris

Catastrophic collisions between proto-satellites have been proposed as a possible origin of Saturn's rings. This argument relies on the concept of the equivalent circular orbit. Here, we re-examine the post-impact dynamical evolution of collision debris using analytical arguments and $N$-body simulations with fragmentation. We focus on the long-term evolution of debris distributed in a broad V-shaped region in the $a$–$e$ plane, with two arms for particles sharing a common collision radius. Because particles on the two arms possess significantly different angular momenta, inter-arm collisions dominate the evolution and drive behavior fundamentally different from the simple circularization assumed in the equivalent circular orbit approach. As a result, the classical equivalent circular orbit concept cannot predict the long-term fate of collision debris. Both our analytical framework and $N$-body simulations show that, although some debris initially passes within the Roche limit on eccentric orbits, successive collisional evolution drives the particles approximately along the original V-shaped constraint curves toward the apex of the V-shape, i.e., the original collision radius. Instead of spreading inward to form a ring, the debris converges and reaccretes near the original collision location. We therefore conclude that catastrophic proto-satellite collisions do not produce massive Saturnian rings. Rather, the debris evolves toward reaccretion into a new generation of satellite-sized bodies near the impact radius. These results fundamentally revise the dynamical interpretation of collision-generated debris and establish a more general framework applicable beyond the Saturnian system, including other planetary ring systems and debris produced during planet formation.
Read the full abstract

Astrophysics of Galaxies

A new sample of Little Red Dots at $z<0.45$ in DESI DR1: Broad Balmer lines, low ionization spectrum and no variability

JWST has unveiled an abundant population of compact broad-line emitters largely at $z\gtrsim4$, the Little Red Dots (LRDs), which might represent a previously unprobed supermassive black hole evolution channel predominant at high redshift. However, the LRDs have remained mostly elusive at lower redshift ($z\lesssim2$) where detailed studies are possible from ground-based observatories. We searched for low-redshift LRDs in the Dark Energy Spectroscopic Instrument (DESI) survey. Our search is primarily based on emission line properties, as opposed to earlier approaches that searched for compact sources with specific photometric spectral energy distributions. We report the discovery of eight LRDs at $z=0.2-0.45$, which show spectral features akin to the high-redshift LRDs in the rest-frame optical. The sources are characterized by broad Balmer lines, steep Balmer decrements, compact morphologies, Balmer absorption features and/or strong He I emission, but weak or absent He II, [Ne V] or other high excitation lines typical of Type I AGN. For 7 out of 8 sources, we retrieve dense-cadence light curves from time-domain surveys and for most sources we find weak to no intrinsic variability ($0.0-0.1$ mag) over $4-17$ years in the rest-frame. We also highlight the identification of a quasar with similar Balmer line profiles as LRDs, but shows differences in Balmer decrement, significant variability, and high-ionisation lines. Given the effective volume $4.9{\rm Gpc^3}$ covered by DESI DR1 at $z<0.45$, our sample corresponds to a number density of $1.6\times10^{-9}$Mpc$^{-3}$, indicating a number density $\sim$10,000 times lower than in the first billion years of cosmic time. We find a dearth of luminous and red LRDs at $z<1$ compared to higher-redshift, which could suggest lower gas feeding rates of LRD activity due to higher metallicities at later cosmic epochs.
Read the full abstract

All authors Kevin Park, Alberto Torralba, Jorryt Matthee, Sara Mascia, Zoltán Haiman, Rohan P. Naidu, Anna de Graaff

Solar and Stellar Astrophysics

The Close Binary V486 Carinae

The hitherto neglected close binary V486 Car is studied with the aid of newly applied satellite photometry (HIPPARCOS and TESS), high dispersion spectrometry (HERCULES) and ground-based B and V photometry. While the sinusoidal light variations are suggestive of a near-contact system, the stars have only shallow eclipse, so highly confident parametrization becomes challenging. We find: $M_1 = 2.1 \pm 0.1$, $M_2 = 0.4 \pm 0.1$; $R_1 = 3.20 \pm 0.02$, $R_2 = 1.48 \pm 0.01$; (${\odot}$); $T_{e1} = 10000 \pm 500$, $T_{e2} = 6200 \pm 200$ (K); distance = 162 $\pm$ 12 (pc). New times of minima for V486 Car have been examined, including recent observations from TESS. The role of the relatively significant O'Connell effect is examined. As well as the conspicuous asymmetry from the main effect of about 0.036 mag (V), a jitter, with amplitude of about 0.005 V mag and quasi-period of order $\sim$ 10 d is noticed. There is a tendency for such photometric excursions at one maximum to precede those at the other. As well, the O – C data indicate the presence of a low mass star $\sim$0.3 M$_{\odot}$ in an orbit separated by a few AU from the close binary. More accurate and plentiful spectroscopic data would be requisite for further investigations. A brief discussion reviews possible approaches to understanding the system in the context of near-contact binary scenarios.
Read the full abstract

All authors Ahmet Erdem, Volkan Bakis, Burcu Ozkardes, Edwin Budding, Mark G. Blackford, Tom Love, Michael D. Rhodes, Timothy S. Banks

Wednesday 13 May

High Energy Astrophysical Phenomena

The Phenomenological Nature of Quasar-type Blazars (BZQ). I. Revisiting the Flat-Spectrum Paradigm

We reevaluate 610 sources classified as Flat-spectrum radio quasars (FSRQs) in the 5th edition of the Roma-BZCAT. Optical spectra from SDSS DR16 confirm broad emission lines within $0.11 \leq z \leq 5.28$. To assess their blazar-like behavior, we combine ZTF optical variability with radio morphologies from FIRST, LOFAR, and VLBI, defining Confirmed, Possible, and Non-Confirmed BZQs. Rest-frame 1.4–10 GHz radio spectra were homogenized and fitted with error-weighted power laws. We show that the scheme of Park et al. (2013) often misclassifies nearly flat spectra as inverted and some prominently steep spectra as flat. Using the individual uncertainty $\sigma_{\alpha,i}$, we classify spectra as flat if $|\alpha_i| \leq 2\sigma_{\alpha,i}$, prominently steep if $\alpha_i > 2\sigma_{\alpha,i}$, and inverted if $\alpha_i < -2\sigma_{\alpha,i}$. This source-by-source criterion, intended as a phenomenological classification for this sample, better reflects the observed spectral shapes and confirms that most BZQs are consistent with being flat within measurement precision, although a non-negligible fraction departs from strict flatness. We also classify full spectral morphologies as power-law, peaked, restarted-peaked, or inverted-peaked, associated with distinct jet processes and activity cycles. About 60% of the sources with at least two decades of frequency coverage exhibit restarted-peaked spectra, suggesting recurrent jet activity. The observed diversity indicates that the label "Flat-spectrum radio quasar" does not fully describe this population, and that the more general term BZQ may better reflect its phenomenological diversity.
Read the full abstract

Solar and Stellar Astrophysics

The Distribution of Blue Straggler Stars in the Color-Magnitude Diagrams of Old Open Clusters

We examine the blue straggler star (BSS) populations of six old ($\geq$4 Gyr) open clusters: M67, NGC 188, NGC 6791, Berkeley 32, Berkeley 39, and Trumpler 19. We find that 50% of BSSs have color-magnitude diagram (CMD) locations corresponding to single stars in the final third of their main-sequence lifetimes. This build-up of BSSs near the terminal-age main sequence (TAMS) is primarily, but not solely, driven by more massive BSSs. Eleven of the BSSs have white dwarf companions with measured cooling ages; their evolution age distributions indicate that more massive BSSs typically form far from the zero-age main sequence, whereas lower mass BSSs can form at every evolutionary age. We show that inferred core helium amounts (above primordial) of late-evolution-age BSSs correspond to the core helium fused by cluster main-sequence stars near the turnoffs. We also find that the masses of asymptotic giant branch (AGB) mass-transfer BSSs require evolved main-sequence accretors and conservative mass transfer. These findings indicate that helium enrichment of progenitor accretors leads to the prevalence of BSSs near the TAMS. We further classify the evolutionary stages of the progenitor donors in M67 and NGC 188 and find mass transfer during the AGB accounts for at least half of the BSSs. We trace how the main-sequence binary population of NGC 188 evolves, and find that only 30-40% of interacting binaries create BSSs and that progenitor orbits must change to match current BSS periods.
Read the full abstract

Instrumentation and Methods for Astrophysics

CIAO: Chandra's Data Analysis System for X-Ray Astronomy and Beyond

The Chandra Interactive Analysis of Observations (CIAO) software, developed by the Chandra X-ray Center, has been the data analysis package for the Chandra X-ray Observatory since its launch in 1999. Over nearly three decades, CIAO has grown from a small software suite into a widely used system for X-ray data analysis and beyond. CIAO provides tools for calibration, spectral, imaging, and timing analysis, together with high-level scripts and the \sherpa\ modeling and fitting application. Its modular design and unified data model allow users to build flexible analysis workflows while maintaining consistency with the Chandra data processing pipeline. Visualization capabilities are provided through integration with SAOImageDS9 and Python-based tools, and simulation components such as ChaRT and MARX extend the analysis environment to include detailed modeling of instrumental effects. In this paper we describe CIAO's design, evolution, and capabilities after 25 years of Chandra operations. We also describe its core architecture, scripting environment, modeling, visualization tools, simulation components, and testing infrastructure, as well as the documentation and user support system that have contributed to its widespread use. CIAO's continued development and broad adoption highlight its important role in X-ray astronomy and its usefulness in multiwavelength astrophysical research.
Read the full abstract

All authors Antonella Fruscione, Jonathan McDowell, Douglas Burke, Mark Cresitello Dittmar, Ian N. Evans, Janet D. Evans, Kenny Glotfelty, Hans Moritz Guenther, David Huenemoerder, William Joye, Nicholas P. Lee, Warren McLaughlin, Joseph B. Miller, Melania Nynka, David A. Principe, Aneta Siemiginowska

Instrumentation and Methods for Astrophysics

Scattered light noise at LIGO Livingston Observatory during O4

Scattered light is one of the most common sources of noise in the LIGO gravitational wave detectors. Light scattering is a highly non-linear process through which motion at low frequencies gets up-converted and creates noise in a higher frequency band in the detector data. From the beginning of the fourth observation run, many glitches appeared in the data of LIGO Livingston detector in the frequency range 10-40 Hz, and the morphology of these glitches suggested that they were produced by scattered light. From our analysis, we identified two different populations of scattered light glitches, one group having higher SNR than the other. The glitches of the high- SNR group were solely modulated by microseismic ground motion (ground motion in 0.1-1.0 Hz) and in this paper, we present models of possible coupling mechanisms for these glitches. We also present results of a statistical correlation analysis based on our models, which indicates that the microseismic ground motion at the corner station along the X direction is the one most correlated with the noise which create these high SNR glitches. After installing baffles very close to the test mass mirrors, we have noticed a significant reduction in the rate and SNR of these glitches. The low-SNR glitches were primarily modulated by high frequency (10-30 Hz) vertical ground motion at the corner station, and this motion was coupling through a specific vacuum chamber at the corner station. After installing an additional seismic isolation platform in that vacuum chamber, these glitches have disappeared.
Read the full abstract

All authors Debasmita Nandi, Anamaria Effler, Siddharth Soni, Tabata Aira Ferreira, Robert Schofield, Huyen Pham, Timothy O'Hanlon, V. V. Frolov, Gabriela González

Solar and Stellar Astrophysics

Toward a Comprehensive Grid of Cepheid Models with MESA. IV. Modest Effects of Rotation on Blue Loops

Evolutionary tracks for $2-8M_\odot$ stars, with metallicities of $Z=0.014$, $0.006$, and $0.002$, including rotation, are computed with Modules for Experiments in Stellar Astrophysics (MESA). We study how rotation impacts the evolutionary properties of classical Cepheids. We examine whether rotation can offer a plausible explanation for the mass discrepancy problem when it is included in the evolutionary code using the fully diffusive approximation for rotationally induced mixing processes. We find that rotation barely influences the appearance and luminosity levels of the blue loops. While luminosity increases with increasing initial rotation rate, the increase does not exceed 0.04 dex, a fraction of the increase resulting from including the main sequence (MS) core overshooting of $0.2H_p$. As a consequence, rotation alone cannot resolve the mass discrepancy problem without simultaneously requiring significant MS core overshooting. Similar to the mass-luminosity relation, the period-radius and period-luminosity relations are barely affected by rotation, while the period-age relation predicts Cepheid ages to be only a few per cent longer compared with models without rotation. The predicted surface rotational velocities are too large compared with observations. These results are in contrast with those obtained with the Geneva code, which implements rotational mixing using the advective-diffusive scheme. In that approach, the luminosity levels of the loops are significantly higher, their luminosity extent increases, and the predicted rotation velocities are lower, compared with MESA models. The differences between the two approaches arise from significantly more efficient rotation-induced mixing during the MS evolution in models computed with the advective-diffusive scheme.
Read the full abstract

Solar and Stellar Astrophysics

Direct Evidence of Non-Ideal Dissipative Dynamics in Solar Wind Magnetic Switchbacks

Magnetic switchbacks, large-amplitude, localized Alfvenic like rotations of the solar wind magnetic field, have been the subject of intensive investigation, with approximately 200 refereed papers published in the last decade. Yet, fundamental controversies persist regarding whether switchbacks can be described with Ideal MHD (magnetohydrodynamic) physics or Hall-MHD physics and whether their origin is at the solar surface or in the solar wind. To settle these controversies, we present Parker Solar Probe electric field measurements between 13 and 40 solar radii, which show that switchbacks have non-zero electric fields in the plasma frame, a finding that definitively settles the physics controversy by proving that switchbacks are Hall-MHD, not Ideal MHD, structures. Along with these electric fields, there are enhanced Poynting vectors having three components with similar magnitudes that exist only inside the switchbacks. These facts contradict the view of switchbacks as simple outward-propagating pulses. Together, they resolve one controversy by showing that switchbacks in the young solar wind are a non-MHD process. They contribute to the second (source) controversy by identifying switchbacks as sites of active, in-situ, evolution. These findings provide a new framework for understanding energy transport and dissipation in astrophysical plasmas.
Read the full abstract

All authors Forrest Mozer, Oleksiy Agapitov, Kyung-Eun Choi, Richard Sydora

Earth and Planetary Astrophysics

Coupled Photochemical-Climate Modeling of Plausible Tenuous Outgassed Atmospheres on the TRAPPIST-1 Planets

Available JWST observations TRAPPIST-1 system have suggested that several of the planets are likely airless, or possess a very tenuous atmosphere. However, the high atmospheric escape rates expected for these planets suggest that any tenuous atmosphere must be replenished by constant outgassing, and past studies on modeling potential atmospheres for the planets have not widely considered surface pressures <1 bar. Here, we show that tenuous atmospheres on the TRAPPIST-1 planets are likely possible, supported by constant plausible rates of water and/or CO$_{2}$ outgassing against assumed high escape rates (up to ~10$^{30}$ s$^{-1}$). We use a coupled photochemical-climate model and sample from a broad phase space of outgassing, surface deposition, and top-of-atmosphere escape rates to test hundreds of atmospheres per planet. Critically, our model also allows surface pressure to vary based on the balance of sources and sinks. We find that 6 different compositional archetypes are generated via H$_{2}$O and/or CO$_{2}$ outgassing across our phase space, and atmospheres commonly fall between 10$^{-4}$ – 1 bar. We find that potentially habitable surface environments are possible for TRAPPIST-1d and e at pressures between 0.05 – 2 bar and 0.5 – 1 bar, respectively. Where possible, we compare our models to JWST observational data for TRAPPIST-1b, c, d, and e; all atmospheres found in this study for these planets match available transmission data to <3$\sigma$. However, emission data are consistent with atmospheric outcomes constrained to thin O$_{2}$-dominated compositions for TRAPPIST-1b ($\lesssim$0.01 bars) and c ($\lesssim$0.2 bars), which may or may not contain trace SO$_{2}$.
Read the full abstract

All authors Megan Gialluca, Victoria Meadows, Andrew Lincowski, Trent Thomas, Parker Hinton, David Brain, David Crisp

Cosmology and Nongalactic Astrophysics

Single field matter bounce with dark energy era: comparison with CMB Planck 2018 data and best fit parameters

In this work, we perform Markov Chain Monte Carlo (MCMC) analyses using the Planck 2018 cosmic microwave background (CMB) datasets, including temperature, polarization, and lensing, in order to compare matter bounce models with observational data. The particular model we considered contains a scalar field with an exponential potential, which behaves as dust in the asymptotic past of the contracting phase, it realizes a quantum bounce, and then behaves as a transient dark energy field at large scales in the expanding phase. The parameter $\lambda$ appearing in the exponential potential is directly related to the model's scalar spectral index, $n_s$, which is set free in the MCMC analyses, as well as the deepness of the bounce, which controls the amplitude of the power spectrum. We provide constraints on the cosmological parameters and compare the model's performance against the standard inflationary $\Lambda$CDM scenario. Our results indicate that Planck data alone cannot favor one model with respect to the other, showing that the model we investigate can be a viable alternative to inflation.
Read the full abstract

Solar and Stellar Astrophysics

SN2023ixf: ultraviolet-to-infrared radiative-transfer modeling of the nebular-phase evolution until 1000 days

We present non-local thermodynamic equilibrium radiative-transfer modeling of SN2023ixf during the nebular phase out to 1000d, using the same ejecta that matched its photospheric evolution, namely a partially stripped red-supergiant star of initially 15Msun whose terminal explosion yielded ejecta with 7-8Msun, kinetic energy of 1.2e51erg, and 56Ni mass of 0.05Msun, augmented with a cold dense shell (CDS) of 0.2Msun at 8000km/s. Interaction with circumstellar material persists at all epochs, powering the ultraviolet (UV) flux at all times, but dominating the optical only after ~600d. Matching the V-band light curve requires invoking both enhanced gamma-ray escape and dust formation after ~200d, first in the CDS and eventually in the inner ejecta as well. Depending on where they form relative to the dust, emission lines are uniformly attenuated or skewed with a blue-red asymmetry. Our models suggest a rising dust mass (chosen as an C-rich and Si-rich mixture) in the CDS and inner ejecta, possibly reaching 1e-4Msun at 700d, while an external cold dust component is required to match the mid-infrared emission. The UV radiation, largely unaffected by dust, is influenced by the emission and absorption from Fe lines, together with strong, blueshifted emission from Lyalpha and MgII2800, both present at >~200d and with a strengthening fractional flux thereafter. Optical-depth effects play a critical role for the UV flux, and most notably on Lyalpha whose strength depends strongly on the CDS structure (mass and extent) and the treatment of power injection. The CDS is continuously slowing down from 8000km/s at 112d to ~6500km/s at 998d, suggesting a growth in mass of several 0.1Msun. SN2023ixf shares many similarities with SN1993J at 1-3yr, but it is eventually fainter due to dust extinction and cooler (i.e., weak [NII] and no [OIII] lines) likely as a result of greater CDS and ejecta masses.
Read the full abstract

All authors Luc Dessart, Wynn V. Jacobson-Galan, K. Azalee Bostroem, Alexei V. Filippenko, WeiKang Zheng, Thomas G. Brink, Lluis Galbany, Claudia Gutierrez, Stefano Valenti

High Energy Astrophysical Phenomena

Planets in Pulsar Winds

Planets around pulsars were the first discovered exoplanets, found thanks to the extremely precise pulsar timing. Here we suggest that they could also be found through the radio emission generated by the pulsar-planet magnetospheric interaction. We present the results of special relativistic numerical simulations of planets in a pulsar wind of velocity $v=0.985~c$, corresponding to a Lorentz factor $\gamma=5.795$. Planets, modeled as a perfectly conducting solid surface in an external magnetic field originating from the pulsar, produce radio emission in the extended magnetic structure on the planet's nightside. We find that the planet around a known pulsar, PSR J0636+5129 b, could be detected via its radio emission. We outline the observational prospects for such objects.
Read the full abstract

All authors T. Kaister, S. Andrés Joya Méndez, P. Marmat, M. Čemeljić, M. Velli, J. Varela, M. Falanga

Solar and Stellar Astrophysics · Cross-listed from physics.plasm-ph

Boris and Exponential Integrators in the Theory of Particles Interacting with Magnetic Turbulence

The interaction of electrically charged particles with magnetic fields is a fundamental problem in several areas of physics. An example is the motion of energetic particles through a magnetized plasma. The most accurate and reliable way to explore theoretically the interactions between particles and fields is via test-particle simulations. In such simulations one creates the turbulent magnetic field and solves the Newton-Lorentz equation numerically by employing an integration scheme. In the current article we discuss exponential integrators and derive systematically from this the Rodrigues scheme as well as the famous Boris integrator. For an approach where one creates the magnetic field anew at each time step, both integrators are overall comparable. In theory the Rodrigues approach should be more accurate due to the fact that the occurring matrix exponential is evaluated without further approximations. Practically, both methods provide very similar results. It is argued in the current article that a Rodrigues based integrator is a very strong alternative because for the specific problem discussed here, it does not require longer computing times.
Read the full abstract

Cosmology and Nongalactic Astrophysics

The Magnetic Origin of Primordial Black Holes: Ultralight PBHs and Secondary GWs

Ultralight primordial black holes (PBHs) provide a compelling window into early-Universe cosmology. Following our earlier work, we explore a mechanism for the formation of ultralight PBHs sourced by primordial inflationary magnetic fields, without invoking an ultra-slow-roll phase of inflation. We propose a magnetogenesis model in which large curvature perturbations are induced at small scales, leading to the efficient production of ultralight PBHs across a broad mass spectrum. We analyze the phenomenological implications of these ultralight PBHs for early-Universe cosmology, particularly during reheating. We compute the resulting stochastic gravitational wave (GW) background generated by both the electromagnetic spectrum and evaporating PBHs, which exhibits distinctive features tied to the underlying magnetogenesis model parameters. Our results demonstrate that inflationary magnetic fields can serve as a viable and testable origin for ultralight PBHs, opening new avenues for probing the interplay between inflation, magnetogenesis, PBHs, and primordial gravitational waves.
Read the full abstract

Solar and Stellar Astrophysics

Spectral Disentangling Reveals Deep CNO-cycle Exposure in ET Cru

Binary stars undergoing mass transfer provide unique laboratories for testing stellar evolution. Here, we present a comprehensive photometric and spectroscopic analysis of the semi-detached system ET Cru. Using spectral disentangling, we independently determined the effective temperatures and chemical abundances of both components with high precision, including nine elements (eleven species). We find masses of $13.41\,M_\odot$ and $6.00\,M_\odot$ for the primary and secondary, respectively, with uncertainties of only $\sim$1.3%. The radii are $5.58\,R_\odot$ and $5.68\,R_\odot$, measured to within 0.4% and 0.5%. Surface gravities are constrained to better than 1%, while effective temperatures are determined to within 3-5%. The secondary exhibits extreme chemical anomalies, with severe carbon depletion and nitrogen enrichment far exceeding those reported in classical Algol systems. Multi-wavelength spectral energy distribution modelling yields a distance of $\sim$2.5 kpc, inconsistent with the $Gaia$ DR3 parallax, suggesting systematic astrometric uncertainties in the parallax distance. Together, these results establish ET Cru as a benchmark Algol-type binary, revealing direct spectroscopic evidence of deep CNO-cycle exposure in the donor and confirming the primary star as a rejuvenated gainer. ET Cru thus provides a chemically and dynamically illustrative case for understanding advanced binary interactions and the late evolutionary stages of massive-star evolution.
Read the full abstract

All authors Gökhan Yücel, Volkan Bak{\i}ş, Christian Nitschelm, Timur Şahin, Ferhat Güney

Astrophysics of Galaxies

Constraining the Galactic Center Dark Cluster with ELT/MICADO Observations

The Galactic Center hosts the densest known stellar environment in the Milky Way, dominated by the massive black hole Sgr A* and the surrounding nuclear star cluster. Theory predicts that this region should also contain a large population of stellar compact objects (SCOs) - black holes, neutron stars, and white dwarfs - forming a "dark cluster" whose distribution and properties remain observationally unconstrained. These unseen stellar remnants are central to questions of mass segregation, cluster dynamics, and the expected rate of extreme mass ratio inspirals (EMRIs) detectable by future gravitational-wave observatories including LISA. Current evidence for SCOs in the Galactic Center is indirect, relying on dynamical mass measurements, X-ray surveys, and a small number of transient sources. Direct detections remain elusive due to crowding, extinction, and the sensitivity limits of existing instruments. We explore how upcoming facilities, in particular the Extremely Large Telescope (ELT) with its first-light imager MICADO, can fundamentally transform this field. MICADO's combination of deep photometry, high spatial resolution, and precise astrometry will enable systematic searches for SCO-star binaries via photometric variability and orbital astrometric signatures, as well as direct detection of isolated accreting black holes interacting with the gas-rich Galactic Center environment. We outline the observational pathways, technical challenges, and expected sensitivities, showing that ELT/MICADO observations can provide the first quantitative constraints on the dark cluster population. Establishing these constraints will be pivotal for understanding the dynamical evolution of the Galactic Center, the role of compact remnants in nuclear star clusters, and the astrophysical context of gravitational-wave sources in galactic nuclei.
Read the full abstract

Solar and Stellar Astrophysics

Analysis of DQZ White Dwarf Evolution through Procyon

Procyon is a great system to probe stellar evolution of non-interacting binaries. We present an extensive grid of MESA (Modules for Experiments in Stellar Astrophysics) evolutionary tracks to constrain the evolution of Procyon A and B. We systematically vary the initial parameters of our grid anchored by precise dynamical masses and spectroscopic determinations of effective temperature (T$_{eff}$) and luminosity ($L$) to match the stars' positions on the H-R Diagram. Our goal is two-fold: (i) to quantify how the inferred system age and the progenitor mass of Procyon B depend on metallicity ($Z$), mixing length ($\alpha$), and core overshoot ($\beta$), and (ii) to determine the best fitting model of Procyon B within a model-based initial-to-final mass relationship (IMFR) for hydrogen-deficient white dwarfs. Our best-fit models reproduce the observed properties for both components, yielding $M_\mathrm{A}=1.487 \pm 0.095$ M$_{\odot}$, $M_\mathrm{B}=0.592\pm 0.082$ M$_{\odot}$, a system age of $2.23 \pm 0.90$ Gyr, and a white dwarf cooling age of $1.20\pm0.49$ Gyr for Procyon B, consistent with independent determinations. Our results point to higher core overshoot than the standard adopted range, with the best fits ranging from $\beta=0.5-1.0$. From our model grid, we map Procyon B to the initial-to-final mass relationship for H-deficient white dwarfs in the $1.9\!-\!2.6$ M$_{\odot}$ progenitor range. Additionally, we implement the accretion of heavy metals onto the surface of the WD and fit our isotopic abundances to spectroscopic observations. We outline the physics used in our analysis.
Read the full abstract

High Energy Astrophysical Phenomena

The Homogeneous MeerKAT and Swift/XRT X-ray Binary Radio:X-ray Plane

During the hard and quiescent spectral states in X-ray binaries, a non-linear correlation is observed between radio and X-ray luminosities, providing a valuable tool to probe the connection between accretion and jet production. This relation was originally thought to define a single 'standard' correlation spanning several orders of magnitude in X-ray luminosity, and was extended to active galactic nuclei by including a mass term. However, subsequent studies revealed a more complex picture, with some sources deviating from the standard correlation and instead populating distinct tracks. To date, all large studies of the radio:X-ray plane have combined data from multiple telescopes, introducing uncertainties due to differing instrument systematics and flux conversions between observing frequencies, thereby complicating comparisons and limiting constraints. ThunderKAT was a five-year programme on the MeerKAT radio telescope that monitored X-ray binaries in outburst, and ran alongside SwiftKAT which provided quasi-simultaneous Swift/XRT X-ray coverage. We present the full set of light curves from these programmes, comprising 948 radio and 1029 X-ray data points. An important finding is the frequent detection of unresolved radio emission during the soft state, likely dominated by previously launched jet ejecta. Using these data, we construct the largest, observationally homogeneous X-ray binary radio:X-ray plane to date. We relate these results to the physical mechanisms proposed to drive inter-source diversity, and outline directions for future observational and theoretical work. This paper is accompanied by a public data release of the ThunderKAT and SwiftKAT measurements and a compiled radio:X-ray plane, available through an interactive website.
Read the full abstract

All authors Justine Crook-Mansour, Rob Fender, Andrew Hughes, Sara Motta, Patrick A. Woudt, Arash Bahramian, Melania Del Santo, Zuobin Zhang, Thomas D. Russell, Jakob van den Eijnden, Joe Bright, David Williams-Baldwin, Francesco Carotenuto, Stéphane Corbel, Fraser J. Cowie, Alex Andersson, Noa Grollimund, James Matthews, Kelebogile Gasealahwe, Itumeleng Monaleng, Lauren Rhodes, Payaswini Saikia, Katie Savard, Evangelia Tremou, Xian Zhang

Astrophysics of Galaxies

Ionization Structure and Metal Enrichment of the Galactic Center Minispiral Observed with JWST

Sagittarius A* (Sgr A*) is the nearest quiescent supermassive black hole, and its proximity offers a unique opportunity to study its surrounding fuel supply. We leverage extensive spatial and spectroscopic information provided by the JWST/MIRI MRS instrument to disentangle mid-infrared ionized gas structures in the central 0.1 parsec of the Galaxy. The Galactic Minispiral's Bar and Northern Arm are revealed by their distinct morphological and kinematic signatures. Several compact (<1 arcsec) gas structures including X7 also appear within ~0.05 parsec of Sgr A* in the plane of the sky, moving with blue-shifted radial velocities >~600 km/s. Fine-structure-line measurements spanning ionization potentials ~7-55 eV are used to constrain the incident radiation field, metal abundances (neon, argon, sulfur, nickel, and iron), and dust depletion/destruction for each identified gas structure. Overall, the Minispiral gas shows mildly super-solar abundances (~1.8*solar), a Wolf–Rayet star-driven radiation field consistent with ionization parameters spanning -1.6<~log U<~-0.6, and significant nickel and iron dust destruction. In particular, the enhanced ionization state and strong dust processing in the high-velocity compact gas structures suggest that these features may experience localized fast shocks.
Read the full abstract

All authors Nicole M Ford, Mayura Balakrishnan, Sebastiano D. von Fellenberg, Daryl Haggard, Joseph M. Michail, Yuzhu Cui, Joseph L. Hora, Joey Neilsen, Giacomo Principe, Tamojeet Roychowdhury, Nadeen B Sabha, Howard A. Smith, Zach Sumners, S. P. Willner

High Energy Astrophysical Phenomena

Infrared Line Diagnostics Fail to Constrain Sgr A*'s UV Output

Sgr A*, the 4 x $10^6 M_{sun}$ supermassive black hole at the Galactic Center, exhibits frequent flaring with X-ray luminosities of $L_X \sim 10^{35}-10^{36}$ erg s$^{-1}$, while its ultraviolet (UV) emission remains unconstrained due to extreme extinction ($A_V \sim 30$ mag). We use JWST/MIRI time-resolved spectroscopy of the central 0.3" region to search for mid-infrared emission line variability driven by Sgr A* flares, comparing the results to CLOUDY photoionization models spanning flare luminosities of $L_{UV} = 10^{32}$-$10^{39}$ erg s$^{-1}$. We detect no continuum-correlated variability in any mid-infrared line, including [Fe II] (5.34 $\mu$m), [Ne II] (12.813 $\mu$m), [Fe II] (17.936 $\mu$m), and [S III] (18.713 $\mu$m) over lags of $\sim1$-8~hr. Despite expectations of a flare-driven response, we show that the lack of variability is consistent with the physical conditions in the spatially extended line-emitting gas, where light-crossing timescales of $\sim0.1$-$10$ days and recombination and cooling timescales much longer than the variability timescales suppress any observable response to individual flares. We further find that the modeled infrared spectra are dominated by continuum emission rather than isolated line emission. The brightest predicted lines are intrinsically weak (lower than 10$^{-5}$ mJy), and their detectability is further reduced by the large kinematic broadening expected at the CLOUDY-predicted emitting radii, reducing their contrast against the continuum. Extending the analysis to higher-ionization mid-infrared and near-infrared lines does not improve sensitivity. These results demonstrate that infrared emission lines trace a steady-state radiation field rather than individual flaring events, and therefore infrared line diagnostics cannot be used to constrain the instantaneous UV flux of Sgr A*.
Read the full abstract

All authors Mayura Balakrishnan, Sebastiano D. von Fellenberg, Daryl Haggard, Joseph M. Michail, Nicole M. Ford, Joseph L. Hora, Laurent Loinard, Sera Markoff, Joey Neilsen, Giacomo Principe, Nadeen B. Sabha, Howard A. Smith, Zach Sumners, Shuo Zhang

Instrumentation and Methods for Astrophysics

Rapid and robust simulation-based inference for kilonovae

With the next generation of both electromagnetic and gravitational wave observatories beginning to come online, rapid analysis methods for kilonova data are becoming increasingly important in astronomy. Traditional Bayesian parameter estimation using Markov chain Monte Carlo (MCMC) is time-consuming and relies on explicit likelihood approximations that can break down when modeling uncertainties are significant. We develop a simulation-based inference (SBI) framework for kilonova parameter estimation using density-estimation likelihood-free inference. The framework uses a Gaussian process emulator trained on $\sim 1300$ POSSIS simulations. We demonstrate that SBI provides a rapid alternative to MCMC that is robust to likelihood misspecification. The standard Gaussian likelihood approximation fails to capture the non-Gaussian, correlated structure of emulator uncertainty; SBI learns this structure directly from forward simulations. Simulation studies show that the SBI method accurately recovers injected parameters, while the MCMC suffers from systematic bias caused by likelihood misspecification. This problem persists when analyzing AT2017gfo, where a subset of the MCMC posteriors pile up at prior boundaries and the SBI posteriors do not. The SBI framework infers a total ejecta mass of $\sim 0.087 M_{\odot}$ dominated by lanthanide-poor ejecta and excludes toroidal and peanut ejecta geometries at the 99th percentile for both components. The SBI framework generates $\sim 2 \times 10^{4}$ posterior samples in seconds.
Read the full abstract

All authors Stephanie M. Brown, Mattia Bulla, Hiranya V. Peiris, Nikhil Sarin, Daniel Mortlock, Stephen Thorp, Gurjeet Jagwani, Stephan Rosswog, Samaya Nissanke

Astrophysics of Galaxies

Lyman-alpha Radiation Pressure in Dense Star Clusters: Implications for Star Formation and Winds at Cosmic Dawn

Observations with the JWST in lensed fields have revealed that galaxies at cosmic dawn may concentrate their star formation in highly dense, compact, star clusters. The high columns and low metallicities encountered in their birth environments suggest that Lyman-alpha (Ly$\alpha$) radiation pressure may be crucial to their formation and evolution. In this study, we address this question by post-processing snapshots from radiation hydrodynamic simulations of dense star cluster-forming clouds ($\Sigma_*\gtrsim10^3{M_\odot{pc}^{-2}}$) with a range of dust abundances ($Z_d=0-0.1Z_{d,\odot}$) using the COLT Monte Carlo code. We infer that Ly$\alpha$ is likely to have mild (~10%) effects on the gas-to-star conversion efficiencies ($\epsilon_*\gtrsim60$%) for $Z_d\gtrsim0.01Z_{d,\odot}$, and even in dust-free environments, $\epsilon_*\gtrsim25$% - much higher than the <10% values typical of star-forming regions in the local Universe. This is because the densest filaments dominating stellar mass assembly ($n\gtrsim10^4{cm}^{-3}$) remain sub-Eddington ($f_{Edd}<1$). On the other hand, the bulk of the gas volume ($n\lesssim10^3{cm}^{-3}$) has $f_{Edd}>1$, with noticeable fractions having $f_{Edd}\gtrsim10$, implying that Ly$\alpha$ can launch dynamically significant winds from these systems rapidly ($\lesssim$4Myr), with possible implications for ionizing photon escape and galactic outflows. The Ly$\alpha$ force multiplier $M_F$ is highly sensitive to $Z_d$, with $M_F\lesssim3$ ($\lesssim 500$) for $0.1Z_{d,\odot}$ (dust-free) environments respectively. Nevertheless, Ly$\alpha$ dominates over UV and IR radiation pressure at all values of $Z_d\lesssim0.1Z_{d,\odot}$, by factors of ~3-500. Our results suggest that Ly$\alpha$ radiation pressure reinforces the emerging picture of locally efficient, bursty star formation accompanied by rapid outflows in galaxies at cosmic dawn.
Read the full abstract