Papers for Tuesday, Apr 15 2025

Fast and accurate waveform models are fundamentally important to modern gravitational wave astrophysics, enabling the study of merging compact objects like black holes and neutron stars. However, generating high-fidelity gravitational waveforms through numerical relativity simulations is computationally intensive, often requiring days to months of computation time on supercomputers. Surrogate models provide a practical solution to dramatically accelerate waveform evaluations (typically tens of milliseconds per evaluation) while retaining the accuracy of computationally expensive simulations. The GWSurrogate Python package provides easy access to these gravitational wave surrogate models through a user-friendly interface. Currently, the package supports 16 surrogate models, each varying in duration, included physical effects (e.g., nonlinear memory, tidal forces, harmonic modes, eccentricity, mass ratio range, precession effects), and underlying solution methods (e.g., Effective One Body, numerical relativity, black hole perturbation theory). GWSurrogate models follow the waveform model conventions used by the LIGO-Virgo-Kagra collaboration, making the package immediately suitable for both theoretical studies and practical gravitational wave data analysis. By enabling rapid and precise waveform generation, GWSurrogate serves as a production-level tool for diverse applications, including parameter estimation, template bank generation, and tests of general relativity.

The galaxy cluster WHL J013719.8-08284 at $z = 0.566$ exhibits a strong-lensing feature known as the Sunrise Arc, which hosts the candidate star Earendel at $z \approx 6.2$, the most distant star candidate observed to date. If this object is a star, or a system of a few stars, its apparent magnitude implies both extreme gravitational lensing magnification and unusually high luminosity. This study revisits Earendel's magnification, which, in previous literature, exhibits significant uncertainty across various lens models ($2\mu = 4{,}000$-$35{,}000$). We present an improved cluster mass reconstruction and a tighter constraint on Earendel's magnification using a joint strong- and weak-lensing analysis with JWST data. Our strong-lensing mass model, incorporating newly identified multiple-image systems from JWST imaging data and modifying the existing multiple-image assignment scheme, produces a root-mean-square (RMS) lens-plane scatter of less than $0.''3$. Additionally, our weak-lensing catalog achieves a source density of $\sim 100$ galaxies arcmin$^{-2}$, providing constraints on the mass profile beyond the strong-lensing regime. In our best-fit model, we estimate the magnification of Earendel to be $\mu = 43$-$67$, significantly lower than previously proposed and thus calling into question its classification as a star.

Very massive stars (VMS) dominate the light of young stellar populations and are sources of intense stellar feedback. Their evolution is mainly driven by strong wind mass loss, yet current evolution models make simplistic assumptions on their atmospheric physics which are incompatible with the nature of VMS. In this work, we aim to understand VMS atmospheres throughout their evolution by supplementing structure models (computed with GENEC) with detailed atmosphere models (computed with PoWR) capable of capturing the physics of a radially-expanding medium in non-LTE. An important aspect is the computation of atmosphere models reaching into deeper layers of the star, notably including the iron-opacity peak as an important source of radiative driving. In this study, we compute atmosphere models at 16 snapshots along the main sequence of a 150 $M_\odot$ star. For each snapshot, we compute two atmosphere models connected to the underlying structure model at different depths (below and above the hot iron bump). We perform a detailed spectroscopic and structural comparison of the two sequences of model atmospheres, and present a generalized method for the correction of the effective temperature in evolution models with strong winds. The choice of connection point between structure and atmosphere models has a severe influence on the predicted spectral appearance, which constitutes a previously unexplored source of uncertainty in quantitative spectroscopy. The simplified atmosphere treatment of current stellar structure codes likely leads to an overestimation of the spatial extension of very massive stars, caused by opacity-induced sub-surface inflation. This inflation does not occur in our deep atmosphere models, resulting in a discrepancy in predicted effective temperatures of up to 20 kK. Future improvements with turbulence and dynamically-consistent models may resolve these discrepancies.

The ALMACAL-22 survey includes over 2700 hours of observations of ALMA phase and amplitude calibrators, spanning frequencies from 84 to 950 GHz across bands 3 to 10. In total, 687 out of the 1,047 calibrators have redshifts confirmed with spectroscopy and we find an additional 50 featureless blazars. The redshift distribution of the ALMACAL-22 sample peaks at $z \approx 1$ and spans a wide range, from the nuclei of nearby galaxies at $z \ll 0.01$ to quasars at $z = 3.742$. 70 new VLT/X-Shooter spectra of these sources covering UV to NIR wavelengths are also presented, which will be used in future stacking experiments to search for cold gas in the circumgalactic medium. Infrared magnitudes from WISE indicate that the majority of the sources are consistent with being quasars or blazars. After fitting the radio spectral energy distributions of the calibrators, we find that most ALMA calibrators exhibit peaked spectra or are re-triggered which is surprising given the large number of blazars in the sample. The peak frequencies span three orders of magnitude from 100 MHz to 170 GHz, corresponding to linear sizes ranging from sub-pc to $>$ 10 kpc. In the future, when combined with high-resolution radio imaging, these results will offer valuable constraints on the molecular gas content of the CGM, as well as the ages and duty cycles of AGN jets. The ever-growing ALMACAL data set will remain an indispensable resource for studying the various aspects of galaxy formation and evolution.

Dynamical friction (DF) may affect the dynamics of stars moving through dense media. This is the case for stars and compact objects (COs) crossing active galactic nuclei (AGN) discs, stellar clusters, and common envelopes (CE), driving stellar migration. DF may decelerate the moving stellar object and may also, under certain conditions, produce an acceleration. In this paper, we study the DF and its effects in the interaction between a star and the ambient gaseous medium through a set of two-dimensional, hydrodynamical numerical simulations using a wind tunnel configuration. Three different stellar wind configurations are considered: isotropic, polar, and equatorial. We confirm that the DF can decelerate and accelerate the star and find the critical value of the normalized velocity ($u_c$) that marks the transition between these regimes, for the three wind profiles. The value of $u_c$ for the isotropic wind differs slightly from that obtained in the thin shell approximation; for an aspherical wind, it may either be larger or smaller. Aspherical winds with small $u$ values produce larger accelerations than isotropic winds, while at high $u$ values, they lead to greater deceleration than the isotropic case. The timescale for DF to substantially affect the velocity of a stellar object is calculated. It is shown to be relevant in AGN discs and CEs.

Fast X-ray transients (FXTs) are a rare and poorly understood population of events. Previously difficult to detect in real time, the launch of the Einstein Probe with its wide field X-ray telescope has led to a rapid expansion in the sample and allowed the exploration of these transients across the electromagnetic spectrum. EP250108a is a recently detected example linked to an optical counterpart, SN 2025kg, or 'the kangaroo'. Together with a companion paper (Rastinejad et al. 2025), we present our observing campaign and analysis of this event. In this letter, we focus on the early evolution of the optical counterpart over the first six days, including our measurement of the redshift of $z=0.17641$. We find that the source is well-modelled by a rapidly expanding cooling blackbody. We show the observed X-ray and radio properties are consistent with a collapsar-powered jet that is low energy ($\lesssim10^{51}$ erg) and/or fails to break out of the dense material surrounding it. The optical emission therefore likely arises from a shocked cocoon resulting from the trapped jet; however, we also examine the possibility that it emerges from the shock produced as the supernova ejecta expand into a dense shell of circumstellar material. We compare to other supernovae and fast transients showing similar features, finding significant similarities with SN 2006aj and SN 2020bvc. This suggests trapped jets could be more common than previously thought and SN 2025kg may herald a larger sample of similar transients.

With a small sample of fast X-ray transients (FXTs) with multi-wavelength counterparts discovered to date, the progenitors of FXTs and their connections to gamma-ray bursts (GRBs) and supernovae (SNe) remain ambiguous. Here, we present photometric and spectroscopic observations of SN 2025kg, the supernova counterpart to the FXT EP 250108a. At $z=0.17641$, this is the closest known SN discovered following an Einstein Probe (EP) FXT. We show that SN 2025kg's optical spectra reveal the hallmark features of a broad-lined Type Ic SN. Its light curve evolution and expansion velocities are also comparable to those of GRB-SNe, including SN 1998bw, and several past FXT SNe. We present JWST/NIRSpec spectroscopy taken around SN 2025kg's maximum light, and find weak absorption due to He I $\lambda 1.0830, \lambda 2.0581$ $\mu$m and a broad, unidentified feature at $\sim$ 4-4.5 $\mu$m. Further, we observe clear evidence for broadened H$\alpha$ in optical data at 42.5 days that is not detected at other epochs, indicating interaction with hydrogen-rich material. From its light curve, we derive a $^{56}$Ni mass of 0.2 - 0.6 $M_{\odot}$. Together with our companion paper (Eyles-Ferris et al. 2025), our broadband data of EP 250108a/SN 2025kg are consistent with a trapped or low energy ($\lesssim 10^{51}$ ergs) jet-driven explosion from a collapsar with a zero-age main sequence mass of 15-30 $M_{\odot}$. Finally, we show that the sample of EP FXT SNe support past rate estimates that low-luminosity jets seen through FXTs are more common than successful (GRB) jets, and that similar FXT-like signatures are likely present in at least a few percent of the brightest Ic-BL SNe.

Observations indicate that early-type galaxies exhibit varying slopes in the relation between their central stellar surface density and stellar mass ($\Sigma_1 - M_{\star}$). Low-mass galaxies tend to follow a steep slope, close to one, while the slope flattens for high-mass early type galaxies. In our study, we investigate the $\Sigma_1 - M_{\star}$ scaling relation and its evolution using the NIHAO suite of cosmological simulations and compare our findings with recent results from the MaNGA survey. Our analysis shows that NIHAO galaxies successfully reproduce the observed scaling relation based on MANGA survey. Our analysis suggests that AGN feedback plays a critical role in flattening the $\Sigma_1$ slope by expelling gas from galactic centers, leading to a decrease in both stellar and dark matter density as the gravitational potential becomes shallower. To further support our findings, we conducted high-resolution N-body simulations, which confirmed that ({\it sudden}) gas removal does substantially alter the stellar density in the central region, consistent with results from NIHAO. Furthermore, our numerical experiments show that even if the same amount of gas is re-accreted on a typical ({\it longer}) free-fall time, it is not able to restore the original stellar density. Our study concludes that AGN feedback assisted gas removal presents a plausible explanation for the decline in central stellar surface density as observed in massive elliptical galaxies.

Soft gamma repeaters (SGRs) are highly magnetised neutron stars (magnetars) notable for their gamma-ray and X-ray outbursts. In this paper, we use near-infrared (NIR) imaging of SGR 0501+4516 in the days, weeks, and years after its 2008 outburst to characterise the multi-wavelength emission, and to obtain a proper motion from our long temporal baseline observations. Unlike most magnetars, the source has only moderate foreground extinction with minimal crowding. Our observations began only 2 hours after the first activation of SGR 0501+4516 in August 2008, and continued for 4 years, including two epochs of Hubble Space Telescope (HST) imaging. The proper motion constraint is improved by a third HST epoch 10 years later. The near-infrared and X-rays faded slowly during the first week, thereafter following a steeper power-law decay. The behaviour is satisfactorily fit by a broken power-law. Three epochs of HST imaging with a 10-year baseline allow us to determine a quiescent level, and to measure a proper motion of 5.4+/-0.6 mas/yr. This corresponds to a low transverse peculiar velocity of 51+/-14 km/s (at 2 kpc). The magnitude and direction of the proper motion rules out supernova remnant HB9 as the birth-site. We can find no other supernova remnants or groups of massive stars within the region traversed by SGR 0501+4516 during its characteristic lifetime (20 kyr). Our observations of SGR 0501+4516 suggest that some magnetars may be either significantly older than expected, that their progenitors produce low supernova ejecta masses, or alternatively that they can be formed through accretion-induced collapse (AIC) or low-mass neutron star mergers. Although the progenitor of SGR 0501+4516 remains unclear, we propose that SGR 0501+4516 is the best Galactic candidate for a magnetar formed through a mechanism other than massive star core-collapse.

The cosmic star-formation rate density, molecular gas density and the AGN activity of the Universe peak at z~ 2-3, showing the Universe is most active at this epoch. The nature of the galaxies at these redshifts and their properties as a function of their environment are particularly interesting to understand the mechanisms driving their star-formation and quenching. At z~ 2.5, a massive (~ 4.8 X 10^15 Msun) proto-supercluster, Hyperion, was identified Cucciati et al. 2018, consisting of 7 groups/peaks and extending over a comoving volume of 60 X 60 X 150 Mpc^3, providing an excellent laboratory to probe the properties and evolution of galaxies as a function of their environments. We use a large compilation of photometric (optical to radio wavelengths, COSMOS2020, COSMOS-Super-deblended, and, A3COSMOS) and spectroscopic (C3VO, HST-Hyperion, VUDS, zCOSMOS, DEIMOS10K, MAGAZ3NE) data to assign membership and study the relation between the local environment and the molecular gas mass, the star-formation rate (SFR), gas depletion timescales, and quenching mechanisms. We find that the depletion timescales and the molecular gas fractions decrease and SFR increases in denser environments at the ~ 2 sigma level, suggesting accelerated evolution in the densest regions of this proto-supercluster resulting from gas stripping, over-consumption, and/or cessation of cold flows. Dedicated observations at sub-millimeter wavelengths enabling further spectroscopic confirmation and better coverage in the sub-millimetric (sub-mm) wavelengths can provide more conclusive results on the environmental implications on gas reservoirs of galaxies in Hyperion.

We report the combined JWST NIRSpec/G395H and NIRISS/SOSS transmission spectrum of the transiting super-Jupiter HAT-P-14 b, from 0.60 $\mu m$ to 5.14 $\mu m$. Initial analysis of these data reported a near-featureless spectrum at NIRSpec wavelengths range (2.87 $\mu m$ to 5.14 $\mu m$) consistent with the small atmospheric scale height of the planet and unexplained bumps and wiggles at NIRISS wavelengths range (0.6 $\mu m$ to 2.8 $\mu m$). Here, we produce a self-consistent spectrum of HAT-P-14 b's atmosphere with an up-to-date reduction. We detect H$_2$O (3.09 $\sigma$) both across NIRISS/SOSS wavelengths range and at the bluest end of NIRSpc/G395H as well as a gray cloud deck (1.90 $\sigma$). We constrain the atmospheric metallicity of HAT-P-14 b to be roughly Solar, with [Fe/H] $= -0.08^{+0.89}_{-0.98}$, consistent with the planet mass-metallicity relationship. The differences compared to previous works are likely due to the improved STScI jwst pipeline, which highlights the need to reanalyze the early NIRISS/SOSS transiting exoplanet targets with the latest methods. As HAT-P-14 b is placed as the 805th best target for transmission spectroscopy according to Transmission Spectroscopy Metrics (TSM), our results showcase JWST's unparalleled photometric precision which can easily characterize a thousand exoplanets' atmospheres through transmission spectroscopy.

