Papers for Monday, Mar 16 2026

Accreted stars in the Milky Way (MW) preserve information about the progenitor galaxies where they formed in their chemical and kinematic properties. In this study, we use the chemo-dynamical signatures in the merger debris to approximate the posterior distribution of disrupted satellite properties at the time of infall. Adopting a simulation-based inference framework, we train an ensemble of normalizing flows using samples of merger debris from the Auriga suite of simulations of MW-like galaxies. Applying this methodology to a local sample of accreted stars in the MW, we infer the lookback times, stellar and halo masses, and halo mass merger ratios of several known accretion events in the Galaxy: Gaia Enceladus-Sausage (GES), Helmi streams, Heracles, I'itoi, LMS-1/Wukong, Sagittarius (Sgr), Sequoia and Thamnos. Our predictions align with the accretion time and mass estimates from the literature, and the expected relation between the progenitor stellar masses and debris metallicities across redshifts. The total stellar mass accreted from these events is predicted to be $2.2^{+1.1}_{-0.6}\times10^{9}~\rm{M_{\odot}}$, with GES and Sgr being the largest contributors. The predicted stellar mass accreted from fully disrupted progenitors in the stellar halo is $1.3^{+1.0}_{-0.5}\times10^{9}~\rm{M_{\odot}}$, which is consistent with previous mass measurements of this component. We provide a prediction for the evolution of the MW halo mass until the accretion of Sgr ($z\approx1$): specifically, we find that the mass growth of the Galaxy from the time of its first merger ($z\approx5$) to $z\approx2$ exceeds the total mass of the known progenitors accreted during that interval, suggesting the presence of unidentified substructures. Our estimate of the Galaxy halo mass after the Sgr merger, but prior to the accretion of the Magellanic Clouds, is $5.9^{+1.4}_{-1.1}\times10^{11}~\rm{M_{\odot}}$.

All observational evidence for dark matter and dark energy is so far exclusively gravitational. Hence, the dark sector may be equivalently described by a theory of the spacetime metric whose dynamics is affected by interactions with an unknown environment. Adapting open-system techniques, we have recently constructed such a general theory of open gravitational dynamics. Here we study a minimal and concrete realization of this theory that describes the late-time acceleration of the Universe. Our model provides a good fit to recent baryon acoustic oscillation measurements by construction, while avoiding violations of the null energy condition. Moreover, it leads to a set of correlated and observationally testable predictions. Studying the modified cosmological perturbation theory and compared to the $\Lambda$CDM model we find: a dissipative suppression of the gravitational-wave luminosity distance relative to the electromagnetic one; a modification in the evolution of the Bardeen potentials with a clear signal in the gravitational slip; and an enhancement of structure formation at low redshift. We present semi-analytical estimates of the magnitude of these effects and show that they lie within the reach of current constraints while providing clear targets for upcoming cosmological surveys.

Active Galactic Nucleus (AGN) feedback plays a critical role in galaxy formation and evolution. AGN-driven winds can significantly influence their host galaxies, although the details of their impact remain unclear. In this study, we investigate the feedback effects of AGN winds on idealized disc galaxies using the SWIFT hydrodynamical code with COLIBRE subgrid physics. We implement a new thermal AGN feedback model in which the energy injection coupling efficiency has a power-law dependence on the Eddington ratio of the black hole (BH) accretion rate, motivated by scaling relations for AGN winds from numerical models and observations. We simulate idealised Milky Way-mass galaxies, incorporating a BH, cold gas disc, stellar disc, and hot circumgalactic medium, within a static dark matter halo. We vary the BH mass and the slope and normalisation of the new coupling efficiency model. For a fixed BH mass, we find that while systematic trends with coupling efficiency exist, most galaxy and BH properties show only modest variations. This likely reflects BH self-regulation in the COLIBRE model, which modulates the effects of changes in the feedback efficiency, provided the BH mass is sufficiently high. Key exceptions are the BH accretion rate and mass growth history, and outflow behaviour, where lower coupling efficiencies lead to faster BH growth and weaker outflows, potentially helping to explain the presence of overmassive BHs at high redshifts. Varying the BH mass, however, has a much larger impact, confirming that BH mass remains the primary factor shaping galaxy and BH evolution in our simulations.

Ultraluminous infrared galaxies are powered by a combination of rapid star formation and active galactic nucleus (AGN) emission, but their relative importance is not always observationally clear. We study the galactic continuum spectrum of a cosmologically simulated $\sim 4 \times 10^{10} M_\odot$ stellar mass starburst galaxy at redshift $z\sim 4.4$ that refines down to resolve beyond the dust sublimation boundary of its super-Eddington-accreting $\sim 10^7 M_\odot$ supermassive black hole. We find that this system resembles the rare class of hot dust-obscured galaxy (Hot DOG), with a roughly flat (in $\nu F_\nu$) IR emission spectrum that sharply drops off at wavelengths $\lesssim 5~\mu\mathrm{m}$. Our system also matches with the observational properties of many Hot DOGs, including undergoing multiple galaxy mergers and being the most massive galaxy within a dense cosmological environment. The distinctive Hot DOG spectral shape in our system is caused by AGN-heated mid-IR warm dust, predominately starburst-heated far-IR cold dust, and a steep near- to mid-IR cutoff caused by strong absorption in the dense ISM of the galactic nucleus, rather than the dust torus itself. This system is lower luminosity ($L_\mathrm{IR} \sim 2 \times 10^{12} L_\odot$) than those detected by the WISE survey at similar redshifts, but will be a prime target for future far-IR surveys such as PRIMA. Our results show that Hot DOGs can naturally result as a transitional phase during rapid AGN accretion, but before significant AGN-driven outflows clear optically thin paths.

Gravitational waves (GWs) exhibit wave-optics effects when their wavelength is comparable to the scale of the gravitational lens. This may occur in lensing from galactic subhalos in GWs emitted by binary black-hole mergers, and is gaining interest as a novel probe of dark matter. Predictions for observables in these cases ultimately rely on evaluating a Fresnel integral that quantifies the effect of lensing on the amplitude of a GW at a given frequency. However, numerical evaluation of this Fresnel integral is tricky, and several algorithms and publicly available codes that implement it have been developed. Here, we show that the dependence of this integral on the lens position can be written as a two-dimensional Fourier transform. Modern FFT techniques then enable rapid evaluation at all-sky positions simultaneously for general lenses without symmetry. Vectorization of FFT routines allows for derivatives with respect to model parameters to be obtained with only incremental additional computational cost. If the lens is axisymmetric, further speedups can be achieved with recently developed techniques for non-uniform fast Hankel transforms. To demonstrate, we make available Fresnel Integral Optimization with Non-uniform trAnsforms (FIONA), an efficient and accurate code that is significantly faster than current methods for dense source grids, reaching 2-3 orders of magnitude speedups for $\sim 10^6$ GW-emitting points. As part of FIONA, we developed code that provides vectorized non-uniform fast Hankel transforms that may have other uses (e.g., calculation of cosmological two-point correlation functions) beyond those considered here.

We investigate the cosmological information in higher-order statistics of the cosmic microwave background (CMB) lensing convergence field for a near-term experiment with noise properties similar to the Simons Observatory (SO). Using a fully field-level forward-modeling pipeline based on ray-traced simulations from the MassiveNuS suite and realistic SO-like CMB lensing reconstruction, we naturally include nonlinear structure formation, post-Born effects, and higher-order reconstruction noise. We measure several non-Gaussian statistics, including Minkowski functionals, peak and minima counts, moments, and wavelet-scattering coefficients. We train Gaussian-process emulators to model each statistic's dependence on the matter density fraction $\Omega_m$, the scalar power spectrum amplitude $A_s$, and the neutrino mass sum $M_\nu$. We quantify the relative information gain these statistics provide beyond the lensing power spectrum and identify which are most robust to reconstruction noise. We find that morphology-based statistics, particularly Minkowski functionals and peak/minima counts, offer significant complementary constraining power: combining all non-Gaussian statistics with the power spectrum yields reductions of 40% and 38% in the marginalized uncertainties on $\Omega_m$ and $A_s$, respectively, and a 70% reduction in the one-sided uncertainty on $M_\nu$. These gains remain non-negligible even when the power spectrum is extended to larger scales and combined with primary CMB and BAO data, with Minkowski functionals providing an additional 11% improvement in $\sigma(M_\nu)$ and 35% in $\sigma(\Omega_m)$ beyond the extended power spectrum. By contrast, moments and wavelet-scattering coefficients provide more limited gains at SO noise levels. Our results highlight the potential of non-Gaussian statistics to enhance cosmological constraints from SO and future CMB surveys.

The Large High Altitude Air Shower Observatory (LHAASO) has revealed a population of Galactic gamma-ray sources radiating beyond 100 TeV, but the nature of several of them is still uncertain. In this contribution, we explore the idea that some of these ultrahigh-energy emitters are not powered by currently active accelerators, but by the fossil remains of microquasars (MQs). We consider systems in which mass transfer onto the stellar-mass black hole has already stopped, so that the central engine and its jets are permanently quenched. During the active phase, powerful transrelativistic jets inflate a hot cocoon whose interior is filled with cosmic rays (CRs) accelerated at the jet termination shocks. Once the jets switch off, the cocoon enters a long afterlife stage in which it behaves as a large reservoir of PeV CRs. If the remnant lies in or near a star-forming region, these relic CRs can still interact with dense clumps and molecular clouds, inside the cocoon or in the surrounding interstellar medium, leading to delayed gamma-ray emission via inelastic pp collisions and the subsequent decay of neutral pions. We present a time-dependent model for the jet-cocoon system, follow the evolution of the CR population during and after the MQ phase, and discuss the conditions under which the resulting microquasar remnants can account for some of the unidentified LHAASO sources.

We study the early stages of an oscillatory reheating phase in an inflaton plus spectator axion-SU(2) system, including both an axion-gauge Chern-Simons coupling $\chi F\tilde{F}$ and a gravitational Chern-Simons coupling $\chi R\tilde{R}$. Assuming an isotropic SU(2) background configuration of chromo-natural type and quadratic potentials, we numerically solve the coupled background and tensor perturbation equations during the first e-fold of reheating. The gravitational Chern-Simons term induces a helicity-dependent modification of the tensor kinetic coefficient, yielding a chiral enhancement of the tensor power spectrum on the order of tens of percent for a representative benchmark. We illustrate how such an early-time enhancement can map to a narrow feature in the present-day stochastic gravitational wave spectrum, potentially relevant for upcoming and proposed space-based detectors, while a fully self-consistent determination of the peak scale requires scanning comoving wavenumbers and specifying the reheating history.

The Neptunian desert is a distinct lack of Neptune-sized planets at short orbital periods, purportedly carved by photoevaporation and tidal circularization following high-eccentricity migration. Constraining these processes and how they vary across different host-star spectral types requires the detailed characterization of planets in the desert and around its boundaries. In this study, we confirm the planetary nature of a massive super-Neptune identified by TESS around the M0 dwarf TOI-672. We analyse photometry from TESS and ExTrA and precise radial velocity measurements taken with the recently commissioned Near-InfraRed Planet Searcher (NIRPS) and HARPS spectrographs. We measure the planetary orbital period, radius, and mass of 3.634 days, 5.31 +0.24 -0.26 Rearth, and 50.9 +4.5 -4.4 Mearth, respectively. Our findings place TOI-672 b within the Neptunian ridge, a pile-up of planets from 3--5 days at the Neptunian desert boundary. We then use a novel approach to determine the desert boundaries in period-radius space and instellation-radius space, and, for the first time, compare the Neptunian desert boundaries for planets orbiting FGK versus M dwarf stars. We determine that the boundary ridge shifts slightly inward from 3.3 +- 1.4 days for FGK host stars to 2.2 +- 1.0 days for M dwarf host stars; these values do not statistically significantly differ from each other, and the shift to shorter periods for M dwarf planets is smaller than theoretical photoevaporation models predict. We also find that TOI-672 b is a single-planet system within the sensitivity limits of our RV and TTV datasets.

The detection of Cosmic Microwave Background primordial $B$-mode polarization would constitute a ``smoking gun" signal of primordial gravitational waves. However, this measurement requires accurate removal of polarized Galactic foregrounds to avoid systematic biases when estimating the tensor-to-scalar ratio. Methods based on Machine Learning techniques (ML), such as Convolutional Neural Networks (CNNs), have recently been proposed as alternative foreground cleaning techniques, but their applicability to real data relies on their ability to generalize beyond the models assumed during training. In this work, we focus on a variety of foreground models (FMs) used for training and conduct a systematic study of the generalization properties of a CNN-based method. We train various CNN architectures on simulations generated from different Galactic FMs, and test their performance on models not used during the training. By characterizing the statistical properties of the FMs using variance, skewness, and Shannon entropy, we define a statistical complexity hierarchy among them. We show that training on the more complex FMs reduces bias and improves precision when testing on unseen FMs, whereas training on the simplest model could introduce systematic errors. These results evidence that a lack of generalization is a relevant source of systematic uncertainty, and emphasize the importance of understanding the impact of the models assumed during training in ML-based methods before applying them to real data.

