Peer Revenue Fully funded

Papers

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

Physics Astrophysics

astro-ph Watch Astrophysics

Tuesday 14 April

Solar and Stellar Astrophysics

KIC 3868420: A high-amplitude $\delta$ Scuti-$\gamma$ Dor hybrid star crossing the Hertzsprung gap

Tao-Zhi Yang and Zhao-Yu Zuo

We report a photometric and asteroseismic analysis of KIC 3868420, a newly identified high-amplitude $\delta$ Scuti-$\gamma$ Doradus hybrid star located in the Hertzsprung gap - a short-lived and rarely observed post-main-sequence phase. Using four years of Kepler long-cadence photometry, we detect 36 significant frequencies, including 11 independent modes spanning both low- and high-frequency regimes. Grid-based modeling with MESA and GYRE, including rotation, shows that five independent frequencies match a combination of radial p-modes and non-radial g-modes, supporting its hybrid nature. The best-fit models yield an evolved post-main-sequence star ($M \sim 2.26 - 2.30 M_\odot$, $R \sim 4.41 - 4.43 R_\odot$, $\tau \sim 5.4 \times 10^{8}$ yr), although degeneracies from rotation and unknown inclination preclude a unique solution. KIC 3868420 thus represents a rare example of a high-amplitude hybrid pulsator in the Hertzsprung gap, demonstrating that high-amplitude pulsation can coexist with mixed p/g-mode behavior in a rapidly evolving star. This finding highlights the value of space-based photometry for identifying such rare objects.
Read the full abstract

Solar and Stellar Astrophysics

Detection and analysis of white-light emission in solar flares through light curve diagnostics

Shuanghong Li, Yongliang Song, Kun Wang et al. (+4 more)

White-light flares (WLFs) are crucial for understanding the energy transport and heating processes in the lower solar atmosphere. Systematic studies are highly necessary. However, most WLFs are very weak and difficult to detect. To address this, we propose a new method of detecting WLFs. Through the observations of SDO/HMI, the light curve of each pixel in the flaring region can be obtained. By subtracting the slowly varying background, we obtained a series of rapidly varying radiative pulses. Pixels for which radiative pulses during flares significantly exceed those occurring before and after the flare were identified as WL emission regions. We applied our method to the detection of the X2.2 flare on September 6, 2017 and validated the method. We found that the WL emission in this flare exhibits two phases, and that different regions show distinct WL emission properties. We also detected the WL emission in all the flares (1 X-class, 2 M-class, and 20 C-class) occurred in active region NOAA 12887. It was found that 15 of the 23 flares are WLFs (1 X-class, 2 M-class, and 12 C-class). The occurrence rate of WLFs in this active region is $\sim65\%$. Surprisingly, the occurrence rate of WLFs in C-class flares even reaches up to $60\%$. It should be noted that most of these C-class WLFs are below C5.0. In addition, a C1.0 WLF was identified; this is the lowest GOES-class event with confirmed WL emission to date. These results demonstrate that WL emission is ubiquitous in most flares, even down to C-class events.
Read the full abstract

All authors Shuanghong Li, Yongliang Song, Kun Wang, Xianyong Bai, Xiao Yang, Nian Liu, Ziyao Hu

Instrumentation and Methods for Astrophysics

FRTSearch: Unified Detection and Parameter Inference of Fast Radio Transients using Instance Segmentation

Bin Zhang, Yabiao Wang, Xiaoyao Xie et al. (+11 more)

The exponential growth of data from modern radio telescopes presents a significant challenge to traditional single-pulse search algorithms, which are computationally intensive and prone to high false-positive rates due to Radio Frequency Interference (RFI). In this work, we introduce FRTSearch, an end-to-end framework unifying the detection and physical characterization of Fast Radio Transients (FRTs). Leveraging the morphological universality of dispersive trajectories in time-frequency dynamic spectra, we reframe FRT detection as a pattern recognition problem governed by the cold plasma dispersion relation. To facilitate this, we constructed CRAFTS-FRT, a pixel-level annotated dataset derived from the Commensal Radio Astronomy FAST Survey (CRAFTS), comprising 2{,}392 instances across diverse source classes. This dataset enables the training of a Mask R-CNN model for precise trajectory segmentation. Coupled with our physics-driven IMPIC algorithm, the framework maps the geometric coordinates of segmented trajectories to directly infer the Dispersion Measure (DM) and Time of Arrival (ToA). Benchmarking on the FAST-FREX dataset shows that FRTSearch achieves a 98.0% recall, competitive with exhaustive search methods, while reducing false positives by over 99.9% compared to PRESTO and delivering a processing speedup of up to $13.9\times$. Furthermore, the framework demonstrates robust cross-facility generalization, detecting all 19 tested FRBs from the ASKAP survey without retraining. By shifting the paradigm from “search-then-identify” to “detect-and-infer,” FRTSearch provides a scalable, high-precision solution for real-time discovery in the era of petabyte-scale radio astronomy.
Read the full abstract

All authors Bin Zhang, Yabiao Wang, Xiaoyao Xie, Shanping You, Xuhong Yu, Qiuhua Li, Hongwei Li, Shaowen Du, Chenchen Miao, Dengke Zhou, Jianhua Fang, Jiafu Wu, Pei Wang, Di Li

Instrumentation and Methods for Astrophysics

In Situ Interferometric Spatial Mapping Of A Microwave Kinetic Inductance Detector Array

Chris Albert, Ritoban Basu Thakur, Farzad Faramarzi et al. (+7 more)

We present a method of spatially mapping microwave kinetic inductance detector (MKID) arrays, in a dark setup. MKIDs are superconducting natively multiplexed resonators which enable kilopixel arrays, such as for the proposed Probe far-Infrared Mission for Astrophysics (PRIMA). In such telescope applications one must map the spatial location of each MKID with their individual resonance frequencies. Traditional LED arrays or beam-mapping methods become increasingly difficult as pixel spacing decreases, e.g., 900 {\mu}m separated MKIDs in the spectrometer module of PRIMA. Our new mapping technique uses a cryogenic interferometer in reflection mode. As on-resonance signals reflect from an MKID, they accrue a phase proportional to the path-length, exactly corresponding to their physical distance on the feedline. Specifically, we use a superconducting transmission line that has nonlinear kinetic inductance. The slow-wave structure of this nonlinear device is designed to have a signal speed of 0.64% the speed of light, enabling a compact system. Current biasing this line allows for varying the wave speed and ensuring that the phase measured is periodic within a nulling interferometric mode. Using this setup, we measure a length ordering that reflects the bimodal MKID distribution of a 44 pixel array of MKIDs designed for PRIMA which contains the same spacing as the final kilopixel array design.
Read the full abstract