Distance ladders which calibrate the luminosity of Type Ia supernovae (SNe Ia) currently provide the strongest constraints on the local value of H0. Recent studies from the Hubble Space Telescope (HST) and James Webb Space Telescope (JWST) show good consistency between measurements of SNe Ia host distances. These are calibrated to NGC 4258 using different primary distance indicators (Cepheids, Tip of the Red Giant Branch (TRGB), J-region Asymptotic Giant Branch, and Miras). However, some sub-samples of calibrated SNe Ia employed to measure H0 yield noteworthy differences due to small sample statistics but also due to differences in sample selection. This issue is particularly important for TRGB-derived calibrations owing to the smaller volume they reach compared to Cepheids, reducing sample size and enhancing the size of statistical fluctuations. To mitigate this issue, we compile the largest and complete (as currently available) sample of HST or JWST measurements of the TRGB in the hosts of normal SNe Ia for a total of N=35, 50% larger than the previous largest. Most are present in the literature, and we compile multiple measures when available. We also add 5 SNe Ia hosts from the HST archive not previously published. The full sample together with the Pantheon+ SN catalog gives H0=72.1-73.3 +/- 1.8 km/s/Mpc (depending on methodology), in good agreement with the value of 72.5 +/- 1.5 km/s/Mpc from HST Cepheids in hosts of 42 SNe Ia calibrated by the same anchor, NGC 4258. We trace the difference in the result of H0=70.4 +/- 1.9 km/s/Mpc from Freedman et al. 2025 to 11 hosts not selected for that CCHP compilation (of N=24) which alone yield H0=74.1 km/s/Mpc, 2$\sigma$ higher than the selected sample. A smaller increase of 0.6 km/s/Mpc comes from a commonly employed correction for peculiar velocities.

We investigate the response of warm absorbers to variations in the ionizing continuum of the Seyfert 1 galaxy NGC 4051 using time-resolved X-ray observations from the \textit{Neutron Star Interior Composition Explorer} (\textit{NICER}). In this work, we have demonstrated we can perform time-resolved spectroscopic studies of warm absorbers of about $\sim 5500$ s time resolution using NICER data. We have extracted 15 spectra for this source, corresponding to 15 different visits to the source, or pointings, each separated by a longer Earth occultation. By modeling the spectral variability of the warm absorber with the \texttt{\sc {\sc warmabs}} analytic model, we detect significant variations in the ionization parameter that correlate with changes in the ionizing flux. A time lag of approximately 5500 seconds is observed between the flux variations and the absorber's ionization response, suggesting that the gas is out of photoionization equilibrium during these periods. Using this time lag, we estimate the lower limit of the gas density $8.91 \times 10^6 \, \text{cm}^{-3}$ and constrain the location of the warm absorber to within $7.02 \times 10^{16} \, \text{cm}$ ($\sim 0.02$ parsec) from the central black hole. This study uses time-resolved spectral analysis to contribute to our understanding of the physical conditions of ionized AGN outflows, such as density and location.

(Abridged) Neutron stars (NSs), the densest known objects composed of matter, provide a unique laboratory to probe whether strange quark matter is the true ground state of matter. We investigate the parameter space of the equation of state of strange stars using a quantum chromodynamics (QCD)-informed model. The parameters - related to the energy density difference between quark matter and the QCD vacuum, the strength of strong interactions, and the gap parameter for color superconductivity - are sampled via quasi-random Latin hypercube sampling to ensure uniform coverage. To constrain them, we incorporate observational data on the maximum mass of NSs (from binary and merger systems), the radii of $1.4$ M$_{\odot}$ NSs (from gravitational wave and electromagnetic observations), and tidal deformabilities (from GW170817). Our results show that quark strong interactions play a key role, requiring at least a $20\%$ deviation from the free-quark limit. We also find that color superconductivity is relevant, with the gap parameter reaching up to $\sim 84$ MeV for a strange quark mass of $100$ MeV. The surface-to-vacuum energy density jump lies in the range $(1.1-1.3)$ $\rho_{\rm{sat}}$, where $\rho_{\rm{sat}} \simeq 2.7 \times 10^{14}$ g cm$^{-3}$. Observational constraints also imply that a $1.4$ M$_{\odot}$ quark star has a radius of $(11.5-12.3)$ km and tidal deformability between $670$ and $970$. These are consistent with the low mass and radius inferred for the compact object XMMU J173203.3-344518. Our results provide useful inputs for future studies on quark and hybrid stars, including their tidal properties, thermal evolution, quasi-normal modes, and ellipticities.

Viewing high-redshift sources at near-opposite directions on the sky can assure, by light-travel-time arguments, acausality between their emitted photons. One utility would be true random-number generation, by sensing these via two independent telescopes that each flip a switch based on those latest-arrived colours; for example, to autonomously control a quantum-mechanical (QM) experiment. Although demonstrated with distant quasars, those were not fully acausal pairs, which are restricted in simultaneous view from the ground at any single observatory. In optical light such faint sources also require large telescope aperture to avoid sampling assumptions when imaged at fast camera framerates: either unsensed intrinsic correlations between them or equivalently-correlated noise may ruin the expectation of pure randomness. One such case which could spoil a QM test is considered. Based on that, allowed geometries and instrumental limits are modelled for any two ground-based sites, and their data simulated. To compare, an analysis of photometry from the Gemini twin 8-m telescopes is presented, using archival data of well-separated bright stars, obtained with the instruments 'Alopeke (on Gemini-North in Hawai'i) and Zorro (on Gemini-South in Chile) simultaneously in two bands (centred at 562 nm and 832 nm) with 17 Hz framerate. No flux correlation is found, calibrating an analytic model, predicting where a search at signal-to-noise over 50 at 50 Hz with the same instrumentation can be made. Finally, the software PDQ (Predict Different Quasars) is presented which searches a large catalogue of known quasars, reporting those with brightness and visibility suitable to verify acausal, uncorrelated photons at those limits.

Metals -- heavy elements synthesized during various phases of stellar evolution or during supernova explosions -- play a fundamental role in shaping galaxy evolution. In fact, their relative abundances, spatial distribution, and scaling with galactic properties reflect the constant interplay between star-formation, nucleosynthesis, and gas flows that drive the cycle of baryons in-and-out of galaxies across the cosmic time. This chapter aims at offering a concise introduction to the methodologies used to measure elemental abundances in galaxies and the basics of chemical evolution modeling. We also provide a brief overview of the current observational framework, including metallicity scaling relations, the study of relative chemical abundances, the distribution of metals within and beyond galaxies, and how these properties evolve with redshift, drawing on both extensive literature and recent developments, and aiming to highlight well-established findings alongside ongoing challenges in this rapidly advancing field.

We use numerical relativity to study the violent preheating era at the end of inflation. This epoch can result in highly nonlinear fluctuations in density and gravitational potential which feed back onto the averaged expansion rate -- an effect known as backreaction. Usually, simulations of preheating use the Friedmann constraint to enforce the Hubble expansion of spacetime during the evolution. In numerical relativity, this is not required and the inhomogeneous spacetime is evolved self-consistently. For a 'vanilla' preheating model, we find a violation of the Friedmann constraint at the level of $0.005\%$ over the entire simulation. This violation increases to $\sim10\%$ as we sample smaller scales in the simulation domain.

The cosmic "Dark Ages" is the period between the last scattering of the Cosmic Microwave Background (CMB) and the appearance of the first luminous sources, spanning redshifts $1100\gtrsim z\gtrsim 30$. The only way to observe this period is by examining the 21 cm hyperfine transition line of neutral hydrogen HI, which -- given the high redshifts (and hence long wavelengths) -- must be observed from outside the Earth's ionosphere. Given the faintness of the signal, concepts for a radio array on the lunar far side (where large collecting areas can be deployed and radio frequency interference is minimal) have been proposed, like FarView or FARSIDE, but designs are still in the preliminary stages. This paper studies multiple aspects of array design to determine the impact of different design decisions on sensitivity to the Dark Ages 21 cm power spectrum. We do so by using the sensitivity package 21cmSense to model and simulate various array configurations. We present a fiducial design based on a modification of the FarView concept, which consists of a collecting area of $\sim2.5\,\rm{km}^{2}$ with 82,944 tightly packed dual-polarization dipoles grouped into 5,184 correlated elements, or subarrays, delivering a $>10\sigma$ detection of the $z=30$ signal with a five year lifetime. We find that, beyond mere collecting area, the most important factor in achieving this sensitivity is the presence of very short baselines that can only be realized with small, closely packed antennas.

Detecting ultrahigh-energy neutrinos can take two complementary approaches with different trade-offs. 1)~Wide and shallow: aim for the largest effective volume, and to be cost-effective, go for wide field-of-view but at the cost of a shallow instantaneous sensitivity -- this is less complex conceptually, and has strong discovery potential for serendipitous events. However, it is unclear if any source can be identified, following detection. And 2)~Deep and narrow: here one uses astrophysical and multi-messenger information to target the most likely sources and populations that could emit neutrinos -- these instruments have deep instantaneous sensitivity albeit a narrow field of view. Such an astrophysically-motivated approach provides higher chances for detection of known/observed source classes, and ensures multi-messenger astronomy. However, it has less potential for serendipitous discoveries. In light of the recent progress in multi-messenger and time-domain astronomy, we assess the power of the deep and narrow instruments, and contrast the strengths and complementarities of the two detection strategies. We update the science goals and associated instrumental performances that envisioned projects can include in their design in order to optimize discovery potential.

Recent spectroscopic observations of the [C\,{\tiny II}] 158$\,\mathrm{\mu m}$ fine-structure line of ionized carbon (C$^+$), using the Stratospheric Observatory for Infrared Astronomy (SOFIA), have revealed expanding [C\,{\tiny II}] shells in Galactic H\,{\tiny II} regions. We report the discovery of a bubble-shaped source (S144 in RCW79), associated with a compact H\,{\tiny II} region, excited by a single O7.5--9.5V/III star, which is consistent with a scenario that the bubble is still mostly ``filled'' with C$^+$. This indicates most likely a very early evolutionary state, in which the stellar wind has not yet blown material away, as it is the case for more evolved H\,{\tiny II} regions. Using the SimLine non-LTE radiative transfer code, the [C\,{\tiny II}] emission can be modeled to originate from three regions. First, a central H\,{\tiny II} region with little C$^+$ in the fully ionized phase, followed by two layers with gas density around $2500\,\mathrm{cm^{-3}}$ of partially photo-dissociated gas. The second layer is a slowly expanding [C\,{\tiny II}] shell with an expansion velocity of $\sim\,$$2.6\,\mathrm{km\,s^{-1}}$. The outermost layer exhibits a temperature and velocity gradient that produces the observed self-absorption features in the optically thick [C\,{\tiny II}] line ($\tau \sim 4$) leading to an apparent deficit in [C\,{\tiny II}] emission and a low ratio of [C\,{\tiny II}] to total far-infrared (FIR) emission. We developed a procedure to approximate the missing [C\,{\tiny II}] flux and find a linear correlation between [C\,{\tiny II}] and FIR without a [C\,{\tiny II}]-deficit. This demonstrates that at least some of the [C\,{\tiny II}]-deficit found in Galactic H\,{\tiny II} bubbles can be attributed to self-absorption.

We present an X-ray spectro-polarimetric study of the weakly magnetized neutron star low-mass X-ray binary GX 9+1, utilizing data from the Imaging X-ray Polarimetry Explorer (IXPE), alongside simultaneous NuSTAR, NICER, and INTEGRAL observations. GX 9+1, located in the Galactic bulge, is a persistently bright Atoll source known for its spectral variability along the color-color diagram. Our spectral analysis during the soft state confirms emission dominated by a soft blackbody and thermal Comptonization components, with no evidence of a hard X-ray tail. These observations suggest a relatively low-inclination system (23 deg < i < 46 deg) with a weak reflection component, consistent with emission from the accretion disk and neutron star boundary layer. Spectro-polarimetric analysis reveals no significant polarization in the 2-8 keV range, with a 3-sigma upper limit for the polarization degree of 1.9%. However, marginal evidence of polarization was detected in the 2-3 keV band at the 95.5% confidence level (2-sigma), suggesting potential contributions from scattering effects in the individual spectral components (disk, reflection, and Comptonization) that could cancel each other out due to the different orientations of their polarization angles. This behavior aligns with other Atoll sources observed by IXPE, which typically exhibit lower and less variable polarization degrees compared to Z-class sources.

We present observations of an eruptive solar flare on 2016 January 6 that occurred behind the solar limb from the perspective of the Earth, but was well observed by STEREO and the Solar Extreme Ultraviolet Monitor on the Mars Atmosphere and {Volatile} EvolutioN (MAVEN) mission. Light curves showing the evolution of the flare's irradiance as a function of time taken by MAVEN are well correlated with the brightness evolution of fan structures observed in the PROBA2 SWAP 174 Å passband, suggesting that the radiance of structures near the flare site was influenced by emission from the flare. Because SWAP did not directly observe the flare itself, this event represents a rare opportunity to study the evolution of emission near a flare without the risk of instrumental scattered light contaminating the observations. We analyze this evolution and implement a simple model to explore the possibility that resonant excitation (or resonant scattering) plays an important role in driving coronal EUV emission during flaring events. Our modeling shows that for a large flare, resonant excitation could increase emission from nearby structures by about 45%, consistent with our findings that the involved structures observed by SWAP increased in brightness by about 60% during the flare. We conclude that resonant excitation may play an important role in driving coronal EUV emission under certain circumstances and should be accounted for in models and emission-based analysis tools.

A white paper submitted to the 2025 NASA Decadal Astrobiology Research and Exploration Strategy (DARES) on the importance of early-career training, support, and retention. The paper identifies two goals for NASA Astrobiology regarding early career researchers (ECRs): (1) Knowledge Retention and Workforce Stability, and (2) Foster Collaboration & Strengthen Community. The paper outlines the challenges of achieving these goals and offers recommendations for actions that NASA Astrobiology can take to further train, support, and retain ECRs in NASA Astrobiology.

