Papers for Wednesday, May 28 2025

In this article, momentum transport generated by the combined effects of pitch-angle diffusion and Background Flow Velocity Inhomogeneities (BFVIs) is proposed to obtain a cosmic rays acceleration mechanism, starting from the well-known focusing equation describing particle diffusion and acceleration. The inhomogeneities of background flow velocity is ubiquitous in astrophysical environment. The isotropic distribution function equation of charged energetic particles is derived, and its solution is obtained, demonstrating the form of momentum power laws of cosmic rays. In addition, if it is assumed that cosmic rays penetrate compressible plasma waves or turbulence, for quasi-steady states, the spectral index $\delta$ of the momentum power law spectrum of cosmic rays is found to be in the range $[-5, -3]$, which includes the observed power law indices of galactic cosmic rays. The results obtained in this article demonstrate that the mechanism proposed in this article, along with shock acceleration, may also contribute to the acceleration of galactic cosmic rays. Furthermore, when momentum convection effect and higher-order momentum derivative terms are considered, the indices of power laws should be smaller than $-5$. This may explain the power laws of solar energetic particle events.

The Solar Maximum Mission of NASA was one of the first satellites with on board digitization of observations. It was launched for the solar maximum of cycle 21 (1980) in order to study the solar activity. It carried many instruments, such as coronagraphs, X and $\gamma$ ray detectors, an Ultra Violet spectrometer and a radiometer. Ground based support was offered by many institutes, such as Paris Meudon observatory under the form of systematic observations or coordinated campaigns with specific instruments. We present here the Meudon Solar Tower (MST) and magnetograph which offered in the eighties a major contribution with observations of velocity and magnetic fields of the photosphere and chromosphere, while SMM was observing the transition region and corona above.

The Cherenkov Telescope Array Observatory (CTAO) will be a ground-based Cherenkov telescope performing wide-sky surveys, ideal for anisotropy studies such as cross-correlations with tracers of the cosmic large-scale structure. Cross-correlations can shed light on high-energy $\gamma$-ray sources and potentially reveal exotic signals from particle dark matter. In this work, we investigate CTAO sensitivity to cross-correlation signals between $\gamma$-ray emission and galaxy distributions. We find that by using dense, low-redshift catalogs like 2MASS, and for integration times around 50 hours, this technique achieves sensitivities to both annihilating and decaying dark matter signals that are competitive with those from dwarf galaxy and cluster analyses.

The formation of stars is governed by the intricate interplay of nonideal magnetohydrodynamic (MHD) effects, gravity, and turbulence. Computational challenges have hindered a comprehensive 3D exploration of this interplay, posing a longstanding challenge in our understanding of clouds and cores. Our objective was to study the spatial features and time evolution of the neutral-ion drift velocity and the mass-to-flux ratio in a 3D nonideal MHD chemo-dynamical simulation of a supercritical turbulent collapsing molecular cloud. The resistivities of the cloud were computed self-consistently from a vast non-equilibrium chemical network containing 115 species. To compute the resistivities we used different mean collisional rates for each charged species in our network. We additionally developed a new generalized method for measuring the true mass-to-flux ratio in 3D simulations. Despite the cloud's turbulent nature, at early times, the neutral-ion drift velocity follows the expected structure from 2D axisymmetric non-ideal MHD simulations with an hourglass magnetic field. At later times, however, the neutral-ion drift velocity becomes increasingly complex, with many vectors pointing outward from the cloud's center. Specifically, we find that the drift velocity above and below the cloud's ``midplane'' is in ``antiphase''. We explain these features on the basis of magnetic helical loops and the correlation of the drift velocity with the magnetic tension force per unit volume. Despite the complex structure of the neutral-ion drift velocity, we demonstrate that, when averaged over a region, the true mass-to-flux ratio monotonically increases as a function of time and decreases as a function of the radius from the center of the cloud. In contrast, the ``observed'' mass-to-flux ratio shows poor correlation with the true mass-to-flux ratio and the density structure of the cloud.

Molecular gas is the key ingredient of the star formation cycle, and tracing its dependencies on other galaxy properties is essential for understanding galaxy evolution. In this work, we explore the relation between the different phases of the interstellar medium (ISM), namely molecular gas, atomic gas, and dust, and galaxy properties using a sample of nearby late-type galaxies. To this goal, we collect CO maps for 121 galaxies from the DustPedia project, ensuring an accurate determination of $M_{H2}$, the global molecular gas mass. We investigate which scaling relations provide the best description of $M_{H2}$, based on the strength of the correlation and its intrinsic dispersion. Commonly used correlations between $M_{H2}$ and star formation rate (SFR) and stellar mass ($M_{\star}$) are affected by large scatter, which accounts for galaxies that are experiencing quenching of their star formation activity. This issue can be partially mitigated by considering a "fundamental plane" of star formation, fitting together $M_{H2}$, $M_{\star}$, and SFR. We confirm previous results from the DustPedia collaboration that the total gas mass has the tightest connection with the dust mass and that the molecular component also establishes a good correlation with dust. Although dust grains are necessary for the formation of hydrogen molecules, the strength of gravitational potential driven by the stellar component plays a key role in driving density enhancements and the atomic-to-molecular phase transition. Eventually, we investigated the correlations between ISM components and monochromatic luminosities at different wavelengths: we proposed mid and far-IR luminosities as reliable proxies of $L^{\prime}_{CO}$ for sources lacking dedicated millimeter observations. Luminosities in mid-IR photometric bands collecting PAH emission can be used to trace molecular gas and dust masses.

We present a spectroscopic sample of 23 broad-line AGNs (BLAGNs) at $3\lesssim z\lesssim 6$ selected using F322W2+F444W NIRCam/WFSS grism spectroscopy of the central 100 ${\rm arcmin^2}$ area of the NEXUS survey. Among these BLAGNs, 15 are classified as Little Red Dots (LRDs) based on their rest-frame UV-optical spectral slopes and compact morphology. The number density of LRDs is $\sim 10^{-5}\,{\rm cMpc^{-3}}$, with a hint of declining towards the lower end of the probed redshift range. These BLAGNs and LRDs span broad H$\alpha$ luminosities of $\sim 10^{42.2}-10^{43.7}\,{\rm erg\,s^{-1}}$, black hole masses of $\sim 10^{6.3}-10^{8.4}\,M_\odot$, and Eddington ratios of $\sim 0.1-1$ (median value 0.4), though the black hole mass and Eddington ratio estimates carry large systematic uncertainties. Half of the LRDs show strong Balmer absorption, suggesting high-density gas surrounding the line-emitting region. We detect extended (hundreds of parsec) rest-frame UV-optical emission from the host galaxy in the majority of these LRDs, which contributes significantly or even dominantly to their total UV emission. This host emission largely accounts for the peculiar UV upturn of the LRD spectral energy distribution. We also measure the small-scale ($\lesssim 1\,{\rm cMpc}$) clustering of these BLAGNs and LRDs by cross-correlating with a photometric galaxy sample. Extrapolating the power-law two-point correlation function model to large linear scales, we infer a linear bias of $3.30_{-2.04}^{+2.88}$ and typical halo masses of a few $\times 10^{11}\,h^{-1}M_\odot$ for BLAGNs at the sample median redshift of $z\sim 4.5$. However, the inferred linear bias and halo masses of LRDs, while formally consistent with those for BLAGNs at $\sim 1.5\sigma$, appear too large to be compatible with their space density, suggesting LRDs may have strong excess clustering on small scales.

Gaussian processes (GPs) have become a common tool in astronomy for analysing time series data, particularly in exoplanet science and stellar astrophysics. However, choosing the appropriate covariance structure for a GP model remains a challenge in many situations, limiting model flexibility and performance. This work provides an introduction to recent advances in GP structure learning methods, which enable the automated discovery of optimal GP kernels directly from the data, with the aim of making these methods more accessible to the astronomical community. We present gallifrey, a JAX-based Python package that implements a sequential Monte Carlo algorithm for Bayesian kernel structure learning. This approach defines a prior distribution over kernel structures and hyperparameters, and efficiently samples the GP posterior distribution using a novel involutive Markov chain Monte Carlo procedure. We applied gallifrey to common astronomical time series tasks, including stellar variability modelling, exoplanet transit modelling, and transmission spectroscopy. We show that this methodology can accurately interpolate and extrapolate stellar variability, recover transit parameters with robust uncertainties, and derive transmission spectra by effectively separating the background from the transit signal. When compared with traditional fixed-kernel approaches, we show that structure learning has advantages in terms of accuracy and uncertainty estimation. Structure learning can enhance the performance of GP regression for astronomical time series modelling. We discuss a road map for algorithmic improvements in terms of scalability to larger datasets, so that the methods presented here can be applied to future stellar and exoplanet missions such as PLATO.

The impact of dark matter-neutrino ($\nu$DM) interactions on cosmological perturbations has regained attention, spurred by indications of non-zero couplings from high-multipole cosmic microwave background data, weak lensing, and Lyman-$\alpha$ observations. We demonstrate that a similar observational preference is obtained if $\nu$DM interactions are primarily enhanced during a specific epoch, $z\sim (10^4-10^5)$, leading to $>3\sigma$ preference for a non-zero interaction in the combined Atacama Cosmology Telescope and cosmic shear data. This redshift-limited enhancement circumvents other cosmological and astrophysical bounds and can be achieved within a neutrino portal dark matter framework incorporating resonantly enhanced scattering rates.

Hydrodynamics is a difficult subject to teach in the classroom because most relevant problems must be solved numerically rather than analytically. While there are numerous public hydrodynamics codes, the complexity of production-level software obscures the underlying physics and can be overwhelming to first-time users. Here we present ULULA, an ultra-lightweight python code to solve hydrodynamics and gravity in 2D. The main goal is for the code to be easy to understand, extend, and experiment with. The simulation framework consists of fewer than 800 active lines of pure python code, but it includes a robust MUSCL-Hancock scheme with exchangeable components such as Riemann solvers, reconstruction schemes, boundary conditions, and equations of state. Numerous well-known hydrodynamics problems are provided and can be run in a few minutes on a laptop. The code is open-source, generously commented, and extensively documented.

Unraveling the multi-phase structure of the diffuse interstellar medium (ISM) as traced by neutral hydrogen (HI) is essential to understanding the lifecycle of the Milky Way. However, HI phase separation is a challenging and under-constrained problem. The neutral gas phase distribution is often inferred from the spectral line structure of HI emission. In this work, we develop a data-driven phase separation method that extracts HI phase structure solely from the spatial morphology of HI emission intensity structures. We combine scattering spectra (SS) statistics with a Gaussian-mixture variational autoencoder (VAE) model to: 1. derive an interpretable statistical model of different HI phases from their multi-scale morphological structures; 2. use this model to decompose the 2D channel maps of GALFA-HI emission in diffuse high latitude ($|b|>30$\degree) regions over narrow velocity channels ($\Delta v=3$ km/s) into cold neutral medium (CNM), warm neutral medium (WNM), and noise components. We integrate our CNM map over velocity channels to compare it to an existing map produced by a spectrum-based method, and find that the two maps are highly correlated, while ours recovers more spatially coherent structures at small scales. Our work illustrates and quantifies a clear physical connection between the HI morphology and HI phase structure, and unlocks a new avenue for improving future phase separation techniques by making use of both HI spectral and spatial information to decompose HI in 3D position-position-velocity (PPV) space. These results are consistent with a physical picture where processes that drive HI phase transitions also shape the morphology of HI gas, imprinting a sparse, filamentary CNM that forms out of a diffuse, extended WNM.

