Peer Revenue Now accepting positions

Papers

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

Physics Astrophysics

astro-ph Watch Astrophysics (sign in to use)

Thursday 13 August

Instrumentation and Methods for Astrophysics

SNIFFLES I: Intended Emission, Unwanted Emission, and Unintended Radiation from Low-Earth Orbiting Satellites Impacting Radio Astronomy from 1-26 GHz

We present the first results of SNIFFLES, an ongoing observational programme to characterise intended emission, unwanted emission, and unintended radiation from NGSO systems across common radio astronomy receiver bands from 1-26 GHz (L, S, C, X, and K bands). Using the Australian 22 metre Mopra radio telescope near Coonabarabran, New South Wales, with followup observations from a single dish of the ATCA interferometer and its new BIGCAT backend, we conducted 4629 tracked observations of satellites from four NGSO constellations (Starlink, OneWeb, Amazon Leo, and Guowang), amounting to 375.9 hours of telescope time. Satellite identification was confirmed by correlating detected Doppler shifts with predicted ephemerides. We detected 2345 instances of intended emission, unwanted emission, and unintended radiation at 300+ unique frequencies from three of the four systems. The detections span all three interference classes affecting radio astronomy: (1) intended emission (including DTD), (2) unwanted emission (OOBE, including up to the fourth harmonic of the Starlink DTD signal at approximately 2.6 GHz, with the fourth harmonic detected near 10.5 GHz), and (3) unintended radiation from satellite platform electronics. Several detections fall within primary radio astronomy allocations, including 1613.19 MHz within the protected OH line band, and at 2690.76 MHz and 2700 MHz. At 2700 MHz, unintended radiation was detected in 76.9 percent of all observations of the relevant satellite version. ATCA followup measurements confirm flux densities up to eleven orders of magnitude brighter than typical astronomical sources, well in excess of levels that saturate radio astronomy receivers.
Read the full abstract

Solar and Stellar Astrophysics · Cross-listed from gr-qc

Temperature effects on white dwarfs in modified gravity

In this article we analyze the effects of a finite temperature equation of state on the equilibrium structure of white dwarfs in massive Brans-Dicke theory as well as the symmetron and dilaton screening mechanisms. We compute and present the numerically obtained mass-radius relation, effective gravitational constant as well as radial profiles of the scalar field, pressure and metric within the star. We show that assuming a non-zero temperature effectively results in a larger radius while leaving the total mass of the star essentially unchanged, and discuss the interplay between the effective gravitational constant, central density, and radius of the star.
Read the full abstract

High Energy Astrophysical Phenomena

Cross-spectral Analysis of the Type-C Quasi-periodic Oscillation Shoulder Component in GX 339-4

We revisit Rossi X-ray Timing Explorer (RXTE) observations of GX~339$-$4 during the rising phase of its 2006/2007 outburst and apply a joint power-density-spectrum (PDS)–cross-spectrum (CS) decomposition to the type-C quasi-periodic oscillation (QPO) region. Within this framework, the QPO region is described by a narrow QPO fundamental and a neighboring high-frequency shoulder, whose amplitudes and phase lags can be measured separately. The shoulder is first detected at MJD~54142.04, mainly through the imaginary part of the CS and a narrow local structure in the phase-lag spectrum, before becoming a resolved high-frequency shoulder in the PDS. It follows the QPO frequency evolution on the high-frequency side, with $R_{\nu}=\nu_{\rm sh}/\nu_{\rm QPO}\simeq1.04$–$1.18$. The QPO lag remains small, typically below $\sim0.17$~rad, whereas the shoulder carries a larger hard lag of $\sim0.5$–$0.8$~rad. Energy-resolved fits show the same separation: the QPO lag is close to zero or only weakly positive across most of the energy band, while the shoulder lag is systematically larger and generally increases with photon energy. The two components have broadly similar rms–energy shapes, although their relative strengths evolve during the observed sequence. Although the shoulder remains broad, with $Q\sim2$–$4$, its lag and rms–energy behavior resemble those of the type-B QPO detected shortly after our observations. This similarity raises the interesting possibility that the shoulder is related to an earlier, broader stage of the variability later seen as the type-B QPO.
Read the full abstract

All authors Haifan Zhu, Mariano Méndez, Pengcheng Yang, Pei Jin, Candela Bellavita, Federico García, Wei Wang, Diego Altamirano, Liang Zhang, Chenxu Gao, Xiao Chen

Astrophysics of Galaxies

Gal3D: Superellipsoid Modeling of Radial 3D Galaxy Structure in IllustrisTNG and EAGLE Simulations

Galaxy morphology and structure are key tracers of galaxy formation and evolution, making accurate measurements of intrinsic three-dimensional (3D) shape essential for linking morphology to galaxy assembly and for comparing numerical simulations. We present Gal3D, a framework that reconstructs smoothed density fields from particle data and quantifies the radial 3D structure of simulated galaxies by fitting superellipsoids to iso-density surfaces. The method recovers axis ratios, orientations, center offsets, and superellipsoid indices ($S_a$, $S_b$, $S_c$), enabling a flexible characterization of diverse galactic structures such as disks, classical bulges, box/peanut bulges, and triaxial components. Applying Gal3D to galaxies in the IllustrisTNG and EAGLE simulations, we find that the radial extent of flattened disk regions increases with stellar mass up to $M_{*,30}\sim10^{11}\,M_\odot$ and then declines sharply, with EAGLE galaxies showing a saturation at $M_{*,30}\sim10^{10.5}\,M_\odot$. The bar-related $ \varepsilon_{ab}\equiv 1-b/a$ strengthens above $M_{*,30}\sim10^{10.5}\,M_\odot$ in both simulations, but remains systematically weaker in EAGLE. In TNG, outer bar regions are commonly associated with elevated $S_a$ and $S_c$, indicating enhanced boxiness and more prominent box/peanut-shaped bulges, whereas such higher-order signatures are weak or absent in EAGLE. At the highest stellar masses, flattened disks become less prominent, while inner prolate or triaxial structures remain common and massive EAGLE galaxies have more prolate or triaxial outer stellar bodies than their TNG counterparts. These results demonstrate that Gal3D provides a practical framework for quantifying intrinsic radial 3D structure and comparing morphology across cosmological simulations.
Read the full abstract