All authors Chris Albert, Ritoban Basu Thakur, Farzad Faramarzi, Byeong Ho Eom, Sumit Dahal, Andrew Bear, Reinier Janssen, Henry LeDuc, Thomas Stevenson, Peter Day

Earth and Planetary Astrophysics

JWST/MIRI Hydrocarbon and Water Absorption in the Wind of a Young Disk: Signatures of Pebble Drift and Carbon Grain Sublimation

María José Colmenares, Edwin A. Bergin, Ke Zhang et al. (+27 more)

We present JWST/MIRI-MRS observations of ISO-Oph 37, a highly inclined flat-spectrum ($\lesssim$1 Myr old) source, to investigate the chemical composition and dynamical origin of its inner-disk gas. The spectrum reveals a rich combination of molecular emission and absorption: H$_2$O, CO, and OH are detected in emission, while strong absorption is observed from CO, H$_2$O, CO$_2$, HCN, C$_2$H$_2$, and CH$_4$, with no detectable ice absorption features. LTE slab modeling of the absorption yields excitation temperatures of $T_{\rm ex}\sim400-600$ K and column densities of $\log N/{\rm cm}^{2}\sim16-19$, characteristic of warm gas located within the inner few au. The absorption lines are significantly blueshifted relative to the systemic velocity, with mid-IR lines exhibiting larger shifts than near-IR CO absorption. This velocity structure points to a velocity- and temperature-stratified molecular disk wind. In this framework, the absorption directly samples disk material lifted from the inner disk surface, preserving the chemical imprint of the wind-launching region. Along the line of sight, ISO-Oph 37 is unusually hydrocarbon-rich compared to other known absorption systems (GV Tau N and IRS 46), exhibiting high (C$_2$H$_2$+CH$_4$)/HCN, (C$_2$H$_2$+CH$_4$)/CO and H$_2$O/CO column density ratios, while the CO and HCN columns remain broadly typical. We find that these molecular ratios are best explained by enhancement of both hydrocarbons and water, driven by inward drift and sublimation of icy pebbles and by thermal processing of carbonaceous grains at the soot line. ISO-Oph 37 thus demonstrates that carbon-rich inner-disk chemistry can be established early in disk evolution and that it can be directly probed through molecular absorption in disk winds.
Read the full abstract

All authors María José Colmenares, Edwin A. Bergin, Ke Zhang, Geoffrey A. Blake, Klaus M. Pontoppidan, Alexa R. Anderson, John Carr, Emma Dahl, Joan Najita, Jonathan P. Williams, Colette Salyk, Till Kaeufer, Mayank Narang, Ilaria Pascucci, Benoît Tabone, Lucas Cieza, Miguel Vioque, Adrien Houge, Sebastiaan Krijt, Aditya M. Arabhavi, Giovanni Rosotti, John Carpenter, Feng Long, Paola Pinilla, Jayatee Kanwar, Eshan Raul, Karina Mauco, James Miley, Abygail Waggoner, the JDISCS collaboration

Earth and Planetary Astrophysics · Cross-listed from cond-mat.mtrl-sci

Fe-H melting curve below 3 GPa: Implications for hydrogen in the lunar core

Jun Takeshita, Kei Hirose, Suyu Fu et al. (+2 more)

It has been assumed that hydrogen is negligibly incorporated into core-forming metals below $\sim$3 GPa, and therefore the presence of hydrogen in iron cores of small terrestrial bodies including the moon has not been considered. Here we performed high-pressure melting experiments on the Fe-H system under H$_2$-saturated conditions, combined with synchrotron X-ray diffraction (XRD) measurements. Results demonstrate substantial depression of the Fe-H melting curve compared to that for Fe at 1.0-3.3 GPa, indicating that hydrogen is incorporated into liquid iron even at low pressures less than 1 GPa and the solubility is enhanced with increasing pressure. Based on the density of liquid Fe-H derived from diffuse scattering signal in XRD data, we found that the solubility of hydrogen in liquid iron is about 0.9 wt% at 3.6 GPa and likely enhanced to 1.2 wt% at 5 GPa corresponding to lunar core conditions. The 1.2 wt% H causes 9 % density reduction, which might fully explain the observed density deficit of the lunar core with respect to iron, depending on the density estimate from seismological data.
Read the full abstract

All authors Jun Takeshita, Kei Hirose, Suyu Fu, Fumiya Sakai, Koutaro Hikosaka

Instrumentation and Methods for Astrophysics

21 cm Power Spectrum Analysis of North Celestial Pole Observations with the Tianlai Dish Pathfinder Array

Guangzhi He, Shifan Zuo, Jixia Li et al. (+11 more)

The Tianlai Dish Pathfinder Array (TDPA) is a radio interferometer designed to test techniques for 21 cm intensity mapping in the post-reionization universe as a means of measuring large-scale cosmic structure. Using 9 nights of observations targeting the North Celestial Pole (NCP) field, totaling approximately 107 hours of integration time, we analyze data in the frequency range 700-800 MHz (corresponding to redshift $z \sim 0.9$). We do the data format conversion, radio frequency interference (RFI) flagging, calibration, imaging and point source subtraction, and foreground removal via Singular Value Decomposition (SVD). The spherically averaged power spectrum $\Delta^2(k)$ is obtained. This work successfully establishes and validates a comprehensive data analysis framework for the TDPA. We identify key improvements including sky model refinement, increased integration time, and pipeline optimization that will enable future detection of the 21 cm signal through auto-correlation and cross-correlation with optical galaxy surveys.
Read the full abstract

All authors Guangzhi He, Shifan Zuo, Jixia Li, Yichao Li, Furen Deng, Shijie Sun, Reza Ansari, Olivier Perdereau, Peter Timbie, Albert Stebbins, Ayodeji Ibitoye, Fengquan Wu, Yougang Wang, Xuelei Chen

Cosmology and Nongalactic Astrophysics

Probing Collapsed Dark Matter Halos with Fast Radio Bursts

Yuxuan He, Weiyang Wang, Chen Zhang et al. (+1 more)