The largest reservoir of carbon in protoplanetary discs is stored in refractory organics, which thermally decompose into the gas-phase at the organics line, well interior to the water iceline. Because this region is so close to the host star, it is often assumed that the released gaseous material is rapidly accreted and plays little role in the evolution of the disc composition. However, laboratory experiments show that the thermal decomposition process is irreversible, breaking macromolecular refractory organics into simpler, volatile carbon-bearing compounds. As a result, unlike the iceline of other volatiles, which traps vapor inwards due to recondensation, the organics line remains permeable, allowing gaseous carbon to diffuse outward without returning to the solid phase. In this paper, we investigate how this process affects the disc composition, particularly the gas-phase C/H and C/O ratios, by incorporating it into a 1D evolution model for gas and solids, and assuming refractory organics dominantly decompose into C$_2$H$_2$. Our results show that this process allows this carbon-rich gas to survive well beyond the organics line (out to $7 \mathrm{~au}$ around a solar-mass star) and for much longer timescales, such that its abundance is increased by an order of magnitude. This has several implications in planet formation, notably by altering how the composition of solids and gas relate, and the fraction of heavy elements available to giant planets. In the framework of our model, refractory organics significantly influence the evolution of the gas-phase C/O ratio, which may help interpreting measurements made with Spitzer and JWST.

We measure the projected two-point correlation functions of emission line galaxies (ELGs) from the Dark Energy Spectroscopic Instrument (DESI) One-Percent Survey and model their dependence on stellar mass and [OII] luminosity. We select $\sim$180,000 ELGs with redshifts of $0.8 < z < 1.6$ and define 27 samples according to cuts in redshift and both galaxy properties. Following a framework that describes the conditional [OII] luminosity-stellar mass distribution as a function of halo mass, we simultaneously model the clustering measurements of all samples at fixed redshift. Based on the modeling result, most ELGs in our samples are classified as central galaxies, residing in halos of a narrow mass range with a typical median of $\sim$10$^{12.2-12.4}$ $h^{-1} M_\odot$. We observe a weak dependence of clustering amplitude on stellar mass, which is reflected in the model constraints and is likely a consequence of the 0.5-dex measurement uncertainty in the stellar mass estimates. The model shows a trend between galaxy bias and [OII] luminosity at high redshift ($1.2 < z < 1.6$) that is otherwise absent at lower redshifts.

Young star clusters can inherit bulk rotation from the molecular clouds from which they have formed. This rotation can affect the long-term evolution of a star cluster and its constituent stellar populations. In this study, we aim to characterize the effects of different degrees of initial rotation on star clusters with primordial binaries. The simulations are performed using NBODY6++GPU. We find that initial rotation strongly affects the early evolution of star clusters. Rapidly rotating clusters show angular momentum transport from the inner parts to the outskirts, resulting in a core collapse. Angular momentum transport is accompanied by a highly elongated bar-like structure morphology. The effects of bulk rotation are reduced on the timescale of two-body relaxation. Rotating and non-rotating clusters experience changes in the direction of angular momentum near the dissolution and early evolution due to the tidal field, respectively. We present synthetic observations of simulated clusters for comparison with future observations in filters of Gaia, CSST, and HST. This work shows the effects of bulk rotation on systems with primordial binaries and could be used for the identification of rotation signatures in observed open clusters.

We present simulated counterparts of the ``Little Red Dot'' (LRD) galaxies observed with JWST, using the large cosmological hydrodynamic simulation, ASTRID. We create mock observations of the galaxies ($5 \leq z \leq 8$) in ASTRID, and find seventeen which fit the color and size criteria of LRDs. These LRDs are galaxies with high stellar masses ($\rm log(M_*/M_{\odot}) \geq 9.7$), and massive black holes ($\rm log(M_{BH}/M_{\odot}) \geq 6.8$). The host galaxies are dense, with stellar half mass radii ($\rm 325\,pc \leq r_{{\rm half},*} \leq 620\,pc$), and dust attenuation in the F444W band above 1.25. Their star formation has been recently quenched. They host relatively bright AGN that are dust-obscured and contribute significantly to the rest-frame optical red slope and have relatively low luminosity in the rest-frame ultraviolet, where the host galaxy's stars are more dominant. These LRDs are in an evolutionary phase of miniquenching that is the result of AGN feedback from their massive black holes. The LRDs in ASTRID are bright with F444W magnitudes of $23.5-25.5$. The less massive and fainter galaxies in ASTRID lack the dust concentration necessary to produce the red slope of an LRD, though this could be an effect of limited resolution. Most of the highest Eddington black holes are not LRDs due to their host galaxies having typical dust levels and relatively high star formation rates accompanying their highly accreting black holes, resulting in their spectra being too flat.

We present an analysis of the hard X-ray emission from the central region of Abell 3667 using deep NuSTAR observations. While previous studies on the nature of the hard X-ray excess have been controversial, our analysis of the central region suggests that the excess is primarily thermal, best described by a two-temperature (2T) model, with the high-temperature component likely arising from merger-induced heating. This interpretation contrasts with some earlier suggestions of non-thermal emission due to inverse Compton scattering of relativistic electrons. Additionally, we set a lower limit on the magnetic field strength of $\sim 0.2 \, \mu$G in the central region, consistent with values found in other dynamically active clusters and compatible with those inferred from equipartition and Faraday rotation measurements. Since our study is focused on the central region of the cluster, further high-resolution observations of the outer regions will be critical to fully disentangle the thermal and non-thermal contributions to the X-ray.

The highly elliptical polar orbit of the Juno mission provides a unique opportunity to simultaneously measure the compression state of Jupiter's magnetosphere and the total power emitted by the planet's ultraviolet aurora, using a single spacecraft. This allows us to study how Jupiter's aurora respond to a compression event. In this paper, we present a case study of an extreme compression event that occurred on December 6-7 2022 when Juno was a distance of 70 R$_{J}$ from Jupiter. This extreme compression was accompanied by a very large increase in the ultraviolet auroral emissions to 12 TW, a factor of six higher than the baseline level. This event coincided with the predicted arrival of a powerful interplanetary shock, which was expected to cause the largest increase in the solar wind dynamic pressure seen thus far during the Juno mission. The simultaneous occurrence of the interplanetary shock, the extreme compression and the bright ultraviolet aurora suggests that in this case, the auroral brightening was caused by the solar wind shock compressing the magnetosphere.

We present a study of equal-mass hyperbolic encounters, embedded in a uniform gaseous medium. Using linear perturbation theory, we calculate the density wakes excited by these perturbers and compute the resulting forces exerted on them by the gas. We compute the changes to orbital energy, orbital angular momentum and apsidal precession across a wide range of eccentrities and pericenter Mach numbers. We identify six distinct classes of hyperbolic orbits, differing through their wake structure and subsequent orbital evolution. We find the gas to always dissipate orbital energy, leading to smaller semi-major axes and higher pericenter Mach numbers. The orbital angular momentum can either increase or decrease, whereas we typically find the orbital eccentricity to be damped, promoting supersonic gas-captures. Additionally, we find that the force exerted by the gas is not strictly frictional -- particularly for asymptotically subsonic trajectories. Therefore, despite the orbit-integrated changes to orbital parameters being similar to those predicted by the \cite{O99} prescription, the time evolution of the density wakes and the instantaneous forces exerted on the perturbers are significantly different.

Understanding the diverse formation and migration pathways that shape exoplanetary systems requires characterizing both their atmospheric properties and their orbital dynamics. A key dynamical diagnostic is the projected spin-orbit angle - the alignment between the stellar spin and the planetary orbit-which provides crucial tests for theoretical models. This angle can be determined using the Rossiter-McLaughlin effect. Although measurements exist for over 200 planets, the overall distribution of these angles is not fully understood, motivating further observations across the full parameter space. We analyze archival HARPS and HARPS-N spectroscopic transit time series of nine gas giant exoplanets on short orbits and one brown dwarf. We derive their projected spin-orbit angle $\lambda$. We find aligned projected orbits for all nine gas giants as well as the brown dwarf. Furthermore, we are able to derive the true spin-orbit angle for the brown dwarf EPIC 219388192b, $\psi = $25$^{+11}_{-14}$ deg. These projected prograde orbits are consistent with quiet disc migration disfavoring violent events exciting the orbits in the history of these systems. Finally, we investigate the current overlap between spin-orbit angle measurements and atmospheric characterization targets. While we find no strong observational biases due to the spin-orbit angle, we note that the majority of planets with atmospheric data still lack spin-orbit measurements. This incompleteness of the dynamical information may limit the interpretation of upcoming atmospheric surveys.

The James Webb Space Telescope (JWST) has started a revolution in exoplanetary science. From studying in exquisite detail the chemical inventories and physical processes in gas giant exoplanets, the structure and chemical diversity of the enigmatic sub-Neptune population to even providing constraints on the atmospheric make-up of rocky exoplanets, the observatory is enabling cutting-edge science that is touching virtually every sub-area in the field. In this review Chapter, we showcase key highlights from exoplanet science being conducted with this state-of-the-art space observatory, which we believe is representative of the transformational science it is producing. One of the key takeaways from these pioneering JWST observations is how they are starting to reshape not only how we think, study and interpret exoplanet observations -- but how they are also reshaping our intuition about our very own Solar System planets.

We present the first asteroseismic analysis of the bright, nearby red giant star, HD145250. We calculate the global seismic quantities of the star from single-sector, 2-minute TESS photometry, and determine its mass and radius to be ~1.4 M$_\odot$ and ~16 R$_\odot$ using asteroseismic scaling relations. Our values agree with published non-seismic mass and radius estimates based on comparisons with stellar evolutionary models.

In the standard model of terrestrial planet formation, planets are formed through giant impacts of planetary embryos after the dispersal of the protoplanetary gas disc. Traditionally, $N$-body simulations have been used to investigate this process. However, they are computationally too expensive to generate sufficient planetary populations for statistical comparisons with observational data. A previous study introduced a semi-analytical model that incorporates the orbital and accretionary evolution of planets due to giant impacts and gravitational scattering. This model succeeded in reproducing the statistical features of planets in $N$-body simulations near 1 au around solar-mass stars. However, this model is not applicable to close-in regions (around 0.1 au) or low-mass stars because the dynamical evolution of planetary systems depends on the orbital radius and stellar mass. This study presents a new semi-analytical model applicable to close-in orbits around stars of various masses, validated through comparison with $N$-body simulations. The model accurately predicts the final distributions of planetary mass, semi-major axis, and eccentricity for the wide ranges of orbital radius, initial planetary mass, and stellar mass, with significantly reduced computation time compared to $N$-body simulations. By integrating this model with other planet-forming processes, a computationally low-cost planetary population synthesis model can be developed.

