Abstract visibility
Abstract text size

Papers for Tuesday, Oct 14 2025

Papers with local authors

Lindsey A. Kwok, Chang Liu, Saurabh W. Jha, Stéphane Blondin, Conor Larison, Adam A. Miller, Mi Dai, Ryan J. Foley, Alexei V. Filippenko, Moira Andrews, Jennifer E. Andrews, Katie Auchettl, Carles Badenes, Thomas G. Brink, Kyle W. Davis, Andreas Flörs, Lluís Galbany, Estefania Padilla Gonzalez, D. Andrew Howell, Sahana Kumar, Réka Könyves-Tóth, Natalie LeBaron, Colin W. Macrie, Keiichi Maeda, Kate Maguire, Curtis McCully, Nicolas E. Meza-Retamal, Rüdiger Pakmor, Jeniveve Pearson, Anthony L. Piro, Abigail Polin, Nabeel Rehemtulla, César Rojas-Bravo, David J. Sand, Chita Sangkachan, Huei Sears, Mridweeka Singh, Bhagya M. Subrayan, Kirsty Taggart, Tea Temim, Jacco H. Terwel, Samaporn Tinyanont, József Vinkó, Xiaofeng Wang, J. Craig Wheeler, Yi Yang
0 votes
Paper 13 — arXiv:2510.09760
0 votes
Paper 13 — arXiv:2510.09760

We present optical + near-infrared (NIR) + mid-infrared (MIR) observations of the normal Type~Ia supernovae (SN Ia) 2022aaiq and 2024gy in the nebular phase, continuously spanning 0.35--28 microns. Medium-resolution JWST spectroscopy reveals novel narrow ($v_{\mathrm{FWHM}}<1500$ km s$^{-1}$) [Ni II] 1.94 and 6.64 micron cores in both events. The MIR [Ni II] 6.64 micron line exhibits a distinct narrow core atop a broader base, indicating a central enhancement of stable Ni. This structure points to high central densities consistent with a near-Chandrasekhar-mass ($M_{Ch}$) progenitor or a high-metallicity sub-$M_{Ch}$ progenitor. From detailed line-profile inversions of SN 2024gy, we derive emissivity profiles for stable iron-group elements (IGEs), radioactive material, and intermediate-mass elements (IMEs), revealing spatially distinct ejecta zones. The [Ni III] 7.35 micron line shows a shallow-to-steep slope transition - a "broken-slope" morphology - that matches predictions for delayed detonation explosions with separated deflagration and detonation ashes. We also reanalyze and compare to archival JWST spectra of SN 2021aefx and the subluminous SN 2022xkq. We estimate a stable $^{58}$Ni mass of $\sim0.1$ M$_\odot$ for SN 2024gy, consistent with delayed detonation models, and $\sim0.01$ M$_\odot$ for SN 2022xkq, favoring sub-$M_{Ch}$ scenarios. These results demonstrate that resolved line profiles, now accessible with JWST, provide powerful diagnostics of explosion geometry, central density, and progenitor mass in SN Ia.

Yingfeng Liu, Shijie Sun, Kaifeng Yu, Furen Deng, Shifan Zuo, Jixia Li, Yougang Wang, Fengquan Wu, Xuelei Chen
0 votes
Paper 39 — arXiv:2510.10067
0 votes
Paper 39 — arXiv:2510.10067

For an interferometric array, an image of the sky can be synthesized from interferometric visibilities, which are the cross-correlations of the received electric voltages of pairs of array elements. However, to search for transient targets such as the fast radio burst (FRB), it is more convenient to use the beam-forming technique, where the real-time voltage outputs of the array elements are used to generate data streams (beams) which are sensitive to a specific direction. This is usually achieved by a weighted sum of the array element voltages, with the complex weight adjusted so that all outputs have the same phase for that direction. Alternatively, beams can also be formed from the weighted sum of the short time averaged correlation (visibility) data. We shall call these two approaches the electric voltage beam forming (EBF) and cross-correlation beam forming (XBF), respectively. All beams formed with the EBF can also be formed by the XBF method, but the latter can also generate beams which can not be generated by the former. We discuss the properties of these two kinds of beams, and the amount of computation required in each case. For an array with large number of elements, the XBF would require much more computation resource, although this is partly compensated by the fact that it allows integration over time. We study the impact of cross-coupling between array elements on the beamforming, first using a toy model, then for the case of the Tianlai Cylinder Pathfinder Array. In both cases, we find that the impact of the cross-coupling on the beam profile is relatively small. The understanding gained in this study is helpful in designing and understanding the beam-forming FRB digital backend for compact arrays such as the Tianlai array.

Marcelo C. Vicentin, Pablo Araya-Araya, Laerte Sodré Jr., Michael A. Strauss
0 votes
Paper 64 — arXiv:2510.10735
0 votes
Paper 64 — arXiv:2510.10735

We present an algorithm designed to identify galaxy (proto)clusters in wide-area photometric surveys by first selecting their dominant galaxy-i.e., the Brightest Cluster Galaxy (BCG) or protoBCG-through the local stellar mass density traced by massive galaxies. We focus on its application to the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Wide Survey to detect candidates up to $\rm z \sim 2$. In this work, we apply the method to mock galaxy catalogs that replicate the observational constraints of the HSC-SSP Wide Survey. We derive functions that describe the probability of a massive galaxy being the dominant galaxy in a structure as a function of its stellar mass density contrast within a given redshift interval. We show that galaxies with probabilities greater than 50\% yield a sample of BCGs/protoBCGs with $\gtrsim 65\%$ purity, where most of the contamination arises from galaxies in massive groups below our cluster threshold. Using the same threshold, the resulting (proto)cluster sample achieves 80\% purity and 50\% completeness for halos with $M_{\rm{halo}} \geq 10^{14} \ M_{\odot}$, reaching nearly 100\% completeness for $M_{\rm{halo}} \geq 10^{14.5} \ M_{\odot}$. We also assign probabilistic membership to surrounding galaxies based on stellar mass and distance to the dominant galaxy, from which we define the cluster richness as the number of galaxies more likely to be true members than contaminants. This allows us to derive a halo mass-richness relation. In a companion paper, we apply the algorithm to the HSC-SSP data and compare our catalog with others based on different cluster-finding techniques and X-ray detections.

Marcelo C. Vicentin, Laerte Sodré Jr., Michael A. Strauss, Erik V. R. de Lima, Pablo Araya-Araya
0 votes
Paper 65 — arXiv:2510.10736
0 votes
Paper 65 — arXiv:2510.10736

We present a selection of candidates of clusters and protoclusters of galaxies identified in the photometric data of the HSC-SSP Wide Public Data Release 3 (PDR3), spanning the redshift range $\rm 0.1 \leq z \leq 2$. The selection method, detailed in Vicentin et al. (2025), involves detecting massive galaxies located in high-density regions of matter, identified as potential central dominant galaxies, i.e., (proto)BCGs. Probabilistic criteria based on proximity to the candidate central galaxy and the expected stellar mass of member galaxies are applied to identify likely members of each structure. We produced updated photometric redshift estimates using deep learning methods trained on a dataset combining spectroscopic redshifts from the HSC-SSP Wide PDR3, high-accuracy photometric redshifts from the COSMOS2020 catalog, and mid-infrared data from the unWISE catalog for matched sources. Our method achieves a predicted purity of $\sim 90\%$ in detecting (proto)clusters, with $\gtrsim 65\%$ correctly identifying the (proto)BCG. A total of 16,007 candidate (proto)clusters were identified over an effective area of $\rm \sim 850 \ deg^{2}$ within the HSC-SSP Wide footprint. Comparisons with other existing catalogs reveal a good level of consistency, while also highlighting that different methods yield complementary discoveries. We further compare richness and halo masses from our optical catalog with those from recent X-ray cluster catalogs (eROSITA and MCXC-II), finding a moderate positive correlation and a scatter of $\rm \sim 0.4$ dex. This catalog provides a valuable new set of targets for the Prime Focus Spectrograph (PFS) instrument.

G. Dransfield, M. Timmermans, D. Sebastian, B. V. Rackham, A. Burgasser, K. Barkaoui, A. H. M. J. Triaud, M. Gillon, J. M. Almenara, S. L. Casewell, K. A. Collins, A. Fukui, C. Jano-Munoz, S. Kanodia, N. Narita, E. Palle, M. G. Scott, A. Soubkiou, A. Stokholm, J. Audenaert, G. Á. Bakos, Y. Beletsky, Z. L. de Beurs, Z. Benkhaldoun, A. Burdanov, R. P. Butler, D. Caldwell, J. D. Crane, Y. T. Davis, B.O. Demory, E. Ducrot, Y. Gómez Maqueo Chew, M. Gachaoui, J. D. Hartman, M. J. Hooton, E. Jehin, S. Mercier, F. Murgas, C. Murray, P. P. Pedersen, F. J. Pozuelos, M. Rice, G. Ross, S. A. Shectman, E. Softich, M. Tala Pinto, A. M. Vanderburg, J. Villasenor, J. de Wit, S. Zúñiga-Fernández
0 votes
Paper 118 — arXiv:2510.11528
0 votes
Paper 118 — arXiv:2510.11528

Giant planets orbiting low-mass stars on short orbits present a conundrum, as in the most extreme cases their existence cannot be reconciled with current models of core accretion. Therefore, surveys dedicated to finding these rare planets have a key role to play by growing the sample to overcome small number statistics. In this work we present MANGOS, a programme dedicated to the search for giant objects (planets, brown dwarfs, and low-mass stars) orbiting M dwarfs. We report on the discovery of five new giant planets (TOI-3288 Ab, TOI-4666 b, TOI-5007 b, TOI-5292 Ab, TOI-5916 b) first detected by TESS, and confirmed using ground-based photometry and spectroscopy. The five planets have radii in the range 0.99-1.12 $\mathrm{R_{Jup}}$, masses between 0.49--1.69~$\mathrm{M_{Jup}}$, and orbital periods between 1.43 and 2.91 days. We reveal that TOI-3288 and TOI-5292 are wide binaries, and in the case of TOI-5292 we are able to characterise both stellar components. We demonstrate that the planets presented are suitable for further characterisation of their obliquities and atmospheres. We detect a small but significant eccentricity for TOI-5007 b, although for this to be more robust, more observations are needed to fully sample the orbit. Finally, we reveal a correlation between stellar metallicity and planet bulk density for giant planets orbiting low-mass stars.

M. S. Lundkvist, J. R. Larsen, Y. Li, M. L. Winther, T. R. Bedding, H. Kjeldsen, T. R. White, M. B. Nielsen, G. Buldgen, C. Guillaume, A. L. Stokholm, D. Huber, J. L. Rørsted, P. Mani, F. Grundahl
0 votes
Paper 119 — arXiv:2510.11532
0 votes
Paper 119 — arXiv:2510.11532

HD 140283 is a well-studied metal-poor subgiant and a Gaia benchmark star, often used for testing stellar models due to its proximity, brightness, and low metallicity ([Fe/H] = -2.3 dex). Here we present the first asteroseismic analysis of HD 140283, providing improved constraints on its fundamental properties. The star was observed by TESS in 20-second cadence during Sector 51. We extracted a custom light curve and performed a frequency analysis, revealing a rich spectrum of solar-like oscillations including mixed modes. These were combined with parameters from the literature to provide constraints on our model inference performed with BASTA. Using a dense grid of models, we find a mass of $0.75 \pm 0.01 \ \mathrm{M}_\odot$, a radius of $2.078 \substack{+0.012\\-0.011} \ \mathrm{R}_\odot$, and an age of $14.2 \pm 0.4$ Gyr, in agreement with the upper limit set by the age of the Universe within $1\sigma$. The observed frequency of maximum power, $\left(\nu_\mathrm{max}\right)_\mathrm{obs} = 611.3 \pm 7.4 \ \mu\mathrm{Hz}$, is significantly higher than predicted from standard scaling relations ($\left(\nu_\mathrm{max}\right)_\mathrm{mod} = 537.2 \substack{+2.9\\-1.8} \ \mu\mathrm{Hz}$), extending known deviations into the metal-poor regime. To our knowledge, the oscillations in HD 140283 have the highest $\nu_\mathrm{max}$ of any metal-poor star to date, which will help to advance our understanding of oscillations in metal-poor stars in general. The results demonstrate the value of asteroseismology for precise age determination in old halo stars and taking custom abundances and opacities into account during the modelling is probably important for further improving models of such stars. In addition, a detailed characterisation of metal-poor stars, such as HD 140283, will also help advance our understanding of Population III stars and their impact on future stellar generations.

C. Y. Tan, W. Cerny, A. B. Pace, J. A. Sharp, K. Overdeck, A. Drlica-Wagner, J. D. Simon, B. Mutlu-Pakdil, D. J. Sand, A. M. Senkevich, D. Erkal, P. S. Ferguson, F. Sobreira, K. R. Atzberger, J. L. Carlin, A. Chiti, D. Crnojević, A. P. Ji, L. C. Johnson, T. S. Li, G. Limberg, C. E. Martínez-Vázquez, G. E. Medina, V. M. Placco, A. H. Riley, E. J. Tollerud, A. K. Vivas, T. M. C. Abbott, M. Aguena, O. Alves, D. Bacon, S. Bocquet, D. Brooks, D. L. Burke, R. Camilleri, J. A. Carballo-Bello, A. Carnero Rosell, J. Carretero, T.-Y. Cheng, Y. Choi, L. N. da Costa, M. E. da Silva Pereira, T. M. Davis, J. De Vicente, S. Desai, P. Doel, S. Everett, B. Flaugher, J. Frieman, J. García-Bellido, D. Gruen, G. Gutierrez, S. R. Hinton, D. L. Hollowood, D. J. James, K. Kuehn, S. Lee, J. L. Marshall, J. Mena-Fernández, F. Menanteau, R. Miquel, J. Myles, M. Navabi, D. L. Nidever, R. L. C. Ogando, A. A. Plazas Malagón, A. Porredon, S. Samuroff, E. Sanchez, D. Sanchez Cid, I. Sevilla-Noarbe, M. Smith, E. Suchyta, M. E. C. Swanson, V. Vikram, A. R. Walker, A. Zenteno
0 votes
Paper 127 — arXiv:2510.11684
0 votes
Paper 127 — arXiv:2510.11684

We report the discovery of three Milky Way satellite candidates: Carina IV, Phoenix III, and DELVE 7, in the third data release of the DECam Local Volume Exploration survey (DELVE). The candidate systems were identified by cross-matching results from two independent search algorithms. All three are extremely faint systems composed of old, metal-poor stellar populations ($\tau \gtrsim 10$ Gyr, [Fe/H] $ \lesssim -1.4$). Carina IV ($M_V = -2.8;\ r_{1/2} = 40 {\rm pc}$) and Phoenix III ($M_V = -1.2;\ r_{1/2} = 19 {\rm pc}$) have half-light radii that are consistent with the known population of dwarf galaxies, while DELVE 7 ($M_V = 1.2;\ r_{1/2} = 2 {\rm pc}$) is very compact and seems more likely to be a star cluster, though its nature remains ambiguous without spectroscopic followup. The Gaia proper motions of stars in Carina IV ($M_* = 2250^{+1180}_{-830} {\rm M_\odot}$) indicate that it is unlikely to be associated with the LMC, while DECam CaHK photometry confirms that its member stars are metal-poor. Phoenix III ($M_* = 520^{+660}_{-290} {\rm M_\odot}$) is the faintest known satellite in the extreme outer stellar halo ($D_{\rm GC} > 100$ kpc), while DELVE 7 ($M_* = 60^{+120}_{-40} {\rm M_\odot}$) is the faintest known satellite with $D_{\rm GC} > 20$ kpc.

All other papers

The Fisher information matrix is used widely in astronomy (and presumably other fields) to forecast the precision of future experiments while they are still in the design phase. Although many sources describe the mathematics of the formalism, few sources offer simple examples to help the beginner. This pedagogical document works through a few simple examples to develop conceptual understanding of the applications.

The most massive galaxies in the Universe also host the largest supermassive black holes (SMBHs), with masses of $10^9 \: \mathrm{M_{\odot}}$ and above. During their hierarchical assembly, these galaxies have experienced only a few major mergers at low redshift, but have accreted many low-mass galaxies across cosmic time, possibly hosting intermediate mass black holes (IMBHs). If some of these IMBHs migrate to the galactic center, they may form compact subsystems around the central SMBH. We investigate the evolution of such subsystems, consisting of ten $10^5 \: \mathrm{M_{\odot}}$ IMBHs at three different concentrations around a $10^9 \: \mathrm{M_{\odot}}$ SMBH. We evolve these systems both in isolation and in the presence of a companion SMBH, using \texttt{MSTAR}, a regularized integration method including relativistic effects up to post-Newtonian order 3.5PN. Our analysis focuses on gravitational--wave--driven intermediate--mass--ratio inspirals (heavy IMRIs) and direct plunges. We show that perturbations from a secondary SMBH enhance the number of IMBH direct plunges by more than a factor of two, making them the dominant merger channel. These plunges and IMRIs with a central $10^9 \: \mathrm{M_{\odot}}$ SMBH will contribute to SMBH growth but will likely evade detection with future gravitational-wave interferometers and pulsar timing arrays (PTAs). However, for galaxies with lower--mass SMBHs ($M_\bullet \lesssim 10^8 \:\mathrm{M_{\odot}}$), heavy IMRIs will be detectable with the Laser Interferometer Space Antenna (LISA) and can provide direct observational constraints on the existence of IMBHs, while the more numerous plunges will still remain hidden.

Genevieve Schroeder (Cornell), Ben Margalit, Brian D. Metzger, Wen-fai Fong, Benjamin P. Gompertz, Kate D. Alexander, Edo Berger, Tanmoy Laskar, Gavin P. Lamb, Andrew Levan, Charles D. Kilpatrick, Jillian C. Rastinejad

In addition to a $\gamma$-ray burst (GRB), the merger of two neutron stars may produce a temporarily or indefinitely stable neutron star remnant with a strong magnetic field (a "magnetar"). As this magnetar remnant spins down, it can deposit its rotational energy into the surrounding kilonova ejecta, producing synchrotron emission that peaks in the radio bands $\sim$months-years after the merger ("boosted kilonova"). The nearby ($z=0.0763$) long-duration GRB 211211A, which has an apparent kilonova counterpart and likely neutron star merger progenitor, may have produced such a remnant. We observed the location of GRB 211211A at 6 GHz with the NSF's Karl G. Jansky Very Large Array (VLA) spanning $\approx 0.54$-$1.7~$years after the burst. We do not detect any radio emission, placing strong limits on the energy deposited into the ejecta by any remnant to $\lesssim 4.4 \times 10^{52}~{\rm erg}$. Due to the proximity of the event, we are also able to place limits on a kilonova afterglow that did not receive any additional energy deposition, though it is possible such emission will be suppressed until $\sim 4~{\rm years}$ after the burst, when the kilonova is expected to overtake the forward shock of the GRB. Future observations with the VLA and next-generation radio facilities will be able to further constrain the magnetar-boosted kilonova and kilonova afterglow scenarios, as well as directly constrain models in the scenario that GRB 211211A was instead produced by a collapsar.

Binaries in which a massive donor star undergoes an extended ($\gtrsim$ kyr) phase of stable mass transfer onto a black hole (BH) accretor offer a promising channel for creating LIGO gravitational wave sources. However, in many systems the mass transfer terminates prematurely in a dynamical instability at orbital periods of a few days, culminating in the BH plunging into the donor and potentially disrupting and accreting its helium core at highly super-Eddington rates. Combining a suite of binary evolution models with analytic estimates and population synthesis, we predict the population of luminous transients from delayed dynamical instability (DDI) and attribute them to the "luminous" class of fast blue optical transients (LFBOTs). The initial plunge of the BH into the partially stripped envelope typically ejects $\sim 10M_{\odot}$ of H/He-enriched material at speeds $\sim 10^{2}-10^{3}$ km s$^{-1}$, generating a compact circumstellar medium (CSM) of radius $\lesssim 1000R_{\odot}$ by the time the BH meets and tidally disrupts the HeC. Rapid BH accretion generates a highly aspherical wind-driven explosion into the environment, powering UV/optical emission via CSM interaction and X-ray reprocessing that rises over a few days to a luminosity $\sim 10^{44}-10^{45}$ erg s$^{-1}$ before fading as the disk spreads outwards and accretion rate drops. Luminous radio/sub-mm emission is generated over several months as the jet collides with the slow quasi-spherical binary outflow, generated by the stable mass transfer preceding DDI, extending to radii $\sim 10^{17}$ cm, in agreement with the inferred CSM environments of LFBOTs. We estimate local rates of DDI merger transients $5-300$ Gpc$^{-3}$ yr$^{-1}$, with a preference for low-metallicities, in agreement with LFBOT demographics. Taken together, our results support LFBOTs as being luminous signposts of "failed" gravitational wave sources.

We characterize the spatial power spectrum of density fluctuations in magnetohydrodynamic flows in a suite of high-resolution, long-time-span general relativistic magnetohydrodynamic (GRMHD) simulations. Extracting the local spatial power spectrum in curved spacetime directly from GRMHD simulations can be challenging for several conceptual and mechanical reasons, including choices of the reference frame, the non-uniform co-ordinate grid of the outputs and limited resolution. Taylor's frozen-in hypothesis describes a mapping between the temporal and spatial power spectrum of turbulence, which we apply to density fluctuations. We explore the validity of the assumptions underlying Taylor's hypothesis and evaluate its applicability in extracting spatial power spectra of density fluctuations of black hole accretion flows. Using outputs from the GRMHD code KORAL, we explore models with strong and ordered magnetic fields (MAD, Magnetically Arrested Disks) as well as weak and disordered magnetic fields (SANE, Standard and Normal Evolution). We explore the effects of black hole spin on the power spectra and characterize their spectral properties as a function of distance from the black hole. The observed power spectra follow a broken power law with two slopes separated by a break frequency. Our analysis shows a decrease in break frequency with increasing radius, with distinct trends between SANE and MAD flows. We also observe the first slope to be steeper for SANE flows and some notable distinctions between prograde and retrograde spins.

The application of deep machine learning methods in astronomy has exploded in the last decade, with new models showing remarkably improved performance on benchmark tasks. Not nearly enough attention is given to understanding the models' robustness, especially when the test data are systematically different from the training data, or "out of domain." Domain shift poses a significant challenge for simulation-based inference, where models are trained on simulated data but applied to real observational data. In this paper, we explore domain shift and test domain adaptation methods for a specific scientific case: simulation-based inference for estimating galaxy cluster masses from X-ray profiles. We build datasets to mimic simulation-based inference: a training set from the Magneticum simulation, a scatter-augmented training set to capture uncertainties in scaling relations, and a test set derived from the IllustrisTNG simulation. We demonstrate that the Test Set is out of domain in subtle ways that would be difficult to detect without careful analysis. We apply three deep learning methods: a standard neural network (NN), a neural network trained on the scatter-augmented input catalogs, and a Deep Reconstruction-Regression Network (DRRN), a semi-supervised deep model engineered to address domain shift. Although the NN improves results by 17% in the Training Data, it performs 40% worse on the out-of-domain Test Set. Surprisingly, the Scatter-Augmented Neural Network (SANN) performs similarly. While the DRRN is successful in mapping the training and Test Data onto the same latent space, it consistently underperforms compared to a straightforward Yx scaling relation. These results serve as a warning that simulation-based inference must be handled with extreme care, as subtle differences between training simulations and observational data can lead to unforeseen biases creeping into the results.

Recent measurements of the mean free path (MFP) of ionizing photons at $z=6$ find that it is significantly shorter than extrapolations from lower $z$. This has a substantial impact on the topology of reionization and thus the prospects of tomography of the 21-cm signal from upcoming radio interferometers. In this work we develop the first analytic model of reionization which explicitly incorporates the MFP as a free parameter, allowing us to transparently explore its effect on the process. Our model is based on the excursion set formalism with an ionization condition which accounts for absorptions parameterized through the MFP. With the goal of observational comparison, we include additional modifications which make our model particularly suitable for predicting one-point statistics of the ionization field (and 21-cm signal), which are among the fundamental quantities for tomography. We find that the effect of the MFP is much more significant during the later stages of reionization, and that including a shorter MFP reduces the size of HII regions by around an order of magnitude towards the end of reionization compared with analytic models which do not account for the MFP. We find that the reported MFP value produces a contrast in the 21-cm signal of $\mathcal{O}$(1 mK) or less at resolutions $\theta \sim $ 15--35 arcmin, an order of magnitude below naive estimates and up to a factor of several smaller than when using a larger MFP value extrapolated from low $z$, requiring significantly more sensitivity for imaging. We compare the contrast to noise estimates for arrays similar in size to HERA and SKA-Low and find that SKA has sufficient sensitivity for direct imaging (at the largest scales considered), while the predicted signal will be challenging for arrays similar in size to HERA. Our model indicates that more detailed sensitivity estimates are warranted in the context of a short MFP.

Jason A. S. Hunt, Michael S. Petersen, Martin D. Weinberg, Kathryn V. Johnston, Marcel Bernet, Kathryne J. Daniel, Sóley Ó. Hyman, Adrian M. Price-Whelan, Arpit Arora, the EXP Collaboration

The Milky Way is known to contain a stellar bar, as are a significant fraction of disc galaxies across the universe. Our understanding of bar evolution, both theoretically and through analysis of simulations indicates that bars both grow in amplitude and slow down over time through interaction and angular momentum exchange with the galaxy's dark matter halo. Understanding the physical mechanisms underlying this coupling requires modelling of the structural deformations to the potential that are mutually induced between components. In this work we use Basis Function Expansion (BFE) in combination with multichannel Singular Spectral Analysis (mSSA) as a non-parametric analysis tool to illustrate the coupling between the bar and the dark halo in a single high-resolution isolated barred disc galaxy simulation. We demonstrate the power of mSSA to extract and quantify explicitly coupled dynamical modes, determining growth rates, pattern speeds and phase lags for different stages of evolution of the stellar bar and the dark matter response. BFE & mSSA together grant us the ability to explore the importance and physical mechanisms of bar-halo coupling, and other dynamically coupled structures across a wide range of dynamical environments.

Dominika Ďurovčíková, Anna-Christina Eilers, Yuzo Ishikawa, Minghao Yue, Marianne Vestergaard, Frederick B. Davies, Jan-Torge Schindler, Xiaohui Fan, Fabrizio Arrigoni Battaia, Marta Volonteri, Robert A. Simcoe, Joseph F. Hennawi, Laura Blecha, Irham T. Andika, Sarah E. I. Bosman, Rebekka Bieri