Earth and Planetary Astrophysics

Chondrule formation in the outer disk from the primary three-dimensional chemical composition of CM chondrules

Chondrules and their associated fine-grained rims record key processes in the early protoplanetary disk, yet the links between chondrule chemistry, morphology, and matrix complementarity remain poorly constrained. We investigate the major, minor, and trace element compositions of 66 chondrules and FGRs from the relatively unaltered CM carbonaceous chondrites Asuka 12236, Paris, and Maribo, together with their 3D morphology, using LA-ICP-MS and X-ray tomography. CM chondrules record systematic metal loss and evaporation of Si-rich mesostasis, driving initially CI-like precursor compositions toward more Mg- and Si-rich bulk compositions along the CI ratio line and toward increasingly Si-poor forsteritic assemblages. GEMS-like materials in pristine CM matrices closely mirror chondrule compositions and likely represent complementary condensates derived from evaporated mesostasis. Dust accreted onto chondrules is predominantly CI-like but contains about 14 wt.% complementary condensate material represented by chondritic amorphous silicates, reconciling Mg/Si complementarity between chondrules and matrix with the preservation of primordial organics and presolar grains. Morphological observations show no significant sectioning bias, consistent with CM chondrules being dominated by agglomerates of ~100 um microspherules. Many display grape-bunch textures produced by welding of smaller chondrules with metal-rich or CI-like rims. This structure may explain the chondrule moderately volatile-element plateau at about 0.3xCI. We propose a "micro-chondrule-first" scenario in which localized heating events produced small molten droplets that subsequently accreted CI-like dust and ice, aggregated, and experienced limited aqueous alteration. These observations place new constraints on chondrule formation in the outer disk and highlight the importance of localized melting and aggregation processes.
Read the full abstract

Solar and Stellar Astrophysics

Elemental Composition Evolution during the 2024 September 30 Solar Eruption: A Comparison of Hot and Cool Plasma Components with Solar Orbiter/SPICE, Hinode/EIS, and Chandrayaan-2/XSM

Solar plasma composition differs between the photosphere and corona over a range of timescales, with preferential enhancement of elements with low first ionization potential (FIP). However, the physical origin of the FIP fractionation remains incompletely understood. Furthermore, during flares, the FIP bias also exhibits rapid changes, associated with fast transport of material with different FIP biases. We present novel observations from Solar Orbiter SPICE and EUI, Hinode/EIS, and Chandrayaan-2 XSM instruments, finding rapid abundance changes in the emitting plasma, on timescales of minutes, during the eruptive M7.6-class solar flare observed on 2024 Sept 30. These instruments have wide temperature coverage and find contrasting abundance-evolution patterns between the hotter and cooler plasma components. 3D reconstruction of the active region and additional observations from the Solar Orbiter STIX X-ray telescope show how the hot and cool plasma components, emitting in different spectral regions and observed with the various instruments, sample the plasma composition evolution in distinct locations within the observed flaring plasma. The bright post-flare loop tops observed by SPICE show coronal FIP bias, while the hot plasma observed with XSM exhibits FIP-bias decreasing from coronal to photospheric during the impulsive phase. We interpret these observations as evidence of the X-ray diagnostics seeing hot coronal reconnection outflows mixing with chromospheric plasma as flare loops sequentially energize and relax, explaining why the FIP bias decreases from coronal to a hybrid; and the cool loop tops seen with SPICE show coronal abundances due to coronal material deposited near the looptops.
Read the full abstract

All authors Momchil Molnar, Joseph Plowman, Amir Caspi, Ritesh Patel, Arpit Shrivastav, Tania Varesano, Don Hassler, Ryan French, Biswajit Mondal, L. P. Chitta

Cosmology and Nongalactic Astrophysics

X-ray signals converted from high-frequency gravitational waves emitted by spinning light primordial black hole dark matter

We investigate the detectability of high-frequency gravitational waves from the superradiance of light primordial black hole dark matter through photons converted in the Galactic magnetic field. We find that the signal is significantly enhanced in the X-ray frequency range around $10^{18}$Hz. For clustered initial conditions, future X-ray observations may detect the converted photons from primordial black holes in the mass range $10^{-15}M_{\odot} \sim 10^{-13}M_{\odot}$. Our results indicate that future X-ray observations could provide a new probe of light primordial black hole dark matter through high-frequency gravitational waves.
Read the full abstract

Instrumentation and Methods for Astrophysics

Tropospheric Ducting Prediction based on GFS Model Data at Inyarrimanha Ilgari Bundara, the CSIRO Murchison Radio-astronomy Observatory

Tropospheric ducting transports terrestrial RFI over hundreds of kilometres into the ARQZWA at Inyarrimanha Ilgari Bundara, the CSIRO Murchison Radio-astronomy Observatory. We present a ducting prediction method based on vertical refractivity profiles derived from GFS data, validated against seven years of continuous spectrum monitoring between 88 MHz and 2.68 GHz. Per-transmitter basic transmission loss to the site is computed with Rec ITU-R P.452-16 above 100 MHz and Rec ITU-R P.1812-6 below it, over Copernicus GLO-90 terrain, driven at each timestep by the GFS-derived refractivity lapse along the path. Skill is scored as a stratified AUC, computed within month and three-hour strata. The median stratified AUC is 0.724 and 24 of 25 confidence intervals exclude chance. A time-percentage calibration derived from four pilot months transfers without retuning to the full record, and skill is flat through three days of forecast lead. A pooled GBM on the same physical features, augmented by a duct-corridor connectivity metric, reaches a mean test AUC of 0.806 over all 25 channels and 0.799 over the 16 on which the analytical predictor is itself skilled. In-situ LTE cell decoding confirms attribution: 86-97% of over-the-horizon cell detections fall within flagged event hours, the range spanning the four monitored LTE channels, and the decoded network identities match the licensing records. Ship AIS receptions at 162 MHz, each a self-located 12 W transmitter at a known sea position, provide VHF ground truth. The FM channels retain only weak skill, far below the ducting channels, and we attribute their events to near-threshold local sources after excluding ducting, sporadic-E and aircraft scatter. The method enables adaptive scheduling of observations away from frequencies affected by forecast ducting, and drives a live forecast service at the observatory. (Abstract abridged for arxiv)
Read the full abstract

