Papers for Monday, Dec 16 2024

Vertical stellar kinematics+density can be used to trace the dark matter distribution [or the equivalent phantom mass in a Modified Newtonian Dynamics (MOND) scenario] through Jeans equations. In this paper, we want to improve this type of analysis by making use of the recent data of the 6D information from the Gaia-DR3 survey in the anticenter and the Galactic poles to obtain the dynamical mass distribution near plane regions, including extended kinematics over a wide region of 8 kpc$<R<$22 kpc, $|z|<3$ kpc. Our conclusions are as follows: (i) the model of the spherical dark matter halos and the MOND model are compatible with the data; (ii) the model of the disky dark matter (with density proportional to the gas density) is excluded; (iii) the total lack of dark matter (there is only visible matter) within Newtonian gravity is compatible with the data; for instance, at solar Galactocentric radius, we obtained $\Sigma=39\pm 18$ M$_\odot$ pc$^{-2}$ for $z=1.05$ kpc, compatible with the expected value for visible matter alone of 44 M$_\odot$ pc$^{-2}$, thus allowing zero dark matter. Similarly, for $R>R_\odot$, $z=1.05$ kpc: $\Sigma=28.7\pm 9.6$, $23.0\pm 5.7$, $16.9\pm 5.8$, $11.4\pm 6.6$ M$_\odot$ pc$^{-2}$, respectively, for $R=10,13,16,19$ kpc, compatible with visible matter alone. Larger error bars in comparison with previous works are not due to worse data or a more awkward technique but to a stricter modeling of the stellar distribution.

We present a comprehensive analysis of a planetary microlensing event OGLE-2015-BLG-1609. The planetary anomaly was detected by two survey telescopes, OGLE and MOA. Each of these surveys collected enough data over the planetary anomaly to allow for an unambiguous planet detection. Such survey detections of planetary anomalies are needed to build a robust sample of planets that could improve studies on the microlensing planetary occurrence rate by reducing biases and statistical uncertainties. In this work, we examined different methods for modeling microlensing events using individual datasets, particularly we incorporated a Galactic model prior to better constrain poorly defined microlensing parallax. Ultimately, we fitted a comprehensive model to all available data, identifying three potential typologies, with two showing comparably high Bayesian evidence. Our analysis indicates that the host of the planet is a brown dwarf with a probability of 34%, or a low-mass stellar object (M-dwarf) with the probability of 66%.

We use high-resolution cosmological zoom-in simulations to model feedback from Seyfert-type supermassive black hole (SMBH) jets onto galaxies with identical dark matter (DM) halos of log(M/M$_\odot$) ~ 11.8. The low mass, ~10$^6$ M$_\odot$, seed SMBHs, have been introduced when the parent DM halos have reached log(M/M$_\odot$) ~ 11, at z ~ 3.7. In a controlled experiment, we vary only the efficiency of the SMBH accretion and focus on galaxies and their immediate environment properties. Our results show that the AGN jet feedback has a substantial effect on the basic properties of Seyfert-type galaxies, such as morphology, gas fraction and distribution, star formation rate and distribution, bulge-to-disk ratio, DM halo baryon fraction, and properties of circumgalactic medium (CGM) and beyond. These have been compared to a galaxy with supernovae only feedback. We focus on the energy deposition by the jet in the ISM and IGM, and follow the expansion of the multiple jet cocoons to 2 Mpc. We find that the jet-ISM interaction gradually pushes the star formation to larger radii with increasing accretion efficiency, which results in increased mass of the outer stellar disk, which is best fit as a double-exponential disk. Furthermore, we compare our galaxies and their properties with the observed nearby Seyfert galaxies, including the scaling relations, and find a close agreement, although statistical analysis of observed Seyferts is currently missing. In a forthcoming paper, we focus on evolution of these objects at z<10 and study the effect of the SMBH seeding redshift on galaxy evolution.

With the advent of time-domain astronomy and the game-changing next generation of telescopes, we have unprecedented opportunities to explore the most energetic events in our Universe through electromagnetic radiation, gravitational waves, and neutrinos. These are elementary particles, which exist in three different flavors and change the latter as they propagate in the dense core of astrophysical sources as well as en route to Earth. To capitalize on existing and upcoming multi-messenger opportunities, it is crucial to understand: 1. the role of neutrinos in explosive transient sources as well as in the synthesis of the elements heavier than iron; 2. the impact of neutrino physics on the multi-messenger observables; 3. the information on the source physics carried by the detectable neutrino signal. In this review, the status of this exciting and fast-moving field is outlined, focusing on astrophysical sources linked to collapsing massive stars and neutron-star mergers. In light of the upcoming plethora of multi-messenger data, outstanding open issues concerning the optimization of multi-messenger detection strategies are discussed.

A better understanding of the ice-ocean couplings is required to better characterise the hydrosphere of the icy moons. Using global numerical simulations in spherical geometry, we have investigated here the interplay between rotating convection and a melting boundary. To do so, we have implemented and validated a phase field formulation in the open-source code \texttt{MagIC}. We have conducted a parameter study varying the influence of rotation, the vigour of the convective forcing and the melting temperature. We have evidenced different regimes akin to those already found in previous monophasic models in which the mean axisymmetric ice crust transits from pole-ward thinning to equator-ward thinning with the increase of the rotational constraint on the flow. The derivation of a perturbative model of heat conduction in the ice layer enabled us to relate those mean topographic changes to the underlying latitudinal heat flux variations at the top of the ocean. The phase change has also been found to yield the formation of sizeable non-axisymmetric topography at the solid-liquid interface with a typical size close to that of the convective columns. We have shown that the typical evolution timescale of the interface increases linearly with the crest-to-trough amplitude and quadratically with the mean melt radius. In the case of the largest topographic changes, the convective flows become quasi locked in the topography due to the constructive coupling between convection and ice melting. The tentative extrapolation to the planetary regimes yields $\mathcal{O}(10^2-10^3)$ meters for the amplitude of non-axisymmetric topography at the base of the ice layer of Enceladus and $\mathcal{O}(10^3-10^4)$ meters for Titan.

The first generation of transiting planet searches in globular clusters yielded no detections, and in hindsight, only placed occurrence rate limits slightly higher than the measured occurrence rate in the higher-metallicity Galactic thick disk. To improve these limits, we present the first results of a new wide field search for transiting hot Jupiters in the globular cluster 47~Tucanae. We have observed 47~Tuc as part of the Multiband Imaging Survey for High-Alpha Planets (MISHAPS). Using 24 partial and full nights of observations taken with the Dark Energy Camera on the 4-m Blanco telescope at CTIO, we perform a search on 19,930 stars in the outer regions of the cluster. Though we find no clear planet detections, by combining our result with the upper limit enabled by Gilliland et al.'s 2000 Hubble search for planets around an independent sample of 34,091 stars in the inner cluster, we place the strongest limit to date on hot Jupiters with periods of $0.8 \leq P \leq 8.3$ days and $0.5~R_{\rm Jup} \leq R_{\rm P} \leq 2.0~R_{\rm Jup}$ of $f_{\rm HJ} < 0.11\%$, a factor of ${\sim}$4 below the occurrence rate in the \textit{Kepler} field. Our search found 35 transiting planet candidates, though we are ultimately able to rule out each without follow-up observations. We also found 4 eclipsing binaries, including 3 previously-uncataloged detached eclipsing binary stars.

Many classical Be stars acquire their very rapid rotation by mass and angular-momentum transfer in massive binaries. Short-lived intermediate-phase objects have only been discovered recently. Data archives and the literature have been searched for additional candidates exhibiting this patterns. Thirteen candidates were identified at various confidence levels. Adding to the two known systems identified as classical Be star+pre-subdwarf binaries (LB-1 and HR6819), two more (V742Cas, HD44637) could be confirmed with interferometry, with V742Cas setting a new record for the smallest visually observed angular semi-major axis, at a=0.663mas. Two further ones (V447Sct, V1362Cyg) are not resolved interferometrically, but other evidence puts them at the same confidence level as LB-1. V2174Cyg is a candidate with very high confidence, but was not observed interferometrically. The remaining ones are either candidates with varying levels of confidence. Of a mostly magnitude complete sample of 328 Be stars, 0.5-1% are found to have recently completed the mass overflow that led to their formation. Another 5% are systems with compact subdwarf companions, i.e., further evolved after a previous overflow, and possibly two more percent harbor white dwarfs. All these systems are of early B-subtypes, however, and if the original sample is restricted to early subtypes (136 objects), these percentages increase by a factor of about 2.5, while dropping to zero for the mid and late subtypes (together 204 objects). This strongly suggests that early- vs. mid- and late-type Be stars have differently weighted channels to acquire their rapid rotation, namely binary interaction vs. evolutionary spin-up.