The relative abundances of exotic environments provides us with (uninformed) bounds on the habitability of those environments relative to our own, on the basis that our presence here is not too atypical. For instance, since red stars outnumber yellow stars 7 to 3, we can infer that red stars must be less than 8.1 times as habitable as yellow, as otherwise our presence around a yellow star would be a statistical outlier at the level of $5\%$. In the multiverse context, the relative abundances of exotic environments can be drastically different from those in our universe, which sometimes allows us to place much stronger bounds on their relative habitability than we would get by restricting our attention to our universe. We apply this reasoning to a variety of different exotic environments: tidally locked planets, binary star systems, icy moons, rogue planets, liquids with properties different from water, and waterworlds. We find that the bounds on the relative habitability of rogue planets and waterworlds are at least an order of magnitude stronger in a multiverse context than from our universe alone. Additionally, the belief that some of water's special properties are essential for life, such as the fact that ice floats and, with some caveats, that it acts as a universal solvent, are incompatible with the multiverse hypothesis. If any of these bounds are found to be violated in the future, the multiverse hypothesis can be falsified to a high degree of confidence.

Future direct imaging space telescopes, such as NASA's Habitable Worlds Observatory (HWO), will be the first capable of both detecting and characterizing terrestrial exoplanets in the habitable zones (HZ) of nearby Sun-like stars. Since this will require a significant amount of time and resources for even a single system or exoplanet, the likelihood that a system will host detectable life should be considered when prioritizing observations. One method of prioritization is to estimate the likelihood that an exoplanet has remained continuously within the HZ long enough for life to emerge and make a detectable impact on the atmosphere. We utilize a Bayesian method to calculate the likelihood that a given orbital radius around a star is currently in the 2 Gyr continuous habitable zone (CHZ$_2$), the approximate time it took life on Earth to significantly oxygenate the atmosphere. We apply this method to the 164 stars in the NASA Exoplanet Exploration Program Mission Star List (EMSL) for HWO, representing a preliminary sample of Sun-like stars with HZs most accessible to a future direct imaging mission. By considering the CHZ$_2$ likelihood at all orbital radii outside a hypothetical inner working angle for HWO, we define a metric for prioritizing targets according to the accessibility and total extent of the CHZ$_2$. We find that the CHZ$_2$ metric peaks between $3-4$ Gyr for late-F and early-G dwarfs, but tentatively determine that stars earlier than $\sim {F3}$ or hotter than $\sim 6600$ K are unlikely to have a CHZ$_2$ at the time of observation.

The detection of helium escaping the atmosphere of exoplanets has revolutionized our understanding of atmospheric escape and exoplanetary evolution. Using high-precision spectroscopic observations from the James Webb Space Telescope (JWST) NIRISS-SOSS mode, we report the detection of significant helium absorption during the pre-transit phase of WASP-107b (17$\sigma$), as well as in the transit and post-transit phases. This unique continuous helium absorption begins approximately 1.5 hours before the planet's ingress and reveals the presence of an extended thermosphere. The observations show a maximum transit depth of 2.395$\% \pm$ 0.01$\%$ near the helium triplet (36$\sigma$; at NIRISS-SOSS resolution $\sim$ 700). Our ellipsoidal model of the planetary thermosphere matches well the measured light curve suggesting an outflow extending to tens of planetary radii. Furthermore, we confidently detect water absorption (log10 H2O=-2.5 $\pm$ 0.6), superimposed with a short-wavelength slope which we attribute to a prominent signature from unocculted stellar spots (5.2$\sigma$), rather than a small-particle haze slope. We place an upper limit on the abundance of K (log10 K$<$-4.86, or K/H$<$ 75$\times$ stellar) at 2$\sigma$, which is consistent with the O/H super-solar metallicity estimate. This study underscores the transformative potential of JWST for tracing atmospheric and mass-loss processes, while offering a benchmark for future studies targeting helium escape and its implications for planetary evolution.

Using high-precision photometric data from TESS, medium resolution spectroscopic data from LAMOST, and long-term eclipse timings, we provide orbital parameters for two early-type detached eclipsing binary systems: TYC 3740-2072-1 and TYC 2888-780-1, and analyze the orbital period variations and evolutionary status of these two targets. TYC 3740-2072-1, with a spectral type of B1V, consists of a $6.914 M_\odot$ subgiant and a $6.233 M_\odot$ main-sequence component. It is expected to evolve into a semi-detached binary, potentially serving as a progenitor of a Type Ia supernova. TYC 2888-780-1 has a spectral type of A3 and consists of two main-sequence components with masses of $1.682 M_\odot$ and $1.673 M_\odot$, respectively. It is undergoing stable evolution. Through eclipse timing analysis, we find that both targets exhibit apsidal motion effects, with observed AM rates of $\dot{\omega}_{obs} = 0.0412~\text{deg}~\text{cycle}^{-1}$ and $\dot{\omega}_{obs} = 0.0205~\text{deg}~\text{cycle}^{-1}$, respectively. Additionally, TYC 2888-780-1 exhibits orbital period variations, which we attribute to the light travel time effect caused by a third body, the minimum mass of this third body is estimated to be $0.598 M_\odot$.

The spectroscopic data from DESI Data Release 1 (DR1) galaxies enables the analysis of 3D clustering by fitting galaxy power spectra and reconstructed correlation functions in redshift space. Given low measurements of the amplitude of structure from cosmic shear at $z\sim1$, redshift space distortions (RSD) + Baryon Acoustic Oscillation (BAO) signals from DESI galaxies combined with weak lensing can break degeneracies and provide a tight alternative constraint on the $z\sim1$ amplitude of structure. In this paper we perform joint analyses that combine full-shape + post-reconstruction information from the DESI DR1 BGS and LRG samples along with angular cross-correlations with Planck PR4 and ACT DR6 CMB lensing maps. We show that adding galaxy-lensing cross-correlations tightens clustering amplitude constraints, improving $\sigma_8$ uncertainties by $\sim 40\%$ over RSD+BAO alone. We also include angular galaxy-galaxy and galaxy-lensing spectra using photometric samples from the DESI Legacy Survey to further improve constraints. Our headline results are $\sigma_8 = 0.803\pm 0.017$, $\Omega_{\rm m} = 0.3037\pm 0.0069$, and $S_8 = 0.808\pm 0.017$. Given DESI's preference for higher $\sigma_8$ compared to lower values from BOSS, we perform a catalog-level comparison of LRG samples from both surveys. We test sensitivity to dark energy assumptions by relaxing our $\Lambda$CDM prior and allowing for evolving dark energy via the $w_0-w_a$ parameterization. We find our $S_8$ constraints to be relatively unchanged despite a $~3.5\sigma$ tension with the cosmological constant model when combining with the Union3 supernova likelihood. Finally we test general relativity (GR) by allowing the gravitational slip parameter ($\gamma$) to vary, and find $\gamma = 1.17\pm0.11$ in mild ($\sim1.5\sigma$) tension with the GR value of $1.0$.

We investigate the star formation and neutral atomic hydrogen (HI) gas properties of galaxies along three large-scale filaments and two galaxy groups in the wide field around the Virgo cluster. Our goal is to understand how galaxies are processed in low-density environments before falling into high-density regions. Combining the spatial distribution of galaxies with multiwavelength colors such as W3-W1, NUV-r, and g-r, we find a predominance of blue galaxies across the structures, indicating normal-to-enhanced star formation, similar to that of isolated galaxies. However, one filament and one group show a significant number of red galaxies (32% and 20%, respectively), suggesting that star formation has been suppressed in low-density environments before reaching high-density regions. Intriguingly, these red galaxies span a wide range of stellar masses, and the presence of red dwarfs support that not only mass but also environment plays an important role in the quenching of star formation in cluster outskirts. One particular filament, potentially connected to Virgo, already has a group of red populations outside Virgo's R_200, making these galaxies good candidates for being "preprocessed" before entering the Virgo cluster. In addition, several galaxies in the filaments and groups possess relatively low HI gas contents, similar to cluster galaxies. However, the overall fraction of HI-deficient galaxies is not as significantly high as the fraction of red galaxies in these structures. This suggests that HI gas properties are less influenced by the environment than star formation properties in low-density regions, possibly due to gas replenishment through accretion.

We present, for the first time, a systematic study of quasar-associated 2175 Å dust absorbers using spectroscopic data from the Sloan Digital Sky Survey (SDSS) Data Release 16 (DR16). By analyzing the optical spectra and multi-band magnitudes of 557,674 quasars in the redshift range of $0.7 \le z \le 2.4$, we identify 843 absorbers that share the same redshifts as quasars and are believed to originate from dust in the quasar nuclei, the host galaxies, or their surrounding environments. These absorbers exhibit weak bump strengths ($A\rm_{bump}=0.49\pm0.15~\mu m^{-1}$) and narrow widths ($\gamma\rm=0.81\pm0.14~\mu m^{-1}$), while their peak positions span a broad range from $x_0 = 4.2$ to $4.84~ \mu m^{-1}$. Their average extinction curves resemble those of the Large Magellanic Cloud (LMC) but exhibit a shallower slope. In broad absorption line (BAL) quasars, the absorption bumps show systematic shifts in peak positions. Although further confirmation is needed, this may suggest environmental differences in dust grain properties. We find a statistically significant negative correlation between bump strength and redshift, suggesting possible evolution in dust properties. These findings highlight the changing composition and physical conditions of dust in quasar environments, likely influenced by factors such as metallicity, radiation fields, and dust processing mechanisms. Future studies incorporating ultraviolet and infrared data will be essential for refining the dust evolution models. Machine learning techniques and high-resolution spectroscopic follow-ups could enhance sample completeness and provide deeper insights into the chemical properties of the dust absorbers.

We introduce a novel mechanism -- Magnetically Arrested Transmutation (MAT) -- to account for the observed absence of ordinary pulsars near the Galactic centre, a longstanding puzzle known as the missing pulsar problem and the over-representation of magnetic white dwarfs in the same region. In this scenario, compact stars capture and accumulate dark matter, eventually forming an endoparasitic black hole (EBH) of initial mass $M_0$ at their core. Although such EBHs generally grow by accreting host matter, we show that sufficiently strong core magnetic fields can establish pressure equilibrium, thereby stalling further accretion and halting the star's transmutation into a black hole. We derive the conditions for MAT to occur, identifying a critical parameter $\beta$, which encapsulates the interplay between the magnetic field strength, host matter density, and EBH mass. For $0 < \beta \leq 4/27$, the growth of the EBH is arrested, limiting its final mass ($M_{\rm f}$) to $M_0 < M_{\rm f} \leq 3/2 M_0$, whereas for $\beta > 4/27$, full transmutation may ensue. This framework offers a plausible unified explanation for the absence of ordinary pulsars and the survival of the magnetar PSR J1745-2900, and the elevated population of magnetic white dwarfs, in the Galactic centre, and hence could be tested and should have implications for understanding dark matter and compact objects.