Cosmology and Nongalactic Astrophysics

Waterfall-modulated $\alpha$-attractors

Hybrid $\alpha$-attractor models \cite{Kallosh:2022ggf} can have significantly greater values of $n_{s}$ and smaller $r$, while preserving the relation $r\cong 3\alpha (1-n_s)^2$, which is valid for exponential T- and E-models at large values of the inflaton field. Here we study single-field $\alpha$-attractors with features inspired by hybrid models: one can uplift the potential, and one can also have a waterfall regime that leads to a premature termination of inflation near the critical point $\varphi_c$. This allows one to increase the effective number of e-foldings $N_c$ in formulas like $n_s\simeq 1-{2\over N_c}$, $r\simeq {12 \alpha\over N_c^2}$. By changing the waterfall's steepness and location, one can continuously move the predictions along the curves with $r\cong 3\alpha (1-n_s)^2$ as $n_s$ increases and $r$ decreases. We also study the effect of waterfall insertions and uplift on $n_s$ in quintessential $\alpha$-attractors that describe inflation and dynamical dark energy.
Read the full abstract

Solar and Stellar Astrophysics

Photospheric Kelvin–Helmholtz Vortices as Possible Drivers of Coronal Heating: Implications of the DKIST Observations

The Daniel K. Inouye Solar Telescope (DKIST) has resolved Kelvin–Helmholtz (KH) vortices at photospheric magnetic-flux boundaries with a characteristic wavelength of 65 km. I estimate whether these vortices can supply the photospheric driver for cross-scale plasma heating through reconnection across different heights from photosphere to low-corona. Using the simulated MURaM shear, density contrast, and 500 km vertical extent, together with a representative photospheric density, gives a shear-energy density of $1.35\times10^{2}$ J m$^{-3}$ and $2.2\times10^{24}$ erg per characteristic vortex. Magnetic fields $1^\circ$–$7^\circ$ from the exact perpendicular orientation ($ {\bf{B}}\perp{\bf{k}}$) remain KH unstable in an idealized calculation and provide an in-plane component that can be wound or compressed into current layers. The limiting case, in which the center-of-momentum shear reservoir becomes new magnetic free energy, gives $b_{\rm cs}=184$ G, identical to the ideal marginal-stability field and equivalent to a $7.6^\circ$ effective twist. This stores at most 135 J m$^{-3}$ in the layers. Using empirical collisionless reconnection heating fractions of 0.28–0.44, the same twist mapped to weakly collisional heights gives ion heating from $\approx$20 eV at the photosphere to $\approx$1.4 keV in the low corona. For an illustrative, snapshot-based KH-active surface fraction of 0.03, quiet-Sun and coronal-hole losses require 5–8% and 14–21%, respectively, of the shear reservoir to become reconnecting magnetic free energy that reaches such heights. Active regions likely require a separate guide-field twist and helicity reservoir. The required upward transport has not been measured by DKIST, but it is directly testable.
Read the full abstract

Earth and Planetary Astrophysics

Are Hot Jupiters Tidally Disrupted During Stellar Main Sequence?

Once hot Jupiters (HJs) reach their very close orbits, they are expected to experience orbital decay due to tidal interactions with their host star. However, the strength of tidal dissipation is highly uncertain, and it remains an open question whether HJs are tidally disrupted during the stellar main sequence. A previous study found that HJ hosts have a smaller Galactic total velocity dispersion than their field star counterparts, which they interpreted as evidence of tidal disruption. We revisit this study and find that, after using the more reliable vertical velocity dispersion ($\sigma_W$) as the age indicator and accounting for the heterogeneity and anisotropy of their HJ sample, the kinematic age difference between their HJ hosts and matched field stars is significantly reduced. As an independent check, we collect HJs newly discovered by TESS and find that their $\sigma_W$ is statistically similar to that of matched field stars. We also find no statistically significant $\sigma_W$ difference between the field stars and the theoretically vulnerable ultra-hot Jupiters with $P<2$ d. Our results suggest that, after accounting for systematics in the age–velocity dispersion relation, there is no statistically strong evidence from the stellar kinematics that a large fraction of hot Jupiters around Sun-like stars are tidally destroyed during the stellar main sequence.
Read the full abstract

All authors Qingru Hu, Wei Zhu, Yang Huang, Bowen Zhang

Earth and Planetary Astrophysics

Are Near Resonant Multiple-planet Systems from Kepler Young?

Recent studies have claimed that Kepler multi-planet systems hosting near-resonant planet pairs—particularly those near second-order mean-motion resonances (MMRs)—exhibit smaller stellar velocity dispersions than the general population of Kepler planet hosts. Interpreting velocity dispersion as an age indicator, these works concluded that near-resonant systems are systematically younger. We revisit this claim, but we explicitly account for contamination by thick disk stars, which are kinematically hotter and follow a different age-velocity dispersion relation (AVR) than thin disk stars. Using the kinematic criterion to separate thin and thick disk stars, we show that systems classified as having plausible second-order resonant pairs are preferentially hosted by brighter, closer stars and are therefore less contaminated by thick disk stars than the overall sample. After applying a cut to remove probable thick disk contaminants (${\rm TD/D}<0.1$), the vertical velocity dispersion of near-resonant systems becomes statistically indistinguishable from that of the overall Kepler multi-planet sample. We conclude that the apparent kinematic youth of near-resonant systems in Kepler may not be due to a genuine age difference, but rather arises from observational selection effects linked to host star properties and planet detectability. We also comment on the kinematic ages of ultra-short-period planets (USPs).
Read the full abstract