The Data Management team of the Vera C. Rubin Observatory has developed a data description language and toolset, Felis, for defining the semantics and metadata of its public-facing data catalogs. Felis uses a rich Pydantic data model for describing and validating catalog metadata, expressed as a human-readable and editable YAML format. Felis also provides a Python library and command line interface for working with these data models. The metadata is used to populate the TAP_SCHEMA tables for the IVOA TAP services utilized by the Rubin Science Platform (RSP). Felis's current capabilities will be discussed along with some future plans.

Gamma-ray bursts (GRBs) are excellent probes of the high-redshift Universe due to their high luminosities and the relatively simple intrinsic spectra of their afterglows. They can be used to estimate the fraction of neutral hydrogen (i.e., neutral fraction) in the intergalactic medium at different redshifts through the examination of their Lyman-alpha damping wing with high quality optical-to-near-infrared spectra. Neutral fraction estimates can help trace the evolution of the Epoch of Reionization, a key era of cosmological history in which the intergalactic medium underwent a phase change from neutral to ionized. We revisit GRB 130606A, a z ~ 5.9 GRB for which multiple analyses, using the same damping wing model and data from different telescopes, found conflicting neutral fraction results. We identify the source of the discrepant results to be differences in assumptions for key damping wing model parameters and data range selections. We perform a new analysis implementing multiple GRB damping wing models and find a 3-sigma neutral fraction upper limit ranging from xHI < 0.20 to xHI < 0.23. We present this result in the context of other neutral fraction estimates and Epoch of Reionization models, discuss the impact of relying on individual GRB lines of sight, and highlight the need for more high-redshift GRBs to effectively constrain the progression of the Epoch of Reionization.

Spectroscopy and direct-imaging of ultra-faint targets such as Earth-like exoplanets and high redshift galaxies are among the primary goals of upcoming large scale astronomy projects like the Habitable World Observatory (HWO). Such objectives pose extreme instrumental challenges, in particular on detectors where dark currents lower than 1 e-/pixel/kilosecond and read noise less than 1 e-/pixel/frame will have to be achieved on large format arrays. Some technologies meet these requirements at optical wavelengths, but none do in the infrared. With this goal in mind, the University of Hawaii has partnered with Leonardo to develop linear-mode avalanche photodiodes (LmAPDs). In this paper, we report recent tests performed on LmAPDs, where we measure a ROIC glow of approximately 0.01 e-/pixel/frame, without which the intrinsic dark current is essentally zero (< 0.1 e- /pixel/kilosecond). We show that at high gain, these devices are capable of detecting single photons

We investigate the properties of active galactic nuclei (AGNs) in the brightest submillimeter galaxies (SMGs) in the COSMOS field. We utilize the bright sample of ALMA/SCUBA-2 COSMOS Survey (AS2COSMOS), which consists of 260 SMGs with $S_{\mathrm{870}\, \mu \mathrm{m}}=0.7\text{--}19.2\,\mathrm{mJy}$ at $z=0\text{--}6$. We perform optical to millimeter spectral energy distribution (SED) modeling for the whole sample. We identify 24 AGN-host galaxies from the SEDs. Supplemented by 23 X-ray detected AGNs (X-ray AGNs), we construct an overall sample of 40 AGN-host galaxies. The X-ray luminosity upper bounds indicate that the X-ray undetected SED-identified AGNs are likely to be nearly Compton thick or have unusually suppressed X-ray emission. From visual classification, we identify $25^{+6}_{-5}$\% of the SMGs without AGNs as major merger candidates. This fraction is almost consistent with the general galaxy population at $z\sim2$, suggesting that major mergers are not necessarily required for the enhanced star formation in SMGs. We also identify $47^{+16}_{-15}$\% of the AGN hosts as major merger candidates, which is about twice as high as that in the SMGs without AGNs. This suggests that major mergers play a key role in triggering AGN activity in bright SMGs.

Constraining the helium enhancement in stars is critical for understanding the formation mechanisms of multiple populations in star clusters. However, measuring helium variations for many stars within a cluster remains observationally challenging. We use Hubble Space Telescope photometry combined with MUSE spectroscopic data for over 7,200 red-giant branch stars in \omc\ to measure helium differences between distinct groups of stars as a function of metallicity separating the impact of helium enhancements from other abundance variations on the pseudo-color (chromosome) diagrams. Our results show that stars at all metallicities have subpopulations with significant helium enhancement ($\Delta Y_{min} \gtrsim$ 0.11). We find a rapid increase in helium enhancement from low metallicities ($\rm{[Fe/H] \simeq -2.05}$ to $\rm{[Fe/H] \simeq -1.92})$, with this enhancement leveling out at \deltay\ $= 0.154$ at higher metallicities. The fraction of helium-enhanced stars steadily increases with metallicity ranging from 10\% at $\rm{[Fe/H] \simeq -2.04}$ to over $90\%$ at $\rm{[Fe/H] \simeq -1.04}$. This study is the first to examine helium enhancement across the full range of metallicities in \omc{}, providing new insight into its formation history and additional constraints on enrichment mechanisms.

The gravitational time delays of macro-lenses can be used to constrain the rest mass of the photon with high accuracy. Assuming a point-mass $+$ external shear lens model, we prove that an upper limit of the photon mass can be derived directly from two observables--the time delay $\Delta t$ and the leading-to-trailing flux ratio $R$ of strongly lensed fast radio bursts (FRBs). Using the observed values of $\Delta t$ and $R$ of a lensed FRB candidate, i.e., FRB 20190308C, as a reference, we obtain a strict upper limit of the photon mass between $m_\gamma < 5.3 \times {10}^{-42}\,\rm kg$, for a given external shear strength of $\gamma' = 0.01$, and $m_{\gamma} < 2.1 \times 10^{-41}-2.4 \times 10^{-42}\,\text{kg}$, within the external shear range of $0<\gamma'<1$. This provides the most stringent limit to date on the photon mass through gravitational lensing time delays, improving by 1 to 2 orders of magnitude the previous results obtained from lensed active galactic nuclei.

We present a qualitative comparison between the host and black hole properties of radio galaxies in the MeerKAT GigaHertz Tiered Extragalactic Exploration~(MIGHTEE) survey with the radio galaxy population in the SIMBA suite of cosmological hydrodynamical simulations. The MIGHTEE data includes a $\sim$1deg$^{2}$ pointing of the COSMOS field observed at 1.28GHz with the MeerKAT radio telescope and cross-matched with multi-wavelength counterparts to provide classifications of high and low excitation radio galaxies (HERGs and LERGs) along with their corresponding host properties. We compare the properties of the MIGHTEE HERGs and LERGs with that predicted by the SIMBA simulations where HERGs and LERGs are defined as radio galaxies dominated by cold or hot mode accretion respectively. We consider stellar masses $\;{M}_{*}$, star formation rates SFR, AGN bolometric luminosity $L_{\rm bol}$, and Eddington fraction $f_{\rm Edd}$, as a function of 1.4GHz radio luminosity and redshift. In both MIGHTEE and SIMBA, the properties of HERGs and LERGs are similar across all properties apart from SFRs due to differences in host cold gas content in SIMBA. We predict a population of HERGs with low $f_{\rm Edd}$ in SIMBA that are confirmed in the MIGHTEE observations and tied to the faint population at low $z$. The predictions from SIMBA with the MIGHTEE observations describe a regime where our understanding of the radio galaxy dichotomy breaks down, challenging our understanding of the role of AGN accretion and feedback in the faint population of radio galaxies.

We investigate the efficacy of a systematic planetary nebula (PN) search in the Census of the Local Universe (CLU) narrowband (H$\alpha$) survey that covers a considerably larger sky region of above declination $-20^\circ$ than most previous surveys. Using PNe observed by the Isaac Newton Telescope Photometric H$\alpha$ Survey (IPHAS) as validation, we are able to visually recover 432 out of 441 cataloged PNe (98\%) within the CLU dataset, with 5 sources having unusable CLU images and 4 missed due to limitations of imaging quality. Moreover, the reference PNe are conventionally divided into three PN classes in decreasing order of identification confidence given their spectra and morphologies. We record consistently high recovery rate across all classes: 95\% of True, 71\% of Likely, and 81\% of Possible sources are readily recovered. To further demonstrate the ability of CLU to find new PNe, we undertake a preliminary search of compact PNe within a sub-region of the validation catalog, mainly utilizing the significance of narrow-band colors ($\Sigma$) as a metric for identification. In a $200\,\rm deg^2$ region, we search the CLU source catalog and find 31 PN candidates after automated and visual scrutiny, of which 12 are new sources not appearing in previous studies. As a demonstration of our ongoing follow-up campaign, we present medium-resolution optical spectra of six candidates and notice that four of them show emission signatures characteristic of confirmed PNe. As we refine our selection methods, CLU promises to provide a systematic catalog of PNe spanning $2/3$ of the sky.