Observations of ultra-dense substructures in strong lensing systems challenge the standard cosmological model at small scales. Self-interacting dark matter (SIDM), as an alternative to the cold and collisionless dark matter (CDM) of the standard cosmological model, provides a natural mechanism for forming such structures via gravothermal core collapse. We show that strong gravitational lensing of fast radio bursts (FRBs) provides an effective approach to detecting these substructures and probing dark matter self-interactions. Core-collapsed SIDM halos exhibit steeper central density profiles than CDM halos, enhancing the lensing cross section and producing longer time delays between FRB images. We compute lensing properties of core-collapsed subhalos and host halos, including maximal impact parameters and time-delay distributions. We demonstrate that future all-sky monitors, such as BURSTT, SKA2-Low, and SKA2-Mid, which are expected to detect $10^{5}$–$10^{7}$ FRBs over a decade, can measure time-delay distributions with high statistical significance. Modeling collapsed halos with a cored power-law density profile with inner slope $\gamma=3$ and assuming no excess beyond the singular isothermal sphere lens model, we show that our strategy can probe self-interaction cross section strengths of $\sigma_{\text{SI}}/m \gtrsim \min\{18,\, 40\lambda_{\text{sub}}\}\,\text{cm}^2/\text{g}$, where $\lambda_{\text{sub}}$ parameterizes the collapse time of a subhalo relative to that of the isolated case.
Read the full abstract

All authors Yuxuan He, Weiyang Wang, Chen Zhang, Yi-Ming Zhong

Astrophysics of Galaxies

Spatially Resolved Kinematics of SLACS Lens Galaxies. II: Breaking Degeneracies with Lensing and Dynamical Models

Shawn Knabel, Tommaso Treu, Michele Cappellari et al. (+4 more)

We model the dynamical mass density profiles of 14 strong gravitational lens galaxies from the Sloan Lens ACS (SLACS) sample using spatially resolved kinematics obtained from Keck KCWI integral-field spectroscopy. We use the Jeans Anisotropic Modeling (JAM) method, combining 2D kinematic maps with joint constraints from lens models from Hubble Space Telescope imaging. We use informative priors on the anisotropy and intrinsic shape from local galaxies to help break the residual mass-anisotropy degeneracy (MAD). We find nearly isothermal power-law total mass density slopes ($\rho_{\rm tot}\propto r^{-\gamma}$) for the sample with a mean of $\gamma = 2.04\pm0.02$ with intrinsic scatter of $0.08^{+0.03}_{-0.02}$. We fit explicitly for deviations from the pure power-law form that are fully sensitive to the mass-sheet degeneracy (MSD) and constrain the value of the mass-sheet parameter $\rm \lambda_{int}$ for each individual galaxy to an average precision of 5.8%. The mean value of $\rm \lambda_{int}$ for the sample is $1.01\pm0.03$, with intrinsic scatter of $0.11\pm0.03$. Values of $\rm \lambda_{int}$ for individual objects and the scatter in the sample are consistent to $1\sigma$ uncertainty with those found by the Time-Delay COSMOgraphy collaboration's 2025 milestone analysis, which used a spherical analysis of the same dataset, but azimuthally averaged. We thus conclude that on average power-law mass profiles are a good first-order description of the SLACS sample and do not introduce measureable bias in time-delay cosmography. However, our analysis indicates that more flexible mass models should be able to reproduce the highly detailed kinematic datasets more accurately.
Read the full abstract

All authors Shawn Knabel, Tommaso Treu, Michele Cappellari, Simon Birrer, Xiang-Yu Huang, Anowar J. Shajib, William Sheu

Monday 13 April

Earth and Planetary Astrophysics

A First Principles Approach to the 100,000-year Problem

Liam Wheen

The 100,000-year problem concerns the dominant period of glacial-interglacial cycles over the past 800,000 years and their correlation with Earth's orbital eccentricity, despite eccentricity's weak influence on solar radiation. Two theories compete: the astronomical theory, in which orbital forcing drives the cycles with amplification from Earth system feedbacks, and the geochemical theory, in which internal dynamics dominate with orbital forcing synchronising oscillations. We investigate these theories using conceptual models. Augmentations to the Budyko energy balance model fail to reproduce the 100,000-year period, revealing formulation limitations. Linearised versions of existing non-linear ice volume models perform comparably to their full counterparts, indicating the data does not necessitate non-linear dynamics. We develop two simple linear models: a feedforward model aligned with the astronomical theory and a feedback model aligned with the geochemical theory. The feedforward model reproduces the ice volume record well and offers a novel explanation for the absence of eccentricity's 400,000-year period, arising from oceanic heat storage and tropospheric energy responding with differing phase lags. Conservative estimates show bulk ocean temperature variation can be explained by eccentricity alone, challenging the geochemical theory's core assumption. We also show that widespread use of Q65 may bias models towards geochemical explanations by underrepresenting eccentricity. The feedback model's improvement is concentrated around Marine Isotope Stage 11, suggesting this anomalous interglacial reflects Earth-based events rather than a general requirement for feedback mechanisms. We conclude that 800,000 years of glacial cycles can be largely reproduced by a linear astronomical model, emphasising the importance of parsimony when interpreting palaeoclimate data.
Read the full abstract

Solar and Stellar Astrophysics

Inspecting Cloudy Substellar Atmospheres with JWST MIRI Synthetic Magnitudes from Spitzer Mid-infrared Spectra

Jolie LHeureux, Genaro Suárez, Johanna M. Vos et al. (+4 more)

We examine the positions of substellar objects in mid-infrared color-magnitude and color-color diagrams to distinguish between cloudy and cloud-free atmospheres. Using Spitzer mid-infrared spectra of 113 M5-T9 ultracool dwarfs, we derive synthetic photometry for the JWST MIRI F560W, E'770W, F1000W, and F 1280W filters, which cover key absorption features including the ~9 um silicate signa-ture. We find that diagrams involving F770W and F1000W best separate L-type objects with silicate clouds in their photospheres. L dwarfs with mE77ow - mF1000w < 0.03 mag are seven times more likely to host cloudy atmospheres. Diagrams using F1000W and F1280W are less informative due to the lower signal of the spectra at long wavelengths. Current model predictions struggle to reproduce the positions of cloudy, warm brown dwarfs, likely because atmospheric models underestimate the ~9 um silicate feature. Cloudy Sonora Diamondback models better match the observed trends, although this may reflect improvements capturing indirect effects of clouds on the 6.25 um water absorption feature rather than accurately modeling the silicate feature itself. Our analysis indicates that JWST MIRI photometry can efficiently identify new cloudy extrasolar atmospheres for targeted spectroscopic follow-up, optimizing the use of telescope time.
Read the full abstract

