Abstract visibility
Abstract text size

Papers for Wednesday, Jun 18 2025

Papers with local authors

We present updated non-adiabatic and inhomogeneous evolution models for Uranus and Neptune, employing an interior composition of methane, ammonia, water, and rocks. Following formation trends of the gas giants, Uranus and Neptune formation models are applied, where both planets begin with layers stable to convection. Both planets are subject to convective mixing throughout their evolution. Uranus undergoes modest convective mixing, preserving much of its primordial internal heat. In contrast, Neptune's interior undergoes extensive mixing, homogenization, and adiabatic cooling of the outer 40\% of its envelope. The subsequent release of internal energy in Neptune, driven by the convective instability of its primordial outer compositional gradient, accounts for its higher luminosity relative to Uranus. Thus, the observed luminosity differences between Uranus and Neptune could be primarily dictated by the convective stability of their outer envelopes. The extensive convective mixing in Neptune leads to a higher metallicity in its outer region compared to Uranus, a feature seen in atmospheric measurements and shown in past interior models of Neptune. Due to Neptune's more pronounced cooling, our models predict favorable conditions for hydrogen-water immiscibility in its envelope.

M. Janssen, C.-k. Chan, J. Davelaar, I. Natarajan, H. Olivares, B. Ripperda, J. Röder, M. Rynge, M. Wielgus
0 votes
Paper 11 — arXiv:2506.13873
0 votes
Paper 11 — arXiv:2506.13873

(abridged) In a series of publications, we describe a comprehensive comparison of Event Horizon Telescope (EHT) data with theoretical models of Sgr A* and M87*. Here, we report on improvements made to our observational data reduction pipeline and present the generation of observables derived from the EHT models. We make use of ray-traced GRMHD simulations that are based on different black hole spacetime metrics and accretion physics parameters. These broad classes of models provide a good representation of the primary targets observed by the EHT. To generate realistic synthetic data from our models, we took the signal path as well as the calibration process, and thereby the aforementioned improvements, into account. We could thus produce synthetic visibilities akin to calibrated EHT data and identify salient features for the discrimination of model parameters. We have produced a library consisting of an unparalleled 962,000 synthetic Sgr A* and M87* datasets. In terms of baseline coverage and noise properties, the library encompasses 2017 EHT measurements as well as future observations with an extended telescope array. We differentiate between robust visibility data products related to model features and data products that are strongly affected by data corruption effects. Parameter inference is mostly limited by intrinsic model variability, which highlights the importance of long-term monitoring observations with the EHT. In later papers in this series, we will show how a Bayesian neural network trained on our synthetic data is capable of dealing with the model variability and extracting physical parameters from EHT observations. With our calibration improvements, our newly reduced EHT datasets have a considerably better quality compared to previously analyzed data.

(abridged) In this second paper in our publication series, we present the open-source Zingularity framework for parameter inference with deep Bayesian artificial neural networks. We carried out out supervised learning with synthetic millimeter very long baseline interferometry observations of the EHT. Our ground-truth models are based on GRMHD simulations of Sgr A* and M87* on horizon scales. We investigated how well Zingularity neural networks are able to infer key model parameters from EHT observations, such as the black hole spin and the magnetic state of the accretion disk, when uncertainties in the data are accurately taken into account. Zingularity makes use of the TensorFlow Probability library and is able to handle large amounts of data with a combination of the efficient TFRecord data format plus the Horovod framework. Our approach is the first analysis of EHT data with Bayesian neural networks, where an unprecedented training data size, under consideration of a closely modeled EHT signal path, and the full information content of the observational data are used. Zingularity infers parameters based on salient features in the data and is containerized. Through parameter surveys and dedicated validation tests, we identified neural network architectures, that are robust against internal stochastic processes and unaffected by noise in the observational and model data. We give examples of how different data properties affect the network training. We show how the Bayesian nature of our networks gives trustworthy uncertainties and uncovers failure modes for uncharacterizable data. It is easy to achieve low validation errors during training on synthetic data with neural networks, particularly when the forward modeling is too simplified. Through careful studies, we demonstrate that our trained networks can generalize well so that reliable results can be obtained from observational data.

(abridged) In the first two papers of this publication series, we present a comprehensive library of synthetic EHT observations and used this library to train and validate Bayesian neural networks for the parameter inference of accreting supermassive black hole systems. The considered models are ray-traced GRMHD simulations of Sgr A* and M87*. In this work, we infer the best-fitting accretion and black hole parameters from 2017 EHT data and predict improvements that will come with future upgrades of the array. Compared to previous EHT analyses, we considered a substantially larger synthetic data library and the most complete set of information from the observational data. We made use of the Bayesian nature of the trained neural networks and apply bootstrapping of known systematics in the observational data to obtain parameter posteriors. Within a wide GRMHD parameter space, we find M87* to be best described by a spin between 0.5 and 0.94 with a retrograde MAD accretion flow and strong synchrotron emission from the jet. Sgr A* has a high spin of $\sim$ 0.8 $-$ 0.9 and a prograde accretion flow beyond the standard MAD/SANE models with a comparatively weak jet emission, seen at a $\sim$ 20$^\circ$ $-$ 40$^\circ$ inclination and $\sim$ 106$^\circ$ $-$ 137$^\circ$ position angle. While previous EHT analyses could rule out specific regions in the model parameter space considered here, we are able to obtain narrow parameter posteriors with our Zingularity framework without being impacted by the unknown foreground Faraday screens and data calibration biases. We further demonstrate that the AMT extension to the EHT will reduce parameter inference errors by a factor of three for non-Kerr models, enabling more robust tests of general relativity. It will be instructive to produce new GRMHD models with the inferred interpolated parameters to study their accretion rate plus jet power.

D. Eckert, F. Gastaldello, L. Lovisari, S. McGee, T. Pasini, M. Brienza, K. Kolokythas, E. O'Sullivan, A. Simionescu, M. Sun, M. Ayromlou, M. A. Bourne, Y. Chen, W. Cui, S. Ettori, A. Finoguenov, G. Gozaliasl, R. Kale, F. Mernier, B. D. Oppenheimer, G. Schellenberger, R. Seppi, E. Tempel
0 votes
Paper 18 — arXiv:2506.13907
0 votes
Paper 18 — arXiv:2506.13907

Supermassive black hole feedback is the currently favoured mechanism to regulate the star formation rate of galaxies and prevent the formation of ultra-massive galaxies ($M_\star>10^{12}M_\odot$). However, the mechanism through which the outflowing energy is transferred to the surrounding medium strongly varies from one galaxy evolution model to another, such that a unified model for AGN feedback does not currently exist. The hot atmospheres of galaxy groups are highly sensitive laboratories of the feedback process, as the injected black hole energy is comparable to the binding energy of halo gas particles. Here we report multi-wavelength observations of the fossil galaxy group SDSSTG 4436. The hot atmosphere of this system exhibits a highly relaxed morphology centred on the giant elliptical galaxy NGC~3298. The X-ray emission from the system features a compact core ($<$10 kpc) and a steep increase in the entropy and cooling time of the gas, with the cooling time reaching the age of the Universe $\sim15$ kpc from the centre of the galaxy. The observed entropy profile implies a total injected energy of $\sim1.5\times10^{61}$ ergs, which given the high level of relaxation could not have been injected by a recent merging event. Star formation in the central galaxy NGC~3298 is strongly quenched and its stellar population is very old ($\sim$10.6 Gyr). The currently detected radio jets have low power and are confined within the central compact core. All the available evidence implies that this system was affected by giant AGN outbursts which excessively heated the neighbouring gas and prevented the formation of a self-regulated feedback cycle. Our findings imply that AGN outbursts can be energetic enough to unbind gas particles and lead to the disruption of cool cores.

Katherine A. Suess, Aliza G. Beverage, Mariska Kriek, Justin S. Spilker, Rachel Bezanson, Vincenzo R. D'Onofrio, Jenny E. Greene, Jamie Lin, Yuanze Luo, Desika Narayanan, Imad Pasha, Sedona H. Price, David J. Setton, Margaret E. Verrico, Yunchong Zhang
0 votes
Paper 48 — arXiv:2506.14361
0 votes
Paper 48 — arXiv:2506.14361

Recent observations at low redshift have revealed that some post-starburst galaxies retain significant molecular gas reservoirs despite low ongoing star formation rates, challenging theoretical predictions for galaxy quenching. To test whether this finding holds during the peak epoch of quenching, here we present ALMA CO(2-1) observations of five spectroscopically confirmed post-starburst galaxies at z ~ 1.4 from the HeavyMetal survey. While four galaxies are undetected in CO emission, we detect M_H2 ~ 10^9.7 Msun of molecular gas in one system. The detected system is a close pair of massive (M* = 10^(11.1-11.2) Msun) post-starburst galaxies with no clear tidal features, likely caught in the early stages of a major merger. These results suggest that mergers may be a key factor in retaining molecular gas while simultaneously suppressing star formation in quenched galaxies at high redshift, possibly by driving increased turbulence that decreases star formation efficiency. Unlike previous studies at z < 1, we find no correlation between molecular gas mass and time since quenching. This may be explained by the fact that -- despite having similar UVJ colors -- all galaxies in our sample have post-burst ages older than typical gas-rich quenched systems at low redshift. Our results highlight the importance of major mergers in shaping the cold gas content of quiescent galaxies during the peak epoch of quenching.

Quasi-periodic pulsations (QPPs) at sub-second periods are frequently detected in the time series of X-rays during stellar flares. However, such rapid pulsations are rarely reported in the hard X-ray (HXR) emission of the small solar flare. We explored the QPP patterns with fast-time variations in HXR and radio emissions produced in a small solar flare on 2025 January 19. By applying the Fast Fourier Transform, the fast-variation pulsations at a quasi-period of about 1 s are identified in the HXR channel of 20-80 keV, which were simultaneously measured by the Hard X-ray Imager and the Konus-Wind. The rapid pulsations with a same quasi-period were also detected in the radio emission at a lower frequency range of about 40-100 MHz. The restructured HXR images show that the QPP patterns mainly locate in footpoint areas that connect by hot plasma loops, and they appear in the flare impulsive phase. Our observations suggest that the fast-variation pulsations could be associated with nonthermal electrons that are periodically accelerated by the intermittent magnetic reconnection, and the 1-s period may be modulated by the coalescence instability between current-carrying loops and magnetic islands.

All other papers

To study the orbits of satellites, a galaxy could be modelled either by means of a static gravitational potential, or by live $N$-body particles. Analytic potentials allow for fast calculations, but are idealized and non-responsive. On the other hand, $N$-body simulations are more realistic, but demand higher computational cost. Our goal is to characterize the regimes in which analytic potentials provide a sufficient approximation, and those where $N$-bodies are necessary. We perform two sets of simulations using both Gala and Gadget, in order to closely compare the orbital evolution of satellites around a Milky Way-like galaxy. Focusing on the periods when the satellite has not yet been severely disrupted by tidal forces, we find that the orbits of satellites up to $10^{8} {\rm M_{\odot}}$ can be reliably computed with analytic potentials to within 5% error, if they are circular or moderately eccentric. If the satellite is as massive as $10^{9} {\rm M_{\odot}}$, errors of 9% are to be expected. However, if the orbital radius is smaller than 30 kpc, the results may not be relied upon with the same accuracy beyond 1--2 Gyr.

Dwarf spheroidal galaxies (dSphs) are known for being strongly dominated by dark matter (DM), which makes them convenient targets for investigating the DM nature and distribution. Recently, renewed interest in the dSph Leo I has resulted from claims suggesting the presence of a central supermassive black hole (BH), with mass estimates that challenge the typical expectations for dSphs, which are generally thought to host intermediate-mass black holes (IMBHs). However, Pascale et al. 2024 presented new upper limits on the BH mass, which are consistent with the range for IMBHs, solving the concerns raised in previous studies. Building on the analysis of Pascale et al. 2024, we examine the DM properties of Leo I inferred from the dynamical models of that paper. Our results indicate that Leo I is the galaxy with the highest DM density among the classical dSphs, with a central DM density (measured at a distance of $150$ pc from the galaxy centre) $\rho_{150}=35.5_{-4.7}^{+3.8}\times10^7\,M_\odot\,$pc$^{-3}$. The DM density profile has logarithmic slope $\gamma_{150}=-0.89_{-0.17}^{+0.21}$ at $150$ pc, in line with literature values. At smaller distances the DM distribution flattens into a core, with a core radius of $r_c=72^{+40}_{-32}$ pc. Combined with the small pericentric distance of Leo I's orbit in the Milky Way, the new estimate of $\rho_{150}$ makes Leo I decisive in the study of the anticorrelation between pericentre and central DM density, and suggests that the anticorrelation could be significantly steeper and more pronounced than previously estimated. Despite its DM dominance, Leo I does not emerge as the most favorable target for indirect DM detection: the inferred DM decay $D$ and annihilation $J$ factors, $\log D(0.5^{\circ})$ [GeV cm$^{-2}$] = $17.94_{-0.25}^{+0.17}$ and $\log J(0.5^{\circ})$ [GeV$^2$ cm$^{-5}$]= $18.13_{-0.18}^{+0.17}$ are consistent with previous estimates.

Despite its small mass fraction typically observed in the interstellar medium, dust plays a significant role as a key component of galaxies, affecting a wide range of properties. This review focuses specifically on how dust grains influence interstellar chemical abundances and on the processes that regulate the evolution of the galactic dust budget. I describe the main physical processes regulating dust evolution, including production by stars and other sources, destruction in supernova shocks and interstellar growth and how they are included in galactic chemical evolution models. I discuss the main effects of interstellar dust on the abundances measured in various high-redshift systems that include Damped Lyman alpha absorbers detected along the lines of sight of distant quasars and in the absorption spectra of Gamma Ray Burst afterglows. I discuss the measure of the dust mass in galaxies and review its global budget, evaluated through the study of the evolution of the comoving dust mass density, for which I present an up-to-date compilation of data chosen from the literature. Interstellar dust growth plays a critical role in regulating the dust budget, for which I present a list of evidence both in favour of it and against. The dust budget at high redshift is one aspect that requires attention to drive significant progress in the future, along with the investigation of the properties of dust in local, low-metallicity systems. Our poor theoretical knowledge of basic aspects related to dust evolution evidences the need for a new high-sensitivity space telescope operating in the far-infrared regime, still awaited by the community since the demise of Herschel.

William McClymont, Sandro Tacchella, Xihan Ji, Rahul Kannan, Roberto Maiolino, Charlotte Simmonds, Aaron Smith, Ewald Puchwein, Enrico Garaldi, Mark Vogelsberger, Francesco D'Eugenio, Laura Keating, Xuejian Shen, Bartolomeo Trefoloni, Oliver Zier