Measurements of quasar lifetimes at high redshift indicate that the earliest billion-solar-mass supermassive black holes (SMBHs) have only been active as luminous quasars for less than a million years. Recently, extended Ly$\alpha$ nebulae around $z\sim6$ quasars have revealed that these short observed lifetimes are unlikely a sightline-dependent effect. However, the interpretation of Ly$\alpha$ emission is not straightforward due to its resonant nature. In this work, we use rest-frame optical emission lines, which more directly trace photoionization by the quasar, to unambiguously validate the short line-of-sight quasar lifetimes observed at early cosmic epochs. We use deep James Webb Space Telescope/NIRSpec IFU observations of five $z\sim 6$ quasars with small proximity zones to search for their extended emission line nebulae in H$\alpha$ and [O III]$5007$, and detect extended emission in both emission lines around four quasars in our sample. We then use the light-crossing time of these nebulae to measure quasar lifetimes along transverse sightlines. Using their H$\alpha$ nebulae, we also confirm that recombination is likely the dominant emission mechanism behind their previously detected Ly$\alpha$ nebulae. Our results confirm the existence of high-redshift quasars with extremely short lifetimes, $t_{\rm Q} \lesssim 10^{5}\ {\rm yr}$, hosting billion-solar-mass black holes, indicating that rapid accretion is likely responsible for the assembly of SMBHs in the early Universe.

Low-mass black holes hosted by dwarf galaxies offer valuable insights into galaxy formation and the growth of the massive black holes found in massive galaxies. Their detection as AGN is challenging due to their low luminosity and compact size. This can be circumvented employing multi-wavelength observational strategies, such as combining optical and radio observations, which enables the detection of AGN features that may be hidden in single-wavelength analyses We aim to detect any jet-like emission indicative of the presence of an AGN in a sample of four dwarf galaxies with AGN signatures based on spatially resolved emission line diagnostic diagrams with SDSS MaNGA. Confirming the presence of an AGN will prove IFU spectroscopy to be a resourceful tool for identifying hidden or switched-off AGN. Using VLA radio observations, we image the radio emission of the four dwarf galaxies and derive their integrated radio flux and luminosity. We compare these to that expected from star formation processes to determine the origin of the radio emission and probe if it is consistent with the results of the emission line diagnostic diagrams. We find that one out of the four galaxies shows AGN radio emission consistent with the analysis of the MaNGA IFU data. The kinetic jet power of this source is Qjet ~ 1e42 erg / s, indicating that dwarf galaxies can host radio jets as powerful as those of massive radio galaxies. This galaxy exhibits an AGN outflow able to escape the gravitational bound produced by the dark matter halo, along with a decrease in the star formation rate of the central region. This suggests the presence of negative feedback from the AGN, which could be suppressing star formation. The other three galaxies exhibit regions of radio emission consistent with a stellar origin and overlapping with the star-forming regions found in the IFU spectroscopy.

Many astrophysical simulations involve extreme dynamic range of timescales around 'special points' in the domain (e.g. black holes, stars, planets, disks, galaxies, shocks, mixing interfaces), where processes on small scales couple strongly to those on large scales. Adaptive resolution, multi-physics, and hybrid numerical methods have enabled tremendous progress on the spatial, physics, and numerical challenges involved. But often the limiter for following the long timescales of global evolution is the extremely short numerical timestep required in some subdomains (which leads to their dominating simulation costs). Recently several approaches have been developed for tackling this in problems where the short timescale solution is sampled and then projected as an effective subgrid model over longer timescales (e.g. 'zooming in and out'). We generalize these to a family of models where time evolution is modulated by a variable but continuous in space-and-time dilation/stretch factor $a({\bf x},\,t)$. This extends previous well-studied approaches (including reduced-speed-of-light and binary orbital dynamics methods), and ensures that the system comes to correct local steady-state solutions, and derive criteria that the dilation factor/timesteps/resolution must obey to ensure good behavior. We present a variety of generalizations to different physics or coupling scales. Compared to previous approaches, this method makes it possible to avoid imprinting arbitrary scales where there is no clear scale-separation, and couples well to Lagrangian or Eulerian methods. It is flexible and easily-implemented and we demonstrate its validity (and limitations) in test problems. We discuss the relationship between these methods and physical time dilation in GRMHD. We demonstrate how this can be used to obtain effective speedup factors exceeding $\gtrsim 10^{4}$ in multiphysics simulations.

Anna Ordog, Rebecca A. Booth, T.L. Landecker, Ettore Carretti, Alex S. Hill, Jo-Anne C. Brown, Artem Davydov, Leonardo Moutinho Caffarello, Luca B. Galler, Jonas Flygare, Jennifer L. West, A.G. Willis, Mehrnoosh Tahani, G.J. Hovey, Dustin Lagoy, Stephen Harrison, Mike Smith, Charl Baard, Rob H. Messing, D. A. Del Rizzo, Benoit Robert, Timothy Robishaw, John M. Dickey, George Morgan, Ian R. Kennedy, Marijke Haverkorn, Andrea Bracco, John Conway

Polarized synchrotron emission at meter to centimeter wavelengths provides an effective tracer of the Galactic magnetic field. Calculating Faraday depth, the most useful parameter for mapping the line-of-sight magnetic field, requires observations covering wide frequency bands with many channels. As part of the Global Magneto-Ionic Medium Survey we have observed polarized emission spanning 350-1030 MHz over the Northern sky, in the declination range ${-20^{\circ}}\leq{\delta}\leq{90^{\circ}}$. We used the 15 m telescope at the Dominion Radio Astrophysical Observatory (DRAO), equipped to receive orthogonal circular polarizations, with the Onsala Space Observatory band-1 feed developed for the Square Kilometre Array. Angular resolution varies across the band from $1.3^{\circ}$ to $3.6^{\circ}$. A digital spectrometer provided 42 kHz frequency resolution. Data were taken with the telescope moving rapidly in azimuth and are absolutely calibrated in intensity. Approximately 20% of the band was obscured by radio-frequency interference. Resolution in Faraday depth is $\sim6$ rad m$^{-2}$, and features as wide as $\sim38$ rad m$^{-2}$ are represented. The median sensitivity of the Faraday depth cube is 11 mK. Approximately 55% of sight-lines in this survey show Faraday complexity. This dataset, called ``DRAO GMIMS of the Northern Sky'' (DRAGONS), is the first to probe Faraday depth of the Northern sky in its frequency range and will support many scientific investigations. The data will be used to calibrate surveys with higher angular resolution, particularly Galactic foreground maps from the Canadian Hydrogen Intensity Mapping Experiment, and to provide information on large structures for aperture-synthesis telescopes, particularly the DRAO Synthesis Telescope. The data are available through the Canadian Astronomy Data Centre.

Mariia Marinichenko, Marcel P. van Daalen, Elena Sellentin, Jeger C. Broxterman, Matthieu Schaller, Joop Schaye

The scattering transform is a wavelet-based statistic capable of capturing non-Gaussian features in weak lensing (WL) convergence maps and has been proven to tighten cosmological parameter constraints by accessing information beyond two-point functions. However, its application in cosmological inference requires a clear understanding of its sensitivity to astrophysical systematics, the most significant of which are baryonic effects. These processes substantially modify the matter distribution on small to intermediate scales ($k\gtrsim 0.1\,h\,\mathrm{Mpc}^{-1}$), leaving scale-dependent imprints on the WL convergence field. We systematically examine the impact of baryonic feedback on scattering coefficients using full-sky WL convergence maps with Stage IV survey characteristics, generated from the FLAMINGO simulation suite. These simulations include a broad range of feedback models, calibrated to match the observed cluster gas fraction and galaxy stellar mass function, including systematically shifted variations, and incorporating either thermal or jet-mode AGN feedback. We characterise baryonic effects using a baryonic transfer function defined as the ratio of hydrodynamical to dark-matter-only scattering coefficients. While the coefficients themselves are sensitive to both cosmology and feedback, the transfer function remains largely insensitive to cosmology and shows a strong response to feedback, with suppression reaching up to $10\%$ on scales of $k\gtrsim 0.1\,h\,\mathrm{Mpc}^{-1}$. We also demonstrate that shape noise significantly diminishes the sensitivity of the scattering coefficients to baryonic effects, reducing the suppression from $\sim 2 - 10 \;\%$ to $\sim 1\;\%$, even with 1.5 arcmin Gaussian smoothing. This highlights the need for noise mitigation strategies and high-resolution data in future WL surveys.