All authors Jolie LHeureux, Genaro Suárez, Johanna M. Vos, Stanimir Metchev, Jacqueline K. Faherty, Sherelyn Alejandro Merchan, Kelle L. Cruz

Solar and Stellar Astrophysics

On the possibility of chemically driven convection in red giants. Implications for the He-core flash and mixing above the Red Giant Branch Bump

M. Miguel Ocampo and Marcelo M. Miller Bertolami

Turbulent mixing remains one of the primary uncertainties in the modeling of stellar interiors. In stellar evolution simulations, regions where mixing occurs are typically identified using instability criteria. A particularly interesting situation arises when nuclear reactions produce inversions in the mean molecular weight within stellar interiors. Under these conditions, the material can become unstable to either thermohaline or a Rayleigh-Taylor instabilities. We demonstrate that the standard criterion adopted in stellar evolution calculations does not accurately distinguish between these two regimes. We derive an alternative criterion and show that chemically driven convection in stellar interiors might be viable under much smaller mean molecular weight inversions than it is normally assumed. We investigate whether inversions in the mean molecular weight can trigger chemically driven convection above the red giant branch bump (RGBB) or during the helium core flash. We find that the inversion at the base of the convective envelope above the RGBB is too weak and short-lived to sustain steady-state convection. In contrast, rapid carbon production at the base of the He-flash-driven convective zone can maintain a steady chemically driven convective region. This process could significantly alter our understanding of the He-core flash and warrants further study.
Read the full abstract

Astrophysics of Galaxies

Identifying Changing-Look AGN Transitions in Light Curve Data with the Zwicky Transient Facility

Margaret E. Verrico, K. Decker French, Vivienne F. Baldassare et al. (+10 more)

Changing-Look AGN (CL-AGN) are AGN which transition between Seyfert types, challenging AGN unification models. Most CL-AGN have been identified via repeat spectroscopy, making it difficult to determine the duration and magnitude of the CL-AGN transition. As such, the physical mechanisms behind this transition are still unknown. We use synthetic photometry in combination with ZTF light curve data to develop a new criterion to identify photometric CL-AGN transitions based on changes in g-band magnitude and g-r color. We find that a CL-AGN criterion of $| \Delta g| > 0.4$ mag and $| \Delta (g-r)| > 0.2$ mag recovers a photometric transition in $9.6^{+4.9}_{-3.4}\%$ of CL-AGN hosts over the six-year ZTF survey, including a candidate repeating changing-look event in SDSS J084957.78+274728.9. Using simulated AGN light curves, we estimate the false positive rate among the simulated Seyferts to be $1.6^{+0.19}_{-0.17}\%$. We find that the rate of similar flares among Type 1 Seyferts is $1.2^{+0.87}_{-0.50}\%$ , and among Type 2 Seyferts is $\leq 0.39\%$ over six years. Photometric CL-AGN transitions last between 21 and 560 days, with a median duration of 360 days, consistent with the thermal or orbital timescales for AGN disks. We do not detect a correlation between black hole mass and transition duration, likely due to the small sample of detected photometric transitions. This method can be applied to the upcoming Legacy Survey of Space and Time to identify CL-AGN candidates and test theories of their origins
Read the full abstract

All authors Margaret E. Verrico, K. Decker French, Vivienne F. Baldassare, Colin J. Burke, Laura Duffy, Nicholas Earl, Megan Harrison, Jason T. Hinkle, Alexander Messick, Samaresh Mondal, Yashasvi Moon, Margaret Shepherd, Zachary Stone

Earth and Planetary Astrophysics

Active moons in our Solar System and beyond – Io, Europa, Enceladus, Triton, and exomoons

Caroline Haslebacher, Emeline Bolmont, Marco Cilibrasi et al. (+19 more)

The outgassing signatures of Io, Europa, Enceladus, Triton, and Io-like exomoons are the focus of this review chapter. The rocky volcanic world of Io is unique in our Solar System, with plumes reaching to hundreds of kilometres in altitude. Io-like exomoons could leave signatures strong enough to be detected with ground-based telescopes. The icy moons Europa and Enceladus, with their subsurface oceans, are currently the best candidates for life. Triton is different in many ways and raises unexplored questions. Our knowledge of these active moons is derived from space- and ground-based observations. To understand their origin, we discuss moon formation in general, before examining evidence and signatures of plumes on these moons. Given the accessibility of subsurface oceanic material through the occurrence of plumes, we expand on possibilities to investigate biosignatures.
Read the full abstract

All authors Caroline Haslebacher, Emeline Bolmont, Marco Cilibrasi, Jonathan Grone, Nico Haslebacher, Ravit Helled, Mathilde Kervazo, Niels F.W. Ligterink, Christophe Lovis, Lucio Mayer, Lorenzo Obersnel, Rafael Ottersberg, Apurva V. Oza, C.H. Lucas Patty, Antoine Pommerol, Ganna Portyankina, Alyssa R. Rhoden, Leander Schlarmann, Yuhito Shibaike, Vishaal Singh, Audrey H. Vorburger, Peter Wurz

High Energy Astrophysical Phenomena

Disk instability model incorporating a variable inner disk radius in SS Cyg and U Gem

Mariko Kimura and Osaki Yoji

Previous theoretical studies indicate that the inner disk in dwarf novae evaporates into a high-temperature, optically thin, and geometrically thick accretion flow during quiescence, with the inner edge moving toward the white dwarf at the onset of an outburst. We incorporate this process into the numerical model developed by Kimura & Osaki (2023) and test the code on two representative dwarf novae, SS Cyg and U Gem. By modeling the inner accretion flow, we calculate the optical, ultraviolet (UV), and X-ray luminosities. Our results show that evaporation suppresses the inside-out outbursts without requiring a radially dependent viscosity parameter in the cold state. The observed time delay between the rise in UV luminosity and the onset of the optical outburst is more than one day, which is successfully reproduced when the inner disk is truncated at several $\times 10^{9}$ cm in the standard evaporation model. However, while the modeled accretion rate at the inner disk edge in U Gem accounts for its quiescent X-ray luminosity, the rate in SS Cyg remains insufficient. This discrepancy in SS Cyg suggests that SS Cyg may require either more efficient evaporation or an additional mass supply into the coronal cavity via gas-stream overflow. By accounting for disk evaporation, our simulations offer a refined version of the disk instability model for dwarf nova outbursts that naturally explains the observed multiwavelength light curves.
Read the full abstract