The Extreme Universe Space Observatory on a Super Pressure Balloon 2 (EUSO-SPB2) is a pathfinder mission toward a space-based observatory such as the Probe of Extreme Multi-Messenger Astrophysics (POEMMA). The aim of POEMMA is the observation of Ultra High Energy COsmic Rays (UHECRs) in order to elucidate their nature and origins and to discover $\gtrsim$ 20 PeV very high energy neutrinos that originate from transient and steady astrophysical sources. EUSO-SPB2 was launched from Wānaka New Zealand on May 13th, 2023 as a NASA Balloon Program Office test flight. The mission goals included making the first near-space altitude observations of the fluorescence emission from UHECR-induced extensive air showers (EASs) and making the first direct Cherenkov light emission from PeV cosmic rays traversing Earth's atmosphere. In addition, a Target of Opportunity program was developed for selecting and scheduling observations of potential neutrino sources as they passed just below the Earth's limb. Although a leaky balloon forced termination over the Pacific Ocean after 37 hours, data was collected to demonstrate the successful commissioning and operation of the instruments. This paper includes a description of the payload and the key instruments, pre-flight instrument characterizations in the lab and in the desert, flight operations and examples of the data collected. The flight was too short to catch a UHECR event via fluorescence, however about 10 candidate EAS events from cosmic rays were recorded via Cherenkov light.

In this manuscript, through applications of TDE (tidal disruption event) expected variability properties, a potential candidate for True type-2 AGN without hidden central broad line regions (=TT2AGN) is reported in the SDSS J233454.07+145712.9 (=SDSS J2334). Through analyzing the 20-years optical light curves of SDSS J2334 from different Sky Survey projects, a TDE is preferred with a $4.7{\rm M_\odot}$ main-sequence star tidally disrupted by the central BH with mass $11.7\times 10^6{\rm M_\odot}$, indicating that central region within distance about 20 light-days to central BH in SDSS J2334 is directly in the line-of-sight. Moreover, AGN activities in SDSS J2334 can be confirmed through applications of BPT diagrams. Meanwhile, comparing virial BH mass determined through assumed broad Balmer emission components and M-sigma expected BH mass by well measured stellar velocity dispersion through stellar absorption features, optical broad emission lines in SDSS J2334 are disfavored with confidence level higher than 6$\sigma$. Therefore, combining the unique properties of the TDE and the spectroscopic results with only narrow emission lines, SDSS J2334 can be well identified as a potential candidate for a TT2AGN. The results indicate the to detect TDE expected flares in normal Type-2 AGN classified by spectroscopic results should be a new practicable method for identifying

We present the discovery of TCP J07222683$+$6220548, a new ultracompact binary system of the AM CVn type. This system was first identified displaying a $\Delta V = 7.6$ mag outburst on 2025-01-20.9416 UTC by the New Milky Way wide-field survey for transients and later independently detected by ASAS-SN and ZTF. The outburst peaked at $V_{\rm max} = 12.45$ and lasted for seven days, followed by a series of rebrightenings. No previous outbursts are found in archival data. Positive superhumps with a period of $0.032546 \pm 0.000084$ d ($46.87 \pm 0.12$ min), barely detectable during the main outburst, became clearly visible during the first rebrightening that lasted from day 18 to day 24 after the initial outburst. No convincing change in the superhump period was detected. Dense time-series photometry follow-up by a pair of 0.5-m INASAN robotic telescopes, together with VSNET and AAVSO observers, was essential for identifying TCP J07222683$+$6220548 as an AM CVn system and triggering confirmation spectroscopy with the 2.5-m CMO SAI telescope. Some outbursting AM CVn systems lacking such detailed follow-up may remain unrecognized among the newly discovered cataclysmic variable candidates.

RR Lyrae stars (RRLs) are old pulsating variables widely used as metallicity tracers due to the correlation between their metal abundances and light curve morphology. With ESA Gaia DR3 providing light curves for about 270,000 RRLs, there is a pressing need for scalable methods to estimate their metallicities from photometric data. We introduce a unified deep learning framework that estimates metallicities for both fundamental-mode (RRab) and first-overtone (RRc) RRLs using Gaia G-band light curves. This approach extends our previous work on RRab stars to include RRc stars, aiming for high predictive accuracy and broad generalization across both pulsation types. The model is based on a Gated Recurrent Unit (GRU) neural network optimized for time-series extrinsic regression. Our pipeline includes preprocessing steps such as phase folding, smoothing, and sample weighting, and uses photometric metallicities from the literature as training targets. The architecture is designed to handle morphological differences between RRab and RRc light curves without requiring separate models. On held-out validation sets, our GRU model achieves strong performance: for RRab stars, MAE = 0.0565 dex, RMSE = 0.0765 dex, R^2 = 0.9401; for RRc stars, MAE = 0.0505 dex, RMSE = 0.0720 dex, R^2 = 0.9625. These results show the effectiveness of deep learning for large-scale photometric metallicity estimation and support its application to studies of stellar populations and Galactic structure.

Context: Studies of the Galactic population of double white dwarfs (DWDs) that would be detectable in gravitational waves by LISA have found differences in the number of predicted detectable DWDs of more than an order of magnitude, depending on the binary stellar evolution model used. Particularly, the binary population synthesis code BPASS predicts 20 to 40 times fewer detectable DWDs than the codes SeBa or BSE, which relates to differing treatments of mass transfer and common-envelope events (CEEs). Aims: We aimed to investigate which of these models are closer to reality by comparing their predictions to the DWDs known from electromagnetic observations. Methods: We compared the DWDs predicted by a BPASS galaxy model and a SeBa galaxy model to a DWD catalogue and the sample of DWDs observed by the Zwicky Transient Facility (ZTF), taking into account the observational limits and biases of the ZTF survey. Results: We found that BPASS underpredicts the number of short-period DWDs by at least an order of magnitude compared to the observations, while the SeBa galaxy model is consistent with the observations for DWDs more distant than 500 pc. These results highlight how LISA's observations of DWDs will provide invaluable information on aspects of stellar evolution such as mass transfer and CEEs, which will allow theoretical models to be better constrained.

Aims. The main goal of this paper is to present an accurate and efficient numerical strategy for solving the radiative transfer problem for polarised radiation in strong resonance lines forming out of local thermodynamic equilibrium, taking angle-dependent (AD) partial frequency redistribution (PRD) effects and J-state interference into account. We consider the polarisation produced both by the Zeeman effect and by the scattering of anisotropic radiation, along with its sensitivity to the Hanle and magneto-optical effects. Methods. We introduce a formalism that allows treating both a two-level and a two-term atom in the presence of arbitrary magnetic and bulk velocity fields. The problem is formulated by treating the population of the lower level/term as a fixed input parameter. This approach makes the problem linear with respect to the radiation field, enabling the application of efficient matrix-free preconditioned iterative methods for its solution. Additionally, the computation of the scattering emissivity in the comoving frame, together with a careful choice of the angular and spectral quadrature nodes, allow us to speed up the calculations by reducing the number of evaluations of the redistribution functions. Results. The proposed solution strategy is applied to synthesise the Stokes profiles of the Mg ii h&k doublet and the H i Ly-{\alpha} line in 1D semi-empirical models. The results demonstrate that the method is both fast and accurate. A comparison with calculations from HanleRT-TIC displays an overall good agreement, thereby validating our solution strategy. Moreover, for the wavelength-integrated polarisation profiles of the H i Ly-{\alpha} line, we find an excellent agreement between the results obtained including PRD effects in their general AD description and those obtained considering the angle-averaged simplifying approximation.

Transformer networks excel in scientific applications. We explore two scenarios in ultra-high-energy cosmic ray simulations to examine what these network architectures learn. First, we investigate the trained positional encodings in air showers which are azimuthally symmetric. Second, we visualize the attention values assigned to cosmic particles originating from a galaxy catalog. In both cases, the Transformers learn plausible, physically meaningful features.

Understanding energy transfer through the chromosphere is paramount to solving the coronal heating problem. We investigated the energy dissipation of acoustic waves in the chromosphere of the quiet Sun using 3D radiative magnetohydrodynamic (rMHD) simulations. We analysed the characteristics of acoustic-wave heating and its dependence on height and magnetic field configuration. We find the typical heights where acoustic waves steepen into shocks and the frequencies and wavenumbers that most efficiently dissipate wave energy through this steepening. We combined a comprehensive large-scale analysis, spanning the entirety of the simulations for several solar hours, with a detailed view of an individual shock. We find that the flux of propagating acoustic waves correlates closely with viscous dissipation in the chromosphere above the temperature minimum. Acoustic waves with frequencies close to the acoustic cut-off frequency can efficiently heat the quiet Sun chromosphere at the plasma-$\beta$ = 1 interface and play an important role in the chromospheric energy balance.

We describe a proposal for the optical design of three dual Fabry-Perot based narrowband filter systems for for the future European Solar Telescope (EST). These are intended to constitute the core elements of three imaging spectropolarimeters, foreseen to become amongst the most important science instruments for EST. The designs proposed here rely heavily on the heritage of CRISP and CHROMIS, developed for the Swedish 1-m Solar Telescope and described in detail in a companion paper (Scharmer et al. 2025, in prep.). The outstanding performance of these systems, and the simplicity of their designs, provide strong support of our proposal to build similar systems for EST. The design concepts involve i) minimising the FPI clear aperture diameter by means of numerical simulations based on constraints on Strehl and spectral resolution set by the EST Science Advisory Group (SAG); ii) a compact telecentric optical design with an optical path length of less than 4.7 m; iii) a straight-through optical system based on lenses and without any folding mirrors; iv) the combination of a high resolution etalon with high reflectivity and a low reflectivity, low resolution etalon, to mitigate the effects of cavity errors (Scharmer 2006, Scharmer et al. 2025); v) flexibility in terms of image scale by simple replacement of the last lens (the camera lens) of the FPI system. We propose to compensate for the focus curve of ESTs Pier Optical Path (POP) by focusing the camera lenses of the FPI systems. s. The proposed systems should offer several advantages over other much more complex systems, including manufacture, alignment, stability, flexibility of changes of image scale, and costs. The underlying design concepts also make the proposed FPI systems robust and highly performing in terms of image quality, overall transmission, and fidelity of the spectral transmission profile.

Compton-thick active galactic nuclei (CT-AGNs), which are defined by column density $\mathrm{N_H} \geqslant 1.5 \times 10^{24} \ \mathrm{cm}^{-2}$, emit feeble X-ray radiation, even undetectable by X-ray instruments. Despite this, the X-ray emissions from CT-AGNs are believed to be a substantial contributor to the cosmic X-ray background (CXB). According to synthesis models of AGNs, CT-AGNs are expected to make up a significant fraction of the AGN population, likely around 30% or more. However, only $\sim$11% of AGNs have been identified as CT-AGNs in the Chandra Deep Field-South (CDFS). To identify hitherto unknown CT-AGNs in the field, we used a Random Forest algorithm for identifying them. First, we build a secure classified subset of 210 AGNs to train and evaluate our algorithm. Our algorithm achieved an accuracy rate of 90% on the test set after training. Then, we applied our algorithm to an additional subset of 254 AGNs, successfully identifying 67 CT-AGNs within this group. This result significantly increased the fraction of CT-AGNs in the CDFS, which is closer to the theoretical predictions of the CXB. Finally, we compared the properties of host galaxies between CT-AGNs and non-CT-AGNs and found that the host galaxies of CT-AGNs exhibit higher levels of star formation activity.