This study explores the possibility of a time-varying dark energy (DE) equation of state (EoS) deviating from -1. We employ a comprehensive dataset of usual astronomical probes (Type Ia supernovae, baryon acoustic oscillations, Big Bang nucleosynthesis, Hubble data, and Planck 2018 CMB) alongside future mock gravitational wave (GW) distance measurements from the Einstein Telescope. We utilize the Pad'e approximation, a versatile framework encompassing well-known DE models like constant EoS, Chevallier-Polarski-Linder parametrization and other time-evolving DE parametrizations. Within Pad'e parametrization, we examine three specific forms (Pad'e-I, SPad'e-I, Pad'e-II) applied to both spatially flat and non-flat universes. Pad'e-II exhibits particularly interesting features in terms of the evidence of dynamical DE at many standard deviations. Our results can be summarized as follows. Flat Universe: When analyzing the combined dataset of standard probes (including CMB) with Pad'e-II in a flat universe, we find a strong preference (6.4{\sigma}) for a dynamical (time-varying) DE EoS. This preference remains significant (4.7{\sigma}) even when incorporating future GW data. Non-Flat Universe: In a non-flat universe, the combined standard datasets (without or with CMB) also indicate dynamical DE EoS at a high confidence level (6.2{\sigma} and 6.4{\sigma}, respectively). The addition of GW data slightly reduces the evidence (3.8{\sigma} and 5.1{\sigma}, respectively), but the preference persists. These results collectively suggest a robust case for dynamical DE in the dark sector. While a non-flat universe is not strongly favored, Pad'e-II hints at a possible closed universe when CMB data is included (with or without GW data).

J. Bätz, M. Mugrauer, K.-U. Michel, J. Reichert, A. Tschirschky, L. Pietsch, F. Edelmann, R. Neuhäuser

We present new radial velocity measurements of 13 selected intermediate mass stars (2 - 6 M$_\odot$). The measurements were performed between 29 April and 6 September 2024 at the University Observatory Jena using the échelle spectrograph FLECHAS. The radial velocity of eight stars was found to be constant during our spectroscopic monitoring, namely: 17 Dra A, HD 148374, HD 169487 A, 57 Cnc, $\gamma$ And, HD 11031, $\kappa$ And, and $\lambda$ Cas. In contrast, the radial velocity of five stars showed significant variability throughout or spectroscopic observation, namely: 7 CrB A, 7 CrB B, HD 214007, $\iota$ Her, and HD 201433 A. In all these cases, Keplerian orbital solutions were fitted to the observational data and the orbital elements of these spectroscopic binary systems were determined. In addition, we searched for wide companions of our targets using the third data release from ESA's Gaia mission, in order to determine the multiplicity status of these stars and contribute to the census of bright, nearby multiple stars.

Precise laser alignment in optical cavities is essential for high-precision laser interferometry. We report on a table-top optical experiment featuring two alignment sensing schemes: the conventional Wavefront Sensing (WFS) scheme which uses quadrant photodetectors (QPDs) to recover optical alignment, and the newly developed Radio Frequency Jitter Alignment Sensing (RFJAS) scheme, which uses an electro-optic beam deflector (EOBD) to apply fast angular modulation. This work evaluates the performance of RFJAS through a direct, side-by-side comparison with WFS. We present a detailed noise budget for both techniques, with particular emphasis on limitations at low frequencies, below 30\,Hz. Our results show that WFS performance is constrained by technical noise arising from beam spot motion (BSM), mainly due to beam miscentering on QPDs. In contrast, RFJAS is primarily limited by residual RF amplitude modulation. A blended scheme that combines both sensing methods may offer the most practical approach for use in gravitational wave detectors such as Advanced LIGO.

Elettra L. Piacentino, Aurelia Balkanski, Jenny Calahan, Anna Fitzsimmons, Mahesh Rajappan, Karin I. Oberg

Aromaticity is a common chemical functionalities in bioactive molecules. In interstellar and circumstellar environments benzene and other small aromatics are considered the precursor for more complex prebiotic molecules and they have shown to potentially have rich ice-phase photochemistry. The availability of small organic molecules in prebiotic networks depends on their photostability in astrophysical environments preceding planet formation, particularly during the protoplanetary disk stage, as the disk composition is linked to the chemical make-up of planets and planetesimals. We study the ultraviolet (UV) photodestruction (120-160 nm) of five aromatic molecules in undiluted ices and, for selected cases, in astrophysically relevant ice matrices (H2O, CO, CO2). For each ice, we measure the destruction cross sections as a function of photon exposure. In undiluted ices, aromatic molecules exhibit substantially lower photodestruction cross sections (sigma < 10-19 cm2) than aliphatic hydrocarbons, including cyclohexane, (sigma = 2.8-4x10-18 cm2). Furthermore, neither substituent nature nor size affects the aromatic stability in pure ices, suggesting that the strong intermolecular interactions among aromatic molecules provide protection against VUV exposure, even with small to mid-sized ring substituents. In mixed ices, the photodestruction and reactivity of aromatic molecules (sigma = 2.5-6.1x10-18 cm2) increases by more than an order of magnitude, but are still lower than in the gas-phase. We attribute this to a weaker cage effect and matrix-specific interactions. We use the experimental photodestruction cross sections to estimate the lifetime of aromatic molecules in protoplanetary disks, denileating the disks regions in which aromatic photochemistry is expected to be the most active.

Gracyn Jewett, Mukremin Kilic, Adam Moss, Alejandro H. Córsico, Matthew J. Green, Murat Uzundag, Pierre Bergeron, Warren R. Brown, Francisco C. De Gerónimo, Alberto Rebassa-Mansergas, Alex J. Brown, Vikram S. Dhillon, Stuart Littlefair

We present time-series photometry of 31 massive DA white dwarfs with $M\gtrsim 0.9~M_\odot$ within the ZZ Ceti instability strip from the Montreal White Dwarf Database 100 pc sample. The majority of the targets had no previous time-series photometry available, though several were classified as non-variable or potential pulsators in the literature. Out of the 31 candidates, we confirm 16 as pulsating. Our observations at three observatories have led us to discover the most massive pulsating white dwarf currently known, J0959$-$1828 ($M=1.32$ or $1.27~M_\odot$ for a CO versus ONe core), which is slightly more massive than the previous record holder J0049$-$2525. We study the sample properties of massive ZZ Ceti white dwarfs, and find several trends with their weighted mean periods. As predicted by theory, we see an increase in the weighted mean periods with decreasing effective temperature, and a decrease in pulsation amplitudes at the red edge of the instability strip. Furthermore, the weighted mean periods decrease with increasing stellar mass. Our observations show that the ZZ Ceti instability strip may not be pure at high masses. This is likely because the non-variable white dwarfs in the middle of the strip may be weakly magnetic, which could escape detection in the available low-resolution spectroscopy data, but may be sufficient to suppress pulsations. Extensive follow-up observations of the most massive white dwarfs in our sample have the potential to probe the interior structures and core-compositions of these white dwarfs with significantly crystallized cores.

Louis Desdoigts, Benjamin Pope, Max Charles, Peter Tuthill, Dori Blakely, Doug Johnstone, Shrishmoy Ray, Anand Sivaramakrishnan, Jens Kammerer, Deepashri Thatte, Rachel Cooper

The James Webb Space Telescope (JWST) hosts a non-redundant Aperture Masking Interferometer (AMI) in its Near Infrared Imager and Slitless Spectrograph (NIRISS) instrument, providing the only dedicated interferometric facility aboard - magnitudes more precise than any interferometric experiment previously flown. However, the performance of AMI (and other high resolution approaches such as kernel phase) in recovery of structure at high contrasts has not met design expectations. A major contributing factor has been the presence of uncorrected detector systematics, notably charge migration effects in the H2RG sensor, and insufficiently accurate mask metrology. Here we present Amigo, a data-driven calibration framework and analysis pipeline that forward-models the full JWST AMI system - including its optics, detector physics, and readout electronics - using an end-to-end differentiable architecture implemented in the Jax framework and in particular exploiting the dLux optical modelling package. Amigo directly models the generation of up-the-ramp detector reads, using an embedded neural sub-module to capture non-linear charge redistribution effects, enabling the optimal extraction of robust observables, for example kernel amplitudes and phases, while mitigating systematics such as the brighter-fatter effect. We demonstrate Amigo's capabilities by recovering the ABDor AC binary from commissioning data with high-precision astrometry, and detecting both HD206893B and the inner substellar companion HD206893c: a benchmark requiring contrasts approaching 10 magnitudes at separations of only 100 mas. These results exceed outcomes from all published pipelines, and re-establish AMI as a viable competitor for imaging at high contrast at the diffraction limit. Amigo is publicly available as open-source software community resource.

Romain Allart, Louis-Philippe Coulombe, Yann Carteret, Jared Splinter, Lisa Dang, Vincent Bourrier, David Lafrenière, Loïc Albert, Étienne Artigau, Björn Benneke, Nicolas B. Cowan, René Doyon, Vigneshwaran Krishnamurthy, Ray Jayawardhana, Doug Johnstone, Adam B. Langeveld, Michael R. Meyer, Stefan Pelletier, Caroline Piaulet-Ghorayeb, Michael Radica, Jake Taylor, Jake D. Turner

Atmospheric escape of planets on short orbital periods, driven by the host star's irradiation, influences their evolution, composition, and atmospheric dynamics. Our main avenue to probe atmospheric escape is through the near-infrared metastable helium triplet, which has enabled mass loss rate measurements for tens of exoplanets. Among them, only a few studies show evidence for out-of-transit absorption, supporting the presence of a hydrodynamic outflow. However, none of these observations precisely identified the physical extent of the outflow, either due to non-continuous or short-duration observations. This limits our measurements of accurate mass loss rates. Here we present the first continuous, full-orbit helium phase curve monitoring of an exoplanet, the ultra-hot Jupiter WASP-121b, obtained with JWST/NIRISS. It reveals helium absorption for nearly 60% of the orbit at >3sigma significance. Our results show that WASP-121b sustains a strong outflow, separating into two tails trailing and leading the planet. The persistent absorption from these tails, together with their measured radial velocity shifts, suggests that they remain in a collisional fluid regime at large distances from the planet and display very different dynamics. The leading trail has a higher density and moves toward the star, whereas the trailing trail is being pushed away from the star, with the latter being blue-shifted due to stellar irradiation pressure. While qualitatively agreeing with theoretical expectations, the observed structure of helium is not self-consistently reproducible by current models, limiting constraints on the mass loss rate. Furthermore, we show that while ground-based observations of the helium triplet are essential to measure the outflow dynamics precisely, they ideally should be combined with continuous JWST phase curves to constrain the absolute level of helium absorption.

Eleonora Parlanti, Giulia Tozzi Natascha M. Förster Schreiber, Claudia Pulsoni, Letizia Scaloni, Stavros Pastras, Pascal Oesch, Capucine Barfety, Francesco Belfiore, Jianhang Chen, Giovanni Cresci, Ric Davies, Frank Eisenhauer, Juan M. Espejo Salcedo, Reinhard Genzel, Rodrigo Herrera-Camus, Jean-Baptiste Jolly. Lilian L. Lee, Minju M. Lee, Daizhong Liu, Dieter Lutz, Filippo Mannucci, Giovanni Mazzolari, Thorsten Naab, Amit Nestor Shachar, Sedona H. Price, Alvio Renzini, Taro T. Shimizu, Amiel Sternberg, Martina Scialpi, Eckhard Sturm, Linda J. Tacconi, Hannah Übler, Stijn Wuyts

Cosmic noon represents the prime epoch of galaxy assembly, and a sweet spot for observations with the James Webb Telescope (JWST) and ground-based near-IR integral-field unit (IFU) spectrographs. This work analyses JWST NIRSpec Micro Shutter Array (MSA), NIRCam Wide Field Slitless Spectroscopy (WFSS) of K20-ID7, a large spiral, star-forming (SF) galaxy at z=2.2, with evidence for radial gas inflows. By exploiting the synergy with ground-based IFU ERIS observations, we conduct a comprehensive and resolved study of the interstellar medium (ISM) and stellar properties, from rest optical to near-IR, via emission-line diagnostics, resolved spectral energy distribution (SED) fitting of high-resolution imaging, and Pa$\beta$ line detection in NIRCam WFSS data. Our analysis reveals massive ($M_{\star}\simeq$(0.67-3.5)$\times$10$^{9}$ $M_{\odot}$) SF clumps with star formation rates (SFRs) ~3-24 $M_{\odot}$/yr, and quite low dust attenuation ($A_V\simeq$0.4), electron density ($n_{e}$<300 cm$^{-3}$), and ionisation (log(U)$\simeq -3.0$). The central bulge turns out to be modestly massive ($M_{\star}$=(7$\pm$3)$\times$10$^{9}$ M$_{\odot}$), heavily obscured ($A_V$=6.43$\pm$0.55), and likely to have formed most of its stellar mass in the past (SFR=82$\pm$42 $M_{\odot}$/yr over the last 100 Myr), yet still forming stars at a lower rate (SFR=12$\pm$8 M$_{\odot}$/yr over the last 10 Myr). We infer a metallicity 12+log(O/H)~8.54 and an apparent enhancement of the N/O abundance (log(N/O)$\simeq -1.0$) in all distinct galaxy regions, a likely consequence of dilution effects due to radial inflows of metal-poor gas. We measure a sub-solar sulfur abundance (log(S/O)$\simeq$-1.9). Finally, the radial stellar age profile reveals older stellar populations in the inner galaxy regions compared to the outskirts, pointing to an inside-out growth of K20-ID7.

Kevin Flaherty, Peter Knowlton, Tasan Smith-Gandy, A. Meredith Hughes, Marina Kounkel, Eric Jensen, James Muzerolle, Kevin Covey

Binary systems are a common site of planet formation, despite the destructive effects of the binary on the disk. While surveys of planet forming material have found diminished disk masses around medium separation ($\sim$10--100 au) binaries, less is known about tight ($<$10 au) binaries, where a significant circumbinary disk may escape the disruptive dynamical effects of the binary. We survey over 100 spectroscopic binaries in the Orion A region with ALMA, detecting significant continuum emission among 21 of them with disk masses ranging from 1--100 M$_{\oplus}$. We find evidence of systematically lower disk masses among the binary sample when compared to single star surveys, which may reflect a diminished planet forming potential around tight binaries. The infrared excess fraction among the binary sample is comparable to single stars, although the tight binaries without significant ALMA emission display tentative evidence of weaker 3-5$\mu$m excesses. The depletion of cold dust is difficult to explain by clearing alone, and the role of additional mechanisms needs to be explored. It may be the result of the formation pathway for these objects, systematic differences in intrinsic properties (e.g., opacity) or a bias in how the sample was constructed.

Steph Campbell, David J. Rosario, Houda Haidar, Enrique López Rodríguez, Dan Delaney, Erin Hicks, Ismael García-Bernete, Miguel Pereira-Santaella, Almudena Alonso Herrero, Anelise Audibert, Enrica Bellocchi, Donaji Esparza-Arredondo, Santiago García-Burillo, Omaira González Martín, Sebastian F. Hönig, Nancy A. Levenson, Chris Packham, Cristina Ramos Almeida, Dimitra Rigopoulou, Lulu Zhang

Broadband mid-infrared (MIR) imaging with high spatial resolution is useful to study extended dust structures in the circumnuclear regions of nearby AGN. However, broadband imaging filters cannot distinguish dust continuum emission from emission lines, and so accounting for the emission line contamination becomes crucial in studying extended dust in these this http URL paper uses Cycle 1 MIR imaging from JWST/MIRI and spectroscopy from JWST/MRS for 11 local Seyfert galaxies, as part of the Galactic Activity, Torus and Outflow Survey (GATOS). Three of the objects (NGC 3081, NGC 5728, and NGC 7172) exist in both datasets, allowing direct measurement of the line emission using the spectroscopy for these objects. We find that extended MIR emission persists on scales of 100s of parsecs after the removal of contamination from emission lines. Further, the line contamination levels vary greatly between objects (from 5% to 30% in the F1000W filter), and across filters, so cannot be generalised across a sample and must be carefully treated for each object and band. We also test methods to estimate the line contamination when only MRS spectroscopy or MIRI imaging is available, using pre-JWST ancillary data. We find that these methods estimate the contamination within 10 percentage points. This paper serves as a useful guide for methods to quantify and mitigate for emission line contamination in MIRI broadband imaging

Samantha N. Hasler, Leonid Pogorelyuk, Riley Fitzgerald, Kerri Cahoy, Rhonda Morgan

Planned and future missions, including the Habitable Worlds Observatory (HWO), will aim to directly image Earth-like exoplanets around Sun-like stars in reflected light. Determining whether an exoplanet is in the habitable zone of its star may be difficult in multi-planet systems when the observer does not know in advance which point source detection corresponds to which planet. This "confusion" problem will be a concern for future missions due to the high occurrence rate of multi-planet systems, and will be exacerbated by lack of prior knowledge about planets' orbital parameters or characteristics, particularly for systems at high inclination with respect to the observer. This work addresses the confusion problem by developing a photometry model and new orbit ranking function to augment the "deconfuser" tool and account for phase variation exhibited by a planet throughout its orbit. We demonstrate the new ranking scheme on a subset of thirty highly confused simulated multi-planet systems among three inclination groupings (low, medium, and high). Results indicate that photometry improves differentiation of previously confused orbits in 7/10 of the low inclination cases, 6/10 of the medium inclination cases, and 6/10 of the high inclination cases. This improvement in handling highly confused systems emphasizes that photometry shows promise for supporting orbit discrimination and deconfusion of directly imaged multi-planet systems, and should be considered when fitting orbits to detections.

Andrey Vayner, Tanio Díaz-Santos, Carl D. Ferkinhoff, Peter R. M. Eisenhardt, Daniel Stern, Lee Armus, Brandon S. Hensley, Daniel Anglés-Alcázar, Roberto J. Assef, Román Fernández Aranda, Andrew W. Blain, Hyunsung D. Jun, Norman W. Murray, Shelley Wright, Chao-Wei Tsai, Thomas Lai, Niranjan Chandra Roy, Drew Brisbin, Manuel Aravena, Jorge González-López, Guodong Li, Mai Liao, Devika Shobhana, Jingwen Wu, Dejene Zewdie

Galaxy-scale outflows are of critical importance for galaxy formation and evolution. Dust grains are the main sites for the formation of molecules needed for star formation but are also important for the acceleration of outflows that can remove the gas reservoir critical for stellar mass growth. Using the MIRI medium-resolution integral field spectrograph aboard the James Webb Space Telescope (JWST), we detect the 3.28 $\mu$m aromatic and the 3.4 $\mu$m aliphatic hydrocarbon dust features in absorption in a redshift 4.601 hot dust-obscured galaxy, blue-shifted by $\Delta$V=$-5250^{+276}_{-339}$ kms$^{-1}$ from the systemic redshift of the galaxy. The extremely high velocity of the dust indicates that the wind was accelerated by radiation pressure from the central quasar. These results pave a novel way for probing the physics of dusty outflows in active galaxies at early cosmic time.

This paper presents a systematic literature review focusing on the application of machine learning techniques for deriving observational constraints in cosmology. The goal is to evaluate and synthesize existing research to identify effective methodologies, highlight gaps, and propose future research directions. Our review identifies several key findings: (1) various machine learning techniques, including Bayesian neural networks, Gaussian processes, and deep learning models, have been applied to cosmological data analysis, improving parameter estimation and handling large datasets. However, models achieving significant computational speedups often exhibit worse confidence regions compared to traditional methods, emphasizing the need for future research to enhance both efficiency and measurement precision. (2) Traditional cosmological methods, such as those using Type Ia Supernovae, baryon acoustic oscillations, and cosmic microwave background data, remain fundamental, but most studies focus narrowly on specific datasets. We recommend broader dataset usage to fully validate alternative cosmological models. (3) The reviewed studies mainly address the $H_0$ tension, leaving other cosmological challenges-such as the cosmological constant problem, warm dark matter, phantom dark energy, and others-unexplored. (4) Hybrid methodologies combining machine learning with Markov chain Monte Carlo offer promising results, particularly when machine learning techniques are used to solve differential equations, such as Einstein Boltzmann solvers, as prior to Markov chain Monte Carlo models, accelerating computations while maintaining precision. (5) There is a significant need for standardized evaluation criteria and methodologies, as variability in training processes and experimental setups complicates result comparability and reproducibility (abridged).

Adalyn Gibson, Meredith A. MacGregor, Ward S. Howard, Ann Marie Cody, Mark Swain, Jennifer A. Burt, Laura Venuti, Evgenya Shkolnik, Neal J. Turner, Alan Didion, Jaime Nastal, David Makowski

We present the TESS discovery of only the second system of transiting exocomets with a sufficient number of events to measure the size distribution in the RZ Psc system, enabling comparisons with the $\beta$ Pictoris and Solar System size distributions. Twenty-four transits with absorption depths (AD) of 1--20\% were observed across three TESS sectors of the 20-50 Myr K0V star, detected as part of our TESS survey of extreme debris disks identified by their IR excess. We discover that the ADs (and hence exocomet radii) follow a broken power-law cumulative frequency distribution not previously seen in extrasolar contexts but similar to that observed in Solar System Kuiper Belt Object sizes, with power-law slopes above and below the break of $\gamma_\mathrm{AD>break}$=2.32$\pm$0.12 and $\gamma_\mathrm{AD<break}$=0.11$\pm$0.04, respectively. We derive size distributions of 1--7~km from two independent lines of evidence. We use the RZ Psc exocomet rate to predict exocomet yields for the Early eVolution Explorer (EVE) NASA astrophysics Small Explorer (SMEX) mission concept to obtain simultaneous photometry of 10$^4$ young stars in NUV, optical, and NIR bands. Assuming occurrence rates scaled from RZ Psc, EVE would detect 590 exocomets from $\approx$70 young systems in the optical band, with $\approx$120 simultaneous 5$\sigma$ detections in all three bands. These data would enable grain sizes of 200--700~nm and graphite--olivine compositions of dozens of events to be distinguished at 2.5--3$\sigma$, as well as a 4$\sigma$ determination of the accuracy of the Herschel-derived M-debris disk fraction.

Kai E. Yang, Xudong Sun, Lucas A. Tarr, Jiayi Liu, Peter Sadowski, S. Curt Dodds, Matthias Rempel, Sarah A. Jaeggli, Thomas A. Schad, Ian Cunnyngham, Yannik Glaser, Linnea Wolniewicz

Inferring the three-dimensional (3D) solar atmospheric structures from observations is a critical task for advancing our understanding of the magnetic fields and electric currents that drive solar activity. In this work, we introduce a novel, Physics-Informed Machine Learning method to reconstruct the 3D structure of the lower solar atmosphere based on the output of optical depth sampled spectropolarimetric inversions, wherein both the fully disambiguated vector magnetic fields and the geometric height associated with each optical depth are returned simultaneously. Traditional techniques typically resolve the 180-degree azimuthal ambiguity assuming a single layer, either ignoring the intrinsic non-planar physical geometry of constant optical-depth surfaces (e.g., the Wilson depression in sunspots), or correcting the effect as a post-processing step. In contrast, our approach simultaneously maps the optical depths to physical heights, and enforces the divergence-free condition for magnetic fields fully in 3D. Tests on magnetohydrodynamic simulations of quiet Sun, plage, and a sunspot demonstrate that our method reliably recovers the horizontal magnetic field orientation in locations with appreciable magnetic field strength. By coupling the resolutions of the azimuthal ambiguity and the geometric heights problems, we achieve a self-consistent reconstruction of the 3D vector magnetic fields and, by extension, the electric current density and Lorentz force. This physics-constrained, label-free training paradigm is a generalizable, physics-anchored framework that extends across solar magnetic environments while improving the understanding of various solar puzzles.

Charlotte Smith-Perez, Aidan Hembruff, Els Peeters, Alexander G. G. M. Tielens, Alessandra Ricca

Polycyclic Aromatic Hydrocarbons (PAHs) constitute a significant fraction of the Universe's carbon budget, playing a key role in the cosmic carbon cycle and dominating the mid-infrared spectra of astrophysical environments in which they reside. Although PAHs are known to form in the circumstellar envelopes of post-AGB stars, their formation and evolution are still not well-understood. We aim to understand how pristine complex hydrocarbons and PAHs in circumstellar environments transition to the PAHs observed in the ISM. The mid-infrared PAH spectra (5-18 micron) of the planetary nebula, NGC 7027, are investigated using spectral cubes from JWST MIRI-MRS. We report the first detection of spatially-resolved variations of the PAH spectral profiles across class A, AB, and B in all major PAH bands (6.2, 7.7, 8.6, and 11.2 micron) within a single source, NGC 7027. These variations are linked to morphological structures within NGC 7027. Clear correlations are revealed between the 6.2, 7.7, and 8.6 micron features, where the red components (6.26, 7.8, 8.65 micron) exhibit a strong correlation and the same is found for the blue components of the 6.2 and 7.7 (6.205 and 7.6 micron). The blue component of the 8.6 (8.56 micron) appears to be independent of the other components. We link this behaviour to differences in molecular structure of their PAH subpopulations. Decomposition of the 11.2 micron band confirms two previously identified components, with the broader 11.25 micron component attributed to emission from very small grains of PAH clusters rather than PAH emission. We show that PAH profile classes generally vary with proximity to the central star's UV radiation field, suggesting class B PAHs represent more processed species while class A PAHs remain relatively pristine, challenging current notions on the spectral evolution of PAHs.

Shi-Min Liang, Jian-Fu Zhang, Na-Na Gao, Nian-Yu Yi

Magnetic reconnection, often accompanied by turbulence interaction, is a ubiquitous phenomenon in astrophysical environments. However, the current understanding of the nature of turbulent magnetic reconnection remains insufficient. We investigate the statistical properties of reconnection turbulence in the framework of the self-driven reconnection. Using the open-source software package AMUN, we first perform numerical simulations of turbulent magnetic reconnection. We then obtain the statistical results of reconnection turbulence by traditional statistical methods such as the power spectrum and structure function. Our numerical results demonstrate: (1) the velocity spectrum of reconnection turbulence follows the classical Kolmogorov type of $E\propto k^{-5/3}$, while the magnetic field spectrum is steeper than the Kolmogorov spectrum, which are independent of limited resistivity, guide field, and isothermal or adiabatic fluid states; (2) most of the simulations show the anisotropy cascade, except that the presence of a guide field leads to an isotropic cascade; (3) reconnection turbulence is incompressible in the adiabatic state, with energy distribution dominated by the velocity solenoidal component; (4) different from pure magnetohydrodynamic (MHD) turbulence, the intermittency of the velocity field is stronger than that of the magnetic field in reconnection turbulence. The steep magnetic field spectrum, together with the velocity spectrum of Kolmogorov type, can characterize the feature of the reconnection turbulence. In the case of the presence of the guide field, the isotropy of the reconnection turbulence cascade is also different from the cascade mode of pure MHD turbulence. Our experimental results provide new insights into the properties of reconnection turbulence, which will contribute to advancing the self-driven reconnection theory.

Zi-Yi Zhou, Long Ji, Ling-Da Kong, Sergey S. Tsygankov, Qing-Cang Shui, Lian Tao

We report the detection of mHz quasi-periodic oscillations (QPOs) in four NuSTAR observations of 4U 1626-67 during its recent spin-down episode. By using a novel method based on the Hilbert-Huang Transform (HHT), we present the first QPO-phase-resolved timing and spectral analysis of accreting X-ray pulsars in low mass X-ray binaries. Broadband QPO waveforms have been reconstructed and exhibit approximately sinusoidal shapes, with fractional amplitudes that vary with energy. In addition, we find that spin pulse profiles exhibit stable shapes between different QPO phases with different instantaneous fluxes, while the fractional root-mean-square (rms) is distinct for different observations. In this source, both QPO-phase-resolved and averaged spectra can be modeled with a negative and positive powerlaws exponential (NPEX) model, and their spectral evolutions show a similar trend, suggesting that the QPO modulation is caused by accretion rate variability instead of a geometric obscuration. These results provide new constraints on accretion physics in strongly magnetized neutron stars and the underlying mechanisms of QPOs.

Type Ia supernovae (SNe Ia) are essential tools for addressing key cosmic questions, including the Hubble tension and the nature of dark energy. Modern surveys are predominantly photometry-based, making the construction of a clean photometric SNe Ia sample crucial. In this study, we investigate whether functional principal component analysis (FPCA) scores derived from photometric light curves, combined with ensemble learning, can reliably distinguish SNe Ia from other transients using the PLAsTiCC dataset. FPCA provides a data-driven, flexible characterization of light curves without relying on rigid theoretical model assumptions. Light curves are fitted by minimizing residuals with penalty terms from clean samples, making the method robust to outliers or poorly sampled bands. The first two FPCA scores and peak magnitudes across the five LSST bands are used as classification features. We implement two complementary binary classifiers: an ensemble boosting model (CatBoost) and a statistical probabilistic method based on Euclidean distances. CatBoost slightly outperforms the statistical method, achieving 98.5% accuracy and 97.8% precision. Performance remains robust (>90%) under typical photometric redshift uncertainties ({\sigma} = 0.1). On the spectroscopic DES Y5 sample, both methods reach approximately 90% accuracy and 95% precision, demonstrating strong out-of-domain generalization compared to state-of-the-art methods with limited cross-survey applicability. Applied to DECam DDF and DESIRT transients, the predictions strongly agree, and their intersection provides a high-confidence SNe Ia sample for cosmological analyses. Overall, this FPCA-based framework offers a powerful, flexible tool for classifying transients in upcoming large-scale surveys such as LSST and Roman.

David P. Huenemoerder (MIT), Sean J. Gunderson (MIT), Richard Ignace (East Tennessee State University), Joy S. Nichols (Harvard and Smithsonian Center for Astrophysics), A. M. T. Pollock (University of Sheffield), Pragati Pradhan (Embry-Riddle Aeronutical University), Norbert S. Schulz (MIT), Dustin K. Swarm (University of Iowa), José M. Torrejón (Universidad de Alicante)

We present 197 ks HETG and 95 ks NuSTAR spectra of the $\gamma\,$Cas-type object V750 Ara. The high-resolution X-ray spectra show that the target is similar to other objects of this class. Data are interpreted under the assumption that the X-rays come from an accreting white dwarf, and our analysis implies an accretion rate of about $3\times10^{-11}M_\odot\mathrm{yr}^{-1}$. Emission lines are weak, and predominantly from hydrogen-like ions: Mg XII, Si XIV, and S XVI. H-like and He-like Fe are both present, but Fe K$\alpha$ fluorescence is weak, being significantly detected only in the NuSTAR spectrum, but was not obviously detected in the HETG dispersed or zeroth-order spectra. The flux was variable above a level expected by Poisson statistics. There were no significant changes in the spectral hardness, though we are limited by lack of soft signal below 1 keV. Emission lines of Mg and Si were strong enough to measure velocity offsets and widths which were found to be marginally inconsistent. The H-like Mg line is consistent with instrumental broadening only, but shows a 300 km/s blueshift. He-like Mg and H-like Si lines have no significant shift in velocity but are broadened by about 1000 km/s. This suggests either different physical origins or velocity structure differing with plasma temperature.

GRS 1915+105 has been well studied since its discovery, and is well-known for its complex light curve variability. Using the full currently available Insight-HXMT dataset from July 2017 to June 2023, we make a comprehensive spectral-timing analysis of this source and report four main findings. First, we uncover a QPO frequency rising branch between MJD 58206 and 58230, where the centroid frequency increases from $\sim$2 Hz to $\sim$6 Hz, consistent with a spectral state transition from the hard to intermediate state. This rising branch completes the full QPO frequency evolution cycle when combined with the subsequent frequency decay phase, and had been missed in prior NICER and Insight-HXMT studies. Second, we identify a previously unreported Flare 3 during the obscured state, which shows distinct spectral and timing properties compared to the earlier flares. Third, we detect sub-Hz QPOs (<1 Hz) in all three flares, specifically at $\sim$0.01 Hz in Flare 1 and $\sim$0.2 Hz in both Flares 2 and 3. In particular, the weak $\sim$0.2 Hz signals observed in Flare 3 indicate ongoing coronal activity despite strong obscuration. Finally, a comparison between QPOs above and below 1 Hz suggests distinct origins, with the former likely arising from Lense-Thirring precession of the inner hot flow and the latter from magnetic perturbations driving a failed disk wind. These findings offer new insights into the unique accretion geometry and variability behaviors of GRS 1915+105.

L. Blot (1 and 2), K. Tanidis (3), G. Cañas-Herrera (4 and 5), P. Carrilho (6 and 7), M. Bonici (8 and 9), S. Camera (10 and 11 and 12), V. F. Cardone (13 and 14), S. Casas (15 and 16), S. Davini (17), S. Di Domizio (18 and 17), S. Farrens (19), L. W. K. Goh (19), S. Gouyou Beauchamps (20 and 21), S. Ilić (22 and 23), S. Joudaki (24), F. Keil (23), A. M. C. Le Brun (2), M. Martinelli (13 and 14), C. Moretti (25 and 26 and 27 and 28), V. Pettorino (4), A. Pezzotta (29), Z. Sakr (30 and 23 and 31), A. G. Sánchez (32), D. Sciotti (13 and 14), I. Tutusaus (21 and 20 and 23), V. Ajani (19 and 33 and 34), M. Crocce (21 and 20), A. Fumagalli (25), C. Giocoli (35 and 36), L. Legrand (37 and 38), M. Lembo (39 and 40 and 41), G. F. Lesci (42 and 35), D. Navarro Girones (5), A. Nouri-Zonoz (43), S. Pamuk (44), A. Pourtsidou (7 and 45), M. Tsedrik (7 and 45), J. Bel (46), C. Carbone (9), C. A. J. Duncan (7), M. Kilbinger (19), D. Sapone (47), E. Sellentin (48 and 5), P. L. Taylor (49 and 50), L. Amendola (30), S. Andreon (29), N. Auricchio (35), C. Baccigalupi (26 and 25 and 27 and 28), M. Baldi (51 and 35 and 36), S. Bardelli (35), P. Battaglia (35), A. Biviano (25 and 26), E. Branchini (18 and 17 and 29), M. Brescia (52 and 53), V. Capobianco (12), J. Carretero (24 and 54), M. Castellano (13), G. Castignani (35), S. Cavuoti (53 and 55), K. C. Chambers (56), A. Cimatti (57), C. Colodro-Conde (58), G. Congedo (7), C. J. Conselice (59), L. Conversi (60 and 61), Y. Copin (62), F. Courbin (63 and 64 and 65), H. M. Courtois (66), M. Cropper (67), A. Da Silva (68 and 69), H. Degaudenzi (70), S. de la Torre (71), G. De Lucia (25), H. Dole (72), M. Douspis (72), F. Dubath (70), X. Dupac (61), S. Dusini (73), S. Escoffier (74), M. Farina (75), F. Faustini (13 and 76), S. Ferriol (62), F. Finelli (35 and 77), M. Frailis (25), E. Franceschi (35), M. Fumana (9), S. Galeotta (25), K. George (78), W. Gillard (74), B. Gillis (7), J. Gracia-Carpio (32), A. Grazian (79), F. Grupp (32 and 80), S. V. H. Haugan (81), H. Hoekstra (5), W. Holmes (82), I. M. Hook (83), F. Hormuth (84), A. Hornstrup

Extracting cosmological information from the Euclid galaxy survey will require modelling numerous systematic effects during the inference process. This implies varying a large number of nuisance parameters, which have to be marginalised over before reporting the constraints on the cosmological parameters. This is a delicate process, especially with such a large parameter space, which could result in biased cosmological results. In this work, we study the impact of different choices for modelling systematic effects and prior distribution of nuisance parameters for the final Euclid Data Release, focusing on the 3$\times$2pt analysis for photometric probes and the galaxy power spectrum multipoles for the spectroscopic probes. We explore the effect of intrinsic alignments, linear galaxy bias, magnification bias, multiplicative cosmic shear bias and shifts in the redshift distribution for the photometric probes, as well as the purity of the spectroscopic sample. We find that intrinsic alignment modelling has the most severe impact with a bias up to $6\,\sigma$ on the Hubble constant $H_0$ if neglected, followed by mis-modelling of the redshift evolution of galaxy bias, yielding up to $1.5\,\sigma$ on the parameter $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m} /0.3}$. Choosing a too optimistic prior for multiplicative bias can also result in biases of the order of $0.7\,\sigma$ on $S_8$. We also find that the precision on the estimate of the purity of the spectroscopic sample will be an important driver for the constraining power of the galaxy clustering full-shape analysis. These results will help prioritise efforts to improve the modelling and calibration of systematic effects in Euclid.

Yiming Ren, Kwan Chuen Chan, Le Zhang, Yin Li, Haolin Zhang, Ruiyu Song, Yan Gong, Xian-Min Meng, Xingchen Zhou

Accurate photometric redshift (photo-$z$) estimation is a key challenge in cosmology, as uncertainties in photo-$z$ directly limit the scientific return of large-scale structure and weak lensing studies, especially in upcoming Stage IV surveys. The problem is particularly severe for faint galaxies with sparse spectroscopic training data. In this work, we introduce nflow-$z$, a novel photo-$z$ estimation method using the powerful machine learning technique of normalizing flow (NF). nflow-$z$ explicitly models the redshift probability distribution conditioned on the observables such as fluxes and colors. We build two nflow-$z$ implementations, dubbed cINN and cNSF, and compare their performance. We demonstrate the effectiveness of nflow-$z$ on several datasets, including a CSST mock, the COSMOS2020 catalog, and samples from DES Y1, SDSS, and DESCaLS. Our evaluation against state-of-the-art algorithms shows that nflow-$z$ performs favorably. For instance, cNSF surpasses Random Forest, Multi-Layer Perceptron, and Convolutional Neutral Network on the CSST mock test. We also achieve a ~30\% improvement over official results for the faint DESCaLS sample and outperform conditional Generative Adversarial Network and Mixture Density Network methods on the DES Y1 dataset test. Furthermore, nflow-$z$ is computationally efficient, requiring only a fraction of the computing time of some of the competing algorithms. Our algorithm is particularly effective for the faint sample with sparse training data, making it highly suitable for upcoming Stage IV surveys.

Nearly every galactic core contains a supermassive compact object, hypothesized to be a Kerr black hole. It was only with the advent of Event Horizon Telescope observations that the predictions of this hypothesis could be observationally tested for our own Galaxy, and the nearby elliptical M87, on spatial scales comparable to the gravitational radius. At the same time it became possible to test whether alternatives such as naked singularities in general relativity, or similar objects in alternative theories of gravity, are excluded by the data. These and other observational developments renewed interest in non-Kerr spacetime metrics, also in the context of active galactic nuclei at cosmological distances. Recently, we have shown that accreting naked singularities in the Reissner-Nordström metric of general relativity tend to produce strong outflows. The geometry and origin of these winds is studied here, and their parameter dependence is investigated. To this end we performed numerical GR hydrodynamical simulations of accretion of electrically neutral matter in the Reissner-Nordström metric and discussed the results in the context of analytic predictions of fluid motion in this spacetime.

The stability of a differentially rotating fluid subject to its own gravity is a problem with applications across wide areas of astrophysics--from protoplanetary discs (PPDs) to entire galaxies. The shearing box formalism offers a conceptually simple framework for studying differential rotation in the local approximation. Aimed at self-gravitating, and importantly, vertically stratified PPDs, we develop two novel methods for solving Poisson's equation in the framework of the shearing box with vertical vacuum boundary conditions (BCs). Both approaches naturally make use of multi-dimensional fast Fourier transforms for computational efficiency. While the first one exploits the linearity properties of the Poisson equation, the second, which is slightly more accurate, consists of finding the adequate discrete Green's function (in Fourier space) adapted to the problem at hand. To this end, we have revisited the method proposed by Vico et al. (2016) and have derived an analytical Green's function satisfying the shear-periodic BCs in the plane as well as vacuum BCs, vertically. Our spectral method demonstrates excellent accuracy, even with a modest number of grid points, and exhibits third-order convergence. It has been implemented in the NIRVANA-III code, where it exhibits good scalability up to 4096 CPU cores, consuming less than 6% of the total runtime. This was achieved through the use of P3DFFT, a fast Fourier Transform library that employs pencil decomposition, overcoming the scalability limitations inherent in libraries using slab decomposition. We have introduced two novel spectral Poisson solvers that guarantees high accuracy, performance, and intrinsically support vertical vacuum boundary conditions in the shearing-box framework. Our solvers enable high-resolution local studies involving self-gravity, such as MHD simulations of gravito-turbulence or gravitational fragmentation.

How planetary systems form and evolve is a key question in astronomy. Revealing how host star properties, such as elemental abundances, age, and mass, differ from those of non-host stars, and how they correlate with planetary characteristics such as radius, provides new insights into the formation and evolutionary pathways of planetary systems. We determine precise ages for 18890 dwarfs and subgiants from the LAMOST-Kepler-Gaia sample with a mean age uncertainty of about 15 percent (median about 10 percent). Within the framework of Galactic chemical evolution, we find that about 86 percent of planet-hosting stars younger than 8 Gyr occupy the upper branch ([Fe/H] > -0.2) of the characteristic V-shaped age-metallicity relation of the Galactic disk. Based on guiding radii (Rg), we further infer that about 19 percent of these young hosts likely originated in the inner disk and subsequently migrated to the solar neighborhood. Among stars older than 10 Gyr, host stars tend to be more metal-rich, with nearly 59 percent having [Fe/H] > -0.2. This suggests that both young and old planet-hosting stars preferentially form in relatively metal-rich environments. However, for host stars with [Fe/H] < -0.2, we find that their metallicities are on average lower by about 0.16 dex compared to non-host stars of similar age and mass, indicating that [Fe/H] is unlikely to be the dominant factor governing planet formation in metal-poor environments. We also identify a systematic depletion of volatile elements, especially carbon, in planet hosts. Moreover, host star [Fe/H] exhibits a weak correlation with planet radius, while [alpha/Fe] primarily support the formation of small planets.

Nikos Sioulas, Marco Velli, Chen Shi, Trevor A. Bowen, Alfred Mallet, Andrea Verdini, B. D. G. Chandran, Anna Tenerani, Jean-Baptiste Dakeyo, Stuart D. Bale, Davin Larson, Jasper S. Halekas, Lorenzo Matteini, Victor Réville, C. H. K. Chen, Orlando M. Romeo, Mingzhe Liu, Roberto Livi, Ali Rahmati, P. L. Whittlesey

We analyze \textit{Parker Solar Probe} and \textit{Solar Orbiter} observations to investigate the propagation and dissipation of Alfvénic fluctuations from the outer corona to 1~AU. Conservation of wave-action flux provides the theoretical baseline for how fluctuation amplitudes scale with the Alfvén Mach number $M_a$, once solar-wind acceleration is accounted for. Departures from this scaling quantify the net balance between energy injection and dissipation. Fluctuation amplitudes follow wave-action conservation for $M_a < M_a^{b}$ but steepen beyond this break point, which typically lies near the Alfvén surface ($M_a \approx 1$) yet varies systematically with normalized cross helicity $\sigma_c$ and fluctuation scale. In slow, quasi-balanced streams, the transition occurs at $M_a \lesssim 1$; in fast, imbalanced wind, WKB-like scaling persists to $M_a \gtrsim 1$. Outer-scale fluctuations maintain wave-action conservation to larger $M_a$ than inertial-range modes. The turbulent heating rate $Q$ is largest below $M_a^{b}$, indicating a preferential heating zone shaped by the degree of imbalance. Despite this, the Alfvénic energy flux $F_a$ remains elevated, and the corresponding damping length $\Lambda_d = F_a/Q$ remains sufficiently large to permit long-range propagation before appreciable damping occurs. Normalized damping lengths $\Lambda_d/H_A$, where $H_A$ is the inverse Alfvén-speed scale height, are near unity for $M_a \lesssim M_a^{b}$ but decline with increasing $M_a$ and decreasing $U$, implying that incompressible reflection-driven turbulence alone cannot account for the observed dissipation. Additional damping mechanisms -- such as compressible effects -- are likely required to account for the observed heating rates across much of the parameter space.

GRB 250702DBE was time-consequently triggered by GBM onboard the Fermi satellite. It is uncertain which celestial catalog is suitable for this special ultra-long event to belong to. In this paper, we comprehensively investigate the lightcurves obtained by Fermi-GBM detectors. In the energy band of 8-1000 keV, no Quasi-Period Oscillation (QPO) signals are found in the lightcurve of the first burst 250702D, a possible QPO signal of 0.046 Hz corresponding to a period of 21.7 s is found in the lightcurve of the second burst 250702B, and a possible QPO signal of 0.024 Hz corresponding to a period of 41.7 s is found in the lightcurve of last burst 250702E. The significance level of the possible QPO signals is comprehensively examined. In addition, we examine the spectral properties of the sources. In general, a broken power law is suitable for modeling the spectral data from 8 keV to 40 MeV. We qualitatively suggest some kinds of celestial object with the periodic characteristic that might be the progenitors of this unique event.