Observations of the intracluster medium (ICM) in the outskirts of galaxy clusters reveal shocks associated with gas accretion from the cosmic web. Previous work based on non-radiative cosmological hydrodynamical simulations have defined the shock radius, $r_\text{shock}$, using the ICM entropy, $K \propto T/{n_\mathrm{e}}^{2/3}$, where $T$ and $n_\text{e}$ are the ICM temperature and electron density respectively; the $r_\text{shock}$ is identified with either the radius at which $K$ is a maximum or at which its logarithmic slope is a minimum. We investigate the relationship between $r_\text{shock}$, which is driven by gravitational hydrodynamics and shocks, and the splashback radius, $r_\text{splash}$, which is driven by the gravitational dynamics of cluster stars and dark matter and is measured from their mass profile. Using 324 clusters from {\small The Three Hundred} project of cosmological galaxy formation simulations, we quantify statistically how $r_\text{shock}$ relates to $r_\text{splash}$. Depending on our definition, we find that the median $r_\text{shock} \simeq 1.38 r_\text{splash} (2.58 R_{200})$ when $K$ reaches its maximum and $r_\text{shock} \simeq 1.91 r_\text{splash} (3.54 R_{200})$ when its logarithmic slope is a minimum; the best-fit linear relation increases as $r_\text{shock} \propto 0.65 r_\text{splash}$. We find that $r_\text{shock}/R_{200}$ and $r_\text{splash}/R_{200}$ anti-correlate with virial mass, $M_{200}$, and recent mass accretion history, and $r_\text{shock}/r_\text{splash}$ tends to be larger for clusters with higher recent accretion rates. We discuss prospects for measuring $r_\text{shock}$ observationally and how the relationship between $r_\text{shock}$ and $r_\text{splash}$ can be used to improve constraints from radio, X-ray, and thermal Sunyaev-Zeldovich surveys that target the interface between the cosmic web and clusters.

Gravitational wave (GW) detection and analysis rely heavily on efficient Signal to Noise Ratio (SNR) computation, a critical parameter for identifying and characterizing astrophysical signals from detectors like LIGO, Virgo, and KAGRA. The gwsnr Python package addresses the computational challenges in GW research by providing a flexible, user-friendly framework for SNR calculations. It supports diverse detector noise models, waveform approximants, detector configurations, and signal parameters, making it highly adaptable to various research needs. The package introduces innovative methods, including a partial-scaling interpolation technique for spin-less binary systems, a noise-weighted inner product approach with multiprocessing capabilities, and an Artificial Neural Network (ANN) model for rapid detection probability estimation of precessing binary black hole systems. Performance is further optimized using numba's Just-In-Time (JIT) compiler, enabling significant reductions in execution time. Seamlessly integrating with other tools, such as the ler package, gwsnr is instrumental in simulating detectable compact binary mergers, estimating GW event rates, and analyzing selection effects in hierarchical Bayesian frameworks. This comprehensive toolkit empowers researchers to conduct efficient, precise, and scalable analyses of GW signals.

Stripped-envelope supernovae (SNe) are H-poor transients produced at the end of the life of massive stars that previously lost their H-rich envelope. Their progenitors are thought to be donor stars in mass-transferring binary systems, which were stripped of their H-rich envelopes some $10^6$yr before core collapse. A subset of the stripped-envelope SNe exhibit spectral and photometric features indicative of interaction between their ejecta and nearby circumstellar material (CSM). We examine whether mass transfer during, or shortly before, core collapse in massive binary systems can produce the CSM inferred from the observations of interacting H-poor SNe. We select 44 models from a comprehensive grid of detailed binary evolution models in which the mass donors are H-free and explode while transferring mass to a main-sequence companion. We find that in these models, mass transfer starts less than $\sim20$kyr before, and often continues until the core collapse of the donor star. Up to $0.8M_\odot$ of H-free material are removed from the donor star during this phase, which may produce a He-rich circumbinary material. We explore plausible assumptions for its spatial distribution at the time of explosion. When assuming that the CSM accumulates in a circumbinary disk, we find qualitative agreement with the supernova and CSM properties inferred from observed Type Ibn SNe, and to a lesser extent with constraints from Type Icn SNe. We find that our mass transferring stripped envelope SN progenitor models may produce up to $\sim$10% of all stripped envelope supernovae. The binary channel proposed in this work can qualitatively account for the observed key properties and rate of interacting H-poor SNe. Models for the evolution of the circumbinary material and the spectral evolution of exploding progenitors from this channel are needed to further test its significance.

Soliton cores at the center of fuzzy dark matter (FDM) halos provide a promising way to distinguish FDM from other dark matter models. However, the relation between solitons and their host halos remains contentious. Here, we rigorously examine this soliton-halo relation (SHR) using a rich set of cosmological simulations across various FDM particle masses, halo masses, and redshifts. We explicitly demonstrate thermal equilibrium between solitons and surrounding halo granules, energy equipartition within halos, and an FDM concentration-mass-nonisothermality relation. For each FDM halo, we confirm that its density profile outside the central soliton matches a collisionless N-body simulation from the same initial condition, serving as stringent numerical convergence tests. Our refined SHR agrees well with virialized halos in simulations, with a $1\sigma$ deviation of less than $30\%$. These findings not only reaffirm the SHR proposed by Schive et al. (2014) but also offer a more comprehensive understanding that extends its applicability. The simulation code GAMER is accessible at this https URL. A Python script for computing the theoretical SHR is available at this https URL.

Vertical shear instability (VSI), driven by a vertical gradient of rotational angular velocity, is a promising source of turbulence in protoplanetary disks. We examine the semi-global stability of thermally stratified disks and find that the VSI consists of surface and body modes: surface modes are confined to regions of strong shear, while body modes extend perturbations across the disk, consistent with the previous findings. In thermally stratified disks, surface modes bifurcate into two branches. The branch associated with the strongest shear at mid-height exhibits a higher growth rate compared to the branch near the surfaces. Surface modes generally grow rapidly and require a high radial wave number $k_R$, whereas body mode growth rates increase as $k_R$ decreases. Thermal stratification enhances the growth rates of both surface and body modes and boosts VSI-driven radial kinetic energy relative to vertical energy. Our results suggest that simulations will initially favor surface modes with large $k_R$, followed by an increase in body modes with smaller $k_R$, with faster progression in more thermal stratified disks.

We conduct three-dimensional hydrodynamic simulations to investigate the nonlinear outcomes and observability of vertical shear instability (VSI) in protoplanetary disks. Our models include both vertically isothermal and thermally stratified disks, with the latter representing realistic conditions featuring a hotter atmosphere above the midplane. We find that the VSI grows more rapidly and becomes stronger in thermally stratified disks due to enhanced shear, resulting in higher levels of turbulence. At saturation, the turbulence stress reaches $\alpha_{R\phi}\gtrsim 10^{-3}$, more than an order of magnitude stronger than the isothermal case. The saturated turbulence is more pronounced near the disk surfaces than at the midplane. On synthetic velocity residual maps, obtained by subtracting the Keplerian rotational velocity, perturbations driven by the VSI manifest as axisymmetric rings in isothermal disks and as ring segments in thermally stratified disks. The latter are visible at disk inclinations as high as $45^\circ$ in thermally stratified disks. The amplitudes of these residual velocities range from $\sim 50$ to $\sim100$ $\mathrm{m\ s}^{-1}$ at a $20^\circ$ inclination, with larger values corresponding to greater thermal stratification. The magnitude of the observed velocity residual increases with the optical depth of the tracer used, as optically thick lines probe the regions near the disk surfaces.