We present the first results from the SMUGGLE-Ring project, a suite of simulations employing the SMUGGLE ISM and stellar feedback model to explore nuclear structures in Milky Way-mass galaxies. We discuss results from three simulations evolved for 5 Gyr in isolation, in which we vary the classical bulge mass, while keeping the disk and halo structures identical. Nuclear stellar disks and rings emerge exclusively in our bulge models, with more massive bulges associated with earlier formation and more extended initial gas reservoirs shortly after bar formation. After gas depletion via active star formation, the nuclear stellar disks bifurcate into pressure-supported nuclear star clusters (NSCs, $v_{\phi}/\sigma_R < 0.7$) and rotationally supported nuclear stellar rings (NSRs, $v_{\phi}/\sigma_R = 1.2$--1.7, radii 0.64--0.76 kpc). The bulgeless model fails to build up and sustain stable nuclear gas disks against feedback disruptions. The enclosed stellar mass of NSCs ($\sim10^{9}\Msun$) dominates over that of NSRs ($\sim10^{8}\Msun$). The star formation rates decline over time due to gas depletion (NSCs 0.1--1 $\Msun$/yr, NSRs 0.01--$0.1 \Msun$/yr). Kinematics reveal outward-shifting rotation peaks with $\sigma$-drops in NSRs, while a fraction of stars in NSCs exhibits radial shift after 3 Gyr. These findings support inside-out NSD formation via secular bar evolution, with NSRs tracing the star-forming outer edge of the nuclear gas disk and NSCs forming the kinematically hotter inner component. The range of nuclear stellar disk sizes (0.25--0.76 kpc) falls within the observationally inferred ranges, but the existence of larger rings would require external gas flow and/or a longer period of evolution. Future SMUGGLE-Ring extensions will incorporate varying gas fractions, tidal/merger effects, and the circumgalactic medium to further elucidate nuclear diversity and outliers.

Precise halo mass ($M_h$) measurements are crucial for cosmology and galaxy formation. HaloFlow introduced a simulation-based inference (SBI) framework that uses state-of-the-art simulated galaxy images to precisely infer $M_h$. However, for HaloFlow to be applied to observations, it must be generalizable even when the underlying galaxy formation physics differ from those in the simulations on which it was trained. Without this generalization, HaloFlow produces biased and overconfident $M_h$ posteriors when applied to simulations with different physics. We introduce HaloFlow$^{\rm DA}$, an extension of HaloFlow that integrates domain adaptation (DA) with SBI to mitigate these cross-simulation shifts. Using synthetic galaxy images forward-modeled from the IllustrisTNG, EAGLE, and SIMBA simulations, we test two DA methods: Domain-Adversarial Neural Networks (DANN) and Maximum Mean Discrepancy (MMD). Incorporating DA significantly reduces bias and improves calibration, with MMD achieving the most stable performance, lowering the normalized residual metric, $\beta$, by an average of 31% and up to 57% when trained and tested on different simulations. Overall, HaloFlow$^{\rm DA}$ produces more robust, less biased with similar precision, $M_h$ constraints than the standard approach using the stellar-to-halo mass relation. HaloFlow$^{\rm DA}$ enables consistent, simulation-trained inference models to generalize across domains, establishing a foundation for robust $M_h$ inference from real HSC-SSP observations.

Paper I in this series introduced a model in which seed radio bursts produced by a hotspot anchored in the magnetosphere of a highly-magnetic neutron star (NS) are greatly amplified by strong gravitational self-lensing and thus give rise to Fast Radio Bursts (FRBs). Key features of the FRB population such as the observed dichotomy between repeating and non-repeating sources, their large luminosities and the high-energy power-law distribution of their bursts naturally arise in the model from the amplification dependence on the relative orientation of the rotation axis with respect to the hotspot and the line of sight. Here we compare the model predictions with Five-hundred-meter Aperture Spherical radio Telescope (FAST) data from repeaters and with the general population of FRBs. We find that the burst energy distribution from FRB 20121102A can be explained by assuming two antipodal hotspots in the NS magnetosphere, both producing seed bursts with the same log-normal energy distribution. This scenario implies a well-aligned system geometry, with the rotation axis, line of sight, and hotspot sites separated by $\lesssim 2$°. Similar constraints are found for FRB 20201124A and FRB 20220912A, and weaker ones for FRB 20190520B, owing to its smaller burst sample. We also show that precession of the NS rotation axis can explain the time evolution of the burst energy distribution from FRB 20121102A as well as its temporary disappearance. In application to a cosmological population of randomly-oriented sources the model predicts distance and fluence distributions of FRBs in good agreement with those from a completeness-selected subsample of the first CHIME/FRB catalogue, provided the energy distribution of seed bursts spans a range of ${\sim10^{35}-10^{38}}$ erg.

We present some of the first infrared spectral maps acquired by SPHEREx. These maps, which to our knowledge are the largest of their type ever compiled in the near-infrared, reveal multiple strong lines due to interstellar ices and polycyclic aromatic hydrocarbons (PAHs) throughout the Cygnus X and North American Nebula regions. The maps emphasize the strongest features arising from the 3 $\mu$m H$_2$O, 4.27 $\mu$m CO$_2$, and 4.67 $\mu$m CO lines and the 3.28 $\mu$m PAH feature, all of which are detected over large areas with complex and filamentary spatial distributions. The ice absorption maps of H$_2$O and CO$_2$ in particular broadly trace dense, cold, and well-shielded regions across Cygnus X, consistent with the established picture of efficient ice formation in dense molecular clouds. The interstellar ice features are also detected abundantly in diffuse absorption over wide areas. The relative strength of the H$_2$O and CO$_2$ features varies among different lines of sight, indicating possible differences in local physical conditions or chemical variations. The 3.28 $\mu$m PAH emission correlates with the emission from the 7.7 and 11.2 $\mu$m features, but shows small differences that may trace the grain size distribution and variations in the ambient UV field. SPHEREx all-sky spectral imaging, of which only a small fraction is showcased in this work, will support numerous science investigations including the structure of the Galaxy, the physics of the interstellar medium, and the chemistry of stars.

With more than 6000 exoplanets discovered so far, about 12 percent are hot Jupiters. Their large sizes and short orbital periods make them valuable targets for studying planetary formation, atmospheres, and orbital evolution. We present a homogeneous analysis of the WASP-19 b system using a 15 year dataset to investigate both its orbital dynamics and atmospheric properties. We test whether the transit times show evidence for tidal orbital decay, apsidal precession, or periodic perturbations from an additional body, and we also construct a photometric transmission spectrum. Multi-wavelength light curves are modeled with PRISM to account for starspots, and linear, quadratic, and cubic ephemeris models are fitted to the transit timing residuals. Our dataset includes 27 new transits and reveals no statistically significant periodic signal. Although none of the tested models fully reproduces the timing scatter, the transit times show systematic deviations from a constant period and are best described by the cubic ephemeris, indicating a slow long-term trend over the full baseline. This behavior is more consistent with gradual apsidal precession than with monotonic tidal decay. A precession model yields a rate of 1.00 +/- 0.12 x 10^-4 rad per orbit and a planetary Love number k2p = 0.107 +/- 0.08. The transmission spectrum shows signatures of Na, K, and H2O, with no strong evidence for TiO or VO. These results suggest that apsidal precession may dominate the long-term orbital evolution of WASP-19 b. Continued high-precision timing and spectroscopic observations are needed to further test this scenario.

The magnetic fields in the solar polar region are important to our understanding of the internal dynamo process, the global coronal structure, and the origin of the solar wind. The inference of polar fields based on spectropolarimetric observation is highly model-dependent and can suffer from various systematic effects. Here we analyze a raster map of the southern polar region taken by the Hinode Spectro-Polarimeter, utilizing the Stokes Inversion based on Response functions code. The inversions provide height-dependent vector magnetic field maps between optical depths $\log_{10}\tau = -2$ and $0$. We examine the impact on the total magnetic flux estimate from adopting (1) 1- vs 2-component atmospheric models via a "filling factor" parameter and (2) different analysis schemes. At $\log_{10}\tau = -1.5$, the polar magnetic flux is estimated to be $(1.84 \pm 0.03) \times 10^{21}$ Mx and $(1.38 \pm 0.02) \times 10^{21}$ Mx under the 1- and 2-component atmosphere assumption, respectively. The magnetic flux is approximately constant or increases slightly with height, respectively. We find that the 2-component (1-component) configuration is preferred for 58.3% (32.3%) of the pixels. Different initial guesses, including the input atmosphere model and the filling factor, as well as different inversion settings, can significantly affect the results, especially for locations with weaker polarization signals. Our work highlights the importance of including unresolved magnetic structures or stray light into consideration. Model degeneracy and the convergence to local minima limit the precision of the polar magnetic flux inference (no better than several tens of percent in this case). Higher-resolution observations and advanced inversion and disambiguation algorithms may alleviate these limitations.

Collisional ring galaxies are produced by the collision of a disk galaxy with a compact galaxy plunging through the disk, forming a ring-shaped expanding density wave, triggering star formation at its wake. The wave expansion is expected to produce negative stellar age gradients in radial profiles of post-collision stellar populations. Integral field spectroscopy combined with stellar population synthesis allows us to spatially resolve the stellar populations, to separate the post-collision and pre-collision components, and to produce the radial profiles. We analyse three candidate galaxies: Arp~143, NGC~2793, and VII~Zw~466. Observations were performed with the Calar Alto 3.5~m~telescope using the PMAS/PPak spectrophotometer. NGC 2793 presents a positive stellar age gradient, dismissing the collision hypothesis. For Arp~143 and VII~Zw~466, we found negative stellar age gradients for the youngest stellar populations, up to the ring radii, consistent with the collision hypothesis. We estimated that the collisions occurred $\sim$300~Myr and $\sim$100~Myr (expansion velocities of 33~$\pm$~10 km s$^{-1}$ and 108~$\pm$~26 km s$^{-1}$), respectively, before the density waves reached the observed ring radii. A spatially resolved analysis of the specific star formation histories (sSFH), reveals an expected star formation enhancement following the collision. The sSFH also allowed to identify the most probable intruder galaxy for VII~Zw~466. We report new redshifts for its group members. Finally, radial profiles of light contributions from pre-collisional and post-collisional stars show that the density wave dragged old pre-collisional stars along, as predicted by simulations.

Interacting dark sector (IDS) models provide a commonly explored extension of the standard $\Lambda$CDM cosmology, allowing for non-gravitational energy--momentum exchange between cold dark matter (CDM) and dark energy (DE). Although such models can be constructed to reproduce the same background expansion history as $\Lambda$CDM, their impact on the growth of cosmic structures is fundamentally different and requires a careful treatment of cosmological perturbations. In this work, we develop the one-loop Standard Perturbation Theory (SPT) formalism for IDS cosmologies without invoking the Einstein--de~Sitter (EdS) approximation. We show that even weak dark sector interactions induce a non-trivial time dependence in the perturbative kernels, leading to a breakdown of the EdS approximation commonly assumed in $\Lambda$CDM analyses. By deriving and numerically solving the evolution equations for the second- and third-order kernels, we compute the corresponding one-loop corrections to the matter power spectrum and find that the resulting deviations can significantly exceed the percent level, even for small interaction strengths. Our results demonstrate that nonlinear corrections are systematically enhanced in IDS models and that neglecting the full time dependence of the kernels can lead to biased predictions on mildly nonlinear scales. These findings establish the necessity of a time-dependent perturbative treatment for IDS scenarios and provide a robust framework for precision tests using nonlinear large-scale structure (LSS) observables.