Astrophysics of Galaxies

ALMA observations of pre-JWST z ~ 10 galaxy candidates: A CO(J = 9-8) line from a ULIRG at z = 2.54 and revisit of the photometric redshifts with JWST photometry

We present Atacama Large Millimetre/submillimetre Array (ALMA) observations targeting the [OIII]$88\,\mu$m line for six $z\sim10$ galaxy candidates selected with the Hubble Space Telescope and the Spitzer Space Telescope. We detect a line ($4.5\sigma$) and dust continuum emission ($30\sigma$) in UDS_18697, while detecting neither robust line nor continuum emission in the remaining five objects. The detected line in UDS_18697 is identified as CO($J=9-8$), because follow-up James Webb Space Telescope (JWST) NIRSpec observations have confirmed the redshift as $z=2.54$. UDS_18697 is classified as an ultra luminous infrared galaxy (ULIRG) with far-infrared (FIR) luminosity of $L_\mathrm{FIR}\approx1.1\times10^{12}\,L_\odot$, assuming a dust temperature of $T_\mathrm{d}\approx42\,$K, estimated using a physically-motivated method. We find that UDS_18697 follows the $L_\mathrm{FIR}-L'_\mathrm{CO}$ relation for local and $z>2$ galaxies, albeit being slightly brighter in CO($J=9-8$). Also, based on the follow-up NIRSpec observations and spectral energy distribution fitting using JWST/NIRCam photometry, we found that most of our targets are suggested to be low-$z$ interlopers. Motivated by these redshift misclassifications, we investigate colour–colour selection criteria for high-$z$ galaxies using JWST spectroscopic survey catalogues. We find that elevating a colour threshold tracing the Lyman break is crucial for constructing a robust high-$z$ sample, particularly for wide field surveys such as Euclid Deep Fields and Roman High-Latitude Wide-Area Survey.
Read the full abstract

All authors Suzuka Arai, Yuma Sugahara, Akio K. Inoue, Takuya Hashimoto, Ken Mawatari, Yi W. Ren, Steven L. Finkelstein, John R. Weaver, Rebecca L. Larson, Seiji Fujimoto, Yuichi Harikane, Takahiro Morishita, Yoichi Tamura, Andreas Faisst, Charles Steinhardt, Nima Chartab, Larry D. Bradley, David B. Sanders

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

Scalar induced gravitational waves as probes of dark QCD

We investigate scalar induced gravitational waves (SIGWs) as probes of a dark QCD crossover. Motivated by twin Higgs and asymmetric twin baryon dark matter scenarios, we consider a dark QCD sector with a confinement scale approximately 5.5 times the Standard Model (SM) QCD scale. We construct the effective energy and entropy degrees of freedom for the SM supplemented by dark QCD sectors containing either three light dark quark flavors or all six dark quark flavors. The resulting equation of state parameter and sound speed are then used to solve the first-order scalar perturbations and the second-order SIGWs through the SM and dark QCD crossover epochs. For a monochromatic primordial curvature power spectrum, we first demonstrate that the realistic SM thermal history modifies the SIGW spectrum relative to the idealized radiation-dominated case. We then show that a dark QCD crossover generates an additional frequency-shifted distortion when the enhanced scalar mode reenters the horizon near the dark confinement scale. This distinctive feature can therefore serve as a characteristic signature of the dark QCD sector. Our results demonstrate that SIGWs provide a complementary cosmological probe of hidden confining sectors, with characteristic spectral features shifted to higher frequencies relative to the SM QCD imprint. The analysis developed in this work can also be extended to other well-motivated theories containing different dark confining sectors.
Read the full abstract

Astrophysics of Galaxies

SPURS: Massive Stars, Dense Gas, and Ly$\alpha$ Escape in GN-z11 at $z = 10.6$

We present ultra-deep JWST spectroscopy of GN-z11 ($z=10.6$) obtained through the SPURS Cycle 4 Large Program, providing the deepest rest-UV view yet obtained of a galaxy at $z>10$. GN-z11 was previously found to be nitrogen-enhanced with detectable Ly$\alpha$. The SPURS spectrum reveals P-Cygni stellar wind features and broad He II emission that are jointly reproduced by stellar population models incorporating very massive stars (VMS; $>100\,M_\odot$) at low metallicity and young ages ($\lesssim3$ Myr). We also detect a broad ($\rm FWHM=1670$ km s$^{-1}$) component to N IV] $\lambda1486$, now seen in several nitrogen emitters, potentially arising from dense WN-like winds or LBV-like outbursts associated with a population of VMS in a dense environment, though an AGN-driven wind cannot be excluded. In either scenario, this broad component may trace the gas producing GN-z11's nitrogen enhancement. Rest-UV absorption lines reveal a fast ($\sim500$~km~s$^{-1}$), highly ionized outflow and a negligible neutral gas covering fraction. We resolve the weak Ly$\alpha$ emission (EW=5.6 Å, $f_{\rm esc,Ly\alpha}=2.7$%), finding a broad red wing (44% of flux at $>500$ km s$^{-1}$) that should experience reduced IGM damping wing suppression and help explain Ly$\alpha$ visibility at $z>10$. Fine-structure O I* $\lambda1304$ emission indicates dense neutral gas near a subset of the ionizing sources, which may also scatter Ly$\alpha$ to the large observed velocities. The weak low-ionization absorption favors a picture in which this dense neutral gas is confined to a compact nuclear region. Together, these results are consistent with a rapid burst of star formation building up the dense nuclear regions and surrounding clusters in GN-z11.
Read the full abstract