We conducted an analysis of the velocity field of dwarf galaxies in the LITTLE THINGS sample, focusing on deriving 2D velocity maps that encompass both the transverse and radial velocity fields. Within the range of radial distances where velocity anisotropies are sufficiently small for the disc to be considered rotationally supported, and where the warped geometry of the disc can be neglected, we reconstructed the rotation curve while taking into account the effect of the asymmetric drift. To fit the rotation curves, we employed the standard halo model and the dark matter disc (DMD) model, which assumes that dark matter is primarily confined to the galactic discs and can be traced by the distribution of \HI{}. Interestingly, our analysis revealed that the fits from the DMD model are statistically comparable to those obtained using the standard halo model, but the inferred masses of the galaxies in the DMD model are approximately 10 to 100 times smaller than the masses inferred in the standard halo model. In the DMD model, the inner slope of the rotation curve is directly related to a linear combination of the surface density profiles of the stellar and gas components, which generally exhibit a flat core. Consequently, the observation of a linear relationship between the rotation curve and the radius in the disc central regions is consistent with the framework of the DMD model.

Anisotropic stochastic gravitational wave background (SGWB) serves as a potential probe of the large-scale structure (LSS) of the universe. In this work, we explore the anisotropic SGWB from local ($z < \sim 0.085$) merging stellar mass compact binaries, specifically focusing on merging stellar binary black holes, merging neutron-star-black-hole binaries, and merging binary neutron stars. The analysis employs seven all-sky mock lightcone gravitational wave event catalogues, which are derived from the Millennium simulation combined with a semi-analytic model of galaxy formation and a binary population synthesis model. We calculate the angular power spectra $\mathrm{C}_\ell$ at multipole moments $\ell$, expressed as $\text{log}_{10} [\ell(\ell+1)\mathrm{C}_\ell/(2\pi)]$, based on the skymaps of the overdensity $\delta_\mathrm{GW}$ in the anisotropic SGWB. The spectra for all three source types exhibit an approximately linear increase with $\text{log}_{10} \ell$ at higher $\ell$ (e.g., $\ell > \sim 30 - 300$) in seven catalogues, with a characteristic slope of $\sim 2$. The spectra of seven catalogues exhibit considerable variations, arising from fluctuations in spatial distribution, primarily in the radial distribution, of nearby sources (e.g., $< 50$ Mpc/h). After subtracting these nearby sources, the variations become much smaller and the spectra for the three source types become closely aligned (within discrepancies of a factor of $\sim 2$ across $\ell = 1 - 1000$ for all catalogues).

This study employs a novel approach for reconstructing the thermal Sunyaev-Zeldovich (tSZ) effect power spectrum from Planck data using the Analytical Blind Separation (ABS) method. The ABS method improves the recovery of weak signals, by applying eigenmode exclusion for low signal-to-noise ratio regimes and introducing a shift parameter to stabilize calculations. Validation through simulated Planck data demonstrates the robustness of ABS in reconstructing the tSZ power spectrum, even under challenging conditions. In the analysis of the Planck PR3 full-mission data, ABS shows lower amplitudes at $\ell \gtrsim 300$ compared to the Planck 2015 band powers using the MILCA and NILC foreground cleaning methods. After marginalizing over residual foreground components, the ABS analysis finds the overall tSZ power spectrum amplitude to be 20\% lower than the Planck best-fit one, suggesting a smaller $S_8$. At $\ell \simeq 200$, the tSZ band power is $10^{12} \ell(\ell+1)C^{yy}_\ell/(2\pi) = 0.126\pm 0.018$, largely independent of the tSZ model choice. These findings highlight the potential of the ABS method as a promising alternative for tSZ power spectrum analysis, offering a robust and independent means of extracting cosmological parameters.

It is well-known that the cores of massive stars sustain a stellar dynamo with a complex magnetic field configuration. However, the same cannot be said for the field's strength and geometry at the convective-radiative boundary, which are crucial when performing asteroseismic inference. In this Letter, we present three-dimensional (3D) magnetohydrodynamic (MHD) simulations of a 7 solar mass mid-main sequence star, with particular attention given to the convective-radiative boundary in the near-core region. Our simulations reveal that the toroidal magnetic field is significantly stronger than the poloidal field in this region, contrary to recent assumptions. Moreover, the rotational shear layer, also important for asteroseismic inference, is specifically confined within the extent of the buoyancy frequency peak. These results, which are based on the inferred properties of HD 43317, have widespread implications for asteroseismic studies of rotation, mixing and magnetism in stars. While we expect our results to be broadly applicable across stars with similar buoyancy frequency profiles and stellar masses, we also expect the MHD parameters and the initial stellar rotation rate to impact the geometry of the field and differential rotation at the convective-radiative interface.

We study the absorption of a gamma-ray burst afterglow in a dense molecular cloud in the X-ray wavelength range. We report the results of numerical simulations of the propagation of the gamma-ray burst radiation in the cloud for various gas densities, metallicities, and distances from the gamma-ray burst progenitor star and the cloud. We consider a sample of 45 gamma-ray bursts with known redshifts in which the isotropic-equivalent gamma-ray energy is approximately equal to the value adopted in our numerical simulations. For these gamma-ray bursts, we have analyzed the Swift/XRT energy spectra of their afterglows at late times, $t \geq 4 \times 10^3$ s. It is shown that the hydrogen column densities estimated using the absorption model in which the ionization of metal ions is not taken into account and the solar metallicity is used are a factor of 1-3 smaller than the actual values - if the molecular cloud is located close to the gamma-ray burst progenitor star. If the gas-dust cloud is located at a distance of $R \geq 10$ pc from the source of the gamma-ray burst or the gas metallicity is $[M/H] \leq -1$, then the effect of the cloud ionization structure on the absorption of the afterglow is minor.

Context. The characterization of exoplanets requires a good description of the host star. Stellar activity acts as a source of noise which can alter planet radii as derived from the transit depth or atmospheric characterization. Aims. Here, we propose PAStar, a model to describe photospheric activity in the form of spots and faculae which could be applied to a wide range of stellar observations, from photometric to spectroscopic time series, to be able to correctly extract planetary and stellar properties. Methods. The adopted stellar atmosphere is a combination of three components, the quiet photosphere, spots and faculae. The model takes into account the effects of star inclination, doppler shifts due to stellar rotation as well as for limb darkening, independent for each component. Several synthetic products have been presented to show the capabilities of the model. Results. The model is able to retrieve the input surface inhomogeneities configuration through photometric or spectroscopic observations. The model has been validated against optical solar data and compared to alternative stellar surface activity models; e.g. SOAP code. The Sun is a unique laboratory to test stellar models because of the possibility to relate unambiguously flux variations to surface inhomogeneities configuration. This validation has been done by analyzing a photometric time series from the VIRGO photometer on board of SOHO mission. Results have been compared to real solar images from the HMI instrument on board of SDO to confirm the goodness of the results in terms of surface inhomogeneities position and dimensions. Conclusions. The description of stellar activity is a fundamental step in several astrophysical contexts and it is covered by the method we have presented. Our model offers a flexible and valuable tool to describe the activity of stars when it is dominated by spots and faculae.

Structure functions, which represent the moments of the increments of a stochastic process, are essential complementary statistics to power spectra for analysing the self-similar behaviour of a time series. However, many real-world environmental datasets, such as those collected by spacecraft monitoring the solar wind, contain gaps, which inevitably corrupt the statistics. The nature of this corruption for structure functions remains poorly understood - indeed, often overlooked. Here we simulate gaps in a large set of magnetic field intervals from Parker Solar Probe in order to characterize the behaviour of the structure function of a sparse time series of solar wind turbulence. We quantify the resultant error with regards to the overall shape of the structure function, and its slope in the inertial range. Noting the consistent underestimation of the true curve when using linear interpolation, we demonstrate the ability of an empirical correction factor to de-bias these estimates. This correction, "learnt" from the data from a single spacecraft, is shown to generalize well to data from a solar wind regime elsewhere in the heliosphere, producing smaller errors, on average, for missing fractions >25%. Given this success, we apply the correction to gap-affected Voyager intervals from the inner heliosheath and local interstellar medium, obtaining spectral indices similar to those from previous studies. This work provides a tool for future studies of fragmented solar wind time series, such as those from Voyager, MAVEN, and OMNI, as well as sparsely-sampled astrophysical and geophysical processes more generally.