A long standing problem in high energy astrophysics is the nature of galactic accelerators of particles with energies above PeV. Such objects are sources of galactic cosmic rays and can produce PeV-regime photons observed by ground-based observatories. Among very likely accelerators are astrospheres of pulsars in gamma-ray binaries. These binaries have long been observed as bright sources of TeV gamma-rays. Recently, 2D relativistic magnetohydrodynamic (rMHD) simulations have shown that the astrospheres can accelerate particles to energies well above PeV, provided that they harbor a Gauss-range magnetic field. Such a strong field is necessary in the region of two colliding winds: the relativistic outflow of the pulsar or accreting black hole and the wind of its stellar companion, a massive early-type star. Here, the wind collision region is explored as the site of PeV protons acceleration. The local structure of colliding flows is illustrated using rMHD simulations of a powerful pulsar wind in 2D and 3D models. The relativistic outflow of a pulsar or black hole, evolving inside the strongly magnetized stellar wind, have an elongated shape and surrounded by a kind of magnetic cocoon providing favorable conditions for acceleration of ultra high energy ions. The simulated spectra of particles, accelerated by intermittent relativistic turbulence in these systems, have piece-wise power-law shape and extend well above PeV energies for powerful outflows. The model indicated that gamma-ray binaries harboring a powerful relativistic outflow, produced either by a pulsar or accreting black hole, can be bright sources of synchrotron MeV-regime photons and multi-PeV regime gamma-rays, as recently detected from galactic microquasars like Cyg X-3. The Gauss-range magnetic field of a massive star wind strongly influences the non-thermal emission of gamma-ray binaries with relativistic companions.

We present results from examining the silicate cloud modeling of four JWST-observed hot Jupiters in the context of brown dwarf theory to further explore signatures of formation in present-day atmospheres. We contextualize our understanding of protoplanetary disk refractory chemistry with empirical evidence from chondritic meteorites to show that giant planets forming and accreting in the outer disk adopt their stellar Mg/Si value. We show that current silicate cloud species determinations of WASP-17 b, WASP-107 b, WASP-39 b and HD 189733 b are in line with predictions laid out in Calamari et al. 2024 based on each system's host star Mg/Si ratio, further supporting this hypothesis. We discuss physical motivations for potential atmospheric scenarios where apparent silicate cloud species is not in agreement with that predicted by its host star chemistry. Additionally, we compare current transit spectroscopy for three of these four exoplanets against brown dwarf spectra to examine molecular absorption trends across the substellar mass temperature regime.

In this article, we investigate the effect of systematics on weak lensing beyond the standard {\Lambda}CDM paradigm. Specifically, we consider the 2- and 3-point statistics of the shear field for the set of cosmological models, including CPL dark energy, interacting dark energy (IDE), and Hu-Sawicki f (R) modified gravity. We consider two major systematics such as photometric redshift uncertainty and intrinsic alignment A IA . Our findings are derived from the Fisher matrix. These results indicate that {\sigma} z and A IA can substantially degrade constraining power, especially for f (R) gravity. Moreover, it also highlights the critical role of higher-order statistics and the need for robust systematic control for future cosmological surveys.

The fate of massive stars during the latest stages of their evolution is highly dependent on their mass-loss rate and geometry. These processes have a significant influence on stars with masses between 25 and 40 Msun, i.e., type II SN progenitors. We aim to investigate the mass-loss history, geometry, and physical conditions of the yellow hypergiant in a post-RSG stage, IRAS 17163-3907. We place it in context with another famous yellow hypergiant, IRC+10420. We combine M-band high-resolution CRIRES+ spectroscopy with VLTI/MATISSE mid-infrared L-band interferometry, and FORS2 optical spectropolarimetry to probe both the small-scale circumstellar structure and the large-scale dusty environment of IRAS 17163. The CRIRES+ spectrum provides the first M-band coverage of IRAS 17163, revealing prominent low-excitation metal lines and hydrogen recombination features, but lacking the pronounced CO absorption seen in IRC+10420. The MATISSE observations reveal the first high angular scales of the source in the L-band and spatially resolve the Br{\alpha} line-emitting region, which hints at a marginally asymmetric and variable ionised wind. FORS2 spectropolarimetry points to deviations from perfect spherical symmetry also on larger scales. The data show no evidence for a binary companion within the explored parameter space, indicating that the observed clumpy and time-variable mass loss is likely intrinsic to the star rather than companion-driven. Our results demonstrate that IRAS 17163 hosts a dense, structured, and time-variable wind, coexisting with extended dusty shells. The comparison with IRC+10420 highlights diversity among post-RSG YHGs. These findings emphasise the role of clumpy and near-symmetric mass-loss in shaping the circumstellar medium of evolved massive stars, with implications for their subsequent evolution and core-collapse supernova progenitor properties.

A handful of planetary systems hosting a Hot Jupiter have been subsequently found to also host long-period giant planets. These ``cold Jupiters,'' giant planets residing beyond the snow line ($\sim$3\,au), play an important role in the dynamical evolution of the system as a whole. In this work, we investigate the detectability of cold Jupiters around a sample of 28 well-studied Hot Jupiter host stars to estimate the occurrence rate of this distinctive system architecture. We perform extensive simulations using the combination of all publicly available radial velocity (RV) data for those stars with synthetic RV data. The synthetic data test observing strategies along three axes: cadence, duration, and measurement precision. For each scenario, we determine detection limits based on the semi-major axis at which a 1 Jupiter mass planet would be recovered 50\% of the time. We find the following: 1) the existing RV data are remarkably insensitive to these Hot Jupiter/Cold Jupiter pairs; 2) the total baseline over which an observational campaign is carried out is the dominant factor in our ability to detect cold Jupiters; and 3) the results are relatively insensitive to the individual RV measurement precision. We conclude that metre-class telescopes with lower RV precision are ideally suited to surveying Hot Jupiter-cold Jupiter systems.

Gravitational waves constitute a powerful probe of the underlying theory of gravity. In extensions of general relativity, additional degrees of freedom, such as scalar fields in the gravitational sector, can modify their propagation through changes in the effective friction term and propagation speed. These modifications may potentially induce resonant phenomena leading to distinctive signatures in the gravitational wave spectrum. One important aspect to be investigated is whether the resonances can be strong enough to enhance the underlying background of primordial tensor modes to levels detectable by upcoming gravitational wave detectors, such as LISA or the Einstein telescope. The characteristic peaks in the SBGW spectrum depend on the parameters of the resonant model as well as on the parameters of the primordial tensor spectrum, such as $r$ and $n_{t}$. Thus these resonance effects open a powerful pathway to explore physics of the very early Universe by amplifying otherwise feeble signals to experimentally detectable levels. Here we analyze how the signals of the primordial Universe can resonate in these scenarios, bringing the early universe physics into the realm of experimental access.

Bars in disk galaxies slow down as they transfer their angular momentum to their dark matter halo via dynamical friction from near-resonant orbits. This bar-halo dynamical friction can become ineffective once phase mixing erases the phase-space gradient around the main resonances. We present fully self-consistent $N$-body simulations of a Milky Way-like disk galaxy with a single dwarf-galaxy fly-by in prograde and retrograde orbits before, during, and after bar formation. In our models, the fly-bys do not trigger a long-lived tidal bar; the bar forms on essentially the same time as in the isolated model. After the encounter, however, all perturbed models develop bars that are stronger and slower than in the isolated one. The final pattern speed depends little on the encounter time, but it does depend on the encounter direction relative to the disk rotation: prograde encounters slow the bar more than retrograde ones. The angular-momentum evolution shows that the disk loses its angular momentum and the halo gains it, consistent with bar-halo friction. By probing the particle distribution of the halo in angle-action space, we demonstrate that the isolated bar enters a metastable, saturated state with a flattened distribution in the phase space around the bar's corotation resonance, whereas a dwarf passage excites long-lived fluctuations in the halo that restore the phase-space gradients near the corotation and thereby sustain the bar-halo friction. This mechanism explains the continued slowdown and growth of bars after fly-bys. It may be relevant to the Milky Way, whose bar formed near the epoch of a major ancient accretion event, suggesting that an early encounter could have influenced the subsequent secular evolution of the bar.

The extent of envelope stripping in the progenitor stars is directly reflected in the diversity of spectral features observed in stripped-envelope supernovae (SESNe). Through extensive spectral observation and analysis, we aim to clarify the statistical differences between the subclasses of SESNe. The Tsinghua Supernova group obtained 249 optical spectra of 62 SESNe during the years from 2010 to 2020, covering phases from $-$16 to over 190 days relative to maximum light. Most spectra were obtained during the photospheric phases after the supernova explosion. For each spectrum, the pseudo-equivalent widths (pEWs) and blueshift velocities of principal lines were measured. We further investigated the common spectral features by analysing their velocity and strength correlations across all subtypes. We identify the feature near 6200~Å in SNe Ib as H$\mathrm{\alpha}$ through comparison with SNe IIb and Ic, which resolves inconsistent literature interpretations. Our finding reveals prevalent residual hydrogen in SNe Ib, further supporting a continuous stripping sequence from SNe IIb to Ib. We observe a trend in increasing velocity among different subtypes of stripped-envelope SNe, with SNe IIb exhibiting the lowest line velocities, followed by Ib, Ic, and Ic-BL. Typically, the O~I lines in SNe Ic/Ic-BL are stronger than those seen in SNe IIb/Ib. In nebular phases, the [Ca II] emission dominates over [O I] in SNe IIb/Ib while [O I] is stronger in SNe Ic, including the He-rich SN 2016coi. This spectral dichotomy implies that progenitors of SNe Ic (BL) have more massive CO cores and hence higher initial masses.

We develop a fully relativistic framework to study the rotational response of gravitationally coupled two-fluid neutron stars within the slow-rotation approximation. Treating the two components as independently conserved perfect fluids interacting only through spacetime curvature, we derive the coupled equilibrium and frame-dragging equations and exploit their linear structure to construct a basis decomposition of the rotational response. This formulation leads to a natural definition of the effective total moment of inertia, which generalizes the single-fluid concept and depends solely on the equilibrium background. It further reveals that the coupled system admits two intrinsic collective rotational eigenmodes, characterized by distinct eigen-moments of inertia, even in the absence of relative rotation between the fluids. Applying this framework to neutron stars containing dark matter, we explore how the presence of an additional gravitationally bound component modifies the global rotational response and its relation to tidal deformability. Our results demonstrate that the persistence or breakdown of rotational-tidal universality in two-fluid neutron stars is governed by dark-sector microphysics rather than by the mere presence of an additional component, and establish a unified framework for interpreting rotational observables, intrinsic mode structure, and universal relations in multi-component relativistic stars.

The cessation of star formation in galaxies, known as 'quenching', is a complex, multi-scale process which has been theorized to be linked to galaxy mergers. In this paper, we investigate the potential role of mergers in quenching galaxies in the IllustrisTNG cosmological hydrodynamical simulation. We track the evolution of over 11,000 central galaxies in the simulation with stellar mass $M_\star \ge 10^9 M_\odot$ at $z = 0$ throughout the entirety of cosmic history. We compare their star formation and merger histories to test whether mergers are necessary or sufficient for inducing quenching in the simulation. Only a very small fraction of mergers (about 3 per cent of major mergers and about 12 per cent of all mergers) lead to quenching within 1 Gyr, indicating that mergers are not sufficient by themselves to cause quenching. Furthermore, the vast majority of quenching events are not preceded by a merger within 1 Gyr. Once random coincidences are accounted for and a stellar mass-matched control sample is applied, no merger excess is observed. Hence, mergers are clearly not necessary for quenching to occur in the simulation. Finally, we perform a series of random forest classification and regression analyses to assess the integrated role of mergers in galaxy quenching and supermassive black hole growth in IllustrisTNG. We determine that secular processes dominate the growth of supermassive black holes and the quenching of central galaxies in this simulation, in stark contrast to prior theoretical expectations from idealized hydrodynamical simulations.

Small scale clumps of ionized gas have been suggested by observations in interstellar medium and circumgalactic medium. The propagation of radio signals can be deflected by these plasma clumps, i.e. plasma lensing. One observable consequence is the magnification and demagnification of background sources. These effects distort the observed luminosity function and potentially introduce bias into population studies. In this work, we investigate these effects on fast radio bursts using Gaussian plasma clumps distributed across multiple lens planes within a small field of view. The central electron density for each clump is sampled from uniform, log-normal, and Gaussian distributions. Two analytical models are employed to mimic the intrinsic luminosity function. Our results show that plasma lensing can modify the observed luminosity functions. On one hand, our model shows that radio sources may be demagnified below the detection threshold, the strength varies between ~1-15% depending on the ionized gas model and the source redshift. On the other hand, magnification can produce anomalously bright sources at the high luminosity end. Both effects introduce potential biases in inferred source properties. The lensing strength correlates with the power spectrum of free electron density. However, scattering effect in the host galaxy or in the Milky Way can suppress the plasma lensing effects.