Solar and Stellar Astrophysics

Astrophysical parameters of LS 437 and the nature of X0726-260

Ignacio Negueruela, Sara R. Berlanas, Lee J. Townsend et al. (+2 more)

Be/X-ray binaries, the most common class of high-mass X-ray binaries, are characterised by OBe companions, but display a rich variety of X-ray behaviours. One of the most atypical systems is X0726-260, which also has the earliest optical counterpart among the whole Milky Way and Magellanic Cloud sample. We intend to improve the characterisation of the optical counterpart, LS 437, and to constrain the physical mechanisms responsible for the anomalous properties of X0726-260. We analyse high-quality, high-resolution optical spectroscopy of LS 437 with standard quantitative methodology to derive stellar parameters. We also make use of archival X-ray monitoring. We derive a moderate projected rotational velocity v sin i $\approx$ 155 km/s and a spectral type O7.5 Ve (Teff = 36 000 K), which makes LS 437 substantially earlier than any other Oe star in an X-ray binary. At this spectral type, the stellar wind likely contributes significantly to mass accretion, and the X-ray light curve is strongly suggestive of an orbitally modulated wind accretor. The source shows marked carbon depletion, whereas nitrogen is only slightly above solar abundance. LS 437 is the earliest Oe star known in the Galaxy, alongside HD 155806. Long-term X-ray lightcurves of X0726-260 strengthen the identification of a persistent 34.5 d periodicity as the orbital period, demonstrating that the X-ray emission is orbitally modulated and no X-ray outbursts have occurred over the past 30 years. Likewise, emission features in the optical spectrum indicate a remarkably stable circumstellar disk, with no sign of major structural changes over the past 40 years. All these characteristics set X0726-260 clearly apart from typical Be/X-ray binaries.
Read the full abstract

All authors Ignacio Negueruela, Sara R. Berlanas, Lee J. Townsend, Javier Lorenzo, Klaus Rübke

Astrophysics of Galaxies

The exponential growth of infinitesimal perturbations in the long-term evolution of simulated galaxies

T. Asano and S. Portegies Zwart

Self-gravitating systems of $N$ particles are chaotic. We study how chaotic the Galaxy is, and what the consequences are. We therefore simulated the dynamical evolution of a galaxy-scale distribution of point masses in order to measure the degree of chaos in such a system. These calculations were performed using the softened gravitational $N$-body tree-code Bonsai, with up to 40 million equal-mass particles. Smaller simulations were performed to establish the scaling of the Lyapunov time $t_L$ with $N$. We established the relations between the degree of chaos, the number of particles, and the softening length in the gravitational force calculation of large-scale $N$-body simulations. The moment in which the bar forms appears to be insensitive to infinitesimal perturbations to the initial realisation. In contrast, the bar strength and its further evolution sensitively depend on these perturbations. Interestingly enough, the maximum in the run-to-run variation in the bar strength is at about the maximum bar strength, and it drops until the bar buckles. The galaxies we simulated are highly chaotic, but the softening in the simulations suppresses chaos. Still, our models show considerable variations in the macroscopic behaviour caused by infinitesimal perturbations to the initial conditions. Real galaxies, however, are expected to be orders of magnitude more chaotic than our simulations, and we are unable to quantify the consequences of this. Smooth galactic potentials for studying individual stellar orbits should be handled with caution on timescales longer than the Lyapunov time. When we extrapolate our results to the number of stars in the Galaxy without planets and other minor bodies, we conclude that Milky Way-size galaxies are chaotic on a timescale of less than a million years.
Read the full abstract

High Energy Astrophysical Phenomena

Mapping the emission and spectral properties of the FRI radio galaxy 3C 449 with LOFAR and the VLA

Luca Ricci, Luisa Ostorero, Raffaella Morganti et al. (+3 more)

The jets and lobes of nearby radio galaxies are ideal laboratories to explore the spectral and dynamical evolution of the radio-emitting plasma that emanates from active galactic nuclei. Here, we present a high-resolution radio continuum study of the low-redshift (z=0.01713), Fanaroff-Riley I (FRI) radio galaxy 3C 449 performed by combining radio data at 145 MHz acquired with the LOw Frequency ARray (LOFAR) and archival Very Large Array (VLA) data at 1365, 1485, 4985, and 8485 MHz. Our LOFAR maps have angular resolutions of 20"x20" and 6.0"x6.0", and show the full extent of the known radio emission at the highest angular resolution to date. Our spectral index maps show the distribution of the spectrum in the 145-8485 MHz frequency range over a region that extends beyond 2.5'. The average 145-8485 MHz spectrum is consistent with a single power law and stays approximately constant over the inner ~50" of both jets. Beyond ~50", on both source sides, the higher-frequency spectrum steepens, indicating no significant downstream sites of particle acceleration. By modelling the spectrum under the assumption of equipartition and of a constant magnetic field across the source, we derive maps of the highest radiative age of the particles. At the outer edges of both the northern tail and southern lobe, the spectral age is ~150 Myr. If the latter age were representative of the dynamical source age, the average expansion speed of both jets during the source lifetime would be supersonic, with Mach numbers of M~4.1 and M~2.8 for the northern and southern jets, respectively. As numerical simulations of FRI jets suggest that the source's current expansion is subsonic, the high Mach numbers might arise either from the source being highly supersonic for a small fraction of its lifetime or from a severe underestimation of the spectral age due to particle acceleration on scales of hundreds of kpc.
Read the full abstract

All authors Luca Ricci, Luisa Ostorero, Raffaella Morganti, Judith H. Croston, Martin J. Hardcastle, Timothy W. Shimwell

Earth and Planetary Astrophysics

Homogeneous accretion of the Earth in the inner Solar System

Paolo A. Sossi and Dan J. Bower