With the progressive release of data from numerous sky surveys, humanity has entered the era of astronomical big data. Multi-wavelength, multi-method research is playing an increasingly crucial role. Binaries account for a substantial fraction of all stellar systems and research into binaries is of fundamental importance. LAMOST J064137.77+045743.8 has not yet been recorded in the SIMBAD astronomical database. We conducted a comprehensive analysis of LAMOST J064137.77+045743.8 using multi-band spectroscopic, astrometric, and photometric data. The low-resolution spectra from Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) suggest that LAMOST J064137.77+045743.8 is a binary consisting of an A7-type subgiant star ($T_{\rm eff}$ $\sim$ 7500\,K and log\,$g$ $\sim$ 3.9) and a cool red dwarf star. Astrometric data from Globe Astrometric Interferometers for Astrophysics support the binary speculation with a Renormalized Unit Weight Error metric value of 1.9. Additional flux observations in the infrared bands further corroborate the presence of a red dwarf companion. The i-band flare detected by the Zwicky Transient Facility (ZTF) photometric observations bolsters the interpretation of an M-type red dwarf companion. The radial velocity variations in the H$\alpha$ lines from LAMOST medium-resolution spectra and the light curves from ZTF both support the classification of the A7 subgiant as a pulsating star. The binary either has a long orbital period, a non-eclipsing binary orbit, or extremely shallow eclipses. Future asteroseismology studies will further probe the internal physics of the A7 subgiants. Research on binaries is incredibly fascinating.

We investigate how primordial magnetic fields (PMFs) affect the formation kinetics of the first molecules, H$_2$, HD, and HeH$^+$, as well as the populations of rovibrational levels and the global signals in the rovibrational transitions of H$_2$ and HD. We show that PMFs can significantly speed up the formation and destruction of the first molecules, leading to an increase in the number density of H$_2$ and HD molecules and a decrease in the number density of HeH$^+$ ion-molecules compared to the case without PMFs. We demonstrate that more frequent collisions of the gas particles in such models alter the ortho-to-para ratio of hydrogen molecules, making it a potential probe of the thermal history of gas in the early Universe. In contrast to the standard cosmological model, where the global signal from the first molecules appears as an absorption feature in the cosmic microwave background spectrum, cosmological models with PMFs can produce an emission signal. Specifically, for non-helical PMFs with $n_B = -2.9$ and a strength of $\sim 1$~nG, the signal transforms into emission with an amplitude of about $\sim 0.5$~Jy/sr. This signal is comparable in magnitude to other known CMB spectral distortions and falls within the detection capabilities of several proposed missions, including Super-PIXIE, Multi-SIMBAD (4 units), and Voyage2050. We show that both the amplitude and the spectral range of the global signals from the first molecules are highly sensitive to the spectral index $n_B$, the strength $B_0$, and the helicity of the PMFs. Therefore, the global signals from the first molecules can serve as a potential probe of PMFs.

Clara C. de la Casa, Kelley M.Hess, Lourdes Verdes-Montenegro, Ralf Kotulla, Hao Chen, Tom H. Jarrett, Michelle E. Cluver, Simon B. De Daniloff, Marie-Lou Gendron-Marsolais, Claude Carignan, John S. Gallagher III, Renée C. Kraan-Korteweg, Roger Ianjamasimanana

We present a new catalog of 196 galaxies of the nearby Hydra I cluster out to $\sim$1.75$\rm r_{200}$, consisting of broad u,g,r,i,z along with narrowband H${\alpha}$ measurements. These deep optical images were obtained with the DECam camera (CTIO) and reach down to a surface brightness limit of $\mu( 3\sigma;10''\times10'')$=26.9 mag $\rm arcsec^2$ in the g band. We also report the HI properties for 89 cluster members detected with MeerKAT. A color magnitude diagram (CMD) shows a bimodal distribution typical of a cluster population, more evolved than those found in isolation. We combine optical H${\alpha}$ and WISE infrared data to compare the star formation history at two distinct timescales. Differences in the star forming activity depicted by both populations manifest as starburst in 24 found members. Of these, 18 starburst galaxies have neutral gas measurements, and show disturbed HI disks that suggest an environmentally-triggered boost in star formation within the last 10$^7$ yrs. Processes such as ram pressure stripping or tidal interactions may underlie their enhanced star-forming activity and asymmetric disks. Since Hydra's dynamical history is unclear, we examine the spatial and velocity distribution of the sample. We reveal a possible link between the large scale structure feeding the Hydra I cluster and the heightened star-forming activity of the starburst galaxies. This feeding pattern matches the few substructure that has been identified in Hydra in previous works, and may explain its origin. Our results portray a picture of a cluster with an evolved nature, plus a population of new infalling galaxies that manifest the impact of their first contact with the cluster environment through star formation, color, morphology and gas content transformations.

The origin of extragalactic fast X-ray transients (EFXTs) remains a fundamental open question in high-energy astrophysics. The Einstein Probe (EP) mission provides a transformative opportunity to investigate their nature. While mounting observations of EP-discovered EFXTs (EP-EFXTs) suggest a possible connection to long gamma-ray bursts (lGRBs), an in-depth comparative analysis between them remains lacking. Here, we present a comparative analysis of their cosmic formation histories, revealing that EP-EFXTs and lGRBs share a similar evolutionary trend-showing a marked decline at $z<1.0$ and a plateau beyond $1.0<z<5$-which clearly distinguishes them from short GRBs. This result is derived from a rigorously selected sample of EP-EFXTs, using Lynden-Bell's $c^{-}$ method to reconstruct, for the first time, the luminosity function and formation rate of EP-EFXTs without any assumptions. Our findings provide independent evidence that EP-EFXTs and lGRBs may originate from a common progenitor channel.

Multiple structures within stellar groups are an intriguing subject for theoretical and observational studies of stellar formation. With the accuracy and completeness of data from Gaia Data Release 3, we now have new opportunities to detect reliable members of stellar groups across a larger field of view than in previous studies. In this work, using machine learning methods and high-accuracy data, we investigate the possibility of detecting multiple structures within 500 arcmin of ASCC 32. We first applied DBSCAN to proper motion and parallax, as multiple structures tend to share similar values for these parameters. Next, we applied GMM to position, proper motion, and parallax for the members detected by DBSCAN. This approach allowed us to identify a filamentary structure among the DBSCAN-detected members. This structure contains all stellar groups previously identified in this region. Subsequently, based on the BIC score, we applied GMM to this filamentary structure. Since multiple structures exhibit distinct positional distributions, GMM was able to effectively separate all groups within the filament. Our methods successfully identified ASCC 32, OC 0395, and HSC 1865 within a 500 arcmin radius. Additionally, we found two distinct substructures within ASCC 32. These four groups exhibit a single main-sequence distribution in the CMD, with proper motion values within three times the standard deviation and slightly differing parallax values, despite having distinct spatial structures. Furthermore, these four groups share the same radial velocity distribution. We provide documentation demonstrating the formation of these stellar groups as a multiple structure, with improved membership identification compared to previous studies.

O. Podladchikova, A. Warmuth, L. Harra, L. Dolla, C. Verbeeck, M. Mierla, L. Rodriguez, S. Parenti, F. Auchere, M. K. Georgoulis, S. J. Hofmeister, N. Engler, M. J. West, R. Aznar Cuadrado, A. M. Veronig, P. Antolin, S. Purkhart, D. M. Long, E. Buchlin, M. Haberreiter, M. Velli, A. N. Zhukov, H. Safari, U. Schuhle, L. Teriaca, A. F. Battaglia, E. Soubrie, V. Buchel, L. P. Chitta, S. Gissot, A. De Groof, M. Gyo, J. P. Halain, B. Inhester, E. Kraaikamp, D. Mueller, D. Pfiffner, P. Rochus, F. Schuller, P. J. Smith, W. Schmutz, K. Stegen

X-ray observations of the Sun led Eugene Parker to propose nanoflares as the basic energy-release units that heat the solar corona. Decades later, Solar Orbiter's Extreme Ultraviolet Imager (HRIEUV), operating halfway between Earth and the Sun, revealed thousands of even smaller brightenings in the quiet corona - tiny "campfires" that are smaller and far more frequent than the fundamental nanoflares observed from 1 AU. We analyze over 12,000 of these events, deriving their thermal energies using multiple geometric models to account for volume uncertainties. Although absolute values vary, all models yield consistent power-law energy distributions and ranges, confirming their flare-like behavior. These picoflares, spanning 10^20-10^24 erg, were detected exactly between the Sun and Earth. They occur up to sixty times more often than nanoflares seen from Earth orbit and supply about 1% of the quiet-Sun coronal heating power. This previously unseen energy source may be a missing component in the solar energy balance. Their discovery extends the flare energy spectrum to smaller scales, and future Solar Orbiter observations at 0.28 AU may reveal the most fundamental flare events that sustain the million-degree solar corona.

Submoons, moons orbiting other moons, may be exotic environments capable of hosting extraterrestrial life. We extend previous studies to revise the maximum lifetime of these objects due to planetary, lunar and sublunar tidal migration. Using the Euler-Lagrange equation with a tidal dissipation process as specified by the Constant Geometric Lag model, we derive and solve the governing equations numerically to map the semi-major axis parameter space for star-planet-moon-submoon systems in which the submoon could be massive enough to host life. We find that Earth could have hosted asteroid-sized submoons ($\sim10^{15}\mathrm{kg}$), whereas a submoon near the previously proposed upper limit ($\sim4.6\cdot10^{17}\mathrm{kg}$) would have driven the Moon $\sim30\%$ farther from Earth than its current orbit. A Warm Jupiter system like Kepler1625 has greater potential of hosting a massive submoon. We found that a submoon of around $10\%M_{\text{Luna}}$ could survive if Kepler1625b's hypothesized moon were $68\%$ farther away then what the best-fit model suggests ($67R_{\mathrm{p}}$ instead of $40R_{\mathrm{p}}$). Giant submoons of mass $1.8M_{\oplus}$ are stable in a Kepler1625-like system. In these cases, the moon orbit is wide ($> 100R_{\mathrm{p}}$). Decreasing the submoon mass to a habitability prerequisite of $0.5M_{\oplus}$, likely needed for a stable atmosphere and plate tectonics, leads to a smaller total number of stable iterations relative to the $m_{sm}=1.8M_{\oplus}$ case. In fact, we identified a minimum number of stable iterations on intermediate submoon mass-scales of around $0.1M_{\oplus}$. This is likely due to an interplay between small tidal forces at small submoon masses and small Roche-Limits at very high submoon masses. If submoon formation pathways in Warm Jupiter systems prefer such intermediate mass-scales, habitable submoons could be a rare phenomenon.

We present a timing analysis of \textit{Insight}-HXMT observations of the black-hole X-ray binary Swift J1727.8$-$1613 across a bright soft X-ray flare on 2023 September 19 (MJD 60206). At the peak of the flare, the source undergoes a brief transition from the hard-intermediate state (HIMS) into the soft-intermediate state (SIMS), marked by the simultaneous appearance of three discrete radio jet ejections, a drop in broadband noise in the 2$-$10 keV band, and the presence of a narrow quasi-periodic oscillation (QPO) with a characteristic ``U''-shaped phase-lag spectrum and a quality factor of $Q \geq 6$, features that robustly identify it as a Type-B QPO. The Type-C QPO, which was clearly detected in the HIMS prior to the flare, is not observed at the flare's peak and only reappears afterward. Most notably, we find that the Type-B QPO is not restricted to the SIMS: it is present throughout all our observations, including those taken in the HIMS, where it appears as a broad shoulder of the Type-C QPO. During the flare, the Type-B and Type-C QPOs exhibit distinct evolutionary trends in frequency, fractional rms amplitude, and phase lag. These results challenge the traditional view that Type-B QPOs are exclusive to the SIMS, a state that is, in fact, defined by their appearance in the power spectrum, and directly linked to discrete jet ejections. Instead, our findings suggest that the physical conditions giving rise to Type-B QPOs occur more broadly within the inner accretion flow.

Rinka Ito, Yusuke Miyamoto, Naomasa Nakai, Aya Yamauchi, Yuichi Terashima

We present the results of very long baseline interferometry (VLBI) observations of water vapor masers in the nucleus of the LINER galaxy NGC 7738. The red- and blue-shifted and newly detected systemic maser features show an almost edge-on disk located at a distance of ${0.031}\mbox{-}{0.222}$ pc from the galactic center and rotating with a velocity of $324\mbox{-}454$ km s$^{-1}$ . The velocity field of the disk indicates sub-Keplerian rotation, suggesting a non-negligible disk mass. The Mestel disk model reveals the central and disk masses to be $(1.2 \pm 0.4) \times 10^6$ $M_{\odot}$ and $(4.7 \pm 1.5) \times10^6$ $M_{\odot}$, respectively. The mean volume density within the inner radius of the disk [$(1.2 \pm 0.5) \times 10^{10}$ $M_{\odot}$ $\mathrm{pc^{-3}}$] strongly suggests the existence of a supermassive black hole at the center.

Michael K Florian, Michael D Gladders, Gourav Khullar, Keren Sharon, Aidan P Cloonan, Eirk Solhaug, Brian Welch, Matthew Bayliss, Hakon Dahle, Taylor A Hutchison, Jane R Rigby

JWST has enabled the discovery of a statistical sample of obscured (type II) active galactic nuclei (AGN) at cosmic noon. Studies comparing those type II AGN with type I AGN at that epoch have reinforced the long-standing idea of an evolutionary link between those classes of objects. Mergers, the idea goes, disturb the morphologies and angular momentum of galaxies. The disruption of angular momentum allows material to be funneled toward galactic cores, sparking AGN activity and potentially also a burst of star-formation. That material enshrouds the galactic nucleus, leading to a type II AGN. Later, AGN feedback clears the circumnuclear dust, leading to a transition into a type I AGN, and also quenches star formation. If this is a common outcome, a class of intermediate objects should exist. Such objects would be somewhat disturbed and dusty and sit below the star-forming galaxy main sequence, and their star-formation histories would show an increase in star-formation at around the time of the suspected merger. We present new JWST observations of SDSSJ2222+2745, a strongly lensed AGN at z=2.801. The lensing magnification enables a detailed study of the host galaxy spanning the rest-ultraviolet through near infrared. JWST and HST photometry, morphological models, and models of the host's spectral energy distribution reveal that SDSSJ2222+2745 is actively transitioning from a type II to type I AGN. Catching a lensed AGN at this special evolutionary phase makes SDSSJ2222+2745 a unique laboratory to study the physical processes involved in the transition and their relationships to the AGN and the host galaxy at incredible spatial-resolution down to about 20pc at z=2.801.

M. Nizovkina, S. S. Larsen, A. G. A. Brown, A. Helmi

The precision of cluster parameter determination has significantly improved with the availability of homogeneous photometric Gaia data, however, challenges such as age-metallicity degeneracy and lack of spectroscopic observations remain. In this paper we investigate whether metallicities derived from low-resolution Gaia XP spectra can be effectively used to break degeneracies and improve the accuracy of OC parameter determinations. We analysed 20 OCs using isochrone fitting methods on Gaia DR3 photometry and metallicity estimates from several Gaia XP-based catalogues. We derived age, distance modulus, and extinction using the Approximate Bayesian Computation (ABC). We compared the parameter estimates to the values obtained in other works through isochrone fitting with spectroscopically constrained metallicities or through neural network techniques applied only to the photometry. We found the systematic difference between Gaia XP derived metallicities and those obtained from high-resolution spectroscopy to be 0.1-0.15 dex. We found a systematic age difference of <0.03 +- 0.13 dex compared to isochrone fitting using high-resolution spectroscopy, and <0.08 +- 0.21 dex compared to neural network-based methods, and a median individual error of ~0.065 dex. Despite their low resolution, Gaia XP metallicities effectively constrain parameters of clusters lacking a well-populated RGB. When used with stringent quality cuts and incorporated as priors, they allow to determine ages comparable in precision to those based on high-resolution spectroscopy and more precise than photometry-only neural network methods. These results highlight the potential of Gaia data for accurate cluster parameter analysis and detailed Galactic studies without relying on traditional spectroscopy.

Zhaoran Liu, Tadayuki Kodama, Brian C. Lemaux, Mariko Kubo, Jose Manuel Pérez-Martínez, Yusei Koyama, Ichi Tanaka, Kazuki Daikuhara, Roy R. Gal, Denise Hung, Masahiro Konishi, Kosuke Kushibiki, Ronaldo Laishram, Lori M. Lubin, Kentaro Motohara, Hidenori Takahashi

We present results from a dual narrow-band imaging survey targeting the CL1604 supercluster at z = 0.9 using the Subaru Telescope. By combining the NB921 filter on HSC and the NB1244 filter on SWIMS, we can detect redshifted H$\alpha$ and H$\beta$ emission lines from the supercluster. This unique technique allows us to measure both star formation rates and dust extinction for a sample of 94 emission-line galaxies across the supercluster. We find that dust extinction, estimated from the Balmer decrement (H$\alpha$/H$\beta$ ratio), increases with stellar mass in star-forming galaxies, whereas relatively quiescent systems exhibit comparatively low extinction. Among galaxies with intermediate masses ($10^{8.5} < M_* < 10^{10.5}\,M_\odot$), the dust-corrected H$\alpha$-based star formation rates align with the main sequence at this epoch. More massive galaxies, however, deviate from this relation, exhibit redder colors, and reside predominantly in higher-density environments. Although stellar mass, SFR, and galaxy color are clearly influenced by environment, we detect no strong, systematic environmental dependence of dust extinction for the whole sample.

Taehyun Kim, Dimitri A. Gadotti, Myeong-gu Park, Yun Hee Lee, Francesca Fragkoudi, Minjin Kim, Woong-Tae Kim

Dark gaps, low surface brightness regions along the bar minor axis, are expected to form as a consequence of secular evolution in barred galaxies. Although several studies have proposed links between dark gap locations and dynamical resonances, the results remain inconclusive. Using DESI Legacy Imaging Survey data, we find that approximately 61% of barred galaxies exhibit pronounced dark gaps. We compare the location of dark gaps with resonance radii derived from the Tremaine-Weinberg method applied to MaNGA data for the same galaxies. Our analysis shows that dark gaps do not preferentially form at specific resonances. Instead, their locations correlate with $\mathcal{R}$ $\equiv$ $R_{CR}/R_{Bar}$: slow bars tend to show shorter dark gap radii, while fast bars show longer ones. This trend reflects a tight relation between bar length and dark gap radius. However, when barred galaxies are classified by their ring morphology, certain types exhibit dark gaps that align with specific resonances. Notably, dark gaps located between the inner and outer rings are closely associated with the corotation radius. In galaxies with two dark gaps along the bar minor axis profile, the inner dark gap typically aligns with the ultraharmonic resonance, and the outer dark gap corresponds to the corotation radius. These findings suggest that some morphological types share similar $\mathcal{R}$ values and exhibit dark gaps near specific resonances. Thus, dark gaps may serve as proxies for dynamical resonances only in certain systems. Our findings may help explain the discrepancies observed in earlier studies.

The measurements of baryon acoustic oscillation by the Dark Energy Spectroscopic Instrument Data Release 2 indicate that dark energy may be a dynamical quantity with a time-varying equation of state. This challenges the core assumptions of the $\Lambda$CDM model and has generated significant interest in dynamical dark energy models. Therefore, studying the parameterization of the equation of state for dynamical dark energy is crucial. Existing work has achieved fruitful results in the dark energy models, exploring various parameterization forms, but it is relatively scattered and lacks systematic parameter constraints based on the latest dataset combinations. We use the $\Lambda$CDM as a baseline model and carry out rigorous statistical constraints on key cosmological parameters for seven representative parameterization models. Planck PR4 and DESI DR2 observations are incorporated into our study. We use three dataset combinations: CMB+BAO+PantheonPlus, CMB+BAO+DES-Y5, and CMB+BAO+Union3. The ${H}_{0}$ and ${\sigma }_{8}$ values of all dynamical dark energy models are lower than the $\Lambda$CDM model, indicating that our results may not effectively alleviate ${H}_{0}$ tension, but can significantly reduce ${\sigma }_{8}$ tension. By comparing the $\chi^2$ and the Akaike Information Criterion obtained for each model, we demonstrate that the linear Chevallier-Polarski-Linder parameterization model is not the optimal choice in all cases. Specifically, when combined with the CMB+BAO+DES-Y5 dataset, the Barboza-Alcaniz, Logarithmic, and Exponential models demonstrate superior statistical fitting performance compared to the $\Lambda$CDM model. The Barboza-Alcaniz model shows a great advantage in fitting performance, leading to the most significant improvement.

Rapidly rotating neutron stars (NSs) are promising targets for continuous gravitational-wave (CGW) searches with current and next-generation ground-based GW detectors. In this Letter, we present the first study of thermal deformations in super-Eddington magnetized NSs with column accretion, where magnetic fields induce anisotropic heat conduction that leads to crustal temperature asymmetries. We compute the resulting mass quadrupole moments and estimate the associated CGW strain amplitudes. Our results show that Galactic magnetized NSs undergoing super-Eddington column accretion can emit detectable CGWs in upcoming observatories. Assuming a 2-yr coherent integration, the Einstein Telescope and Cosmic Explorer could detect such CGW signals from rapidly spinning NSs with spin periods $P \lesssim 20\,\rm ms$, while the LIGO O5 run may detect systems with $P \lesssim 6 \,{\rm ms}$. These findings suggest that super-Eddington magnetized NSs could represent a new class of CGW sources, providing a unique opportunity to probe the NS crust and bridge accretion physics with GW astronomy.

Yanee Tangjai, Agnieszka Leszczynska, Tatphicha Promfu, Achara Seripienlert, Serap Tilav (for the IceCube Collaboration)

This study evaluates the response of the IceTop tanks to low-energy air showers in the GeV to TeV energy range based on simulated and measured count rates. Correlating this response with primary cosmic rays provides a tool to study Galactic and solar cosmic-ray flux modulations, particularly for solar particle events. We present long-term behavior of the IceTop scaler rates for a range of discriminator thresholds to better understand and calibrate the detector's response to changing environmental conditions.

Jiaxin Liu, Haoyu Yuan, Xiangli Lei, Wenlong Xu, Jumpei Takata, Weihua Lei

We investigate the optical linear polarization caused by Thomson scattering of the stellar radiation for gamma-ray binary \lsi61, which likely contains a young pulsar. Based on the pulsar binary scenario, we model the interaction between the pulsar wind and stellar wind from the massive companion star, which creates a shock. To accurately compute the resulting polarization of the stellar wind, we develop a method for the Thomson scattering that accounts for the finite size of the companion star. By fitting the optical polarization data, we constrain the system parameters, such as eccentricity, the momentum ratio of the two winds, and mass-loss rate from the companion star. We find that (i) the predicted eccentricity $e\sim 0.1$ is smaller than the values derived from the radial velocity curve and (ii) the orbital phase of the periastron is $\nu_{\rm p}=0.5-0.6$, which is consistent with the previous polarization study of Kravtsov et al. Additionally, we estimate the mass-loss rate from the companion star and the momentum ratio of two winds as $\dot{M}\sim 2\times 10^{-6}\rm M_{\odot}~{\rm year^{-1}}$ and $\eta>0.1$, respectively. Assuming that the pulsar wind carries the spin-down energy, the spin-down magnetic field of the putative pulsar inferred from these parameters is of the order of $B\sim 10^{14}\mathrm{G}$, which may support the highly-B pulsar or magnetar scenario for the compact object of $\rm{LS\ I} +61^{\circ}303$. We also discuss the dispersion measure under the predicted orbital geometry and provide a corresponding interpretation of the pulsed radio signal detected by FAST.

We present a refined estimation of the stochastic gravitational wave background (SGWB) based on observed dual active galactic nuclei (AGNs) together with AGN X-ray luminosity functions, in light of recent Pulsar Timing Array detections of an nHz SGWB. We identify a characteristic luminosity dependence in dual AGN fractions by compiling recent observational datasets, providing crucial constraints on supermassive black hole binary (SMBHB) populations. Our AGN-based model reproduces the current SGWB measurements within PTA observational uncertainties of $2 - 4 \sigma$ uncertainties, demonstrating consistency between electromagnetic and gravitational wave observations. These findings establish SMBHBs as the dominant source of the nHz gravitational wave signal, providing crucial insights into their demographics and evolution.

Jaime E. Pineda, Rachel K. Friesen, Erik Rosolowsky, Ana Chacón-Tanarro, Michael Chun-Yuan Chen, James Di Francesco, Helen Kirk, Anna Punanova, Youngmin Seo, Yancy Shirley, Adam Ginsburg, Stella S. R. Offner, Ayush Pandhi, Ayushi Singh, Feiyu Quan, Héctor G. Arce, Paola Caselli, Spandan Choudhury, Alyssa A. Goodman, Fabian Heitsch, Peter G. Martin, Christopher D. Matzner, Philip C. Myers, Elena Redaelli, Samantha Scibelli