We derive bolometric correction functions for oscillation mode amplitudes observed by the different cameras of the ESA PLATO mission. Such corrections between bolometric (full light) and mission instrument-specific amplitudes enable comparisons to theoretical expectations and amplitude conversion between different photometric missions, which is essential for proper detectability yields and target selection. Bolometric correction functions were calculated assuming a Planck function approximation for the stellar spectral flux distribution. The calculations follow the procedures applied in earlier analyses for the NASA Kepler and TESS missions. We derived power-law and polynomial parametrisations of the bolometric corrections with $T_{\rm eff}$. We find that on average, oscillation mode amplitudes from PLATO's normal cameras (N-CAMs) are expected to be ~6.7% lower compared to Kepler, and ~12.5% higher compared to TESS. A significant average amplitude ratio of ~25% is expected for amplitudes measured using the blue PLATO fast camera (F-CAM) compared to TESS. We find that observations of bright solar-like oscillators, especially with PLATO's F-CAMs, would provide an important test of the predicted corrections.

We report the spectroscopic confirmation and lens modeling of the complex strong-lens system J0233-0205, in which the deflector consists of a pair of galaxies at zd = 0.790 +/- 0.022, the background source lies at zs = 2.160 +/- 0.002, and the circularized Einstein radius is thetaE = 1.680 arcsec +/- 0.003 arcsec. Our lens modeling requires two lens components, namely two elliptical galaxies with Einstein radii of 0.669 arcsec +/- 0.002 arcsec and 0.735 arcsec +/- 0.002 arcsec, respectively, and a projected separation of 0.513 arcsec (about 3.832 kpc), as well as three source components: two disk galaxies separated by 0.4965 arcsec (about 3.712 kpc), plus a point-like component closely aligned with one of the disks. From a joint lensing and stellar-population analysis, we infer a total stellar mass within the critical curve of the lens pair of (1.956 +/- 0.418) x 10^11 solar masses and a total enclosed mass of (1.107 +/- 0.008) x 10^12 solar masses, corresponding to a projected dark-matter fraction of 82 +/- 4%. The stellar masses of the two lens galaxies are (8.548 +/- 2.128) x 10^10 solar masses and (1.525 +/- 0.295) x 10^11 solar masses, implying dark-matter fractions within the z-band effective radius of 57 +/- 11% and 70 +/- 6%, respectively. The small separation of the lens pair, together with its relatively high deflector redshift, makes J0233-0205 a potentially ideal laboratory for probing the mass distribution and dark-matter content of close galaxy pairs. In addition, the two disk galaxies and the associated point-like source make this system valuable for investigating the merger process in the source plane.

Neutrinos are the most abundant fermions in the Universe and influence the formation of large-scale structure through both their non-zero masses and a possible chemical potential which can be described by a single asymmetry parameter. While most previous studies have focused on the impact of the neutrino mass, the role of neutrino asymmetry remains comparatively unexplored. In this work, we investigate how massive neutrinos ($M_{\nu}=0-0.24\,\mathrm{eV}$) with a non-zero asymmetry parameter ($\eta^{2}=0-0.8$) modify the halo mass function (HMF) and halo bias using cosmological N-body simulations with cosmological parameters consistently refitted to CMB observations. We find that at all redshifts, neutrino mass suppresses the abundance of massive halos, whereas neutrino asymmetry enhances the HMF over a broad mass range. At z=0, the abundance of the most massive halos is reduced by up to ~30% in the largest-mass case ($M_{\nu}=0.24\,\mathrm{eV}$), while neutrino asymmetry ($\eta^{2}=0.8$) produces a maximum ~5% enhancement. These effects become increasingly pronounced at higher redshifts: by z=4 and z=9, the enhancement induced by neutrino asymmetry reaches ~25% and ~75%, respectively, while the corresponding suppression due to neutrino mass deepens to below ~40% and ~70% of the massless case. For halo bias, we find that halos with masses above $10^{13.4}\,\mathrm{M_\odot}$ exhibit an enhanced large-scale bias due to neutrino mass, reaching up to ~5% at z=0, while neutrino asymmetry reduces the bias by a few percent on linear scales. These trends strengthen with redshift, with the enhancement and suppression growing to ~15% and ~10% at z=2, respectively. Linear bias models provide an adequate, though not exact, description of halo bias in massive-neutrino cosmologies. Our results demonstrate that halo abundance and clustering offer sensitive probes of both neutrino mass and asymmetry.

We present a census of local active galactic nuclei (AGN) at a redshift of $z\leq0.025$ selected using the high-ionization [Ne v] $\lambda14.32\,\mu$m emission line from the Infrared Database of Extragalactic Observables from Spitzer (IDEOS). We identify 103 sources with detected [Ne v] emission, which we regard as AGN within the volume. This sample represents $\sim18\%$ of the galaxy population within this redshift range, consistent with AGN fractions derived using other selection techniques. We investigate the biases and properties of this [Ne v]-selected AGN sample by comparing it with traditional AGN selection methods based on hard X-ray, optical, and mid-infrared colors. We find that our selection significantly misses AGN with underdeveloped narrow line regions (NLRs), which account for approximately half of the AGN identified by NLR-independent methods. However, approximately $\sim10\%$ of our sample are undetected in optical diagnostics, while $\sim40\%$ are missed by hard X-rays and $\sim70\%$ by infrared continuum. Notably, $\sim15\%$ of our AGN are missed by all classical methods, constituting a population of previously unidentified AGN revealed solely by the [Ne v] emission line. Based on our analysis, we show that this line can efficiently select heavily Compton-thick and host-dominated AGN systems. Our analysis also yields mean bolometric luminosities of $\log(L_{\rm bol}/{\rm erg~s^{-1}})=44.5\pm0.7$, black hole masses of $\log(M_{\rm BH}/M_{\odot})=7.3\pm0.6$, and Eddington ratio of $\lambda_{\rm Edd}=0.15\pm0.11$. Our sample harbors AGN with comparable luminosities but systematically lower-mass black holes accreting at higher Eddington ratios than those in the hard X-ray-selected sample. This suggests that our AGN may represent local analogs of the rapidly growing SMBH population prevalent at cosmic noon.

Determining the Hubble constant tension requires alternative strategies, and multiply imaged quasars, with their intermediate redshifts, can potentially be used in this regard. We provide a currently complete catalogue of spectroscopically confirmed lensed quasars with ESA/{\it Gaia} astrometry and photometry, as well as redshifts and time delays when available. In addition to the improved astrometry, the catalogue increases the number of lensed quasars by a factor of 1.5 (now 364, of which 277 are doubles and 87 are quads or triples) and significantly increases the number of lensing galaxies detected (now 218), which represents a major step forward. Redshifts are provided for 347 quasars and 188 deflectors. A completely new table of time delays, required for estimates of $H_0$, is presented, with 195 time delays from 73 systems. {\it Gaia} absolute astrometry is sub-milliarcsecond and covers the entire sky. Future {\it Gaia} data releases will provide long-term photometry, which should provide many more time delays. The catalogues as presented here enable machine-learning techniques to be trained and tested and subsequently applied to the {\it Gaia} data releases. Finally, we derive simple but homogeneous models of the 18 quadruply imaged quasars for which images of all four components are presented in {\it Gaia} DR3.}

Stars form in cold, dense clouds embedded in galactic discs, but whether their formation is primarily regulated by gravitational collapse, turbulence, or stellar feedback remains unclear. Using four high-resolution dwarf galaxy simulations with and without supernova (SN) feedback and magnetic fields, we test how feedback regulates the supply of dense gas and, consequently, the star formation rate (SFR). Although the SFR does increase when SNe are turned off, this increase is only by a factor of a few. Instead, across all models, the theoretical maximum SFR originally proposed by Zuckerman and Palmer, defined as the ratio of the total dense gas mass to its mean free-fall time (${M_{\rm dense}}/{\tff}$), always exceeds the measured SFR by nearly two orders of magnitude. Moreover, the increase of the SFR in the case without SNe is accompanied by a nearly corresponding increase of the total dense gas mass ($M_{\rm dense}$), such that the dense-gas depletion time, $\tau \equiv {\rm SFR}/M_{\rm dense}$, decreases by only $\sim 33\%$ in the hydrodynamical case and by about 55\% in the magnetohydrodynamical models. This indicates that SN feedback does not primarily act by slowing the collapse of dense gas, but instead by limiting how much diffuse gas can be converted into dense gas. Our results suggest that the main contribution to the regulation of the SFR, at least in dwarf galaxies, may arise from stabilization by galactic rotation, rather than by SN feedback.

The blazar PKS 0903-57 has recently been classified as a flat spectrum radio quasar at a redshift of $z=0.2621$. In March and April 2020, Fermi-LAT and AGILE reported tremendous activity in high-energy $\gamma$ rays with the flux increasing by $\sim$2 orders of magnitude compared to quiescence. The flare was observed with H.E.S.S. in very-high-energy $\gamma$ rays for six nights with a total observation time of 13.1 h, resulting in the discovery of PKS 0903-57 in this energy band with an average flux of $1.5\times 10^{-10}$ ph cm$^{-2}$s$^{-1}$ above an energy threshold of $\sim 180$ GeV corresponding to $60\%$ of the Crab Nebula flux above the same threshold. The very-high-energy $\gamma$-ray flux was strongly variable. X-ray and optical data were collected with Swift and ATOM, and also indicate significant variability. The observed multiwavelength flux and spectral variability during the H.E.S.S. observation window suggest variability time scales on the order of a few hours and reveal complex correlation patterns. The lack of absorption beyond that of the extragalactic background light in the $\gamma$-ray domain suggests that the emission region was located outside of the broad-line region. A leptonic one-zone modeling of the six H.E.S.S. observation nights using the dusty torus as seed photons for the inverse-Compton scattering, results in a low magnetization of the emission region. This implies that shock acceleration is likely the main driver during the event.

Vertical gas flows, such as winds and meridional circulations, are natural outcomes of protoplanetary disk processes and play a critical role in the earliest stages of planet formation. We analyze vertical gas motions in 14 disks as part of the exoALMA Large Program, focusing on the 12CO J=3-2 and 13CO J=3-2 emission lines. Using discminer to model the Keplerian velocity field, we extract line-of-sight velocity residuals and measure the radial and vertical components of the gas motion. Vertical motions are detected in most disks. Two types of patterns emerge: (1) oscillatory up/down flows, likely linked to instabilities, and (2) transitions from downward to upward motions that we interpret as the base of a disk wind. In most cases, the velocity amplitudes are of a few tens of m/s. Two disks, however, MWC758 and CQ Tau, show two spiral velocity features in their residual maps, red- and blue-shifted, which we interpret as vertical velocities reaching up to 350 m/s (0.7 Cs), consistent with gas motion in eccentric disks. Fast upward motions (up to 500 m/s; 1.8 Cs) is also detected in the outer disk of MWC758. Synthetic observations from (magneto)hydrodynamic simulations validate the reliability of our method. Although strong molecular winds appear to be relatively rare in 12CO and 13CO, our study shows that, when traced by deep high spectral resolution line data, protoplanetary disks exhibit ubiquitous vertical flows. However, their overall velocity structure is highly complex, preventing to identify a coherent, dominant physical mechanism driving the vertical motions across all disks, thus requiring further theoretical investigation.

Filaments are ubiquitous throughout the Galaxy. Massive star formation is often observed in hub-filament systems, where multiple filaments appear to be interconnected and merging. Filament-filament collisions are therefore a likely triggering mechanism for massive star formation. We derive basic physical properties of filament-filament collisions, such as the collision cross section (CCS), the hub mass, and its mass function, based on a simple cylindrical filament model. We assume a cylindrical filament with length $2p$, full width $2q$, and line-mass $\lambda_0$, and consider the CCS between two identical filaments. The collision is specified by three vectors: the directions of the colliding filaments ($n_1$ and $n_2$) and the direction of the relative velocity between the two filaments ($n_v=v/|v|$). For the thin filament, $p\gg q$, the CCS is expressed as $S=4p^2|n'_1\times n'_2|$, where $n'_1$ and $n'_2$ represent the directional vectors projected onto a plane perpendicular to the relative velocity $n_v$. As the angle between $n'_1$ and $n'_2$ becomes smaller, the cross section proportional to $p\cdot q$ becomes relatively important. We propose a simple model in which the hub mass is estimated by the overlapping portion of the two colliding filaments. The hub mass function is derived using the CCSs and the geometrically estimated overlapping mass. When the directions and relative velocities of the filaments are isotropically distributed, the mass function expected from a single species of filaments fits well to a power law and the power exponent is $\gamma_M\simeq -2.96$ ~ $-3.78$. The power exponent of the global hub mass function is the same as that of the line-mass distribution function, $\gamma_\lambda\simeq -1.5$. This means that a massive hub is formed by the collision of two massive filaments.