JWST has revealed the apparent evolution of the black hole (BH)-stellar mass ($M_\mathrm{BH}$-$M_\rm{\ast}$) relation in the early Universe, while remaining consistent the BH-dynamical mass ($M_\mathrm{BH}$-$M_\mathrm{dyn}$) relation. We predict BH masses for $z>3$ galaxies in the high-resolution THESAN-ZOOM simulations by assuming the $M_\mathrm{BH}$-$M_\mathrm{dyn}$ relation is fundamental. Even without live BH modelling, our approach reproduces the JWST-observed $M_\mathrm{BH}$ distribution, including overmassive BHs relative to the local $M_\mathrm{BH}$-$M_\mathrm{\ast}$ relation. We find that $M_\mathrm{BH}/M_\mathrm{\ast}$ declines with $M_\mathrm{\ast}$, evolving from $\sim$0.1 at $M_\mathrm{\ast}=10^6\,\mathrm{M_\odot}$ to $\sim$0.01 at $M_\mathrm{\ast}=10^{10.5}\,\mathrm{M_\odot}$. This trend reflects the dark matter ($f_\mathrm{DM}$) and gas fractions ($f_\mathrm{gas}$), which decrease with $M_\mathrm{\ast}$ but show little redshift evolution down to $z=3$, resulting in small $M_\mathrm{\ast}/M_\mathrm{dyn}$ ratios and thus overmassive BHs in low-mass galaxies. We use $\texttt{Prospector}$-derived stellar masses and star-formation rates to infer $f_\mathrm{gas}$ across 48,022 galaxies in JADES at $3<z<9$, finding excellent agreement with our simulation. Our results demonstrate that overmassive BHs would naturally result from a fundamental $M_\mathrm{BH}$-$M_\mathrm{dyn}$ relation and be typical of the gas-rich, dark matter-dominated nature of low-mass, high-redshift galaxies. Such overmassive, rapidly growing BHs may strongly influence the earliest stages of galaxy formation.

We present spectroscopic Ly$\alpha$ luminosity functions (LFs) at z = 5.7 and z = 6.6 based on a large 209 source sample of Ly$\alpha$ emitter (LAE) candidates identified in Subaru/Hyper Suprime-Cam narrowband imaging and confirmed with Keck II/DEIMOS spectroscopy over a multi-year observing campaign. After applying photometric and spectroscopic cuts to produce homogeneous samples, we use the resulting samples of 49 z = 5.7 and 56 z = 6.6 LAEs to compute spectroscopic LAE LFs at each redshift. We correct our LFs for incompleteness using a source-injection simulation. We find excellent agreement with current spectroscopic and photometric LAE LFs from the literature. We look for evolution over the redshift range z = 5.7-6.6. We find a strong convergence of the LFs at L(Ly$\alpha$) > $10^{43.4}$ erg s$^{-1}$. This convergence (noted in previous literature) provides strong evidence that the most luminous LAEs at z = 6.6 form ionized bubbles around themselves, allowing for greater Ly$\alpha$ transmission through the neutral intergalactic medium, which is measured as an increase in the bright end of the z = 6.6 LF over the expected evolution in the LF based on the faint end. We infer that ultraluminous LAEs may play a significant role in reionization.

Wenzer Qin, Soubhik Kumar, Priyamvada Natarajan, Neal Weiner

We propose a new mechanism for the formation of seeds of supermassive black holes at early cosmic epochs. Enhanced density fluctuations with amplitudes that are not large enough to form primordial black holes post-inflation can still lead to collapsed dark matter halos at very early times. For halos forming prior to $1+z \approx 200$, the Cosmic Microwave Background (CMB) is energetic enough to suppress the formation of molecular hydrogen, hence preventing cooling and fragmentation, as a consequence of which baryons falling into the potential well of the halo may undergo "direct collapse" into a black hole. We show using a few illustrative models how this mechanism may account for the abundance of high-redshift black holes inferred from observations by the James Webb Space Telescope while remaining consistent with current limits from CMB spectral distortions. Limits on the primordial power spectrum are also derived by requiring that the universe not reionize too early.

Line intensity mapping (LIM) is an emerging technique for probing the aggregate emission of a spectral line from all sources, without requiring individual detections. Through the wavelength-redshift relation, one can map the line-of-sight evolution of the line emission that traces the underlying large-scale structure in a spectral-imaging survey. In this work, we present a new technique -- feature intensity mapping -- as an extension of the LIM formalism to map broad spectral features in 3D, rather than the narrow emission lines typically targeted by LIM. By accounting for the convolution of spectral features with the instrument's spectral response across redshift, our technique enables simultaneous constraints on the redshift-dependent emission from multiple features. This approach enables 3D intensity mapping with some of the brightest features in galaxies' infrared spectra: the polycyclic aromatic hydrocarbon (PAH) emission bands. We forecast the detectability of PAH signals using feature intensity mapping with the ongoing SPHEREx mission in the near-infrared and the proposed PRIMA mission in the far-infrared. We find that $S/N$ of $\gtrsim 10$ per redshift bin of widths $\Delta z = 0.1$ and $0.5$ can be achieved at $z < 0.5$ and $1 < z < 5$ with SPHEREx and PRIMA, respectively, for multiple PAH features, suggesting a promising prospect for mapping the aggregate PAH emission at cosmological distances with upcoming datasets.

Pulsar Timing Array (PTA) searches for nHz gravitational-wave backgrounds (GWBs) typically model time-correlated noise by assuming a diagonal covariance in Fourier space, neglecting inter-frequency correlations introduced by the finite observation window. We show that this diagonal approximation can lead to biased estimates of spectral parameters, especially for the common red process that represents the GWB. To address these limitations, we present a method that (i) computes the time-domain autocorrelation on a coarse grid using a fast Fourier transform (FFT), (ii) interpolates it accurately to the unevenly sampled observation times, and (iii) incorporates it into a low-rank likelihood via the Sherman--Morrison--Woodbury identity. Using both analytic covariance comparisons and end-to-end simulations inspired by the NANOGrav 15-year dataset, we demonstrate that our method captures frequency correlations faithfully, avoids Gibbs ringing, and recovers unbiased spectral parameters with modest computational cost. As PTA datasets increase in sensitivity and complexity, our approach offers a practical and scalable path to fully accurate covariance modeling for current and future analyses.

Benedetta Mestichelli, Michela Mapelli, Filippo Santoliquido, Manuel Arca Sedda, Marica Branchesi, Lavinia Paiella, Guglielmo Costa, Giuliano Iorio, Matthew Mould, Veronika Lipatova, Boyuan Liu, Ralf S. Klessen

Population III (Pop.$~$III) stars are expected to be massive and to undergo minimal mass loss due to their lack of metals, making them ideal progenitors of black holes and neutron stars. Here, we investigate the formation and properties of binary neutron star (BNS) and black hole-neutron star (BHNS) mergers originating from Pop.$~$III stars, and compare them to their metal-enriched Population II (Pop.$~$II) counterparts, focusing on their merger rate densities (MRDs), primary masses and delay times. We find that, despite the high merger efficiency of Pop.$~$III BNSs and BHNSs, their low star formation rate results in a MRD at least one order of magnitude lower than that of Pop.$~$II stars. The MRD of Pop.$~$III BNSs peaks at redshift $z\sim15$, attaining a value $\mathcal{R}_{\rm BNS}(z\sim15) \sim 15\,\rm Gpc^{-3}\,yr^{-1}$, while the MRD of Pop.$~$III BHNSs is maximum at $z\sim13$, reaching a value $\mathcal{R}_{\rm BHNS}(z\sim13) \sim 2\,\rm Gpc^{-3}\,yr^{-1}$. Finally, we observe that the black hole masses of Pop.$~$III BHNS mergers have a nearly flat distribution with a peak at $\sim 20\,\rm M_{\odot}$ and extending up to $\sim 50\,\rm M_{\odot}$. Black holes in Pop.$~$II BHNS mergers show instead a peak at $\lesssim 15\,\rm M_{\odot}$. We consider these predictions in light of recent gravitational-wave observations in the local Universe, finding that a Pop.$~$III origin is preferred relative to Pop.$~$II for some events.

Modeling energy-dependent X-ray pulse profiles from rotation-powered millisecond pulsars observed with NICER has emerged as a promising avenue for measuring neutron star radii and probing the equation of state of cold, ultra-dense matter. However, pulse profile models have often required an unwieldy number of parameters to account for complex surface emission geometries, introducing the risk of overfitting and degeneracies. To explore the number of model parameters that can be inferred uniquely, we perform a quantitative assessment of the information content in X-ray pulse profiles by applying Fourier methods. We determine the number of independent observables that can be reliably extracted from the pulse shapes, as well as from complementary X-ray spectral data obtained with XMM-Newton, for key NICER targets. Our analysis provides a framework for evaluating the match between model complexity and data constraints. It also demonstrates the importance of incorporating in the model the pulsed components of the magnetospheric non-thermal emission, which often contributes significantly to the observed spectra. Our results highlight limitations in previous inferences of neutron-star radii from NICER observations, which have incorporated model complexity not supported by the data.

Wilfried Mercier, Boris Sindhu Kalita, Marko Shuntov, Rafael C. Arango-Toro, Olivier Ilbert, Laurence Tresse, Yohan Dubois, Clotilde Laigle, Hossein Hatamnia, Nicolas McMahon, Andreas Faisst, Isa Cox, Maxime Trebitsch, Leo Michel-Dansac, Si-Yue Yu, Michaela Hirschmann, Marc Huertas-Company, Arianna Long, Anton Koekemoer, Grégoire Aufort, Joseph Lewis, Ghassem Gozaliasl, R. Michael Rich, Jason Rhodes, Henry Joy McCracken, Caitlin Casey, Jeyhan Kartaltepe, Brant Robertson, Maximilien Franco, Daizhong Liu, Hollis Akins, Natalie Allen, Sune Toft

Clumps in the rest-frame UV emission of galaxies have been observed for decades. Since the launch of the James Webb Space Telescope (JWST), a large population is detected in the rest-frame near-infrared (NIR), raising questions about their formation mechanism. We investigate the presence and properties of NIR over-densities (hereafter substructures) in star-forming and quiescent galaxies at 1 < z < 4 to understand their link to the evolution of their host galaxy. We identify substructures in JWST/NIRCam F277W and F444W residual images at a rest-frame wavelength of 1 um. The fraction of galaxies with substructures with M* > 10^9 Msun has been steadily decreasing with cosmic time from 40% at z = 4 to 10% at z = 1. Clumps, the main small substructures in the rest-frame NIR, are the most common type and are much fainter (2% of the flux) than similar UV clumps in the literature. Nearly all galaxies at the high-mass end of the main sequence (MS), starburst, and green valley regions have substructures. However, we do not find substructures in low-mass galaxies in the green valley and red sequence. Although massive galaxies on the MS and in the green valley have a 40% probability of hosting multiple clumps, the majority of clumpy galaxies host only a single clump. The fraction of clumpy galaxies in the rest-frame NIR is determined by the stellar mass and SFR of the host galaxies. Its evolution with redshift is due to galaxies moving towards lower SFRs at z < 2 and the build-up of low-mass galaxies in the green valley and red sequence. Based on their spatial distribution in edge-on galaxies, we infer that most of substructures are produced in-situ via disk fragmentation. Galaxy mergers may still play an important role at high stellar masses, especially at low SFR.

Rory Smith, Stephanie Tonnesen, Katarina Kraljic, Paula Calderon-Castillo, Antonino Marasco, Yara Jaffe, Benedetta Vulcani, Bianca M. Poggianti

Context: In dense environments, disk galaxies can be subjected to tidal interactions with other galaxies and/or ram pressure stripping. Some morphological features are clearly associated with one or the other interaction (e.g. tidal bridges vs long one-sided linear gas tails). But, under certain circumstances, both mechanisms can result in morphological features that could be confused, such as lopsided or asymmetric disks and unwinding spiral arms. Aims: Our aim is to develop new measures for application to asymmetric galaxies of this type that distinguish gravitational-only tidal interactions from ram pressure stripping, and that can be applied directly to simulations, and potentially to observations. Methods: We define a new measure for galaxies called the Size-Shape Difference (SSD) measure. This measure is sensitive to differences in the size and shape of a younger stellar population (<200 Myr) compared to that of an intermediate age stellar population (200-400 Myr). We use numerical simulations of galaxies undergoing gravitational-only tidal interactions and/or undergoing ram pressure stripping to test the measure. Results: Because ram pressure tends to directly alter the gas distribution, the younger stellar population (which best traces out the gas distribution) tends to change shape and morphology with respect to the intermediate age population. The SSD measure is sensitive to this change, and we find it can effectively distinguish between ram pressure and gravitational-only tidal encounters. In fact, we find it is even more effective when a combination of a tidal interaction and ram pressure has occurred together, as may arise in dense environments. As tidal interactions tend to enhance the spiral structure in disk galaxies, the effectiveness of the SSD measure is further enhanced when combined with a measure of the strength of the spiral arms.

R. Anirudh, Melanie Kaasinen, Gergö Popping, Desika Narayanan, Karolina Garcia, Dariannette Valentin-Martinez

(Abridged) Resolved observations of the CO emission from $z=1-3$ star-forming galaxies are becoming increasingly common, with new high-resolution surveys on the horizon. We aim to inform the interpretation of this resolved CO emission by creating synthetic observations and testing to what extent routinely observed CO transitions can be used to trace H$_2$ across galaxy disks. To this end, we extract $z=1-3$ massive star-forming galaxies (on and above the main sequence) from the SIMBA cosmological simulation and predict their spatially resolved CO(1$-$0)-to-CO(5$-$4) emission using the $\texttt{SLICK}$ pipeline, which combines sub-resolution modeling of the cloud population with the DESPOTIC spectral line calculation code. We find that the CO(1$-$0)-to-H$_2$ ratio ($\alpha_{\rm CO}$) varies significantly within these galaxy disks$-$from values of $\sim1-5$ $\mathrm{M_\odot}$ (K km s$^{-1}$ pc$^2$)$^{-1}$ in the central 1-3 kpc of the most massive galaxies to $>100$ $\mathrm{M_\odot}$ (K km s$^{-1}$ pc$^2$)$^{-1}$ at $\sim$ 15 kpc. Thus, the use of a single $\alpha_{\rm CO}$ to derive the H$_2$ surface density leads to severe underestimates of the H$_2$ contribution in its outskirts. As expected, higher-$J$ CO transitions trace molecular gas in the centers at higher densities, whereas CO(1$-$0) better traces the more diffuse, extended molecular gas. We see significant variations in the CO excitation, with CO(3$-$2)/CO(1$-$0) line luminosity ratios of the most massive galaxies at $z\sim2$ declining from $\sim$ 0.9 in the galaxy centers to $\sim$ 0.1 in the outskirts. On average, line ratios increase substantially toward higher redshifts and lower galaxy stellar masses. We predict that tracing molecular gas with CO beyond 3-5 kpc of cosmic noon galaxies will be challenging with current facilities due to the drastic increase in $\alpha_{\rm CO}$.

Jack Stethem (1), Christopher S. Kochanek (1, 2), Anya Phillips (3), Lyra Cao (4), Marc Pinsonneault (1, 2) ((1) Department of Astronomy, The Ohio State University, (2) Center for Cosmology and AstroParticle Physics, The Ohio State University, (3) Center for Astrophysics - Harvard and Smithsonian, (4) Department of Physics and Astronomy, Vanderbilt University)