We present an overview of the final data release (DR2) from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_\mathrm{V} \gtrsim 7$~mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This final release includes the data for all the regions observed: Heiles Cloud 2 and B18 in Taurus; Barnard 1, Barnard 1-E, IC348, NGC 1333, L1448, L1451, and Per7/34 in Perseus; L1688 and L1689 in Ophiuchus; Orion A (North and South) and Orion B in Orion; Cepheus, B59 in Pipe; Corona Australis (CrA) East and West; IC5146; and Serpens Aquila and MWC297 in Serpens. Similar to what was presented in GAS DR1, we find that the NH$_3$ emission and dust continuum emission from Herschel correspond closely. We find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores, and we find that the transition between coherent core and turbulent cloud is a common result. This shows that the regions of coherence are common throughout different star forming regions, with a substantial fraction of the high column density regions displaying subsonic non-thermal velocity dispersions. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, $\sigma_v$, and gas kinetic temperature, $T_{\rm kin}$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of $\sigma_v$ and $T_{\rm kin}$ with increasing star formation activity. The data presented in this paper are publicly available via \dataset[DOI: https://doi.org/10.11570/24.0091]{this https URL}.

Deep learning has generated diverse perspectives in astronomy, with ongoing discussions between proponents and skeptics motivating this review. We examine how neural networks complement classical statistics, extending our data analytical toolkit for modern surveys. Astronomy offers unique opportunities through encoding physical symmetries, conservation laws, and differential equations directly into architectures, creating models that generalize beyond training data. Yet challenges persist as unlabeled observations number in billions while confirmed examples with known properties remain scarce and expensive. This review demonstrates how deep learning incorporates domain knowledge through architectural design, with built-in assumptions guiding models toward physically meaningful solutions. We evaluate where these methods offer genuine advances versus claims requiring careful scrutiny. - Neural architectures overcome trade-offs between scalability, expressivity, and data efficiency by encoding physical symmetries and conservation laws into network structure, enabling learning from limited labeled data. - Simulation-based inference and anomaly detection extract information from complex, non-Gaussian distributions where analytical likelihoods fail, enabling field-level cosmological analysis and systematic discovery of rare phenomena. - Multi-scale neural modeling bridges resolution gaps in astronomical simulations, learning effective subgrid physics from expensive high-fidelity runs to enhance large-volume calculations where direct computation remains prohibitive. - Emerging paradigms-reinforcement learning for telescope operations, foundation models learning from minimal examples, and large language model agents for research automation-show promise though are still developing in astronomical applications.

Kshitij Duraphe, Kartik Mandar, Chooda Khanal, Abha Pareek, Tejaswi Kondhiya, V Sree Suswara, Deeksha Dinesh, Vidyasagar Bhat, Gopal Bhatta

We present a comprehensive timing analysis of the black hole X-ray binary Cygnus X-1 using 26 NuSTAR observations spanning 2012-2024, providing the most detailed characterization to date of its accretion flow variability across spectral states. Our analysis reveals fundamental insights into the physics governing state transitions in stellar-mass black holes. We discover distinct bimodal flux distributions in the 8-79 keV band with well-separated peaks, contrasting with overlapping distributions in the 3-8 keV band. This energy-dependent bimodality establishes hard X-rays as the optimal diagnostic for state classification, directly tracing the geometric transformation between corona-dominated and disk-dominated configurations. Power spectral analysis uncovers state-dependent characteristic frequencies shifting from 0.050 Hz (hard) to 0.074 Hz (intermediate), with featureless red noise in soft states. These frequencies correspond to disk truncation radii evolving from $\sim$5.5 $R_g$ to $\sim$2 $R_g$, providing direct observational evidence for the inward progression of the accretion disk during state transitions. Frequency-dependent time lags evolve systematically from $\sim$50 ms hard lags at 0.1 Hz in hard states to near-zero in soft states, quantifying the collapse of the Comptonizing corona. Linear rms-flux relations persist across all states with parameters that precisely track the relative contributions of thermal versus non-thermal emission components. Most remarkably, we identify a failed state transition (observation 30302019006) exhibiting anticorrelated band behavior, suppressed variability ($F_{var}$ < 1.38\%), and apparent sub-ISCO truncation. This discovery challenges standard transition models and suggests new pathways for accretion flow evolution in wind-fed systems.

Lucía Bravo Ferres, Francisco Nogueras-Lara, Rainer Schödel, Rubén Fedriani, Adam Ginsburg, Samuel Crowe, Jonathan C. Tan, Morten Andersen, Joseph Armstrong, Yu Cheng, Zhi-Yun Li

Determining the infrared extinction curve towards the Galactic centre is crucial for accurately correcting observed data and deriving the underlying stellar populations. However, extinction curves reported in the literature often show discrepancies. We aim to derive the infrared extinction curve towards the Galactic centre based on JWST-NIRCam data for the first time, using observations of the Sagittarius C region in the 1-5 $\mu$m range. We determined extinction ratios using two different methods, both based on measuring the reddening vector using the slope of red clump stars, whose intrinsic properties are well known, in observed colour-magnitude diagrams. The extinction curve derived in this work is in good agreement with previous results in the literature. We obtained the following extinction ratios relative to F162M: $A_\mathrm{F115W} : A_\mathrm{F162M} : A_\mathrm{F182M} : A_\mathrm{F212N} : A_\mathrm{F360M} : A_\mathrm{F405N} : A_\mathrm{F470N} : A_\mathrm{F480M} = 1.84 \pm 0.03 : 1.00 : 0.789 \pm 0.005 : 0.607 \pm 0.014 : 0.306 \pm 0.011 : 0.248 \pm 0.017 : 0.240 \pm 0.019 : 0.21 \pm 0.03$. Besides, we found different values of the extinction index for the short- ($\lambda \sim 1-2.5\,\mu$m, $\alpha \sim 2$) and long-wavelength ($\lambda \sim 2.5-5\,\mu$m, $\alpha \sim 1.4$) regimes, with the extinction curve flattening at longer wavelengths. Comparison with extinction curves derived both inside and outside the Galactic centre suggests that the infrared extinction curve does not significantly vary in the central regions, and shows no significant evidence for variations between different lines of sight beyond the inner Galaxy within the uncertainties.

Scattering between dark matter (DM) and protons leads to suppressed small-scale fluctuations, with implications for a variety of cosmological observables. In this work, we search for evidence of DM-proton scattering with an interaction cross section $\sigma\!=\!\sigma_0 (\frac{v}{c})^n$ for $n=0,2$ and $4$, corresponding e.g. to velocity-independent contact interactions from heavy mediators, velocity-dependent pseudoscalar-mediated scattering, and higher-order dipole interactions, respectively, using high-redshift ($z \sim4-10$) ultraviolet galaxy luminosity functions (UVLFs) observed by Hubble Space Telescope (HST). We employ an adjusted implementation of GALLUMI combined with the modified Boltzmann solver CLASS DMeff that accounts for interacting DM, and incorporate UVLF data from both blank and lensed HST fields, alongside Planck CMB data and the Pantheon supernova catalog in a Bayesian analysis framework to set constraints on $\sigma_0$. Our results show that including lensed UVLF data, which probe fainter galaxies than the blank HST fields and thus smaller scales, leads to a substantial improvement in the constraints on $\sigma_0$ for $n>0$, surpassing existing bounds from Milky-Way (MW) satellite abundance and CMB anisotropies. For $m_{\chi} = 1\,\rm MeV $, for example, we set the upper bounds at $8.3\times 10^{-26} \, \rm cm^2$ for $n=2$ and $1.2\times 10^{-22} \, \rm cm^2$ for $n=4$. For $n=0$, our bound is within an order of magnitude of those from the Lyman-$\alpha$ forest and MW satellites.

Greta Zucchi, Xihan Ji, Piero Madau, Roberto Maiolino, Ignas Joudzbalis, Francesco D'Eugenio, Sophia Geris, Yuki Isobe

Observations with the James Webb Space Telescope (JWST) have uncovered a substantial population of high-redshift, broad-line active galactic nuclei (AGNs), whose properties challenge standard models of black hole growth and AGN emission. We analyze a spectroscopic sample of 34 Type 1 AGNs from the JWST Advanced Deep Survey (JADES) survey, spanning redshifts 1.7 < z < 9, to constrain the physical nature of the accretion flows powering these sources with broad-line diagnostics statistically for the first time. At z > 5, we find a marked suppression of high-ionization emission lines (HeII, CIV, NV) relative to prominent broad Halpha and narrow [OIII] features. This contrast places strong constraints on the shape of the ionizing spectral energy distribution (SED) and on the physical conditions in the broad-line region (BLR). By comparing the observations to photoionization models based on SEDs of black holes accreting at sub-Eddington ratios, we show that standard AGN continua struggle to reproduce the observed broad line ratios and equivalent widths across a wide ionization parameter range. These results suggest the need for modified SEDs -- either intrinsically softened due to super-Eddington accretion or radiative inefficiencies in the innermost disk, or externally filtered by intervening optically thick gas that absorbs or scatters the highest-energy photons before they reach the BLR.

We address the challenge of reconstructing the energy of three ultra-high-energy cosmic rays registered with a small fluorescence telescope EUSO-TA that operated in 2015 at the site of the Telescope Array experiment in Utah, USA. Each of these events was recorded within one time frame. Conventional methods of energy reconstruction are not applicable in this case because the events do not have light curves but a single data point. As an alternative, we consider a number of approaches based on artificial neural networks. We demonstrate that a signal recorded by a fluorescence telescope within one time frame might be enough to reconstruct the energy of a primary particle with reasonable accuracy using an ensemble of simple convolutional neural networks. Contrary to the conventional approach, reconstruction of the shower geometry is not needed for this. More than this, preliminary estimates can be obtained even without recognizing the shower track. However, there remain some problems that do not allow us to claim that the suggested method is universal and always works. We discuss difficulties that we faced and possible ways of improving the method.

Cheukyu Edward Tong, Keara Carter, Paul Grimes, Eugene Lauria, Dan Marrone, Gabriella Montano, Matthew Morgan, Yoshirnori Uzawa, Lingzhen Zeng

A dual band receiver has been designed for the Black Hole Explorer (BHEX) mission, which is a space Very-Long Baseline Interferometry (VLBI) mission concept, aimed at unveiling the photon ring of black holes. The cryogenic receiver comprises a 228-320 GHz Superconductor-Insulator Superconductor (SIS) receiver, paired with a 76-106.7 GHz HEMT receiver. The details of the design are described in this talk. A novel comb generator, which will be used for delay tracking, has been designed and tested.

Recently a new meteoritic mineral, Zolenskyite (Fe0.99Mn0.04Ca0.01Cr1.99S3.98), was discovered from the Indarch meteorite. Zolenskyite was structurally indexed as the monoclinic C2/m CrNb2Se4 - Cr3S4 type structure of synthetic FeCr2S4, with unit cell parameters a = 12.84(1) Å, b = 3.44(1) Å, c = 5.94(1) Å and \b{eta} = 117(1)°. Zolenskyite was reported as high-pressure phase formed from Daubréelite at high pressures and temperatures in highly shocked regions of the EH parent asteroid. Although this discovery provides valuable information about the origin of meteoritic mineral assemblages, the results and conclusions raise controversies with those reported in previous articles where the synthetic FeCr2S4 was described. In this review, an alternative analysis of the supplementary X-ray data from Zolenskyite was made and yields to the monoclinic I2/m Cr3S4 type structure of synthetic FeCr2S4, with unit cell parameters a = 5.940 Å, b = 3.440 Å, c = 11.441 Å and \b{eta} = 90.55°, in agreement with previous results in the literature and differ from those reported above. Regarding the genesis of Zolenskyite and according to solid-state laboratory synthesis, the transformation of cubic FeCr2S4 phase (ideal composition of Daubréelite) to monoclinic FeCr2S4 (ideal composition of Zolenskyite) requires a process whose replication in Shock metamorphism stages is not well established and remains an open issue to address, together with another open issue related to the real composition of meteoritic minerals, with minor and trace metals, which is often overlooked when compared with synthetic ideal compositions. Whatever the results to be obtained by addressing the above open issues, they will shed more light on the genesis of Zolenskyite. In this context, it is worth to promote an open debate about the hypothesis of artificial origin.

We explore extreme mass-ratio inspirals (EMRIs) in the co-evolution of massive black holes (MBHs) and nuclear star clusters (NSCs), which host diverse stellar populations across a wide range of masses. The dynamics are simulated self-consistently with GNC, which we have updated to incorporate gravitational wave orbital decay, the loss cone of a spinning MBH, and stellar evolution. Over $12$ Gyr, we investigate the evolution of the NSC with a mass-growing MBH, as well as the EMRIs of stellar black holes, neutron stars, white dwarfs, brown dwarfs (BDs), and low-mass main-sequence stars (MSs), along with tidal disruption events (TDEs) involving MSs, BDs, and post-MSs. The mass growth of the MBH contributed by TDEs is typically $\sim 10^7\,M_{\odot}$, $\sim 10^6\,M_{\odot}$, and $\sim 5\times10^4\,M_{\odot}$ for massive, Milky-Way-like, and smaller NSCs, respectively. Between $40\%$ and $70\%$ of the stellar mass is lost during stellar evolution, which dominates the mass growth of the MBH if a significant fraction of the lost mass is accreted. The evolution of EMRI rates is generally affected by the cluster's size expansion or contraction, stellar population evolution, MBH mass growth, and the stellar initial mass function. The EMRI rates for compact objects peak at early epochs ($\lesssim 1$ Gyr) and then gradually decline over cosmic time. LISA-band ($0.1$ mHz) EMRIs involving compact objects around Milky-Way-like MBHs tend to have high eccentricities, while those around spinning MBHs preferentially occupy low-inclination (prograde) orbits. In contrast, MS- and BD-EMRIs usually have eccentricity and inclination distributions that are distinct from those of compact objects.

Emeline Bolmont, Edward Galantay, Sergi Blanco-Cuaresma, Apurva V. Oza, Christoph Mordasini

We investigate the origin and stability of extrasolar satellites orbiting close-in gas giants, focusing on whether these satellites can survive planetary migration within a protoplanetary disk. To address this question, we used Posidonius, an N-Body code with an integrated tidal model, which we expanded to account for the migration of a gas giant within a disk. Our simulations include tidal interactions between a $1M_\odot$ star and a $1M_{Jup}$ planet, as well as between the planet and its satellite, while neglecting tides raised by the star on the satellite. We adopt a standard equilibrium tide model for the satellite, planet, and star, and additionally explore the impact of dynamical tides in the convective regions of both the star and planet on satellite survival. We examine key parameters, including the initial satellite-planet distance, disk lifetime (proxy for the planet's final orbital distance), satellite mass, and satellite tidal dissipation. For simulations incorporating dynamical tides, we explore three different initial stellar rotation periods. We find that satellite survival is rare if the satellite has nonzero tidal dissipation. Survival is only possible for initial orbital distances of at least 0.6 times the Jupiter-Io separation and for planets orbiting beyond about 0.1 AU. Satellites that fail to survive are either 1) tidally disrupted, as they experience orbital decay and cross the Roche limit, or 2) dynamically disrupted, where eccentricity excitation drives their periastron within the Roche limit. Satellite survival is more likely for low tidal dissipation and higher satellite mass. Given that satellites around close-in planets appear unlikely to survive planetary migration, our findings suggest that if such satellites do exist, another process should be invoked. In that context, we also discuss the claim of the existence of a putative satellite around WASP-49 A b.

Hillary Davis, Thor Tepper-Garcia, Naomi McClure-Griffiths, Joss Bland-Hawthorn, Oscar Agertz

The development of an N-body/hydrodynamic `surrogate' model of the Milky Way (MW) - a model that resembles the MW in several key aspects after many Gyrs of evolution - would be extremely beneficial for Galactic Archaeology. Here we present four new `surrogate' models, all built with the Nexus framework. The simulations contain stars, dark matter and gas. Our most sophisticated model allows gas to evolve thermodynamically, and includes star formation, metal production, and stellar feedback. The other three models in this work have an isothermal gas disc. We examine these new simulations in the context of cold gas observations of the Galaxy. Our focus is the so-called `HI terminal velocity curve' - a heliocentric measurement of the maximum Vlos as a function of Galactic longitude, which dates back to the early days of radio astronomy. It is a powerful approach to indirectly estimating the gas dynamics because it does not require knowledge about the distance to individual gas clouds, which is difficult to estimate. A comparison of the terminal velocities and recovered rotation curve values in the simulations against observations suggests that our models are in need of further refinement. The gravitational torques associated with our synthetic bars are too strong, driving excessive streaming motion in the inner gas disc. This causes the simulated terminal velocity curves in the Galactic Quadrant I and IV to deviate substantially from each other, unlike what is seen in observed HI terminal velocities of the MW. We suggest possible ways forward for future models.

Xiaohan Chen, Man I Lam, Yingying Zhou, Hongrui Gu, Jinzhi Lai, Zhou Fan, Jing Li, Xin Zhang, Hao Tian

Slitless spectroscopy eliminates the need for slits, allowing light to pass directly through a prism or grism to generate a spectral dispersion image that encompasses all celestial objects within a specified area. This technique enables highly efficient spectral acquisition. However, when processing CSST slitless spectroscopy data, the unique design of its focal plane introduces a challenge: photometric and slitless spectroscopic images do not have a one-to-one correspondence. As a result, it becomes essential to first identify and count the sources in the slitless spectroscopic images before extracting spectra. To address this challenge, we employed the You Only Look Once (YOLO) object detection algorithm to develop a model for detecting targets in slitless spectroscopy images. This model was trained on 1,560 simulated CSST slitless spectroscopic images. These simulations were generated from the CSST Cycle 6 and Cycle 9 main survey data products, representing the Galactic and nearby galaxy regions and the high galactic latitude regions, respectively. On the validation set, the model achieved a precision of 88.6% and recall of 90.4% for spectral lines, and 87.0% and 80.8% for zeroth-order images. In testing, it maintained a detection rate >80% for targets brighter than 21 mag (medium-density regions) and 20 mag (low-density regions) in the Galactic and nearby galaxies regions, and >70% for targets brighter than 18 mag in high galactic latitude regions.

Max Charles, Louis Desdoigts, Benjamin Pope, Peter Tuthill, Dori Blakely, Doug Johnstone, Shrishmoy Ray, K. E. Saavik Ford, Barry McKernan, Anand Sivaramakrishnan

Flying on board the James Webb Space Telescope (JWST) above Earth's turbulent atmosphere, the Aperture Masking Interferometer (AMI) on the NIRISS instrument is the highest-resolution infrared interferometer ever placed in space. However, its performance was found to be limited by non-linear detector systematics, particularly charge migration - or the Brighter-Fatter Effect. Conventional interferometric Fourier observables are degraded by non-linear transformations in the image plane, with the consequence that the inner working angle and contrast limits of AMI were seriously compromised. Building on the end-to-end differentiable model & calibration code amigo, we here present a regularised maximum-likelihood image reconstruction framework dorito which can deconvolve AMI images either in the image plane or from calibrated Fourier observables, achieving high angular resolution and contrast over a wider field of view than conventional interferometric limits. This modular code by default includes regularisation by maximum entropy, and total variation defined with $l_1$ or $l_2$ metrics. We present imaging results from dorito for three benchmark imaging datasets: the volcanoes of Jupiter's moon Io, the colliding-wind binary dust nebula WR 137 and the archetypal Seyfert 2 active galactic nucleus NGC 1068. In all three cases we recover images consistent with the literature at diffraction-limited resolutions. The performance, limitations, and future opportunities enabled by amigo for AMI imaging (and beyond) are discussed.

Atmospheric scintillation is one of the largest sources of error in ground-based spectrophotometry, reducing the precision of astrophysical signals extracted from the time-series of bright objects to that of much fainter objects. Relative to the fundamental Poisson noise, scintillation is not effectively reduced by observing with larger telescopes, and alternative solutions are needed to maximize the spectrophotometric precision of large telescopes. If the chromatic covariance of the scintillation is known, it can be used to reduce the scintillation noise in spectrophotometry. This paper derives analytical solutions for the chromatic covariance of stellar scintillation on a large telescope for a given atmospheric turbulence profile, wind speed, wind direction, and airmass at optical/near-infrared wavelengths. To demonstrate how scintillation noise is isolated, scintillation-limited exoplanet transit spectroscopy is simulated. Then, a procedure is developed to remove scintillation noise and produce Poisson-noise limited light curves. The efficacy and limits of this technique will be tested with on sky observations of a new, high spectrophotometric precision, low resolution spectrograph.

We show that core-collapsed self-interacting dark matter halos of mass $\sim 10^6\,{M_\odot}$, originally simulated to explain the dense perturber of the GD-1 stellar stream, also reproduce the structural properties inferred for the dense perturber detected in the strong lensing system JVAS B1938+666 from radio observations. Furthermore, these halos are sufficiently compact and dense to gravitationally capture field stars in satellite galaxies of the Milky Way, providing a natural explanation for the origin of Fornax 6, a stellar cluster in the Fornax dwarf spheroidal galaxy. Our results demonstrate that observations of halos with similar masses but residing in different cosmic environments offer a powerful and complementary probe of self-interacting dark matter.

Nina Bonaventura, Jianwei Lyu, George H. Rieke, Andrew J. Bunker, Chris J. Willott, Christopher N. A. Willmer

In Paper I, we exploited the unsurpassed resolution and depth of JWST/NIRCam imagery to investigate the relationship between AGN and host-galaxy properties in the JWST era, finding a correlation between the level of spatial disturbance (as measured by shape asymmetry, $A_S$) and obscuration ($N_H$). Here in Paper II, we report an expansion of our X-ray and infrared analysis of Seyfert-luminosity host galaxies with four additional metrics to the single-metric morphology analysis of Paper I, as well as new samples of inactive control galaxies. This expanded study of one of the largest and most complete, multi-wavelength samples of AGN detected at $0.6<z<2.4$ in the GOODS-South and North fields, confirms that mergers surprisingly play a significant role in obscured, sub-quasar AGN host galaxies. Additionally, the pattern of morphological disturbances observed amongst the X-ray- and mid-IR-selected AGN suggests that these represent different phases of AGN evolution tied to a major-merger timeline, as opposed to distinct populations of AGN. These results indicate that mergers are important in triggering sub-quasar AGN at these redshifts.

Clarissa R. Do Ó, Jaehan Bae, Quinn M. Konopacky, Jayke S. Nguyen, Patrick Diamond, Krzysztof Goździewski, Dawid Jankowski

Direct imaging has revealed exoplanet systems hosting multiple wide-orbit Super-Jupiters, where planet-planet interactions can shape their long-term dynamical evolution. These strong perturbations may lead to orbital instability, raising questions about the long-term survival of such systems. Shortly after formation, planet-disk interactions can shepherd planets into mean-motion resonances, which may promote long-term stability as seen in HR 8799. However, early-stage processes such as disk photoevaporation and viscosity can influence these outcomes. The $\sim$5 Myr-old PDS 70 system offers a unique laboratory to investigate these processes: its two massive ($>$4 $M_{Jup}$), wide-orbit ($>$20 AU) giants are still embedded in their natal disk. We perform 2D hydrodynamic simulations of the system, allowing the disk to disperse via photoevaporation. Once the disk dissipates, we continue to track the planets' orbital evolution over Gyr timescales using N-body simulations. We find that the system is likely to remain stable for $>$ 1 Gyr. To assess the importance of disk-driven evolution, we compare these results with disk-free N-body simulations using orbital parameters constrained by orbit fits that include recent relative astrometry and radial velocities from the literature. In this case, we find that only $\lesssim 4\%$ of posterior is stable for 100 Myr, highlighting the importance of considering disk-driven evolution for long-term dynamics stability of exoplanetary systems. We also simulate two three-planet configurations including the proposed inner candidate "PDS 70 d", finding that a higher photoevaporation leads the system to become unstable in $<$ 10 Myr.

A key measure of gravity is the relation between the Weyl potential $\Psi+\Phi$ and the matter overdensity $\delta_m$, capsulized as an effective gravitational constant $G_{\rm light}$ for light motion. Its value, together with the possible spatial and temporal variation, is essential in probing physics beyond Einstein gravity. However, the lack of an unbiased proxy of $\delta_m$ prohibits direct measurement of $G_{\rm light}$. We point out that the equivalence principle ensures the dispersion measure (DM) of localized fast radio bursts (FRBs) as a good proxy of $\delta_m$. We further propose a FRB-based method $F_G$ to directly measure $G_{\rm light}$, combining galaxy-DM of localized FRBs and galaxy-weak lensing cross-correlations. The measurement, with a conservative cut $k\leq 0.1h$/Mpc, can achieve a precision of $\lesssim 10\% \sqrt{10^5/N_{\rm FRB}}$ over 10 equal-width redshift bins at $z\lesssim 1$. The major systematic error, arising from the clustering bias of electrons traced by the FRB DM, is subdominant ($\sim 5\%$). It can be further mitigated to the $\lesssim 1\%$ level, based on the gastrophysics-agnostic behavior that the bias of total baryonic matter (ionized diffuse gas, stars, neutral hydrogen, etc) approaches unity at sufficiently large scales. Therefore, FRBs shed light on gravitational physics across spatial and temporal scales spanning over 20 orders of magnitude.

Blazars, particularly Flat Spectrum Radio Quasars (FSRQs), are well-known for their ability to accelerate a substantial population of electrons and positrons, as inferred from multiwavelength radiation observations. Therefore, these astrophysical objects are promising candidates for studying high-energy electron--positron interactions, such as the production of $W^{\pm}$ and $Z$ bosons. In this work, we explore the implications of electron--positron annihilation processes in the jet environments of FSRQs, focusing on the resonant production of electroweak bosons and their potential contribution to the diffuse neutrino flux. By modeling the electron distribution in the jet of the FSRQ 3C~279 during a flaring state, we calculate the reaction rates for $W^{\pm}$ and $Z$ bosons and estimate the resulting diffuse fluxes from the cosmological population of this http URL incorporate the FSRQ luminosity function and its redshift evolution to account for the population distribution across cosmic time, finding that the differential flux contribution exhibits a pronounced peak at redshift $z \sim 1$. While the expected fluxes remain well below the detection thresholds of current neutrino observatories such as IceCube, KM3NeT, or Baikal-GVD, the expected flux from the $Z$ boson production could account for approximately $10^{-3}$ of the total diffuse astrophysical neutrino flux. These results provide a theoretical benchmark for the role of Standard Model electroweak processes in extreme astrophysical environments and emphasize the interplay between particle physics and astrophysics, illustrating that even rare high-energy interactions can leave a subtle but quantifiable imprint on the diffuse astrophysical neutrinos.

Kuldeep J. Purohit, Jitesh R. Bhatt, Subhendra Mohanty, Prashant K. Mehta

We perform a linear mode analysis of a uniformly distributed cloud of axion-like particles (ALPs) embedded in a magnetized intergalactic medium, in order to investigate the stability of axion stars under realistic astrophysical conditions. We find that when the frequency $\omega$ of transverse waves is much smaller than the collision frequency $\nu_c$ of the intergalactic plasma, the conversion of ALPs into photons occurs on timescales far longer than the age of the Universe, ensuring stability of the star. In the opposite regime, $\omega \gg \nu_c$, significant axion-to-photon conversion may occur if the condition $\tfrac{\beta^2}{m_a^2-\omega_p^2} < 1$ is satisfied, where $\beta$ depends on the ALP--photon coupling and the magnetic field, $m_a$ is the ALP mass, and $\omega_p$ is the plasma frequency. We have calculated up to second order in perturbations to compute the effect of an ALP star. Since the calculated value of parameter $\beta ^2$ is extremely small in comparison with $\omega^2_p$, we argue that the direct detection of an axion star is highly unlikely in experiments like NCLE. However, since the calculated $\beta$ is extremely small compared to $\omega_p$, this requires an unrealistically fine-tuned coincidence between $m_a$ and $\omega_p$. As a consequence we argue that that detection of Our results therefore suggest that axion stars remain stable in typical intergalactic environments, though extreme magnetic fields (e.g.\ near magnetars) may lead to different outcomes.

Argus is a high-performance Python package for detecting and characterising nanohertz gravitational waves in pulsar timing array data. The package provides a complete Bayesian inference framework based on state-space models, using Kalman filtering for efficient likelihood evaluation. Argus leverages JAX for just-in-time compilation, GPU acceleration, and automatic differentiation, facilitating rapid Bayesian inference with gradient-based samplers. The state-space approach provides a computationally efficient alternative to traditional frequency-domain methods, offering linear scaling with the number of pulse times-of-arrival, and natural handling of non-stationary processes.

The correlations between the positions and shapes of galaxies, i.e. intrinsic alignments, have been measured in many observational studies and hydrodynamical simulations. The alignments of disk galaxies in hydrodynamical simulations have been measured to be positive, null and negative with varying methodologies, samples and hydrodynamical simulations. This work compares the correlations of disks and ellipticals around all galaxies and disks around ellipticals at $z=0$ and $z=1$ for simple and reduced shapes in TNG300, Horizon-AGN and EAGLE for multiple morphological definitions in a consistent way. All types of signals are positive and robust in TNG300 and EAGLE and positive or null in Horizon-AGN, except for the disks around ellipticals correlation for reduced shapes at $z=1$ when defined by $|v/\sigma|$, which is negative. A re-weighting of the ellipticals around galaxies signals in TNG300, according to the underlying stellar mass distributions of the samples, highlights the importance of the influence of (sub-grid) physics at these non-linear scales.

G. Kaur, M. Bilicki, S. Bellstedt, E. Tempel, W. A. Hellwing, I. Baldry, B. Bandi, S. Barsanti, S. Driver, N. Guerra-Varas, B. Holwerda, C. Lagos, J. Loveday, A. Robotham

The Wide-Area VISTA Extragalactic Survey (WAVES) on the 4-metre Multi-Object Spectroscopic Telescope (4MOST) includes two flux-limited subsurveys with very high (95\%) completeness requirements: Wide over $\sim\!1200$ deg$^2$ and Deep over $\sim\!65$ deg$^2$. Both are $Z$-band selected, respectively as $Z<21.1$ and $Z<21.25$ mag, and additionally redshift-limited, while the true redshifts are not known a priori but will be only measured by 4MOST. Here, we present a classification-based method to select the targets for WAVES-Wide. Rather than estimating individual redshifts for the input photometric objects, we assign probabilities of them being below $z=0.2$, the redshift limit of the subsurvey. This is done with the supervised machine learning approach of eXtreme Gradient Boosting (XGB), trained on a comprehensive spectroscopic sample overlapping with WAVES fields. Our feature space is composed of nine VST+VISTA magnitudes from $u$ to $K_s$ and all the possible colors, but most relevant for the classification are the $g$-band and the $u-g$, $g-r$ and $J-K_s$ colors. We check the performance of our classifier both for the fiducial WAVES-Wide limits, as well as for a range of neighboring redshift and magnitude thresholds, consistently finding purity and completeness at the level of 94-95\%. We note, however, that this performance deteriorates for sources close to the selection limits, due to deficiencies of the current spectroscopic training sample and the decreasing signal-to-noise of the photometry. We apply the classifier trained on the full spectroscopic sample to 14 million photometric galaxies from the WAVES input catalog, which have all 9 bands measured. Our work demonstrates that a machine-learning classifier could be used to select a flux- and redshift-limited sample from deep photometric data.

I report on the discovery of 34 new quasi-periodic oscillations (QPOs) in the prompt light curves of long gamma-ray bursts (GRBs) from the Swift/BAT catalog: with one or more constant leading periods, as well as several chirping signals. This is the largest homogenously identified sample or GRB QPOs to date. The presence of QPOs suggests the existence of characteristic time scales that at least in some GRBs might be related to the dynamical properties of plasma trajectories in the accretion disks powering the relativistic jets. Several scenarios for their origin were examined. We identify non-planar orbits around Kerr black holes, the Lense-Thirring effect, and shock oscillations as plausible mechanisms for the QPO generation.

The rapid growth of large-scale radio surveys, generating over 100 petabytes of data annually, has created a pressing need for automated data analysis methods. Recent research has explored the application of machine learning techniques to address the challenges associated with detecting and classifying radio galaxies, as well as discovering peculiar radio sources. This paper provides an overview of our investigations with the Evolutionary Map of the Universe (EMU) survey, detailing the methodologies employed-including supervised, unsupervised, self-supervised, and weakly supervised learning approaches -- and their implications for ongoing and future radio astronomical surveys.

Gan Luo, Arshia M. Jacob, Marco Padovani, Daniele Galli, Ana López-Sepulcre, Ningyu Tang, Di Li, Jing Zhou, Pei Zuo

Methylidyne (CH) has long been considered a reliable tracer of molecular gas in the low-to-intermediate extinction range. Although extended CH 3.3 GHz emission is commonly observed in diffuse and translucent clouds, observations in cold, dense clumps are rare. In this work, we conducted high-sensitivity CH observations toward 27 PGCCs with the Arecibo 305m telescope. Toward each source, the CH data were analyzed in conjunction with $^{13}$CO (1--0), HINSA, and H$_2$ column densities. Our results revealed ubiquitous subsonic velocity dispersions of CH, in contrast to $^{13}$CO, which is predominantly supersonic. The findings suggest that subsonic CH emissions may trace dense, low-turbulent gas structures in PGCCs. To investigate environmental effects, particularly the cosmic-ray ionization rate (CRIR), we estimated CRIR upper limits from HINSA, yielding values from $(8.1\pm4.7)\times10^{-18}$ to $(2.0\pm0.8)\times10^{-16}$ s$^{-1}$ ($N_{H_2}$ from $(1.7\pm0.2)\times10^{21}$ to $(3.6\pm0.4)\times10^{22}$~cm$^{-2}$). This result favors theoretical predictions of a cosmic-ray attenuation model, in which the interstellar spectra of low-energy CR protons and electrons match {\it Voyager} measurements, although alternative models cannot yet be ruled out. The abundance of CH decreases with increasing column density, while showing a positive dependence on the CRIR, which requires atomic oxygen not heavily depleted to dominate CH destruction in PGCCs. By fitting the abundance of CH with an analytic formula, we place constraints on atomic O abundance ($2.4\pm0.4\times10^{-4}$ with respect to total H) and C$^+$ abundance ($7.4\pm0.7\times10^{13}\zeta_2/n_{\rm H_2}$). These findings indicate that CH formation is closely linked to the C$^+$ abundance, regulated by cosmic-ray ionization, while other processes, such as turbulent diffusive transport, might also contribute a non-negligible effect.

We measure galaxy structural properties and colour gradients using HST images to trace the evolution of galaxy components. We jointly fit 3D bulge and disk models to 2505 galaxies in GOODS-South across seven bands (bvizYJH) to IAB = 25.5, accounting for different component ellipticities and inclination-dependent dust extinction. Extinction strongly affects structural parameters and colour gradients in ~26% of the sample - primarily edge-on galaxies with central obscuration (B-band face-on optical depth tau ~ 4) that reveal clear bulge components in the near-infrared. Despite irregular morphologies, the model captures observed colour gradients well. Bulges at z ~ 1 differ markedly from z ~ 0, with typical Sersic index n ~ 1.0 and bulge-to-disc size ratio Re/hd ~ 0.15, suggesting most galaxies host pseudo-bulges formed via secular evolution. Galaxy ellipticity correlates strongly with disk scale-length and absolute magnitude, partly driven by dust extinction variations. We trace bulge and disk evolution from z ~ 0 to z ~ 2.5: bulges are redder than disks (observed-frame) at z < 1.4, but colours converge at higher redshifts and fainter magnitudes. Redder galaxies show redder cores relative to their outskirts, and brighter galaxies have redder cores.

Jiwon Shin (1), C. Y. Hui (1), Sangin Kim (1), Kwangmin Oh (2), Ellis R. Owen (3 and 4) ((1) Chungnam National University, (2) Michigan State University, (3) RIKEN, (4) Osaka University)

Using 16 years of data collected by Fermi Large Area Telescope and 1523 days of survey data from High Altitude Water Cherenkov (HAWC) Observatory, we discovered the long-sought second GeV-TeV connection towards the globular cluster (GC) UKS 1 (Shin et al. 2025). Gamma-ray spectroscopy suggests that the GeV emission can be attributed to both the pulsar magnetosphere and inverse Compton scattering (ICS) by the pulsar wind. In particular, the TeV peak is displaced from the cluster center by several tidal radii in the trailing direction of the proper motion of UKS 1. This alignment supports a scenario in which relativistic leptons, likely driven by a millisecond pulsar population, produce very-high-energy (VHE) gamma-rays via ICS within a bow shock tail. Our findings not only highlights GCs as potential sources of VHE gamma-rays, but also offers a rare opportunity to probe cosmic ray transport in the Milky Way by studying particle propagation and anisotropic gamma-ray production associated with the extended, offset TeV feature of UKS 1.

Judhajeet Basu, G.C. Anupama, Jan-Uwe Ness, Kulinder Pal Singh, Sudhanshu Barway, Shatakshi Chamoli

We report on UV and X-ray observations of the 2024 eruption of the recurrent nova LMCN 1968-12a, a rapidly recurring extragalactic system with a $\sim$4.3 year recurrence period and a massive white dwarf (WD). The eruption was discovered on 2024 August 1.8 by \textit{Swift}, and subsequently monitored using \textit{AstroSat}'s UVIT and SXT, along with Swift/UVOT and XRT. The multi-wavelength light curves reveal a rapid UV-optical decline, followed by a plateau phase exhibiting 1.26-day modulations consistent with the orbital period. The Supersoft (SSS) X-ray emission, that emerged by day 5, exhibited a double peak, suggesting variable obscuration that could be due to an inhomogeneous nova ejecta or due to a nova super-remnant along the line of sight. Time-resolved X-ray spectroscopy shows a blackbody component with T $\approx 10^6$ K. The SEDs obtained concurrently in the UV, peaking at T $\approx$ 20,000 K and with a source radius $\sim$2-3 R$_\odot$, are inconsistent with emission from the secondary star or nova photosphere alone. Instead, the UV emission is attributed to an irradiated accretion disk that survived the eruption. The persistent UV plateau and its temperature suggest that the accretion disk was not completely disrupted and resumed activity within days, consistent with recent findings in other rapidly recurring novae such as U~Sco and M31N~2008-12a.

In the era of mega-constellations, the need for accurate and publicly available information has become fundamental for satellite operators to guarantee the safety of spacecrafts and the Low Earth Orbit (LEO) space environment. This study critically evaluates the accuracy and reliability of publicly available ephemeris data for a LEO mega-constellation - Starlink. The goal of this work is twofold: (i) compare and analyze the quality of the data against high-precision numerical propagation. (ii) Leverage Physics-Informed Machine Learning to extract relevant satellite quantities, such as non-conservative forces, during the decay process. By analyzing two months of real orbital data for approximately 1500 Starlink satellites, we identify discrepancies between high precision numerical algorithms and the published ephemerides, recognizing the use of simplified dynamics at fixed thresholds, planned maneuvers, and limitations in uncertainty propagations. Furthermore, we compare data obtained from multiple sources to track and analyze deorbiting satellites over the same period. Empirically, we extract the acceleration profile of satellites during deorbiting and provide insights relating to the effects of non-conservative forces during reentry. For non-deorbiting satellites, the position Root Mean Square Error (RMSE) was approximately 300 m, while for deorbiting satellites it increased to about 600 m. Through this in-depth analysis, we highlight potential limitations in publicly available data for accurate and robust Space Situational Awareness (SSA), and importantly, we propose a data-driven model of satellite decay in mega-constellations.

Marion Zannese, Jacques Le Bourlot, Evelyne Roueff, Emeric Bron, Franck Le Petit, Dries Van De Putte, Maryvonne Gerin, Naslim Neelamkodan, Javier R. Goicoechea, John Black, Ryan Chown, Ameek Sidhu, Emilie Habart, Els Peeters, Olivier Berné

The James Webb Space Telescope enabled the first detection of several rovibrational emission lines of HD in the Orion Bar, a prototypical photodissociation region. This provides an incentive to examine the physics of HD in dense and strong PDRs. Using the latest data available on HD excitation by collisional, radiative and chemical processes, our goal is to unveil HD formation and excitation processes in PDRs by comparing our state-of-the-art PDR model with observations made in the Orion Bar and discuss if and how HD can be used as a complementary tracer of physical parameters in the emitting region. We compute detailed PDR models, using an upgraded version of the Meudon PDR code, which are compared to NIRSpec data using excitation diagrams and synthetic emission spectra. The models predict that HD is mainly produced in the gas phase via the reaction D + H2 = H + HD at the front edge of the PDR and that the D/HD transition is located slightly closer to the edge than the H/H2 transition. Rovibrational levels are excited by UV pumping. In the observations, HD rovibrational emission is detected close to the H/H2 dissociation fronts of the Orion Bar and peaks where vibrationally excited H2 peaks, rather than at the maximum emission of pure rotational H2 levels. We derive an excitation temperature around Tex ~ 480 - 710 K. Due to high continuum in the Orion Bar, fringes lead to high noise levels beyond 15 $\mu$m, no pure rotational lines of HD are detected. The comparison to PDR models shows that a range of thermal pressure P = (3-9)x10$^7$ K cm$^{-3}$ with no strong constraints on the intensity of the UV field are compatible with HD observations. This range of pressure is compatible with previous estimates from H2 observations with JWST. This is the first time that observations of HD emission lines in the near-infrared are used to put constraints on the thermal pressure in the PDR.

Three ring systems have been discovered to date around small irregular objects of the solar system (Chariklo, Haumea and Quaoar). For the three bodies, material is observed near the second-order 1/3 Spin-Orbit Resonance (SOR) with the central object, and in the case of Quaoar, a ring is also observed near the second-order resonance 5/7 SOR. This suggests that second-order SORs may play a central role in ring confinement. This paper aims at better understanding this role from a theoretical point of view. It also provides a basis to better interpret the results obtained from N-body simulations and presented in a companion paper. A Hamiltonian approach yields the topological structure of phase portraits for SORs of orders from one to five. Two cases of non-axisymmetric potentials are examined: a triaxial ellipsoid characterized by an elongation parameter C22 and a body with mass anomaly mu, a dimensionless parameter that measures the dipole component of the body's gravitational field. The estimated triaxial shape of Chariklo shows that its corotation points are marginally unstable, those of Haumea are largely unstable, while those of Quaoar are safely stable. The topologies of the phase portraits show that only first- (aka Lindblad) and second-order SORs can significantly perturb a dissipative collisional ring. We calculate the widths, the maximum eccentricities and excitation time scales associated with first- and second-order SORs, as a function of C22 and mu. Applications to Chariklo, Haumea and Quaoar using mu ~ 0.001 show that the first- and second-order SORs caused by their triaxial shapes excite large (>~ 0.1) orbital eccentricities on the particles, making the regions inside the 1/2 SOR inhospitable for rings. Conversely, the 1/3 and 5/7 SORs caused by mass anomalies excite moderate eccentricities (<~ 0.01), and are thus a more favorable place for the presence of a ring.

Sven Buder, Tobias Buck, Ása Skúladóttir, Melissa Ness, Madeleine McKenzie, Stephanie Monty

Understanding how the Milky Way's present-day structure was shaped by past major mergers is a key goal of Galactic archaeology. The chemical and dynamical structure of the Galaxy retains the imprint of such events, including a major accretion episode around 8-10 Gyr ago. Recent findings suggest that present-day orbital energy correlates with stellar chemistry and birth location within the merging progenitor galaxy. Using a high-resolution NIHAO-UHD cosmological zoom-in simulation of a Milky Way analogue, we trace the birth positions, ages, and present-day orbits of stars accreted in its last major merger. We show that stars born in the progenitor's core are more tightly bound to the Milky Way and chemically enriched, while those from the outskirts are less bound and metal-poor. This supports the scenario proposed by Skúladóttir et al. (2025) using idealised simulations, now demonstrated in a cosmological context. The preserved chemodynamical memory is also evident in elemental planes such as [Al/Fe] vs. [Mg/Mn], reflecting gradients in star formation efficiency. The median abundance trends for different orbital energies are broadly similar, though less pronounced and more non-linear than the linear relations reported by Skúladóttir et al. (2025), with relatively good agreement at the highest metallicities. Our results reveal that spatial and temporal memory is retained across the merger -- connecting birth locations to present-day properties like a golden thread. We demonstrate that selection methods in integrals-of-motion space systematically bias progenitor reconstructions by missing its enriched core, and we outline strategies to obtain a more representative census of accreted stars.

Namrata Roy, Alaina Henry, Tucker Jones, Ivana Barisic, Ryan L. Sanders, Kevin Bundy, Matthew A. Malkan, Themiya Nanayakkara, Karl Glazebrook, Timothy Heckman, Juan M. Espejo Salcedo, Xin Wang, Danail Obreschkow, Tommaso Treu

We present spatially resolved rest-optical spectroscopy of 38 star-forming galaxies at 0.5 < z < 1.7 from the JWST/NIRSpec MSA-3D survey, which uses slit-stepping to build IFU-like datacubes at 0.1'' resolution. We map emission-line morphology, excitation, and kinematics of the warm ionized gas using [N II]/H$\alpha$, [S II]/H$\alpha$, and [O III]/H$\beta$. Relative to z$\sim$0 galaxies at fixed stellar mass, our sources show systematically lower [N II]/H$\alpha$ and [S II]/H$\alpha$ and elevated [O III]/H$\beta$, consistent with harder radiation fields and lower metallicities. Radially, [O III]/H$\beta$ profiles are typically flat or mildly positive, whereas [N II]/H$\alpha$ also remains flat or declines outward, mirroring metallicity trends. On kpc scales, we find a strong positive correlation between [N II]/H$\alpha$ and velocity dispersion ($\sigma$), linking local excitation to turbulent or shock-driven kinematics. Six galaxies ($\sim$ 16% of the sample) host spatially localized regions with elevated [N II]/H$\alpha$, high EW(H$\alpha$), and V_RMS = $\sqrt{V^2 + \sigma^2} > 200$ km/s, indicative of weak AGN activity, shocks, or outflows. For these candidates we infer modest warm-ionized outflow rates of 1-4 Msun/yr and kinetic powers $\sim$ 0.1-1% of the AGN bolometric luminosity (from central [O III] or H$\alpha$). These values place our sample at the low-energy tail of known AGN-driven outflows yet in continuity with $\dot{M_{out}}-L_{AGN}$ scaling relations across 0 < z < 6. A completeness assessment shows MSA-3D is sensitive to AGN with $L_{AGN} \geq 10^{43}$ erg/s, underscoring both the promise and current limitations of detecting weak AGN activity in distant galaxies with resolved spectroscopy.

Lucie E. Rowland, Kasper E. Heintz, Hiddo Algera, Mauro Stefanon, Jacqueline Hodge, Rychard Bouwens, Manuel Aravena, Elisabete da Cunha, Pratika Dayal, Andrea Ferrara, Rebecca Fisher, Valentino González, Hanae Inami, Olena Komarova, Ilse de Looze, Themiya Nanayakkara, Katherine Ormerod, Andrea Pallottini, Clara L. Pollock, Renske Smit, Paul van der Werf, Joris Witstok

Neutral gas in galaxies during the Epoch of Reionisation regulates star formation, dust growth, and the escape of ionising photons, making it a key ingredient in understanding both galaxy assembly and reionisation. Yet, direct constraints on the HI content of galaxies at z>6 have been scarce. With JWST, Ly$\alpha$ damping wings in galaxy spectra can now provide a direct probe of this neutral component. We analyse JWST/NIRSpec prism spectra of 12 UV-luminous galaxies from the REBELS-IFU program at z~6.5-7.7, deriving HI column densities by modelling Ly$\alpha$ damping wings. Significant damped Ly$\alpha$ absorption is detected in eight galaxies, with $N_{\mathrm{HI}}\gtrsim10^{21}$ cm$^{-2}$. We use the column densities and sizes derived for these sources to estimate their HI mass and compare with $L_{\mathrm{[CII]}}$-$M_{\mathrm{HI}}$ calibrations. The resulting HI masses show a tentative correlation with those inferred from [CII], although the [CII]-based estimates are systematically larger, suggesting that the HI reservoirs may extend beyond the [CII]-emitting gas. We also combine the DLA-based measurements with FIR-derived dust-to-gas ratios, dust attenuation, and gas-phase metallicities. No correlation is found between DLA-based and FIR-based dust-to-gas ratios, but combining the REBELS-IFU sample with literature samples at lower metallicities reveals a strong correlation between $A_{\mathrm{V}}/N_{\mathrm{HI}}$ and metallicity. These findings suggest that by $z\sim7$ massive galaxies can already host substantial, enriched reservoirs of neutral gas and dust, consistent with $A_{\mathrm{V}}$/$N_{\mathrm{HI}}$-metallicity trends at lower redshift. At the highest redshifts ($z>8$), however, we see tentative evidence for systematically lower $A_{\mathrm{V}}$/$N_{\mathrm{HI}}$ at fixed metallicity, which may point to pristine gas accretion or more efficient dust destruction/expulsion.

M.P. Snelders, J.W.T. Hessels, J. Huang, N. Sridhar, B. Marcote, A.M. Moroianu, O.S. Ould-Boukattine, F. Kirsten, S. Bhandari, D.M. Hewitt, D. Pelliciari, L. Rhodes, R. Anna-Thomas, U. Bach, E.K. Bempong-Manful, V. Bezrukovs, J.D. Bray, S. Buttaccio, I. Cognard, A. Corongiu, R. Feiler, M.P. Gawroński, M. Giroletti, L. Guillemot, R. Karuppusamy, M. Lindqvist, K. Nimmo, A. Possenti, W. Puchalska, D. Williams-Baldwin

FRB 20121102A is the original repeating fast radio burst (FRB) source and also the first to be localised to milliarcsecond precision using very-long-baseline interferometry (VLBI). It has been active for over 13 years and resides in an extreme magneto-ionic environment in a dwarf host galaxy at a distance of ~1 Gpc. In this work, we use the European VLBI Network (EVN) to (re-)localise FRB 20121102A and its associated persistent radio source (PRS). We confirm that the two are co-located -- improving on previous results by a factor of ~4 and constraining the FRB and PRS co-location to ~12 pc transverse offset. Over a decade, the PRS luminosity on milliarcsecond scales remains consistent with measurements on larger angular scales, showing that the PRS is still compact. We also present the detection of 18 bursts with the Nancay Radio Telescope (NRT) as part of our ÉCLAT monitoring program. These bursts, together with previously published results, show that the observed dispersion measure (DM) of FRB 20121102A has dropped by ~25 pc/cc in the past five years, highlighting a fractional decrease in the local DM contribution of >15%. We discuss potential physical scenarios and highlight possible future observations that will help reveal the nature of FRB 20121102A, which is one of only a few known FRBs with a luminous PRS.

For any elliptical potential with an external parallel shear, Witt has proven that the gravitational center lies on a rectangular hyperbola derived from the image positions of a single quadruply lensed object. Moreover, it is predicted that for an isothermal elliptical potential the source position both lies on Witt's Hyperbola and coincides with the center of Wynne's Ellipse (fitted through the four images). Thus, by fitting Witt's Hyperbolae to several quartets of images - ten are known in Abell 1689 - the points of intersection provide an estimate for the center for the assumed isothermal elliptical potential. We introduce a new figure of merit defined by the offset of the center of Wynne's Ellipse from Witt's Hyperbola. This offset quantifies deviations from an ideal elliptical isothermal potential and serves as a discriminant to exclude poorly fitted quadruples and assign greater weight to intersections of hyperbolae of better fitting systems. Applying the method to 10 quads (after excluding 7 poorly fitted quads) in Abell 1689, we find the potential is centered within 11" of the BCG, X-ray center, flexion-based center and the center found from a total strong lensing analysis. The Wynne-Witt framework thus delivers a fast, analytic, and self-consistency-checked estimator for centers in clusters with multiple quads.

Laurine Martinien, Gaspard Duchêne, François Ménard, Karl R. Stapelfeldt, Ryo Tazaki, Jennifer B. Bergner, Emmanuel Dartois, Jennifer A. Noble, William Thompson

The James Webb Space Telescope provides unprecedented information to study ices in protoplanetary disks. However, the saturation of ice bands in highly inclined disks hinders the measurement of ice abundances using classical spectroscopy. This is unfortunate as the presence and more importantly abundance of ices plays a key role in, e.g., the evolution of dust (because it modifies the sticking properties) and the composition of planetesimals and exoplanetary atmospheres. To overcome this issue and quantify the ice abundance within disks, we introduce a new method based on measuring the changes in the apparent disk thickness as a function of wavelength, which is directly and quantitatively related to the grain opacity. Specifically, we expect i) that the increased opacity within ice bands should result in a thicker disk than in the adjacent continuum, and ii) the thickness variations to be proportional to the abundance of ice. We extracted the disk thickness in model images of edge-on disks containing different abundances of water ice, as well as in James Webb Space Telescope spectral imaging of four edge-on disks. For both models and observations, the disk thickness decreases toward longer wavelengths except across the positions of ice absorption features where the thickness is enhanced across the band. In the model images, we demonstrate that this effect increases with ice abundance without any hint of saturation. This definitely demonstrates the presence of the ice species within each disk and confirms our expectation that this method can be applied to estimate ice abundances. Thanks to this method, it will thus be possible to constrain the ice abundance in highly inclined disks with disks model fitting. Unlike spectroscopic analysis, this method is not subject to saturation and should therefore be more robust and applicable to all disks for which the two surfaces can be resolved.