We present the full sample of 76 galaxies in 39 galaxy cluster fields at z=0.04-0.07 observed with VLT/MUSE by the GASP survey. Most of them (64) were observed as possible ram pressure stripped galaxies (stripping candidates) based on optical B-band images, while the remaining 12 were a control sample of both star-forming and passive galaxies. Based on spatially resolved ionized gas and stellar kinematics, we assess the physical origin of the gas asymmetries and find that 89% of the stripping candidates are confirmed by the VLT/MUSE data. In addition, also 3 of the 4 star-forming galaxies in the control sample show signs of ram pressure. These control galaxies display a ring of unusual emission line ratios, which we see also in field galaxies, possibly originating from the interaction with a hotter surrounding medium. The stripped galaxies are classified into various classes corresponding to different degrees of stripping, from weakest stripping to strong and extreme (jellyfish galaxies) stripping, as well as truncated gas disks with gas left only in the galaxy center. Our results show that selecting cluster stripping candidates based on optical imaging yields a sample that is indeed largely dominated by galaxies affected by ram pressure at different stages and stripping strength, though some contamination is present, mostly by tidal processes. Strong ram pressure cases are found in galaxies over the whole range of stellar masses studied (10^9-10^11.5 Msun) both in low-mass and high-mass clusters (cluster velocity dispersions sigma = 500-1100 km/s). We examine the possible connection between the progressive stages of stripping, up to the phase of a truncated gas disk, and the subsequent complete stripping of gas. We discuss the incompleteness intrinsic to this and other methods of selection to obtain a complete census of ram pressure stripping in clusters.

We examine the impact of the magnetic field on Population III star formation by varying the magnetic field strength. We perform simulations with magnetic field strengths ranging from $10^{-20}$ G to $10^{-4}$ G, in addition to a model without a magnetic field. The simulations are run for $>1000-1400$ yr after the first protostar forms. In weak-field models, the surrounding disk fragments, forming multiple protostars, and the magnetic field is amplified by the orbital motion and rotation of these protostars. In the model without a magnetic field, frequent fragmentation occurs, and the most massive protostar reaches $\sim200 M_\odot$. However, in models with a magnetic field, once the magnetic field is amplified, the protostars merge to form a single massive protostar, and no further fragmentation occurs except in the model with the strongest magnetic field. Even after the formation of the single protostar, the magnetic field continues to amplify, leading to the formation of a thick disk supported by magnetic pressure and a global spiral pattern. In models with moderate or strong magnetic fields, a rotating disk can form, but fragmentation does not occur, and a strong magnetic field drives an outflow. However, the range of parameters for both disk formation and outflow driving is very narrow, making their appearance under realistic conditions unlikely. Given the weak magnetic field in the early universe, Population III stars are expected to form as single stars, surrounded by a thick disk with a spiral pattern. Thus, the magnetic field, regardless of its strength, plays a crucial role in Population III star formation.

We apply population synthesis techniques to analyze TYPHOON long slit spectra of the starburst barred spiral galaxy M83. The analysis covers a central square of 5 arcmin side length. We determine the spatial distribution of dust through the analysis of reddening and extinction, together with star formation rates, ages, and metallicities of young and old stellar populations. For the first time, a spatial one-to-one comparison of metallicities derived from full-spectral fitting techniques with those obtained from individual young stellar probes has been carried out. The comparison with blue supergiant stars, young massive star clusters, and super star clusters shows a high degree of concordance when wavelength coverage in the $B$-band is available. The metallicity of the young population is supersolar and does not show a radial metallicity gradient along the investigated part of the disk, in agreement with our chemical evolution model. However, a notable decrease in metallicity is observed in a tightly confined region at the galaxy center, coinciding with circumnuclear orbits. We attribute this to matter infall either from the circumgalactic medium or a dwarf galaxy interloper or, alternatively, to AGN-interrupted chemical evolution. We confirm the presence of a dust cavity with a diameter of 260~pc close to the galaxy center. Dust absorption and molecular CO emission are spatially well correlated. We find an anticorrelation between R$_V$, the ratio of dust attenuation to reddening, and the emission strength of molecular species present in photo-dissociation regions. We confirm our results by using alternative fitting algorithms and stellar libraries.

Transitional millisecond pulsars (tMSPs) in tight binary systems represent an important evolutionary link between low-mass X-ray binaries and radio millisecond pulsars. To date, only three confirmed tMSPs and a few candidates have been discovered. Most of them are gamma-ray sources. For this reason, searching for multiwavelength counterparts to unassociated Fermi gamma-ray sources can help to find new tMSPs. Here we investigate whether the unassociated gamma-ray source 4FGL J1824.2+1231 belongs to the tMSP family. To find the counterpart to 4FGL J1824.2+1231, we used data from SRG/eROSITA and Swift X-ray catalogues, and from different optical catalogues. We also performed time-series photometric optical observations of the source with the 2.1-m telescope of the Observatorio Astronomico Nacional San Pedro Martir, the 1.5-m telescope of the Maidanak Astronomical Observatory and the 1.5-m Russian-Turkish telescope. In addition, we carried out optical spectroscopic observations with the Russian-Turkish telescope and used archival spectroscopic data obtained with the Gemini-North telescope. Within the position error ellipse of 4FGL J1824.2+1231, we found only one X-ray source which coincides with an optical object. We consider it as a likely multiwavelength counterpart to 4FGL J1824.2+1231. The source shows strong optical variability and significant proper motion. The latter strongly implies that this is a Galactic source. Double-peaked H and He emission lines are detected in its spectrum with a flat continuum, as often observed in accretion disks of compact binary systems. The X-ray spectrum is well fitted by a power law with the photon index 1.7. The derived intrinsic X-ray-to-gamma-ray flux ratio is about 0.2. If the X-ray/optical source is the true counterpart to 4FGL J1824.2+1231, then all its properties suggest that it is a tMSP in the subluminous disk state.

We search for signatures of magnetic flux cancellation in a 3D resistive MHD flux-emergence simulation of coronal jets and eruptions in a coronal-hole-like environment. To do this, we analysed the output from a 3D MHD simulation of an emerging twisted horizontal flux tube from the convection zone into the solar atmosphere. The simulation considered the impact of neutral hydrogen on the magnetic induction equation, that is, it employed partially ionised plasma. Standard and blowout jets as well as eruptions were observed during the simulation. We observe clear evidence of magnetic flux cancellation in a short segment along the internal polarity-inversion line (iPIL) of the photospheric Bz during an extended period of the simulation characterised by eruptions and blowout jets. Converging magnetic footpoint motions at ~ 1 km/s carried sheared fields within the magnetic tails of the emerging flux tube towards the iPIL. These fields reconnect at the iPIL and generate concave-upward and slowly rising field lines causing a flux decrease that is associated with magnetic flux cancellation. We show evidence of magnetic flux cancellation in 3D MHD simulations of coronal hole eruptions and jets associated with an emerging twisted flux tube. The magnetic flux cancellation can be traced up to about 520 km above the photosphere and might contribute to the formation of pre-eruptive magnetic flux rope seeds. Although our results are consistent with several basic aspects of magnetic flux-cancellation observations associated with coronal jets, the observations nevertheless also suggest that cancellation involves much larger fractions of the available flux than our numerical simulation. We supply avenues to address this discrepancy in future work.

Simulation-based inference (SBI) enables cosmological parameter estimation when closed-form likelihoods or models are unavailable. However, SBI relies on machine learning for neural compression and density estimation. This requires large training datasets which are prohibitively expensive for high-quality simulations. We overcome this limitation with multifidelity transfer learning, combining less expensive, lower-fidelity simulations with a limited number of high-fidelity simulations. We demonstrate our methodology on dark matter density maps from two separate simulation suites in the hydrodynamical CAMELS Multifield Dataset. Pre-training on dark-matter-only $N$-body simulations reduces the required number of high-fidelity hydrodynamical simulations by a factor between $8$ and $15$, depending on the model complexity, posterior dimensionality, and performance metrics used. By leveraging cheaper simulations, our approach enables performant and accurate inference on high-fidelity models while substantially reducing computational costs.

Extragalactic foregrounds in cosmic microwave background (CMB) observations are both a source of cosmological and astrophysical information and a nuisance to the CMB. Effective field-level modeling that captures their non-Gaussian statistical distributions is increasingly important for optimal information extraction, particularly given the precise and low-noise observations from current and upcoming experiments. We explore the use of Wavelet Flow (WF) models to tackle the novel task of modeling the field-level probability distributions of multi-component CMB secondaries. Specifically, we jointly train correlated CMB lensing convergence ($\kappa$) and cosmic infrared background (CIB) maps with a WF model and obtain a network that statistically recovers the input to high accuracy -- the trained network generates samples of $\kappa$ and CIB fields whose average power spectra are within a few percent of the inputs across all scales, and whose Minkowski functionals are similarly accurate compared to the inputs. Leveraging the multiscale architecture of these models, we fine-tune both the model parameters and the priors at each scale independently, optimizing performance across different resolutions. These results demonstrate that WF models can accurately simulate correlated components of CMB secondaries, supporting improved analysis of cosmological data. Our code and trained models can be found here (this https URL).

We set up and perform collision rate simulations between dark matter in the form of asteroid-mass primordial black holes (PBHs) and white dwarf stars. These encounters trigger prompt detonations and could be the key to solving the ignition mystery of type Ia supernovae. Our framework is flexible enough to cover the full range of progenitor white dwarf masses, host galaxy stellar masses, galactocentric radial offsets, and cosmic time. The rate distribution pattern is consistent with exhaustive literature observational determinations for a slightly extended log-normal PBH mass spectrum. Most strikingly, the so far unexplained brightness distribution comes out without finetuning. We find no severe contradictions, except that the inferred PBH mass scale is unpredicted from first principles.

The dark matter (DM) conundrum is one of the most intriguing due to its resistance in direct detection experiments. In recent years, attempts to identify non-gravitational signatures as the result of DM traversing or accumulating within stars have attracted a lot of attention. These calculations are usually evaluated at the order-of-magnitude level for stellar populations where the DM density is highest, such as galactic centers. However, if the signature implies the destruction of the host star, their population could have been diminished over a Hubble time in the most DM-dense regions, unless replenished by star formation. This circumstance exemplifies the need for galactic star formation history profiles when deriving DM-induced transient rates, in particular for predicting the host-offset distribution. Here, we combine theoretical and empirical scaling relations of galaxy structure, star formation, and stellar initial mass function to construct a simple and efficient framework that permits us to estimate the target population formation rate and mass function within galactocentric radial zones across galaxy stellar masses and cosmic time. In a companion paper, we apply the framework to the hypothesis that DM in the form of primordial black holes accounts for the ignition of normal type Ia supernovae when colliding with white dwarf stars.