Rotational variables are stars that vary in brightness due to star spots modulated by rotation. They are probes of stellar magnetism, binarity, and evolution. Phillips et al. (2023) explored distinct populations of ~50,000 high-amplitude rotational variables from the All-Sky Automated Survey for Supernovae (ASAS-SN), examining correlations between stellar rotation, binarity, and activity. Here, we carry out a similar analysis of ~50,000 much lower amplitude Kepler rotational variables. The Kepler population is dominated by slowly rotating, single, main sequence stars, with a striking absence of the rapidly rotating main sequence group in the ASAS-SN sample. The binary fractions of the Kepler rotators are significantly lower than for the ASAS-SN systems and they are significantly less spotted, as expected from their lower amplitudes. The scope of these statistical surveys will dramatically increase in the near future.

G. Yu. Mozgunov (1), A. S. Pozanenko (1, 2), P. Yu. Minaev (1), I. V. Chelovekov (1), S. A. Grebenev (1), A. G. Demin (3), A. V. Ridnaya (3), D. S. Svinkin (3), Yu. R. Temiraev (4), D. D. Frederiks (3) ((1) Space Research Institute, Moscow, (2) National Research University `Higher School of Economics', Moscow, (3) Ioffe Institute, St. Petersburg, (4) GlowByte Company, Moscow)

We have searched for ultra-long (> 100 s) gamma-ray transients in the data from the anticoincidence shield (ACS) of the SPI gamma-ray spectrometer onboard the INTEGRAL orbital observatory and classified them by machine learning methods. We have found about 4364 candidates for such events in the SPI-ACS data by the `blind' threshold search method. We have developed an algorithm for automatic processing of their light curves that distinguishes a candidate for transients on various time scales and allows its duration and fluence to be determined. The algorithm has been applied to calculate (and compare) the fluxes in the light curves recorded by various INTEGRAL detectors: IREM, SPI-ACS, SPI, ISGRI, and PICsIT. These fluxes have been used to train the classifier based on gradient boosting. Subsequently, we have performed a cluster analysis of the candidates found by the dimensionality reduction and clustering methods. In conclusion we have compared the remaining candidates with the data from the Konus-WIND gamma-ray detectors. Thus, we have confirmed 16 candidates for astrophysical transients, including four candidates for ultra-long gamma-ray bursts from the events detected by the SPI-ACS detector. Out of the probable events, but unconfirmed by other experiments, up to 270 events can be classified as real gamma-ray bursts.

Mid-infrared (mid-IR) observations of Near-Earth Objects (NEOs) have historically been a valuable tool for understanding their physical properties. However, the current state of mid-IR instruments on ground-based telescopes places several limitations on performing thermal characterization of NEOs. The complexity of maintaining these instruments in operational conditions on telescopes has led to their decommissioning. Here, we present the first science commissioning observations out to 12.5 microns from the upgraded Mid-Infrared Spectrograph and Imager (MIRSI) at the NASA-IRTF. We obtained 42 observations of 31 NEOs and derived their diameters and albedos. Since MIRSI allows simultaneous optical observations with its MIRSI Optical Camera (MOC), we were able to determine the absolute magnitude for most of the targets at the time of the thermal acquisition. We present ejecta characterization for the Didymos system from observations made 11 hours and 9 days after the Double Asteroid Redirection Test (DART) impact. We present albedo and size measurements for (98943) Torifune 2001 CC21, the fly-by target of the Japanese Extended Hayabusa2 Mission. We also highlight several applications that the MIRSI system will provide for future airless body characterization, such as constraining thermal inertia from simultaneous optical and thermal lightcurves. This work also demonstrates the importance of having MIRSI as an available rapid-response instrument for planetary defense purposes.

E. Bañados (1), V. Le Brun (2), S. Belladitta (1 and 3), I. Momcheva (1), D. Stern (4), J. Wolf (1), M. Ezziati (2), D. J. Mortlock (5 and 6), A. Humphrey (7 and 8), R. L. Smart (9 and 10), S. L. Casewell (11), A. Pérez-Garrido (12), B. Goldman (13 and 14), E. L. Martín (15 and 16), A. Mohandasan (17), C. Reylé (17), C. Dominguez-Tagle (15 and 16), Y. Copin (18), E. Lusso (19 and 20), Y. Matsuoka (21), K. McCarthy (4), F. Ricci (22 and 23), H.-W. Rix (1), H. J. A. Rottgering (24), J.-T. Schindler (25), J. R. Weaver (26), A. Allaoui (2), T. Bedrine (2), M. Castellano (23), P.-Y. Chabaud (2), G. Daste (2), F. Dufresne (2), J. Gracia-Carpio (27), M. Kümmel (28), M. Moresco (29 and 3), M. Scodeggio (30), C. Surace (2), D. Vibert (2), A. Balestra (31), A. Bonnefoi (2), A. Caillat (2), F. Cogato (29 and 3), A. Costille (2), S. Dusini (32), S. Ferriol (18), E. Franceschi (3), W. Gillard (33), K. Jahnke (1), D. Le Mignant (2), S. Ligori (9), E. Medinaceli (3), G. Morgante (3), F. Passalacqua (34 and 32), K. Paterson (1), S. Pires (35), C. Sirignano (34 and 32), I. T. Andika (36 and 37), H. Atek (38), D. Barrado (39), S. Bisogni (30), C. J. Conselice (40), R. Decarli (3), H. Dole (41), T. Dupuy (42), A. Feltre (19), S. Fotopoulou (43), B. Gillis (42), X. Lopez Lopez (29 and 3), M. Onoue (44 and 45), G. Rodighiero (34 and 31), N. Sedighi (15 and 16), F. Shankar (46), M. Siudek (47 and 48), L. Spinoglio (49), D. Vergani (3), G. Vietri (30), F. Walter (1), G. Zamorani (3), M. R. Zapatero Osorio (39), J.-Y. Zhang (16 and 15), M. Bethermin (14), N. Aghanim (41), B. Altieri (50), A. Amara (51), S. Andreon (52), C. Baccigalupi (53 and 54 and 55 and 56), M. Baldi (57 and 3 and 58), S. Bardelli (3), A. Basset (59), P. Battaglia (3), A. Biviano (54 and 53), A. Bonchi (60), D. Bonino (9), E. Branchini (61 and 62 and 52), M. Brescia (63 and 64), J. Brinchmann (7 and 65), S. Camera (66 and 67 and 9), V. Capobianco (9), C. Carbone (30), J. Carretero

We demonstrate the potential of Euclid's slitless spectroscopy to discover high-redshift (z>5) quasars and their main photometric contaminant, ultracool dwarfs. Sensitive infrared spectroscopy from space is able to efficiently identify both populations, as demonstrated by Euclid Near-Infrared Spectrometer and Photometer Red Grism (NISP RGE) spectra of the newly discovered z=5.404 quasar EUCL J181530.01+652054.0, as well as several ultracool dwarfs in the Euclid Deep Field North and the Euclid Early Release Observation field Abell 2764. The ultracool dwarfs were identified by cross-correlating their spectra with templates. The quasar was identified by its strong and broad CIII] and MgII emission lines in the NISP RGE 1206-1892 nm spectrum, and confirmed through optical spectroscopy from the Large Binocular Telescope. The NISP Blue Grism (NISP BGE) 926-1366 nm spectrum confirms CIV and CIII] emission. NISP RGE can find bright quasars at z~5.5 and z>7, redshift ranges that are challenging for photometric selection due to contamination from ultracool dwarfs. EUCL J181530.01+652054.0 is a high-excitation, broad absorption line quasar detected at 144 MHz by the LOw-Frequency Array (L144=4e25 W/Hz). The quasar has a bolometric luminosity of 3e12 Lsun and is powered by a 3.4e9 Msun black hole. The discovery of this bright quasar is noteworthy as fewer than one such object was expected in the ~20 deg2 surveyed. This finding highlights the potential and effectiveness of NISP spectroscopy in identifying rare, luminous high-redshift quasars, previewing the census of these sources that Euclid's slitless spectroscopy will deliver over about 14,000 deg2 of the sky.

T. Joseph W. Lazio (Jet Propulsion Laboratory, California Institute of Technology)

The Deep Space Network (DSN) is the primary means of commanding, tracking, and receiving data from all of NASA's deep space missions, as well as a number of deep space missions operated by other international space agencies. The current number of missions enabled by the DSN is approximately 40 missions, but there has been concern about the level of "over-subscription" of the DSN, namely that the number of missions currently using the DSN is larger than can be enabled reasonably. This manuscript assesses the maximum number of missions that could be enabled, based on recent performance and with the constraint that the total number of hours used per week does not exceed the available number of DSN antenna-hours. Three different models are considered, and the maximal number of missions that could be enabled ranges between approximately 40 missions and 70 missions, assuming that there continues to be approximately six Mars missions and that those Mars missions continue to make use of the DSN's multiple spacecraft per antenna (MSPA) capability. Crucially, the conclusion that an approximately 50% growth in the DSN mission suite rests on the assumption that the DSN antennas are "interchangeable," but they are not, with some spacecraft able to use only certain antennas. Efforts to make the DSN antennas more "interchangeable," primarily in their transmitter and receiver suites, would be an effective means of ensuring expanded capability. Additional findings from this work are that, while additional use of the MSPA capability might appear to be a promising means for increasing the mission suite, there appear to be no locations in the Solar System, other than Mars, for which it would be effective.

This work proposes a multiple machine learning method (MMLM) aiming to improve the accuracy and robustness in the analysis of star clusters. The MMLM performance is evaluated by applying it to the reanalysis of the old binary cluster candidate - NGC 1605a and NGC 1605b - found by Camargo (2021) (hereafter C21). The binary cluster candidate is analyzed by employing a set of well established machine learning algorithms applied to the Gaia-EDR3 data. Membership probabilities and open clusters (OCs) parameters are determined by using the clustering algorithms pyUPMASK, ASteCA, Kmeans, GMM, and HDBSCAN. In addition, a KNN smoothing algorithm is implemented to enhances the visualization of features like overdensities in the 5D space and intrinsic stellar sequences on the color-magnitude diagrams (CMDs). The method validates the clusters' parameters previously derived, however, suggests that their probable members-stars are distributed over a wider overlapping area. Finally, a combination of the elbow method, t-SNE, kmeans, and GMM algorithms group the normalized data into 6 clusters, following C21. In short, these results confirm NGC1605a and NGC1605b as genuine OCs and reinforce the previous suggestion that they form an old binary cluster in an advanced stage of merging after a tidal capture during a close encounter. Thus, MMLM has proven to be a powerful tool that helps to obtain more accurate and reliable clusters parameters and its application in future studies may contribute to a better characterization of the Galaxy's star cluster system.

Indrajit Paul, R. Vasanth Kashyap, Tuhin Ghosh, Rainer Beck, Luke Chamandy, Srijita Sinha, Anvar Shukurov

To explore the spatial variations of the regular (mean) magnetic field of the Andromeda galaxy (M31), we use Fourier analysis in azimuthal angle along four rings in the galaxy's plane. Earlier analyses indicated that the axisymmetric magnetic field (azimuthal Fourier mode $m=0$) is sufficient to fit the observed polarization angles in a wide range of galactocentric distances. We apply a Bayesian inference approach to new, more sensitive radio continuum data at $\lambda \lambda3.59$, $6.18$, and $11.33$ cm and the earlier data at $\lambda 20.46$ cm to reveal sub-dominant contributions from the modes $m=1$, 2, and 3 along with a dominant axisymmetric mode. Magnetic lines of the axisymmetric mode are close to trailing logarithmic spirals which are significantly more open than the spiral arms detectable in the interstellar dust and neutral hydrogen. The form of the $m=0$ mode is consistent with galactic dynamo theory. Both the amplitudes and the pitch angles of the higher azimuthal modes ($m>1$) vary irregularly with $r$ reflecting local variations in the magnetic field structure. The maximum strength of the mean magnetic field of $1.8\text{--}2.7\mu$G (for the axisymmetric part of the field) occurs at $10\text{--}14$ kpc but we find that its strength varies strongly along the azimuth; this variation gives rise to the $m=1$ mode. We suggest a procedure of Bayesian inference which is independent of the specific nature of the depolarization and applies when the magneto-ionic layer observable in polarized emission is not symmetric along the line of sight because emission from its far side is completely depolarized.

C. Mac Cormack (1, 2), S. B. Shaik (3), P. Hess (4), R. Colaninno (4), T. Nieves-Chinchilla (1) ((1) Heliospheric Physics Laboratory, Heliophysics Science Division, NASA Goddard Space Flight Center, MD, USA, (2) The Catholic University of America, Washington, DC, USA, (3) George Mason University, Fairfax, VA, USA, (4) U.S. Naval Research Laboratory, Washington, D.C., USA)

With the combination of observations from in-situ and high-resolution remote sensing instruments, the Solar Orbiter (SolO) mission has become a particularly valuable mission for studying the inner heliosphere. With a field of view (FOV) of 40$^\circ$ to the east of the Sun, the Solar Orbiter Heliospheric Imager (SoloHI) is one of the six remote sensing instruments on board the Solar Orbiter (SolO) spacecraft. SoloHI's high-resolution imaging observations of the heliosphere and its higher cadence than previous generations of heliospheric imagers make it a perfect candidate to perform and complement studies of CMEs evolution through the heliosphere. In this work, we present the first prominence material detected by SoloHI along with other remote sensing instruments. We report that as the CME propagates out in the heliosphere, the associated filament material reaches a heliocentric height of $\sim$122.5$\,$R$_{\odot}$ ($\sim$0.57$\,$au).

Kara L. Kniezewski, Emily I. Mason, Daniel J. Emmons, Kyle E. Fitch, Seth H. Garland

Solar flares are intense bursts of electromagnetic radiation, which occur due to a rapid destabilization and reconnection of the magnetic field. While pre-flare signatures and trends have been investigated from magnetic observations prior to flares for decades, analysis which characterizes the variability of the magnetic field in the hours prior to flare onset has not been included in the literature. Here, the 3D magnetic field is modeled using a Non-Linear Force Free Field extrapolation for 6 hours before and 1 hour after 18 on-disk solar flares and flare quiet windows for each active region. Parameters are calculated directly from the magnetic field from two field isolation methods: the "Active Region Field", which isolates field lines where the photospheric field magnitude is $\geq$ 200 Gauss, and the "High Current Region", which isolates field lines in the 3D field where the current, non-potential field, twist, and shear exceed pre-defined thresholds. For this small pool of clean events, there is a significant increase in variation starting 2-4 hours before flare onset for the current, twist, shear, and free energy, and the variation continues to increase through the flare start time. The current, twist, shear, and free energy are also significantly stronger through the lower corona, and their separation from flare quiet height curves scales with flare strength. Methods are proposed to combine variation of the magnetic fields with variation of other data products prior to flare onset, suggesting a new potential flare prediction capability.

Chandreyee Sengupta, Tom C. Scott, Hao Chen, Hyein Yoon, Yogesh Chandola, Mengtian Li, Gyula I. G. Józsa, O. Ivy Wong, Yin-Zhe Ma, Patricio Lagos, Ruta Kale, Denis Tramonte

A large catalogue of low surface brightness galaxies (LSBGs) from the Dark Energy Survey showed significant clustering around nearby galaxy groups and clusters. Using the HIPASS survey, we tried to determine the redshift of a sub-sample of these LSBGs and determine whether they were members of the groups they were projected near, but this was hampered by HIPASS's high spectral rms. This letter reports on MeerKAT H I observations to determine the redshifts of 52 LSBG candidates projected in the vicinity of two groups from our previous HIPASS study. The main goal is to investigate and ascertain whether these LSBGs are genuine group members. H I was detected with MeerKAT and redshifts were determined for only five of the 52 candidates within a velocity range of $\pm$ 2500 km/s of their respective group velocities. All five H I detections were blue LSBGs and two of them were confirmed to be ultradiffuse galaxies (UDGs). Both these UDGs were group members, while the other three detections were either foreground or background galaxies. In this letter we explore scenarios that can explain the 90% non-detection. MeerKAT's excellent sensitivity allows us to conclude that the majority of the non-detected candidates, particularly the blue galaxies, are not group members but lie at higher redshifts. However, this still leaves the open question as why Tanoglidis LSBG candidates, in particular the red ones, appear to be clustered in projection around nearby groups.

Where in the present-day Milky Way should we search for the remnants of its earliest stars? We address this question using the DORCHA (Gaelic for Dark; DUR-uh-khuh) suite: a set of 25 high-resolution, dark-matter-only cosmological zoom-in simulations of Milky Way analogue (MWA) haloes evolved to $z=0$. Of these, 15 are isolated and the rest are in pairs, similar to the MW and M31. By identifying and tagging the most bound material in high-redshift ($z\geq5$) progenitor haloes -- those likely to host early star formation -- we track the present-day phase-space distribution of this ancient component. We find that this material is highly centrally concentrated at $z=0$, with 90 - 100 per cent residing within $r \lesssim 15\,h^{-1}\,\mathrm{kpc}$. It exhibits steep density profiles ($\rho\propto\,r^{-4}$), low velocity dispersions ($\sigma_r / \sigma_{\rm max} \lesssim 0.6$), and radially biased orbits ($\beta \gtrsim 0.5$ for $r \gtrsim 0.1\,R_{200}$), consistent with a relaxed, centrally embedded population. These results hold across haloes with diverse formation histories and environments, suggesting that the dynamical signature of early progenitors is robust to later mergers and interactions. Our findings imply that the fossil record of the first generations of stars -- including Population III and extremely metal-poor stars -- should be sought in the innermost regions of the Milky Way, where they retain distinctive kinematic imprints. While these stellar populations may overlap, we caution that low metallicity does not uniquely identify ancient stars, nor vice versa. The DORCHA suite thus provides a physically motivated baseline for interpreting observations from Galactic Archaeology surveys targeting the bulge and inner halo.

Primordial black holes (PBHs) arise from the collapse of density perturbations in the early universe and serve as a dark matter (DM) candidate and a probe of fundamental physics. There remains an unconstrained ``asteroid-mass'' window where PBHs of masses $10^{17} {\rm g} \lesssim M \lesssim 10^{23} {\rm g}$ could comprise up to $100\%$ of the dark matter. Current $e^{\pm}$ Hawking radiation constraints on the DM fraction of PBHs are set by comparing observed spatial- and time-integrated cosmic ray flux measurements with predicted Hawking emission fluxes from the galactic DM halo. These constraints depend on cosmic ray production and propagation models, the galactic DM density distribution, and the PBH mass function. We propose to mitigate these model dependencies by developing a new local, time-dependent Hawking radiation signature to detect low-mass PBHs transiting through the inner Solar System. We calculate transit rates for PBHs that form with initial masses $M \lesssim 5\times10^{17}\text{g}$. We then simulate time-dependent positron signals from individual PBH flybys as measured by the Alpha Magnetic Spectrometer (AMS) experiment in low-Earth orbit. We find that AMS is sensitive to PBHs with masses $M\lesssim 2\times10^{14} \, {\rm g}$ due to its lower energy threshold of $500 \, {\rm MeV}$. We demonstrate that a dataset of daily positron fluxes over the energy range $5-500 \, {\rm MeV}$, with similar levels of precision to the existing AMS data, would enable detection of PBHs drawn from present-day distributions that peak within the asteroid-mass window. Our simulations yield ${\cal O} (1)$ detectable PBH transits per year across wide regions of parameter space, which may be used to constrain PBH mass functions. This technique could be extended to detect $\gamma$-ray and X-ray Hawking emission to probe further into the asteroid-mass window.

Quiet Sun magnetism plays an important role in the global energy balance of the solar atmosphere. The study of the magnetism of the quiet Sun has been a major effort in the last decades, and, as a result, very important advances in our knowledge have been achieved. The recent commissioning of the Daniel K. Inouye Solar Telescope, the largest ground-based solar telescope in the world today, provides us with data of high spectropolarimetric quality and high spatial resolution. The combination of data acquired with the ViSP instrument on the DKIST telescope around 630.1 nm and multiline inversions enables the achievement of exceptionally high spatial resolution in the solar atmosphere, both at the surface and in depth. In this paper we determine the log(gf) values of the spectral lines around 630.1 nm observed with the ViSP instrument and search for the best multi-inversion strategy, by means of MHD simulations, to analyze quiet Sun magnetism data within such spectral range. We present results of quiet Sun magnetism with special emphasis on its evolution with optical depth. In the internetwork, we find the decrease with height of the averaged magnetic field strength and the zone from which the fields become more horizontal. Additionally, we explore the depth-dependent evolution of the networks canopy, uncovering intriguing insights into its thermodynamic and magnetic properties. Notably, we detect a temperature enhancement in very localized areas of the vicinity of the network and distinguish mass motions with opposing velocities inside the flux tubes.

Ian D. Roberts, Michael L. Balogh, Visal Sok, Adam Muzzin, Michael J. Hudson, Pascale Jablonka

We report the discovery of COSMOS2020-635829 as a likely jellyfish galaxy undergoing to ram pressure stripping in a (proto)cluster at $z > 1$. High-resolution imaging from the James Webb Space Telescope reveals a symmetric stellar disk coupled to a unilateral tail of star-forming knots to the south. We show that these extra-planar continuum sources are embedded within an ionized gas tail that is kinematically connected to the disk of COSMOS2020-635829. Likely representing the highest-redshift discovery of a ram pressure stripped ionized gas tail. The tail sources are characterized by extremely young stellar populations ($\lesssim 100\,\mathrm{Myr}$), have stellar masses of ${\sim}10^8\,\mathrm{M_\odot}$, and star formation rates of $0.1\text{--}1\,\mathrm{M_\odot\,yr^{-1}}$. This work reinforces the notion that ram pressure stripping can perturb group and cluster galaxies at $z > 1$ and likely contributes to environmental quenching even near Cosmic Noon.

Koji Sugitani, Jan G. A. Wouterloot, Harriet Parsons, Sarah Graves, Takayoshi Kusune, Archana Soam

We carried out polarimetric observations of the bright-rimmed cloud IC 1396E/SFO 38 with SCUBA-2/POL-2 to study the effect of ultraviolet (UV) light on its structure and magnetic field. This bright-rimmed cloud appears optically to be one single cloud illuminated by the UV light from the excited star of IC 1396, however our Stokes I image and 13CO(J=3-2) archival data suggest that this cloud is not a simple, single structure, but appears to be composed of two parts on first glance; a head part with wings and a tail, and a north-west extension this http URL molecular clouds are generally filamentary and it seems likely that the initial structures of bright-rimmed clouds are expected to be also generally elongated, we examined the possibility that the structure was created from a single elongated cloud by the UV impact. We compared the cloud structure with a simulation study that investigated the evolution of prolate clouds exposed to the UV radiation from various directions and found that this apparent two-part structure could be reproduced in a situation where a single filamentary cloud is obliquely illuminated by UV light. The magnetic field directions of the cloud are different from the ambient field direction, demonstrating the field reconfiguration due to the UV impact. The distortion or pinch of the magnetic field is seen toward the cloud head, where an intermediate-mass star cluster is located, suggesting gravitational contraction. We roughly estimated the magnetic strength and stability in three parts of the cloud and found that the cloud head is most likely to be supercritical.

We study the thermal and spectral properties of irradiated circumbinary disks (CBDs) around binary black holes (BBHs), using analytic, hydrogen-based opacity models that capture key dependencies on temperature, density, and ionization. We solve the vertically hydrostatic energy balance equations with Rosseland mean opacities from free-free absorption, bound-free absorption, and electron scattering processes, combined with ionization fractions derived by the Saha equation. Four opacity models are considered, including a reference model with no physical opacity, constructed by Lee et al. (2024), and three physically motivated alternatives. The midplane temperature profiles show significant variation across models, while the surface temperature remains largely unchanged in regions dominated by viscous heating. Opacity effects become pronounced in the outer disk, where irradiation reprocessing shapes the IR-optical continuum. Inclusion of bound-free opacity introduces a noticeable flattening and a mid-frequency peak in the spectral energy distribution. We also compute spectra of a triple disk system including the CBD and two accreting minidisks. The high-frequency peak arises from the hot minidisks, while the low-frequency excess originates from irradiated outer CBD layers. By comparing model spectra with detection limits of Subaru, JWST, and Swift, we find that systems within $\sim10$ Mpc may show detectable IR excess from the CBD. Our results highlight the need for accurate opacity modeling to interpret electromagnetic signatures of black hole mergers and support future integration of opacity tables with metallicity.

Gravitational waves (GWs) accompanied by electromagnetic (EM) counterparts provide a novel methodology to measure the Hubble constant ($H_0$), known as bright sirens. However, the rarity of such multi-messenger events limits the precision of the $H_0$ constraint. Recently, the newly-discovered quasi-periodic eruptions (QPEs) show intriguing evidence of a stellar-mass companion captured by a supermassive black hole (SMBH) in an extreme mass-ratio inspiral (EMRI), which is the most promising sources of the space-based GW detectors, such as Laser Interferometer Space Antenna (LISA). Here we study the secular orbital evolution of QPE systems under different theoretical frameworks, and assess their GWs detectability by LISA. We find that LISA can detect the EMRI signals of two or three known QPE systems with a four-year observation. One EMRI event can measure the $H_0$ with an uncertainty of $3\%-8\%$ at a 68.3 percent confidence level, while a few events will reduce it below $2\%$. The EMRI events surpass current bright siren limitations and offer an independent pathway to resolve the Hubble constant tension.

A. Waszewski, J. S. Morgan, R. Ekers, M. Johnston-Hollitt, M. C. M. Cheung, N. D. R. Bhat, R. Chhetri, S. C. Fu

We present a study of the solar wind over different periods of the solar cycle, specifically focussing on the minimum between solar cycles 24 and 25, and the active, ascending phase of solar cycle 25. With the use of interplanetary scintillation (IPS) data taken by the Murchison Widefield Array (MWA) from mid-2019 and early 2023, we have sampled over a wide range of solar latitudes and elongations, probing a large section of the surrounding heliosphere from $\sim$90 to 140R$_\odot$. The MWA observations provide the highest density of sampled IPS radio sources to date, allowing for an investigation into the latitudinal dependence of the scattering effect caused by the solar wind on a radio source as observed throughout the solar cycle. We find our measurements during periods of heightened solar activity are consistent with a spherically symmetric solar wind. On the other hand, with a reduction in solar activity we find the solar wind density inherits a latitudinal dependence. As is consistent with prior studies, an elliptical function better represents the transition from poles to equator, although we find a more exaggerated sigmoid shaped curve is required to represent the low- to mid-latitude transition region during the minimum of solar cycle 24. We find for a heliospheric distance range of 108 - 123R$_{\odot}$ the reduction ratio between the equator and the southern pole is 1.62$\pm$0.02.

We present a physically motivated model of dark energy (DE) rooted in the topological structure of the Quantum ChromoDynamic (QCD) vacuum. In this framework, DE emerges as the Universe expands, from the difference in vacuum energy between the expanding Universe and Minkowski spacetime, driven by QCD vacuum topological sectors. This leads to a modification of the dark energy term in the Friedmann equation, which then scales with the Hubble parameter as $\rho_{\rm DE}(t) \propto H(t)$ when dark energy dominates the expansion of the Universe. The QCD scale, $\Lambda_{\rm QCD} \sim 100~{\rm MeV}$, naturally sets the energy density of DE and provides a compelling explanation for why its impact on cosmic expansion becomes significant only in the recent cosmological epoch. Importantly, the entire framework is grounded in Standard Model (SM) physics, involving no new fields or coupling constants. Key predictions of the model include: (a) A present-day DE equation of state $w_{\rm DE} > -1$, asymptotically approaching the de Sitter limit $w_{\rm DE} = -1$ in the future, with the corresponding asymptotic Hubble constant $\overline{H}$ set by $\Lambda_{\rm QCD}$. (b) At redshifts $z \ge 0$, $w_{\rm DE}$ can lie above or below $-1$ and may cross this boundary multiple times, behavior qualitatively consistent with recent DESI observations. (c) The solution to the Friedmann equation in this framework can deviate from the canonical $\Lambda$CDM form at $z \ge 0$. (d) If such deviation occurs, it can be tested using cosmological observations, CMB anisotropies, BAO, SNIa, and large-scale structure, which we propose to explore in future work. (e) Finally, we point out a potential connection between our framework and the observed $H_0$ tension.

Polarimetry has the capacity to provide a unique probe of the surface properties of asteroids. Trends in polarization behavior as a function of wavelength trace asteroid regolith mineral properties that are difficult to probe without measurements in situ or on returned samples. We present recent results from our ongoing survey of near-infrared polarimetric properties of asteroids. Our data reveal a mineralogical link between asteroids in the broader M- and K- spectral classes. In particular, M-type objects (16) Psyche, (55) Pandora, (135) Hertha, and (216) Kleopatra show the same polarimetric-phase behavior as K-type objects (89) Julia, (221) Eos, and (233) Asterope from visible through near-infrared light. The near-infrared behavior for these objects is distinct from other classes observed to date, and shows a good match to the polarimetric properties of M-type asteroid (21) Lutetia from the visible to the near-infrared. The best link for these objects from laboratory polarimetric phase curve measurements is to a troilite-rich fine-grained regolith. Our observations indicate that the M- and K-type spectral classes are most likely part of a continuum, with the observed spectral differences due to heterogeneity from partial differentiation, shock darkening of the surface material, or other later evolution of the original parent population. We also provide incidental J- and H-band polarimetric observations of other Main Belt asteroids obtained during our survey.

R. Arun, Blesson Mathew, B. Shridharan, K. Ujjwal, Akhil Krishna R, G. Maheswar, Sreeja S. Kartha

We present a comprehensive mid-infrared spectroscopic survey of 124 Herbig Ae/Be stars using newly processed Spitzer/IRS spectra from the newly released CASSISjuice database. Based on prominent dust and molecular signatures (polycyclic aromatic hydrocarbons, silicates, and hydrogenated amorphous carbons), we classify the stars into five groups. Our analysis reveals that 64% of the spectra show PAH emission, with detections peaking in the stellar effective temperature range 7000-11000 K (B9-A5). Silicate features appear in 50% of the sample and likewise diminish at higher temperatures. Additionally, we find that future PAH studies can focus on Herbig Ae/Be stars with a spectral index ($n_{2-24}$ > -1) and flared morphologies to maximize PAH detections. The 6.2 $\mu$m PAH band is the most frequently observed in our sample, shifting blueward with increasing stellar temperature, and this is the largest sample yet used to test that peak shift. The weaker 6.0 $\mu$m feature does not shift with 6.2 $\mu$m, implying a distinct origin of C=O (carbonyl) or olefinic C=C stretching relative to C--C vibrations. We examined the 11.0/11.2 $\mu$m PAH ratio using high-resolution Spitzer spectra for the first time in a sample of Herbig Ae/Be stars, finding a range of ionization conditions. This study provides a strong foundation for future JWST observations of intermediate-mass pre-main-sequence stars.

Purpose: We aim to investigate an accretion disk (AD) temperature profile of a doubly lensed quasar SBS 1520+530. Temperature profiles of accretion disks are known to be diagnostic of many of the physical properties of AGNs and quasars. Methods: Our approach involves application of the photometric reverberation mapping to the light curves of SBS 1520+530 obtained in the Johnson-Cousins V and R filters. The RM method implies that the time shift between the light curves taken in different spectral ranges determines the light travel time between the AD zones with different physical conditions. In determining the time shifts, we applied a method based on some useful properties of the orthogonal polynomials. Results: The variations in filter R lag those in filter V , with 1.25$\pm$ 0.63 days for the inter-band time shift averaged between the two image components and over the three seasons. The obtained time lag noticeably exceeds the value following for SBS 1520+530 from the classical model of optically thick geometrically thin AD. We considered two possible ways to resolve the discrepancy between the theory and observations. The first assumes an AD temperature profile flatter than the classical one. The second way is to consider an extended optically thick scattering envelope originated due to matter outflow from an AD interior. Conclusion: Both explanations may be consequences of the super-Eddington accretion in SBS 1520+530. Using bolometric luminosity estimates available for SBS 1520+530 from the literature, we obtained a value of $\approx$ 3.4 for the Eddington ratio, which indeed indicates a moderately super-Eddington accretion regime.

Gijs Vermariën, Thomas G. Bisbas, Serena Viti, Yue Zhao, Xuefei Tang, Rahul Ravichandran

Computational astrochemical models are essential for helping us interpret and understand the observations of different astrophysical environments. In the age of high-resolution telescopes such as JWST and ALMA, the substructure of many objects can be resolved, raising the need for astrochemical modeling at these smaller scales, meaning that the simulations of these objects need to include both the physics and chemistry to accurately model the observations. The computational cost of the simulations coupling both the three-dimensional hydrodynamics and chemistry is enormous, creating an opportunity for surrogate models that can effectively substitute the chemical solver. In this work we present surrogate models that can replace the original chemical code, namely Latent Augmented Neural Ordinary Differential Equations. We train these surrogate architectures on three datasets of increasing physical complexity, with the last dataset derived directly from a three-dimensional simulation of a molecular cloud using a Photodissociation Region (PDR) code, 3D-PDR. We show that these surrogate models can provide speedup and reproduce the original observable column density maps of the dataset. This enables the rapid inference of the chemistry (on the GPU), allowing for the faster statistical inference of observations or increasing the resolution in hydrodynamical simulations of astrophysical environments.

The dense cores of Milky Way globular clusters (GCs) play host to a variety of dynamical encounters between stellar objects, which can accelerate stars to velocities high enough to escape the GC. The most extreme examples of these encounters are interactions between single GC stars and binaries including at least one compact object. These interactions can result in ejection velocities of up to several hundred $\mathrm{km \ s^{-1}}$, approaching or even exceeding the escape velocity of the Galaxy itself. In order to study whether these interactions contribute to the Galactic population of hypervelocity stars (stars moving faster than the Galactic escape speed), we combine Monte Carlo $N$-body GC simulations, observations of Galactic GCs, and a particle spray code to generate realistic populations of stars which have escaped from Milky Way GCs following star + compact object binary (S+COB) interactions. We find that over the last 500 Myr, S+COB interactions have likely ejected $\sim$6300 stars from Galactic GCs, of which $839_{-67}^{+70}$ have present-day velocities exceeding $500 \; \mathrm{km \ s^{-1}}$. Using mock photometric observations, we find that $290_{-23}^{+28}$ ejected stars are detectable in Gaia Data Release 3, however, only $1_{-1}^{+2}$ stars faster than $500 \; \mathrm{km \ s^{-1}}$ are detectable. Even so, we show that observational prospects in the upcoming Legacy Survey of Space and Time are more optimistic, and future detected fast extratidal GC stars will serve as a useful probe of GC cores.

The Baryon Acoustic Oscillation (BAO) feature, imprinted in the transverse and radial clustering of dark matter tracers, enables the simultaneous measurement of the angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$ at a given redshift. Further, measuring the redshift space anisotropy (RSD) allows us to measure the combination $f_8\equiv f\sigma_8(z)$. Motivated by this, we simultaneously study the dynamics of background evolution and structure formation in an abstract phase space of dynamical quantities $x=H_0 D_A/c$, $p = dx/dz$, and $f_8$. We adopt a semi-cosmographic approach, whereby we do not pre-assume any specific dark energy model to integrate the dynamical system. The Luminosity distance is expanded as a Padé rational approximation in the variable $(1+z)^{1/2}$. The dynamical system is solved by using a semi-cosmographic equation of state, which incorporates the dark matter parameters along with the parameters of the Padé expansion. The semi-cosmographic $D_A(z), H(z)$ and $f\sigma_8(z)$, thus obtained, are fitted with BAO and RSD data from the SDSS IV. The reconstructed phase trajectories in the $3D$ $(x,p,f_8)$ space are used to reconstruct some diagnostics of background cosmology and structure formation. At low redshifts, a discernible departure from the $\Lambda$CDM model is observed. The geometry of the phase trajectories in the projected spaces allows us to identify three key redshifts where future observations may be directed for a better understanding of cosmic tensions and anomalies.

Suman Chatterjee, Shouvik Sarkar, Samir Choudhuri, Khandakar Md Asif Elahi, Somnath Bharadwaj, Shiv Sethi, Akash Kumar Patwa

Studying the diffuse Galactic synchrotron emission (hereafter, DGSE) at arc-minute angular scale is important to remove the foregrounds for the cosmological 21-cm observations. Statistical measurements of the large-scale DGSE can also be used to constrain the magnetic field and the cosmic ray electron density of our Galaxy's interstellar medium (ISM). Here, we have used the Murchison Widefield Array (MWA) drift scan observations at $154.2 \, {\rm MHz}$ to measure the angular power spectrum $({\cal C}_{\ell})$ of the DGSE of a region of the sky from right ascension (RA) $349^{\circ}$ to $70.3^{\circ}$ at the fixed declination $-26.7^{\circ}$. In this RA range, we have chosen 24 pointing centers (PCs), for which we have removed all the bright point sources above $\sim430 \, {\rm mJy}\,(3\sigma)$, and applied the Tapered Gridded Estimator (TGE) on residual data to estimate the ${\cal C}_{\ell}$. We use the angular multipole range $65 \le \ell \le 650$ to fit the data with a model, ${\cal C}^M_{\ell}=A\times \left(\frac{1000}{\ell}\right)^{\beta}+C$, where we interpret the model as the combination of a power law $(\propto \ell^{-\beta})$ nature of the DGSE and a constant part due to the Poisson fluctuations of the residual point sources. We are able to fit the model ${\cal C}^M_{\ell}$ for six PCs centered at $\alpha=352.5^{\circ}, 353^{\circ}, 357^{\circ}, 4.5^{\circ}, 4^{\circ}$ and $1^{\circ}$. We run the Markov Chain Monte Carlo (MCMC) ensemble sampler to get the best-fit values of the parameters $A, \beta$ and $C$ for these PCs. We see that the values of $A$ vary in the range $155$ to $400$ mK$^{2}$, whereas the $\beta$ varies in the range $0.9$ to $1.7$. We find that the value of $\beta$ is consistent at $2-\sigma$ level with the earlier measurement of the DGSE at similar frequency and angular scales.

The stability of the heliopause, the tangential discontinuity separating the solar wind from the interstellar medium, is influenced by various processes, including the Kelvin-Helmholtz instability. This study investigates the role of charge exchange collisions between protons and hydrogen (H) atoms in reduction of the Kelvin-Helmholtz growth rate at the heliopause. Using a two-dimensional gasdynamic model with the inclusion of H atoms, we perform numerical simulations of the plasma flow near the heliopause flanks. We conduct a parametric study by varying the Knudsen number. Our results indicate that charge exchange collisions play a crucial role in suppressing the Kelvin-Helmholtz instability. As the Knudsen number decreases, the flow transitions from an unstable regime to a smoother state.

Diffuse Interstellar Bands (DIBs) are crucial tracers of the interstellar medium (ISM), yet their carriers remain poorly understood. While large-scale surveys have advanced DIB studies in cool stellar spectra, measurements in hot stellar spectra are still limited. Using 287 277 high signal-to-noise (S/N $>$ 50) hot stellar spectra from the tenth data release of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope low-resolution spectroscopic survey (LAMOST LRS DR10), we systematically measured the three prominent optical DIBs at 5780, 5797, and 6614 Å. We published three catalogs containing 285 103, 279 195, and 281 146 valid measurements for the DIBs at 5780, 5797, and 6614 Å, respectively. Among them, 112 479, 25 232, and 71 048 are high-quality samples after rigorous quality control. To our knowledge, these are the largest hot-star DIB datasets in the northern sky. The catalogs provide spectral metadata, added astrometeric information, DIB profiles, and quality metrics. Our methodology and open-source pipeline ensure reproducibility, while the scale and precision of the data support future statistical studies. We anticipate that these catalogs will highlight the LAMOST's role in advancing DIB research and deepening our understanding of the ISM.

We have analyzed a sample of 2MIG isolated galaxies hosting AGNs (isolated AGNs) to assess whether their nuclear activity differs from that in a denser environment. The isolation criteria rule out interactions with other galaxies of similar size for at least $\sim$3 Gyr. We systemized the available Swift, NuSTAR, XMM-Newton, Chandra, and INTEGRAL X-ray data for isolated AGN at $z<0.05$, determining their general X-ray properties and peculiarities, spectral models, SMBH masses, and possible parameter correlations. We investigated the best spectral models of 25 isolated AGNs, including 10 recently fitted in this work. The intrinsic luminosities $\log L_{2-10,\mathrm{keV}}$ do not exceed $\sim43$. Our results indicate that the isolation of galaxies in the nearby Universe does not significantly affect nuclear activity. This is confirmed by the diversity of accretion types and the absence of any advantage of one particular basic/compound spectral model. We note an interesting case discovered for the ESO499-041 spectrum, where the data-to-model ratio shows significant changes after 6.4 keV. These deviations from the continuum may indicate the presence of a relativistic Iron line in the Chandra spectrum of ESO 499-041 compared to the Swift spectrum. Most isolated AGNs contain SMBH with $M_{\mathrm{SMBH}} \leq 10^7,M_{\odot}$. For the first time, we found a linear correlation between $L_{2-10,\mathrm{keV}}$ and $\log M_{\mathrm{SMBH}}$ that is not observed in other AGN samples. This trend may be affected by the limited sample size and requires further confirmation with larger datasets. Finally, we identified UGC10120, NG6300, and CGCG243-024 as promising candidates for Milky Way analogues.

The inspirals of supermaissive black hole binaries provide a convinced gravitational wave background in the nHz band, serving as the fiducial model of the recent gravitational wave signal reported by the PTA experiments. The uncertainties of the number of binaries contributing to each frequency bin introduce a foreground noise in the nHz and $\mu$Hz bands against the observation of the underlying gravitational wave backgrounds of the cosmological origin. In this work, we investigate a new method to constrain the cosmological gravitational wave strength under the astrophysical foreground. The energy density fluctuations from cosmological gravitational-wave sources can generally trigger the formation of compact subhalos of dark matter, and the upcoming Square Kilometer Array has the ability to constrain the abundance of the subhalos at the $\mathcal{O}(1)$ level. The cosmological gravitational wave energy spectra from various sources are expected to be constrained several orders of magnitude below the astrophysical foreground, providing more strict constraints on the parameter spaces of corresponding new physics models.

Very long baseline interferometry (VLBI) astrometry is used to determine the three-dimensional position and proper motion of astronomical objects. A typical VLBI astrometric campaign generally includes around ten observations, making it challenging to characterise systematic uncertainties. Our study on two bright pulsars, B0329+54 and B1133+16, involves analysis of broadband Very Long Baseline Array (VLBA) data over $\sim30$ epochs (spanning approximately $3.5\, {\rm years}$). This extended dataset has significantly improved the precision of the astrometric estimates of these pulsars. Our broadband study suggests that, as expected, the primary contribution to systematic uncertainties in L-band VLBI astrometry originates from the ionosphere. We have also assessed the effectiveness of the modified TEC (total electron content) mapping function, which converts vertical TEC to slant TEC, in correcting ionospheric dispersive delays using global TEC maps. The parallax and proper motion obtained from the multiple data sets, calibrated using the traditional and the modified TEC mapping functions, are consistent. However, the reduced chi-square values from least-squares fitting and precision of the fitted astrometric parameters show no significant improvement, and hence, the effectiveness of the new TEC mapping function on astrometry is unclear. For B0329+54, the refined parallax estimate is $0.611^{+0.013}_{-0.013}\, {\rm mas}$, with best-fit proper motion of $\mu_{\alpha} = 16.960^{+0.011}_{-0.010}\, {\rm mas\, yr^{-1}}$ in R.A. and and $\mu_{\delta} = -10.382^{+0.022}_{-0.022}\, {\rm mas\, yr^{-1}}$ in Dec. For B1133+16, the new estimated parallax is $2.705^{+0.009}_{-0.009}\, {\rm mas}$, with proper motions of $\mu_{\alpha} = -73.777^{+0.008}_{-0.008}\, {\rm mas\, yr^{-1}}$ and $\mu_{\delta} = 366.573^{+0.019}_{-0.019}\, {\rm mas\, yr^{-1}}$.

H. Bouy, G. Duchêne, G. Strampelli, J. Aguilar, J. Olivares, D. Barrado, S. N. Raymond, N. Huélamo, M. Tamura, E. Bertin, W. Brandner, J.-C. Cuillandre, P. A. B. Galli, N. Miret-Roig

Free-floating planetary mass objects-worlds that roam interstellar space untethered to a parent star-challenge conventional notions of planetary formation and migration, but also of star and brown dwarf formation. We focus on the multiplicity among free-floating planets. By virtue of their low binding energy (compared to other objects formed in these environments), these low-mass substellar binaries represent a most sensitive probe of the mechanisms at play during the star formation process. We use the Hubble Space Telescope and its Wide Field Camera 3 and the Very Large Telescope and its ERIS adaptive optics facility to search for visual companions among a sample of 77 objects members of the Upper Scorpius and Taurus young nearby associations with estimated masses in the range between approximately 5-33 M$_{\rm Jup}$. We report the discovery of one companion candidate around a Taurus member with a separation of 111.9$\pm$0.4 mas, or $\sim$18 au assuming a distance of 160pc, with an estimated primary mass in the range between 3-6 M$_{\rm Jup}$ and a secondary mass between 2.6-5.2 M$_{\rm Jup}$, depending on the assumed age. This corresponds to an overall binary fraction of 1.8$^{+2.6}_{-1.3}$% among free-floating planetary mass objects over the separation range $\ge$7 au. Despite the limitations of small-number statistics and variations in spatial resolution and sensitivity, our results, combined with previous high-spatial-resolution surveys, suggest a notable difference in the multiplicity properties of objects below $\sim$25 M$_{\rm Jup}$ between Upper Sco and Taurus. In Taurus, five companions were identified among 78 observed objects (4.9$^{+2.8}_{-2.0}$%), whereas none were found among 97 objects in Upper Sco ($\le$1.2%).}