Meteorites are classified as either non-carbonaceous- (NC) or carbonaceous (CC), representing bodies that likely formed in the inner- or outer solar system, respectively. Despite its location in the inner solar system, the Earth is thought to contain either minor- (~6 %) or substantial amounts (~40 %) of outer solar system material. However, because neither interpretation leverages variations among multiple isotopic systems simultaneously, Earth's provenance remains equivocal. Here, we examine variations in 10 nucleosynthetic isotope anomalies among planetary- and meteorite parent bodies to show that the linear extension of an array defined by NC bodies in any two isotopic anomalies always intersects the observed isotopic composition of the bulk silicate Earth to within 1 standard deviation. The Earth therefore formed exclusively from inner solar system material whose composition did not vary over the course of accretion and was, on average, unlike that of any chondrite. Extension of the NC array yields isotopic compositions for Mercury and Venus that are more extreme than for Earth, implying a spatial or temporal gradient during the formation of the terrestrial planets.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Constraints on Coupled Dark Energy in the DESI Era

Adrià Gómez-Valent, Ziyang Zheng and Luca Amendola

We investigate the current viability of a well-known coupled dark energy scenario in which cold dark matter (DM) interacts with a spin-0 dark energy component through a non-trivial field dependence of the DM mass. This ultra-light scalar mediates a fifth force between DM particles, which can leave signatures on cosmological scales. We use state-of-the-art data on the cosmic microwave background from Planck's CamSpec likelihood, baryon acoustic oscillations from the second DESI data release as well as the supernovae of Type Ia (SNIa) from Pantheon+ and DES-Dovekie. We perform the analysis considering both a flat potential and a Peebles-Ratra (PR) potential for the scalar field in order to assess the impact of the potential slope on the fitting performance of the model. While for a constant potential the scalar field dynamics is insensitive to the sign of the coupling parameter $\beta$, the PR potential breaks the existing symmetry in the solutions at late times and could induce a difference at the phenomenological level between positive and negative values. We study for the first time if it is actually the case, finding no important asymmetry in the fitting results. In the light of the aforesaid datasets, we find in all cases a peak at $|\beta|\sim 0.03$ - less pronounced than reported in some recent works -, excluding the no-coupling scenario at $\sim 95\%$ CL at most. The model is able to explain an effective crossing of the phantom divide, with the equation-of-state parameter lying within the $2\sigma$ bands of model-agnostic reconstructions. Our results are very robust under changes in the SNIa sample used in the analysis and is not significantly altered when we replace a constant potential with the PR one, although the latter is crucial to produce the aforesaid crossing. In passing, we also provide constraints obtained with the PR potential in the uncoupled case.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Beyond the Standard Model of Cosmology: Testing new paradigms with a Multiprobe Exploration of the Dark Universe

Juan Garcia-Bellido

Cosmology is living through fascinating times, where new observations from ground and space telescopes are questioning the established paradigm, the so-called Lambda Cold Dark Matter model. The particle nature of Dark Matter is severely constrained by underground experiments, while recent observations by galaxy surveys indicate that the cosmological constant (Lambda) may not be constant after all. Furthermore, observations at high redshift of fully-formed galaxies with massive black holes at their centers by the James Webb Space Telescope, as well as black holes with unexpected properties observed by gravitational wave detectors LIGO-Virgo, are driving an in-depth revision of our assumptions in models of structure formation and the evolution of the universe. I propose to explore two new paradigms to account for Dark Matter and Dark Energy, based on known physics without the need for new particles nor new degrees of freedom. I will extend the primordial spectrum of fluctuations to small scales with new statistical properties to provide a viable Primordial Black Hole scenario for Dark Matter, and will include non-equilibrium thermodynamics in the expanding universe, in the form of General Relativistic Entropic Acceleration, to explain Dark Energy. My proposal could provide a unified explanation for a plethora of interrelated multi-epoch, multi-scale and multi-probe observations from present and future Gravitational Wave detectors, Large Scale Structure observatories and Cosmic Microwave Background experiments. It emphasizes the need to develop new theoretical ideas hand-in-hand with observations to acquire a deeper understanding of our universe. If these ideas are correct, they will open a new window into the early universe and a new fundamental understanding of gravity in the late universe.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Effective field theory of a single scalar pion field for large scale structure in the Universe

Lara Celik, Bart Horn, Bhavya Mishra et al. (+1 more)

We discuss the effective field theory of large scale structure in terms of a single scalar degree of freedom, corresponding to the velocity potential of the matter fluid in a $\Lambda$CDM universe. This cosmic “pion” field is nonlinearly related to the overdensity and the gravitational potential, and corresponds to the Goldstone boson of spontaneously broken spacetime symmetry, allowing us to organize perturbation theory in a systematic way while keeping the symmetries manifest. We develop the effective field theory of the pion field to next-to-leading order, and we use it to calculate the corrections to the power spectrum and to check that these are consistent with the consistency relations of spontaneously broken spacetime symmetry. We compare our results against computer simulations for the evolution of large scale structure in the pion field picture, and we make use of N-body simulations to measure EFT coefficients and analyze the growth of additional degrees of freedom in the deep nonlinear regime. We conclude with a discussion of how the pion field picture may help suggest new variables for analyzing simulations and experimental surveys of large scale structure.
Read the full abstract

All authors Lara Celik, Bart Horn, Bhavya Mishra, David Muqattash

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

MMS Insights into CME Driven Sub-Alfvénic Solar Wind at 1 AU

Harsha Gurram, Li-Jen Chen, Matthew R. Argall et al. (+4 more)

We report the properties of electron distributions and turbulence during a Coronal Mass Ejection (CME) in April 2023 observed by Magnetospheric Multiscale (MMS). The CME exhibits a clear sheath and magnetic cloud (MC), and within the MC, the solar wind becomes sub-Alfvénic for two hours. We investigate plasma and turbulence properties of the sub-Alfvénic CME wind and compare them with those in the super-Alfvénic solar wind in the MC and CME sheath. Electrons within the sub-Alfvénic MC show significantly higher temperatures than those in the CME sheath and the super-Alfvénic MC, with their one-dimensional distributions revealing super-thermal tail and a depletion in electron populations between 15-50 eV. Within the CME sheath, isolated regions of electron heating are observed, where parallel energy flux is enhanced up to ~1 keV. Magnetic field fluctuations within the sub-Alfvénic MC interval exhibit negligible cross helicity and steeper-than-Kolmogorov scaling in the inertial range, with no clear spectral break. These fluctuations also show reduced intermittency at ion and sub-ion scales, emerging intermittency at electron scales, and weak magnetic compressibility. Together, these observations point to the presence of weak magnetohydrodynamic (MHD) turbulence within the sub-Alfvénic MC, resembling conditions commonly observed in planetary magnetospheres such as Jupiter's.
Read the full abstract