We conduct a search for stellar-mass binary black hole mergers in gravitational-wave data collected by the LIGO detectors during the LIGO-Virgo-KAGRA (LVK) third observing run (O3). Our search uses a machine learning (ML) based method, Aframe, an alternative to traditional matched filtering search techniques. The O3 observing run has been analyzed by the LVK collaboration, producing GWTC-3, the most recent catalog installment which has been made publicly available in 2021. Various groups outside the LVK have re-analyzed O3 data using both traditional and ML-based approaches. Here, we identify 38 candidates with probability of astrophysical origin ($p_\mathrm{astro}$) greater than 0.5, which were previously reported in GWTC-3. This is comparable to the number of candidates reported by individual matched-filter searches. In addition, we compare Aframe candidates with catalogs from research groups outside of the LVK, identifying three candidates with $p_\mathrm{astro} > 0.5$. No previously un-reported candidates are identified by Aframe. This work demonstrates that Aframe, and ML based searches more generally, are useful companions to matched filtering pipelines.

High rates of stable mass transfer likely occur for some binary star systems, but the resulting flow of mass and angular momentum (AM) is unclear. We perform hydrodynamical simulations of a polytropic donor star and a point mass secondary to determine the mass, AM, and velocity of gas that escapes the system, and the dependence on binary parameters such as mass ratio. The simulations use an adiabatic equation of state and do not include any radiative cooling or irradiation of the outflow. Mass transfer is initiated by injecting heat into the stellar envelope, causing it to gradually inflate and overflow its Roche lobe. The transferred mass flows into an accretion disk, but soon begins to escape through the outer Lagrange point (L2), with a lesser amount escaping through the L3 point. This creates an equatorially concentrated circumbinary outflow with an opening angle of 10 to 30 degrees with a wind-like density profile $\rho \propto r^{-2}$. We find that the ratios of the specific AM of the outflowing gas over that of the L2 point are approximately {0.95, 0.9, 0.8, 0.65} for binary mass ratios $q$ (accretor/donor) of {0.25, 0.5, 1, 2}. The asymptotic radial velocity of the outflowing gas, in units of the binary orbital velocity, is approximately 0.1 to 0.2 for the same mass ratios, except for $q=0.25$ where it might be higher. This outflow, if ultimately unbound from the binary, may be a source of circumstellar material that will interact with ejecta from a subsequent supernova or stellar merger.

Transit timing variation (TTV) is a useful tool for studying the orbital properties of transiting objects. However, few TTV studies have been done on transiting brown dwarfs (BDs) around solar-type stars. Here we study the long-term TTV of a population of close BD companions around solar-type stars using TESS data. We use the measured orbital period change rate to constrain the tidal interaction strength between the host star and the BD companion and put limits on the destruction timescale of these transiting BDs. However, we find no statistically significant evidence of orbital decay or expansion in our sample based on the current data. This may be due to either poor observational data or inherently weak tidal dissipation. We then perform simulations to investigate future observation strategies for detecting orbital decay of transiting BDs, which show NGTS-7A b, TOI-263~b and LP 261-75 b are the most promising targets in the next few years. Our study demonstrates the potential of TTV technique to probe the formation and evolution of close BD companions around solar-type stars.

The mantle-inner core gravitational (MICG) mode is the free mode axial oscillation between the mantle and inner core sustained by the gravitational torque between their degree 2 order 2 density structures. Here, we investigate how the MICG mode is affected by oscillations of cylindrical surfaces in the fluid outer core in the form of Alfvén waves. The latter are triggered by oscillations of the tangent cylinder (TC) moving jointly with the inner core and propagate away from the rotation axis. We show that the MICG mode remains a distinct normal mode of oscillation of the core-mantle system only when the triggered Alfvén waves are attenuated before they traverse the width of the fluid core. For an internal magnetic field strength of a few mT, as we expect in Earth's core, Alfvén waves can readily traverse the width of the core, and the MICG mode is absorbed into the spectrum of torsional oscillation (TO) modes. The MICG period retains a dynamical influence, acting as a point of resonance for TO modes, and marking the transition from a TO mode in which the motion of the TC (including the inner core) is weakly impacted by gravitational coupling to one in which the oscillating motion of the TC is strongly restricted. Our results imply that the observed 6-year periodic signal in the length of day cannot be interpreted as the signature of the MICG mode and must instead be caused by TO modes, or more generally, by the propagation of Alfvén waves.

Context. Recent spectroscopic measurements have revealed absorption from higher rotational levels in C$_2$ than previous observations. These improvements are accompanied by the availability of updated radiative and collisional data. Aims. We revisit the density and radiation field intensity diagnostics provided by the observations of many rotational levels of inter- stellar C$_2$ and extensive molecular information. Methods. We built an excitation model of C2 without spatial structure, including levels up to J= 34 where updated radiative and collisional excitation data are introduced as well as excitation by chemical formation. Results. We confirm the importance of the recent collisional excitation rate coefficients of C$_2$ by molecular H$_2$. We show that the new higher level observations cannot be explained by the standard balance between collisional excitation and radiative transitions. We propose that chemical excitation at formation provides a plausible mechanism to explain the observed high excitation of C$_2$. In addition, it allows us to lift the degeneracy of the density over radiation field strength parameter in the excitation model. Conclusions. A 0D model remains limited and it is highly desirable to use a full Photon Dominated Region (PDR) model, which includes all excitation processes introduced here and full chemical and thermal balance.

The 12C/13C isotope ratio has been derived towards numerous cold clouds (20-50 K) and a couple protoplanetary disks and exoplanet atmospheres. However, direct measurements of this ratio in the warm gas (>100 K) around young low-mass protostars remain scarce, but are required to study its evolution during star and planet formation. We derived 12C/13C ratios from the isotopologues of the complex organic molecules (COMs) CH3OH and CH3CN in the warm gas towards seven Class 0/I protostellar systems to improve our understanding of the evolution of the 12C/13C ratios during star and planet formation. We used the data that were taken as part of the PRODIGE large program with the NOEMA at 1mm. The emission of CH3OH and CH3CN is spatially unresolved in the PRODIGE data (300au scale). Derived rotational temperatures exceed 100K, telling us that they trace the gas of the hot corino, where CH3CN probes hotter regions than CH3OH on average (290 K versus 180 K). The column density ratios between the 12C and 13C isotopologues, derived from LTE analysis, range from 4 to 30, thus, are significantly lower than the expected local ISM isotope ratio of about 68. Assuming that CH3CN and CH3OH may inherit the 12C/13C ratio from their precursor species, astrochemical models were conducted for the latter and compared with our observational results. We conclude that an enrichment in 13C in COMs at the earliest protostellar stages is likely inherited from the COMs' precursor species, whose 12C/13C ratios are set during the prestellar stage via isotopic exchange reactions. This also implies that low 12C/13C ratios observed at later evolutionary stages could at least partially be inherited. A final conclusion on 12C/13C ratios in protostellar environments requires improved observations to tackle current observational limitations and additional modelling efforts.

We present an algorithmic method for efficiently planning a long-term, large-scale multi-object spectroscopy program. The Sloan Digital Sky Survey V (SDSS-V) Focal Plane System performs multi-object spectroscopy using 500 robotic positioners to place fibers feeding optical and infrared spectrographs across a wide field. SDSS-V uses this system to observe targets throughout the year at two observatories in support of the science goals of its Milky Way Mapper and Black Hole Mapper programs. These science goals require observations of objects over time with preferred temporal spacinges (referred to as "cadences"), which can differ from object to object even in the same area of sky. robostrategy is the software we use to construct our planned observations so that they can best achieve the desired goals given the time available as a function of sky brightness and local sidereal time, and to assign fibers to targets during specific observations. We use linear programming techniques to seek optimal allocations of time under the constraints given. We present the methods and example results obtained with this software.

We present the discovery of a peculiar central stellar structure in the collisional ring galaxy AM0644-741 using HST imaging and MUSE integral field unit (IFU) data. We identified two Sérsic components with a Sérsic index of 1.72 (inner part) and 1.11 (outer part) in the HST F814W band optical image using \textsc{Galfit}. We utilized the MUSE data cube to construct stellar line of sight velocity (V$_{\rm LOS}$), velocity dispersion ($\sigma_{\rm LOS}$), h$_3$ \& h$_4$ velocity moments, and stellar population age maps using the \textsc{GIST} pipeline for further investigating both Sérsic components, which have a difference of $\sim$ 60 degrees in their position angle. The inner component, with an effective radius $\sim$1 kpc, shows a strong anticorrelation between V$_{\rm LOS}$/$\sigma_{\rm LOS}$ and h$_3$, indicating the presence of a rotating stellar structure. In addition, the inner component also shows a relatively higher velocity dispersion (average values reaching up to $\sim$240 km sec$^{-1}$) along with disky isophotes and stronger Mg~$b$ line strength, which all together highlight a peculiar dynamical state of AM0644-741's central region. Our analysis suggests that the recent encounter has had a smaller impact on the stellar orbits within the inner component. In contrast, it has specifically affected the stellar orbits of the progenitor's outer disk when forming the star-forming ring. The BPT analysis of the unresolved nuclear source shows a LINER-type ionization, hinting at AGN activity in the galaxy. Our study projects the dynamical evolution of collisional systems and provides scope for simulations to explore the central region in greater detail.

Kepler-51 is a 500 Myr G dwarf hosting three "super-puffs" and one low-mass non-transiting planet. Kepler-51d, the coolest (T_eq ~ 350 K) transiting planet in this system, is also one of the lowest density super-puffs known to date (rho_p = 0.038 +/- 0.009 g/cm^3). With a planetary mass of Mp = 5.6 +/- 1.2 Earth masses and a radius of Rp = 9.32 +/- 0.18 Earth radii, the observed properties of this planet are not readily explained by most planet formation theories. Hypotheses explaining Kepler-51d's low density range from a substantial H/He envelope comprising more than 30% of its mass, to a high-altitude haze layer, to a tilted ring system. To test these hypotheses, we present the NIRSpec-PRISM 0.6-5.3 micron transmission spectrum of Kepler-51d observed by the James Webb Space Telescope. We find a spectrum best fit by a sloped line covering the entire wavelength range. Based on forward modeling and atmosphere retrievals, Kepler-51d likely possesses a low-metallicity atmosphere with high-altitude hazes of submicron particle sizes spanning pressures of 1-100 microbars. However, the spectrum could also be explained by a tilted ring with an estimated lifetime on the order of ~0.1 Myr. We also investigate the stellar activity of this young Sun-like star, extracting a spot temperature significantly hotter than sunspots and spot covering fractions on the order of 0.1-10%, depending on the assumed spot parameters.