In this work we present the results of our analysis of medium-resolution LAMOST spectra of late-type candidate members of the Pleiades. We have used the code ROTFIT to determine the atmospheric parameters (Teff, logg, and [Fe/H]), radial velocity (RV), and projected rotation velocity (vsini) for 1581 spectra of 283 stars. Moreover, for late-type stars (Teff<6500 K), we also calculated the H$\alpha$ and LiI-6708 net equivalent width by means of the subtraction of inactive photospheric templates. We have also used rotation periods available in the literature and we have determined them for 89 stars by analyzing TESS photometry. The RV distribution of the members peaks at 5.0 km/s with a dispersion of 1.4 km/s, while the average metallicity is [Fe/H]=-0.03$\pm$0.06, in line with previous determinations. Fitting empirical Lithium isochrones, we obtain a reliable age for the Pleiades of 118$\pm$6 Myr, in agreement with the recent literature. The activity indicators H$\alpha$ line flux (Fha) and luminosity ratio (R'ha) show the hottest stars to be the less active ones. The R'ha values display the typical activity-rotation trend with the Rossby number with a steep decay for $R_{\rm O}\geq$0.2 and a nearly flat (saturated) activity level for smaller values. However, we still see a slight dependence on $R_{\rm O}$ in the saturated regime which is well fitted by a power law with a slope of $-0.18\pm0.02$, in agreement with some previous work. For three sources we have found LAMOST spectra acquired during flares showing strong and broad H$\alpha$ profiles and the presence of the HeI-6678 emission line. We identify 39 possible SB1 and ten SB2 systems. We have also shown the potential of the LAMOST-MRS spectra, which allowed us to refine the orbital solution of some binary and to discover a new double-lined binary.

The presence of a shallow heat source of unknown origin in accreting neutron star crusts has been inferred by analyzing their cooling behavior in quiescence. To investigate a diverse bursting history for KS 1731-260 during accretion outbursts, we use realistic crust compositions and nuclear heating and cooling sources from detailed nuclear reaction network calculations to interpret observed cooling curves. We find that the required strength of the shallow heat source is reduced by more than a factor of 3 compared to previous analysis, and obtain constraints on the most likely dominant surface burning modes of KS 1731-260 over its history. Our analysis suggests an impure nuclear pasta layer in the inner crust, though future observations will provide more stringent constraints.

In the framework of general relativity, dark energy was proposed to explain the cosmic acceleration. A pivotal inquiry in cosmology is to determine whether dark energy is the cosmological constant, and if not, the challenge lies in constraining how it evolves with time. In this paper, we utilize the latest observational data to constrain some typical dark energy models, and make a comparison for them according to their capabilities of fitting the current data. Our study is confined to late-universe observations, including the baryon acoustic oscillation, type Ia supernova, cosmic chronometer, and strong gravitational lensing time delay data. We employ the Akaike information criterion (AIC) and Bayesian information criterion (BIC) to assess the worth of models. The AIC analysis indicates that all dark energy models outperform the $\Lambda$CDM model. However, the BIC analysis leaves room for $\Lambda$CDM due to its heavier penalty on the model complexity. Compared to $\Lambda$CDM, most dark energy models are robustly supported by AIC while being explicitly disfavored by BIC. The models that are robustly favored by AIC and not explicitly disfavored by BIC include the $w$CDM, interacting dark energy, and Ricci dark energy models. Furthermore, we observe that an alternative modified gravity model exhibits superior performance when compared with $\Lambda$CDM from both the AIC and BIC perspectives.

We present the long-term spectral evolution of eight black hole ultra-luminous X-ray sources (BH-ULXs), namely NGC1313 X-1, NGC5408 X-1, NGC6946 X-1, IC342 X-1, NGC55 ULX1, NGC4395 ULX1, NGC5204 X-1 and NGC4190 ULX1 using {\it XMM-Newton} monitoring data spanning over a decade or more. An in-depth spectral modeling with thermal Comptonization ({\it nthComp}) and standard disc ({\it diskbb}) components reveals NGC5204 X-1, IC342 X-1, NGC4190 ULX1 and NGC1313 X-1 exhibiting harder spectral characteristics with dominant effect of Comptonization ($F_{nth}>F_{disc}$, $\Gamma_{nth}\lesssim2$). However, NGC6946 X-1 and NGC55 ULX1 remain in a disc-dominated state ($F_{disc}\sim2F_{nth}$, $\Gamma_{nth}\gtrsim2$), while NGC5408 X-1 shows intermediate spectral characteristics. The spectral analyses indicate an anti-correlation between disc luminosity ($L_{disc}$) and temperature ($T_{col}$) for all sources except NGC5204 X-1. These anti-correlations follow a relation $L_{disc}\propto T_{col}^{\alpha}$ with steeper exponents of $\alpha=-6.01\pm0.25$ (NGC55 ULX1), $-8.93\pm0.11$ (NGC6946 X-1), and $-10.31\pm0.10$ (NGC5408 X-1) for sources with softer or intermediate spectral characteristics. For harder sources, NGC1313 X-1 and IC342 X-1, the combined results provide $\alpha=-3.58\pm0.04$. However, for NGC5204 X-1, a positive correlation is found, yielding $\alpha=1.4\pm0.1$, suggesting that the emission mechanism is associated with the transition from the `standard disc' to the `slim disc' scenario. These findings suggest that the observed $L_{disc}-T_{col}$ correlations, along with the overall spectro-temporal properties of BH-ULXs, seems to be governed by disc-corona-wind driven accretion processes at various inclinations. Finally, we report a QPO-like feature ($\sim20$ mHz) with $rms\%\sim6.6$, Q-factor $\sim6.7$ and significant $2.8\sigma$ in NGC55 ULX1.

From a long-term Doppler monitoring campaign of 373 giant stars, we have identified ten giants with periodic radial velocity variations that are challenging to associate with planets. Similar cases in the literature are attributed to poorly understood intrinsic processes. Our goal is to confirm or refute the presence of planets around these ten evolved stars. Additionally, we evaluate the reliability and sensitivity of planet-confirmation metrics when applied to giant stars and present cases of intrinsically induced radial velocity variations, aiming to enhance the physical understanding of the phenomenon. We combined 25 years of radial velocity data from the Hamilton/Lick, SONG, and CARMENES spectrographs. To assess consistency with Keplerian models, we examined the residuals and tracked changes in statistical significance as new data were incorporated. Additionally, we compared radial velocity amplitudes across optical and infrared wavelengths, searched for periodic variations of activity indicators, and examined their correlations with radial velocities. Seven of the ten giants exhibit intrinsically induced radial velocity variations. The strongest arguments against planets orbiting the giants are guided by long-term radial velocity monitoring that detects changing periodicity on long timescales or detects systematics close to the original period in the radial velocity residuals. While activity indicators offer some support, their signals are generally weak. Comparing optical and infrared radial velocity amplitudes also proves insufficient for confirming or refuting planets. We find HIP 64823 remains a promising candidate for hosting a giant exoplanet with orbital period $P\sim$ 7.75 yr. For two stars, the evidence remains inconclusive. Long-term radial velocity monitoring is essential for distinguishing planetary companions from intrinsic variations in evolved stars.

Radio Frequency Interference (RFI) presents a significant challenge for carrying out precision measurements in radio astronomy. In particular, RFI can be a show stopper when looking for faint cosmological signals such as the red-shifted 21-cm line from cosmic dawn (CD) and epoch of reionization (EoR). As wireless communications, satellite transmissions, and other RF technologies proliferate globally, understanding the RFI landscape has become essential for site selection and data integrity. We present findings from RFI surveys conducted at four distinct locations: three locations in India, the Gauribidanur Radio Observatory in Karnataka, Twin Lakes in Ladakh, Kalpong Dam in the Andaman Islands, and the Gruvebadet Atmosphere Laboratory in Ny-Alesund, Svalbard, Norway. These sites, selected based on their geographical diversity and varying levels of human activity, were studied to assess RFI presence in 30-300 MHz bands-critical for low-frequency observations and experiments targeting the 21-cm CD/EoR signal. Using an automated RFI detection approach via the Hampel filter and singular value decomposition, the surveys identified both persistent and transient interference, which varies with location and time. The results provide a comprehensive view of the RFI environment at each site, informing the feasibility of long-term cosmological observations and aiding in the mitigation of RFI in radio astronomical data. The methods developed to characterize RFI can be easily generalized to any location and experiment.

We study the cosmological constraints of the time variation of the electron mass $m_e$ and the fine-structure constant $\alpha$, using data of cosmic microwave background, supernovae light curve and baryon acoustic oscillation (BAO) data including the recent DESI BAO DR2 measurements. The results are slightly depending on the BAO data set included in the analysis. The latest DESI BAO DR2 data strongly indicates that $m_e$ or $\alpha$ is slightly larger than the previous data from 6DF+SDSS and DESI BAO DR1. We also compare the varying $m_e$ model, the varying $\alpha$ model, and the simultaneous variation of $m_e$ and $\alpha$. When considering the Hubble tension, a larger electron mass is the most promising option and the variation of the fine-structure constants does not alleviate the tension.

A statistical analysis of magnetic energies of the nonlinear force-free and potential fields, and their difference (a proxy for the free magnetic energy) in active regions (ARs) on the Sun of different Hale (Mount Wilson) and McIntosh classes for the period from May 1, 2010 to June 12, 2024 is presented. The magnetic fields in ARs are calculated using the GX Simulator based on the information about ARs contained in the daily Solar Region Summary (SRS) files provided by the NOAA SWPC and vector magnetograms by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). Total unsigned and signed magnetic fluxes and vertical electric currents on the photosphere are also calculated. For the parameters considered, distributions have been determined in total for all ARs and separately for each Hale and McIntosh class. Minimum, maximum, mean values of the parameters and standard deviations were calculated for each class. The information about the parameters is presented in the form of graphs and tables. The magnetic energies, unsigned magnetic flux, unsigned vertical current, as well as the integral number of sunspots, number of ARs, and area of sunspots, integrated over ARs visible per day on the solar disk, exhibit similar approximately 11.6-year cyclicity. On average, magnetic energies of ARs increase with increasing Hale and McIntosh class, while the average fraction of the free magnetic energy in ARs of different classes differs weakly. We also found that the Poisson Flare Probabilities (PFPs) correlate with the parameters, and the Pearson correlation coefficient is up to 0.89. The results reveal relationships between various parameters of ARs and may be used in developing prediction of space weather effects.

In protoplanetary disks, the formation of planetesimals via streaming and/or gravitational instabilities requires regions with a locally enhanced dust-to-gas mass ratio. Conventionally, gas pressure maxima sustained by gas surface density maxima have been considered as the primary cause of such dust accumulation. However, the disk's pressure structure depends not only on gas density but also on the temperature structure, which itself is influenced by the distribution of dust. In this study, we propose a novel mechanism for dust accumulation, which is driven by the coevolution of dust and disk temperature. In the inner disk region where the midplane temperature is primarily determined by the balance between viscous heating and radiative cooling, a perturbation in dust surface density distribution may affect radiative cooling efficiency, potentially producing a local maximum in the temperature and pressure profiles. To test this hypothesis, we perform coupled calculations of dust and disk temperature evolution, incorporating the advection, diffusion, coagulation, and fragmentation of dust particles along with viscous heating, radiative cooling, and radial thermal diffusion. Our results demonstrate that a pressure maximum formed by a perturbation in the dust surface density can spontaneously induce dust accumulation, even in the absence of a gas surface density maximum, under conditions where dust drift is significantly faster than diffusion and the thermal evolution occurs faster than the inward migration of dust. This mechanism requires viscous heating to dominate disk heating, and typically occurs interior to the snow line. In this spontaneous dust trap, the dust-to-gas density ratio at the midplane can exceed unity, suggesting the potential for rocky planetesimal formation via streaming and gravitational instabilities.

We present a practical implementation of the perturbation theory derived by Lynden-Bell (2015) for describing, to arbitrary precision, the orbit of a particle in an arbitrary spherically-symmetric potential. Our implementation corrects minor but important errors in the initial derivation, and extends the formalism in two ways. First, a numerical method is developed to efficiently and precisely solve the analogue to the Kepler problem, and second, a method is introduced to track the particle's vertical oscillations about an axisymmetric disk, even when the vertical oscillation frequency varies with radius. While not as flexible as numerical integration, this method guarantees conservation of energy, angular momentum, and related quantities, and may be used to evaluate a particle's position and velocity in constant time. Our implementation is written in Python and is pip installable as the package lbparticles.

Dust-obscured galaxies (DOGs) are considered to be in a co-evolution phase, with the associated active galactic nuclei (AGN) obscured by dust and gas. Although the DOGs are thought to harbor rapidly growing SMBHs, their X-ray statistical properties, crucial for understanding the properties of obscuring gas as well as the accretion disk state and the hot electron corona around the SMBHs, remain unexplored due to the combination of the low number density of DOGs and the lack of X-ray surveys achieving both of the wide-area and uniformly high-sensitivity observations. We construct a sample of X-ray-detected DOGs in the eROSITA Final Equatorial Depth Survey (eFEDS) field and examine their X-ray statistical properties. By using Subaru/HSC SSP, VIKING, and WISE all-sky surveys, our results reveal the discovery of 5738 IR-bright DOGs in the footprint covered by both of the eFEDS and VIKING surveys (60 deg^2), with 65 objects identified as X-ray-detected DOGs (eFEDS-DOGs). Among them, 41 eFEDS-DOGs show a power-law slope in the near to mid-IR bands (power-law DOGs), indicating dust-obscured AGN. The hydrogen column density (N_H) suggests that eFEDS-DOGs cover even unobscured AGN, spanning 10^20 < N_H <= 10^23. On the other hand, the majority of IR-bright DOGs are not detected by eROSITA, suggesting that most IR-bright DOGs are heavily obscured by dust and gas with N_H > 10^23. Therefore, eFEDS-DOGs, discovered thanks to the wide-area survey by eROSITA, are newly found populations showing less obscured phases among the lifetime of DOGs. Additionally, some eFEDS-DOGs exhibit deviations, down to nearly 1.0 dex below the monochromatic luminosity at 6 micron versus absorption-corrected intrinsic X-ray luminosity between 0.5-2 keV relation, suggesting that it may signal high Eddington ratios reaching the Eddington limit.

The origin of lambda Boo stars is currently unknown. After several efforts by many authors, no bona fide lambda Boo stars have been confirmed as members of open clusters. Their detection could provide an important test bed for a detailed study of lambda Boo stars. Results. For the first time, we present the surprising finding of two lambda Boo stars as members of open clusters: HD 28548 belongs to the cluster HSC 1640 and HD 36726 belongs to the cluster Theia 139. This was confirmed using a detailed abundance analysis, while the cluster membership was independently analyzed using Gaia DR3 data and radial velocities. We compared the lambda Boo star HD 36726 with other cluster members and showed that the lambda Boo star was originally born with a near-solar composition. This also implies one of the highest chemical differences detected between two cluster members (0.5 dex). In addition, we suggest that the lambda Boo peculiarity strongly depletes heavier metals, but could also slightly modify lighter abundances such as C and O. We also found that both lambda Boo stars belong to the periphery of their respective clusters. This would suggest that lambda Boo stars avoid the strong photoevaporation by UV radiation from massive stars in the central regions of the cluster. We preliminarily suggest that peripheral location appears to be a necessary, though not sufficient, condition for the development of lambda Boo peculiarity. We also obtained a precise age determination for the lambda Boo stars HD 28548 (26.3 Myr) and HD 36726 (33.1 Myr). Conclusions. We have confirmed, for the first time, that two lambda Boo stars belong to open clusters. This remarkable finding could make open clusters excellent laboratories for studying the origin of lambda Boo stars.

Stellar spectral interpolation is critical technique employed by fitting software to derive the physical parameters of stars. This approach is necessary because on-the-go generation of synthetic stellar spectra is not possible due to the complex and high cost of computation. The goal of this study is to develop a spectral interpolator for a synthetic spectral library using artificial neural networks (ANNs). The study aims to test the accuracy of the trained interpolator through self-inversion and, subsequently, to utilize the interpolator to derive the physical parameters of stars in the globular cluster NGC 6397 using spectra obtained from the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT). In this study, ANNs were trained to function as spectral interpolators. The ULySS full-spectrum fitting package, integrated with the trained interpolators, was then used to extract the physical parameters of 1587 spectra of 1063 stars in NGC 6397. The trained ANN interpolator achieved precise determination of stellar parameters with a mean difference of 31 K for $T_{\rm eff}$ and 0.01 dex for [Fe/H] compared to previous studies. This study demonstrates the efficacy of ANN-based spectral interpolation in stellar parameter determination, offering faster and more accurate analysis.