All authors Harsha Gurram, Li-Jen Chen, Matthew R. Argall, Subash Adhikari, Lynn B. Wilson, Jason R. Shuster, Victoria D. Wilder

Solar and Stellar Astrophysics

OB runaway stars originating in the Vel OB1 association

N. Azatyan, L. Kaper, A. Samsonyan et al. (+4 more)

OB runaway stars are massive stars moving through interstellar space at high velocities (up to 200 km/s), produced by dynamical ejections in young massive clusters or supernova explosions in massive binaries. They can travel several hundred parsec before exploding as supernovae, affecting the dynamical and chemical evolution of the Galaxy. The Vel OB1 association, one of the largest OB associations, hosts about 20 O-type and more than 50 B-type stars. We aimed to identify OB runaways in this region, quantify their number, identify their parent clusters, and understand their production channels and impact on the surrounding medium. Using Gaia DR3 coordinates, parallaxes, and proper motions, we identified OB runaways by measuring their peculiar velocity. We inspected infrared WISE images to identify wind bow shocks and reconstructed runaway trajectories to locate parent clusters and estimate travel times. We identified six young stellar clusters hosting most of the massive-star population in Vel OB1 (distance 1.6-2.1 kpc; age 1-10 Myr) and derived a threshold velocity of 15 km/s to classify runaways. We identified 25 OB runaways (including HMXB VelaX-1) and one F-type runaway. We detected 16 arc-like features, six associated with runaways selected by peculiar velocity, and ten bow shocks aligned with runaway proper motions. Parent clusters are identified for seven runaways, most likely ejected dynamically. The runaway fraction is about 30%. Wind bow shocks from OB runaways reveal valuable information on local ISM conditions.
Read the full abstract

All authors N. Azatyan, L. Kaper, A. Samsonyan, M. Stoop, D. Andreasyan, J. van den Eijnden, E. Nikoghosyan

High Energy Astrophysical Phenomena

Relativistic Effects on Circumbinary Orbit Stability

Gonzalo C. de Elía, Macarena Zanardi and Rebecca G. Martin

With n-body simulations and analytic approximations we study the dynamics and stability of low eccentricity misaligned test particles around binary systems with varying mass fraction and eccentricity. General relativity (GR) plays a primary role in determining the motion of an outer particle since it drives apsidal precession of the binary orbit. The effects of GR can drive particle instability close to the binary orbit, depending upon the binary parameters and the initial inclination of the particle. For the binary parameters we consider, we find instability up to a semimajor axis of about 8 ab, where ab is the binary semimajor axis. In particular, we identify and analyse three different regions of instability that are driven by GR in the phase plane of the initial semimajor axis and the initial inclination of the particle. The results have implications for circumbinary orbits and circumbinary disks on all scales, but are particularly important around supermassive black hole binaries where the effects of GR can be strong.
Read the full abstract

Astrophysics of Galaxies

Reconstructing chemical enrichment pathways in disc galaxies: A phylogenetic approach

Brian Tapia-Contreras, Patricia B. Tissera, Emanuel Sillero et al. (+11 more)

Phylogenetic methods, traditionally used in biology to trace the evolutionary relationships among species, are emerging as a powerful framework to reconstruct evolutionary processes in galaxies from chemical information. We apply galactic phylogenetics to study the chemical evolution of stellar populations in distinct regions of a simulated disc galaxy, assessing its capability to unveil assembly histories. We used a high-resolution simulation that follows the chemical enrichment of an isolated disc galaxy, by different stellar progenitors. We track gas particles as they turn into stars and inherit their parent gas chemical composition. Target particles are selected to store the chemical history of each chemical element considered in the simulation. Two regions were analysed: an inner ring, influenced by early bar-driven inflows, and an outer ring, shaped by spiral arms. We built phylogenetic trees for stellar populations in each region and quantified their structure using the Corrected Colless index, a standard metric of tree balance used in biology. The inner ring tree reveals a compact clade of old stars enriched by rapid SNII feedback, followed by a hierarchical sequence with increasing SNIa and AGB contributions. In contrast, the outer ring exhibits more symmetric, caterpillar-like trees with smoother abundance gradients, consistent with more prolonged star formation and efficient local mixing. Chemical enrichment rates corroborate these trends, showing fast early enrichment in the inner ring and gradual, spatially extended enrichment in the outer disc. The structural indices differ significantly between the two regions and converge robustly even for modest stellar samples (NSSP = 100). Galactic phylogenetics provides a novel and complementary tool to decode the fossil record of galaxies.
Read the full abstract

All authors Brian Tapia-Contreras, Patricia B. Tissera, Emanuel Sillero, Paula Jofré, Keaghan Yaxley, Xia Hua, Robert M. Yates, Álvaro Márquez S., Theosamuele Signor, Payel Das, Álvaro Rojas-Arriagada, Claudia Aguilera-Gómez, Francisco Jara-Ferreira, Robert A. Foley

Solar and Stellar Astrophysics

No country for old stars -Spectroscopic confirmation of the first intermediate-age RR Lyrae in the open cluster Trumpler 5

Valentina D'Orazi, Cecilia Mateu, Giuliano Iorio et al. (+13 more)

RR Lyrae stars are widely considered tracers of ancient (greater than 10 Gyr), metal-poor stellar populations. However, recent kinematic and photometric studies suggest the existence of a metal-rich RRL sub-population associated with the thin disc and intermediate ages (approximately 2-5 Gyr), challenging canonical evolutionary models. We aim to provide the first spectroscopic confirmation of a member of this elusive population. Specifically, we target a metal-rich RRL candidate recently identified photometrically as a member of the intermediate-age open cluster Trumpler 5. We obtained high-resolution spectroscopy using PEPSI at the LBT and GHOST at Gemini-South Telescope. We measured radial velocities from multiple epochs to constrain cluster membership and derived detailed chemical abundances (Mg, Ca, Sc, Ti, Mn, Fe, Cu, Zn, Y, and Ba) to compare the RRL's composition with that of red clump stars in the cluster. The RRL's systemic velocity Vgamma = 50.57 +0.78/-0.36 km/s is in excellent agreement with the cluster mean V = 50.76 +/- 0.49 km/s. Combining RVs, proper motions, and parallax, the probability of the star being a background interloper is negligible (approximately 0.002%, better than 4 sigma). We derived a metallicity of [Fe/H] = -0.40 +/- 0.05, matching the cluster value. While most abundance ratios (Mg, Ti, Mn, Cu, and Zn) align with cluster members, the RRL exhibits significant depletion in Ca, Sc, Y, and Ba. Notably, [Sc/Fe] is under-abundant by approximately 0.6 dex relative to the cluster stars, following trends seen in field metal-rich RRLs. We provide strong constraints on the membership status between an RRL variable and an intermediate-age open cluster [...]
Read the full abstract