Recent studies indicate that the physical properties of eclipsing binaries can be extracted from the derivatives of their light curves. A classification scheme for the derivatives of light curves would be helpful for identifying key characteristics of eclipsing binaries. In this study, we propose a new classification method for the light curves of overcontact eclipsing binaries by using their derivatives. We synthesized 89,670 sample light curves of overcontact binaries and categorized them into five types on the basis of their first to fourth derivatives. For each type, we examined the statistical distributions of four parameters: the mass ratio, orbital inclination, fill-out factor, and eclipse obscuration. Their distributions demonstrated that parameter values exhibit certain trends depending on the classified types. With the proposed classification method, general properties of overcontact binaries can be understood, providing a foundation for further detailed analysis.

The gravitational waves emitted by massive black hole binaries can be affected by a variety of environmental effects, which, if detected, could inform astrophysics and cosmology. We here study how gravitational waves emitted by black holes in quasi-circular orbits are affected by the presence of an ultra-light, vector-field, dark-matter environment that is minimally coupled to the binary. This dark-matter environment induces oscillatory gravitational potentials that perturb the orbit of the binary, leaving an imprint in the binary's binding energy, and thus, on the gravitational waves emitted. We here compute the effect of this environment on the gravitational-wave phase using the stationary-phase approximation within the post-Newtonian formalism. We then perform a Fisher analysis to estimate the detectability of this environmental effect with a four-year LISA observation, focusing on vector fields with ultra-light masses in the $(10^{-19}, 10^{-16}) \; \rm{eV}$ range. We conclude that the observation of such gravitational waves with space-borne interferometers, like LISA, could yield a measurement or constraint on local, vector dark-matter environments, provided the dark-matter density is larger than roughly $10^{14} \rm{M}_\odot/{\rm{pc}}^3$.

In the study of optical transients, parameter inference is the process of extracting physical information, i.e. constraints on the source's characteristics, by comparing the observed lightcurves to the predictions of different models and finding the model and parameter combination that make the closest match. In the developing field of the study of kilonovae (KNe), systematic uncertainties in modelling are still very large, and many models struggle to fit satisfactorily the whole multi-wavelength dataset of the AT2017gfo kilonova, associated to the Binary Neutron Star (BNS) merger GW170817. In a multi-messenger context, we sometime observe tensions between KN-only inference results and constraints from other messengers. In order to discuss the compatibility of KN models with observations and with the information derived from other messengers, we detail the process of Bayesian parameter inference, identifying the many sources of uncertainty embedded in KN analyses. We highlight the systematic error margin hyperparameter $\sigma_{\rm sys}$, which can be exploited as a metric for a model's goodness-of-fit. We then discuss how to assess the performance of parameter inference analyses by quantifying the information gain using the Kullback-Leibler divergence between prior and posterior. Using the example of the Bu2019lm model with the NMMA Bayesian inference framework, we showcase the expected performance that dedicated KN follow-ups with telescope networks could reasonably reach, highlighting the different factors (observational cadence, error margins) that influence such inference performances. We finally apply our KN analysis to the dataset of AT2017gfo to validate our performance predictions and discuss the complementarity of multi-messenger approaches.

Non-minimal couplings between the electromagnetic field strength and the spacetime curvature are part of the effective field theory of gravity and matter. They alter the local propagation of light in a significant way if the ratio of spacetime curvature to the non-minimal coupling is of order one. Spacetime curvature can become appreciable around black holes, and yet the effect of non-minimal couplings on electromagnetic observations of black holes remains underexplored. A particular feature of the non-minimal coupling between the electromagnetic field-strength and the Riemann tensor is that it generates two distinct photon rings for different polarizations. Working within the paradigm of lensing bands, and focusing on the $n=1$ photon ring, we investigate which values of the non-minimal coupling can be robustly excluded.

The chemical evolution of the inner regions of protoplanetary discs is a complex process. Several factors influence it, one being the inward drift and evaporation of volatile-rich pebbles. During the disc's evolution, its inner part is first enriched with evaporating water-ice, resulting in a low C/O ratio. Afterwards, C-rich gas from the outer disc is transported inwards. Consequently, the C/O ratio of the inner disc increases again after 2 Myr. Previously, we studied how internal photoevaporation influences these processes in discs around Sun-like stars. We now extend our study to lower-mass stars, where the time evolution of the disc's C/O ratio is different due to the closer-in position of the evaporation fronts and differences in disc mass, size and structure. Our simulations are carried out with the semi-analytical 1D disc model chemcomp, which includes viscous evolution and heating, pebble growth and drift, pebble evaporation and condensation, as well as a simple chemical partitioning model. We show that internal photoevaporation plays a major role in the evolution of protoplanetary discs: As for Sun-like stars, photoevaporation opens a gap, which stops inward drifting pebbles. In addition, volatile-rich gas from the outer disc is carried away by the photoevaporative winds. Consequently, the C/O ratio in the inner disc remains low, contradicting observations of discs around low-mass stars. Our model implies that young inner discs (< 2 Myr) should be O-rich and C-poor, while older discs (> 2 Myr) should be C-rich. The survival of discs to this age can be attributed to lower photoevaporation rates, which either originate from a large spread of observed X-ray luminosities or from the photoevaporation model used here, which likely overestimates the photoevaporation efficiency. A reduction of the latter brings the calculated elemental abundances into better agreement with observations.

PSR B1259-63/LS 2883 is a well-studied gamma-ray binary hosting a pulsar in a 3.4-year eccentric orbit around a Be-type star. Its non-thermal emission spans from radio to TeV energies, exhibiting a significant increase near the periastron passage. This paper is dedicated to the analysis of INTEGRAL observations of the system following its last periastron passage in June 2024. We aim to study the spectral evolution of this gamma-ray binary in the soft (0.3-10 keV) and hard (30-300 keV) X-ray energy bands. We performed a joint analysis of the data taken by INTEGRAL/ISGRI in July-August 2024 and quasi-simultaneous Swift/XRT observations. The spectrum of the system in the 0.3-300 keV band is well described by an absorbed power law with a photon index of $\Gamma=1.42\pm 0.03$. We place constraints on potential spectral curvature, limiting the break energy $E_\mathrm{b}>30$ keV for $\Delta\Gamma>0.3$ and cutoff energy $E_\mathrm{cutoff}>150$ keV at 95% confidence level. For one-zone leptonic emission models, these values correspond to electron distribution spectral parameters of $E_\mathrm{b,e}>0.8$ TeV and $E_\mathrm{cutoff,e}>1.7$ TeV, consistent with previous constraints derived by H.E.S.S.

Understanding the evolution of dust in galaxies is crucial because it affects the dynamics and cooling of gas, star formation, and chemical evolution. Recent work on dust removal in galaxies indicates timescales of gigayears, with old stellar populations and AGNs as the primary drivers of this process. However, most statistically significant studies are focused on low redshifts $z < 0.4$. Here, we determine the dust removal timescale in galaxies over a wide range of redshifts, up to $z \sim 5$. We use publicly available catalogue data of infrared-selected galaxies, observed by \textit{Herschel}. Using the inferred dust masses, stellar masses, and stellar ages, we calculate the dust removal timescale in a sample of more than 120,000 galaxies. We find that, with increasing redshift, the dust removal timescale decreases from 1.8 Gyr at redshift $z \sim 0.05$ to less than 500\,Myr at $z > 3$. Galaxies at higher redshifts undergo more efficient dust removal than galaxies at lower redshift, likely driven by AGN activity, supernova shocks, and astration. These findings indicate that dust removal evolves over cosmic time, reflecting the changing mechanisms regulating dust content of galaxies as the Universe evolves.

Detailed understanding and suppression of backgrounds are among the key challenges faced by Coherent Elastic Neutrino-Nucleus Scattering (CE\ensuremath{\nu}NS) experiments. The sensitivity of these experiments is largely determined by the background levels arising from various sources. Above-ground and shallow-overburden neutrino experiments typically employ passive shielding, primarily composed of lead (Pb), to suppress environmental $\gamma$ background. However, such shielding can introduce additional backgrounds that are particularly challenging for CE\ensuremath{\nu}NS experiments. These backgrounds arise mainly from $\gamma$ and neutrons produced by cosmic muon interactions in the shielding, and their contribution can become significant depending on the amount of Pb shielding used. In the current work, we measure the yield of secondary particles originating from Pb as a result of high-energy cosmic muon interaction, using a high-purity germanium (HPGe) detector and plastic scintillators. A time-coincidence technique is used to identify and reject these secondary background events from the experimental data. The obtained mean characteristic time of these residual background events is 11 $\pm$ 4 $\mu$s, which is consistent with the Geant4-based MC simulation result of 11 $\pm$ 1 $\mu$s. The measured efficiency-corrected rate of muon-induced events in the HPGe detector is 34 $\pm$ 1 (stat.) $\pm$ 3 (sys.) day$^{-1}$kg$^{-1}$ within the energy range of 30 keV to 2000 keV. The yield of muon-induced secondary backgrounds in 10 cm thick Pb shielding is evaluated to be $(11 \pm 1 (\text{stat.}) \pm 1 (\text{sys.}))$ secondary events$\thinspace\text{kg}^{-1}\thinspace\mathrm{m^{-2}}\thinspace\text{muon}^{-1}$ at sea level.

We revisit perturbative unitarity in scalar field inflation with a nonminimal coupling, with Higgs inflation serving as the most prominent example. Although such models are phenomenologically successful, it is critical to examine whether or not unitarity violations spoil their theoretical self-consistency. The analysis of these issues has so far typically relied on order-of-magnitude estimates of scattering amplitudes, which are appropriate for generic parameters. It is not evident that these methods apply to scenarios relying on a near-critical inflationary potential, for which an interplay of both small scalar self-couplings and nonminimal couplings could partially alleviate the unitarity issues. To allow for an exploration of this possibility, we consider the full $S$-matrix for the relevant scattering processes, taking into account important phase space volume factors, leading to a precise evaluation of the cut-off scale. In the single-field case, we demonstrate that near-criticality raises the cut-off scale considerably, compared to previous estimates. In the multifield case, momentum-dependent self-interactions in the kinetic sector lower the cut-off compared to the single-field case to a value comparable to but slightly larger than previous estimates. We carefully study both the single-field and multifield cases in metric and metric-affine (Palatini) formulations of gravity, as well as introduce a new phenomenologically viable model with a canonical kinetic term and a significantly raised cut-off, and discuss the importance of background field effects.

In neutrino-dense astrophysical environments, these particles exchange flavor through a coherent weak field, forming a collisionless neutrino plasma with collective flavor dynamics. Instabilities, which grow and affect the environment, may arise from neutrino-neutrino refraction alone (fast limit), vacuum energy splittings caused by masses (slow limit), or neutrino-matter scattering (collisional limit). We present a comprehensive analytical description of the dispersion relation governing these unstable modes. Treating vacuum energy splittings and collision rates as small perturbations, we construct a unified framework for fast, slow, and collisional instabilities. We classify modes into gapped, where collective excitations are already present in the fast limit but rendered unstable by slow or collisional effects, and gapless, which are purely generated by these effects. For each class, we derive approximate dispersion relations for generic energy and angle distributions, which reveal the order of magnitude of the growth rates and the nature of the instabilities without solving directly the dispersion relation. This approach confirms that slow and collisionally unstable waves generally grow much more slowly than they oscillate. Consequently, the common fast-mode approximation of local evolution within small boxes is unjustified. Even for fast modes, neglecting large-distance propagation of growing waves, as usually done, may be a poor approximation. Our unified framework provides an intuitive understanding of the linear phase of flavor evolution across all regimes and paves the way for a quasi-linear treatment of the instability's nonlinear development.