Recent DESI baryon acoustic oscillation (BAO) measurements, combined with Planck cosmic microwave background (CMB) data and DESY5 type Ia supernova (SN) data, indicate a significant deviation from $\Lambda$CDM, which seems to suggest that this deviation can be explained by an interaction between dark energy and dark matter. In this work, we perform a comprehensive analysis by utilizing the latest DESI DR2 BAO data in conjunction with CMB data from ACT, SPT, Planck, and WMAP, along with SN data from PantheonPlus and DESY5. We consider four interacting dark energy (IDE) models with different forms of the interaction term $Q$. Our analysis indicates that CMB experiments other than Planck enhance the evidence for an interaction in the IDE models with $Q \propto \rho_{\rm de}$. In particular, when using the SPT+DESI+DESY5 data, the IDE model with $Q = \beta H_0 \rho_{\rm de}$ gives $\beta = -0.4170 \pm 0.1220$, with a deviation from zero reaching $3.4\sigma$ level. When replacing DESY5 with PantheonPlus, this deviation weakens to $2.1\sigma$ level, but remains relatively significant. Furthermore, the Bayes factors of the IDE model with $Q = \beta H_0 \rho_{\rm de}$ are positive in all cases, providing a moderate-to-strong preference over $\Lambda$CDM. Overall, our comprehensive analysis clearly suggests that the IDE models with $Q \propto \rho_{\rm de}$ (especially, $Q = \beta H_0 \rho_{\rm de}$) provide strong evidence supporting the existence of interaction and are more preferred by the current cosmological data.