The acceleration process of massive particles as well as the production of very high energy (VHE) photons and neutrinos remains a fundamental challenge in astrophysics. We investigate the parsec-scale properties of the blazar PKS 1424+240, that has been selected on the basis of both gamma-ray and neutrino VHE radiation. We analyze VLBA observations of this BL Lac object, stacking 42 polarization-sensitive images collected in 2009-2025 to enhance the signal and reveal persistent parsec-scale structure. Our observations indicate that this object is viewed inside the jet cone, very close to the axis of its relativistic jet, with a viewing angle of <0.6 deg. This effectively maximizes Doppler boosting to values ~30 and enhances both electromagnetic and neutrino emission in the observer's direction. Based on polarimetric observations, we unambiguously detect a net toroidal component of the jet's magnetic field. Blazars with very small jet viewing angles offer a solution to the longstanding mismatch between Doppler factors inferred from VLBI and those derived from VHE observations -- the so-called `Doppler factor crisis'. We show that relativistic beaming plays the critical role in the gamma-ray and neutrino emission of blazars, with direct implications for models of their multi-messenger emission.

We have constructed and analysed the secular light curve of BP Tau, a classical T Tauri-type star. Wave-like variations in the average brightness were detected, with an amplitude of $\Delta B\approx 0.2$ and characteristic time-scales of several decades. We argue that three deep dimming events $(\Delta B \sim 1.5)$, lasting from 1 hour to several days, are caused by the eclipse of a hot (accretion) spot by dust falling onto the star together with gas. Such eclipses, albeit with smaller amplitudes, may explain the absence of a strictly defined periodicity in the brightness variations of BP Tau associated with axial rotation. We also show that within the distance range of 0.1 to 200 AU, BP Tau does not have a companion with a mass exceeding 0.2M$_\odot.$ The causes of brightness and colour index variations on different time-scales are discussed.

The wide-area component of the LOFAR Two-Metre Sky Survey (LoTSS) is currently the largest radio survey ever carried out, and a large fraction of the 4.5 million radio sources it contains have been optically identified with galaxies or quasars with spectroscopic or photometric redshifts. Identification of radio-luminous AGN from this LoTSS source catalogue is not only important from the point of view of understanding the accretion history of the universe, but also enables a wide range of other science. However, at present the vast majority of the optical identifications lack spectroscopic information or well-sampled spectral energy distributions. We show that colour and absolute magnitude information from the Wide-Field Infrared Survey Explorer (WISE) allows for the robust and efficient selection of radio AGN candidates, generating a radio AGN candidate sample of around 600,000 objects with flux density $> 1.1$ mJy, spanning 144-MHz luminosities between $10^{21}$ and $10^{29}$ W Hz$^{-1}$. We use the catalogue to constrain the total sky density of radio-luminous AGN and the evolution of their luminosity function between $z=0$ and $z\approx 1$, and show that the typical mass of their host galaxies, around $10^{11} M_\odot$, is essentially independent of radio luminosity above around $L_{144} \approx 10^{24}$ W Hz$^{-1}$. Combining with Very Large Array Sky Survey (VLASS) data, we show that the core prominences, radio spectral indices and variability of extended sources from the sample are qualitatively consistent with the expectations from unified models. A catalogue of the radio AGN candidates is released with this paper.

We present a new analytical approach to the longitudinal development of electromagnetic air showers, offering improvements to the classical Greisen formalism. We introduce a novel profile for the slope function $\lambda_1(s)$ that achieves an agreement less than $0.75\%$ with the original $\lambda_1$ for shower age parameter $s$ between $0.3 < s < 1.4$, where $s$ represents the stage of shower development. Our new formalism provides an improved representation of shower evolution, particularly near and beyond the shower maximum. In addition, we derive a complete expression for the number of particles $N(t)$. Our implementation includes the zenithal angle dependence on the number of particles at the detector level at high altitudes, making it particularly useful for high-altitude observatories. This expression is suitable for implementing air shower simulation tool fitting procedures over a wide range of energies and geometries. Our analysis suggests that the proposed new formalism may provide better agreement with the expected evolution of particle numbers compared to the traditional Greisen formulation.

Kepler and TESS observations led to the discovery of many close-in super Earths, including some with ultra-short orbital periods ($\lesssim 1$ day). During and shortly after their multi-Myr formation epoch, their GKM host stars generally have kilogauss magnetic fields which can exert torques on the orbits of nearby super- Earths. In this work, we examine one aspect of this interaction: the magnetic torque resulting from Alfvén-wing drag on non-corotating, non-magnetized planets engulfed by the host stars' stellar wind. We compute the magnitude of this torque for a range of stellar magnetic field strengths, and planetary orbital velocities. We also model the planets' orbital evolution, taking into account for stellar spin down and magnetic field decay, and derive the boundaries within which ultra-short-period super-Earths can survive.

High-velocity clouds (HVCs) are composed of neutral hydrogen (HI) moving at velocities that deviate from the general rotational motion of the Milky Way. Currently, the origins of the HVCs remain poorly known due to the difficulty in determining their distance and the lack of any other suitable identification. Here we report the detection of a compact gas clump in HVC AC-I, which displays characteristics typical of a disk galaxy, named AC G185.0-11.5, using the HI observations. We estimated the distance of AC G185.0-11.5 to be about 277.7 kpc using the Baryonic Tully-Fisher relation and constrained its HI gas mass to be between 3.0*10^7 and 4.7*10^8 solar masses. The distance determination indicates that the HVC AC-I hosting AC G185.0-11.5 is an extragalactic object in the Local Group. The absence of molecular gas and an optical counterpart for AC G185.0-11.5 implies that it may be a rare dark galaxy.

A discrete jet component (blob) ejection and its subsequent deceleration was observed in the 2019/2020 outburst in the low-mass X-ray binary MAXI J1348-630. A first kinematic analysis of the deceleration due to an abrupt transition from an evacuated cavity to the interstellar medium suggested a kinetic energy exceeding 1046 erg, surpassing estimates of the available total ejection energy. However, incorporating a transition layer with exponential density growth between the cavity and interstellar medium recently enabled a kinematic analysis with much more realistic energy requirements of approximately $10^{44}$ erg. Here, we study the expected radiative signatures of electrons accelerated within the decelerating blob by introducing a model akin to the relativistic blast wave model for gamma-ray bursts, considering radiative energy losses and radiation drag, to simulate the deceleration of a relativistically moving plasmoid. This model yields snap-shot spectral energy distributions and multi-wavelength light curves from synchrotron and Synchrotron-Self-Compton (SSC) emission. Notably, the synchrotron emission peaks in the X-rays, and the predicted Radio and X-ray light curves closely resemble the observed ones during the jet decleration phase following the outburst in 2019/2020.

This work presents a new detailed study on the energy-dependent variation in the X-ray polarisation of the Crab Pulsar Wind Nebula (PWN), obtained using data from the Imaging X-ray Polarimetry Explorer (IXPE). For the entire PWN, we observed a linear variation in polarisation degree (PD), and detected the rotation of the polarisation angle (PA) with the energy at higher than 99.9999\% of the confidence level. This energy-dependent polarisation variation is in line with the indication found in Vela PWN by IXPE, and it can be interpreted as the emitting region of the polarised photons shrinks with increasing energy, leading to higher PD because they are less influenced by the turbulence of the magnetic field. We compared the IXPE polarisation results with those of other hard X-ray/gamma observatories (PoGO+, Intregral, AstroSat) for the PWN, finding the same trend from soft-X to hard-X with the PD increasing with the energy and the PA approaching the pulsar's spin axis. In fact, in this wide energy band, the fitting results show an energy trend for the PA compatible with the estimated pulsar's spin axis within 3$\sigma$ of confidence level.

The Gaia experiment recently reported the observation of a binary system composed of a Sun-like star orbiting a dark compact object, known as Gaia BH1. The nature of the compact object remains uncertain. While the Gaia mission identifies it as a black hole candidate, the absence of X-ray or radio detections challenges that interpretation, and alternative exotic compact objects such as boson stars have also been suggested. In this paper, we study whether a boson star could account for the observed properties of the source. To do so we compute the X-ray luminosity of the central dark object as a result of spherically symmetric (Bondi-Michel) accretion of matter, comparing our results for the cases in which the dark object is a Schwarzschild black hole or a non rotating boson star. Our model incorporates realistic interstellar medium properties, ranging from hot ionized gas to dense molecular clouds. By solving the governing equations numerically, we calculate mass accretion rates and derive the resulting Bremsstrahlung X-ray luminosities. Black holes and boson stars fundamentally differ by the absence of an event horizon in the latter, which directly impacts accretion dynamics as there is an accumulation of mass in regions closer to the boson star, which will significantly change the observed X-ray emission. For the Gaia BH1 system we find that accretion onto a black hole yields luminosities of $\sim10^{27} \ \text{erg}\, \text{cm}^{-2}\, \text{s}^{-1}$ which corresponds to an X-ray flux undetectable by Chandra sensitivity. On the other hand, boson star accretion can produce observable luminosities in the order of $10^{27} \ \text{to} \ 10^{41} \ \text{erg}\, \text{cm}^{-2}\, \text{s}^{-1}$.

We discuss an interacting decaying vacuum energy and dark matter empowered by gravitationally induced matter creation model, and its impact on structure formation by analysing the growth rate perturbations. Our work is motivated by the possibility that the decaying vacuum is due to quantum field theory and dark matter originates from gravitationally induced matter creation. We delve deeper into our investigation and explore both theoretical and statistical analysis of the cosmological model to test its ability to describe the evolution of the Universe. To achieve this, we use three distinct combinations of datasets from CC, Pantheon SNIa sample, BAO, CMB distance priors and $f(z)\sigma_8(z)$ datapoints to constrain the model parameters. Our statistical analysis employs Markov Chain Monte Carlo (MCMC) methods. The deceleration parameter shows that the model transitions from a decelerated phase to an accelerated phase of expansion. The current Hubble parameter values are estimated to be $H_0=67.517 \pm 0.869$ km/s/Mpc, $H_0=67.534\pm 0.874$ km/s/Mpc, and $H_0=67.533 \pm 0.884$ km/s/Mpc using DS1, DS2 and DS3 datasets, respectively. These values of $H_0$ are very close to those derived from the Planck data. The effective equation of state parameter indicates an accelerating phase, with density parameter for vacuum energy exhibiting expected values. We analyse the stability characteristics through the selection information criteria. We also perform thermodynamic analysis by studying the evolution of entropy in the Universe for the model and find it to be in agreement with the generalized second law of thermodynamics. These findings support that the proposed model effectively describes the evolutionary features of the Universe at both theoretical and observational levels.

We investigate the formation of primordial black holes in curvaton models of inflation, where the curvature perturbation is not only generated by the inflaton but also by a light scalar field (the curvaton) that decays after inflation. During inflation, both fields are subject to quantum diffusion, owing to small-scale vacuum fluctuations crossing out the Hubble radius. After inflation, whether the curvaton dominates the universe or not depends on its field value when inflation ends. Since that value is stochastic, different regions of the universe undergo different post-inflationary histories. In practice, we show that this results in a double-peaked distribution for the number of e-folds realised in these models. Since that number of e-folds is related to the curvature perturbation by the delta-N formalism, the presence of a second peak has important consequences for primordial black holes that we discuss.

We examined the spatially resolved polarization variations in the Crab Nebula over 2 yr, using observational data from the Imaging X-ray Polarimetry Explorer, and offer key insights into its magnetic field structures and evolution. The results show significant temporal changes in the polarization degree (PD) across three regions of interest in the 2-8 keV energy band. Regions (a) and (b), located in the northern and the southwestern parts of the study area, exhibit PD variations with significance levels greater than 4 sigma and 3 sigma , respectively. Region (c), located in the southwest,shows a notable decrease in PD with a significance greater than 5 sigma. However, no significant variation in the polarization angle was observed. Meanwhile, notable flux variations were detected, likely influenced by dynamic processes such as magnetized turbulence within the nebula.

It is important to explore the potential existence of lepton asymmetry in the neutrino sector. Conducting a joint analysis of DESI DR2 BAO data and \emph{Planck} 2018 CMB data, we obtain the upper limits on the neutrino degeneracy parameter, i.e., $\xi_{3}<0.56$ for the normal mass hierarchy while $\xi_{3}<0.62$ for the inverted mass hierarchy, at 95\% confidence level. Considering the influence of the dynamical dark energy, we find that these upper limits remain to be robust. This work may provide helpful implications for model buildings of the matter-antimatter asymmetry in the universe.

This work discusses the influence of galaxy mergers in the evolution of a parabolic Lema\^ıtre-Tolman-Bondi (LTB) cosmology with simultaneous big bang endowed with two consecutive single fractal galaxy distributions systems possessing fractal dimension $D$. Based on recent empirical findings, it is assumed that the resulting galaxy mass from mergers can be expressed by a redshift dependent decaying power law. The proposed cosmological model modifies the relativistic fractal number counts distribution by including a merger rate evolution that estimates the model's radial density. Numerical solutions for the first order small-merger-rate approximation (SMRA) are found and the results show that a fractal galaxy distribution having $D=1.5$ in the range $0.1<z<1.0$, and $D=0.5$ for $1<z<6$, as suggested by recent empirical findings, the SMRA allows consistent description of the model for a merger rate power law exponent up to $q=0.2$ considering a fractal galaxy distribution starting from the Local Group. Consistent values were also found up to $q=2.5$ and $z=7$ from a scale smaller than the Local Supercluster. These results show that galaxy mergers can be successfully incorporated into the dynamics of a parabolic LTB fractal cosmology.

This study focuses on high-redshift, z > 3, quasars where resolved X-ray jets remain underexplored in comparison to nearby sources. Building upon previous work, we identify and confirm extended kpc-scale jets emission in two quasars (J1405+0415, z = 3.215 ; J1610+1811, z = 3.122) through meticulous analysis of Chandra X-ray data. To deepen our understanding, high-resolution radio follow- up observations were conducted to constrain relativistic parameters, providing valuable insights into the enthalpy flux of these high-redshift AGN jets. The investigation specifically aims to test the X-ray emission mechanism in these quasars by exploring the inverse Compton scattering of cosmic microwave background (IC/CMB) photons by synchrotron-emitting electrons. Our novel method uses a prior to make a Bayesian estimate of the unknown angle of the jet to our line of sight, thus breaking the usual degeneracy of the bulk Lorentz factor and the Doppler beaming factor.

Massive early-type galaxies (ETGs) are thought to form in two phases: an initial phase of rapid star formation and a later phase of mergers. A small fraction of these galaxies, called red nuggets, formed during the first phase may have survived frozen until today, experiencing no massive mergers since z~2. Nearby massive compact ETGs are considered candidates for such relic galaxies. We study the internal dynamical structures of 15 compact ETGs with existing integral field unit (IFU) observations and 79 compact ETGs from the TNG50 simulation. We dynamically decompose each galaxy into a disk, bulge, and hot inner stellar halo, for both observations and simulations. In TNG50, the luminosity fraction of the hot inner stellar halo (or the size of the spheroid, which includes the bulge and halo) strongly correlates with the galaxy merger history. The true merger-free galaxies show an extremely low fraction of a hot inner stellar halo (or an extremely compact spheroid). Although such compactness could result from the tidal stripping of satellites, tidal forces would also destroy the dynamically cold disk (if one exists) when the halo is removed. Thus, a galaxy is guaranteed to be merger-free if it has a very low fraction of the hot inner stellar halo and retains a dynamically cold disk. Comparing observed galaxies with TNG50, we identify 7 of the 15 compact ETGs, PGC 11179, UGC 3816, NGC 2767, NGC 1277, PGC 32873, PGC 12562, and PGC 70520, as true merger-free galaxies. These galaxies have compact, massive bulges, likely formed through secular heating, as supported by their TNG50 analogues.

The DESI DR2 BAO data exclude the flat $\Lambda$CDM model at more than 2.5$\sigma$, depending on different data combinations when analyzed through the $w_0w_a$CDM parametrization. This simple parametrization may introduce bias in the results. We use null tests that probe for deviations from flat $\Lambda$CDM at late times, independent of any specific dark energy parametrization. We provide several diagnostics for null tests and discuss their advantages and disadvantages. In particular, we derive diagnostics that improve on previous ones, such as the popular $O_{\rm m}$ diagnostic. The diagnostics are derived from both background and perturbed quantities. Using the combination of DESI DR2 BAO and supernova data, with or without CMB data, we find that deviations from flat $\Lambda$CDM are at $\sim1\sigma$ confidence level in most of the redshift range (more than 1$\sigma$ for a few small redshift intervals in a few cases). These deviations are minor for other non-DESI SDSS-IV BAO data combined with Pantheon+, with or without CMB data. Since spatial curvature can potentially modify the results, we also test for curvature in the general $\Lambda$CDM model and the general FLRW model. While there is slight evidence for nonzero cosmic curvature at lower redshifts in a general $\Lambda$CDM model, there is no statistically significant evidence in a general FLRW model.

Faraday rotation describes the change of the linear polarization angle of radiation passing through a magnetized plasma and it is quantified by the rotation measure (RM), which is related to the line-of-sight (LOS) magnetic field component and the thermal electron density traversed by light along its path toward the observer. However, it is challenging to disentangle the signal from different LOS portions and separate the contribution from the local ISM. This is particularly relevant since the Sun is located within the Local Bubble (LB), a low-density and hot cavity formed by past SN events, making it essential to investigate how this environment may influence the observed RM values. The present study investigates the imprint of the local environment on the synthetic RM signal, as measured by an observer within a LB-like cavity. The RM derived from diffuse polarized synchrotron radiation produced by CR electrons at decimeter wavelengths is also analyzed. We produce synthetic RM maps for an observer placed inside a LB candidate, selected from a MHD simulation that resembles the properties of the ISM in the Solar vicinity. Using the capabilities of the radiative transfer code POLARIS, we study the imprint of the cavity walls on the RM signal. As the MHD simulation does not account for CR diffusion, we develop a CR toy-model to study the Faraday rotation of the diffuse polarized synchrotron radiation. We find that (i) the imprint of local structures, such as the walls of the LB candidate and the edges of other supernovae blown cavities, is of fundamental importance for interpreting the global Faraday sky; (ii) the LB has a non negligible contribution to the sinusoidal patterns of RM as a function of Galactic longitude seen in observations; and (iii) the RM signal from diffuse synchrotron emission shows a strong correspondence with the RM signal generated by the LB candidate walls.

Measuring the magnetic field of the Milky Way reveals the structure and evolution of the galaxy. Pulsar rotation measures (RMs) provide a means to probe this Galactic magnetic field (GMF) in three dimensions. We use the largest single-origin data set of pulsar measurements, from the MeerKAT Thousand-Pulsar-Array, to map out GMF components parallel to pulsar lines of sight. We also present these measurements for easy integration into the consolidated RM catalogue, RMTable. Focusing on the Galactic disk, we investigate competing theories of how the GMF relates to the spiral arms, comparing our observational map with five analytic models of magnetic field structure. We also analyse RMs to extragalactic radio sources, to help build up a three-dimensional picture of the magnetic structure of the galaxy. In particular, our large number of measurements allows us to investigate differing magnetic field behaviour in the upper and lower halves of the Galactic plane. We find that the GMF is best explained as following the spiral arms in a roughly bisymmetric structure, with antisymmetric parity with respect to the Galactic plane. This picture is complicated by variations in parity on different spiral arms, and the parity change location appears to be shifted by a distance of 0.15 kpc perpendicular to the Galactic plane. This indicates a complex relationship between the large-scale distributions of matter and magnetic fields in our galaxy. Future pulsar discoveries will help reveal the origins of this relationship with greater precision, as well as probing the locations of local magnetic field inhomogenities.

Observations of the Epoch of Reionization (EoR) have the potential to answer long-standing questions of astrophysical interest regarding the nature of the first luminous sources and their effects on the intergalactic medium (IGM). We present astrophysical constraints from a Neural Density Estimation-Accelerated Bayesian joint analysis of constraints deriving from Cosmic Microwave Background power spectrum measurements from Planck and SPT, IGM neutral fraction measurements from Lyman-line-based data sets and 21-cm power spectrum upper limits from HERA, LOFAR and the MWA. In the context of the model employed, the data is found to be consistent with galaxies forming from predominantly atomic-cooled hydrogen gas in dark matter halos, with masses $M_\mathrm{min} \gtrsim 2.6 \times 10^{9}~M_{\odot} ((1+z)/10)^{\frac{1}{2}}$ at 95% credibility ($V_\mathrm{c} \gtrsim 50~\mathrm{km~s^{-1}}$) being the dominant galactic population driving reionization. These galaxies reionize the neutral hydrogen in the IGM over a narrow redshift interval ($\Delta z_\mathrm{re} < 1.8$ at 95% credibility), with the midpoint of reionization (when the sky-averaged IGM neutral fraction is 50%) constrained to $z_{50} = 7.16^{+0.15}_{-0.12}$. Given the parameter posteriors from our joint analysis, we find that the posterior predictive distribution of the global 21-cm signal is reduced in amplitude and shifted to lower redshifts relative to the model prior. We caution, however, that our inferences are model-dependent. Future work incorporating updated, mass-dependent star formation efficiencies in atomic cooling halos, informed by the latest UV luminosity function constraints from the James Webb Space Telescope, promises to refine these inferences further and enhance our understanding of cosmic reionization.

We present a new approach for understanding how galaxies lose or retain baryons by utilizing a pipeline of two machine learning methods applied the IllustrisTNG100 simulation. We employed a Random Forest Regressor and Explainable Boosting Machine (EBM) model to connect the retained baryon fraction of $\approx10^5$ simulated galaxies to their properties. We employed Random Forest models to filter and used the five most significant properties to train an EBM. Interaction functions identified by the EBM highlight the relationship between baryon fraction and three different galactic mass measurements, the location of the rotation curve peak, and the velocity dispersion. This interpretable machine learning-based approach provides a promising pathway for understanding the baryon cycle in galaxies.

Photometric and spectroscopic observations of GV Leo were performed from 2017 to 2024. The light curves show a flat bottom at the primary eclipse and the conventional O'Connell effect. The echelle spectra reveal that the effective temperature and rotation velocity of the more massive secondary are $T_{\rm eff,2}$ = 5220$\pm$120 K and $v_2 \sin i$ = 223$\pm$40 km s$^{-1}$, respectively. Our binary modeling indicates that the program target is a W-subclass contact binary with a mass ratio of $q$ = 5.48, an inclination angle of $i$ = 81$^\circ$.68, a temperature difference of ($T_{\rm eff,1}-T_{\rm eff,2}$) = 154 K, and a filling factor of $f$ = 36 \%. The light asymmetries were reasonably modeled by a dark starspot on the secondary's photosphere. Including our 26 minimum epochs, 84 times of minimum light were used to investigate the orbital period of the system. We found that the eclipse times of GV Leo have varied by a sinusoid with a period of 14.9 years and a semi-amplitude of 0.0076 days superimposed on a downward parabola. The periodic modulation is interpreted as a light time effect produced by an unseen outer tertiary with a minimum mass of 0.26 M$_\odot$, while the parabolic component is thought to be a combination of mass transfer (secondary to primary) and angular momentum loss driven by magnetic braking. The circumbinary tertiary would have caused the eclipsing pair of GV Leo to evolve into its current short-period contact state by removing angular momentum from the primordial widish binary.

The nature of sub-Neptunes is one of the hottest topics in exoplanetary science. Temperate sub-Neptunes are of special interest because some could be habitable. Here, we consider whether these planets might instead be rocky worlds with thick, hot atmospheres. Can recent JWST observations of TOI-270 d be understood in terms of such a model? We perform thermochemical equilibrium calculations to infer conditions of quenching of C-H-O-N species. Our results indicate apparent CO$_2$-CH$_4$ equilibrium between ~900 and ~1100 K. The CO abundance should be quenched higher in the atmosphere where the equilibrium CO/CO$_2$ ratio is lower, potentially explaining a lack of CO. N$_2$ is predicted to dominate the nitrogen budget. We confirm that the atmosphere of TOI-270 d is strongly enriched in both C and O$_g$$_a$$_s$ relative to protosolar H, whereas N is likely to be less enriched or even depleted. We attempt to reproduce these enrichments by modeling the atmosphere as nebular gas that extracted heavy elements from accreted solids. This type of model can explain the C/H and O$_g$$_a$$_s$/H ratios, but despite supersolar C/N ratios provided by solids, the NH$_3$ abundance will probably be too high unless there is a nitrogen sink in addition to N$_2$. A magma ocean may be implied, and indeed the oxygen fugacity of the deep atmosphere seems sufficiently low to support the sequestration of reduced N in silicate melt. The evaluation presented here demonstrates that exoplanetary geochemistry now approaches a level of sophistication comparable to that achieved within our own solar system.

Scaling relations are fundamental tools for exploring the morphological properties of galaxies and understanding their formation and evolution. Typically, galaxies follow a scaling relation between mass and size, measured by effective radius. However, a compact class of galaxies exists as outliers from this relation, and the origin of these compact galaxies in the local universe remains unclear. In this study, we investigate the compact dwarf galaxy SDSS J134313.15+364457.5 (J1343+3644), which is the result of a merger. Our analysis reveals that J1343+3644 has a half-light radius of 482~pc, significantly smaller than typical galaxies with the same brightness ($M_\text{r} = -19.17$ mag). With a high star-formation rate (SFR) of 0.87~M$_{\sun}$ year$^{-1}$, J1343+3644 is expected to evolve into a compact elliptical galaxy in a few million years. J1343+3644 could, therefore, be a progenitor of a compact elliptical galaxy. The phenomenon happened in early universe, where compact galaxies were common.

Mass-loss influences stellar evolution, especially for massive stars with strong winds. Stellar wind bow shock nebulae driven by Galactic OB stars can be used to measure mass-loss rates ($\dot{M}$). The standoff distance ($R_{0}$) between the star and the bow shock is set by momentum flux balance between the stellar wind and the surrounding interstellar medium (ISM). We created the Milky Way Project: MOBStIRS (Mass-loss rates for OB Stars driving IR bow Shocks) using the online Zooniverse citizen science platform. We enlisted several hundred students to measure $R_0$ and two other projected shape parameters for 764 cataloged IR bow shocks. MOBStIRS incorporated 1528 JPEG cutout images produced from Spitzer GLIMPSE and MIPSGAL survey data. Measurements were aggregated to compute shape parameters for each bow shock image deemed high-quality by participants. The average statistical uncertainty on $R_0$ is $12.5\%$ but varies from ${<}5\%$ to ${\sim}40\%$ among individual bow shocks, contributing significantly to the total error budget of $\dot{M}$. The derived nebular morphologies agree well with (magneto)hydrodynamic simulations of bow shocks driven by the winds of OB stars moving at $V_a = 10-40~km~s^{-1}$ with respect to the ambient interstellar medium (ISM). A systematic correction to $R_0$ to account for viewing angle appears unnecessary for computing $\dot{M}$. Slightly more than half of MOBStIRS bow shocks are asymmetric, which could indicate anisotropic stellar winds, ISM clumping on sub-pc scales, time-dependent instabilities, and/or misalignments between the local ISM magnetic field and the star-bow shock axis.

Towering storms, swirling clouds, and vortices are the cloud tops manifestation of complex weather systems shaping the atmosphere of Jupiter. We use observations from Juno's MicroWave Radiometer (MWR), the Very Large Array (VLA) and the Hubble Space Telescope (HST) to probe for the first time the depth and impact of weather on Jupiter. We use ammonia, the main source of opacity at radio wavelengths on Jupiter, as the tracer for the weather by fitting ammonia anomalies to the MWR brightness temperature variations. We show that the majority of the weather on Jupiter is confined to regions where the clouds are forming. Both the South Equatorial Belt and the Equatorial Zone have surprisingly shallow weather systems (P < 2 bar), and even in the North Equatorial Belt most of the ammonia variations is above the water condensation level (P ~ 6 bar). This confirms that the water condensation layer plays a crucial role in controlling the dynamics and the weather on Jupiter. However, the shallow nature of the weather cannot explain the deep-seated depletion down to 30 bar that the Juno mission has revealed. We do find three features, however, that extend below the water condensation layer: a vortex in the northern hemisphere reaching down to 30 bar, an ammonia plume down to 20-30 bars, and the signature of precipitation down to 20 bar. This work highlights the interplay of large-scale processes (vortices, plumes) and small-scale processes (storms) are responsible for shaping the atmospheric makeup of Jupiter.

We utilize the DESI Legacy Imaging Surveys DR10 to investigate the previously undetected faint extension of the Palomar 5 stellar stream. By applying the HDBSCAN clustering algorithm, we identify stream members and successfully extend the leading arm of the stream to approximately $\mathrm{DEC} \sim -15^\circ$. Combining the fully detected stream with a suite of mock stream simulations, we conduct a detailed comparison to constrain both the intrinsic properties of the stream and the dynamical parameters of the Milky Way (MW) halo. Our analysis yields a best-fit model characterized by eight parameters: $M_{\mathrm{halo}} = 5.67\times10^{11}\ M_{\odot}$, $r_{s,\mathrm{halo}} = 28.94\ \mathrm{kpc}$, $q_z = 0.93$, $M_{\mathrm{gc}} = 4.31\times10^{3}\ M_{\odot}$, $dM_{\mathrm{gc}}/dt = 1.81\ M_{\odot}\ \mathrm{Myr}^{-1}$, $\mu_{\alpha}\cos\delta = -2.28\ \mathrm{mas\ yr}^{-1}$, $\mu_{\delta} = -2.26\ \mathrm{mas\ yr}^{-1}$, and $D = 23.25\ \mathrm{kpc}$. Notably, our constraints on the halo shape indicate that the MW's dark matter halo exhibits a flattened potential, with a minor-to-major axis ratio of $q_z = 0.93$. This finding aligns well with theoretical expectations and previous observational estimates. Additionally, the best-fit model accurately reproduces the observed stream morphology and dynamics, providing a more precise understanding of both the evolution of the stream and the overall structure of the Galactic halo.

We present a complete treatment of the bispectrum of intrinsic alignments, both in three spatial dimensions and in projection in the flat-sky approximation. Since intrinsic alignment is a spin-2 observable, the bispectrum of intrinsic alignments contains a parity-even and a parity-odd part, the latter being nonzero even in the absence of parity violation. Moreover, all possible combinations of scalar, E- and B-mode bispectra are nonzero in the absence of parity violation. In analogy to the galaxy bispectrum in redshift space, we construct a complete set of multipoles for anisotropic bispectra of projected spin-2 fields. We then construct separable bispectrum estimators, both for parity-even and parity-odd bispectra, which can be computed by means of Fast Fourier Transforms (FFTs). We compare several different choices of angular weighting in terms of signal-to-noise ratios (SNR) for a Stage IV setup using luminous red galaxies (LRGs) from the Dark Energy Spectroscopic Instrument (DESI) with galaxy shapes measured by the Legacy Survey of Space and Time (LSST). Assuming an overlapping area of $\sim 4,000$ square degrees (yielding $\sim 1.3$ million LRGs) and including scales up to $k_\text{max} = 0.14\,h$/Mpc, we find that the position-position-E-mode bispectrum $B_{DDE}$ (which is parity-even) can be strongly detected at SNR $\sim 30$, while detecting parity-odd bispectra (such as $B_{DDB}$, SNR $\sim 5$) or bispectra with more than one shape field (such as $B_{DEE}$, SNR $\sim 5$) may also be possible.

This study evaluates the impact of Starlink satellites on low-frequency radio astronomy below 100 MHz, focusing on challenges on data processing and scientific goals. We conducted 40 hours of imaging observations using NenuFAR, in the 30.8-78.3 MHz range. Observations included both targeted tracking of specific satellites based on orbital predictions and untargeted searches focused on high-elevation regions of the sky. Images in total intensity and polarimetry were obtained, and full Stokes dynamic spectra were generated for several hundred directions within the Field of View. Detected signals were cross-matched with satellite orbital data to confirm satellite associations. Detailed analyses of the observed spectra, polarization, and temporal characteristics were performed to investigate the origin and properties of the detected emissions. We detected broadband emissions from Starlink satellites, predominantly between 54-66 MHz, with flux densities exceeding 500 Jy. These signals are highly polarized and unlikely to originate from ground-based RFI or reflected astronomical sources. Instead, they are likely intrinsic to the satellites, with distinct differences in emission properties observed between satellite generations. These findings highlight significant challenges to data processing and scientific discoveries at these low frequencies, emphasizing the need for effective mitigation strategies, particularly through collaboration between astronomers and satellite operators.

Spider pulsars are systems in which a millisecond pulsar (MSP) tightly orbits (Pb $\lesssim$ 1 day) a low mass (mc $\lesssim$ 0.5 M$_\odot$) semi-degenerate star. Spider often display eclipses around superior conjunction. This eclipse phenomenon is currently poorly understood. We analyzed eclipses via pulsar timing. The eclipses were fit with a phenomenological model which gives a measurement of the duration and asymmetry of the eclipses. These parameters were then compared to other eclipse and system measurements to discuss the potential link between the presence of eclipses and orbital inclination, eclipsing systems being known to have higher mass functions than non-eclipsing ones. We present here a comprehensive review of the NRT NUPPI backend spider pulsars dataset. We also present the first review and systematic analysis of a large sample of eclipsers, monitored with the NRT over several years. The phenomenological fit allowed us to compare the eclipsers with each other, which led to the categorization of eclipsers depending on the shape of their eclipses. We present the polarimetric properties of the 19 spiders in the sample alongside their profiles, which were previously unpublished in some cases. For the eclipsing systems, we found evidence for a positive correlation between eclipse duration and mass function, as expected if more eclipsing material crosses the line-of-sight in higher inclination systems. For the entire sample, we found marginal evidence for increasing pulse profile width with decreasing mass function. We finally conducted a comprehensive literature review of the published inclination measurements for the pulsars in the sample and compared the inclinations to eclipse parameters. Nevertheless, the small number of available orbital inclination constraints, contradicting each other in some cases, hinders such searches for correlations

The rapid advancement of observational capabilities in astronomy has led to an exponential growth in the volume of light curve (LC) data, presenting both opportunities and challenges for time-domain astronomy. Traditional analytical methods often struggle to fully extract the scientific value of these vast datasets, especially as their complexity increases. Machine learning (ML) algorithms have become indispensable tools for analyzing light curves, offering the ability to classify, predict, discover patterns, and detect anomalies. Despite the growing adoption of ML techniques, challenges remain in LC classification, including class imbalance, noisy data, and interpretability of models. These challenges emphasize the importance of conducting a systematic review of ML algorithms specifically tailored for LC analysis. This survey provides a comprehensive overview of the latest ML techniques, summarizing their principles and applications in key astronomical tasks such as exoplanet detection, variable star classification, and supernova identification. It also discusses strategies to address the existing challenges and advance LC analysis in the near future. As astronomical datasets continue to grow, the integration of ML and deep learning (DL) techniques will be essential for unlocking the full scientific potential of LC data in the era of astronomical big data.

Ground-level particle detection has recently emerged as an extremely powerful approach to TeV-PeV gamma-ray astronomy. The most successful observatories of this type, HAWC and LHAASO, utilise water-Cherenkov based detector units, housed in tanks or buildings. Here we explore the possibility of deploying water-Cherenkov detector units directly in to a natural or artificial lake. Possible advantages include reduced cost and improved performance due to better shielding. The lake concept has been developed as an option for the future Southern Wide-view Gamma-ray Observatory, and is now under consideration for a possible future extension of the observatory, beyond its recently selected land site. We present results from prototypes operated in a custom built facility, and concepts for full-scale array deployment and long-term operation.

Previous studies have observed significant photometric differences between non-emission B-type and classical Be stars, however the precise mechanism responsible for these differences is unclear. This study combines the Bright Star Catalogue with Tycho and Gaia photometry to create a homogeneous sample of 1015 of the closest and brightest B and Be-type field stars with 90 per cent of objects at distances < 500pc. Due to their proximity, the extinction towards these objects is very low, ensuring we minimise any obfuscation in the reddening correction and final photometry. We present our findings in both Tycho and Gaia photometry through colour magnitude diagrams and present intrinsic colours and absolute magnitudes for each spectral type. We find Be stars are on average ~0.5 magnitudes brighter in both Gaia $G$ and Tycho V$_T$ compared to non-emission B stars of the same spectral type. Additionally, we find tentative evidence that Be stars are redder in Gaia B$_P$$-$R$_P$, particularly for the earlier types, but have similar Tycho B$_T$$-$V$_T$ colours. We test the effects of gravitational darkening due to rapid rotation and binarity on the photometry of our sample and find both to be insufficient to explain the observed photometric differences between B and Be stars. We conclude that the most likely mechanism responsible for the observed photometric differences is the combined effect of the circumstellar disc and stellar evolution up the Main Sequence, with the disc dominating early-types and evolution dominating late type stars.

We investigated the potential of polarimetric observations in the optical wavelength range for the detection of exomoons and the characterization of exoplanet-exomoon systems. Using the three-dimensional Monte Carlo radiative transfer code POLARIS, we calculated flux and polarization phase curves of Earth-like exoplanets with a satellite similar to Earth's moon. Of particular interest are mutual events, when one of the two bodies casts a shadow on the other or transits in front of it. We find that the signatures of mutual events in the polarization phase curve show significant variations depending on the inclination of the lunar orbit. If the planet-satellite pair is spatially resolved from the star but the satellite is spatially unresolved, the increase in the degree of polarization during a transit of the exomoon in front of the center of the exoplanet reaches $2.7\%$ in our model system near quadrature. However, the change is less than $0.5\%$ if the orbit of the exomoon is inclined such that it transits the planet noncentrally at the same phase angles. The influence of an exomoon on the polarization phase curve of an exoplanet-exomoon system is dependent on the lunar polarization phase curve. Observations of full eclipses and occultations of the exomoon allow the determination of separate polarization phase curves for the two bodies. Information about the lunar orbital inclination can be obtained with polarimetric observations of shadows or transits. Measuring the influence of large satellites not only on the total flux, but also on the polarization of the reflected stellar radiation during mutual events thus facilitates the prediction of future mutual events and the verification of exomoon candidates.

Studying chemically rich protostellar outflows and their jet provides an important insight into the low-mass star formation process and its related chemistry. Whilst well-known shock tracers such as SiO can be used to study the jet properties and give information about the dynamics of the system, interstellar complex organic molecules (iCOMs) have been useful in constraining the age of shocked gas, for example. Yet, the number of outflows mapped in iCOMs is still limited. In this work, we study the outflow driven by the protostar FIR6c-a (HOPS 409) located in the OMC-2/3 filament. We report the detection of the red-shifted jet, left undetected in previous studies, as well as the detection of the iCOMs methanol (CH$_3$OH) and methyl cyanide (CH$_3$CN) for the first time towards this outflow. Using SiO, we derived some jet properties (i.e., collimation and dynamical time). We found a clear dichotomy between the blue- and red-shifted jets, likely due to the density of the medium in which the jets propagate. In addition, we identified two bow shocks within the blue-shifted part of the outflow, which we attribute to two different ejection events. Finally, using the CH$_3$OH} and \ce{CH$_3$CN} abundance ratio and chemical modelling, we constrained the outflow age to be $\geq 1000$ yr old and, surprisingly, found that a cosmic-ray ionization rate of $10^{-14}$ s$^{-1}$ is needed to reproduce the observed ratio towards the source.

We aim to develop a new method to infer the sub-beam probability density function (PDF) of H2 column densities and the dense gas mass within molecular clouds using spatially unresolved observations of molecular emission lines in the 3 mm band. We model spatially unresolved line integrated intensity measurements as the average of an emission function weighted by the sub-beam column density PDF. The emission function, which expresses the line integrated intensity as a function of the gas column density, is an empirical fit to high resolution (< 0.05 pc) multi-line observations of the Orion B molecular cloud. The column density PDF is assumed to be parametric, composed of a lognormal distribution at moderate column densities and a power law distribution at higher column densities. To estimate the sub-beam column density PDF, the emission model is combined with a Bayesian inversion algorithm (the Beetroots code), which takes account of thermal noise and calibration errors. We validate our method by demonstrating that it recovers the true column density PDF of the Orion B cloud, reproducing the observed emission line integrated intensities. We apply the method to 12CO(J=1-0), 13CO(J=1-0), C18O(J=1-0), HCN(J=1-0), HCO+(J=1-0) and N2H+(J=1-0) observations of a 700 x 700 pc2 field of view (FoV) in the nearby galaxy M51. On average, the model reproduces the observed intensities within 30%. The column density PDFs obtained for the spiral arm region within our test FoV are dominated by a power-law tail at high column densities, with slopes that are consistent with gravitational collapse. Outside the spiral arm, the column density PDFs are predominantly lognormal, consistent with supersonic isothermal turbulence. We calculate the mass associated with the powerlaw tail of the column density PDFs and observe a strong, linear correlation between this mass and the 24$\mu$m surface brightness.

Gamma-ray bursts (GRBs) are the most powerful explosions in the Universe; their energy release reache s us from the end of the re-ionization era, making them invaluable cosmological probes. GRB 230307A i s the second-brightest GRB ever observed in the 56 years of observations since the discovery of the phenomenon in 1967. Follow-up observations of the event at longer wavelengths revealed a lanthanide-ri ch kilonova with long-lasting X-ray emission immediately following the prompt gamma-rays. Moreover, t he gamma-ray light curve of GRB 230307A collected with INTEGRAL's SPectrometer of INTEGRAL AntiCoincidence Shield (SPI-ACS) and Fermi's Gamma-Ray Burst Monitor (GBM). We use Fourier analysis, wavelets and Gaussian Processes to search for periodic and quasi-periodic oscillations (QPOs) in the prompt gamma-ray emission of GRB 230307A. We critically assess all three methods in terms of their robustness for detections of QPOs in fast transients such as GRBs. Our analyses reveal QPOs at a frequency of $\sim 1.2$ Hz (0.82s period) near the burst's peak emission phase, consistent across instruments and detection methods. We also identify a second, less significant QPO at $\sim 2.9$ Hz (0.34s) nearly simultaneously. We hypothesise that the two QPOs originate from the transition epoch at the end of the jet acceleration phase. These QPOs re present plasma circulation periods in vorticity about the jet axis carried outwards to the prompt radiation zone at much larger radii. They are sampled by colliding structures (e.g., shocks) in the spinning jet, possibly marking the evolution of plasma rotation during the final stages of the progenitor neutron star coalescence event.

Alfven waves propagating in a vertically stratified plasma, such as those travelling from the solar photosphere to the corona, are partially reflected due to the gradient in the Alfven speed. Wave reflection naturally results in the superposition of upward- and downward-propagating waves. A simple analytic model demonstrates that this superposition leads to the non-equipartition of kinetic and magnetic energies in the Alfven wave perturbations and slows down the net energy transport. A numerical model of Alfven wave propagation in the lower solar atmosphere reveals significant wave reflection below the transition region, leading to highly variable kinetic and magnetic energy content in the lower chromosphere. At higher altitudes, the kinetic energy eventually dominates, depending on the wave frequency. The velocity at which net energy propagates upward is significantly smaller than the local Alfven speed throughout the chromosphere. Consequently, the commonly used expression for unidirectional Alfven waves in a uniform plasma, which relates the energy flux to the kinetic energy density, is not generally applicable in the stratified lower solar atmosphere and cannot be reliably used to estimate the energy content of observed waves. A generalized expression is given, incorporating correction factors that account for wave reflection and energy non-equipartition. The applicability of the expression is discussed.

Accretion of pre-main sequence stars (PMS) is a key process in stellar formation, governing mass assembly, influencing angular momentum conservation and stellar internal structure, and shaping disc evolution, which serves as the birthplace of exoplanets. Classical T Tauri stars (cTTSs), low-mass PMS stars actively accreting from a disc, hold a well-described magnetospheric accretion model. Their strong, inclined dipole magnetic fields truncate the disc at a few stellar radii, channelling material along magnetic field lines to fall onto the stellar surface near the dipole pole. However, this paradigm assumes the presence of a single star, and a complete description of the accretion process in multiple systems remains to be achieved. Building on our previous work on DQ Tau and AK Sco, we aim to describe the accretion processes in cTTS binaries, accounting for the influence of stellar magnetic fields. Specifically, we sought to explore how the magnetospheric accretion model of cTTSs can be applied to V4046 Sgr, a spectroscopic binary composed of equal-mass and coeval cTTSs in a circular orbit with synchronous rotation, surrounded by a circumbinary disc. We analysed a time series of ESPaDOnS spectra covering several orbital cycles. A variability analysis was performed on the radial velocities and on the Balmer, He I D3, and Ca II emission lines, which are associated with the accretion process. We identified the secondary as the system's main accretor, operating in an unstable regime. Additionally, we detected an accretion funnel flow connecting the dipole pole of the primary star with a nearby bulk of gas. We concluded that the two components exhibit dissimilar accretion patterns. The primary operates in an "ordered chaotic" regime, where accretion funnel flows and accretion tongues coexist. Conversely, the secondary appears to be in a chaotic regime, with accretion tongues dominating.

This work continues the analysis of the model for calculating the thermal structure of an axisymmetric protoplanetary disk, initiated in the paper by Pavlyuchenkov (2024). The model is based on the well-known Flux-Limited Diffusion (FLD) approximation with separate calculation of heating by direct stellar radiation (hereinafter referred to as the FLD$^{\rm s}$ method). In addition to the previously described FLD$^{\rm s}$ model with wavelength-averaged opacities, we present a multiband model mFLD$^{\rm s}$, where the spectrum of thermal radiation is divided into several frequency bands. The model is based on an implicit finite-difference scheme for the equations of thermal radiation diffusion, which reduces to a system of linear algebraic equations written in hypermatrix form. A modified Gauss method for inverting the sparse hypermatrix of the original system of linear equations is proposed. The simulation results described in the article show that the midplane radial temperature profile obtained with the mFLD$^{\rm s}$ method has a variable slope in accordance with the reference Monte Carlo radiative transfer simulations. The mFLD$^{\rm s}$ model also qualitatively reproduces the non-isothermality of the temperature distribution along the angular coordinate near the midplane, which is not provided by the FLD$^{\rm s}$ method. However, quantitative differences remain between the reference temperature values and the results of mFLD$^{\rm s}$. These differences are likely due to the diffusive nature of the FLD approximation. It is also shown that the characteristic times for the disk to reach thermal equilibrium within the mFLD$^{\rm s}$ model can be significantly shorter than in FLD$^{\rm s}$. This property should be taken into account when modeling non-stationary processes in protoplanetary disks within FLD-based models.

The abundance of galaxy clusters as a function of mass and redshift is a well-established and powerful cosmological probe. Cosmological analyses based on galaxy cluster number counts have traditionally relied on explicitly computed likelihoods, which are often challenging to develop with the required accuracy and expensive to evaluate. In this work, we implement an alternative approach based on simulation-based inference (SBI) methods that relies solely on synthetic galaxy cluster catalogues generated under a given model. These catalogues are much easier to produce than it is to develop and validate a likelihood. We validate this approach in the context of the galaxy cluster survey of the upcoming Simons Observatory for a setup in which we can also evaluate an exact explicit likelihood. We find that our SBI-based approach yields cosmological parameter posterior means that are within $0.2\,\sigma$ of those obtained with the explicit likelihood and with biases smaller than $0.1\,\sigma$. We also introduce and validate a procedure to assess the goodness of fit using only synthetic catalogues similar to those used for training. This demonstrates, for the first time, that a galaxy cluster number count cosmological analysis can be performed fully without resorting to a likelihood at any stage. Finally, we apply our SBI-based approach to the real Planck MMF3 cosmology sample, obtaining cosmological parameter constraints that are within $0.1\,\sigma$ of their likelihood-based counterparts. This constitutes the first SBI-based number count cosmological analysis of a real galaxy cluster catalogue.

This study presents highly precise measurements of the cross-correlation between volume-limited galaxy samples from the DESI legacy survey catalogue and X-ray selected galaxy clusters from eROSITA, allowing for detailed analysis across redshift and color. Two key findings emerge. First, the cluster-galaxy cross-correlation, when split into quiescent and star-forming galaxies, contains significant information about the infall, feedback, and quenching processes of blue cloud galaxies in massive environments. These results align well with existing galaxy evolution models for higher stellar masses ($\log_{10}(M^*[M_\odot]) > 10.75$), though the red fraction may be slightly underestimated in the intermediate mass range ($10.25 < \log_{10}(M^*[M_\odot])< 10.75$). Second, the integral of the cross-correlation within 500 kpc enables a model-independent measurement of the red sequence and its scatter in clusters, providing a robust alternative to existing red-sequence calibration methods without requiring spectroscopic redshifts or classifications of galaxies. Similar analyses on upcoming photometric surveys as Euclid and LSST together with spectroscopic samples like 4MOST and DESI should lead to a significant increase in the signal-to-noise ratio and in particular at small separations.

Modeling the internal structure of self-gravitating solid and liquid bodies presents a challenge, as existing approaches are often limited to either overly simplistic constant-density approximations or more complex numerical equations of state. We present a detailed analysis of a tractable and physically motivated model for perfectly elastic, spherically symmetric self-gravitating bodies in hydrostatic equilibrium. The model employs a logarithmic equation of state (logotropic EOS) with a non-zero initial density and constant bulk modulus. Importantly, scaling properties of the model allow all solutions to be derived from a single, universal solution of an ordinary differential equation, resembling the Lane-Emden and Chandrasekhar models. The model provides new insights into stability issues and reveals oscillatory asymptotic behavior in the mass-radius relation, including the existence of both a maximum mass and a maximum radius. We derive useful, simple analytical approximations for key properties, such as central overdensity, moment of inertia, binding energy, and gravitational potential, applicable to small, metallic bodies like asteroids and moons. These new approximations could aid future research, including space mining and the scientific characterization of small Solar System bodies.

This study explores the detection of Quasi-Periodic Oscillations (QPOs) in blazars as a method to identify kink events within their jets, utilizing both $\gamma$-ray and polarized light observations. Focusing on a sample of 9 blazars, we analyze $\gamma$-ray light curves to identify significant QPOs. In addition to $\gamma$-ray data, we incorporated polarized light data corresponding to the same temporal segments to cross-validate the presence of QPOs. However, the limited availability of comprehensive polarized data restricted our ability to perform a thorough analysis across all datasets. Despite these limitations, our analysis reveals a segment where QPOs in polarized light coincided with those observed in $\gamma$-rays, providing preliminary evidence supporting the kink origin of these oscillations.

The Pierre Auger Observatory advances the study of ultra-high-energy cosmic rays through a hybrid system of surface and fluorescence detectors. This paper presents recent results, including refined spectrum measurements, anisotropy evidence, and new insights into cosmic-ray composition. Studies at energies beyond terrestrial accelerators reveal implications for particle physics. The AugerPrime upgrade will further enhance particle identification and extend the sensitivity to photons and neutrinos, broadening the Observatory's capability to explore cosmic-ray sources and propagation, paving the way for new discoveries.

Black holes can acquire magnetic flux from their magnetized progenitor or via prolonged accretion. We study the evolution of black hole magnetospheres by means of axisymmetric general relativistic magnetohydrodynamic simulations. We show that all simulated initial magnetic field geometries of varying complexity ultimately evolve into a split monopole magnetosphere. The magnetospheric evolution consists of two phases. In the first phase, the magnetosphere evolves toward pressure equilibrium accompanied by a large magnetic flux decrease on the event horizon on a fast Alfvénic timescale of $\sim 60$ light-crossing times of the gravitational radius. The second phase proceeds in a pressure balance in which the magnetic flux decays and current sheets shift in polar angle over the event horizon on slower resistive timescales. We present an analytic model for the second phase. Furthermore, we show that in a split monopole magnetosphere the magnetic flux on the event horizon decays exponentially with a timescale that depends on the black hole spin, where higher spin results in slower decay. Our results can have an implications for the timescales of reconnection-powered flares and for multimessenger counterparts to gravitational wave events.

Cosmogenic neutrinos (CNs) are produced by ultra-high energy cosmic rays (UHECRs) interacting with cosmic background radiation. We investigated the properties of CN point/extended sources, i.e, the neutrino spectrum, and angular profile as functions of time, by assuming that UHECR sources are transient events, such as gamma-ray bursts. The properties depend much on the intergalactic magnetic field (IGMF), but the angular extent is in general sub-degree, within which the CN flux can overshoot the diffuse CN flux in early time. The nearby CN point sources could be detected for the low IGMF case by future neutrino telescopes. The recent KM3-230213A event is possible to account for by a nearby transient CN source, rather than diffuse CN emission. Observations of CN point sources will provide a chance to search for UHECR sources.

Subhaloes are critical in distinguishing dark matter models, yet their evolution within galactic haloes, particularly in the Fuzzy Dark Matter (FDM) model, remains challenging to fully investigate in numerical simulations. In this work, we employ the fluid-wave hybrid scheme recently implemented in the GAMER-2 code to perform a cosmological zoom-in simulation of a Milky Way-sized halo with an FDM particle mass of m = 2 x 10^(-23) eV. It simultaneously resolves the solitonic core of the host halo and tracks the complex tidal evolution of subhaloes down to redshift z = 0. We examine the internal structure of subhaloes by analyzing their density profiles, velocity dispersions, and density power spectra across various redshifts. Our findings show that partially tidally stripped subhaloes deviate from the core-halo mass relation; their solitons remain intact and are enveloped by smaller granules predominantly from the host halo. Furthermore, our simulation unravels a complex tidal evolution of FDM subhaloes. On the one hand, we observe a subhalo core undergoing complete tidal disruption at z ~ 0.14, which later reemerges near the outskirts of the host halo around z ~ 0. This disruption event, characterized by a core contaminated with interference fringes from the host halo's wave function, occurs earlier than previously predicted. On the other hand, FDM subhaloes have denser cores before infall due to the presence of central solitons, making them more resilient to tidal disruption than their N-body counterparts. Our results demonstrate GAMER-2's capability to resolve non-linear FDM substructure down to z = 0, paving the way for future studies of larger FDM subhalo samples with heavier particle masses.

Context. One of the primary goals of Galactic Archaeology is to reconstruct the Milky Way's accretion history. To achieve this, significant efforts have been dedicated to identifying signatures of past accretion events. In particular, the study of integrals-of-motion (IoM) space has proven to be highly insightful for uncovering these ancient mergers and understanding their impact on the Galaxy's evolution. Aims. This paper evaluates the effectiveness of a state-of-the-art method for detecting debris from accreted galaxies, by testing it on four Milky Way-like galaxies from the Auriga suite of cosmological magneto-hydrodynamical simulations. Methods. We employ the innovative method from Lövdal et al. (2022) to identify substructures in the integrals-of-motion space within the local stellar halos of the four simulated galaxies. This approach enables us to evaluate the method's performance by comparing the properties of the identified clusters with the known populations of accreted galaxies in the simulations. Additionally, we investigate whether incorporating chemical abundances and stellar age information can help to link distinct structures originating from the same accretion event. Results. This method is very effective in detecting debris from accretion events that occur less than 6-7 Gyr ago but struggles to detect most of the debris from older accretion. Furthermore, most of the detected structures suffer from significant contamination by in-situ stars. Our results also show that the method may also generate artificial detections. Conclusions. Our work show that the Milky Way's accretion history remains uncertain, and question the reality of some detected structures in the Solar vicinity.

In the presence of spacetime torsion, any generic $f(R)$ model of gravity is conformally dual to a scalar-tensor theory augmented with a second rank antisymmetric massless degree of freedom. We investigate the stochastic gravitational wave background (SGWB) that may be sourced directly at the second order by such a torsional field, treated perturbatively during an epoch of canonical, single-field, slow-roll inflation. The resulting second-order induced SGWB, which dominates over the primary inflationary GW background at all scales, peaks only at ultra-low frequencies, and is found to be extremely red-tilted with an effective tensor spectral index $\alpha_{\rm T}\sim-6$ on matter-dominated scales. The signal is potentially within the reach of upcoming indirect GW probes on very large scales $k\lesssim10^{-2}\:\textrm{Mpc}^{-1}$, i.e., next-generation CMB experiments like the LiteBIRD. In the near future, observation of such a markedly red-tilted SGWB on CMB scales could hence provide a novel and unique clue in favour of torsional gravity during the inflationary era.

We present an analysis of DESI Data Release 1 (DR1) that incorporates Halo Occupation Distribution (HOD)-informed priors into Full-Shape (FS) modeling of the power spectrum based on cosmological perturbation theory (PT). By leveraging physical insights from the galaxy-halo connection, these HOD-informed priors on nuisance parameters substantially mitigate projection effects in extended cosmological models that allow for dynamical dark energy. The resulting credible intervals now encompass the posterior maximum from the baseline analysis using gaussian priors, eliminating a significant posterior shift observed in baseline studies. In the $\Lambda$CDM framework, a combined DESI DR1 FS information and constraints from the DESI DR1 baryon acoustic oscillations (BAO)-including Big Bang Nucleosynthesis (BBN) constraints and a weak prior on the scalar spectral index-yields $\Omega_{\rm m} = 0.2994\pm 0.0090$ and $\sigma_8 = 0.836^{+0.024}_{-0.027}$, representing improvements of approximately 4% and 23% over the baseline analysis, respectively. For the $w_0w_a$CDM model, our results from various data combinations are highly consistent, with all configurations converging to a region with $w_0 > -1$ and $w_a < 0$. This convergence not only suggests intriguing hints of dynamical dark energy but also underscores the robustness of our HOD-informed prior approach in delivering reliable cosmological constraints.

We develop a novel approach to constrain the Hubble parameter $H_0$ and the primordial power spectrum amplitude $A_\mathrm{s}$ using supernovae type Ia (SNIa) data. By considering SNIa as tracers of the peculiar velocity field, we can model their distance and their covariance as a function of cosmological parameters without the need of calibrators like Cepheids; this yields a new independent probe of the large-scale structure based on SNIa data without distance anchors. Crucially, we implement a differentiable pipeline in JAX, including efficient emulators and affine sampling, reducing inference time from years to hours on a single GPU. We first validate our method on mock datasets, demonstrating that we can constrain $H_0$ and $\log 10^{10}A_\mathrm{s}$ within $\sim10\%$ using $\sim10^3$ SNIa. We then test our pipeline with SNIa from an $N$-body simulation, obtaining $7\%$-level unbiased constraints on $H_0$ with a moderate noise level. We finally apply our method to Pantheon+ data, constraining $H_0$ at the $10\%$ level without Cepheids when fixing $A_\mathrm{s}$ to its $\it{Planck}$ value. On the other hand, we obtain $15\%$-level constraints on $\log 10^{10}A_\mathrm{s}$ in agreement with $\it{Planck}$ when including Cepheids in the analysis. In light of upcoming observations of low redshift SNIa from the Zwicky Transient Facility and the Vera Rubin Legacy Survey of Space and Time, surveys for which our method will develop its full potential, we make our code publicly available.

We investigate whether a violation of the distance duality relation (DDR), $D_L(z) = (1+z)^2 D_A(z)$, connecting the angular diameter and luminosity distances, can explain the Hubble tension and alter the evidence for dynamical dark energy in recent cosmological observations. We constrain five phenomenological parameterisations of DDR violation using Baryon Acoustic Oscillation measurements from the DESI survey calibrated with the sound horizon derived from \textit{Planck} Cosmic Microwave Background data and the Pantheon+ Type Ia supernova (SNIa) catalogue calibrated with the supernova absolute magnitude from S$H_0$ES. We find that two toy models can resolve the tension: a constant offset in the DDR (equivalent to a shift in the calibration of the SNIa data), $D_L(z)/D_A(z)\simeq 0.925(1+z)^2$, which leaves the hint for evolving dark energy unaffected; or a change in the power-law redshift-dependence of the DDR, restricted to $z\lesssim 1$, $D_L(z)/D_A(z)\simeq(1+z)^{1.866}$, together with a {\it constant} phantom dark energy equation of state $w\sim -1.155$. The Bayesian evidence slightly favours the latter model. Our phenomenological approach motivates the investigation of physical models of DDR violation as a novel way to explain the Hubble tension.

In this work, we develop a generic formalism for the study of tensor perturbations induced at second order by first-order vector metric perturbations, dubbing these induced tensor modes $\textit{vector-induced gravitational waves}$ (VIGWs). Notably, considering an inflation-inspired power-law type magnetic field power spectrum of the form $P_B(k)\propto k^{n_\mathrm{B}}$ (where $n_{\rm B}$ is the magnetic spectral index), we show that the VIGW signal is enhanced for stiff post-inflationary EoS, with the maximum enhancement happening for $w=1$. We explicitly demonstrate this contribution is dominant over the first-order magnetically-sourced GWs. The VIGW spectrum exhibits a maximum at around the scale crossing the cosmological horizon at the end of reheating, $k_\mathrm{reh}$, with its present day peak amplitude scaling as $\Omega_{\rm GW}(k_{\rm reh},\eta_0)\propto \Delta N_{\rm reh}\times(H_{\rm inf}/M_{\rm Pl})^{8}$, where $H_{\rm inf}$ is the Hubble parameter at the end of inflation and $\Delta N_{\rm reh}$ the duration of the post-inflationary era in $e$-folds. For $w=1$ (kination) and $n_{\rm B}>-3/2$, one further obtains a nearly $n_{\rm B}$-independent frequency scaling of the GW spectrum of the form $\Omega_{\rm GW}(f,\eta_0)\propto \left(\frac{f}{f_{\rm reh}}\right)^{-2.8}$ for $f>f_\mathrm{reh}\equiv k_\mathrm{reh}/(2\pi)$. Finally, we need to highlight that the VIGW signal can be well within the detection bands of several next-generation interferometric GW missions at small scales. Indicatively, for $H_{\rm inf} \sim O(10^{7})\:\mathrm{GeV}$ and $O(10^{14})\:\mathrm{GeV}$, and $\Delta N_{\rm reh} \sim 15$ and $10$, the VIGW signal is found to be detectable by LISA and ET respectively.

We study the angular power spectrum of gravitational-wave and galaxy catalogs in tomographic redshift and distance bins as a probe of late-time cosmology, focusing specifically on next-generation ground-based interferometers in combination with the Euclid photometric survey. We assess the potential of this technique to constrain the Hubble constant and the matter energy density. Our analysis incorporates realistic gravitational-wave source populations, error modelling calibrated on recent detector designs, and accounts for nuisance parameters. We show that the tomographic angular cross-correlation could determine the Hubble constant to percent or sub-percent precision depending on the binning choice, configuration and operation time of gravitational-wave observatories. This conclusion holds even when marginalising over the unknown tracer biases, primordial power-spectrum parameters and baryon density. In particular, we show that the combination of the galaxy auto-correlation spectra and the cross-correlation of gravitational waves and galaxy surveys can lead to an improvement of up to a factor ${\sim}10$ in constraining power over either of the two probes taken individually. However, this prospect crucially relies on the presence of multiple gravitational-wave interferometers able to yield precise sky localisation. We also discuss the use of a spectroscopic redshift catalog, as well as the detectability of the clustering bias of gravitational-wave sources.

We revisit scale separation for compactifications of ten- and eleven-dimensional supergravity. For cosmological solutions rolling down flux-generated potentials, we observe that scale separation is achieved as time flows, and is fairly generic. This is realized without the need of orientifolds nor corrections to the classical supergravity approximation. We then confront scale separation with the Covariant Entropy Bound (CEB) and the CKN bound. We show that a naive application of these bounds to vacua hints at the existence of at least two extra dimensions. For rolling solutions, we observe that the CEB is not always respected, but since these examples lack a cosmic horizon, the application of entropy bounds remains delicate.

Ehlers' Frame Theory is a class of geometric theories parameterized by $\lambda := 1/c^2$ and identical to the General Theory of Relativity for $\lambda \neq 0$. The limit $\lambda \to 0$ does not recover Newtonian gravity, as one might expect, but yields the so-called Newton-Cartan theory of gravity, which is characterized by a second gravitational field $\boldsymbol{\omega}$, called the Coriolis field. Such a field encodes at a non-relativistic level the dragging feature of general spacetimes, as we show explicitly for the case of the $(\eta,H)$ geometries. Taking advantage of the Coriolis field, we apply Ehlers' theory to an axially symmetric distribution of matter, mimicking, for example, a disc galaxy, and show how its dynamics might reproduce a flattish rotation curve. In the same setting, we further exploit the formal simplicity of Ehlers' formalism in addressing non-stationary cases, which are remarkably difficult to be treated in the General Theory of Relativity. We show that the time derivative of the Coriolis field gives rise to a tangential acceleration which allows to study a possible formation in time of the rotation curve's flattish feature.

Paleodetection has been proposed as a competitive method for detecting dark matter and other new physics interactions, complementing conventional direct detection experiments. In this work, we utilise TRIM simulations to improve the modelling of track length distributions. Our findings suggest that previous studies have overestimated the number of tracks caused by weakly interacting particles, and that the lowest observable dark matter mass should be higher than previously predicted. These differences are mainly attributed to the fact that (a) the recoil energy-track length relation is not one-to-one, (b) at low recoil energies, a substantial fraction of recoils do not yield any tracks, and (c) at high energies, electronic stopping becomes dominant, resulting in a track length barrier at $\sim200$ nm. In addition to WIMPs, we also modelled tracks from generalised coherent elastic neutrino nucleus scattering (CE$\nu$NS) via new light mediators and estimated the projected sensitivity for these interactions.

We investigate the cosmological dynamics of interacting dark energy within the framework of $\alpha$-attractor models. Specifically, we analyze the associated autonomous system, focusing on its fixed points that represent dark energy and scaling solutions, along with their stability conditions. We employ center manifold theory to address cases where some fixed points display eigenvalues with zero and negative real parts. The model reveals attractors describing dark energy, enabling a smooth transition from the radiation-dominated era to the matter-dominated era, and ultimately into the dark-energy-dominated phase. Additionally, we identify a scaling matter solution capable of modifying the growth rate of matter perturbations during the matter-dominated epoch. Consequently, we study the evolution of matter perturbations by obtaining both analytical and numerical solutions to the density contrast evolution equation. Based on these results, we compute numerical solutions for the weighted growth rate $f\sigma_{8}$, indicating that interacting $\alpha$-attractor dark energy models may provide a better fit to structure formation data than the standard $\Lambda$CDM scenario.

The spontaneous breaking of $SO(10)$ grand unified symmetry to $SU(3)_c \times SU(2)_L \times U(1)_Y \times U(1)_\chi$ yields the GUT monopole as well as a comparably heavy $U(1)_\chi$ monopole which also carries $U(1)_Y$ flux. A metastable string scenario in this case requires that the $U(1)_\chi$ symmetry is necessarily broken close to the GUT scale, thus resulting in a dimensionless string tension $G \mu \sim 10^{-6}$. We show that the $\chi$ monopole does not carry any unconfined flux following the electroweak symmetry breaking. Coupled with $G \mu \sim 10^{-6}$, this metastable string network appears to provide a good fit to the recent Pulsar Timing Array data on the stochastic gravitational background. Gauge coupling unification, especially in the presence of low scale supersymmetry, determines the GUT scale and, in combination with constraints from proton decay experiments, one is able to constrain some of the key parameters in this setup. The breaking of $SO(10)$ via $SU(5) \times U(1)_\chi$ also yields superheavy metastable strings with no unconfined flux associated with the monopoles. Finally, we consider $SO(10)$ breaking via $SU(4)_c \times SU(2)_L \times U(1)_R$, $SU(3)_c \times SU(2)_L \times SU(2)_R \times U(1)_{B-L}$ and flipped $SU(5)$ that yield metastable strings where the associated monopoles carry unconfined flux after the electroweak breaking.

The AION project has built a tabletop prototype of a single-photon long-baseline atom interferometer using the 87Sr clock transition - a type of quantum sensor designed to search for dark matter and gravitational waves. Our prototype detector operates at the Standard Quantum Limit (SQL), producing a signal with no unexpected noise beyond atom shot noise. Importantly, the detector remains at the SQL even when additional laser phase noise is introduced, emulating conditions in a long-baseline detector such as AION or AEDGE where significant laser phase deviations will accumulate during long atom interrogation times. Our results mark a key milestone in extending atom interferometers to long baselines. Such interferometers can complement laser-interferometer gravitational wave detectors by accessing the mid-frequency gravitational wave band around 1 Hz, and can search for physics beyond the Standard Model.

We propose a pseudo-scalar quantity, which is an analogue of the Chern-Simons invariant, in the framework of non-metricity gravity. By considering the coupling between the pseudo-scalar quantity and the axion, we give scenarios which may solve the problems of the axion misalignment, the $S_8$ problem, and the beginning of inflation. When the phase transition associated with the spontaneous breaking of the gauge symmetry of the electroweak theory or grand unified theories (GUTs) occurs, the pseudo-scalar quantity has a non-trivial value, which induces the misalignment of the axion field and axion particles are produced. If the gradient of the potential is small, the $S_8$ problem might be solved. We also propose a mechanism which induces inflation by the misalignment of the axion field generated by the phase transition of the GUTs.

We present a very detailed derivation of solutions describing hairy black holes within the gravity-coupled Weinberg-Salam theory, which were previously reported in \href{this https URL}{this http URL. 133 (2024) 171402}. These black holes support a strong magnetic field that polarizes the electroweak vacuum and creates a condensate of massive fields carrying superconducting currents along the black hole horizon. The currents, in turn, generate a ``corona'' of magnetic vortex segments attached to the horizon at both ends. The condensate and corona together constitute the black hole hair. The extremal solutions approach, in the far field, the magnetic Reissner-Nordström configuration, with a total mass that is {\it lower} than the total charge, $M<|Q|$, due to the negative Zeeman energy of the condensate. This makes the removal of the hair energetically unfavorable. The maximally hairy black holes exhibit masses comparable to terrestrial values, with approximately 11\% of their total mass stored in the hair. Given that these solutions arise within a well-tested theoretical framework, they are likely to have physical relevance.

Radio Frequency Interference (RFI) from anthropogenic radio sources poses significant challenges to current and future radio telescopes. Contemporary approaches to detecting RFI treat the task as a semantic segmentation problem on radio telescope spectrograms. Typically, complex heuristic algorithms handle this task of `flagging' in combination with manual labeling (in the most difficult cases). While recent machine-learning approaches have demonstrated high accuracy, they often fail to meet the stringent operational requirements of modern radio observatories. Owing to their inherently time-varying nature, spiking neural networks (SNNs) are a promising alternative method to RFI-detection by utilizing the time-varying nature of the spectrographic source data. In this work, we apply Liquid State Machines (LSMs), a class of spiking neural networks, to RFI-detection. We employ second-order Leaky Integrate-and-Fire (LiF) neurons, marking the first use of this architecture and neuron type for RFI-detection. We test three encoding methods and three increasingly complex readout layers, including a transformer decoder head, providing a hybrid of SNN and ANN techniques. Our methods extend LSMs beyond conventional classification tasks to fine-grained spatio-temporal segmentation. We train LSMs on simulated data derived from the Hyrogen Epoch of Reionization Array (HERA), a known benchmark for RFI-detection. Our model achieves a per-pixel accuracy of 98% and an F1-score of 0.743, demonstrating competitive performance on this highly challenging task. This work expands the sophistication of SNN techniques and architectures applied to RFI-detection, and highlights the effectiveness of LSMs in handling fine-grained, complex, spatio-temporal signal-processing tasks.