Alexander C. Mayer, Thorsten Naab, Paola Caselli, Alexei V. Ivlev, Tommaso Grassi, Oliver Zier, Rüdiger Pakmor, Stefanie Walch, Volker Springel

We present simulations of the supernova-driven turbulent interstellar medium (ISM) in a simulation domain of volume $(256\,{\rm pc})^3$ within which we resolve the formation of protostellar accretion disks and their stellar cores to spatial scales of $\sim 10^{-4}$ au, using the moving-mesh code {\small AREPO}. We perform simulations with no magnetic fields, ideal magnetohydrodynamics (MHD) and ambipolar diffusion, and compare the resulting first Larson cores and their associated structures, including the accretion disks, their location within the larger-scale structure and the streamers connecting these. We find that disks of sizes $10-100\,{\rm au}$ form early in the simulations without magnetic fields, while there are no disks larger than 10 au with ideal MHD. Ambipolar diffusion causes large disks to form in a subset of cases (two out of six cores), and generally reduces the strength of outflows, which are seen to play a central role. When they are able to carry away significant angular momentum, they prevent the formation of a rotationally supported disk. Magnetic fields strengths grow from $0.1 - 1$ mG in the protostellar core to more than 10 G in the first Larson core in all simulations with ideal MHD. The rotationally supported disks which form can have rotation speeds $> 1$ km s$^{-1}$ even out to further than 100 au from the centre, become gravitationally unstable and form complex spiral substructures with Toomre $Q < 1$. We conclude that the impact of magnetic fields and non-ideal MHD on the formation of protostellar disks is substantial in realistic formation scenarios from the turbulent ISM.