Vedant Dhruv, Ben Prather, Mani Chandra, Abhishek V. Joshi, Charles F. Gammie

The black holes in the Event Horizon Telescope sources Messier 87* and Sagittarius A* (SgrA*) are embedded in a hot, collisionless plasma that is fully described in kinetic theory yet is usually modeled as an ideal, magnetized fluid. In this Letter, we present results from a new set of weakly collisional fluid simulations in which leading order kinetic effects are modeled as viscosity and heat conduction. Consistent with earlier, lower-resolution studies, we find that overall flow dynamics remain very similar between ideal and non-ideal models. For the first time, we synthesize images and spectra of SgrA* from weakly collisional models -- assuming an isotropic, thermal population of electrons -- and find that these remain largely indistinguishable from ideal fluid predictions. However, most weakly collisional models exhibit lower light curve variability, with all magnetically dominated models showing a small but systematic decrease in variability.

Romain A. Meyer, Feige Wang, Koki Kakiichi, Gabe Brammer, Jackie Champagne, Katharina Jurk, Zihao Li, Zijian Li, Marat Musin, Sindhu Satyavolu, Jan-Torge Schindler, Marko Shuntov, Yi Xu, Siwei Zou, Fuyan Bian, Caitlin Casey, Eiichi Egami, Xiaohui Fan, Danyang Jiang, Nicolas Laporte, Weizhe Liu, Pascal Oesch, Lidia Tasca, Jinyi Yang, Zijian Zhang, Hollis Akins, Zheng Cai, Dave A. Coulter, Jiamu Huang, Mingyu Li, Weizhe Liu, Yongming Liang, Xiangyu Jin, Jeyhan Kartaltepe, Jasleen Matharu, Maria Pudoka, Wei-Leong Tee, Callum Witten, Haowen Zhang, Yongda Zhu

We present a spectroscopically-selected [OIII]+Hb emitters catalogue at 6.75<z<9.05 and the resulting [OIII] 5008 ÅLuminosity Function (LF) in the COSMOS field. We leverage the 0.3 deg$^{2}$ covered to date by COSMOS-3D using NIRCam/WFSS F444W (90% of the survey) to perform the largest spectroscopic search for [OIII] emitters at 6.75<z<9.05. We present our catalogue of 237 [OIII] emitters and their associated completeness function. The inferred constraints on the [OIII] LF enable us to characterise the knee of the [OIII] LF, resulting in improved [OIII] LF constraints at z~7,8. Notably, we find evidence for an accelerated decline of the [OIII] luminosity density between z~7 and z~8, which could be expected if the metallicity of [OIII] emitters, as well as the cosmic star-formation rate density, is declining at these redshifts. We find that theoretical models that reproduce the z~7,8 [OIII] LF do not reproduce well the [OIII] equivalent width distribution, pointing to potential challenges in the modelling of[OIII] and other nebular lines in the early Universe. Finally, we provide the first constraints on the cosmic variance of [OIII] emitters, estimating at 15% the relative uncertainty for the z~7,8 [OIII] LF in the 0.3 deg$^2$ field. This estimate is in good agreement with that inferred from clustering, and shows that the [OIII] LF derived from smaller extragalactic legacy fields is strongly affected by cosmic variance. Our results highlight the fundamental role that wide-area JWST slitless surveys play to map the galaxy large-scale structure down into the reionisation era, serving as a springboard for a variety of science cases.

In this work, we present an extensive review and detailed analysis of sunspot measurements, drawings, and engravings made by John Flamsteed and, mainly, by Philippe de La Hire during the Maunder minimum. All available information and contemporary knowledge about the sunspot nature are shown. The coordinates, areas, and numbers of sunspots and sunspot groups are reconstructed. Based on these observations, La Hire, Jean-Dominique Cassini, and his son Jacques Cassini regularly published results that shed light on the purpose of sunspot measurements and the scientific paradigm of that time. In particular, astronomers believed that sunspots were recurrent over decades. We compare the reconstructed time-latitude diagram with those obtained by Spoerer (Ueber die periodicitat der sonnenflecken seit dem Jahre 1618..., 1889) and Ribes and Nesme-Ribes (Astron. Astrophys. 276, 549, 1993). The sidereal differential rotation rate is estimated, and its latitudinal profile is reconstructed. We also evaluate the fraction of sunspot groups that obey or violate Joy's law.

The local Hubble flow offers a powerful laboratory to study the interplay between cosmic expansion and gravitational dynamics. On large scales, galaxy velocities follow Hubble's law, but within groups and clusters local gravitational effects introduce significant departures from linearity. Using the IllustrisTNG cosmological simulations, we investigate whether dark energy leaves detectable imprints on the local velocity-radius relation. We model the kinematics with extensions of the Lemaitre-Tolman framework and apply Bayesian inference to recover halo masses and the Hubble constant H0. The fits reveal systematic biases: halo masses are underestimated with a median ratio $M_{fit}/M_{true} = 0.95 \pm 0.28$, while the inferred Hubble constant clusters around $H_0 = 64 \pm 16 km/s/Mpc$, compared to the simulation input of 67.74. This corresponds to an average 25\% uncertainty in H0 recovery from the local flow method. While the mass and expansion rate can be constrained, different model variants whether including angular momentum, friction, or altered radial scaling-remain statistically indistinguishable. Our results highlight both the promise and the limitations of using local kinematics as a precision probe of dark energy.

Yue Wang, Chen-Wei Wang, Shaolin Xiong, Xiao Xiao, Yanqiu Zhang, Sheng-Lun Xie, Lin Lin, Yuan-Pei Yang, Haoxuan Guo, Ce Cai, Yue Huang, Cheng-Kui Li, Bing Li, Xiaobo Li, Jiacong Liu, Xiang Ma, Liming Song, Wen-Jun Tan, Ping Wang, Wang-Chen Xue, Shu-Xu Yi, Yun-Wei Yu, Zheng-Hang Yu, Jin-Peng Zhang, Peng Zhang, Shuang-Nan Zhang, Wen-Long Zhang, Zhen Zhang, Xiao-Yun Zhao, Chao Zheng, S.J. Zheng

Magnetar X-ray Burst (MXB) is usually composed of a single pulse or multiple pulses with rapid rise and brief duration mostly observed in hard X-ray (soft gamma-ray) band. Previous work studied the temporal behavior of some magnetar bursts and employed the Fast Rise Exponential Decay (FRED) model to fit pulses of MXB. However, whether there is other kind of pulse shape has not been explored. In this study, we systematically examined light curve of MXBs from SGR J1935+2154 detected by GECAM between 2021 and 2022. We find that there are different light curve morphologies. Especially, we discover a peculiar and new pattern, Exponential Rise and Cut-Off Decay (ERCOD), which is significantly different from FRED and could be well described by a mathematical function we proposed. We find that MXBs with ERCOD shape are generally longer in duration, brighter in the peak flux, and harder in spectrum. We note that the ERCOD shape is not unique to SGR J1935+2154 but also present in other magnetars. This new light curve pattern may imply a special burst and radiation mechanism of magnetar.

Alexandra Rochon, Étienne Artigau, Drew Weisserman, Lisa Dang, René Doyon, Charles Cadieux, Ryan Cloutier

We present the reanalysis of three 15 micron JWST/MIRI secondary eclipses of LHS 1140 c, a warm super-Earth (R$_{\rm{p}}$ = 1.272 R$_{\oplus}$) in a 3.78-day orbit around an M4.5 dwarf. We present a novel method for data reduction that leverages spatial derivatives of the point-spread function and compare it to widely used aperture photometry. Both methods yield eclipse depth consistent within 1 sigma of the values reported in the literature. We measure an eclipse depth of 271$^{+31}_{-30}$ ppm corresponding to a brightness temperature of $T_B=595^{+33}_{-34}$ K, consistent with a bare rock. The secondary eclipse occurs 4.1$\pm$0.8 minutes before the circular-orbit predicted time. We explore the implications of our results on the internal structure of LHS 1140 c, the orbital architecture of the system and the possibility of future observations with JWST. We find a core-mass fraction (CMF) informed by the stellar abundances of refractory elements of 0.34$\pm$0.11, inflated compared to the CMF from radius and mass measurements, suggesting the possible presence of bulk volatiles in the interior.

Guotang Wu, Xiaoli Yan, Zhike Xue, Jincheng Wang, Zhe Xu, Liheng Yang, Yian Zhou, Liping Yang, Xinsheng Zhang, Qifan Dong, Zongyin Wu

To better understand the characteristics, driving mechanisms, and potential heating contributions of chromospheric jets, we analyze two contrasting types: one originating from within the sunspot penumbra (inside jets), and the other originating from outside the penumbra (outside jets). Statistical analysis of 100 jets (50 inside jets and 50 outside jets) reveals that inside jets have a projected velocity range of 4--14~km\,s$^{-1}$, a length range of 1--4~Mm, a width range of 0.2--0.6~Mm, and a lifetime range of 135--450~s, with mean values of 7.90~km\,s$^{-1}$, 2.61~Mm, 0.41~Mm, and 260~s, respectively. About 52\% of inside jets are associated with brightenings in H$\alpha$ blue wing images, and some show high-temperature signatures, suggesting a connection with localized energy release. In contrast, outside jets have higher velocities (8--50~km\,s$^{-1}$, average 19.04~km\,s$^{-1}$), greater lengths (average 6.26~Mm, up to 27.27~Mm), slightly larger widths (average 0.46~Mm), and longer lifetimes (135--630~s, average 327~s). They typically originate from regions of opposite magnetic polarities and are associated with magnetic flux emergence and EUV brightenings. Some outside jets correspond to coronal jets with inverted Y-shaped structures and temperatures exceeding one million Kelvin. Our results suggest that both jet types are driven by magnetic reconnection occurring in distinct magnetic field configurations and contribute to chromospheric and coronal heating.

We present comprehensive spectral and timing results of 14 Chandra, 6 XMM-Newton and 19 Swift-XRT observations of the ultraluminous X-ray source NGC 4490 ULX-8, spanning from 2000 to 2024. We model the source spectra using absorbed power-law and absorbed multicolour disc blackbody models. The best-fit photon indices span 0.92-2.68, with typical uncertainties ranging from $\pm$0.1 to $\pm$1 depending on data quality. The inner disk temperature range from 0.97 to 1.69 keV, consistent with blackbody emission from an accretion disk. Our results reveal significant long-term variability in intrinsic X-ray source fluxes while the source remains relatively stable within individual observations. A Hardness-Intensity Diagram of the source shows no clear transition between hard and soft states, but an increase in brightness during two recent observations taken on 2022 December 1 and 2024 May 4. We find a positive correlation of X-ray luminosity and photon index that persists even when the hydrogen column density is tied across observations, suggesting a physical origin. The X-ray luminosity-inner disk temperature relation yields a weakly constrained slope owing to large temperature uncertainties, but a simpler fixed-slope test indicates consistency with a standard thin-disk. Using the derived disk parameters, we estimate the black hole mass to lie in the range of 16-75 $M_{\odot}$, under the assumption of a geometrically thin accretion flow, where the lower and upper bounds correspond to a Schwarzchild and a Kerr black hole respectively. Alternatively, we consider the scenario of ULX-8 hosting an accreting neutron star and estimate the corresponding magnetic field strength required to explain the observed properties.

Maja Lujan Niemeyer, Eiichiro Komatsu, José Luis Bernal, Chris Byrohl, Robin Ciardullo, Olivia Curtis, Daniel J. Farrow, Steven L. Finkelstein, Karl Gebhardt, Caryl Gronwall, Gary J. Hill, Matt J. Jarvis, Donghui Jeong, Erin Mentuch Cooper, Deeshani Mitra, Shiro Mukae, Julian B. Muñoz, Masami Ouchi, Shun Saito, Donald P. Schneider, Lutz Wisotzki

We present a measurement of the Lyman-$\alpha$ (Ly$\alpha$) intensity mapping power spectrum from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We measure the cross-power spectrum of the Ly$\alpha$ intensity and Ly$\alpha$-emitting galaxies (LAEs) in a redshift range of $1.9 < z < 3.5$. We calculate the intensity from HETDEX spectra that do not contain any detected LAEs above a signal-to-noise ratio of $5.5$. To produce a power spectrum model and its covariance matrix, we simulate the data using lognormal mocks for the LAE catalog and Ly$\alpha$ intensity in redshift space. The simulations include the HETDEX sensitivity, selection function, and mask. The measurements yield the product of the LAE bias, the intensity bias, the mean intensity of undetected sources, and the ratio of the actual and fiducial redshift-space distortion parameters, $b_\mathrm{g} b_I \langle I \rangle \bar{F}_{\rm RSD} / \bar{F}^{\rm fid}_{\rm RSD}= (6.7 \pm 3.1)$, $(11.7 \pm 1.4)$, and $(8.3 \pm 1.5) \times 10^{-22} \, \text{erg}\, \text{s}^{-1} \, \text{cm}^{-2} \, \text{arcsec}^{-2} \, \text{Å}^{-1}$ in three redshift bins centered at $\bar z=2.1$, 2.6, and 3.2, respectively. The results are reasonably consistent with cosmological hydrodynamical simulations that include Ly$\alpha$ radiative transfer. They are, however, significantly smaller than previous results from cross-correlations of quasars with Ly$\alpha$ intensity. These results demonstrate the statistical power of HETDEX for Ly$\alpha$ intensity mapping and pave the way for a more comprehensive analysis. They will also be useful for constraining models of Ly$\alpha$ emission from galaxies used in modern cosmological simulations of galaxy formation and evolution.

Davide Bevacqua, Danilo Marchesini, Paolo Saracco, Francesco La Barbera, Richard Pan, Sirio Belli, Gabriel Brammer, Guido De Marchi, Fabio R. Ditrani, Giovanna Giardino, Karl Glazebrook, Valentina La Torre, Jamie Lin, Adam Muzzin, Namrata Roy, Paola Santini, Benedetta Vulcani, Peter J. Watson, Xin Wang

We report the spectroscopic detection of neutral gas inflow into a massive ($M_* \simeq 4\times 10^{10} M_\odot$) quiescent galaxy observed at $z_{\rm{spec}} = 2.6576$ with JWST. From the redshifted absorption of the NaI doublet at $\lambda \lambda 5890, 5896 $ Ang, we estimate an inflow velocity $v=278^{+79}_{-79}$ km s$^{-1}$ and a column density $\log(N_{NaI}/\rm{cm^2}) = 13.02^{+0.03}_{-0.03}$. We derive the inflowing mass of the gas $M_{in} = 1.6^{+0.1}_{-0.1} \times 10^8 M_\odot$ and rate $\dot{M}_{in} = 19^{+6}_{-7} \, M_\odot \, \rm{yr}^{-1}$. The presence of several surrounding galaxies suggests that the galaxy may be accreting gas from nearby companions. However, we cannot confirm it with current data and the intergalactic medium or cosmic filaments are also viable sources of the inflowing gas. Despite the ongoing inflow, the galaxy remains quiescent, with an upper limit to the star formation rate of $0.2 \, M_\odot \, \rm{yr}^{-1}$. Moreover, its star formation history suggests that the galaxy has remained quiescent during the past $\sim1$ Gyr, with half of its stars formed by redshift $z_{50}=11^{+18}_{-3}$. We discuss that the inflow is not massive, dense, or long-lived enough to ignite significant star formation, or it is fueling low-level AGN activity instead. This is direct evidence that quiescent galaxies can accrete cold gas after their quenching while keeping their star formation subdued. Follow-up observations with JWST and ALMA will be needed to constraint the nature of the inflowing gas.

Xilu Wang, Amol V. Patwardhan, Yangming Lin, Junbo Zheng, Michael J. Cervia, Yanwen Deng, A. Baha Balantekin, Haining Li, Ian U. Roederer, Rebecca Surman

The astrophysical origin of the lanthanides is an open question in nuclear astrophysics. Besides the widely studied $s$, $i$, and $r$ processes in moderately-to-strongly neutron-rich environments, an intriguing alternative site for lanthanide production could in fact be robustly $\textit{proton-rich}$ matter outflows from core-collapse supernovae under specific conditions -- in particular, high-entropy winds with enhanced neutrino luminosity and fast dynamical timescales. In this environment, excess protons present after charged particle reactions have ceased can continue to be converted to neutrons by (anti-)neutrino interactions, producing a neutron capture reaction flow up to A~200. This scenario, christened the $\nu i$ process in a recent paper, has previously been discussed as a possibility. Here, we examine the prospects for $\nu i$ process through the lens of stellar abundance patterns, bolometric lightcurves, and galactic chemical evolution models, with a particular focus on hypernovae as candidate sites. We identify specific lanthanide signatures for which the $\nu i$ process can provide a credible alternative to $r$/$i$ processes.

The population of ultra-hot Jupiters (UHJs) provide unique opportunities to probe the extreme formation and evolutionary pathways in exoplanets. Owing to their very high temperatures and inflated atmospheres, UHJs are among the most favorable targets for both transmission and emission spectroscopy, enabling detailed characterization of their atmospheric properties. Here, we present a reanalysis of the JWST NIRSpec/G395H transmission spectra of the extreme ultra-hot Jupiter (EUHJ) WASP-178b, aimed at precisely characterizing its atmospheric composition. Our approach combines data reduction using two independent pipelines, lightcurve modeling with robust detrending techniques, and rigorous atmospheric retrievals. We report statistically significant detections of CO (7.24 $\sigma$) and CO$_2$ (7.22 $\sigma$), along with marginal evidence for C$_2$H$_2$ (1.34 $\sigma$), but no clear evidence for H$_2$O, suggesting depletion. From these retrieved abundances, we constrain the C/O ratio to a precise super-solar value of 0.954$\pm$0.033, consistent with an emerging trend in other UHJs. We also infer a very high atmospheric metallicity for a Jupiter-sized gas giant$\unicode{x2014}$11.44$_{-6.94}^{+12.54}$ $\times$solar$\unicode{x2014}$indicating unique atmospheric evolutions. These findings provide a critical benchmark for an extreme exoplanet atmosphere, offering a testbed for developing next-generation atmospheric evolution models and enabling comparative population-level studies across the UHJ population.

A. S. Bonomo, L. Naponiello, A. Sozzetti, S. Benatti, I. Carleo, K. Biazzo, P. E. Cubillos, M. Damasso, C. Di Maio, C. Dorn, N. Hara, D. Polychroni, M.-L. Steinmeyer, K. A. Collins, S. Desidera, X. Dumusque, A. F. Lanza, B. S. Safonov, C. Stockdale, D. Turrini, C. Ziegler, L. Affer, M. D'Arpa, V. Fardella, A. Harutyunyan, V. Lorenzi, L. Malavolta, L. Mancini, G. Mantovan, G. Micela, F. Murgas, D. Nardiello, I. Pagano, E. Pallé, M. Pedani, M. Pinamonti, M. Rainer, G. Scandariato, R. Spinelli, T. Zingales

Sub-Neptunes with planetary radii of $R_{p} \simeq 2-4 R_{\oplus}$ are the most common planets around solar-type stars in short-period ($P<100$ d) orbits. It is still unclear, however, what their most likely composition is, that is whether they are predominantly gas dwarfs or water worlds. The sub-Neptunes orbiting bright host stars are very valuable because they are suitable for atmospheric characterization, which can break the well-known degeneracy in planet composition from the planet bulk density, when combined with a precise and accurate mass measurement. Here we report on the characterization of the sub-Neptune TOI-5789 c, which transits in front of the bright ($V=7.3$ mag and $K_{s}=5.35$ mag) and magnetically inactive K1V dwarf HIP 99452 every 12.93 days, thanks to TESS photometry and 141 high-precision radial velocities obtained with the HARPS-N spectrograph. We find that its radius, mass, and bulk density are $R_{c}=2.86^{+0.18}_{-0.15} R_\oplus$, $M_{c}=5.00 \pm 0.50 M_\oplus$, and $\rho_{c}=1.16 \pm 0.23$ g cm$^{-3}$, and we show that TOI-5789 c is a promising target for atmospheric characterization with both JWST and, in the future, Ariel. By analyzing the HARPS-N radial velocities with different tools, we also detect three additional non-transiting planets, namely TOI-5789 b, d, and e, with orbital periods and minimum masses of $P_{b}=2.76$ d, $M_{b}\sin{i}=2.12 \pm 0.28 M_\oplus$, $P_{d}=29.6$ d, $M_{d}\sin{i}=4.29 \pm 0.68 M_\oplus$, and $P_{e}=63.0$ d, $M_{e}\sin{i}=11.61 \pm 0.97 M_\oplus$. The mutual orbital inclination between planets b and c must be higher than $\sim4$ deg, which points to a dynamically hot system. Nevertheless, from sensitivity studies based on both the HARPS-N and archival HIRES radial-velocity measurements, we can exclude that such high mutual inclinations are due to the perturbation by an outer gaseous giant planet.

C. Moutou, P. Petit, P. Charpentier, P. Cristofari, C. Baruteau, P. Thébault, L. Arnold, E. Artigau, A. Carmona, N.J. Cook, F. Debras, X. Delfosse, J.-F. Donati, L. Malo, M. Ould-Elhkim

One of the first exoplanet hosts discovered thirty years ago, the star 55 Cnc has been constantly observed ever since. It is now known to host at least five planets with orbital periods ranging from 17 hours to 15 years. It is also one of the most extreme metal rich stars in the neighbourhood and it has a low-mass secondary star. In this article, we present data obtained at the Canada-France-Hawai'i Telescope with the SPIRou spectropolarimeter on both components of the 55 Cnc stellar system. We revisit the long-period radial-velocity signals of 55 Cnc A, with a focus on the role of the magnetic cycle, and propose the existence of a sixth planet candidate, whose period falls close to that of the magnetic cycle, or half of it. The other massive outer planet has a revised period of 13.15 years and a minimum mass of 3.8 MJup. Although some uncertainty remains on these outer planets, the characterization of the four inner planets is very robust through the combination of many different data sets, and all signals are consistent in the nIR and optical domains. In addition, the magnetic topology of the solar-type primary component of the system is observed by SPIRou at the minimum of its activity cycle, characterized by an amplitude ten times smaller than observed during its maximum in 2017. For the low-mass component 55 Cnc B, we report the discovery of two exoplanets in the system, with a period of 6.799+-0.0014 and 33.75+-0.04 days and a minimum mass of 3.5+-0.8 and 5.3+-1.4 MEarth, respectively. The secondary magnetic field is very weak and the current data set does not allow its precise characterization, setting an upper limit of 10 G. The system 55 Cnc stands out as the sixth binary system with planetary systems around both components, and the first one with non equal-mass stellar components.

The Alcock-Paczynski (AP) parameter $F_{AP}$ is independent of the sound horizon $r_d$, making the Dark Energy Spectroscopic Instrument (DESI) baryon acoustic oscillation (BAO) AP measurements particularly well suited for cosmological applications. We propose a novel null test of cosmic curvature tailored to DESI BAO data that combines $F_{AP}$ with the ratios $D_V'/D_V$ or $D_M'/D_M$. This null test can also be performed using a joint dataset of DESI BAO and type Ia supernova (SNe Ia) observations. Additionally, we use the test to assess the internal consistency and mutual compatibility of these datasets. We find that the data are compatible. Although the results show that a spatially flat universe is inconsistent with the data at low redshift $z\lesssim 0.5$, we cannot draw the conclusion that the observational data prefers $\Omega_k\neq 0$ because there is no observational data in that region.

Sean M. O'Brien, Megan E. Schwamb, Christopher A. Watson, Louise D. Nielsen, Edward M. Bryant, Sarah L. Casewell, Matthew R. Burleigh, Lucy Fortson, Samuel Gill, Chris J. Lintott, Katlyn L. Hobbs, Ioannis Apergis, Daniel Bayliss, Jorge Fernández Fernández, Maximilian N. Günther, Faith Hawthorn, James S. Jenkins, Alicia Kendall, James McCormac, Ernst J. W. de Mooij, Toby Rodel, Suman Saha, Laura Trouille, Richard G. West, Peter J. Wheatley, Marius Constantin Agafitei, Deniz Rüzgar Apaydın, Elisabeth Baeten, Bruce Baller, Jeff Carabott, Sallyann Chesson, Sebastián Alejandro Freigeiro, Virgilio Gonano, Matthias Hanke, Pete Hermes, Avery Hildebrand, John S. Langley, See Min Lim, Leo Ryan McCarthy, Graham Mitchell, Ken O'Neill, Charles R. Pearson, Nolan Reket, Jeanne Riethmiller, Juergen Saeftel, Arttu Sainio, Charlie Steiner, Amanda Strickland, Christopher Tanner, Ivan A. Terentev, Ernest Jude P. Tiu, Sergey Y. Tumanov, Marciniak Urszula, Pia Vahlenkamp, Femke de Vroome, Paweł Wantuch, Timothy Woodruff

We report the identification and characterization of a new binary system composed of two near-equal mass M-dwarfs. The binary NGTS-EB-8 was identified as a planet candidate in data from the Next Generation Transit Survey (NGTS) by citizen scientists participating in the Planet Hunters NGTS project. High-resolution spectroscopic observations reveal the system to be a double-lined binary. By modeling the photometric and radial velocity observations, we determine an orbital period of 4.2 days and the masses and radii of both stars to be $M_A=0.250^{+0.005}_{-0.004}$ M$_{\odot}$, $M_B=0.208^{+0.005}_{-0.004}$ M$_{\odot}$, $R_A=0.255^{+0.004}_{-0.005}$ R$_{\odot}$, $R_B=0.233^{+0.006}_{-0.005}$ R$_{\odot}$. We detect Balmer line emission from at least one of the stars but no significant flare activity. We note that both components lie in the fully convective regime of low-mass stars ($<0.35$ M$_{\odot}$), therefore can be a valuable test for stellar evolutionary models. We demonstrate that the photometric observations, speckle imaging and initial radial velocity measurements were unable to identify the true nature of this system and highlight that high-resolution spectroscopic observations are crucial in determining whether systems such as this are in fact binaries.