It is predicted that half or more of all cataclysmic variables (CVs) should have evolved past the period minimum and now exist as so-called "period bouncers" where a white dwarf should be accreting from a Roche-lobe filling substellar companion. However, this prediction stands in stark contrast to observations, where only a few per cent of CVs are found in this evolutionary phase. A potential solution to this discrepancy is that a magnetic field emerges from within the white dwarf after the system has reached the period minimum. The transfer of angular momentum from the spin of the white dwarf into the orbit then pushes the two stars apart, detaching them for potentially billions of years. Here we present the discovery of ZTF J021804.16+071152.93, a detached $0.69\pm0.01 M_{\odot}$, 19 MG magnetic white dwarf plus $37\pm5 M_\mathrm{Jup}$ brown dwarf binary with an orbital period of 1.7 hours. The kinematics of the system indicate that it is a high probability member of the galactic thick disk. However, this strongly disagrees with the much younger age of the system obtained from the white dwarf parameters, implying that the system may have been accreting in the past. This system is therefore consistent with having detached as a result of the emergence of the magnetic field of the white dwarf when the system was still mass transferring, and may represent the ultimate fate for many (perhaps even most) CVs.

Supermassive black hole (SMBH) masses can be measured using molecular gas kinematics. Here we present high angular resolution ($0.12$ arcsec or $\approx11$ pc) Atacama Large Millimeter/submillimeter Array observations of the $^{12}$CO(2-1) line emission of the early-type galaxy NGC 1387. The observations reveal a face-on, regularly-rotating central molecular gas disc with a diameter of $\approx18$ arcsec ($\approx1.7$ kpc) and a central depression slightly larger than the SMBH sphere of influence. We forward model the CO data cube in a Bayesian framework with the \textsc{Kinematic Molecular Simulation} code, and use \textit{Hubble Space Telescope} data to constrain the stellar gravitational potential contribution to the molecular gas kinematics. We infer a SMBH mass of $1.10^{+1.71}_{-0.95}[\text{stat},3\sigma]^{+2.45}_{-1.09}[\text{sys}]\times10^8$ M$_\odot$ and a F160W-filter stellar mass-to-light ratio of $0.90^{+0.44}_{-0.35}[\text{stat}, 3\sigma]^{+0.46}_{-0.36}[\text{sys}]$ M$_\odot$/L$_{\odot,\text{F160W}}$. This SMBH mass is consistent with the SMBH mass -- stellar velocity dispersion relation.

Our study seeks to address the debate over the spin of MAXI J1820+070 through broadband spectral modeling of NuSTAR observations obtained during the soft state. We further compare our results with previous spin estimates and examine the source variability across the soft state. In addition, we investigate the origin of the soft X-ray excess, which we argue does not originate from the plunge region as previously suggested. To further investigate the origin of this excess, we calculate spin-dependent radial disk temperature profiles across all epochs. Our results indicate that the black hole in MAXI J1820+070 is rapidly spinning, with spin $a$ > 0.75, potentially powering the relativistic jets. Our analysis reveals a significant decline in the inner disk temperature midway through the soft state, accompanied by a modest increase in the inferred inner disk radius up to 3.5Rg. This behavior is consistent with slight disk truncation, possibly associated with a reduction in gas ionization and nonthermal processes. Furthermore, the soft excess emission below 10 keV is well described by a blackbody component with kT=0.5 keV, approximately 38% cooler than the inner disk. This suggests that the emission may originate from a warm corona layer located beyond 10Rg, analogous to warm Comptonization models proposed to explain the soft X-ray excess in active galactic nuclei.

Photosphere radius expansion (PRE) bursts provide a crucial tool for constraining the mass and radius of neutron stars. In this study, we analyze time-resolved spectroscopic data from XTE J1810-189 in 2008, which exhibit evidence of a PRE event. We report here the possibility of a small-size and low-mass neutron star in XTE J1810-189 with use of the advantage of the direct cooling tail method. We obtained three sets of results, which can be broadly divided into high metal abundance (20 $\rm{Z}_{\odot}$ and 40 $\rm{Z}_{\odot}$), low metal abundance and hydrogen-rich (pure hydrogen, $\rm{Z}_{\odot}$, 0.3 $\rm{Z}_{\odot}$, 0.1 $\rm{Z}_{\odot}$, 0.01 $\rm{Z}_{\odot}$), and pure helium. In the high-metallicity scenario, the inferred neutron star mass is $<1.3\,M_{\odot}$ with a radius $<8\,\rm{km}$. In the low-metallicity, hydrogen-rich case, the mass ranges from 0.3 to 2.1 $M_{\odot}$ with radii of 7-13 km. For a pure-helium composition, we find two mass solutions: $1.08_{-0.22}^{+1.32}M_{\odot}$ (with $R>14\,\rm{km}$) and $2.5-2.9\,M_{\odot}$ (above the highest observed neutron star masses). Additionally, we applied the touchdown method combined with an MCMC analysis, the results are consistent with those from the direct cooling tail method, but with a broader range. Our analysis of the time-resolved spectrum of burst suggests a high-metallicity atmosphere, but new observations are required to confirm this result.

Context. Multi-element phased-array radio telescopes use digital beamforming to widen their field-of-view with numerous tied-array beams (TABs). These beams share bandpass variations and radio frequency interference (RFI). Yet, most pulsar and transient pipelines process each beam independently, ignoring shared spatial information. This leads to many RFI-dominated false positives that require extensive later sifting. Aims. We exploit multi-beam spatial information to stabilize bandpasses, suppress red noise and broad-band RFI, and drastically reduce false positives without degrading genuine astrophysical signals. Methods. We derive tied-array gain against residual phase dispersion, showing off-beam sources converge to the incoherent limit. Using chi-squared statistics, we analyze dividing a TAB by a beam-averaged reference and quantify the necessary smoothing. We test these predictions using LOFAR high-band antenna voltages (PSR B0329+54), simulations, and LOTAAS survey data (PSR J0250+5854). Results. Off-beam sources contribute nearly uniform power across beams once primary-beam effects are handled. Dividing by a smoothed multi-beam reference yields flatter dynamic spectra and equal or higher pulse signal-to-noise ratios compared to incoherent subtraction. Applied to LOTAAS data, this "beam flatfielding" cuts single-pulse false triggers by a factor of ~200 while preserving profile morphology and peak S/N. Conclusions. Beam flatfielding is a computationally cheap, simple post-beamforming step. For current and future multi-beam facilities, it provides stable bandpasses, closer-to-Gaussian noise statistics, and drastically fewer false positives, easing downstream classification without sacrificing sensitivity.

The Nancy Grace Roman Space Telescope currently plans to survey nearly 6000 square degrees of the sky, mainly in the High-Latitude Wide-Area Survey (HLWAS) and Galactic Plane Survey (GPS). Although these surveys are optimized for other science, they are also a treasure trove for studying the nearby universe. The foreground of the HLWAS includes 59 known stellar streams, 14 known satellite galaxies, and 9 globular clusters in the Milky Way, and an additional 63 galaxies within 10 Mpc spanning several orders of magnitude in stellar mass. The GPS includes an additional 38 globular clusters in its footprint. We summarize and visualize these populations and discuss some of the relevant characteristics of the planned Roman observations. We also examine the expected astrometric performance of the core surveys based on the anticipated time-baselines between observations, and point out the substantial improvement provided by longer time intervals between repeat observations. In particular, the plan for a 6-month revisit timescale in the HLWAS is a missed opportunity from the perspective of proper motions. These data will nonetheless be a powerful new resource for studying the Milky Way and its neighborhood.

The eruption of Nova Persei 1901 (GK Per) occurred 125 yrs ago; remarkably it still holds major surprises. Using data from the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer (SPHEREx), we find it has a bipolar molecular hydrogen shell. This shell, which has dimensions 18'x10', is co-spatial with the Halpha nebulosity surrounding the nova, which is purported to be an ancient planetary nebula (PN). The shell is detected most strongly in the 0--0 S(9) 4.6947 micron line. A filament of emission in the S(9) 4.6947 micron line is seen 45" SW of GKPer. This coincides, over much of its length, with the site of X-ray and non-thermal radio emission where the 1901 nova ejecta impinges on the ambient medium. We propose that the H_2 emission from the filament arises from the predicted neutral zone between the forward and reverse shocks. Since it is common for bipolar PNe to be accompanied by H_2 envelopes, it ostensibly suggests that the 18'x10' nebulosity is a conventional PN with a luminous, ionizing central source. We show this is not the case, and that the H$\alpha$ nebulosity may be surrounding gas belonging to pre-existing material that was ionized during the 1901 eruption. The ionized gas is presently undergoing recombination on a timescale of ~3000 years, explaining why the nebulosity is still visible.

Turbulence plays an important role in shaping the interstellar medium, and strongly influences star formation. We aim to identify the physical processes capable of sustaining HI turbulence in the solar neighborhood. We compare recent HI line-of-sight velocity observations within a volume of radius 70-500 pc centered on the Sun with a suite of 1 kpc numerical simulations that include two distinct turbulent drivers: (i) supernova (SN) feedback and (ii) imposed large-scale turbulent forcing. For each simulation, we construct synthetic sky maps that closely mimic the observational one, allowing for a consistent comparison between the simulations and the observational data. HI observations show a median velocity dispersion of 11.1 km s-1 in the solar neighborhood. SN-driven simulations systematically underpredict this value, yielding dispersions in the range 4.9-6.7 km s-1. Simulations with strong enough large-scale forcing can reproduce not only the median observed velocity dispersion, but also the observed velocity distribution.

The time since a galaxy first became a satellite is central to understanding how environment drives galaxy evolution, yet it cannot be measured directly. Using the TNG300 and TNG-Cluster simulations, we track satellites from $z=1$ to $z=0$ and derive a simple, redshift-dependent prescription for ${T}_{\rm{inf}}$ based on position in projected phase space and stellar mass, via symbolic regression. The resulting calibration provides continuous, observation-ready estimates of infall time across projected phase space. In projected phase space, ${T}_{\rm{inf}}$ is often well described by two components, and we provide analytic expressions for the corresponding characteristic timescales. This framework can be applied directly to spectroscopic samples to infer environmental histories in galaxy groups and clusters.

We perform for the first time full simulation-based inference on the Lyman-$\alpha$ forest 1D power spectrum. In particular, we consider the prediction of the Lyman-$\alpha$ forest $P_{\rm 1D}(k)$ at $2.0<z<3.5$ from the CAMELS cosmological hydrodynamic simulations run with the IllustrisTNG and SIMBA galaxy formation models. We train a normalizing flow to perform neural posterior estimation of two cosmological parameters ($\Omega_m$ and $\sigma_8$) and four astrophysical parameters parametrizing supernova and AGN feedback. When training and testing the neural network on the same baryon physics model, the posterior distributions of the cosmological parameters are found to be in excellent agreement with the true parameters values (within $10\%$ deviations in $\gtrsim 75\%$ and $\gtrsim 90\%$ of the cases for $\Omega_m$ and $\sigma_8$, and a precision better than $10\%$ in both), while the astrophysical parameters are generally unconstrained due to the limited probed volume. When training on one model and testing on the other (e.g., training on IllustrisTNG and testing on SIMBA, or viceversa), the performance is significantly worse, both in accuracy and in precision, resulting in a $\sim 10\%$ positive bias on the predicted values for $\sigma_8$. We show that a multi-domain training based on the combination of simulations from both models recovers unbiased constraints, offering an effective solution to cope with the complex problem of the lack of convergence in the predictions from different galaxy formation models. This study represents a promising way forward to constrain cosmology and fundamental physics with the Lyman-$\alpha$ forest with artificial intelligence.

Protoplanetary disks are the birthplaces of planetary systems, and deviations from Keplerian rotation imprinted in disk gas kinematics serve as key tracers of physical processes and the presence of protoplanets within disks. Using the the CO (J=3-2) data from the exoALMA Large Program encompassing 15 disks, we constructed two-dimensional (2D) maps of centroid velocity, line width, and peak intensity, and extracted non-Keplerian deviations by subtracting smooth Keplerian models. This paper provides the first systematic and uniform overview of 2D gas substructures across the entire exoALMA sample. We find that all targets exhibit large-scale deviations from smooth Keplerian disks, displaying a variety of morphologies including spiral-like structures, arc- or ring-like features, and patterns indicative of variations in the emitting surface height. Non-axisymmetric spiral-arm features are detected or suggested in five disks (CQ Tau, MWC 758, HD 135344B, HD 34282, and SY Cha), and are preferentially found in Herbig Ae/Fe systems. In contrast, some other sources (J1852, PDS 66, and V4046 Sgr), despite exhibiting noticeable deviations, appear to be dynamically quieter. This 2D atlas suggests that kinematic substructures are ubiquitous in large ($\gtrsim$ 100 au) protoplanetary disks with ages of a few million years, based on the observations obtained with sufficient sensitivity at moderate-to-high spatial resolution of $\sim$20 au and high velocity resolution of $\sim$0.1 km s$^{-1}$.