Photometric data from the Transiting Exoplanet Survey Satellite (TESS) mission and radial velocities from the Gaia mission and ground-based observations were used to model the light curves and calculate the physical parameters of the eccentric eclipsing systems CH Ind, V577 Oph, CX Phe, and TIC 35481236. The components of these systems have temperatures between 6450 and 7500 K, masses between 1.4 and 1.85 solar masses, and radii between 1.49 and 3.05 solar radii. The residuals of these models were further analyzed using the Fourier method to reveal the pulsational frequencies of their oscillating components. Due to the similarity of the components of each system, the eclipses were used as spatial filters in order to determine which member is the pulsating star. CH Ind was found to pulsate in 46 frequencies; its primary component is a $\gamma$ Dor star and the secondary a $\delta$ Sct star. The primary component of V577 Oph oscillates in both the regimes of $\gamma$ Dor and $\delta$ Sct stars. Moreover, using past timings of minima, an eclipse timing variation analysis was also performed for V577 Oph, resulting in the calculation of the apsidal motion parameters and the existence of a third body around the system. Both components of CX Phe were found to be $\delta$ Sct stars; its primary has three independent frequencies in the range of 14.5-17.4 d$^{-1}$ and its secondary has two main modes of 5.19 and 7.22 d$^{-1}$. The analysis of TIC 35481236 indicates the hybrid $\delta$ Sct-$\gamma$ Dor nature of its secondary component. The physical and pulsational properties of the $\delta$ Sct stars of these systems were compared with those of other $\delta$ Sct stars-members of binaries in evolutionary diagrams.