Ram pressure stripped galaxies are rare cases of environmental evolution in action. However, our ability to understand these galaxies is limited by the small number of identified galaxies experiencing ram pressure stripping (RPS). Our aim is to explore the efficacy of citizen science classifications in identifying ram pressure stripped galaxies, and use this to aid in motivating new samples of ram pressure stripped candidates. We compile a sample of over 200 known ram pressure stripped galaxies from existing literature, with morphological classifications obtained from Galaxy Zoo. We compare these galaxies with magnitude and redshift-matched comparison cluster and field galaxies. Additionally, we create a sample of SDSS cluster galaxies, with morphological classifications similar to known ram pressure stripped galaxies, and compare the fraction of potential new RPS candidates against control samples. We find that ram pressure stripped galaxies exhibit a higher proportion of odd and irregular morphological classifications compared to field and cluster comparison samples. This trend is particularly pronounced in galaxies displaying strong optical ram pressure stripping features, but absent from galaxies with only radio tails. We find that SDSS galaxies with Galaxy Zoo classifications consistent with the known RPS galaxies have a higher fraction of visible ram pressure stripping features ($19\%$) compared with other cluster galaxies ($12\%$) when classified by experts. We identify 101 new ram pressure stripping candidate galaxies through these expert classifications. We demonstrate that indirect morphological classifications from citizen science projects can increase the efficiency in which new stripping candidates are found. Projects such as Galaxy Zoo can aid in the identification of ram pressure stripped galaxies that are key to understanding galaxy evolution in clusters.

Calibration of radio interferometer data ought to be a solved problem; it has been an integral part of data reduction for some time. However, as larger, more sensitive radio interferometers are conceived and built, the calibration problem grows in both size and difficulty. The increasing size can be attributed to the fact that the data volume scales quadratically with the number of antennas in an array. Additionally, new instruments may have up to two orders of magnitude more channels than their predecessors. Simultaneously, increasing sensitivity is making calibration more challenging: low-level RFI and calibration artefacts (in the resulting images) which would previously have been subsumed by the noise may now limit dynamic range and, ultimately, the derived science. It is against this backdrop that we introduce QuartiCal: a new Python package implementing radio interferometric calibration routines. QuartiCal improves upon its predecessor, CubiCal, in terms of both flexibility and performance. Whilst the same mathematical framework - complex optimization using Wirtinger derivatives - is in use, the approach has been refined to support arbitrary length chains of parameterized gain terms. QuartiCal utilizes Dask, a library for parallel computing in Python, to express calibration as an embarrassingly parallel task graph. These task graphs can (with some constraints) be mapped onto a number of different hardware configurations, allowing QuartiCal to scale from running locally on consumer hardware to a distributed, cloud-based cluster. QuartiCal's qualitative behaviour is demonstrated using MeerKAT observations of PSR J2009-2026. These qualitative results are followed by an analysis of QuartiCal's performance in terms of wall time and memory footprint for a number of calibration scenarios and hardware configurations.

The popularity of the CLEAN algorithm in radio interferometric imaging stems from its maturity, speed, and robustness. While many alternatives have been proposed in the literature, none have achieved mainstream adoption by astronomers working with data from interferometric arrays operating in the big data regime. This lack of adoption is largely due to increased computational complexity, absence of mature implementations, and the need for astronomers to tune obscure algorithmic parameters. This work introduces pfb-imaging: a flexible library that implements the scaffolding required to develop and accelerate general radio interferometric imaging algorithms. We demonstrate how the framework can be used to implement a sparsity-based image reconstruction technique known as (unconstrained) SARA in a way that scales with image size rather than data volume and features interpretable algorithmic parameters. The implementation is validated on terabyte-sized data from the MeerKAT telescope, using both a single compute node and Amazon Web Services computing instances.

Stimela2 is a new-generation framework for developing data reduction workflows. It is designed for radio astronomy data but can be adapted for other data processing applications. Stimela2 aims at the middle ground between ease of development, human readability, and enabling robust, scalable and reproducible workflows. It represents workflows by linear, concise and intuitive YAML-format "recipes". Atomic data reduction tasks (binary executables, Python functions and code, and CASA tasks) are described by YAML-format "cab definitions" detailing each task's "schema" (inputs and outputs). Stimela2 provides a rich syntax for chaining tasks together, and encourages a high degree of modularity: recipes may be nested into other recipes, and configuration is cleanly separated from recipe logic. Tasks can be executed natively or in isolated environments using containerization technologies such as Apptainer. The container images are open-source and maintained through a companion package called cult-cargo. This enables the development of system-agnostic and fully reproducible workflows. Stimela2 facilitates the deployment of scalable, distributed workflows by interfacing with the Slurm scheduler and the Kubernetes API. The latter allows workflows to be readily deployed in the cloud. Previous papers in this series used Stimela2 as the underlying technology to run workflows on the AWS cloud. This paper presents an overview of Stimela2's design, architecture and use in the radio astronomy context.

Most of the baryonic matter of the Universe resides in a highly-ionized gaseous intergalactic medium. This gas flows along dark-matter filaments toward galaxy superclusters, clusters, and groups until it pools around the galaxies into a circumgalactic medium. Eventually, the gas settles into the interstellar medium of the galaxies, where it fuels the successive generations of star formation that ultimately produce the stars and heavy elements that make up galaxies today. The gas has been studied for decades using absorption lines produced by Hydrogen and various ions of heavy elements in the spectra of background quasi-stellar objects (QSOs). But directly imaging the extremely faint glow of this "cosmic web" of intergalactic and circumgalactic gas has remained an elusive goal of observational cosmology. Some recent progress has been made by using imaging spectrographs to record high-redshift Ly$\alpha$ emission, although over only very narrow fields of view. Here we report direct images of intergalactic and circumgalactic gas in the distant Universe obtained using the Condor Array Telescope that we purposely built to reach extremely low-surface-brightness sensitivities over very wide fields of view. We show that these images directly detect and characterize the imprint of Ly$\alpha$ emission from the cosmic web at an overwhelming statistical significance. By stacking portions of the images centered on tens of thousands of galaxies of known redshift, we show that they also reveal extremely faint emission from H$^0$, C$^{3+}$, and Mg$^+$ and absorption from cosmic dust in the tenuous outskirts of the galaxies. Our results demonstrate that sensitive imaging observations can now detect and characterize emission (and absorption) from the cosmic web of intergalactic and circumgalactic gas (and dust).

Artificial intelligence (AI) is revolutionizing research by enabling the efficient analysis of large datasets and the discovery of hidden patterns. In astrophysics, AI has become essential, transforming the classification of celestial sources, data modeling, and the interpretation of observations. In this review, I highlight examples of AI applications in astrophysics, including source classification, spectral energy distribution modeling, and discuss the advancements achievable through generative AI. However, the use of AI introduces challenges, including biases, errors, and the "black box" nature of AI models, which must be resolved before their application. These issues can be addressed through the concept of Human-Guided AI (HG-AI), which integrates human expertise and domain-specific knowledge into AI applications. This approach aims to ensure that AI is applied in a robust, interpretable, and ethical manner, leading to deeper insights and fostering scientific excellence.

The pCam6060 photodetecting system was developed at SAO RAS and is based on the GSENSE6060BSI photodetector manufactured by GPixel (China) with a frame format of 6144x6144 active pixels and a pixel size of 10 mkm. The readout speed reached 11 fps. The back-illuminated detector has a wide spectral range of 200-1040 nm with a minimum quantum efficiency (QE) of 10% and maximum sensitivity of 95% at 580 nm. The quantum efficiency in the near-infrared range was 58% at 850 nm. The pCam6060 system controller implements a mode of simultaneous image readout via two 12-bit video channels with different gain and their subsequent combination in the controller into a single frame with an extended 16-bit dynamic range. This method simultaneously achieves a low readout noise level (about 3 e-) in the high-gain channel and a large dynamic range (full well capacity of about 100000 e-) in the low-gain channel. Back-illuminated CMOS detectors, unlike front-illuminated devices, do not have the effect of long-term preserving the residual charge from previous exposures, which makes them suitable for recording faint objects in photometric long-exposure observation methods. Communication between the host computer and the camera was carried out via a fiber optic line at distances of up to 50 m. Video data are recorded on the computer hard drive in real-time. The pCam6060 photodetecting system is designed for astronomical applications and has a moisture-proof design.

Bayesian analysis has begun to be more widely adopted in X-ray spectroscopy, but it has largely been constrained to relatively simple physical models due to limitations in X-ray modelling software and computation time. As a result, Bayesian analysis of numerical models with high physics complexity have remained out of reach. This is a challenge, for example when modelling the X-ray emission of accreting black hole X-ray binaries, where the slow model computations severely limit explorations of parameter space and may bias the inference of astrophysical parameters. Here, we present RTFAST-Spectra: a neural network emulator that acts as a drop in replacement for the spectral portion of the black hole X-ray reverberation model RTDIST. This is the first emulator for the reltrans model suite and the first emulator for a state-of-the-art x-ray reflection model incorporating relativistic effects with 17 physically meaningful model parameters. We use Principal Component Analysis to create a light-weight neural network that is able to preserve correlations between complex atomic lines and simple continuum, enabling consistent modelling of key parameters of scientific interest. We achieve a $\mathcal{O}(10^2)$ times speed up over the original model in the most conservative conditions with $\mathcal{O}(1\%)$ precision over all 17 free parameters in the original numerical model, taking full posterior fits from months to hours. We employ Markov Chain Monte Carlo sampling to show how we can better explore the posteriors of model parameters in simulated data and discuss the complexities in interpreting the model when fitting real data.