We propose a novel modification to the optical benches of space-based gravitational wave detectors (SGWDs) to enable the detection of axion-like dark matter (ALDM)-induced birefringence without altering the polarization of inter-spacecraft laser links. Our design introduces an auxiliary interferometer to convert polarization modulation into measurable phase shifts. Analytical expressions for sensitivity to the ALDM-photon coupling are derived for various time-delay interferometry (TDI) combinations. Projected sensitivity curves demonstrate complementary coverage across the ALDM mass range $10^{-19}\sim10^{-14}\mathrm{eV}$. This approach preserves the original interferometric stability while enabling new physics capabilities for SGWDs.

Highly eccentric binary neutron star mergers exhibit unique dynamical and observational signatures compared to quasi-circular ones in terms of their gravitational wave signal and the ejection of matter, leading to different electromagnetic counterparts. In this article, we present general relativistic magneto-hydrodynamic simulations of binary neutron star systems on highly eccentric orbits. While in quasi-circular binaries, the influence of the magnetic field is too weak to affect the general pre-merger dynamics, the close encounters in eccentric systems could potentially trigger magneto-hydrodynamic instabilities. Therefore, we investigate possible effects before, during, and after the merger for a total of three different systems with varying initial eccentricity. We study the f-mode oscillations excited by tidal interaction in close encounters and find good agreement with predicted f-mode frequency estimates. However, our simulations reveal no significant differences compared to results neglecting the magnetic field. Although we observe a rearrangement of the poloidal structure of the magnetic field inside the stars, there is no relevant increase in the magnetic energy during the encounters. Also, during the merger, the amplification of the magnetic field seems to be largely independent of the eccentricity in our systems. Consistent with studies of merging non-magnetized binary neutron stars, we find a correlation between eccentricity and mass ejection, with a higher impact parameter leading to a larger amount of unbound material.