A. J. Goodwin, R. Arcodia, G. Miniutti, J.C.A Miller-Jones, S. van Velzen

Quasi-periodic X-ray eruptions (QPEs) are a new class of repeating nuclear transient in which repeating X-ray flares are observed coming from the nuclei of generally low mass galaxies. Here we present a comprehensive summary of the radio properties of 12 bona-fide quasi-periodic eruption sources, including a mix of known tidal disruption events (TDEs) and AGN-like hosts. We include a combination of new dedicated radio observations and archival/previously published radio observations to compile a catalogue of radio observations of each source in the sample. We examine the overall radio properties of the sample and compare to the radio properties of known TDEs, given the apparent link between QPEs and TDEs. Overall we find compact, weak radio sources associated with 5/12 of the QPE sources and no signatures of strong AGN activity via a luminous radio jet. We find no radio variability or correlation between radio emission and the X-ray QPE properties, implying that the mechanism that produces the X-ray flares does not generate strong radio-emitting outflows. The compactness of the radio sources and lack of correlation between radio luminosity and SMBH mass is very unusual for AGN, but the radio spectra and luminosities are consistent with outflows produced by a recent TDE (or accretion event), in both the known TDE sources as well as the AGN-like sources in the sample.

Théo Lebeau, Stefano Ettori, Jenny G. Sorce, Nabila Aghanim, Jade Paste

Within the deep gravitational potential of galaxy clusters lies the intracluster medium (ICM). At the first order, it is considered to be at hydrostatic equilibrium within the potential well. However, evidence is growing that the ICM dynamics is non-negligible, is mostly turbulent in origin, and provides a non-thermal pressure support to the equilibrium. In this work, we intend to characterise the properties of the velocity field in the ICM of a simulated replica of the Virgo cluster. We first study the 3D and projected properties of the ICM velocity field by computing its probability density functions (PDFs) and its statistical moments. We then estimate the non-thermal pressure fraction from an effective turbulent Mach number, including the velocity dispersion. We finally compute the velocity structure function (VSF) from projected maps of the sightline velocity. We first show that the components of the 3D velocity field and the projected quantities along equivalent sightlines are anisotropic and affected by the accretion of gas from filaments. Then, we compare the mean statistical moments of the 3D velocity field to the mean properties of a hundred random projections. We show, in particular, an almost linear relation between the standard deviation estimated from direct simulation outputs and sightline velocity dispersion projections, comparable to the line broadening of X-ray atomic lines. However, this linear relation does not hold between the direct simulation outputs and the standard deviation of the sightline velocity projections, comparable to the line shift of X-ray atomic lines. We find a non-thermal pressure fraction around $6\%$ within $R_{500}$ and $9\%$ within $R_{vir}$ from sightline velocity dispersion, which is in good agreement with direct simulation outputs. Finally, we show that the VSF might probe the Active Galactic Nuclei (AGN) feedback turbulent injection scale.

Esin G. Gulbahar, Camille M. Diez, Aitor Ibarra, Ivan Valtchanov, Richard Saxton, Ignacio de la Calle Perez, Jose Lopez-Miralles, Alejandro Gonzalez Ganzabal, Peter Kretschmar

XMM-Newton is a cornerstone mission of the European Space Agency (ESA) for X-ray astronomy, providing high-quality X-ray data for astrophysical research since the start of the century. Its Science Analysis System (SAS) has been a reliable data reduction and analysis software, evolving throughout the years to meet changing user needs, while incorporating new methods. This paper presents the XMM-SAS Datalab, a tool within the cloud-based ESA Datalabs platform, designed to enhance the interactivity and collaborative potential of SAS. By integrating SAS with a modern, Python-based JupyterLab interface, it enables shared analysis workspaces, removes the need for local software setup, and provides faster access through containerised environments and preconfigured libraries. Moving SAS to the cloud preserves a consistent software setup while eliminating installation complexities, saving time and effort. A case study of the X-ray binary Vela X-1 demonstrates that the Datalabs platform reliably replicates local SAS outputs, with minimal deviations attributed to calibration file versions. The XMM-SAS Datalab allows straightforward X-ray data analysis with collaborative process, setting the way for future adaptations in e-science platforms and multi-wavelength astronomy, while offering traceability and reproducibility of scientific results.

Recently, Sodini et al. (2024) presented a sample of OI damped Lyman-$\alpha$ absorption system (DLA) analogs at $z\sim6$ that contain possible chemical signatures of Population III (Pop III) stars. In this paper, we use an N-body simulation-based semi-analytic model of the first stars and galaxies to predict the impact of Pop III stars on high-redshift DLAs. These Pop III DLA predictions are the first to include a number of important physical effects such as Lyman-Werner (LW) feedback, reionization, and external metal enrichment (all of which account for three-dimensional spatial fluctuations caused by halo clustering). We predict the abundance of DLAs as a function of their carbon-to-oxygen ratios ([C/O]). We find that our fiducial model is strongly ruled out by the data as it contains too few high-[C/O] DLAs, which have metals primarily from Pop III stars. However, increasing the delay time between Pop III and metal-enriched star formation due to supernovae feedback leads to better agreement with the data. Our results suggest that DLA analogs at $z\sim6$ are a promising probe of Pop III star formation for two key reasons. First, for reasonable parameter choices there are significant numbers of DLAs with metals primarily originating from Pop III stars. Second, we find that the number of DLAs with substantial Pop III contributions depends strongly on the Pop III star formation efficiency and the delay time between Pop III and metal-enriched star formation.

Gautam Nagaraj, Robin Ciardullo, Caryl Gronwall, Vandana Ramakrishnan, Kyoung-Soo Lee, Eric Gawiser, Nicole M. Firestone, Govind Ramgopal, J. Aguilar, Steven Ahlen, Davide Bianchi, David Brooks, Francisco Javier Castander, Todd Claybaugh, Andrei Cuceu, Axel de la Macorra, Arjun Dey, Biprateep Dey, Peter Doel, Jaime Forero-Romero, Enrique Gaztanaga, Satya Gontcho A Gontcho, Gaston Gutierrez, Hiram K. Herrera-Alcantar, Klaus Honscheid, Mustapha Ishak, Robert Kehoe, David Kirkby, T. Kisner, Anthony Kremin, Martin Landriau, Laurent Le Guillou, Michael Levi, Christophe Magneville, Marc Manera, Paul Martini, Aaron M. Meisner, Ramon Miquel, John Moustakas, Nathalie Palanque-Delabrouille, Francisco Prada, Ignasi Perez-Rafols, Graziano Rossi, Lado Samushia, Eusebio Sanchez, David J. Schlegel, Michael F. Schubnell, Hee-Jong Seo, Joseph H. Silber, David Sprayberry, Gregory Tarle, Francisco Valdes, Benjamin A. Weaver, Martin White, Rongpu Zhou, Hu Zou

The ubiquity and relative ease of discovery make $2\lesssim z\lesssim 5$ Ly$\alpha$ emitting galaxies (LAEs) ideal tracers for cosmology. In addition, because Ly$\alpha$ is a resonance line, but frequently observed at large equivalent width, it is potentially a probe of galaxy evolution. The LAE Ly$\alpha$ luminosity function (LF) is an essential measurement for making progress on both of these aspects. Although several studies have computed the LAE LF, very few have delved into how the function varies with environment. The large area and depth of the One-hundred-deg2 DECam Imaging in Narrowbands (ODIN) survey makes such measurements possible at the cosmic noon redshifts of z~2.4, ~3.1, and ~4.5. In this initial work, we present algorithms to rigorously compute the LAE LF and test our methods on the ~16,000 ODIN LAEs found in the extended COSMOS field. Using these limited samples, we find slight evidence that protocluster environments either suppress the numbers of very faint and very bright LAEs or enhance medium-bright LAEs in comparison to the field. We also find that the LF decreases in number density and evolves towards a steeper faint-end slope over cosmic time from z~4.5 to z~2.4.

Gael Buldgen, Gloria. Canocchi, Arthur. Le Saux, Vladimir A. Baturin, Regner Trampedach, Anna V. Oreshina, Sergey V. Ayukov, Anil Pradhan, Jean-Christophe Pain, Masanobu Kunitomo, Thierry Appourchaux, Rafael A. Garcia, Morgan Deal, Nicolas Grevesse, Arlette Noels, Joergen Christensen-Dalsgaard, Tristan Guillot, Devesh Nandal, Jérôme Bétrisey, Christophe Blancard, James Colgan, Philippe Cossé Christopher, J. Fontes, Ludovic Petitdemange, Charly Pincon