Macarena G. del Valle-Espinosa, Matilde Mingozzi, Bethan James, Ruben Sanchez-Janssen, Juan Antonio Fernandez-Ontiveros, Ryan J. Rickards Vaught, Ricardo O. Amorin, Leslie Hunt, Alessandra Aloisi, Karla Z. Arellano-Cordova, Danielle A. Berg, John Chisholm, Matthew Hayes, Svea Hernandez, Alec Hirschauer, Logan Jones, Crystal L. Martin, Livia Vallini, Xinfeng Xu

Polycyclic Aromatic Hydrocarbons (PAHs) are key diagnostics of the physical conditions in the interstellar medium and are widely used to trace star formation in the mid-infrared (mid-IR). The relative strengths of mid-IR PAH emission features (e.g., 6.2, 7.7, 11.3 um) are sensitive to both the size and ionization state of the molecules and can be strongly influenced by the local radiation field. However, at low metallicities ( Z < 0.2 Zsun), detecting PAHs remains notoriously difficult, likely reflecting a combination of suppressed formation and enhanced destruction mechanisms. We present new JWST/MIRI MRS observations of the metal-poor (Z = 0.1 Zsun) dwarf galaxy CGCG 007-025. We confirm the tentative PAH detection previously reported from Spitzer data and, for the first time, identify a compact (approx. 50 pc) PAH-emitting region nearly co-spatial with the newly detected [NeV](I.P. = 97 eV) emission and the galaxy's most metal-poor, strongly star-forming region. The 11.3 um PAH feature is clearly detected, while no emission is found from the other typically brighter features, suggesting a PAH population dominated by large, neutral molecules resilient to hard ionizing fields. When compared with models, mid-IR line ratios involving [NeIII], [OIV], and [NeV] can only be reproduced by a combination of star formation and AGN ionization, with the latter contributing 4--8%. The [OIV] and [NeV] luminosities exceed what massive stars or shocks can produce, highlighting a puzzling scenario in line with recent JWST observations of similar galaxies. This work provides a crucial reference for studying the physical conditions of the dust and star formation in low-metallicity starburst regions, environments typical of the early universe.

The light elements beryllium (Be; $Z=4$) and boron (B; $Z=5$) are mainly produced by spallation reactions between cosmic rays and carbon (C; $Z=6$), nitrogen (N; $Z=7$), and oxygen (O; $Z=8$) nuclei. Only traces of Be or B would have been produced in the Big Bang, but there could be a contribution from the $\nu$-process in type II supernovae. Their abundances at very low metallicities have been debated in the literature, with the aim of understanding their origin. Our aim is to derive the boron abundance in a sample of metal-poor stars based for the first time on observations with the STIS spectrograph on board the Hubble Space Telescope, using clean B lines measured in space ultraviolet. We identified a measurable line of B I at 2089.6 A. In our sample of metal-poor warm stars, this line is practically free from blending lines, and for this reason the precision of the presently derived boron abundances is unprecedented. We find that in the interval -2.6<[Fe/H]<-1.0, the slope of the relation A(B) versus [Fe/H] is significantly larger than 1, and thus steeper than that obtained with Be abundances. As a consequence, we find in this interval of metallicity a B/Be ratio that slightly increases with [Fe/H]. Since at [Fe/H]=-1 the abundance of B is already close to the solar abundance, there should be a break in the B enrichment at a metallicity of about [Fe/H]=-1.

Primordial black holes (PBHs) are a compelling dark matter candidate and a unique probe of small-scale cosmological fluctuations. Their formation is usually attributed to large positive curvature perturbations, which collapse upon Hubble re-entry during radiation domination. In this work we investigate instead the role of negative curvature perturbations, corresponding to the growth of primordial void (PV) like regions. Using numerical relativity simulations, we show that sufficiently deep PV can undergo a nonlinear rebounce at the center, generating an effective overdensity that eventually collapses into a PBH. We determine the critical threshold for this process for a variety of equations of state, and demonstrate that the resulting black holes obey a scaling relation analogous to the standard overdensity case. These results establish primordial voids as a novel channel for PBH formation and highlight their potential impact on PBH abundances and cosmological signatures.

Francesco D'Eugenio, Erica J. Nelson, Daniel J. Eisenstein, Roberto Maiolino, Stefano Carniani, Jan Scholtz, Mirko Curti, Christopher N. A. Willmer, Andrew J. Bunker, Jakob M. Helton, Ignas Juodžbalis, Fengwu Sun, Sandro Tacchella, Santiago Arribas, Alex J. Cameron, Stéphane Charlot, Emma Curtis-Lake, Kevin Hainline, Benjamin D. Johnson, Brant Robertson, Christina C. Williams, Chris Willott, William M. Baker, Jacopo Chevallard, A. Lola Danhaive, Yuki Isobe, Xihan Ji, Zhiyuan Ji, Gareth C. Jones, Nimisha Kumari, Tobias J. Looser, Jianwei Lyu, Eleonora Parlanti, Michele Perna, Dávid Puskás, Pierluigi Rinaldi, Charlotte Simmonds, Yang Sun, Giacomo Venturi, Joris Witstok, Zihao Wu, Yongda Zhu

We present JWST/NIRSpec dense-shutter spectroscopy (DSS). This novel observing strategy with the NIRSpec micro-shutter assembly (MSA) deliberately permits a high number of controlled spectral overlaps to reach extreme multiplex while retaining the low background of slit spectroscopy. In a single configuration over the JADES Origins Field we opened shutters on all faint (F444W<30 mag) z$_{\rm phot}$>3 candidates in the MSA, prioritising emission-line science and rejecting only bright continuum sources. Using 33.6 and 35.8 ks on-source in G235M and G395M, we observed a single mask with ~850 sources, obtaining secure spectroscopic redshifts for ~540 galaxies over 2.5<z<8.9. The per-configuration target density in DSS mode is 4-5x higher than standard no- and low-overlap MSA strategies (<200 sources), with no loss in redshift precision or accuracy. Line-flux sensitivities are 30 percent lower at fixed exposure time, matching the expected increase in background noise, but the gain in survey speed is 5x in our setup, more than justifying the penalty. The measured line sensitivity exceeds NIRCam WFSS by a minimum factor of ~5 (i.e. ~25 in exposure time) at $\lambda$~4 $\mu$m, demonstrating that controlled overlap is a compelling method to gain deep, wide-band spectra for large samples. Crucially, we envisage the MSA could deliver even higher target allocation densities than what used here. We derive Balmer-line based SFRs, gas-phase metallicities (including a large sample suitable for strong-line calibrations), and identify rare sources (mini-quenched systems and broad-line AGN). This approach is immediately applicable wherever deep imaging enables robust pre-selection and astrometry, providing an efficient method to obtain large samples of faint emission-line galaxies, a compelling middle ground between the completeness of slitless surveys and the sensitivity and bandwidth of NIRSpec/MSA.

Ryan Hazlett, Jennifer Mead, Eli Visbal, Greg L. Bryan, Mordecai-Mark Mac Low, Mihir Kulkarni, Eric P. Andersson, Kaley Brauer, John H. Wise

We present an extension of our semi-analytic model that follows the formation of Population III stars and their metal-enriched descendants, incorporating dark matter halo merger trees from cosmological $N$-body simulations and feedback from reionization. Our extended model is calibrated using two complementary cosmological hydrodynamical simulations: Aeos, which resolves individual Population III and II stars to $z\sim14.6$, and Renaissance, which is lower resolution but follows large-scale metal-enriched star formation to $z \sim 11$. With a combined calibration, we capture small-scale physics of primordial star formation over a large range in halo mass. We find good agreement between our calibrated model and Aeos, reproducing the evolution in number of star-forming halos and total stellar mass. Achieving this agreement requires increasing the normalization of, flattening the redshift dependence of, and adding scatter to the commonly used critical mass threshold $M_{\mathrm{crit}}$. Our treatment of the delay between Pop III stellar death and subsequent Pop II star formation emphasizes the need to account for halos that have yet to transition to Pop II, since incomplete sampling of this delay in simulations limits physically motivated calibrations. Finally, we apply our model to larger-volume dark matter only simulations and predict $\sim10$ active Pop III sources at $z = 10$ lie within the area strongly lensed by galaxy cluster MACS J0416 with a magnification exceeding $\mu > 30$. These results demonstrate that semi-analytic approaches, when calibrated to hydrodynamical simulations, can provide accurate, computationally efficient predictions for the earliest stages of cosmic star formation.

Compact binary millisecond pulsars (also known as spiders) allow us to probe pulsar winds in their innermost regions, between the light cylinder (radius $\sim10^{7}$ cm) and the companion star (at $\sim10^{11}$ cm). Their flux is known to vary along the orbit, from radio to X-rays. During the past decade, gamma-ray orbital modulation (GOM) has been discovered in a handful of spiders, but its origin remains largely unknown. We present the results of a systematic search for GOM among 43 systems, selecting pulsed 0.1-1 GeV photons and using spin and orbital ephemeris from Fermi's Third Pulsar Catalog. We discover GOM from three spiders - PSR J1124-3653, PSR J1946-5403 and PSR J2215+5135 - and confirm four previous detections. In all seven cases so far, the GOM peaks near the pulsar's superior conjunction. The X-ray orbital light curves are usually in anti-phase, peaking when the pulsar is at inferior conjunction, but we find one case where both gamma-rays and X-rays peak around superior conjunction: PSR J1946-5403. We measure the modulated fractions of the GOM and find consistent values for all seven spiders, with an average $22.0\pm2.6\%$. Including eclipsing systems seen edge-on, we find no clear dependence of the modulated fraction on the orbital inclination (within $\simeq$45-90$^\circ$). Our results challenge previous models proposed to explain GOM in spiders, based on inverse Compton and synchrotron emission close to the companion, since these predict a clear dependence with orbital inclination (stronger modulation at high inclinations). We nearly double the number of GOM detections in spiders, showing that it is more common than previously thought.

T. Emil Rivera-Thorsen, Brian Welch, Taylor Hutchison, Matthew J. Hayes, Jane R. Rigby, Keunho Kim, Suhyeon Choe, Michael Florian, Matthew B. Bayliss, Gourav Khullar, Keren Sharon, Håkon Dahle, John Chisholm, Erik Solhaug, M. Riley Owens, Michael D. Gladders

At present, the best opportunity for detailed Lyman Continuum escape studies is in gravitationally lensed galaxies at z >= 2. Only one such galaxy currently exists in the literature with sufficient spatial magnification: The Sunburst Arc at redshift z = 2.37. Here, we present rest-frame optical JWST NIRSpec integral field observations of the Sunburst Arc that cover a large fraction of the source plane. From this dataset, we generate precise maps of ISM kinematics, dust geometry, ionization, and chemical enrichment. We find that the galaxy rotates but also shows strong, possibly dominant, signatures of turbulence, which are indicative of recent or ongoing major interaction. The cluster that leaks ionizing photons shows little variation in kinematics or dust coverage, but dramatically elevated ionization, indicating that photoionization is the predominant mechanism that creates paths for LyC escape. We conjecture that tidal stripping of H I gas due to an interaction could have removed a large portion of the neutral ISM around the LyC emitting cluster, making it easier for the cluster to completely ionize the rest.

Yolanda G. C. Frensch, François Bouchy, Gaspare Lo Curto, Alexandrine L'Heureux, Roseane de Lima Gomes, João Faria, Xavier Dumusque, Lison Malo, Marion Cointepas, Avidaan Srivastava, Xavier Bonfils, Elisa Delgado-Mena, Nicola Nari, Khaled Al Moulla, Romain Allart, Jose M. Almenara, Étienne Artigau, Khalid Barkaoui, Frédérique Baron, Susana C. C. Barros, Björn Benneke, Marta Bryan, Charles Cadieux, Bruno L. Canto Martins, Izan de Castro Leão, Amadeo Castro-González, Ryan Cloutier, Karen A. Collins, Nicolas B. Cowan, Eduardo Cristo, Jose R. De Medeiros, Xavier Delfosse, René Doyon, David Ehrenreich, Sergio B. Fajardo-Acosta, Thierry Forveille, Tianjun Gan, João Gomes da Silva, Jonay I. González Hernández, Nolan Grieves, Steve Howell, David Lafrenière, Christophe Lovis, Claudio Melo, Lina Messamah, Lucile Mignon, Christoph Mordasini, Louise D. Nielsen, Ares Osborn, Léna Parc, Francesco Pepe, Caroline Piaulet-Ghorayeb, Rafael Rebolo, Jason Rowe, Nuno C. Santos, Damien Ségransan, Keivan G. Stassun, Stephanie Striegel, Alejandro Suárez Mascareño, Stéphane Udry, Solène Ulmer-Moll, Diana Valencia, Valentina Vaulato, Gregg Wade, Cristilyn N. Watkins

Gas giant planets orbiting low-mass stars are uncommon outcomes of planet formation. Increasing the sample of well-characterised giants around early M dwarfs will enable population-level studies of their properties, offering valuable insights into their formation and evolutionary histories. We aim to characterise giant exoplanets transiting M dwarfs identified by TESS. High-resolution spectroscopic data are obtained in the optical and nIR, combining HARPS and NIRPS. We derive RVs via the cross-correlation function and implement a novel post-processing procedure to further mitigate telluric contamination in the nIR. The resulting RVs are jointly fit with TESS and ground-based photometry to derive the orbital and physical parameters of the systems. We confirm two gas giants transiting the low-mass stars TOI-3288 A (K9V) and TOI-4666 (M2.5V). TOI-3288 A hosts a Hot Jupiter with a mass of $2.11\pm0.08~M_{\rm Jup}$ and a radius of $1.00 \pm 0.03~R_{\rm Jup}$, with an orbital period of 1.43 days ($T_{\rm eq} = 1059 \pm 20~{\rm K}$). TOI-4666 hosts a $0.70_{-0.06}^{+0.05}~M_{\rm Jup}$ warm Jupiter ($T_{\rm eq} = 713 \pm 14~{\rm K}$) with a radius of $1.11 \pm 0.04~R_{\rm Jup}$, and an orbital period of 2.91 days. We identify a decrease in planetary mass with spectral type, where late M dwarfs host less massive giant planets than early M dwarfs. More massive gas giants that deviate from this trend are preferentially hosted by more metal-rich stars. Furthermore, we find an increased binarity fraction among low-mass stars hosting gas giants, which may play a role in enhancing giant planet formation around low-mass stars. The observed population trends agree with theoretical expectations, where higher metallicity can compensate for lower disk masses, and wide binary systems may influence planet formation and migration through Kozai-Lidov cycles or disk instabilities.

M. Grayling, S. Thorp, K. S. Mandel, M. Pascale, J. D. R, Pierel, E. E. Hayes, C. Larison, A. Agrawal, G. Narayan

We present BayeSN-TD, an enhanced implementation of the probabilistic type Ia supernova (SN Ia) BayeSN SED model, designed for fitting multiply-imaged, gravitationally lensed type Ia supernovae (glSNe Ia). BayeSN-TD fits for magnifications and time-delays across multiple images while marginalising over an achromatic, Gaussian process-based treatment of microlensing, to allow for time-dependent deviations from a typical SN Ia SED caused by gravitational lensing by stars in the lensing system. BayeSN-TD is able to robustly infer time delays and produce well-calibrated uncertainties, even when applied to simulations based on a different SED model and incorporating chromatic microlensing, strongly validating its suitability for time-delay cosmography. We then apply BayeSN-TD to publicly available photometry of the glSN Ia SN H0pe, inferring time delays between images BA and BC of $\Delta T_{BA}=121.9^{+9.5}_{-7.5}$ days and $\Delta T_{BC}=63.2^{+3.2}_{-3.3}$ days along with absolute magnifications $\beta$ for each image, $\beta_A = 2.38^{+0.72}_{-0.54}$, $\beta_B=5.27^{+1.25}_{-1.02}$ and $\beta_C=3.93^{+1.00}_{-0.75}$. Combining our constraints on time-delays and magnifications with existing lens models of this system, we infer $H_0=69.3^{+12.6}_{-7.8}$ km s$^{-1}$ Mpc$^{-1}$, consistent with previous analysis of this system; incorporating additional constraints based on spectroscopy yields $H_0=66.8^{+13.4}_{-5.4}$ km s$^{-1}$ Mpc$^{-1}$. While this is not yet precise enough to draw a meaningful conclusion with regard to the `Hubble tension', upcoming analysis of SN H0pe with more accurate photometry enabled by template images, and other glSNe, will provide stronger constraints on $H_0$; BayeSN-TD will be a valuable tool for these analyses.

We study the gravitational microlensing of various static and spherically symmetric non-singular black holes (and horizonless, non-singular compact objects of similar size). For pointlike sources we extend the parametrized post-Newtonian lensing framework to fourth order, whereas for extended sources we develop a ray tracing approach via a simple radiative transfer model. Modelling non-relativistic proper motion of the lens in front of a background star we record the apparent brightness as a function of time, resulting in a photometric lightcurve. Taking the star radius to smaller values, our numerical results approach the theoretical predictions for point-like sources. Compared to the Schwarzschild metric in an otherwise unmodified lensing geometry, we find that non-singular black hole models (and their horizonless, non-singular counterparts) at finite size tend to feature larger magnifications in microlensing lightcurves, contrary to the point-source prediction.

Building on initial work on the Thermodynamic Split Conjecture (TSC), which posits that black hole and cosmological horizon thermodynamics are generically inequivalent, we examine the consequences of that split for the Gibbons Hawking temperature and its role across cosmology. We consider many key results in both early and late universe cosmology and show that many important results such as those governing eternal inflation, vacuum tunneling, quantum breaking and primordial black holes can change. The analysis further reveals that small, TSC motivated corrections to horizon thermodynamics can subtly modify Friedmann dynamics, potentially helping to address the $H_0$ and $S_8$ tensions. The work thus provides a unified route from quantum gravity motivated thermodynamics to observational cosmology and motivates dedicated tests of the thermal laws governing the Universe itself.

In this work we consider the scattering between non-relativistic particles with different finite sizes. We first calculate their interaction potential and apply the partial wave method to obtain their scattering cross section. Our findings show that the particle size can significantly affect the scattering between non-relativistic particles. Then we apply such a study to direct detection of puffy dark matter. We find that the finite size of the target nucleus may introduce non-perturbative effects that differ from the scenario of point-like dark matter. For large-size dark matter particles, this non-perturbative regime in the dark matter nucleus scattering cross section effectively disappears; while for small values of the size-to-range ratio in the scattering process, a significant non-perturbative regime can maintain. Finally, for the direct detection of nugget-type puffy dark matter with a small number of constituent particles, we find that the stability conditions for the formation of bound-state dark matter can provide constraints on the dark matter nucleus scattering cross section.

We consider the implications of the modified dispersion relations, due to the noncommutativity of the spacetime, for a photon gas filling the early Universe in the framework of the Big Bang Nucleosynthesis (BBN) processes, during the period of light elements formation. We consider three types of deformations present in the dispersion relations for the radiation gas, from which we obtain the low temperature corrections to the energy density and pressure. The cosmological implications of the modified equations of state in the BBN era are explored in detail for all radiation models. The effects induced on the nucleosynthesis process by spacetime noncommutativity are investigated by evaluating the abundances of relic nuclei (Hydrogen, Deuterium, Helium-3, Helium-4, and Lithium-7). The primordial mass fraction estimates and their deviations due to changes in the freezing temperature impose an upper limit on the energy density of the deformed photon gas, which follows from the modified Friedmann equations. The deviations from the standard energy density of the radiative plasma are therefore constrained by the abundances of the Helium-4 nuclei. Upper limits on the free parameters of the spacetime noncommutativity are obtained via a numerical analysis performed using the \texttt{PRyMordial} software package. The primordial abundances of the light elements are obtained by evaluating the thermonuclear reaction rates for the considered noncommutative spacetime models. An MCMC (Markov Chain Monte Carlo) analysis allows to obtain restrictions on the free parameters of the modified dispersion relations. The numerical and statistical approach is implemented in the python code \texttt{PRyNCe}, available on GitHub.

We show that axions can be produced from Abelian-Higgs cosmic strings due to the axion-gauge coupling. The strong magnetic field is confined in the string, and the electric field is induced around the moving string, allowing axion productions from the dynamics of cosmic strings. Our numerical analysis on the string collision shows that a sizable number of axions can be produced at the reconnection, and further emissions occur from moving kinks afterward. Furthermore, the simulation on the string network shows multimodal axion emissions in the sense that axions are produced in both the low-energy and high-energy regimes. The former can contribute to the cold dark matter and the latter can be regarded as dark radiation. We found that the axion with sub-GeV mass can explain the current relic dark matter abundance and simultaneously predicts a sizable amount of dark radiation which can be probed by future observations.

Jiamin Liang, Mingqiu Li, Yu Gao, Wei Ji, Sichun Sun, Qi-Shu Yan

The observation of gravitational waves has opened a new window into the Universe through gravitational-wave astronomy. However, high-frequency gravitational waves remain undetected. In this work, we propose that spin systems can be employed to detect gravitational waves in this unexplored frequency regime. We derive the spin's response to gravitational waves and identify three distinct effects: the well-known Gertsenshtein effect, a metric-induced interaction, and the gravitational spin Hall effect. We focus on nuclear spins and utilize nuclear magnetic resonance to enhance the gravitational response, leveraging the advantages of long coherence time, high polarization, and a small gyromagnetic ratio. The proposed experimental scheme is capable of probing gravitational waves in the kilohertz to gigahertz range, with projected sensitivities reaching $\sqrt{S_h}\approx10^{-20}~\mathrm{Hz}^{-1/2}$.

The observability framework Kieker provides a range of analysis capabilities, but it is currently only able to instrument a smaller selection of languages and technologies, including Java, C, Fortran, and Python. The OpenTelemetry standard aims for providing reference implementations for most programming languages, including C# and JavaScript, that are currently not supported by Kieker. In this work, we describe how to transform OpenTelemetry tracing data into the Kieker framework. Thereby, it becomes possible to create for example call trees from OpenTelemetry instrumentations. We demonstrate the usability of our approach by visualizing trace data of the Astronomy Shop, which is an OpenTelemetry demo application.

The LIGO, Virgo, and KAGRA (LVK) gravitational-wave observatories have opened new scientific research in astrophysics, fundamental physics, and cosmology. The collaborations that build and operate these observatories release the interferometric strain data as well as a catalogue of observed signals with accompanying Bayesian posterior distributions. These posteriors, in the form of equally-weighted samples, form a dataset that is used by a multitude of further analyses seeking to constrain the population of merging black holes, identify lensed pairs of signals, and much more. However, many of these analyses rely, often implicitly, on the ability to reconstruct the likelihood and prior from the inputs to the analysis and apply resampling (a statistical technique to generate new samples varying the underlying analysis assumptions). In this work, we first provide a guide on how to reconstruct and modify the posterior density accurately from the inputs for analyses performed with the Bilby inference library. We then demonstrate and compare resampling techniques to produce new posterior sample sets and discuss Pareto-smoothing to improve the efficiency. Finally, we provide examples of how to use resampling to study observed gravitational-wave signals. We hope this guide provides a useful resource for those wishing to use open data products from the LVK for gravitational-wave astronomy.

The crucial role of hydrodynamic instabilities in soliton field theory is revealed. We demonstrate that the essential of soliton formation mechanism is the sound mode instability induced by thermodynamic instability. This instability triggers phase separation, where new thermal phases are generated to produce solitons. These solitons can be regarded as a coexistence state composed of a matter phase and a vacuum phase, with an interface proving surface tension to maintain dynamical equilibrium. The phase separation mechanism naturally allows the existence of vacuum bubbles, characterized by a vacuum phase surrounded by a matter phase with negative pressure. Furthermore, we show that the soliton interface resemble a fluid membrane, whose interface pressure satisfies a Young-Laplace-type relation, resulting in the emergence of the membrane instability induced by surface tension. In the thin-wall limit, the dispersion relation is analytically derived. This instability triggers topological transition of the interface, splitting a cylindrical interface into multiple spheres with a smaller total surface area. Such results highlight the duality between solitons and fluids.

Three black-hole low-mass X-ray binaries (LMXBs) in the Milky Way show rates of period decay which cannot be easily explained by standard mechanisms. Recently, it has been claimed that the anomalous period decays in two of these systems may be explained by dynamical friction due to very high dark matter (DM) densities around the black holes. We critically assess these claims by performing $N$-body simulations of binaries embedded in dense DM ``spikes". We simulate the previously-studied systems XTE J1118+480 and A0620--00, as well as studying the third binary Nova Muscae 1991 for the first time in this context. These simulations show that feedback on the DM distribution plays a crucial role and we rule out previously-claimed shallow DM spikes. We set lower limits on the steepness $\gamma$ of DM density profiles required to explain the period decay in these LMXBs, requiring $\gamma \gtrsim 2.15-2.20$ in XTE J1118+480 and A0620--00 and $\gamma \gtrsim 2.3$ in Nova Muscae 1991. Improved modeling of the long-term evolution of binaries embedded in DM spikes may allow us to exclude even larger densities in future.

We revisit the upper bound on the annihilation cross-section, $\langle\sigma v\rangle$ of a stable dark matter (DM) of mass $5\times10^2-10^{14}$ GeV by considering five different channels: $W^+W^-$, $b\bar{b}$, $\mu^+\mu^-$, $\tau^+\tau^-$, and $e^+e^-$. We use the observed electron and positron fluxes from CALET, DAMPE, HESS, positron flux from AMS-02, and gamma-ray flux from HAWC, GRAPES-3, CASA-MIA to constrain the annihilation cross-section. We also consider unstable DM of mass $10^3-10^{16}$~GeV decaying to $W^+W^-$, $b\bar{b}$, $\mu^+\mu^-$, $\tau^+\tau^-$, and $e^+e^-$ and derive the corresponding lower bound on the DM lifetime, $\tau_{\rm DM}$. We find that the latest data from CALET gives a stringent constraint on $\langle\sigma v\rangle$ in the low DM mass regime. For a typical DM mass of 1 TeV, we show that $\langle\sigma v\rangle_{{\rm DM~DM}\rightarrow\mu^+\mu^-}\gtrsim\mathcal{O}(10^{-24})~\rm cm^3/s$ is disfavored. On the other hand in the low mass regime, the AMS-02 gives a much stringent limit on the DM lifetime, excluding $\tau_{\rm DM\rightarrow\mu^+\mu^-}\lesssim\mathcal{O}(10^{27})$ s for a 1 TeV mass of DM. In the high mass regime, typically $M_{\rm DM}\gtrsim\mathcal{O}(10^5)$ GeV, HAWC and CASA-MIA give the strongest constraints on $\langle\sigma v\rangle$ and $\tau_{\rm DM}$.