We put into test the idea of replacing dark energy by a vector field against the cosmic microwave background (CMB) observation using the simplest vector-tensor theory, where a massive vector field couples to the Ricci scalar and the Ricci tensor quadratically. First, a remarkable Friedmann-Lemaître-Robertson-Walker (FLRW) metric solution that is completely independent of the matter-energy compositions of the universe is found. Second, based on the FLRW solution as well as the perturbation equations, a numerical code calculating the CMB temperature power spectrum is built. We find that though the FLRW solution can mimic the evolution of the universe in the standard $\Lambda$CDM model, the calculated CMB temperature power spectrum shows unavoidable discrepancies from the CMB power spectrum measurements.

We derive analytic dispersion relations for cold, orbitally constrained systems governed by the Vlasov equation. For magnetized plasmas, we obtain the first explicit relation for two-dimensional anisotropic BGK modes with finite magnetic field, showing that only a finite number of angular modes can become unstable and identifying a magnetic-field threshold for stabilization. In the gravitational case, we establish a bound on the growth rate of core perturbations, set by the potential's curvature. These results clarify how orbital constraints shape the spectrum and growth of kinetic instabilities in cold, collisionless media.

We compute a holographic entanglement entropy via Ryu--Takayanagi prescription in the three-dimensional Friedmann--Lemaître--Robertson--Walker universe. We consider two types of holographic scenarios analogous to the static patch holography and the half de Sitter holography, in which the holographic boundary is timelike and placed in the bulk. We find in general that the strong subadditivity can be satisfied only in the former type and in addition the holographic boundary has to fit inside the apparent horizon. Also, for the universe filled with an ideal fluid of constant equation of state $w<-1$, the condition is sharpened as that the holographic boundary has to fit inside the event horizon instead. These conditions provide a necessary condition for the dual quantum field theory to be standard and compatible with the strong subadditivity.