All authors Zuyi Chen, Daniel P. Stark, Charlotte A. Mason, Peter Senchyna, Mengtao Tang, Keerthi Vasan G. C., Lily Whitler, Adele Plat, Viola Gelli

High Energy Astrophysical Phenomena

An Evolving Leptonic Jet Model for Delayed Radio Flares in Neutrino Blazars

The jets of blazar active galactic nuclei (AGN) are promising sites of hadron acceleration and subsequent neutrino production, owing to their extreme intrinsic power and high radiation density. Potential associations of IceCube neutrinos with blazars displaying delayed radio flares, such as TXS~0506+056 and PKS~1424+240, may support this scenario. However, the mechanisms and location of particle acceleration in the jet remain unclear. The 2017 IceCube event associated with TXS 0506+056 was concurrent with the initial peak of a flare in gamma rays, followed by a longer flare at radio frequencies peaking three years later. State-of-the-art single-zone radiative source models focus on the compact high-energy-emitting region, and thus fail to describe this subsequent radio enhancement. In this work, we model the time-domain multi-frequency evolution of TXS 0506+056 by considering electron emission along the physically extended jet. We couple a numerical particle interaction framework with a dynamic description of the temporal and spatial jet evolution down to the parsec scale. We fit the model to multi-wavelength data, including multi-frequency radio light curves. The results suggest that the parsec-scale jet meets a population of older, cooled electrons near the edge of the observable radio core, establishing a causal relation between the initial high-energy emission and trailing radio activity. This spatially and temporally resolved modeling approach can offer new insights into the physical conditions along the jet of flaring neutrino candidate blazars.
Read the full abstract

Solar and Stellar Astrophysics

Titanium Oxide Absorption as a Proxy to Detect Long term Variation and Activity Cycle in Proxima Centauri

Stellar activity cycles on magnetically active stars can be estimated by molecular absorption bands. We have previously introduced a molecular index which compares absorptional line strength of the TiO567nm with its nearby continuum has previously been introduced. In this work we use this indicator to evaluation long-term activity variations for Proxima Centauri star, using spectroscopic data from HARPS. The results indicate periodicity with an activity period of 2873(+47.4-53.9) days, which is similar to the previous measurements from other indicators.
Read the full abstract

Wednesday 12 August

Cosmology and Nongalactic Astrophysics

CosmoFit: A Graphical User Interface for Bayesian Cosmological Parameter Estimation with User-Defined H(z) Models

We present CosmoFit, a desktop graphical interface for Bayesian parameter inference in background cosmology. The application is designed to simplify cosmological analyses by reducing the need to manually program each step, while maintaining the transparency and reproducibility of the results. Users can define an expansion history H(z), specify fixed and sampled parameters, assign prior ranges and proposal scales, select the available observational likelihoods, and run Markov chain Monte Carlo analyses from a unified graphical interface. CosmoFit uses Cobaya for Bayesian sampling and GetDist for posterior analysis and visualization. The current version supports cosmic-chronometer data and the official Pantheon+, Pantheon+SH0ES, and Union3 supernova likelihoods available through Cobaya. For each run, the application stores the configuration files, metadata, execution logs, chains, and numerical summaries required to inspect and reproduce the analysis. The performance of CosmoFit was evaluated through unit tests, a parameter-recovery test using synthetic data, and an independent comparison with an analysis performed directly in Cobaya. The posterior mean values obtained with both procedures agree to within less than one percent. In addition, the synthetic-data test recovers the fiducial parameters within the expected credible regions. As an illustrative example, we perform a spatially flat wCDM analysis using Pantheon+SH0ES and obtain marginalized constraints and triangle plots directly from the interface. Within its current scope, CosmoFit provides an accessible, transparent, and reproducible tool for research and academic applications involving explicit models of the background cosmological expansion.
Read the full abstract

High Energy Astrophysical Phenomena · Cross-listed from hep-ph

Massive cold hybrid stars in a modified Polyakov-Nambu-Jona-Lasinio model

We propose a modified Polyakov-loop Nambu–Jona-Lasinio (mPNJL) model in which the Polyakov potential is given by an explicit dependence on the quark chemical potential, allowing it to remain finite at zero temperature and thus to describe the confinement-deconfinement transition in cold dense matter. Combining this modified quark sector with hadronic equations of state via a Maxwell construction, we find that, depending on the model parameters, the equation of state can exhibit either two phase transitions, from hadronic matter to confined quark matter and subsequently to deconfined quark matter, or a single transition directly from hadronic to deconfined quark matter or from hadronic to confined quark matter. Stable massive cold hybrid stars with only confined and/or deconfined quark phase are obtained. We systematically examine how the parameters of the modified Polyakov potential and the quark vector interactions control the location of these transitions, and find that repulsive vector interactions are essential to obtain a stable quark core. Hybrid stars with confined and/or a deconfined core can reach maximum masses above $2M_\odot$, provided a sufficiently stiff hadronic equation of state is used at low density. In the core of the maximum-mass configurations, the speed of sound shows variations at finite baryon densities, with c$_s^2(\mu)$ departing from the asymptotic conformal value c$_s^2$ = 1/3 in confined core stars. These variations serve as a diagnostic of the equation-of-state stiffness while remaining fully consistent with causality and thermodynamic stability. This work establishes the qualitative role of each model parameter in shaping hybrid-star structure.
Read the full abstract

All authors Sk Md Adil Imam, Pedro Costa, Mariana Dutra, Odilon Lourenço, Renan Pereira, Constança Providência

High Energy Astrophysical Phenomena

How Neutron Star Radii Encode the Dense-Matter Equation of State and Hadron-Quark Transition