We examine ejecta generated by ultra low velocity impacts under asteroid conditions. In an environment of precisely controlled milligravity and under vacuum, impacts with velocities in the range of centimeters/second are performed with irregularly shaped impactors onto granular beds. The resulting ejecta velocities are compared to existing literature values and extend the observed systematic trends towards lower impact energies, broadening the parameter range. Simulations are performed to reason the systematics and the absence thereof for measurements performed at earth gravity. We find, that the cutoff induced by gravity dependent minimal observable velocities plays a crucial role in the values obtained for mean ejecta velocities.

Sympathetic solar flares are eruptions that occur nearby in space and time, driven by an apparent interaction between the active regions in which they are triggered. Their statistical existence on the Sun has yet to be firmly established. The main goal of this paper is to identify a statistical signature of sympathetic flares, characterize their properties and determine a potential mechanism driving their interaction. We perform a statistical analysis of a large number of flares observed by the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO), the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and the Spectrometer Telescope for Imaging X-rays (STIX) on Solar Orbiter during solar cycle 24 and 25. We examine the spatiotemporal distribution of consecutive flare pairs across solar cycle phases and hemispheres, along with the propagation velocity of potential causal interactions and the relationship between flare magnitudes. We observe an excess of hemispheric flares separated by about 30 degrees of longitude and triggered in less than 1.5 hours from each other. This peak in angular separation varies with the solar cycle phase and hemisphere. Moreover, we identify a deficit of transequatorial events separated by 25-30 degrees in latitude and less than 5 degrees in longitude, a phenomenon we term unsympathetic flares. We provide strong statistical evidence for the existence of sympathetic flares on the Sun, demonstrating that their occurrence rate reaches approximately 5% across the three instruments used in this study. Additionally, we propose an interpretation of the observed angular scale of the sympathetic phenomenon, based on the separation between magnetic field line footpoints derived from potential field source surface extrapolations.

Context. The study of open cluster morphology is pivotal for exploring their formation and evolutionary processes. Aims. We manage to assess the morphological stability of 105 nearby open clusters within tidal radii on the X-Y, X-Z, and Y-Z planes of the heliocentric Cartesian coordinate system, utilizing member catalogs from the literature. Meanwhile, we also delve into factors potentially impacting the clusters' morphological stability on these projection planes. Methods. We used the rose diagram method by constructing the 3D projected morphology of sample clusters to quantify the morphological stability of their 3D projected morphology. Results. Our analysis indicates there is a demonstrated linear positive correlation between the number of sample clusters' member stars within tidal radii and their morphological stability in different 3D projection planes. This may suggest that the more member stars there are within the tidal radius of a cluster, the stronger its own gravitational binding capacity is, resulting in strong morphological stability. We find a direct link between the clusters' morphological stabilities in the X-Z plane within tidal radii and their velocity dispersion in the same plane, suggesting that the morphological stabilities in the X-Z plane may be dependent on internal dynamics. Moreover, the morphological stability of the clusters' 3D projection is influenced by their spatial positions along the Y axis, but not linearly, indicating that the environmental changes at the clusters' location may affect their morphological stability. Likewise, specific external forces can have an effect on their morphological stability. Conclusions. This research introduces a novel perspective for understanding the morphological stability of open clusters, with a particular focus on their 3D projected morphologies.

The braking index, $n$, of a pulsar is a measure of its angular momentum loss and the value it takes corresponds to different spin-down mechanisms. For a pulsar spinning down due to gravitational wave emission from the principal mass quadrupole mode alone, the braking index would equal exactly 5. Unfortunately, for millisecond pulsars, it can be hard to measure observationally due to the extremely small second time derivative of the rotation frequency, $\ddot{f}$. This paper aims to examine whether it could be possible to extract the distribution of $n$ for a whole population of pulsars rather than measuring the values individually. We use simulated data with an injected $n=5$ signal for 47 millisecond pulsars and extract the distribution using hierarchical Bayesian inference methods. We find that while possible, observation times of over 20 years and RMS noise of the order of $10^{-5}$ ms are needed, which can be compared to the mean noise value of $3\times10^{-4}$ ms for the recent wideband 12.5-year NANOGrav sample, which provided the pulsar timing data used in this paper.

Next generations of radio surveys are expected to identify tens of millions of new sources, and identifying and classifying their morphologies will require novel and more efficient methods. Self-Organising Maps (SOMs), a type of unsupervised machine learning, can be used to address this problem. We map 251,259 multi-Gaussian sources from Rapid ASKAP Continuum Survey (RACS) onto a SOM with discrete neurons. Similarity metrics, such as Euclidean distances, can be used to identify the best-matching neuron or unit (BMU) for each input image. We establish a reliability threshold by visually inspecting a subset of input images and their corresponding BMU. We label the individual neurons based on observed morphologies and these labels are included in our value-added catalogue of RACS sources. Sources for which the Euclidean distance to their BMU is $\lesssim$ 5 (accounting for approximately 79$\%$ of sources) have an estimated $>90\%$ reliability for their SOM-derived morphological labels. This reliability falls to less than 70$\%$ at Euclidean distances $\gtrsim$ 7. Beyond this threshold it is unlikely that the morphological label will accurately describe a given source. Our catalogue of complex radio sources from RACS with their SOM-derived morphological labels from this work will be made publicly available.

We present results from commissioning observations of the mid-IR instrument, MIRAC-5, on the 6.5-m MMT telescope. MIRAC-5 is a novel ground-based instrument that utilizes a state-of-the-art GeoSnap (2 - 13 microns) HgCdTe detector with adaptive optics support from MAPS to study protoplanetary disks, wide-orbit brown dwarfs, planetary companions in the contrast-limit, and a wide range of other astrophysical objects. We have used MIRAC-5 on six engineering observing runs, improving its performance and defining operating procedures. We characterize key aspects of MIRAC-5's performance, including verification that the total telescope, atmosphere, instrument, and detector throughput is approximately 10%. Following a planned dichroic upgrade, the system will have a throughput of 20% and background limiting magnitudes (for SNR = 5 and 8 hour exposure times) of 18.0, 15.6, and 12.6 for the L', M', and N' filters, respectively. The detector pixels experience 1/f noise but, if the astrophysical scene is properly modulated via chopping and nodding sequences, it is less than 10% the Poisson noise from the observed background in an 85 Hz frame. We achieve close to diffraction-limited performance in the N-band and all bands are expected to reach diffraction-limited performance following the adaptive optics system commissioning. We also present an exposure time calculator calibrated to the on-sky results. In its current state, MIRAC-5 will be capable of achieving several scientific objectives including the observation of warm wide-orbit companions. Once the adaptive optics is commissioned and a coronagraph installed in 2025, MIRAC-5 will have contrast-limited performance comparable to JWST, opening new and complementary science investigations for close-in companions.

We investigate the impact of CO$_2$ on TRAPPIST-1 e, f and g during the magma ocean stage. These potentially habitable rocky planets are currently the most accessible for astronomical observations. A constraint on the volatile budget during the magma ocean stage is a link to planet formation and also needed to judge their habitability. We perform simulations with 1-100 terrestrial oceans (TO) of H$_2$O with and without CO$_2$ and for albedos 0 and 0.75. The CO$_2$ mass is scaled with initial H$_2$O by a constant factor between 0.1 and 1. The magma ocean state of rocky planets begins with a CO$_2$-dominated atmosphere but can evolve into a H$_2$O dominated state, depending on initial conditions. For less than 10 TO initial H$_2$O, the atmosphere tends to desiccate and the evolution may end with a CO$_2$ dominated atmosphere. Otherwise, the final state is a thick (>1000 bar) H$_2$O-CO$_2$ atmosphere. Complete atmosphere desiccation with less than 10 TO initial H$_2$O can be significantly delayed for TRAPPIST-1e and f, when H$_2$O has to diffuse through a CO$_2$ atmosphere to reach the upper atmosphere, where XUV photolysis occurs. As a consequence of CO$_2$ diffusion-limited water loss, the time of mantle solidification for TRAPPIST-1 e, f, and g can be significantly extended compared to a pure H$_2$O evolution by up to 40 Myrs for albedo 0.75 and by up to 200 Mrys for albedo 0. The addition of CO$_2$ further results in a higher water content in the melt during the magma ocean stage. Our compositional model adjusted for the measured metallicity of TRAPPIST-1 yields for the dry inner planets (b, c, d) an iron fraction of 27 wt-%. For TRAPPIST-1 e, this iron fraction would be compatible with a (partly) desiccated evolution scenario and a CO$_2$ atmosphere with surface pressures of a few 100 bar. A comparative study between TRAPPIST-1 e and the inner planets may yield the most insights about formation and evolution scenarios.