Helioseismology and solar modelling have enjoyed a golden era thanks to decades-long surveys from ground-based networks such as for example GONG, BiSON, IRIS and the SOHO and SDO space missions which have provided high-quality helioseismic observations that supplemented photometric, gravitational, size and shape, limb-darkening and spectroscopic constraints as well as measurements of neutrino fluxes. However, the success of solar models is also deeply rooted in progress in fundamental physics (equation of state of the solar plasma, high-quality atomic physics computations and opacities, description of convection and the role of macroscopic transport processes of angular momentum and chemicals, such as for example meridional circulation, internal gravity waves, shear-induced turbulence or even convection. In this paper, we briefly outline some key areas of research that deserve particular attention in solar modelling. We discuss the current uncertainties that need to be addressed, how these limit our predictions from solar models and their impact on stellar evolution in general. We outline potential strategies to mitigate them and how multidisciplinary approaches will be needed in the future to tackle them.

W. C. Jing, J. L. Han, C. Wang, P. F. Wang, T. Wang, N. N. Cai, J. Xu, Z. L. Yang, D. J. Zhou, Yi Yan, W. Q. Su, X. Y. Gao, L. Xie

The turbulent ionized interstellar medium diffracts radio waves and makes them propagate in multiple paths. The pulse-broadening observed at low frequencies results from the scattering effect of interstellar clouds of ionized gas. During the Galactic Plane Pulsar Snapshot (GPPS) survey and other projects by using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we detect the pulse-broadening for 122 pulsars in the radio frequency band between 1.0 and 1.5 GHz, including 60 newly discovered pulsars in the GPPS survey and 62 previously known pulsars. We find that a more accurate dispersion measure can be obtained from aligning the front edge of the scattered subband pulses at the 1/4 or 1/2 peak level for most pulsars with one dominant component in the intrinsic profile, and the best DM values from aligning the intrinsic profile components from the model-fitting. From the pulse profiles at a few subbands we derive the pulse-broadening timescale and the scattering spectral index. These scattering parameters are measured for the first time for 93 pulsars. For 29 pulsars with previously detected scattering features, our measurements of the pulse-broadening timescale are consistent with results in the literature. We find that pulsars behind spiral arms show a stronger scattering effect due to greater density fluctuations in the arm regions. With a properly derived dispersion measure and careful calibration, we also present polarization profiles for 41 pulsars in three subbands of FAST observations.

This paper deals with the effect of errors in the B and V magnitudes, or measurements in any other color system, on the width of the main sequence in a color-magnitude (Hertzsprung-Russell) diagram. The width is defined as the dispersion in apparent (or absolute) magnitude at a fixed, measured photometric color. I find that the dispersion is larger than might be thought, a priori. A statistical analysis is presented which demonstrates that the error in the magnitude residual from a linear approximation to the main sequence is Gaussian, but with a standard deviation which is much larger, in general, than the errors in the individual B and V magnitudes. This result is confirmed by a Monte Carlo simulation of a main sequence population with specified errors in B and V magnitudes, and can be explained on the basis of simple algebraic arguments.

Rahima Mokeddem, Bruno B. Bizarria, Jiajun Zhang, W.S. Hipólito-Ricaldi, Carlos Alexandre Wuensche, Elcio Abdalla, Filipe B. Abdalla, Amilcar R. Queiroz, Thyrso Villela, Bin Wang, Chang Feng, Edmar C. Gurjão, Alessandro Marins

We present a novel application of cosmological rescaling, or "remapping," to generate 21 cm intensity mapping mocks for different cosmologies. The remapping method allows for computationally efficient generation of N-body catalogs by rescaling existing simulations. In this work, we employ the remapping method to construct dark matter halo catalogs, starting from the Horizon Run 4 simulation with WMAP5 cosmology, and apply it to different target cosmologies, including WMAP7, Planck18 and Chevallier-Polarski-Linder (CPL) models. These catalogs are then used to simulate 21 cm intensity maps. We use the halo occupation distribution (HOD) method to populate halos with neutral hydrogen (HI) and derive 21 cm brightness temperature maps. Our results demonstrate the effectiveness of the remapping approach in generating cosmological simulations for large-scale structure studies, offering an alternative for testing observational data pipelines and performing cosmological parameter forecasts without the need for computationally expensive full N-body simulations. We also analyze the precision and limitations of the remapping, in light of the rescaling parameters $s$ and $s_m$, as well as the effects of the halo mass and box size thresholds.

The heating of a molecular cloud by photons emitted by a primary black hole (PBH) located inside the cloud is considered. For graphite and silicate dust particles, the dependence of dust temperature on the distance from PBH is derived, along with the emission spectrum of dust particles. The obtained spectrum is compared with the sensitivity of the Millimetron Space Observatory for various values of concentration and size of dust particles and different PBH masses.

L. Testi, A. Natta, S. Gozzi, C.F. Manara, J.P. Williams, R. Claes, U. Lebreuilly, P. Hennebelle, R. Klessen, S. Molinari

The value of the accretion luminosity during the early phases of star formation is a crucial information which helps us understand how stars form, yet it is still very difficult to obtain. We develop a new methodology to measure accretion luminosity using mid-infrared hydrogen recombination lines, and apply it to a limited sample of Class~I protostars in the Taurus and Ophiuchus star forming regions. We adopt the commonly used assumption that the properties of disk-protostar accretion in Class I objects is similar to the disk-star accretion in Class II objects. Using simultaneous observations of three hydrogen recombination lines Brg, Pfg, and Bra, we derive the mean intrinsic line ratios, and we verified that these are constant across the probed range of photospheric and accretion properties. We establish correlations between the line luminosities and accretion luminosity. We measure the extinction towards the line emission regions in Class I protostars comparing the observed line ratios to the Class II mean values. We then derive the Class I accretion luminosities from the established Class II correlations. We find that the accretion luminosity dominates the bolometric luminosity for the more embedded protostars, corresponding to lower values of the bolometric temperature. As the bolometric temperature increases above ~700K, there is a sharp drop of the contribution of the accretion from the bolometric luminosity. Our finding are in qualitative agreement with numerical simulations of star formation. We suggest that this methodology should be applied to larger and more statistically significant samples of Class I objects, for a more detailed comparison. Our results also suggest that by combining multiple infrared line ratios, it will be possible to derive a more detailed description of the dust extinction law in protostellar envelopes.

Thomas A. Baycroft, Alexandre Santerne, Amaury H.M.J. Triaud, Neda Heidari, Daniel Sebastian, Yasmin T. Davis, Alexandre C.M. Correia, Lalitha Sairam, Alix V. Freckelton, Aleyna Adamson, Isabelle Boisse, Gavin A.L. Coleman, Georgina Dransfield, João Faria, Salomé Grouffal, Nathan Hara, Guillaume Hébrard, Vedad Kunovac, David V. Martin, Pierre F.L. Maxted, Richard P. Nelson, Madison G. Scott, Owen J. Scutt, Matthew R. Standing

Planetary systems orbiting close binaries are valuable testing grounds for planet formation and migration models. More detections with good mass measurements are needed. We present a new planet discovered during the BEBOP survey for circumbinary exoplanets using radial velocities. We use data taken with the SOPHIE spectrograph at the Observatoire de Haute-Provence, and perform a spectroscopic analysis to obtain high precision radial velocities. This planet is the first radial velocity detection of a previously unknown circumbinary system. The planet has a mass of $0.56$ $M_{Jup}$ and orbits its host binary in 550 days with an eccentricity of 0.25. Compared to most of the previously known circumbinary planets, BEBOP-3b has a long period (relative to the binary) and a high eccentricity. There also is a candidate outer planet with a $\sim1400$ day orbital period. We test the stability of potential further candidate signals inside the orbit of BEBOP-3b, and demonstrate that there are stable orbital solutions for planets near the instability region which is where the Kepler circumbinary planets are located. We also use our data to obtain independent dynamical masses for the two stellar components of the eclipsing binary using High Resolution Cross-Correlation Spectroscopy (HRCCS), and compare those results to a more traditional approach, finding them compatible with one another.

Adam C. Schneider, Frederick J. Vrba, Justice Bruursema, Jeffrey A. Munn, Mike Irwin, Mike Read, Watson Varricatt, Tom Kerr, Klaus Hodapp, Simon Dye, Stephen J. Williams, Andrew T. Cenko, Trudy M. Tilleman, Marc A. Murison, Barry Rothberg, Scott Dahm, Bryan Dorland, Andy Lawrence, Kenneth C. Chambers

This paper describes the second and third data releases (DR2 and DR3, respectively) from the ongoing United Kingdom Infrared Telescope (UKIRT) Hemisphere Survey (UHS). DR2 is primarily comprised of the $K$-band portion of the UHS survey, and was released to the public on June 1, 2023. DR3 mainly includes the $H$-band portion of the survey, with a public release scheduled for September 2025. The $H$- and $K$-band data releases complement the previous $J$-band data release (DR1) from 2018. The survey covers approximately 12,700 square degrees between declinations of 0 degrees and $+$60 degrees and achieves median 5$\sigma$ point source sensitivities of 19.0 mag and 18.0 mag (Vega) for $H$ and $K$, respectively. The data releases include images and source catalogs which include $\sim$581 million $H$-band detections and $\sim$461 million $K$-band detections. DR2 and DR3 also include merged catalogs, created by combining $J$- and $K$-band detections (DR2) and $J$-, $H$-, and $K$-band detections (DR3). The DR2 merged catalog has a total of $\sim$513 million sources, while the DR3 merged catalog contains $\sim$560 million sources.

S. Belov, T. Parmenter, T. Arber, D. Kolotkov, F. Reale, T. Goffrey

Hot solar coronal loops, such as flaring loops, reach temperatures where the thermal transport becomes non-local. This occurs when the mean-free-path of electrons can no longer be assumed to be small. Using a modified version of the Lare2d code, we study the evolution of flare-heated coronal loops under three thermal transport models: classical Spitzer-Harm (SH), a flux-limited local model (FL), and the non-local Schurtz-Nicolai-Basquet (SNB) model. The SNB model is used extensively in laser-plasma studies. It has been benchmarked against accurate non-local Vlasov-Fokker-Planck models and proven to be the most accurate non-local model which can be applied on fluid time-scales. Analysis of the density-temperature evolution cycles near the loop apex reveals a distinct evolutionary path for the SNB model, with higher temperatures and lower densities than local models. During energy deposition, the SNB model produces a more localised and intense temperature peak at the apex due to heat flux suppression, which also reduces chromospheric evaporation and results in lower post-flare densities. EUV emission synthesis shows that the SNB model yields flare light curves with lower peak amplitudes and smoother decay phases. We also find that non-local transport affects equilibrium loop conditions, producing hotter and more rarefied apexes. These findings emphasise the need to account for non-local conduction in dynamic solar phenomena and highlight the potential of the SNB model for improving the realism of flare simulations. Flux-limited conduction models cannot reproduce the results of non-local transport covered by the SNB model.

Chengcai Shen, Xiaocan Li, Yuan-Kuen Ko, John C. Raymond, Fan Guo, Vanessa Polito, Viviane Pierrard

In astronomical environments, the high-temperature emission of plasma mainly depends on ion charge states, which requires accurate analysis of the ionization and recombination processes. For various phenomena involving energetic particles, the non-Maxwellian distributions of electrons exhibiting high-energy tails can significantly enhance the ionization process. Therefore, accurately computing ionization and recombination rates with non-Maxwellian electron distributions is essential for emission diagnostic analysis. In this work, we report two methods for fitting various non-Maxwellian distributions by using the Maxwellian decomposition strategy. For standard \{kappa} distributions, the calculated ionization and recombination rate coefficients show comparable accuracy to other public packages. We apply the above methods to two specific non-Maxwellian distribution scenarios: (I) accelerated electron distributions due to magnetic reconnection revealed in a combined MHD-particle simulation; (II) the high-energy truncated \{kappa} distribution predicted by the exospheric model of the solar wind. During the electron acceleration process, ionization rates of high-temperature iron ions increase significantly compared to their initial Maxwellian distribution, while the recombination rates may decrease due to the electron distribution changes in low-energy ranges. This can potentially lead to an overestimation of the plasma temperature when analyzing the Fe emission lines under the Maxwellian distribution assumption. For the truncated \{kappa} distribution in the solar wind, the ionization rates are lower than those for the standard \{kappa} distribution, while the recombination rates remain similar. This leads to an overestimation of plasma temperature when assuming a \{kappa} distribution.

Zachary K. Garvin (1), Anaïs Roussel (1), Luoth Chou (2), Marco E. Castillo (2 and 3), Xiang Li (2), William B. Brinckerhoff (2), Sarah Stewart Johnson (1) ((1) Georgetown University, Washington, D.C., USA, (2) NASA Goddard Space Flight Center, Greenbelt, MD, USA, (3) Aerodyne Industries, Cape Canaveral, FL, USA)

The ESA/NASA Rosalind Franklin rover, planned for launch in 2028, will carry the first laser desorption ionization mass spectrometer (LDI-MS) to Mars as part of the Mars Organic Molecule Analyzer (MOMA) instrument. MOMA will contribute to the astrobiology goals of the mission through the analysis of potential organic biosignatures. Due to the minimal availability of comparable equipment, laboratory analyses using similar techniques and instrumentation have been limited. In this study, we present a modified commercial benchtop LDI-MS designed to replicate MOMA functionality and to enable rapid testing of samples for MOMA validation experiments. We demonstrate that our instrument can detect organic standards in mineral matrices, with MS/MS enabling structural identification even in complex mixtures. Performance was additionally validated against an existing LDI-MS prototype through the comparison of spectra derived from natural samples from a Mars analog site in the Atacama Desert. Lastly, analysis of Mars analog synthetic mineral mixes highlights the capacity of the instrument to characterize both the mineralogical and organic signals in mission-relevant samples. This modified benchtop instrument will serve as a platform for collaborative research to prepare for MOMA operations, test LDI parameters, and generate pre-flight reference data in support of the mission science and astrobiology specific goals.

Eric J. Korpela, David P. Anderson, Jeff Cobb, Matt Lebofsky, Wei Liu, Dan Werthimer (Space Sciences Laboratory, University of California, Berkeley)

SETI@home is a radio Search for Extraterrestrial Intelligence (SETI) project, looking for technosignatures in data recorded at multiple observatories from 1998 to 2020. Most radio SETI projects analyze data using dedicated processing hardware. SETI@home uses a different approach: time-domain data is distributed over the Internet to $\gt 10^{5}$ volunteered home computers, which analyze it. The large amount of computing power this affords ($\sim 10^{15}$ floating-point operations per second (FPOP/s)) allows us to increase the sensitivity and generality of our search in three ways. We use coherent integration, a technique in which data is transformed so that the power of drifting signals is confined to a single discrete Fourier transform (DFT) bin. We perform this coherent search over 123 000 Doppler drift rates in the range ($\pm$100 Hz s$^{-1}$). Second, we search for a variety of signal types, such as pulsed signals and arbitrary repeated waveforms. The analysis uses a range of DFT sizes, with frequency resolutions ranging from 0.075 Hz to 1221 Hz. The front end of SETI@home produces a set of detections that exceed thresholds in power and goodness of fit. We accumulated $\sim 1.2\times 10^{10}$ such detections. The back end of SETI@home takes these detections, identifies and removes radio frequency interference (RFI), and looks for groups of detections that are consistent with extraterrestrial origin and that persist over long timescales. This paper describes the front end of SETI@home and provides parameters for the primary data source, the Arecibo Observatory; the back end and its results are described in a companion paper.

Enrique Pérez-Montero, Borja Pérez-Díaz, José M. Vílchez, Igor A. Zinchenko, Asier Castrillo, Marta Gavilán, Sandra Zamora, Ángeles I. Díaz

The derivation of sulphur chemical abundances in the gas-phase of star-forming galaxies is explored in this work, using the emission lines produced in these regions in the optical part of the spectrum and by means of photoionization models. We adapted the code HII-CHI-mistry to account for these abundances by implementing additional grids of models that assume a variable sulphur-to-oxygen abundance ratio, beyond the commonly assumed solar ratio. The addition of these models, and their use in a new iteration of the code allows us to use sulphur lines to precisely estimate the sulphur abundance, even in the absence of auroral lines. This approach aligns with the results from the direct method, and no additional assumptions about the ionization correction factor are needed, as the models directly predict the total sulphur abundance. We applied this new methodology to a large sample of star-forming regions from the MaNGA survey, and we explored the variation of the S/O ratio as a function of metallicity, making corrections for the significant contribution from diffuse ionized gas, which particularly affects the [SII] emission. Our results indicate no significant deviations from the solar S/O value in the range 8.0 < 12+log(O/H) < 8.7, where the bulk of the MaNGA sample stays, , but also with possible enhancements of sulphur production both at low and high metallicity regimes. The latter may be linked to the depletion of oxygen in the gas-phase due to its incorporation onto dust grains

David P. Anderson (1), Eric J. Korpela (1), Dan Werthimer (1), Jeff Cobb (1), Bruce Allen (2) ((1) Space Sciences Laboratory, University of California, Berkeley, (2) Max Planck Institut für Gravitationsphysik (Albert Einstein Institut), Hanover)

SETI@home is a radio Search for Extraterrestrial Intelligence (SETI) project that looks for technosignatures in data recorded at the Arecibo Observatory. The data were collected over a period of 14 years and cover almost the entire sky visible to the telescope. The first stage of data analysis found billions of detections: brief excesses of continuous or pulsed narrowband power. The second stage removed detections that were likely radio frequency interference (RFI), then identified and ranked signal candidates: groups of detections, possibly spread over the 14 years, that plausibly originate from a single cosmic source. We manually examined the top-ranking signal candidates and selected a few hundred. In the third and final stage we are reobserving the corresponding sky locations and frequency ranges using the Five-hundred-meter Aperture Spherical Telescope (FAST) radio telescope. This paper covers SETI@home's second stage of data analysis. We describe the algorithms used to remove RFI and to identify and rank signal candidates. To guide the development of these algorithms, we used artificial candidate birdies that model persistent ET signals with a range of power, bandwidth, and planetary motion parameters. This approach also allowed us to estimate the sensitivity of our detection system to these signals.

Eleanor M. Gallay, James R. A. Davenport, Steve Croft

We present an exploration of technosignature research that is possible using real-time alert brokers from surveys such as the Zwicky Transient Facility (ZTF) and the upcoming Legacy Survey of Space and Time (LSST). Nine alert brokers currently stream up to 1 million alerts each night from ZTF, and LSST is projected to increase this volume by an order of magnitude. While these brokers are primarily designed to facilitate real-time follow-up of explosive transients such as supernovae, they offer a unique platform to discover rare forms of variability from nearby stars in real time, which is crucial for follow-up and characterization. We evaluate the capability for both spatial and temporal searches for extraterrestrial intelligence (SETI) methods using the currently available brokers, and present examples of technosignature searches using ZTF alert and archival data. We have deployed optical SETI techniques, such as planetary transit zone geometries and the SETI Ellipsoid. We have also developed a search for novel high-amplitude stellar dippers, and present a workflow that integrates features available directly through the brokers, as well as post-processing steps that build on the existing capabilities. Though the SETI methods that alert brokers can execute are currently limited, we provide suggestions that may enhance future technosignature and anomaly searches in the era of the Vera C. Rubin Observatory.

Caroline Heneka, Florian Nieser, Ayodele Ore, Tilman Plehn, Daniel Schiller

For many fundamental physics applications, transformers, as the state of the art in learning complex correlations, benefit from pretraining on quasi-out-of-domain data. The obvious question is whether we can exploit Large Language Models, requiring proper out-of-domain transfer learning. We show how the Qwen2.5 LLM can be used to analyze and generate SKA data, specifically 3D maps of the cosmological large-scale structure for a large part of the observable Universe. We combine the LLM with connector networks and show, for cosmological parameter regression and lightcone generation, that this Lightcone LLM (L3M) with Qwen2.5 weights outperforms standard initialization and compares favorably with dedicated networks of matching size.

Ismael Ayuso, Mariam Bouhmadi-López, Che-Yu Chen, Xiao Yan Chew, Konstantinos Dialektopoulos, Yen Chin Ong

We explore the role of the affine connection in $f(Q)$ gravity, a modified theory where gravity is governed by non-metricity within the symmetric teleparallel framework. Although the connection is constrained to be flat and torsionless, it is not uniquely determined by the metric, allowing for multiple physically distinct formulations. We analyze three such connections compatible with a homogeneous and isotropic universe to show that they yield markedly different cosmological dynamics, even under the same functional form of $f(Q)$. Using both analytical and numerical methods, including a Born-Infeld type model of $f(Q)$, we demonstrate that specific connections can resolve cosmological singularities like the Big Bang and Big Rip, replacing them with smooth de Sitter phases. Others retain singularities but with notable modifications in their behavior. These findings highlight the physical relevance of connection choice in $f(Q)$ gravity and its potential to address fundamental cosmological questions.

Inspired by the well-studied $\Lambda_{\rm s}$CDM model, we propose and investigate a class of dynamical dark energy models that evolve from a negative cosmological constant, transitioning to a positive value at low redshifts. Specifically, we introduce a dark energy framework based on a generalised ladder-step function for the cosmological constant. Furthermore, we present two additional models in which the cosmological constant undergoes a smooth sign change at a specified redshift. We provide a detailed discussion of the construction and theoretical properties of these models, analysing their background cosmological evolution using a cosmographic approach and the statefinder hierarchy parameters. Our results are compared with those obtained for the standard $\Lambda$CDM and $\Lambda_{\rm s}$CDM models. We also perform a careful analysis of the types of singularities that may arise in these models due to a sign change in the cosmological constant. Notably, we show that a continuous sign change removes the sudden singularity present in the $\Lambda_{\rm s}$CDM model, replacing it with a milder $w$-singularity.

Third-generation ground-based gravitational wave detectors are expected to observe $\mathcal{O}(10^5)$ of overlapping signals per year from a multitude of astrophysical sources that will be computationally challenging to resolve individually. On the other hand, the stochastic background resulting from the entire population of sources encodes information about the underlying population, allowing for population parameter inference independent and complementary to that obtained with individually resolved events. Parameter estimation in this case is still computationally challenging, as computing the power spectrum involves sampling $\sim 10^5$ sources for each set of hyperparameters describing the binary population. In this work, we build on recently developed importance sampling techniques to compute the SGWB efficiently and train neural networks to interpolate the resulting background. We show that a multi-layer perceptron can encode the model information, allowing for significantly faster inference. We test the network assuming an observing setup with CE and ET sensitivities, where for the first time we include the intrinsic variance of the SGWB in the inference, as in this setup it presents a dominant source of measurement noise.

The idea that coherent oscillations of a scalar field, oscillating over a time period that is much shorter than the cosmological timescale, can exhibit cold dark matter (CDM) like behavior was previously established. In our work we first show that this equivalence between the oscillating scalar field model and the CDM sector is exact only in a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime in the absence of cosmological constant and any other possible matter components in the universe when the mass of the scalar field is very large compared to the Hubble parameter. Then we show how to generalize the equivalence between the coherently oscillating scalar field model and the CDM sector in a spatially curved universe with multiple matter components. Using our general method, we will show how a coherently oscillating scalar field model can represent the CDM sector in the presence of non-minimal coupling of the CDM sector with radiation. Our method is powerful enough to work out the dynamics of gravitational collapse in a closed FLRW spacetime where the coherently oscillating scalar field model represents the CDM sector. We have, for the first time, presented a consistent method which specifies how a coherently oscillating scalar field model, where the scalar field is ultralight, acts like the CDM sector in a multicomponent universe.

We study fluctuations of rotating viscous stars, using the causal relativistic hydrodynamics of Bemfica, Disconzi, Kovtun, and Noronha. We derive, in a slow-rotation approximation, a coupled system of equations describing the propagation of axial gravitational waves through the star, which couple to internal viscous modes. We show that rotating viscous stars amplify incoming low-frequency gravitational waves, a phenomenon which we argue to be universal. Superradiant amplification does not seem to trigger an instability for uniformly rotating stars, even if the object is compact enough to have light rings.

Debasish Borah, Nayan Das, Sahabub Jahedi, Dipankar Pradhan

We propose a multi-messenger frontier probe of non-thermal or freeze-in massive particle (FIMP) dark matter (DM) by considering an effective field theory (EFT) setup. Assuming leptophilic operators connecting DM with the standard model (SM) bath, we consider DM mass $m_{\rm DM}$ and the reheat temperature of the Universe $T_{\rm rh}$ in a regime which prevents DM-SM thermalisation. Low $T_{\rm rh}$ allows sizeable DM-SM interactions even for non-thermal DM allowing the latter to be probed at direct, indirect detection frontiers as well as future electron-positron and muon colliders. An extended reheating period governed by monomial inflaton potential after its slow-roll phase not only generates the required abundance of non-thermal DM via ultraviolet (UV) freeze-in but also brings the scale-invariant primordial gravitational waves (GW) within reach of near future experiments across a wide range of frequencies. While particle physics experiments can probe sub-GeV scale $T_{\rm rh}$ and FIMP with mass $m_{\rm DM} \lesssim O(10)$ GeV, future GW detectors are sensitive to a much wider parameter space.

The so-called memory-burden effect implies that evaporating Primordial Black Holes (PBHs) inevitably stabilize before complete decay. This stabilization opens a new mass window for PBH Dark Matter below $10^{15}\,$g. The transition to the memory-burdened phase is not instantaneous but unfolds over cosmological timescales, with some PBHs entering this phase in the present epoch. Additionally, a fraction of PBHs undergo mergers today, forming ''young'' semiclassical black holes that evaporate at unsuppressed rates. Both processes generate fluxes of stable astrophysical particles, which are constrained by current measurements of high-energy $\gamma$-rays and neutrinos. Moreover, the steep increase in energy injection at higher redshifts perturbs the ionization history of the Universe, leading to complementary bounds from observations of the CMB temperature and polarization anisotropies. We find that the reopened window enabled by the memory-burden effect is largely within reach of detection, both locally and across cosmological distances. We further describe how our findings restrict the values of the critical exponent characterizing the memory burden phenomenon.

The suppression of cosmological structure at small scales is a key signature of dark matter (DM) produced via freeze-in in the low-mass regime. We present a comprehensive analysis of its impact, incorporating recent constraints from Milky Way satellite counts, strong gravitational lensing with JWST data, and the Lyman-$\alpha$ forest. We adopt a general strategy to translate existing warm dark matter (WDM) bounds into lower mass limits for a broad class of DM candidates characterized by quasi-thermal phase space distributions. The benefits of this approach include computational efficiency and the ability to explore a wide range of models. We derive model-independent bounds for DM produced via two-body decays, scatterings, and three-body decays, and apply the framework to concrete scenarios such as the Higgs portal, sterile neutrinos, axion-like particles, and the dark photon portal. Results from specific models confirm the validity of the model-independent analysis.

Dark matter (DM) particles gravitationally captured by the Sun can accumulate in its core and subsequently annihilate, producing neutrino fluxes that may be detectable on Earth. The intensity of these fluxes is highly sensitive to the properties of the underlying DM model, especially when the DM candidate is a scalar particle originating from spontaneous or non-linear symmetry breaking mechanisms. In this work, we explore the potential of solar neutrino fluxes to distinguish between the Standard Model extended by a scalar singlet and the non-linear Higgs portal scenarios in the context of a future DM discovery. We compute the expected neutrino fluxes within the regions of parameter space consistent with both relic density and current direct detection limits. Our results show that the non-linear model predicts neutrino fluxes that are systematically larger than those of the linear case, typically by at least one order of magnitude, and up to six orders of magnitude for DM masses around 1 TeV. These findings suggest that solar neutrino observations could provide a valuable probe to discriminate between these competing dark matter frameworks.

While the third LIGO--Virgo gravitational-wave transient catalog includes 90 signals, it is believed that ${\cal O}(10^5)$ binary black holes merge somewhere in the Universe every year. Although these signals are too weak to be detected individually with current observatories, they combine to create a stochastic background, which is potentially detectable in the near future. LIGO--Virgo searches for the gravitational-wave background using cross-correlation have so far yielded upper limits. However, Smith \& Thrane (2017) showed that a vastly more sensitive ``coherent'' search can be carried out by incorporating information about the phase evolution of binary black hole signals. This improved sensitivity comes at a cost; the coherent method is computationally expensive and requires a far more detailed understanding of systematic errors than is required for the cross-correlation search. In this work, we demonstrate the coherent approach with realistic data, paving the way for a gravitational-wave background search with unprecedented sensitivity.

Black holes gradually settle into their static configuration by emitting gravitational waves, whose amplitude diminish over time according to a power-law decay at fixed spatial locations. We show that the nonlinear tails in the presence of a quadratic source, which have been recently found to potentially dominate over the linear ones, can be simply derived from the AdS$_2$$\times$S$^2$ spacetime perspective with their amplitudes being related to the Aretakis constants.

In a recent study, arXiv:2401.13043 found evidence for a 6% flux contribution from Jupiter to the total flux rate time series data from the BOREXINO solar neutrino experiment, specifically during the time intervals 2019-2021 and 2011-2013. The significance of this detection was estimated to be around $2\sigma$. We reanalyze the BOREXINO data and independently confirm the Jovian signal with the same amplitude and significance as that obtained in arXiv:2401.13043. However, using the same technique, we also find a spurious flux contribution from Venus and Saturn (at $\sim 2\sigma$ significance), whereas prima facie one should not expect any signal from any other planet. We then implement Bayesian model comparison to ascertain whether the BOREXINO data contain an additional contribution from Jupiter, Venus or Saturn. We find Bayes factors of less than five for an additional contribution from Jupiter, and less than or close to one for Venus and Saturn. This implies that the evidence for an additional contribution from Jupiter is very marginal.

High-energy neutrino observation from the Seyfert galaxy NGC 1068 offers new insights into the non-thermal processes of active galactic nuclei. Simultaneous gamma-rays emitted by such sources can possibly oscillate into axion-like particles (ALPs) when propagating through astrophysical magnetic fields, potentially modifying the observed spectrum. To probe for ALP-induced signals, a robust understanding of the emission processes at the source is necessary. In this work, we perform a dedicated multi-messenger analysis by modeling a jet in the innermost vicinity of the central supermassive black hole of NGC 1068. We model in particular the neutrino and gamma-ray emission originating in lepto-hadronic collisions between jet accelerated particles and background particles from the corona, reproducing both the Fermi-LAT and IceCube data. These source models serve as a baseline for ALP searches, and we derive limits on the ALP-photon coupling by marginalizing over motivated ranges of astrophysical parameters. We find $g_{a\gamma} \lesssim 7 \times 10^{-11}$GeV$^{-1}$ for $m_a \lesssim 10^{-9}$ eV. These limits may be weaker than existing constraints, but they demonstrate the potential of multi-messenger observations to probe new physics. We conclude by discussing how additional upcoming multi-messenger sources and improved observational precision can enhance ALP sensitivity.

Quadratic gravity theories emerge as the low-energy limit of many grand-unified and quantum-gravity theories. We report the first gravitational-wave ringdown constraints on such theories obtained purely from the quasinormal-mode spectra of rapidly-spinning black-hole remnants. We find no quadratic-gravity signatures from ringdown signals measured by the LIGO-Virgo-KAGRA detectors, thereby constraining the coupling length scale of axi-dilaton gravity below 34 km, dynamical Chern-Simons gravity below 49 km, and scalar Gauss-Bonnet gravity below 47 km.

Extracting the faint gravitational-wave background (GWB) signal from dominant detector noise and disentangling its %diverse astrophysical and cosmological components remain significant challenges for traditional methods like cross-correlation analysis. We propose a novel hybrid approach that combines deep learning with Bayesian inference to identify and characterize the GWB more rapidly than current techniques. Our method utilizes a custom-designed multi-scale multi-headed autoencoder (MSMHAutoencoder) architecture to separate GWB signals from detector noise, and subsequently Marcov Chain Monte Carlo parameter estimation to disentangle the GWB components. Using simulated data representative of the LIGO-Virgo-KAGRA network at design sensitivity, we show that our MSMHAutoencoder can detect with high confidence (log noise Bayes factor of 3) a GWB from binary black hole mergers with fractional energy density $\Omega_{\text{BBH}} \approx 10^{-9}$ at 25 Hz. In the presence of such an astrophysical GWB, we can simultaneously measure a cosmological component as faint as $\Omega_{\text{Cosmo}} \approx 1.3 \times 10^{-10}$ using 47.4 days of training data.