We investigate how future high-precision neutron star (NS) radius measurements encode microscopic information about the dense-matter equation of state (EOS), focusing on a possible first-order hadron–quark phase transition and the resulting mass–radius topology. Within a Bayesian framework using meta-model EOSs with nine microscopic parameters, we analyze mock radius measurements $R_{1.4}=11.9\pm\sigma_R$ km with $\sigma_R=0.9$ and $0.1$ km for canonical NSs. We introduce inverse EOS–radius mappings by studying the conditional mean of each EOS parameter as a function of the model-predicted canonical radius $R_{1.4}$. Their slope measures the sensitivity of an EOS parameter to the radius, while their curvature determines the leading precision dependence of its posterior mean through the Jensen expansion. The symmetry-energy parameters $L$ (slope) and $K_{\rm sym}$ (curvature) are strongly encoded in $R_{1.4}$ and their posterior means shift appreciably with improved radius precision, whereas the higher-order hadronic parameters show stronger topology dependence. Among the transition parameters, the transition density $\rho_t$ is the most strongly encoded in $R_{1.4}$, while the energy-density jump and quark-matter sound speed are more strongly associated with the topology of the full mass–radius sequence. Since the different topologies have strongly overlapping $R_{1.4}$ distributions, even precise radius measurements cannot by themselves identify the topology or uniquely determine the high-density transition properties. These results provide a parameter-dependent hierarchy for assessing the scientific return of future high-precision radius measurements and complementary probes of high-density NS EOS.
Read the full abstract

Astrophysics of Galaxies

Gaia parallax bias via spherical harmonics: A Python tool and discussion of possible causes