Traditional methods for determining the radius of a 1.4 $M_{\odot}$ neutron star ($R_{1.4}$) rely on specific equations of state (EOS) models that describe various types of dense nuclear matter. This dependence on EOS models can introduce substantial systematic uncertainties, which may exceed the measurement uncertainties when constraining $R_{1.4}$. In this study, we explore a novel approach to constraining $R_{1.4}$ using data from NICER observations of PSR J0030+0451 (J0030) and PSR J0437-4715 (J0437). However, this work presents a more data-driven analysis framework, substantially decreasing the need for EOS assumptions. By analyzing the Mass-Radius measurements of these two neutron stars, we infer $R_{1.4}$ using statistical methods based mostly on observational data. We examine various hotspot configurations for J0030, along with new J0437 observations, and their effects on the inferred radius. Our results are consistent with X-ray timing, gravitational wave, and nuclear physics constraints, while avoiding EOS-related biases. The same method has also been applied to a simulated mass-radius dataset, based on our knowledge of future X-ray telescopes, demonstrating the model's ability to recover the injected $R_{1.4}$ value in certain cases. This method provides a data-driven pathway for extracting neutron star properties and offers a new approach for future observational efforts in neutron star astrophysics.

The Near-Earth Asteroid (NEA) 2024 PT5 is on an Earth-like orbit which remained in Earth's immediate vicinity for several months at the end of 2024. PT5's orbit is challenging to populate with asteroids originating from the Main Belt and is more commonly associated with rocket bodies mistakenly identified as natural objects or with debris ejected from impacts on the Moon. We obtained visible and near-infrared reflectance spectra of PT5 with the Lowell Discovery Telescope and NASA Infrared Telescope Facility on 2024 August 16. The combined reflectance spectrum matches lunar samples but does not match any known asteroid types -- it is pyroxene-rich while asteroids of comparable spectral redness are olivine-rich. Moreover, the amount of solar radiation pressure observed on the PT5 trajectory is orders of magnitude lower than what would be expected for an artificial object. We therefore conclude that 2024 PT5 is ejecta from an impact on the Moon, thus making PT5 the second NEA suggested to be sourced from the surface of the Moon. While one object might be an outlier, two suggest that there is an underlying population to be characterized. Long-term predictions of the position of 2024 PT5 are challenging due to the slow Earth encounters characteristic of objects in these orbits. A population of near-Earth objects which are sourced by the Moon would be important to characterize for understanding how impacts work on our nearest neighbor and for identifying the source regions of asteroids and meteorites from this under-studied population of objects on very Earth-like orbits.

We derive the characteristic scales for physical quantities of galaxies, such as mass, size, acceleration, and angular momentum, within the self-interacting ultralight dark matter (ULDM) model. Due to the small mass of ULDM, even minor self-interactions can drastically alter these scales in the Thomas-Fermi limit. Even in this limit, the characteristic mass can be determined by quantum pressure rather than by repulsive forces. We suggest that these characteristic scales are connected to certain mysteries of observed galaxies such as the universal acceleration scale or the constant surface density of galaxies. Oscillation of ULDM field can explain the current cosmological density of dark matter. Many cosmological constraints imply the energy scale $\tilde{m}$ of order of $10~eV$ and a GUT scale phase transition related to ULDM.

The redshifted 21cm signal from Cosmic Dawn promises to open a new window into the early history of our universe and enable the probing of an unprecedented comoving survey volume. In this work, we revisit the imprint of Warm Dark Matter (WDM) on the 21cm signal power spectrum using an updated implementation of the WDM effect in the public code $\texttt{21cmFast}$ and considering a single population of cosmic dawn galaxies. By focusing on inferring the WDM mass, we analyze the degeneracies between the latter and the astrophysics parameters characterizing star formation and X-ray heating and we emphasize the role of the threshold mass for star-forming galaxies, $M_{\rm turn}$. We study the capability of the recently built HERA telescope to reconstruct the WDM mass by adopting the statistical approach of simulation-based inference. We include a comparison of the per-parameter reconstruction quality for different number of simulations used in the training of the algorithm. Our results indicate that HERA could surpass current Lyman-$\alpha$ forest constraints if Cosmic Dawn galaxies exhibit a threshold mass $M_{\rm turn}\lesssim 10^{8}\, M_\odot$. The X-ray source properties considered in this study may also influence the strength of the WDM constraint for lower threshold masses.

For many pulsars, the scattering structures responsible for scintillation are typically dominated by a single, thin screen along the line of sight, which persists for years or decades. In recent years, an increasing number of doubly-lensed events have been observed, where a secondary lens crosses the line of sight. This causes additional or distorted scintillation arcs over time scales ranging from days to months. In this work we report such a transient event for pulsar B1737+13 and propose a possible lensing geometry including the distance to both lenses, and the orientation of the main screen. Using phase retrieval techniques to separate the two lenses in the wavefield, we report a curvature and rate of motion of features associated with the secondary lens as it passed through the line of sight. By fitting the annual variation of the curvature, we report a possible distance and orientation for the main screen. The distance of the secondary lens is found by mapping the secondary feature onto the sky and tracking its position over time for different distances. We validate this method using B0834+06, for which the screen solutions are known through VLBI, and successfully recover the correct solution for the secondary feature. With the identified lensing geometry, we are able to estimate the size of the secondary lens, 1 - 3 au. Although this an appropriate size for a structure that could cause an extreme scattering event, we do not have conclusive evidence for or against that possibility.

We study the variability of the thermal (accretion disc) and non-thermal (jet) emission of thirteen flat spectrum radio quasars in the optical and near infrared (OIR) regime using light curves spanning years with an average sampling of three observations per week. We fit a combination of a blackbody and a power-law function to the OIR data, in the blazar rest frame, to extract the corresponding thermal (disc) and non-thermal (jet) components from the total flux. We carry out this analysis for the entire duration of the light curves to obtain the variation of the disc and jet components over years. Reliability of our fits have been affirmed by successfully retrieving accurate parameters by employing our method to simulated data and by comparing our results with published disc luminosity obtained by other methods for a few well-observed blazars. In blazars, the thermal (disc) emission is difficult to extract because the relativistically beamed radiation of the jet dominates at all wavelengths. By employing this method, the disc emission in blazars may be estimated directly from photometric data at OIR bands instead of indirect methods, such as, inferring it from the emission line luminosities. We find that the variability of the disc and jet emission obtained by the above method are strongly correlated in most cases.

Photographic spectra of Plaskett's Star (PS; HR2420, HD47129, V640 Mon) that were recorded at the Dominion Astrophysical Observatory (DAO) have been digitized with a flatbed scanner. Many of the spectra were recorded during campaigns in 1922 and 1937, and sample wavelengths between 0.39 and 0.50um. Spectra of poor quality are identified. Mean spectra near orbital phases 0.25 and 0.75 match many characteristics of synthetic spectra, although Hgamma and HeI 4388 are exceptions. Evidence is presented that Hgamma was affected by transient activity in 1937, but not in 1922. Emission lines of NIII and HeII move with wavelength in a manner that is consistent with them tracking the motion of the secondary, indicating that an 'f' spectral type designation should be assigned to the secondary. The location of the peak that is associated with the secondary in cross-correlation functions changes with time near phase 0.75, although the mean amplitude of the radial velocity curve of the secondary did not change between the two campaigns. There is also an offset in velocities of the primary measured from Hgamma and HeI 4472 near phase 0.25. The velocity curves of the components suggest a mass ratio that is larger than previous estimates, although uncertainties associated with the spectroscopic features attributed to the secondary, coupled with the wavelength resolution of the spectra, complicate efforts to determine robust masses. We conclude that peculiarities in the radial velocity curves of the components have thus been in place for over a thousand orbital cycles.