All authors Valentina D'Orazi, Cecilia Mateu, Giuliano Iorio, Alexey Bobrick, Zdenek Prudil, Ricardo Salinas, Angela Bragaglia, Lyudmila Mashonkina, Raffaele Gratton, Ilya Ilyn, Natalia Alvarez Baena, Vittorio F. Braga, Antonino Nunnari, Venu Kalari, Felice Cusano, Silvia Tosi

Astrophysics of Galaxies

The role of small-scale environments in the quenching of massive galaxies at $1<z<5$

Takumi Kakimoto, Masayuki Tanaka, Makoto Ando et al. (+10 more)

Massive quiescent galaxies (QGs) at high redshifts are likely progenitors of massive elliptical galaxies in the local Universe. Recent discoveries of high-redshift QGs in overdensities (galaxy groups and proto-clusters) have highlighted the relationship between massive galaxy quenching and their surrounding environment. We spectroscopically confirm a galaxy group associated with a massive QG at $z_\mathrm{spec}=4.53$ from the Lyman break feature using Subaru/FOCAS. This group consists of at least three star-forming galaxies within 150 pkpc of the QG, suggesting a physical association necessary for galaxy quenching. To investigate the role of the surrounding environment, we also perform a statistical analysis to characterize the typical environment of QGs at high redshifts. By selecting QGs using an SFR-based selection in the COSMOS field, we find that the quiescent fraction is higher in overdensities than in the field at $1<z<5$. In particular, the elevated quiescent fraction within small-scale overdensities $(\sim 200\,\mathrm{pkpc})$ demonstrates that environmental quenching, possibly driven by galaxy mergers and interactions, plays a major role in the evolution of massive QGs at high redshifts. By extending the low-redshift mass and environmental quenching model up to $z \sim 5$, we confirm similar trends to those inferred from the SDSS red galaxy fraction at $z<1$. While environmental dependence weakens at $z \gtrsim 3$ for the overall population, a clear correlation persists specifically for more massive galaxies. This mass-dependent environmental effect suggests a downsizing trend, highlighting that the accelerated evolution of massive galaxies in high-density regions is already in place in the early Universe.
Read the full abstract

All authors Takumi Kakimoto, Masayuki Tanaka, Makoto Ando, Kei Ito, Ken Mawatari, Masato Onodera, Rhythm Shimakawa, Francesco Valentino, Po-Feng Wu, Kiyoto Yabe, Shuowen Jin, Mariko Kubo, Sune Toft

Earth and Planetary Astrophysics

How leaky? A large parameter study of leaky dust traps to quantify the transport of pebbles and ice in protoplanetary discs

Adrien Houge, Anders Johansen, Andrea Banzatti et al. (+1 more)

In protoplanetary discs, the presence of dust traps can significantly alter the transport of solids from the outer to the inner regions, and hence they are often invoked as an explanation for the chemical diversity of inner discs observed with JWST (e.g., varying oxygen abundances and C/O ratios). As a detailed treatment of dust transport around dust traps is computationally expensive, earlier works investigating the impact of outer traps on the inner disc composition have often used simplified dust models representing the size distribution with a single effective size and drift speed. In this paper, we revisit the impact of outer traps on dust transport using the state-of-the-art one-dimensional dust evolution code DustPy, which simulates the transport and evolution of dust particles including detailed coagulation and fragmentation. We quantify and map the leakiness of dust traps across a broad parameter space, performing over 300 simulations while varying the disc viscosity, turbulence strength, planet mass and location, and dust fragmentation velocity. We find that dust traps are leakier than previously thought, on a broader parameter space, such that most outer traps (r > 5 au) will result in a long-lived O-rich inner disc with gas-phase C/O < 1. In similar conditions (e.g., carved by the same planet mass), we find inner traps are much leakier than outer traps, though their relative efficiency in reducing the pebble flux is time-dependent. Highly blocking traps altering the inner disc composition dramatically (leading, e.g., to C/O > 1) are possible to set up but necessitate low viscosity and weak turbulence, along with efficient planetesimal formation by the streaming instability. In that case, we find that is the formation of planetesimals, rather than the dust traps themselves, that is capable of significantly altering the inner disc composition.
Read the full abstract

All authors Adrien Houge, Anders Johansen, Andrea Banzatti, Sierra Grant

Earth and Planetary Astrophysics

Sub-Neptunes as Soot Factories: Deep Atmosphere Hydrocarbon Formation and Quenching as the Origin of Sub-Neptune Aerosol Trends

Jeehyun Yang, Eliza M.-R. Kempton and Arjun B. Savel

Recent population-level studies of sub-Neptune atmospheres have identified a tentative parabolic trend in transmission spectrum amplitude for planets with Teq ~ 500-800 K. While the trend has been commonly attributed to hydrocarbon aerosols, we lack a first-principles explanation of its underlying chemical mechanism. Previous work has focused on the role of methane photolysis and subsequent polymerization, but with limited reaction networks that truncated at C2-species and couldn't reproduce the observed parabolic trend. In this work, enabled by a computer-automated, rate-based chemical network generator, we construct the most comprehensive carbon reaction network for exoplanet atmospheres to date. We explicitly model the formation of polycyclic aromatic hydrocarbons (PAHs), which are well established as soot precursors in combustion chemistry. We calculate the chemical timescales of hydrocarbon species through an eigenvalue timescale method and model their quenched abundances across a range of C/O, metallicities, and Teq. In this framework, the deep atmosphere acts as a "soot factory" analogous to a combustion engine, transporting the ingredients for hydrocarbon aerosol formation to the JWST-observable region of the atmosphere, where it may be further augmented by photochemistry. We find that the predicted abundances of PAHs peak near 600 K, and fall off toward higher and lower Teq, consistent with the observed muted-spectra regime suggested in observational studies by HST and JWST. We also show that PAH abundances are expected to vary with C/O and metallicity, thus providing a natural explanation for observed diversity among planets with similar Teq.
Read the full abstract