Parallaxes in Gaia DR3 are known to suffer from a complex set of sky-correlated and magnitude-dependent offsets or biases at the level of a few tens of $\mu$as. Estimated from a sample of one million distant quasars and AGNs from the CRF catalog, the average offset is negative, but the actual distribution of this important parameter shows significant variations on the sky. We propose a practical method to evaluate the parallax correction as a function of sky position and, optionally, of $G$ magnitude using a spherical harmonic series, and supply a tested Python tool varpi3.py available on Zenodo\footnote{ https://zenodo.org/records/21708614}. We find that only the constant $Y_{00}$ term is significantly dependent on magnitude, while the other 80 harmonic terms are either close to zero or flat with magnitude. The directions of the smallest and largest parallax offsets are $(l,b)\simeq(220\degr,+43\degr)$ and $(l,b)\simeq(45\degr,-45\degr)$, which are close to the orientation of the quasar density dipole reported in recent publications. Motivated by this curious coincidence, we review possible physical effects resulting in a negative bias of measured parallaxes, including an anisotropic universe with a positive curvature and an orbital aberration component. The proposed method of parallax correction is tested using independent asteroseismology data for four different areas on the sphere. Finally, we show that the parallax zero-point propagates into the CRF proper-motion field through the parallax–proper-motion covariance, biasing the vector spherical harmonic determination of the secular-aberration glide, and hence the Galactocentric acceleration, at the microarcsecond-per-year level.
Read the full abstract

Solar and Stellar Astrophysics

The Wave-Regulated Precursor of a Near-Parallel Interplanetary Shock Observed by Parker Solar Probe

Diffusive shock acceleration, at shocks from coronal mass ejections to supernova-remnant blast waves, presupposes a scattering wave field that the accelerated particles themselves maintain. This self-regulation has not been resolved in situ. We report Parker Solar Probe observations of a fast (~2800 km/s), near-parallel interplanetary shock at 0.24 AU on 2023 March 13 and separate its upstream wave field into four families, a classification not made before at a fast shock near the Sun. Right-hand and left-hand circularly polarized families over a common wavenumber band, with a field-aligned linearly polarized family, are cyclotron-resonant with the suprathermal-to-MeV protons streaming from the shock: the beam drives the field that scatters it, and the measured mean free path, half the precursor scale, leaves the beam anisotropic enough to sustain the drive. Outside this loop lies a weak, oblique, linearly polarized component, a few per cent of the wave power, resolved here for the first time at an in situ foreshock. Its in-phase density and field-magnitude fluctuations identify the compressive part as fast magnetosonic and shift the cyclotron-resonance energies of the resonant families by up to 13 % along the precursor. Acceleration at shocks inside 0.3 AU is governed upstream, in a foreshock the shock builds for itself.
Read the full abstract

All authors Immanuel Christopher Jebaraj, Lucas Colomban, Oleksiy Agapitov, Michael Gedalin, Mikhail Malkov, Athanasios Kouloumvakos, Edin Husidic, Shiladittya Mondal, Sunil Yadav, Nicolas Wijsen

Instrumentation and Methods for Astrophysics

Laser-micromachined silicon-platelet feedhorns for large-scale submillimeter and millimeter-wave focal planes

We present the fabrication and characterization of the first silicon-platelet feedhorn arrays produced using laser micromachining. First, we present a demonstration of the technology for the millimeter-wave band of 80~GHz to 170~GHz, i.e. covering the 90/150~GHz bands typical of CMB experiments. Next, we expand the technology to large-scale production on 150~mm wafers and demonstrate operation at submillimeter wavelengths. This feedhorn array is optimized for operation in a band centered at 350~GHz (330~GHz to 370~GHz) and is being deployed as one of the focal plane elements of the CCAT 350~GHz module of Prime-Cam. We present the design and fabrication processes for these feedhorn arrays and compare the optical performance directly to simulation and to feedhorns of identical design but produced using traditional deep reactive-ion etching (DRIE). We conclude with a discussion of future expansions of this technology, including the potential of sidewall control and using thicker (and thus fewer) wafers, which could significantly reduce production costs and labor.
Read the full abstract

All authors Jason E. Austermann, James Beall, Andrew Forsman, Haibo Huang, Johannes Hubmayr, Matthew A. Koc, Jeff van Lanen, Pavel Lapa, Josh Raimist, Sara M. Simon, Jordan Stutz, Anatolios Tambazidis, Joel N. Ullom

Astrophysics of Galaxies

Observing Co-Located Neutral and Ionized Gas-Phase Iron Depletion in the Magellanic Clouds

Depletion is the observed phenomenon where gas-phase elemental abundances are reduced through accretion onto dust grains. We measure neutral gas-phase elemental abundances (S, Fe) in the Magellanic Clouds along 33 sightlines using high-resolution UV spectroscopy (HST/COS and HST/STIS), and compare them to ionized gas-phase abundances (S, Fe) adopted from the literature for six co-located H\,ii regions (with the furthest separation of $\lesssim3'$, 50 pc). Comparing S abundances show that S is minimally depleted in the H\,ii regions and surrounding diffuse ISM. However, we find that the gas-phase Fe abundances in H\,ii regions can be lower than those of the neighboring neutral ISM by 0.3 to 2 dex. This difference is likely an offset in the amount of Fe depleted into dust grains. As accretion of gas-phase Fe is likely not effective at the temperatures of the H\,ii regions, Fe depletion into solid form would have occurred in the dense atomic or molecular clouds prior to star formation. Stronger depletion in the H\,ii regions shows that Fe-bearing grains survive destruction in the first few million years following ionization. Our observations highlight that Fe depletion in H\,ii regions can be a useful tracer of Fe depletion in dense molecular clouds, which are challenging to observe directly via UV absorption.
Read the full abstract

All authors Yun Qi Li, Jessica K. Werk, Caleb R. Choban, Julia Roman-Duval, Kirill Tchernyshyov, J. Xavier Prochaska, Doyeon A. Kim, Arianna S. Long

Solar and Stellar Astrophysics

Large-Scale Dynamos Driven by Shear-Flow-Induced Jets

At every scale they occupy, magnetic fields affect various phenomena, including star formation, cosmic ray transport, charged particle acceleration, space weather, transport in planetary atmospheres, and laboratory plasmas. These fields are often generated and sustained by turbulent flows in a process called the dynamo. In 1955, E. N. Parker parameterized the effects of small-scale turbulence to propose a mean-field dynamo theory. The widely used theory reproduces observed large-scale fields but suffers from difficulty in tuning parameters as they are not justified from first principles: Studies of turbulent flows show tangled magnetic fields, which are folded and fragmented into small-scale structures due to shear-flow straining. Here, considering a shear flow that is unstable and driven, we develop analytic theory and perform three-dimensional (3D), advanced computer simulations of turbulence with up to 4096 x 4096 x 8192 grid points, showing ab initio generation of quasi-periodic, large-scale magnetic fields. The generation occurs via the mean-vorticity effect—an additional mean-field dynamo process postulated in 1990. Crucial to this dynamo is the prior generation of large-scale 3D jets, robustly produced as topologically protected and exact nonlinear solutions of the magnetohydrodynamic equations. The jet-driven dynamo applies to shear-driven laboratory and astrophysical systems. These include binary neutron star mergers, where the reported dynamo likely operates on microsecond timescales to produce in milliseconds some of the strongest magnetic fields in the Universe, providing signals for multimessenger astronomy.
Read the full abstract

All authors B. Tripathi, A. E. Fraser, P. W. Terry, E. G. Zweibel, M. J. Pueschel, R. Fan

High Energy Astrophysical Phenomena

3D Radiative Transfer of Lyman-series Lines with SKIRT

Context. High-resolution X-ray spectroscopy and polarimetry provided by XRISM and IXPE offer new diagnostics of the geometry and kinematics of photo-ionised plasmas around compact objects. Interpreting reprocessed X-ray emission in such systems requires full three-dimensional radiative transfer (3D RT) including photon-ion interactions. Aims. We extend the Monte Carlo (MC) RT code SKIRT by implementing the Lyman-series lines of H-like ions, enabling self-consistent modelling of resonance scattering, radiative recombination, and polarisation of these lines in X-ray photo-ionised plasmas. Methods. We implemented Lyman-series transitions (up to $n=10$) for ions with $Z=$1–30, including fine-structure splitting and linear polarisation in resonance scattering. Two channels for the production of the Lyman-series lines (resonance scattering and radiative recombination) are considered. Results. The implementation reproduces analytical expectations and shows good agreement with Cloudy. The SKIRT simulations naturally capture RT effects such as P Cygni profiles and line-profile distortion in optically thick media, which are inaccessible to the conventional 1D RT codes commonly used in X-rays. In 3D geometries, we find that anisotropic illumination and velocity fields significantly modify the Lyman series line ratios and profiles, all of which are observable with XRISM. Conclusions. The extended version of SKIRT provides a powerful framework for interpreting X-ray line spectra and polarisation from photo-ionised plasmas. It is particularly suited for constraining the geometry and velocity structure in the vicinity of compact objects in the XRISM and IXPE era.
Read the full abstract

All authors N. Sameshima, A. Lauwers, B. Vander Meulen, M. Tsujimoto, L.-Y. Gu, M. Baes, P. Camps

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

Spontaneous Scoto-leptogenesis

We propose a low-scale spontaneous leptogenesis scenario within the dynamical minimal scotogenic model for accommodating neutrino masses and inert scalar dark matter simultaneously. Thus, we dub the mechanism Spontaneous Scoto-leptogenesis. In this setup, a rolling Majoron arising from the global $U(1)_{B-L}$ symmetry breaking induces an effective chemical potential for the $B-L$ charge in the presence of $B-L$ violating interactions that allow for the efficient decays and inverse decays of right handed neutrinos (RHN), so it gives rise to the observed baryon asymmetry of the Universe through the electroweak sphaleron conversion. The mechanism becomes effective in the strong washout regime and successfully lowers the viable mass scale of the lightest RHN to the range of TeV scales, thereby making the thermal scotogenic leptogenesis with two hierarchical RHNs accessible to direct tests. We identify the roles of the $\lambda_5$ coupling for spontaneous leptogenesis and inert scalar dark matter through the efficient erasure of the inert scalar asymmetry. We also explore the regime for Majoron dark matter from the kinetic misalignment, showing that a multicomponent dark sector comprising the inert scalar and the Majoron can be realized in the model. The resulting framework provides a unified origin for low-scale baryogenesis, neutrino masses, and multicomponent dark sector, so it can be tested by complementary experimental probes through direct detection experiments, collider searches for inert scalars, and future detection of Majoron dark matter or dark radiation.
Read the full abstract

Astrophysics of Galaxies

Ripples in the OCEANS: Broad Line Variability of Little Red Dots

Little Red Dots (LRDs) are a unique class of compact, red sources discovered in the JWST extragalactic deep fields. Determining if they are indeed powered by accreting supermassive black holes (SMBHs) is one of the main drivers of the intense study of these objects. Evidence for variability in these objects provides a direct test for the active galactic nucleus (AGN) nature of their central engine. In this study, we present a variability analysis of 6 LRDs observed by the $R \sim 2700$ OCEANS survey and leverage archival $R \sim 1000$ spectroscopic data from the CEERS and RUBIES surveys. We report marginal detections of $\rm H\alpha$ broad-line (BL) variability in the LRDs OCEANS-100424/RUBIES-42232 (27% variability at 2.1$\sigma$ significance) and OCEANS-35829/RUBIES-49140 (GlimmIr/Irony; 50% variability at 1.5$\sigma$ significance). The other 4 LRDs in our sample do not show evidence for BL variability, with a 1$\sigma$ upper limit of $4.8 \% - 30\%$ variability between their epochs of observations. We also find no evidence ($<1\sigma$) for continuum variability in our LRD sample. We compare our results to a sample of SDSS-RM quasars to determine the probability of our broad $\rm H\alpha$ variability detections. We find that the probability of reproducing 2 variable and 4 nonvariable quasars is $4.71\%$, corresponding to $\sim 2 \sigma$ departure from typical quasar variability. The detection of BL $\rm H\alpha$ variability in 2 LRDs provides some evidence for the AGN nature of these objects as opposed to pure scattering models.
Read the full abstract

All authors Madisyn Brooks, Kelcey Davis, Jonathan R. Trump, Raymond C. Simons, Erini Lambrides, Pablo Arrabal Haro, Bren E. Backhaus, Nikko J. Cleri, Steven L. Finkelstein, Mauro Giavalisco, Norman A. Grogin, Michaela Hirschmann, Dale D. Kocevski, Anton M. Koekemoer, Rebecca L. Larson, Ray A. Lucas, Stephen M. Wilkins, Stijn Wuyts, Jorge A. Zavala

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

Singlet-Doublet fermion origin of dark matter, neutrino mass and inverse first-order electroweak phase transition

We study the possibility of an inverse first-order electroweak phase transition (IFOEWPT) and observable gravitational waves (GW) in a radiative neutrino mass model of scotogenic type where singlet-doublet (SD) fermions, the lightest of whom is the dark matter (DM) candidate, generate the necessary seesaw at one-loop level. Considering the possibility of light neutrinos being Dirac for simplicity and additional detection prospects, we extend the standard model (SM) with two generations of $SU(2)_L$ singlet and doublet fermions, one singlet scalar, and three right-handed neutrinos (RHNs). While RHNs provide the right chiral parts of light Dirac neutrinos, the SD fermions and the scalar singlet facilitate the one-loop neutrino mass diagram. The neutral component of the lighter SD fermion, stabilized under a residual $Z_2$ symmetry plays the role of DM while the heavier SD fermions strongly couple to the Higgs leading to an IFOEWPT where the Universe undergoes two different first-order phase transition as it goes from the symmetric to the final broken Higgs phase. We constrain the parameter space from the requirements of generating the correct neutrino mass, DM relic as well as IFOEWPT while incorporating the existing constraints from different experiments. The final allowed parameter space of the model can be probed at collider, direct-detection, GW and cosmic microwave background (CMB) experiments in near future.
Read the full abstract

All authors Debasish Borah, Indrajit Saha, Sujit Kumar Sahoo, Narendra Sahu, Shashwat Sharma

Astrophysics of Galaxies

CLASSY. XV. Kinematics and Spatial Distributions of Outflows in Local Highly Star-Forming Galaxies

Star-forming galaxies drive massive outflows that play an important role in galaxy evolution by regulating feedback and influencing the dynamics of surrounding media. Measuring galactic outflow rates is essential for quantifying feedback efficiency and the amount of mass, momentum, and energy deposited into the circumgalactic medium. In this paper, we examine 17 galactic outflows from the CLASSY survey with radiative transfer modeling of UV absorption lines presented in M. Huberty et al. (2024), to study their spatial distributions and kinematic properties. We study the SiII, SiIII, and SiIV ionization states that trace the cool and warm phases of the outflows and find that SiII traced winds generally behave differently than the warmer SiIII and SiIV traced winds. We derive the mass, momentum, and energy loading factors, which we find scale inversely proportional to stellar mass. We find that our measurements of the mass and momentum loading factors are in agreement with the hydrodynamic FIRE-2 simulations. We model the velocity profiles of the winds, with profiles reaching a maximum velocity of 620 km/s on average, in agreement with hydrodynamic simulations from CGOLS. We also investigate the relationship between outflow properties and the age of the stellar population from SED fitting. We find that outflows associated with young star forming regions are more likely to have a column density dominated by cooler gas and have mass outflow rates which decrease with radius.
Read the full abstract

All authors Mason S. Huberty, Cody A. Carr, Claudia Scarlata, Alaina Henry, Matthew Hayes, Xinfeng Xu, Timothy Heckman, Karla Z. Arellano-Cordova, Danielle A. Berg, R. Michael Jennings, Crystal L. Martin, Kaelee S. Parker, Stephane Charlot, John Chisholm, Simon Gazagnes, Weida Hu, Bethan L. James, Claus Leitherer, Matilde Mingozzi, Evan D. Skillman