Previous Hubble Space Telescope (HST) observations of the star-forming cluster NGC 346 in the Small Magellanic Cloud (SMC) had revealed a large population of pre-main sequence (PMS) candidates, characterised by Halpha excess emission in their photometry. However, without access to spectroscopy, the nature of these objects remained unclear. Using the NIRSpec instrument on board JWST, we studied a sample of these stars, with masses in the range ~0.9-1.8 Msun, effective temperatures in the range 4,500-8,000 K, and PMS ages between ~0.1 and 30 Myr. Here we present the first spectra of solar-mass PMS stars in the metal-poor SMC (Z=1/8 Zsun) and discuss the physical properties of ten representative sources with good signal-to-noise ratio. The observations indicate that even the oldest of these PMS candidates are still accreting gas with typical rates of ~10^{-8} Msun/yr for stars older than ~10 Myr, confirming their PMS nature. The spectra also reveal near-infrared excess and molecular hydrogen excitation lines consistent with the presence of discs around these stars. These findings suggest that in a low-metallicity environment circumstellar discs can live longer than previously thought.

We study realistic models predicting primordial black hole (PBH) formation from density fluctuations generated in a first-order phase transition. We show that the second-order correction in the expansion of the bubble nucleation rate is necessary for accurate predictions and quantify its impact on the abundance of PBHs and gravitational waves (GWs). We find that the distribution of the fluctuations becomes more Gaussian as the second-order term increases. Consequently, models that predict the same PBH abundances can produce different GW spectra.

Considering inflationary magnetogenesis induced by time-dependent kinetic and axial couplings of a massless Abelian vector boson field breaking the conformal invariance we show in this Letter that, surprisingly, the spectral shape of the primordial magnetic field power spectrum is insensitive to the post-inflationary history, namely the barotropic parameter ($w$) and the gauge coupling functions of the post-inflationary era.

The electron spectrum exhibits a complex structure and has controversially proposed this http URL work reproduce the evolution of the electron spectrum based on a spatially dependent propagation (SDP) model. The key point is that our SPD model features two diffusion regions leading to two diffusion timescales, competing with the cooling timescale. This results in a three-segment power-law electron spectrum: (1) The spectrum below tens of GeV is primarily influenced by cooling effects from distant sources. (2)the spectrum dominated by diffusion effects from nearby sources from tens of GeV to TeV; (3) the spectrum above TeV, which is predominantly governed by cooling effects from nearby sources. This evolution is unique to the SDP model, and we offer a comprehensive and clear depiction of electron evolution under a single propagation scenario for the first time.

At sub-Kelvin temperatures, two-level systems (TLS) present in amorphous dielectrics source a permittivity noise, degrading the performance of a wide range of devices using superconductive resonators such as qubits or kinetic inductance detectors. We report here on measurements of TLS noise in hydrogenated amorphous silicon (a-Si:H) films deposited by plasma-enhanced chemical vapor deposition (PECVD) in superconductive lumped-element resonators using parallel-plate capacitors (PPCs). The TLS noise results presented in this article for two recipes of a-Si:H improve on the best achieved in the literature by a factor >5 for a-Si:H and other amorphous dielectrics and are comparable to those observed for resonators deposited on crystalline dielectrics.

We study an inflation model with nonminimal derivative coupling that features a coupling between the derivative of the inflaton field and the Einstein tensor. This model naturally amplifies curvature perturbations at small scales via gravitationally enhanced friction, a mechanism critical for the formation of primordial black holes and the associated production of potentially detectable scalar-induced gravitational waves. We derive analytical expressions for the primordial power spectrum, enabling efficient exploration of the model parameter space without requiring computationally intensive numerical solutions of the Mukhanov-Sasaki equation. Using the third data release of the Parkes Pulsar Timing Array (PPTA DR3), we constrain the model parameters characterizing the coupling function: $\phi_c = 3.7^{+0.3}_{-0.5} M_\mathrm{P}$, $\log_{10} \omega_L = 7.1^{+0.6}_{-0.3}$, and $\log_{10} \sigma = -8.3^{+0.3}_{-0.6}$ at 90\% confidence level. Our results demonstrate the growing capability of pulsar timing arrays to probe early Universe physics, complementing traditional cosmic microwave background observations by providing unique constraints on inflationary dynamics at small scales.

Gravitational-wave observatories like LIGO are large-scale, terrestrial instruments housed in infrastructure that spans a multi-kilometer geographic area and which must be actively controlled to maintain operational stability for long observation periods. Despite exquisite seismic isolation, they remain susceptible to seismic noise and other terrestrial disturbances that can couple undesirable vibrations into the instrumental infrastructure, potentially leading to control instabilities or noise artifacts in the detector output. It is, therefore, critical to characterize the seismic state of these observatories to identify a set of temporal patterns that can inform the detector operators in day-to-day monitoring and diagnostics. On a day-to-day basis, the operators monitor several seismically relevant data streams to diagnose operational instabilities and sources of noise using some simple empirically-determined thresholds. It can be untenable for a human operator to monitor multiple data streams in this manual fashion and thus a distillation of these data-streams into a more human-friendly format is sought. In this paper, we present an end-to-end machine learning pipeline for features-based multivariate time series clustering to achieve this goal and to provide actionable insights to the detector operators by correlating found clusters with events of interest in the detector.

We investigate the detection of boosted dark matter in a two component dark matter model with the hidden gauge $U(1)_D$ symmetry. The model introduces heavy and light fermionic dark matter components, where the heavy dark matter annihilation in the Galactic center boosts the light dark matter. Directional direct detection experiments, such as NEWSdm, are suitable for testing the boosted scenario by focusing on signals from the Galactic center. We found that the expected event rate could reach up to O(10) events per year per 5 kg of the detector material.

Magnetic reconnection, a fundamental plasma process, is pivotal in understanding energy conversion and particle acceleration in astrophysical systems. While extensively studied in two-dimensional (2D) configurations, the dynamics of reconnection in three-dimensional (3D) systems remain under-explored. In this work, we extend the classical tearing mode instability to 3D by introducing a modulation along the otherwise uniform direction in a 2D equilibrium, given by $g(y)$, mimicking a flux tube-like configuration. We perform linear stability analysis (both analytically and numerically) and direct numerical simulations to investigate the effects of three-dimensionality. Our findings reveal that the 3D tearing instability exhibits reduced growth rates compared to 2D by a factor of $\int g(y)^{1/2} dy~/\int dy$, with the dispersion relation maintaining similar scaling characteristics. We show that the modulation introduces spatially varying resistive layer properties, which influence the reconnection dynamics. Remarkably, we find that Sweet-Parker scaling for the reconnection rate persists even in the absence of a guide field.

We investigate the level-crossing phenomenon in two-axion systems, where the mass eigenvalues intersect as the mass of one axion increases with the cooling of the universe. This phenomenon can significantly alter the abundance of axions in the early universe. Our study focuses on its impact on the QCD axion and an axion-like particle, identifying viable regions of axion mass and decay constant that explain the observed dark matter. We demonstrate the equivalence of two different bases for describing the axion system in the existing literature. Furthermore, we derive an improved expression for the adiabatic condition that overcomes limitations in earlier formulations. This new formulation is basis-independent, and we numerically validate its effectiveness. Our analysis reveals specific relations between axion masses and axion-photon couplings within the viable region. These relations could potentially serve as a smoking gun signal for this scenario if confirmed experimentally. We also find that, using the chiral perturbation model, the thermal friction on the QCD axion might be significantly larger than previously estimated. Additionally, we show that a simple model with axion mixing can naturally realize either a heavier or lighter QCD axion.

This paper presents a systematic exploration of exact solutions for electrically charged wormholes, black holes, and black bounces within the hybrid metric-Palatini gravity (HMPG) framework. HMPG combines features of the metric and Palatini formulations of modified gravity, offering a powerful approach to address challenges in General Relativity, particularly those related to cosmic acceleration and dark matter. We examine configurations characterized by a zero scalar potential under spherical symmetry, and present solutions in both the Jordan and Einstein conformal frames. A diverse set of solutions emerges, including traversable wormholes, black holes with extremal horizons, and "black universe" models in which spacetime beyond the horizon leads to an expanding cosmological solution rather than a singularity. Each configuration is categorized according to the properties of the scalar field, with an in-depth analysis of the horizon and throat structures, asymptotic behaviour, and singularities. These findings underscore the versatility of HMPG in capturing complex gravitational phenomena and broadens the scope of the theory, offering a robust framework for modelling gravitational phenomena across a range of astrophysical contexts. Future work will benefit from extending these solutions to include scalar potentials, addressing both early-universe inflation and late-time acceleration, and applying observational data, such as gravitational lensing and gravitational wave measurements.