Extracting cosmological information from Stage IV weak lensing surveys requires non-linear modelling of the matter power spectrum that is accurate across a broad range of scales and redshifts and robust to baryonic feedback. We forecast the application of the two-loop effective field theory of large-scale structure (EFTofLSS) to Roman Space Telescope, carefully considering parameterization, scale cuts, and priors. We develop neural network emulators for the two-loop integrals, allowing rapid evaluation of the likelihood. Weak lensing demands a continuous-in-redshift description of the EFT, potentially introducing tens of nuisance parameters. We address this by calibrating the counterterm redshift evolution against the Euclid Emulator 2 and accounting for the residual freedom in redshift with spline functions. A principal component analysis of the free parameters reduces the dimensionality to a few degrees of freedom that the data can constrain. Next, we calibrate the priors on those degrees of freedom by using a suite of hydrodynamical simulations. We forecast the $S_8$ constraints as a function of scale cuts, showing that the two-loop EFT with Roman cosmic shear provides unbiased $S_8=\sigma_8\sqrt{\Omega_{\rm m}/0.3}$ constraints with relative errors of about $0.9\%$ and $1.4\%$ when allowing for $5\%$ and $1\%$ contamination from ultraviolet modes, respectively. The two-loop EFT improves the scale reach beyond the one-loop EFT and non-linear dark matter-only models when baryonic effects are included. This framework provides a robust path for extracting small-scale information from future cosmic shear data.

Transmission spectroscopy has enabled the analysis of exoplanet atmospheres. However, a major challenge is the noise from host stars, caused by stellar activity such as dark spots and bright plages. This noise can mimic or obscure signals in transmission spectra, complicating the study of exoplanetary atmospheres. We aim to characterize how unocculted stellar spots impact planetary absorption line profiles during transit by analyzing planet-occulted line distortions (POLDs). We used the SOAPv4 tool to simulate transits of a hot Jupiter orbiting a Sun-like star under different spot configurations. We analyze the induced POLDs in the Ca II K, the Na I doublet, and H-alpha lines. Our simulations show that POLDs vary with spot size, position, and stellar rotation. The Na I and Ca II K lines exhibit the strongest distortions, while H-alpha is comparatively less affected. Low-latitude spots and higher values of v sin i enhance both the amplitude and asymmetry of distortions, whereas high-latitude spots have a weaker impact. Larger spots generally lead to more pronounced modifications of line profiles, although their relative effect can decrease due to rotational broadening. Our results show that non-occulted stellar spots imprint structured and line-dependent distortions in high-resolution transmission spectra, with amplitudes and velocity shifts shaped by the combined effects of activity level, stellar rotation, and spot geometry. The projected spot area emerges as the dominant factor controlling the strength of these signatures, while the line response varies, with Ca II K being the most sensitive and H-alpha displaying distinctive asymmetric features. These findings demonstrate that stellar surface heterogeneities can mimic or alter planetary signals, highlighting the importance of detailed modeling for the reliable interpretation of upcoming observations.

We report on a search using the IceCube Neutrino Observatory for MeV neutrinos from compact binary mergers detected through gravitational waves during the LIGO-Virgo-KAGRA (LVK) O1, O2, and O3 observing runs. The search focuses on events involving at least one candidate neutron star, such as binary neutron star (BNS) and neutron star--black hole (NSBH) mergers, which may produce a burst of thermal neutrinos due to the hot and dense conditions created during the merger. We looked for short-time increases in IceCube's detector activity around each gravitational-wave event, using four time windows centered on the merger time. We also performed a binomial test for two populations, those with and without at least one neutron star. No significant excess of neutrinos was found. We set upper limits on the MeV neutrino flux for each event, and we place constraints on MeV neutrino emission from mergers that have at least one neutron star. We showcase upper limits for GW170817, the first confirmed BNS merger, providing one of the strongest limits to date on MeV neutrino emission from such sources.

Gravitational-wave (GW) dark sirens provide an independent probe of the cosmic expansion history. Their cosmological constraining power, however, depends critically on precise luminosity-distance measurements and sky localizations for cross-matching with galaxy catalogs. Multiband GW observations can track GW events across different frequency bands and thus improve both. Motivated by this, we forecast the cosmological potential of intermediate-mass black hole binaries (IMBHBs) observed by a three-band GW detector network composed of Taiji (TJ), the Lunar Gravitational-wave Antenna (LGWA), and the Einstein Telescope (ET). We simulate detectable IMBHB populations and analyze them with a hierarchical Bayesian dark-siren framework that includes galaxy-catalog completeness and redshift uncertainties. We find that the TJ-LGWA-ET network outperforms all two-detector configurations considered here. In the $\Lambda$CDM model, it constrains the Hubble constant and matter density to $\sim 0.12\%$ and $\sim 0.6\%$, respectively. In the $w$CDM model, a 4-year dark-siren sample alone constrains the dark-energy equation-of-state parameter $w$ to $\sim 2.7\%$. Adding baryon acoustic oscillation (BAO) and Type Ia supernova (SNe Ia) data improves the $w$ constraint to $\sim 2.1\%$, slightly better than that from the current CMB+BAO+SNe Ia combination. We also show that the final constraints remain sensitive to IMBHB population assumptions and galaxy-catalog limitations, which highlights the need for deep galaxy surveys with precise redshift measurements.

The presence of asymmetries and substructures in protoplanetary disks, revealed by both dust and gas emission, highlights the potential interplay and the broader connection between chemistry and dynamics in disk evolution. We explore multiple relationships using the nonparametric Kendall-$\tau$ correlation to examine formaldehyde (H$_2$CO) emission with relation to stellar and disk properties for a subset of disks from the exoALMA sample. We also retrieve the H$_2$CO column density and excitation temperature using four transitions, measured in radial bins of 100 au, and quantify the level of asymmetry in the resolved peak intensity of the H$_2$CO emission. From our correlation analysis, we find no correlations with sufficient statistical significance. However, we identify tentative relationships that can be tested with larger samples. In particular, we report a proposed correlation ($2.1\sigma$) between stellar effective temperature and the formaldehyde excitation conditions, suggesting that, to first order, the central star dominates the nature of the H$_2$CO emission over possible dynamical asymmetries traced by dust. Although a correlation with the stellar luminosity was also expected, a larger sample is required to confirm or refute this trend. A possible correlation with spectral type, together with the broad range of H$_2$CO excitation temperatures within the inner 100 au of the studied disks, hint at possible multiple chemical formation pathways for H$_2$CO, including both gas-phase reactions and ice-surface chemistry on dust grains.