Scalar-tensor theories with a scalar field coupled to the Gauss-Bonnet invariant can evade no-hair theorems and allow for non-trivial scalar profiles around black holes. This coupling is characterized by a length scale $\lambda$, which, in an effective field theory perspective, sets the threshold below which deviations from General Relativity become significant. LIGO/VIRGO constraints indicate $\lambda$ is small, implying supermassive black holes should not scalarize. However, recent work suggests that scalarization can occur within a narrow window of masses, allowing supermassive black holes to scalarize, while leaving LIGO/VIRGO sources unaffected. We explore the impact of this scenario on the stochastic gravitational wave background recently observed by Pulsar Timing Arrays. We find that scalarization can alter the characteristic strain produced by circularly inspiralling SMBH binaries and that current data shows a marginal preference for a non-zero $\lambda$. However, similar signatures could arise from astrophysical effects such as orbital eccentricity or environmental interactions, emphasizing the need for improved modeling and longer observations to discriminate among the different scenarios.

Large Language Models (LLMs) are being explored for applications in scientific research, including their capabilities to synthesize literature, answer research questions, generate research ideas, and even conduct computational experiments. Ultimately, our goal is for these to help scientists derive novel scientific insights. In many areas of science, such insights often arise from processing and visualizing data to understand its patterns. However, evaluating whether an LLM-mediated scientific workflow produces outputs conveying the correct scientific insights is challenging to evaluate and has not been addressed in past work. We introduce AstroVisBench, the first benchmark for both scientific computing and visualization in the astronomy domain. AstroVisBench judges a language model's ability to both (1) create astronomy-specific workflows to process and analyze data and (2) visualize the results of these workflows through complex plots. Our evaluation of visualizations uses a novel LLM-as-a-judge workflow, which is validated against annotation by five professional astronomers. Using AstroVisBench we present an evaluation of state-of-the-art language models, showing a significant gap in their ability to engage in astronomy research as useful assistants. This evaluation provides a strong end-to-end evaluation for AI scientists that offers a path forward for the development of visualization-based workflows, which are central to a broad range of domains from physics to biology.

We investigate anisotropic compact stars comprising two non-interacting fluids: quark matter and condensed dark matter. Using the MIT Bag model equation of state for quark matter and Bose-Einstein Condensate equation of state for dark matter, we numerically compute interior solutions for those two-fluid component spherical configurations. Varying the initial central density ratio of dark matter to quark matter, we examine how different proportions of these components influence the mass-radius profile, the factor of compactness as well as the quark mass fraction. Recent studies suggest that quark matter may exist in the cores of massive neutron stars, significantly affecting their structure and stability. We calculate the factor of compactness for both negative and positive anisotropy cases explored in this article. Our findings demonstrate that dark matter-admixed quark stars are more compact yet less massive compared to pure quark matter stars, aligning with recent theoretical predictions and gravitational wave observations.

Four-dimensional gravitational theories derived from an infinite sum of Lovelock curvature invariants, combined with a conformal rescaling of the metric, are equivalent to a subclass of shift-symmetric Horndeski theories that possess a single scalar degree of freedom. Under the assumption of a homogeneous and isotropic cosmological background, the theory admits an inflationary solution that replaces the Big Bang singularity. This can be achieved by a solution where the Hubble expansion rate $H$ is equal to the time derivative of the scalar field $\dot{\phi}$. We show that the solution $H=\dot{\phi}$ suffers from a strong coupling problem, characterized by the vanishing kinetic term of linear scalar perturbations at all times. Consequently, nonlinear scalar perturbations remain uncontrolled from the onset of inflation throughout the subsequent cosmological evolution. Moreover, tensor perturbations are generally subject to Laplacian instabilities during inflation. This instability in the tensor sector also persists under background initial conditions where $H \neq \dot{\phi}$. In the latter case, both the coefficient of the kinetic term for scalar perturbations and the scalar sound speed diverge at the onset of inflation. Thus, the dominance of inhomogeneities in this theory renders the homogeneous background solution illegitimate.

The generalized Lanczos algorithm can provide a universal method for constructing the wave function under the group structure of Hamiltonian. Based on this fact, we obtain an open two-mode squeezed state as the quantum origin for the curvature perturbation. In light of this wave function in the open system, we successfully develop a new method to calculate its corresponding power spectrum by using the Bogoliubov transformation. Unlike traditional approaches, we explicitly retain the Bogoliubov coefficients in terms of the squeezing amplitude \( r_k \) and the squeezing rotation angle \( \phi_k \). As a result, the power spectrum of the open two-mode squeezed state will match that of the Bunch-Davies vacuum numerically. Furthermore, the derivation of the open two-mode squeezed state relies on the second kind Meixner polynomial (equivalent to the generalized Lanczos algorithm) and the symmetry of the Hamiltonian. Therefore, our research may offer a new insight into the calculation of the correlation functions through a group-theoretic perspective.

We study the orbital dynamics and relativistic precession effects in the spacetime of rotating braneworld black holes within the Randall-Sundrum framework. For test particles on spherical orbits, we analyze three conserved quantities-energy, angular momentum, and Carter constant-and examine how the innermost stable spherical orbit depends on the tidal charge and orbital inclination. Compared to Kerr black holes, braneworld corrections significantly modify both nodal and periastron precession frequencies: positive tidal charges suppress precession rates, while negative charges enhance them. For stationary gyroscopes, we calculate the Lense-Thirring precession frequency and demonstrate its sensitivity to the tidal charge, black hole spin, and gyroscope orientation. Our results show that a positive tidal charge weakens frame-dragging effects even as it enhances gravitational attraction-offering a distinctive signature of extra-dimensional gravity. These results have important implications for astrophysical observations, including accretion disk behavior, stellar orbital dynamics, and gravitational wave detection. The modified orbital and gyroscopic precession provide new ways to test braneworld gravity in strong-field regimes.

Background: Constraining the nuclear matter equation of state (EoS) from neutron star observations is one of the main subjects in nuclear physics today. In general, neutron stars rotate rapidly and structure of neutron stars can be affected, especially in millisecond pulsars. To better constrain the nuclear EoS, it is important to describe neutron star structure taking into account the effects of rotation in a fully relativistic manner. Purpose: In this study, we investigate the internal structure of neutron stars under the influence of rotation. We explore correlations between rotational effects and EoS parameters, based on realistic calculations of rapidly rotating neutron stars based on the KEH method, which provides stable solutions for axially symmetric rotating equilibrium configurations. Results: Using 5 different Skyrme EoS parameter sets, we find that the maximum angular frequency achievable by rotating neutron stars, as calculated via the KEH method, varies depending on the stiffness of the equation of state. We confirm that an increase in the rotating frequency leads to an overall increase in both the mass and radius along the M-R curve. By performing calculations at two frequently referenced neutron stars, we further examine how the changes in mass and radius correlate with the nuclear matter properties at saturation density. Our results suggest that the 716Hz rotational constraint may require a more conservative interpretation when accounting for realistic stellar deformation effects. Conclusions: To place stringent constraints on the nuclear EoS based on observational data, it is sometimes essential to account for the effects of rotation in neutron star models. In particular, the influence of rotation becomes increasingly significant at higher spin frequencies and cannot be neglected in rapidly rotating systems with $>$400Hz.

When gravitational waves (GWs) propagate near massive objects, they undergo gravitational lensing that imprints lens model dependent modulations on the waveform. This effect provides a powerful tool for cosmological and astrophysical studies. However, conventional Bayesian parameter inference methods for GWs are computationally expensive, especially for lensed events with additional lens parameters, necessitating more efficient approaches. In this work, we explore the use of neural spline flows (NSFs) for posterior inference of microlensed GWs, and successfully apply NSFs to the inference of 13-dimensional lens parameters. Our results demonstrate that compared with traditional methods like Bilby dynesty that rely on Bayesian inference, the NSF network we built not only achieves inference accuracy comparable to traditional methods for the main parameters, but also can reduce the inference time from approximately 3 days to 0.8 s on average. Additionally, the network exhibits strong generalization for the spin parameters of GW sources. It is anticipated to become a powerful tool for future low-latency searches for lensed GW signals.

Dark matter (DM) remains one of the most compelling unresolved problems in fundamental physics, motivating the search for new detection approaches. We propose a network-based quantum sensor architecture to enhance sensitivity to ultralight DM fields. Each node in the network is a superconducting qubit, interconnected via controlled-Z gates in symmetric topologies such as line, ring, star, and fully connected graphs. We investigate four- and nine-qubit systems, optimizing both state preparation and measurement using a variational quantum metrology framework. This approach minimizes the quantum and classical Cramer-Rao bounds to identify optimal configurations. Bayesian inference is employed to extract the DM-induced phase shift from measurement outcomes. Our results show that optimized network configurations significantly outperform conventional GHZ-based protocols while maintaining shallow circuit depths compatible with noisy intermediate-scale quantum hardware. Sensitivity remains robust under local dephasing noise. These findings highlight the importance of network structure in quantum sensing and point toward scalable strategies for quantum-enhanced DM detection.

In this paper, we investigate the cosmological dynamics of teleparallel dark energy in the presence of nonzero spatial geometry. Extending previous analyses of nonminimal scalar-tensor theories in the torsion-based framework, we consider different scalar field potentials and examine the resulting background evolution and linear perturbations. Adopting a dynamical systems approach, we reformulate the field equations and constrain the model parameters via a Markov chain Monte Carlo analysis combining updated datasets from Pantheon+SH0ES supernovae, cosmic chronometers, and growth rate measurements. Our results suggest a mild preference for an open geometry, although all models remain consistent with a flat universe at the $1\sigma$ level. Notably, Bayesian information criteria indicate that the nonflat teleparallel scenario with a vanishing potential is strongly favored over the standard $\Lambda$CDM model. Furthermore, all teleparallel scenarios are compatible with local determinations of the Hubble constant and exhibit better agreement with low-redshift structure formation data compared to $\Lambda$CDM. These findings highlight the potential of nonflat teleparallel gravity to address current observational tensions and motivate its further investigation as a viable alternative to standard cosmology.

In this Letter, we study scalar wave perturbations of arbitrary frequency to the 5D Schwarzschild-Tangherlini black hole (STBH) within general relativity. For the first time, we derive a closed formula for the 5D partial wave gravitational Raman scattering amplitude applicable to a broad class of boundary conditions, expressed in terms of the Nekrasov-Shatashvili (NS) function for the reduced confluent Heun problem. Furthermore, up to $O(G^2)$ we compute the dynamical $\ell=0$, and the static $\ell=1$, scalar tidal Love numbers of the STBH by matching an effective field theory description for a scalar wave scattering off the black hole, to our novel ultraviolet-NS solutions. The matched Love numbers do not vanish and present renormalization group running behavior.