We study quantum decoherence of curvature perturbations at superhorizon scales caused by the gravitational nonlinearities. We show that cubic gravitational couplings, constrained by the spatial diffeomorphism invariance, lead to infrared (IR) and ultraviolet (UV) divergences in the decoherence rate at one loop. These divergences arise from fluctuations of deep IR modes which look like a background mode for a local observer and violent zero-point fluctuations in the deep UV, respectively. We argue that these divergences are unobservable, as they vanish when considering proper observables. We consider correlators defined using the geodesic distance for IR divergences and time-averaged correlators for UV divergences. To account for these observer's perspectives, we propose to consider an effective quantum state, defined in terms of actual observables, as a more appropriate probe of the quantum coherence of the system measured by an observer. We then evaluate the finite decoherence rate induced by superhorizon environment during inflation and at late universe.

The Atacama Cosmology Telescope (ATC) has recently released new measurements and constraints on inflationary observables. In this paper it is shown that a component of a dynamical affine connection, which is independent of the metric, can easily drive inflation in agreement with these observations. Such geometrical explanation of inflation is analysed in detail here in the minimal model, including the predictions for the scalar spectral index $n_s$ and its running $\alpha_s$, the amplitude of the scalar perturbations and the tensor-to-scalar ratio $r$. Furthermore, this minimal model is shown to provide an inflationary attractor: arbitrary initial values of the kinetic energy density are dynamically attracted down to negligible values compared to the potential energy density in homogeneous and isotropic metrics. The role of the Higgs boson during and after inflation is also briefly discussed.