The extent to which galactic-scale astrophysical processes conspire with the underlying cosmological model to expel baryons from haloes remains a central question in galaxy formation. We present an analytical model for the gas distribution within and beyond haloes, based on the balance between gravitational collapse, hydrostatic pressure, and cosmic expansion. Our model predicts, from first principles, the halo-centric distance enclosing a baryon mass fraction equal to the cosmic value $f_{\rm b} = \Omega_{\rm b}/\Omega_{\rm m}$ (`closure radius') in an arbitrary $\Lambda$CDM cosmology. We compare the predictions with the results of six variants of the EAGLE cosmological, hydrodynamical simulation, encompassing values of the cosmological constant ranging from 0 to 100 times its observed value in our Universe, $\Lambda_0$. Despite its simplicity, our model exhibits excellent agreement with the simulations for haloes with mass $M_{\rm 200c} > 10^{11} M_\odot$ in the redshift range $0<z<3$, suggesting that it captures the key astrophysical processes and highlighting its robustness to the cosmological parameters. Thus, it provides the first physical explanation for the empirical closure radius-halo mass relation previously observed in simulations. Furthermore, we find that dark energy plays a non-negligible role in baryon evacuation: the simulations reveal that in the fiducial cosmological model, the closure radius at $z<2$ is $\sim 30\%$ larger than in an Einstein-de Sitter universe. In cosmologies with $\Lambda \geq 10 \Lambda_0$, dark energy emerges as the dominant factor in this process -- suggesting that, as our Universe transitions towards $\Lambda$-domination, dark energy eventually becomes the primary driver of baryon evacuation from massive haloes.

Galaxy surveys demand fast large-scale structure forward models that preserve large-scale phases while providing realistic nonlinear morphology at fixed force resolution. Single-step Lagrangian Perturbation Theory (LPT) solvers are efficient, but they typically yield overly diffuse filaments and knots and underpredict small-scale clustering. We introduce Ridged Lagrangian Perturbation Theory (RLPT), a modular two-step scheme: a standard long-range LPT/ALPT transport is followed by a single post-processing Eulerian {ridging} update that reconstructs a short-range, curl-free displacement from the realised density field through a smooth scale separation and a Poisson inversion. This explicit completion layer is inexpensive, preserves the large-scale solution, and provides a small set of transparent parameters to tune the short-range response. We test RLPT against particle-mesh and $N$-body references and find that one additional ridging step systematically improves both nonlinear power and field-level agreement relative to 2LPT/ALPT baselines. Finally, we demonstrate that ridging can be repurposed as a deterministic subgrid relocation model: even when the underlying dark-matter field is only ``good enough'' on the mesh, ridging enables controlled tuning of tracer clustering beyond the nominal resolution, which is particularly relevant for mock-galaxy production and observational systematics sensitive to close pairs.

Context. Cyclotron resonant scattering features (CRSFs) in accreting X-ray pulsars are often difficult to detect, especially when shallow or variable. Recent studies have shown that combining spectral and timing analyses enhances their detectability. Aims. We investigated the evolution of energy-resolved pulse profiles of the X-ray pulsar 4U 1901+03 during its 2019 giant outburst, focusing on the 30-40 keV range where there have been disputed claims of a cyclotron line detection. Methods. We analysed four NuSTAR observations of 4U 1901+03 at different luminosities. We studied energy-resolved pulse profiles using harmonic decomposition, cross-correlation analysis, energy-phase maps, and pulsed-fraction spectra. We also used Bayesian spectral modelling to assess the presence and properties of a cyclotron line. Results. We detected significant spectral-timing variability in the 30-40 keV range, which becomes stronger at lower luminosities. We found a pronounced drop in the pulsed fraction near 35 keV only in the lowest accretion state and in the first harmonic of one intermediate-luminosity observation. Adopting a Bayesian informative approach, we find evidence for a cyclotron line in all examined energy spectra, with an average centroid energy of E_cyc approx 32 keV (varying by only 1.6%), and an anti-correlation between line depth and luminosity. Conclusions. We show that a combined spectral-timing approach is more sensitive than phase-averaged spectroscopy to shallow cyclotron features. The luminosity-dependent evolution of pulse profiles and cyclotron line depth point to a drastic change in the emission geometry and accretion flow structure.

We investigate a specific emergent dark energy scenario, known as critically emergent dark energy (CEDE), in which dark energy is effectively absent in the early Universe and becomes dynamically relevant only after a critical cosmic epoch through a phase transition. We constrain this model using recent cosmological observations, including cosmic microwave background (CMB) data from \emph{Planck} 2018, baryon acoustic oscillation (BAO) measurements from SDSS and DESI DR2, and two independent Type Ia supernova compilations, PantheonPlus and Union3. Our results show that within the CEDE framework a dark energy phase transition is not ruled out. In particular, CMB-only, CMB+SDSS, and CMB+DESI datasets provide evidence for a nonzero transition scale factor and, according to standard statistical indicators such as $\Delta\chi^2$ and Bayesian evidence, can favor CEDE over the $\Lambda$CDM model. At the same time, we find that CEDE does not fully resolve the Hubble constant tension. Overall, our analysis indicates that dark energy models featuring a phase transition remain a viable and phenomenologically interesting extension of the standard cosmological framework. Upcoming high-precision cosmological surveys will be essential to further assess whether such emergent dark energy scenarios represent a genuine departure from $\Lambda$CDM or an effective description of current data.

We present a detailed analysis of the structural parameters of the globular cluster M 15 using g- and i-band photometric data from Pan-STARRS1 DR2. The central coordinates ($X_{\rm{C}}, Y_{\rm{C}}$), ellipticity ($\epsilon$), and position angle (PA) are derived via two independent methods: ellipse fitting of the two-dimensional stellar number isodensity distribution and Markov Chain Monte Carlo (MCMC) sampling. Our analysis of 38 stellar density intervals reveals a cluster center offset by only $4\,.\!\!^{\prime\prime}1\pm9\,.\!\!^{\prime\prime}6$ from the commonly accepted in the literature value, a good agreement on the order of our map resolution. We find a radial variation in the ellipticity, with a mean value of $\epsilon=0.09\pm0.02$ for the inner region ($R\leq4\,.\!\!^{\prime}5$) and $\epsilon=0.04\pm0.02$ for the outer region ($R>4\,.\!\!^{\prime}5$), where the errors correspond to $1\sigma$. The MCMC analysis of 75 datasets yields a mean $\epsilon=0.022\pm0.005$ for the entire cluster. The PA remains constant with increasing distance from the cluster center, $\rm{PA}=44\,.\!\!{\rm{^\circ}}4\pm16\,.\!\!{\rm{^\circ}}2$, and the MCMC method providing a consistent value of $\rm{PA}=46\,.\!\!{\rm{^\circ}}6\pm7\,.\!\!{\rm{^\circ}}1$. Our results are in agreement with some recent studies but challenge others, suggesting that a single $\epsilon$ value may be insufficient to fully characterize the overall oblateness of M 15 due to incompleteness and crowding effects in its core.

The exoALMA large program provided an unprecedented view of the morphology and kinematics of 15 circumstellar disks, offering a biased but homogenous and well-characterized sample for population-level analysis. Continuum observations revealed numerous dust substructures, known to be potential signatures of embedded planets. We analyze the observed dust morphologies with the simulation-based inference tool DBNets2.0, assuming these are due to embedded planets at fixed locations, to infer the system properties. We estimate the putative planet mass, the disk $\alpha$-viscosity, scale-height, and dust Stokes number that would reproduce 19 substructures in 13 of the 15 exoALMA disks. We compare our results with literature estimates derived with different methods, and find good agreement in most cases. We further explore the implications of the inferred disk properties for accretion, showing that for the Herbig stars in our sample, the implied viscous accretion timescales are too long to account for their observed stellar accretion rates. Regarding planet migration, our results favor inward migration, with only three putative planets expected to migrate outward. Finally, we check for correlations of the inferred disk and planet properties with the disks' gas-to-dust mass ratio, non-axisymmetry index, and masses of the gas, dust, and host stars, finding no remarkable trend.

The exoALMA Large Program has revealed a wealth of substructures in the dust and molecular line emission of several protoplanetary discs, suggesting that planet formation may unfold within highly dynamic environments. Using synthetic observations of planet-disc interactions and disc instabilities, we demonstrate how the origin of these substructures can be investigated through a tomographic study of molecular lines, extending the scope of the analysis beyond line-centroid kinematics alone. Our results indicate that with only a few hours of ALMA integration at moderate angular resolution ($0.15''-0.30''$), it is possible to identify the key signatures driven by planets more massive than 0.1% of the stellar mass. These signatures manifest not only as deviations from Keplerian motion but also as localized line broadening, enabling accurate constraints on the orbital radius and azimuthal location of the planets. We further show that a diagnostic based on line skewness in spectrally resolved observations can help distinguish between planetary and instability-driven signatures, owing to the distinct degrees of velocity coherence associated with each mechanism. Finally, we apply this tomographic analysis to exoALMA CO line data for the discs of HD 135344B and MWC 758. In HD 135344B, we identify strongly localized velocity and line-width perturbations, suggesting the possibility of three massive planets embedded in the disc: one at $R=95$ au, exterior to the continuum substructures, and two within dust gaps at $R=41$ au and $R=73$ au. For MWC 758, the dominance of vertical-velocity spirals over localized signatures is consistent with predictions from models of moderate disc eccentricities or warps, potentially induced by a substellar companion in the inner regions of the system.

The unprecedented accuracy of JWST has led to the detection of silicate clouds in exoplanet atmospheres, allowing for the first time to probe cloud formation in extreme environments. While parametrized cloud descriptions can fit these observations, the results do not fully agree with microphysical models. To bridge this gap, we developed Nimbus, a fast microphysical cloud model that can constrain cloud formation processes from observations and utilize Virga, an equilibrium condensation model balancing gravitational settling and diffusion. Using both models, we investigate WASP-107 b, WASP-17 b, VHS-1256 b, and YSES-1 c to determine their cloud structure and constrain cloud formation processes. Our results show that all four planets have cluster-sized silicate particles (r ~ 1 nm) at high altitudes. Within Nimbus and Virga, these particles can only be explained by highly inefficient cloud particle settling (fsed < 0.1) or by inefficient growth rates due to low sticking coefficients (s < 0.0001). Our results also show that the sticking coefficient is directly linked to the vertical extent of clouds and can therefore be constrained using the broad shape of the spectral energy distribution. The sticking coefficients found for VHS-1256 b and YSES-1 c are in agreement with expectations from laboratory experiments under Earth-like conditions (0.01 < s < 0.3). Panchromatic observations were crucial to achieve these constraints. Future cloud studies should therefore aim to combine observational data from 1 micron to 10 micron whenever possible.

The spatial distribution of the chemical reservoirs in protoplanetary disks is key to elucidate the composition of planets, especially habitable ones. However, the partitioning of the main elements among the refractory and volatile phases is still elusive. Key parameters such as the carbon-to-oxygen C/O elemental ratio and the ionization fraction remain poorly constrained, with the latter potentially orders of magnitude lower than in the interstellar medium. Moreover, the thermal structure of the gas is also poorly known, despite its deep influence on gas-phase chemistry. In this context, ortho-to-para ratios could provide selective and sensitive probes. Recent ALMA observations have measured the spatially resolved column density of ortho-and para-H2CO in the transition disk orbiting TW Hya and derived the radial profile of the ortho-to-para ratio. Yet, current disk models do not include the nuclear-spin-resolved chemistry required to interpret these observations. The present work aims to fill this gap, by combining a parametric disk physical model of TW Hya with the UGAN network, updated to include a comprehensive description of the nuclear-spin-resolved chemistry of formaldehyde. This new model reproduces the observed column density of H2CO to within a factor of 2, as well as the measured ortho-to-para ratio which varies from 1.5 in the outer disk to 3 inside 90au. In particular the low value of this ratio beyond 90au is well explained by our model. However, the statistical value of 3 measured below 70au cannot be reproduced, suggesting that additional processes involving ices may be involved. Our parameter space exploration shows that the abundance of H2CO is highly sensitive to the C/O elemental ratio and to the cosmic-ray ionization rate. Future observations of ortho-and para-H2CO, based on well selected rotational transitions, in a large sample of disks, appear highly desirable.

The dynamics of Saturn's satellite system offer a rich framework for studying orbital stability and resonance interactions. Traditional methods for analysing such systems, including Fourier analysis and stability metrics, struggle with the scale and complexity of modern datasets. This study introduces a machine learning-based pipeline for clustering approximately 22,300 simulated satellite orbits, addressing these challenges with advanced feature extraction and dimensionality reduction techniques. The key to this approach is using MiniRocket, which efficiently transforms 400 timesteps into a 9,996-dimensional feature space, capturing intricate temporal patterns. Additional automated feature extraction and dimensionality reduction techniques refine the data, enabling robust clustering analysis. This pipeline reveals stability regions, resonance structures, and other key behaviours in Saturn's satellite system, providing new insights into their long-term dynamical evolution. By integrating computational tools with traditional celestial mechanics techniques, this study offers a scalable and interpretable methodology for analysing large-scale orbital datasets and advancing the exploration of planetary dynamics.

The global 21-cm signal provides a powerful probe of early-Universe astrophysics, but its detection is hindered by Galactic foregrounds that are orders of magnitude brighter than the signal and distortions introduced by beam chromaticity. These challenges require accurate foreground modelling, rigorous Bayesian model comparison, and robust validation frameworks. In this work, we substantially accelerate global 21-cm inference by exploiting GPU architectures, enabling likelihood evaluations to achieve near-constant wall-clock time across a wide range of model dimensionalities and data volumes. Combined with algorithmic parallelisation of Nested Sampling, this reduces the total inference runtime of this work from hundreds of CPU-years to approximately two GPU-days, corresponding to a cost reduction of over two orders of magnitude. Leveraging this capability, we advance the physically motivated forward-modelling approach, in which foregrounds are represented by a discrete set of sky regions by introducing a novel, observation-dependent sky-partitioning scheme that defines regions using the antenna beam-convolved sky power of a given observing window. We show that this scheme improves modelling performance in three ways: firstly, by enforcing a strictly nested region hierarchy that enables clear identification of the Occam penalty in the Bayesian evidence, facilitating principled optimisation of model complexity; secondly, by enabling more accurate recovery of spatially varying spectral indices, with posterior estimates centred within physically plausible ranges; and thirdly, by allowing complex foregrounds to be modelled for robust global 21-cm signal inference using substantially fewer parameters. Overall, this approach achieves validated recovery at lower region counts, corresponding to an approximate 40% reduction in foreground-model dimensionality.

The brief (10 nanoseconds) transient radio emission from cosmic ray air showers carries key information about the energy and mass composition of high energy cosmic rays, but anthropogenic radio frequency interference has historically prevented radio-based cosmic ray studies from being carried out independently from other types of detectors. We describe a cosmic ray detection system for the Owens Valley Radio Observatory Long Wavelength Array that searches for radio emission from cosmic ray air showers without relying on an external trigger, and runs alongside the other observing modes of the array. The OVRO-LWA, located in Eastern California, recently completed an expansion to 352 dual-polarization antennas and new signal processing infrastructure. In order to detect cosmic rays in the presence of radio frequency interference (RFI), initial event classification and RFI rejection is performed on Field Programmable Gate Array boards, which each process a sampled voltage timeseries from both polarizations of a subarray of 32 antennas. Each board uses dedicated RFI veto antennas outside the air shower radio footprint to reject RFI events. We present the trigger design, RFI flagging strategy, and candidate cosmic rays.

We investigate the impact of spatial curvature, $\Omega_k$, and dynamical dark energy on the cosmological constraints of the neutrino mass sum, $\sum m_\nu$. Using a joint analysis of the latest CMB (Planck and ACT DR6), BAO (DESI DR2) and SNe Ia (DESY5 and DES-Dovekie) datasets, we perform an exploration of the neutrino mass parameter space. To mitigate prior-driven biases near the physical boundary, we implement a symmetric extension wrapper that allows for effective negative masses. We find that the inclusion of spatial curvature significantly modifies the posterior distributions, exhibiting a smooth transition across the $\sum m_\nu = 0$ threshold. In the $\Lambda$CDM + $\Omega_k$ + $\sum m_{\nu,\mathrm{eff}}$ framework, we obtain $\sum m_{\nu,\mathrm{eff}} = -0.011^{+0.052}_{-0.050}$, reducing the tension with the terrestrial lower limit of 0.06 eV from $2.59\sigma$ for the $\Lambda$CDM + $\sum m_{\nu,\mathrm{eff}}$ model to $1.17\sigma$. For the most flexible scenario $w_0 w_a$CDM + $\Omega_k$ + $\sum m_{\nu,\mathrm{eff}}$, we find $\sum m_{\nu,\mathrm{eff}} = -0.07 \pm 0.11$ with a tension of $1.13\sigma$, illustrating how the increased parameter freedom notably degrades the precision of the mass estimate compared to simpler extensions. Our results demonstrate that current cosmological bounds on $\sum m_\nu$ are heavily influenced by boundary effects and geometric degeneracies.

Line intensity mapping (LIM) is a technique for producing 3D maps of the Universe by scanning the sky with a spectrometer sensitive to a range of wavelengths corresponding to the redshifted spectral lines of atoms or molecules, such as hydrogen or carbon, commonly found in galaxies and the diffuse media around them. While LIM experiments have successfully detected the 21 cm line of neutral hydrogen, other lines that reveal large-scale structure or astrophysical processes remain undetected. Many LIM experiments are in development or are underway to fill this gap, but will likely suffer from contamination from systematics, like Galactic foregrounds, or noise. Cross-correlation techniques offer the smoothest route for making detections and constraining astrophysical processes in this regime. In this work, we apply three cross-correlation techniques (stacking, the conditional voxel intensity distribution (CVID), and the cross power spectrum) to simulated LIM maps produced using [CII] luminosity models for a pathfinder LIM experiment (EXCLAIM). We find that these cross-correlation techniques allow for mean detection of the target signal line ([CII]) at redshifts 2.5-3.5 at the 4.5$\sigma$, 3.9$\sigma$, and 8.4$\sigma$ level, respectively, and offer moderate constraints on the line emission model. Under a futuristic scenario with reduced noise, the techniques improve substantially, with detections at the 44.0$\sigma$, 24.6$\sigma$, and 34.3$\sigma$ levels and percent-level constraints. Each technique offers unique information, with the strongest constraints achieved by using the three techniques in combination.

We present a new constraint on the effective number of relativistic species in the early universe, $N_{\rm eff}$, by combining recent primordial helium abundance measurements from the Large Binocular Telescope $Y_p$ Project with primordial deuterium abundance data, cosmic microwave background (CMB) observations from $\it{Planck}$, the Atacama Cosmology Telescope, and the South Pole Telescope, and baryon acoustic oscillation (BAO) data from the Dark Energy Spectroscopic Instrument, yielding $N_{\rm eff}=2.990\pm0.070$ (68% C.L.). This is the tightest constraint on $N_{\rm eff}$ to date, and is in excellent agreement with the standard model prediction of $N_{\rm eff}=3.044$. Furthermore, we constrain excess contributions to $N_{\rm eff}$ beyond the three neutrino species, finding $\Delta N_{\rm eff}<0.107$ (95% C.L.). This bound nearly approaches the minimum contribution to $\Delta N_{\rm eff}$ from a light spin-3/2 particle that decoupled at any time after inflation ended. Our baseline analysis does not include large-scale $\it{Planck}$ polarization information, enabling a fully consistent combination of state-of-the-art CMB and BAO measurements. As a byproduct, we show that current $N_{\rm eff}$ bounds are essentially insensitive to the inclusion or exclusion of optical depth constraints inferred from large-scale CMB polarization data, making $N_{\rm eff}$ highly robust in this regard. Our constraints place stringent limits on light particles in the early Universe and on a broad range of models aimed at increasing the CMB-inferred value of the Hubble constant.

Pulsar Timing Arrays have recently reported strong evidence for a stochastic gravitational wave background. In standard analyses, it is modeled through pulsar-dependent Fourier coefficients assumed to follow gaussian statistics, so that the signal is fully characterized by its two-point function. However, if the background arises from a finite population of inspiralling supermassive black hole binaries, non-gaussian features may emerge, making the determination of higher-order correlators essential. In this work, we compute the complete four-point correlator of the stochastic gravitational wave background Fourier coefficients for four arbitrary pulsar positions, identifying it as the leading probe of non-gaussianity. The result separates into a gaussian contribution, proportional to the square of the two-point function, and a genuinely non-gaussian connected component, whose non-trivial angular dependence generalizes the Hellings and Downs correlation to four pulsars. This angular structure depends only on averages of products of antenna pattern functions, and is therefore expected to be independent of the specific physical origin of the background. We further propose to incorporate the four-point correlator into the parameter-estimation pipeline by deriving a marginalized likelihood that perturbatively accounts for non-gaussian effects. Our results provide the theoretical framework to search for non-gaussian features in pulsar timing array data, opening the way to a more complete characterization of gravitational-wave backgrounds.

Antideuterons and antihelium nuclei in the cosmic-ray spectrum have long been considered a smoking gun signature of dark matter annihilation, making the tentative observation of several such events by AMS highly intriguing. Conventional dark matter models, however, can produce only up to O(1) antideuteron events at AMS and are not capable of generating observable fluxes of antihelium. In this letter, we propose a class of models in which dark matter annihilates into particles carrying baryon and lepton number, whose subsequent decays produce enhanced fluxes of antinucleons and antinuclei. Such scenarios are motivated by Grand Unified Theories and can lead to an order-of-magnitude or larger enhancement in the resulting antideuteron and antihelium-3 fluxes, providing a means by which to potentially explain the events reported by the AMS Collaboration.

Accurate predictions for dark matter-electron scattering in solids require an all-electron treatment together with a faithful description of dielectric screening beyond simple approximations. In particular, local field effects, arising from microscopic inhomogeneities of the electronic response, can significantly modify scattering rates across relevant momentum and energy scales. We present an all-electron framework for computing dark matter-electron scattering rates that incorporates dielectric screening at the random-phase approximation (RPA) level, including local field effects. Using crystalline silicon as a benchmark, we show that local field effects play an important role both at large momentum transfers, spanning multiple Brillouin zones, and at low momentum near the plasmon resonance. We compute electron recoil spectra and projected sensitivities for non-relativistic halo dark matter and for boosted dark matter or other dark-sector particles, which are sensitive to the impact of local field effects in these high and low momentum regimes, respectively. We further present RPA dielectric functions including local field effects for Ge, GaAs, SiC, and diamond, enabling a systematic comparison across target materials. These developments are implemented in the open source code QCDark2.

We investigate spherically symmetric and static traversable wormholes supported by exotic matter, focusing on solutions sourced by physically motivated dark matter energy density profiles. Considering the Thomas-Fermi-type distribution, we construct explicit forms of the shape function $b(r)$ and analyze the resulting radial and tangential pressures, carefully addressing the requirements of the flare-out condition at the throat and the absence of horizons. We explore zero-tidal-force configurations as well as inhomogeneous equations of state, demonstrating how appropriate choices of the radial pressure allow for finite and well-behaved redshift functions throughout the spacetime. Boundary conditions at a finite radius are implemented to ensure vanishing energy density and pressures, and asymptotic expansions are derived to characterize the behavior of the metric and matter content near the edge of the dark matter halo. Additionally, we reformulate the Einstein field equations entirely in terms of the energy density, radial and tangential pressures, and their derivatives, providing a framework to analyze the matter distribution independently of the explicit metric functions. Our results offer a systematic methodology to construct physically consistent wormhole geometries supported by realistic dark matter halos, highlighting the intricate interplay between matter profiles, equations of state, and geometric constraints.

We present a cold-fluid analysis of the purely transverse Weibel (current-filamentation) instability across four regimes: non-relativistic (NR) single-species, NR multi-species, relativistic single-species, and relativistic multi-species (electron--positron and electron--proton). Beginning from linearized fluid equations, we derive the dispersion relations in each regime and extract scaling laws for the maximum growth rate $\gamma_{\rm max}$ and characteristic unstable wavenumber $k_{\rm max} = \omega_{pi}/c$. Relativistic corrections suppress $\gamma_{\rm max}$ by up to 40 per cent above $v_0 \approx 0.2c$, peaking near $v_0 \approx 0.9c$. Multi-species effects are significant only for $m_e/m_i \gtrsim 1/500$. For the tabletop laser experiment of Bai et al., this http URL., 16, 3770 (2025), the cold-fluid prediction gives $d_i = c/\omega_{pi} \approx 31.7\,\mu{\rm m}$, within 2 per cent of the measured filament spacing $\lambda_F \approx 31\,\mu{\rm m}$. The saturation field estimate $B_{\rm sat} \approx 2.3\times10^4$ T is an upper bound, consistent with the measured $\approx 5000$ T under kinetic suppression. Two MMS burst-mode bow shock crossings (October 16, 2015 and November 25, 2017) confirm $k_{\rm max} d_i = 1$ from FGM/FPI data. A multi-environment scatter plot spans 21 orders of magnitude in $n_i$, with all points within a factor of 3 of the 1:1 line.

An axion-like inflaton coupled to non-Abelian gauge bosons provides a compelling microphysical framework for warm inflation. Starting even from cold initial conditions, in these systems, sphaleron heating may generate thermal friction sufficient to sustain finite temperatures throughout the inflationary epoch. Insisting on shift-symmetric potentials, in this work we revisit the viability of these scenarios under the designation of Minimal Warm Inflation. We examine both observational constraints and model-building limitations on models with a hierarchy between the decay constants appearing in the friction rate and in the inflaton potential. We conclude that the popular clockwork mechanism cannot generate the required hierarchy; however, partial-wave unitarity bounds admit effective descriptions that remain consistent with observations.

A purely transmissive Huygens' metasurface model under plane-wave illumination is used to derive circuit parameters describing a constituent unit cell, such that diverse refraction angles are attained at two distinct frequency bands. Various levels of accuracy of the circuit description approaching the analytical are possible by constraining certain numbers of parameters. This theoretical study is then tested by calculating the exact formulas of the two representations for the various strategies proposed. By using simulations of a candidate unit-cell, we then examine whether such circuit parameters correspond to rudimentary versions of the geometry of a so-called parallel 'dogbone' structure. A device of this type is intended as dual-band (dichroic), dual-angle beam refractor diverting an incoming beam at different directions in two different bands without reflections.

A Monte Carlo generator of high energy cosmic ray interactions, relying on a very basic and transparent theoretical formalism, in the framework of the Reggeon Field Theory, is presented. The main motivation for our work is to provide a new cosmic ray interaction model characterized by relatively transparent physics, sufficient parameter freedom, and a high computational efficiency, which can be easily managed by external users, including a re-tuning of the model parameters. Such a model can be used for studying potential modifications of the interaction treatment, necessary for describing particular sets of data on extensive air showers initiated by high energy cosmic rays, at a microscopic level, thereby keeping a consistency with general restrictions, like the unitarity, energy-momentum and charge conservation, Lorentz and isospin invariance. Importantly, this should allow one to study a compatibility of such modifications with relevant accelerator data. The model results for particle production and for basic extensive air shower characteristics are presented and discussed.

The Migdal effect in semiconductors, prompt ionization from a primary nuclear scattering event, can be described across all kinematic regimes using an effective field theory that encodes the complex vibrational and electronic degrees of freedom of the crystal in measurable structure factors. Simultaneously, general dark matter-nucleus interactions can be systematically described using non-relativistic effective field theory operators. We combine these two effective field theory frameworks to calculate the Migdal effect in semiconductors for all ten dimension-six non-relativistic operators. From the effective Hamiltonian, we find that DM-nucleus scattering factorizes from the ionization and vibrational excitation signal as it does in the free-atom case. Using data from EDELWEISS that was taken with a germanium detector, we derive new experimental bounds on each operator and compare these limits to other direct-detection constraints in the literature. We find the accessible parameter space to be disfavored by bounds on heavy mediators contained in simple UV completions that generate the effective operators.

Axion dark matter (DM) is predicted to convert into radio waves in neutron star magnetospheres. We assess the detectability of this signal using a 5 m radio telescope to be installed at the Fan Mountain Observatory, operating in the UHF, L- and S-bands from 0.5 to 4~GHz. We demonstrate that such a telescope can search new parameter space for axion-like particles over a broad range from $2\,\mu\text{eV}<m_a<17\,\mu\text{eV}$ for axion-photon couplings $g_{a\gamma\gamma} \gtrsim 2\times 10^{-12}\text{ GeV}^{-1}$ with a three year observing period assuming the standard halo model -- improving neutron star observations by more than an order of magnitude. The search is broadband and is thus complementary to other techniques in the same frequency range. We describe in detail our neutron star population model, noise model, and proposed observing strategy. Most constraining power comes from neutron stars at the Galactic centre, where the smooth DM halo is densest. If a DM spike exists at the Galactic centre, the search is sensitive in the QCD axion model band. UHF and L-band observations (0.5 to 2~GHz) represent the pathfinder phase of a wider program we call ``Axion Search with Telescope for Radio Astronomy'' (ASTRA). Future higher mass searches aimed at discovery potential for the post-inflation axion require further hardware development to cover S, C, X and Ku bands (2 to 18~GHz).

Penrose's weak cosmic censorship conjecture asserts that spacetime singularities produced by gravitational collapse are generically hidden behind event horizons, thus preventing them from causally influencing distant observers and preserving the predictability of the exterior region. In this work, we probe this conjecture in a setup that deliberately violates one of its central assumptions - the dominant energy condition - by considering the spherical collapse of a phantom scalar field with negative energy density. In principle, such a field could produce a Schwarzschild geometry with negative mass and therefore no event horizon. Our aim is to assess whether, once the dominant energy condition is abandoned, the fully coupled evolution of matter and geometry can dynamically generate or expose naked singularities, thereby probing the robustness of cosmic censorship. To this end, we perform high-accuracy numerical relativity simulations based on fourth-order finite-difference schemes. Starting from smooth, asymptotically flat initial data representing regular phantom scalar wave packets, we follow their fully nonlinear evolution through collapse or dispersion. While an ordinary (positive-energy) scalar field exhibits the standard Choptuik critical behavior at the threshold of black-hole formation, the phantom field displays qualitatively different dynamics. For all amplitudes considered, we find no evidence for trapped surfaces, naked singularities, or alternative stationary end states. Instead, the phantom scalar field always disperses, suggesting that cosmic censorship remains dynamically preserved even in the presence of negative-energy matter.