Papers for Friday, Apr 18 2025

We perform simulations of magnetohydrodynamic accretion onto equal-mass, non-spinning binary black holes in 3+1 full general relativity addressing the effects of orbital eccentricity. We find that binary black holes with non-negligible eccentricity accrete matter with periodicity that matches the binary orbital period, whereas quasi-circular binaries exhibit accretion rate modulation at approximately $\sim 0.7\times$ their binary orbital period. Additionally, we find that the total jet luminosity is modulated at the orbital period for eccentric binaries, while quasi-circular binaries only exhibit long-term modulations. We perform a radiative transfer calculation of the dual jet synchrotron emission and demonstrate that the optically thin synchrotron emission varies on the binary orbital period for eccentric binaries. Moreover, eccentric binaries spend more time in a low state, where the synchrotron emission is minimum, than in a high state, where the synchrotron emission peaks. The quasi-circular binary also exhibits variability in its optically thin synchrotron emission but the exact frequency of variability does not appear robust against different parameters. Our suite of simulations is an essential step towards providing a comprehensive catalog of multi-messenger theoretical models that will enable studies of supermassive binary black holes detectable across the electromagnetic and gravitational wave spectra.

A variety of dust attenuation/extinction curves have been observed in high-redshift galaxies, with mixed results regarding their correlations with global galaxy properties. These variations are likely driven by factors such as intrinsic dust properties, total dust content, and the dust-star geometry. In this work, we explore how the shape of dust attenuation curves-quantified by the UV-optical slope (S) and UV bump strength (B)-correlates with galaxy properties. Our goal is to identify the key physical mechanisms shaping attenuation curves through cosmic time. We build on arXiv:2402.05996, analyzing 173 dusty galaxies at z ~ 2-11.5, with attenuation curves inferred via SED fitting of JWST data using a modified version of BAGPIPES (arXiv:2304.11178). We investigate trends between S, B, and properties inferred from SED fitting: AV, SFR, stellar mass (M*), specific SFR (sSFR), mass-weighted stellar age (a*), ionization parameter (U), and metallicity (Z). For a subset, we also consider oxygen abundance (12 + log(O/H)), derived via Te-based methods. We find that lower AV galaxies tend to have steeper slopes and stronger UV bumps, consistent with radiative transfer predictions involving dust geometry and content. S also correlates with a* and sSFR, suggesting that strong radiation fields in young, bursty galaxies may destroy small grains, flattening the slope. B correlates with 12 + log(O/H), possibly due to metallicity-driven dust composition changes. Overall, attenuation curve shapes appear most strongly linked to: (1) redshift (dust evolution), (2) AV (RT effects), (3) a* or sSFR (radiation field), and (4) oxygen abundance (dust composition). Disentangling these effects requires spatially resolved data and theoretical models including dust evolution.

The IRAS Vela Shell (IVS) is a structure of enhanced FIR emission located towards the Gum Nebula, a prominent region of $\rm H\alpha$ emission in the local Milky Way shaped by various galactic stellar feedback over the past several million years. We constrain the 3D spatial geometry of the IVS using a parsec-resolution 3D dust map and contextualize it within the broader Gum Nebula. Our analysis reveals a dense, bowl-like IVS structure below the Galactic plane, with a more diffuse component above. We obtain a total shell mass of $5.1_{-2.4}^{+2.4}\times 10^{4}\;\rm M_{\odot}$ and, incorporating previous studies on shell expansion, a momentum of $6.0_{-3.4}^{+4.7}\times 10^{5}\;\rm M_{\odot}\;km\; s^{-1}$. We find a spatial correlation between the morphology of the dust-traced IVS and the Gum Nebula's $\rm H\alpha$ emission when projected onto the sky. We quantify contributions of feedback from stellar winds, an expanding HII region, and supernovae to the IVS formation, finding that stellar winds are subdominant. Our momentum analysis shows that both an HII region and supernova feedback could drive the shell's expansion. Using astrometric constraints from Gaia and Hipparcos, we trace back nearby feedback sources and find that the massive stars $\gamma2$ Velorum and $\zeta$ Puppis are currently within the IVS, producing enough ionizing luminosity to form an HII region of comparable size. Alternatively, if the IVS' momentum is primarily driven by supernovae, $1-2$ events would be required. We also identify several young massive clusters that could have hosted supernovae within the past 3 Myr.

We present near-infrared follow-up observations of the International Gravitational Wave Network (IGWN) event S250206dm with the Wide-Field Infrared Transient Explorer (WINTER). WINTER is a near-infrared time-domain survey designed for electromagnetic follow-up of gravitational-wave sources localized to $\leq$300 deg$^{2}$. The instrument's wide field of view (1.2 deg$^2$), dedicated 1-m robotic telescope, and near-infrared coverage (0.9-1.7 microns) are optimized for searching for kilonovae, which are expected to exhibit a relatively long-lived near-infrared component. S250206dm is the only neutron star merger in the fourth observing run (to date) localized to $\leq$300 deg$^{2}$ with a False Alarm Rate below one per year. It has a $55\%$ probability of being a neutron star-black hole (NSBH) merger and a $37\%$ probability of being a binary neutron star (BNS) merger, with a $50\%$ credible region spanning 38 deg$^2$, an estimated distance of 373 Mpc, and an overall false alarm rate of approximately one in 25 years. WINTER covered $43\%$ of the probability area at least once and $35\%$ at least three times. Through automated and human candidate vetting, all transient candidates found in WINTER coverage were rejected as kilonova candidates. Unsurprisingly, given the large estimated distance of 373 Mpc, the WINTER upper limits do not constrain kilonova models. This study highlights the promise of systematic infrared searches and the need for future wider and deeper infrared surveys.

Dark matter (DM) particles can interact with particles of the Standard Model. Although there exist constraints from direct and indirect detection experiments, the dynamical evolution of astrophysical objects could provide a promising probe for these interactions. Obtaining astrophysical predictions is challenging and limited by our ability to simulate scatterings between DM and baryonic particles within N-body and hydrodynamics simulations. We develop a novel scheme that allows simulating these interacting dark matter (IDM) models and accurately accounts for their angular and velocity dependence, as well as the mass ratio between the DM and baryonic scattering partners. To describe DM-baryon interactions, we use an N-body code together with its implementation of smoothed-particle hydrodynamics and meshless finite mass. The interaction itself is realised in a pairwise fashion by creating a virtual scattering partner from the baryonic particle and letting it interact with a DM particle using a scattering routine initially developed for self-interacting dark matter. After the interaction, the virtual particle is rejoined with the baryonic particle, fulfilling energy and momentum conservation. Through several test problems, we demonstrate that we can reproduce their analytic solutions with our IDM scheme. We comment on various numerical aspects and challenges as well as describe the limitations of our numerical scheme. Furthermore, we study the impact of IDM on halo formation with a collapsing overdensity. Overall, it is possible to accurately model IDM within N-body and hydrodynamics simulations, commonly used in astrophysics. In consequence, our scheme allows for making novel predictions and obtaining new constraints of DM-baryon scattering.

Understanding the processes that transform star-forming galaxies into quiescent ones is key to unraveling the role of environment in galaxy evolution. We present measurements of the luminosity functions (LFs) and stellar mass functions (SMFs) of passive red-sequence galaxies in four galaxy clusters at $0.8 < z < 1.3$, selected using deep VLT observations complemented with data from the GCLASS and GOGREEN surveys. We find a significant enhancement in the abundance of faint/low-mass passive galaxies in both the LFs and SMFs of all four clusters compared to the field. This is further evidenced by a shallower low-mass slope in the composite passive cluster SMF, which yields a Schechter parameter $\alpha = -0.54^{+\,0.03}_{-0.03}$, compared to $\alpha = 0.12^{+\,0.01}_{-0.01}$ for the field. Our findings indicate that quenching processes that act in clusters are enhanced compared to the field, suggesting that environmental quenching mechanisms may already be active by $z\sim1$. To reproduce the observed passive cluster SMF, we estimate that $25\pm5\%$ of the star-forming field population that falls into the cluster must have been quenched. Our results largely support traditional quenching models but highlight the need for deeper studies of larger cluster samples to better understand the role of environmental quenching in the distant Universe.

Approximately half of all disc galaxies exhibit appreciable warps in both their stellar and HI discs. The typical warp amplitude is small (a few degrees), becoming noticeable only at the periphery of the galaxy disc. As a result, warps remain a relatively poorly studied phenomenon. In this study, we investigate a large sample of distant edge-on galaxies (approximately 1,000 objects) to examine the frequency and characteristics of stellar disc warps up to a redshift of $z \sim 2$. For the selected galaxies, we utilize HST data from the Cosmic Evolution Survey field and JWST observations from the Cosmic Dawn Center Archive. We measure the properties of disc warps and investigate their evolution as a function of redshift. Our results indicate a potential evolution in the observed frequency of strong S-shaped warps (with an amplitude greater than 4$^\circ$) in stellar discs as a function of redshift. At $z \approx 2$, the frequency of strong warps reaches approximately 50%, while at $z \approx 0$, this fraction decreases to around 10-15%. We attribute the observed evolution in the occurrence of strong warps to the changing frequency of galaxy interactions and mergers. If galaxy interactions represent one of the primary mechanisms responsible for the formation of warps, then the prevalence of vertical disc deformations should increase in tandem with the rising interaction and merger rate.

Multi-field inflation can be inherently non-predictive, with the exception of models with strong attractors. In this work, we focus on models with multiple scalar fields that are non-minimally coupled to the space-time Ricci curvature scalar, motivated by the expectation of a rich particle spectrum at high energies. We show that in this family of models, the single- and multi-field predictions for CMB observables are identical, as long as there exists at least one non-minimal coupling with $\xi \gg 1$. We provide simple expressions for the Hubble scale, the number of $e$-folds, and the turn rate for systems with an arbitrary number of fields and explore the statistical properties for different priors.

Supernova 2023ixf is a type IIP supernova that was observed in May 2023 in the spiral galaxy Messier 101. This was the closest supernova observed of the decade, making this an exciting discovery. Combining the observed brightness and duration with theoretical scaling relations, we model the lightcurve of this supernova in order to reveal the properties of the progenitor star at the time of explosion, including its mass, radius, and explosion energy. We simulate these explosions using the stellar evolution and radiation-hydrodynamics codes MESA+STELLA. We find that SN2023ixf is not easily explained with "normal" stellar evolution, and only models with a small mass of H-rich ejecta can fit the lightcurve. We also find that the late time properties of the lightcurve are better fit by a higher initial-mass star with substantial mass loss during its lifetime, as compared to models with lower initial mass and less mass loss.

The next generation of weak gravitational lensing surveys has the potential to place stringent constraints on cosmological parameters. However, their analysis is limited by systematics such as the intrinsic alignments of galaxies, which alter weak lensing convergence and can lead to biases in cosmological parameter estimations. In this work, we investigate the impact of intrinsic alignments on non-Gaussian statistics of the weak lensing field using galaxy shapes derived from the IllustrisTNG hydrodynamical simulation. We create two catalogs of ray-traced convergence maps: one that includes the measured intrinsic shape of each galaxy and another where all galaxies are randomly rotated to eliminate intrinsic alignments. We compare an exhaustive list of weak lensing statistics between the two catalogs, including the shear-shear correlation function, the map-level angular power spectrum, one-point, peak count, minimum distribution functions, and Minkowski functionals. For each statistic, we assess the level of statistical distinguishability between catalogs for a set of future survey angular areas. Our results reveal strong small-scale correlation in the alignment of galaxies and statistically significant boosts in weak lensing convergence in both positive and negative directions for high-significance peaks and minimums, respectively. Weak lensing analyses utilizing non-Gaussian statistics must account for intrinsic alignments to avoid significantly compromised cosmological inferences.

We present the first results from the CAPERS survey, utilizing PRISM observations with the JWST/NIRSpec MSA in the PRIMER-UDS field. With just 14 % of the total planned data volume, we spectroscopically confirm two new bright galaxies ($M_{\rm UV}\sim -20.4$) at redshifts $z = 10.562\pm0.034$ and $z = 11.013\pm0.028$. We examine their physical properties, morphologies, and star formation histories, finding evidence for recent bursting star formation in at least one galaxy thanks to the detection of strong (EW$_0\sim70$ A) H$\gamma$ emission. Combining our findings with previous studies of similarly bright objects at high-$z$, we further assess the role of stochastic star formation processes in shaping early galaxy populations. Our analysis finds that the majority of bright ($M_{\rm UV}\lesssim -20$) spectroscopically-confirmed galaxies at $z>10$ were likely observed during a starburst episode, characterized by a median SFR$_{10}$/SFR$_{100}\sim2$, although with substantial scatter. Our work also finds tentative evidence that $z>10$ galaxies are more preferentially in a bursting phase than similarly bright $z\sim6$ galaxies. We finally discuss the prospects of deeper spectroscopic observations of a statistically significant number of bright galaxies to quantify the true impact of bursting star formation on the evolution of the bright end of the ultraviolet luminosity function at these early epochs.

The intergalactic medium (IGM) comprises all the matter that lies between galaxies. Hosting the vast majority ($\gtrsim 90\%$) of the baryons in the Universe, the IGM is a critical reservoir and probe for cosmology and astrophysics, providing insights into large-scale structure formation and galaxy evolution. In this Chapter, we present an overview of the general properties of the IGM, focusing on their dependence on cosmic environment and cosmic time. Emphasis is given to the basic physical principles that allow us to model the density, temperature, and ionization state of the IGM, supported by results from cosmological hydrodynamical simulations. We also cover the foundational principles of quasar spectroscopy used to probe the IGM in absorption, with a particular focus on HI absorption lines. Finally, we briefly discuss future prospects and complementary observational techniques to enhance our understanding of the IGM.

Rocky planets orbiting M-dwarf stars are prime targets for characterizing terrestrial atmospheres, yet their long-term evolution under intense stellar winds and high-energy radiation remains poorly understood. The Kepler-1649 system, which hosts two terrestrial exoplanets orbiting an M5V star, presents a valuable opportunity to explore atmospheric evolution in the extreme environments characteristic of M-dwarf stellar systems. In this Letter we show that both planets could have retained atmospheres over gigayear timescales. Using a multi-species magnetohydrodynamic model, we simulate atmospheric ion escape driven by stellar winds and extreme ultraviolet radiation from 0.7 to 4.8 Gyrs. The results show that total ion escape rates follow a power-law decline ($\propto \tau^{-1.6}$ for Kepler-1649 b, $\propto \tau^{-1.5}$ for Kepler-1649 c$\,$), with O$^{+}$ dominating atmospheric loss (76.8%-98.7%). The escape rates at 4.8 Gyrs are two orders of magnitude lower than those during the early epochs ($1.9\times10^{27}$ s$^{-1}$ at 0.7 Gyr vs. $3.0\times10^{25}$ s$^{-1}$ at 4.8 Gyrs for planet b$\,$), while planet b consistently exhibits 1.1-1.9$\times$ higher O$^{+}$ escape rates than planet c due to its closer orbit (0.051 AU vs. 0.088 AU). Despite substantial early atmospheric erosion, both planets may still retain significant atmospheres, suggesting the potential for long-term habitability. These findings offer predictive insight into atmospheric retention in M-dwarf systems and inform future JWST observations aimed at refining habitability assessments.

We investigate the unusual emission line luminosity ratios observed in the JADES NIRSpec spectroscopy of GN-z11, which reveal exceptionally strong emission lines and a significant detection of the rarely observed N III] $\lambda1748-1753$Å multiplet. These features suggest an elevated N/O abundance, challenging existing models of stellar populations and nebular emission. To assess whether Wolf-Rayet (WR) stars can account for the observed line ratios, we construct a suite of stellar and nebular models incorporating high-resolution stellar spectral libraries, enabling a more accurate treatment of WR evolution and its influence on the ionising radiation field. We find that the inclusion of WR stars is essential for reproducing the observed position of GN-z11 in the C III]/He II versus C III]/C iv diagnostic plane, resolving discrepancies from previous studies. The model-derived metallicity (0.07$\lesssim$Z/Z$_{\odot}\lesssim$0.15), ionisation parameter ($\log\,U$$\approx$-2) and stellar ages are consistent with the literature estimates. However, our models under-predict the N III/O III] ratio, suggesting that WR stars alone cannot fully explain the nitrogen enrichment. This suggests that additional mechanisms, such as rapid chemical enrichment in a young, metal-poor environment, may be necessary to explain the nitrogen excess. While our models successfully reproduce most observed line ratios, further refinements to the models are needed to fully characterise the stellar populations and the enrichment processes of high-redshift galaxies like GN-z11.

Type-I disk migration can form a chain of planets engaged in first-order mean-motion resonances (MMRs) parked at the disk inner edge. However, while second- or even third-order resonances were deemed unlikely due to their weaker strength, they have been observed in some planetary systems (e.g. TOI-178 bc: 5:3, TOI-1136 ef: 7:5, TRAPPIST-1 bcd: 8:5-5:3). We performed $>6,000$ Type-I simulations of multi-planet systems that mimic the observed {\it Kepler} sample in terms of stellar mass, planet size, multiplicity, and intra-system uniformity over a parameter space encompassing transitional and truncated disks. We found that Type-I migration coupled with a disk inner edge can indeed produce second- and third-order resonances (in a state of libration) in $\sim 10\%$ and 2\% of resonant-chain systems, respectively. Moreover, the fraction of individual resonances in our simulations reproduced that of the observed sample (notably, 5:3 is the most common second-order MMR). The formation of higher-order MMRs favors slower disk migration and a smaller outer planet mass. Higher-order resonances do not have to form with the help of a Laplace-like three-body resonance as was proposed for TRAPPIST-1. Instead, the formation of higher-order resonance is assisted by breaking a pre-existing first-order resonance, which generates small but non-zero initial eccentricities ($e\approx10^{-3}$ to 10$^{-2}$). We predict that 1) librating higher-order resonances have higher equilibrium $e$ ($\sim 0.1$); 2) be more likely found as an isolated pair in an otherwise first-order chain; 3) more likely emerge in the inner pairs of a chain.

Metric radio bursts are often said to be valuable diagnostic tools for studying the near-sun kinematics and energetics of the Interplanetary Coronal Mass Ejections (ICMEs). Radio observations also serve as an indirect tool to estimate the coronal magnetic fields. However, how these estimated coronal magnetic fields are related to the magnetic field strength in the ICME at 1 AU has rarely been explored. We aim to establish a relation between the coronal magnetic fields obtained from the radio observations very close to the Sun and the magnetic field measured at 1 AU when the ICME arrives at the Earth. We performed statistical analysis of all metric type II radio bursts in solar cycles 23 and 24, which were found to be associated with ICMEs. We estimated the coronal magnetic field associated with the corresponding CME near the Sun (middle corona) using a split-band radio technique and compared those with the magnetic fields recorded at 1 AU with in-situ observations. We found that the estimated magnetic fields near the Sun using radio techniques are not well correlated with the magnetic fields measured at 1 AU using in-situ observations. This could be due to the complex evolution of the magnetic field as it propagates through the heliosphere. Our results suggest that while metric radio observations can serve as effective proxies for estimating magnetic fields near the Sun, they may not be as effective close to the Earth. At least, no linear relation could be established using metric radio emissions to estimate the magnetic fields at 1 AU with acceptable error margins.

Classical Cepheids can be used as age indicators due to well-established period-age and period-age-color relations. \citet{Desomma2021} refined these relations by including a metallicity term and different Mass-Luminosity assumptions. In this study, we apply the period-age-metallicity relation for the first time to samples of Classical Cepheids in M31 and M33. For both galaxies, we consider Cepheid coordinates and spatial distributions, along with the metallicity gradients by \citet{Zaritsky1994} and \citet{Magrini2007}, to provide a metallicity estimate for each pulsator. By applying the period-age-metallicity relation, we derive individual ages for each Cepheid. Combining these ages and spatial distributions, we construct detailed age maps for both galaxies. Our analysis confirms a radial age gradient, with younger Cepheids preferentially found toward the galactic centers. In M31, we confirm an outer ring at $\sim 11$ kpc, consistent with previous studies, and identify for the first time an inner ring at $\sim 7$ kpc, possibly associated with star formation episodes. Comparing age gradients at different angles, we find a consistent general trend of ages increasing monotonically with radius. At the same time, we observe smaller-scale differences, particularly in the $90^\circ$-$180^\circ$ quadrant, suggesting asymmetric star formation and possible dynamical influences. In contrast, M33 displays a steeper global age gradient, indicating a higher concentration of young stars toward its center. This study highlights the utility of Cepheids as stellar population tracers, providing insights into the star formation and dynamical evolution of spiral galaxies. Future works will extend this methodology to additional galaxies.

M-dwarfs frequently produce flares, and their associated coronal mass ejections (CMEs) may threaten the habitability of close-in exoplanets. M-dwarf flares sometimes show prominence eruption signatures, observed as blue/red asymmetries in the H$\alpha$ line. In Paper I, we reported four candidates of prominence eruptions, which shows large diversity in their durations and velocities. In this study, we statistically investigate how blue/red asymmetries are related with their flare and starspot properties, using the dataset from 27 H$\alpha$ flares in Paper I and previously reported 8 H$\alpha$ flares on an M-dwarf YZ Canis Minoris. We found that these asymmetry events tend to show larger H$\alpha$ flare energies compared to non-asymmetry events. In particular, 5 out of 6 blue asymmetry events are not associated with white-light flares, whereas all 7 red asymmetry events are associated with white-light flares. Furthermore, their starspot distributions estimated from the TESS light curve show that all prominence eruption candidates occurred when starspots were located on the stellar disk center as well as on the stellar limb. These results suggest that flares with lower heating rates may have a higher association rate with prominence eruptions and/or the possibility that prominence eruptions are more detectable on the limb than on the disk center on M-dwarfs. These results provide significant insights into CMEs that can affect the habitable world around M-dwarfs.

We study the nonlinear evolution of matter overdensities using the spherical collapse model in degenerate higher-order scalar-tensor (DHOST) theories beyond Horndeski, employing the effective field theory (EFT) of dark energy approach. We investigate the impact of the EFT parameters characterising DHOST theories on the formation of large-scale structure. We identify the parameter space in which the collapse of the spherical overdensity is prevented by the scalar field turning imaginary at some moment, which allows us to place constraints on the model parameters. We show how the collapse time and the critical density contrast depend on the EFT parameters. To assess the observational implications, we compute the halo mass function using the Press-Schechter formalism. We find that the number density of halos is suppressed compared to the $\Lambda$CDM model due to ``beyond Horndeski'' effects, upon imposing the stability of linear perturbations.

A recent study using weak gravitational lensing reveals that there are some isolated galaxies having almost flat rotation curves at very large distance from the galactic centres. According to the authors of the study this provides a strong challenge the standard cold dark matter model, since the dark haloes are too small to explain their observations, especially for small stellar masses. In this article, we show that improving their model, the virial radius is larger than their estimates. The NFW rotational curve, and especially the pseudo Isothermal one, are in agreement with their flat rotational curves, especially for the larger baryonic mass bins used by the authors.

Magnetic flux tubes in the presence of background rotational flows are abundant throughout the solar atmosphere and may act as conduits for MHD waves to transport energy throughout the solar atmosphere. Here we investigate the contribution from MHD waves to the Poynting flux in a 3D numerical simulation of a realistic solar atmosphere, modelling a structure resembling a solar vortex tube, using the PLUTO code in the presence of different plasma flow configurations. These simulations feature a closed magnetic loop system where a rotational flow is imposed at one foot-point in addition to photospheric perturbations acting as a wave driver mimicking those of p-modes. We find that a variety of MHD waves exist within the vortex tube, including sausage, kink and torsional Alfvén waves, owing to the photospheric wave driver and the nature of the rotational flow itself. We demonstrate how the visual interpretation of different MHD modes becomes non-trivial when a background rotational flow is present compared to a static flux tube. By conducting a simulation both with and without the rotational plasma flow, we demonstrate how the perturbed Poynting flux increases in the presence of the rotational flow as the waves transport increased magnetic energy. We attribute this increase to the dynamical pressure from the rotational flow increasing the plasma density at the tube boundary, which acts to trap the wave energy more effectively inside the vortex. Moreover, we demonstrate how the Poynting flux is always directed upwards in weakly twisted magnetic flux tubes.

The spot evolution on the Sun and solar-type stars is important for understanding the nature of consequential flaring activity. This study statistically investigates the variance of flare occurrence rate through the time evolution of spots on the Sun and solar-type stars. We have compiled the 28-year catalogs of solar flares and their source sunspots obtained from solar surface observations by NOAA and GOES for the Sun. Also, we combined the cataloged stellar flares with the time evolution of starspots estimated by light curves obtained by the 4-year Kepler mission for solar-type stars. For the obtained 24124 solar flares and 180 stellar flares, we calculate the flare occurrence distribution with respect to $t_\mathrm{flare}-t_\mathrm{max}$, which represents the timing of flare through the spot evolution, where $t_\mathrm{flare}$ is the flare occurrence time, and $t_\mathrm{max}$ is the time when the source spot takes its maximum area. When normalized by the spot lifetime, we found that the flare occurrence distribution for $t_\mathrm{flare}-t_\mathrm{max}$ shows a similar distribution regardless of spot size or flare energy, suggesting that the Sun and the solar-type star share the same physical process in the spot-to-flare activity. On this basis, we propose a formula for the time variation of the flare occurrence rate per spot. Also, the correlation between the temporal variation of flare occurrence rate and the time evolution of spot area and the lack of difference in flare occurrence rate between the emergence and decaying phases provide a milestone for the nature of flare-productive spots.

Quasi-periodic eruptions (QPEs), the repeated outbursts observed in soft X-ray bands, have attracted broad interest, but their physical origin is under debate. One of the popular models, the star-disk collision model, suggests that QPEs can be produced through periodic collisions of an orbiting star with the accretion disk of a central black hole (BH). However, previous tests of the star-disk collision model mainly focus on the timing analysis. Other observed properties, such as peak luminosities $L_{\rm{p}}$, durations $t_{\rm{e}}$, and radiation temperatures $T_{\rm{p}}$ of the eruptions, are not systematically investigated. For a sample of six QPE sources and two QPE-like sources, we test the star-disk collision model by using these observables to derive the constraints on the stellar radius $R_*$. We find that, except for two sources (eRo-QPE3 and eRo-QPE4), the rest of the sample either has no allowed $R_*$ to simultaneously reproduce the observed $L_{\rm{p}}$ and $t_{\rm{e}}$, or the required $R_*$ is too large to avoid being disrupted by the central BH. For the two exceptions, a stellar radius of the order of $1\ R_{\rm{\odot}}$ is necessary to satisfy all the constraints. Another issue with the simplest version of this model is that it predicts $k T_{\rm{p}} \sim 10\ \rm{eV}$, one order of magnitude lower than the observed value.

We study the counts-in-cells reduced skewness $s_3$ for dark matter, halo, and galaxy distributions in both real and redshift space, using the ELEPHANT ($\textit{Extended LEnsing PHysics with ANalytical ray Tracing}$) suite of $N$-body simulations. We compare General Relativity (GR) with two extended (EG) gravity models: $f(R)$ gravity with chameleon screening and the normal-branch Dvali-Gabadadze-Porrati (nDGP) model with Vainshtein screening. We quantify the suppression of $s_3$ by redshift-space distortions (RSD), finding that while small-scale skewness is strongly reduced, the $F5$ model retains a $\sim 4\%$ deviation from GR in galaxy samples, corresponding to a $2\sigma$ significance. We show that the ratio $s_3^{\mathrm{RSD}}/s_3^{\mathrm{real}}$ is approximately independent of the gravity model across tracers and redshifts. Our results demonstrate that real-space predictions can help reliably infer redshift-space skewness in both GR and extended gravity, providing a new tool for testing gravity with current and forthcoming galaxy redshift surveys.

We carry out a numerical calculation of magnetar-powered shock break-outs (SBOs) and supernova (SN) light-curves. In particular, we investigate the impact of gravitational wave (GW) emission by the magnetar central engine on its electromagnetic (EM) counterparts in the ULTRASAT band. Our results show that GW emission by the magnetar has only a minor effect on the SBO light-curve. However, we find that SN light-curves can carry a direct signature of GW emission, which becomes more evident at late times (> 20-30 days).~Our results demonstrate that future ULTRASAT observations will provide crucial insights into the magnetar formation process, and unique information for direct searches of long-transient signals with current and future generation GW detectors. In particular, we estimate a rate of multi-messenger (UV+GW) detections of newly formed magnetars $>$ 1 every two years with ULTRASAT and the Einstein Telescope.

Ram pressure stripping (RPS) has a non-negligible impact on the gas content of cluster galaxies. We use the semi-analytic model GAEA and the hydro-simulation TNG to investigate whether cluster galaxies suffer a strong RPS that is sufficient to remove a significant fraction of their gas during the first pericentric passage. We estimate that a ram pressure of $10^{-10.5}$, $10^{-12} $, $10^{-13.5} $ $g cm^{-1} s^{-2}$ can remove at most $90\%$, $50\%$, and $20\%$ of the cold gas reservoir from low-mass galaxies with $9<\log M_{\star}/{\rm M}_{\odot} <9.5$, assuming the gas can be stripped instantaneously. We then use this information to divide the phase space diagram into `strong', `moderate', `weak', and `no' RPS zones. By tracing the orbit of galaxies since $2.5R_{vir}$, we find in both GAEA and TNG that about half of the galaxies in Virgo-like halos ($\log M_h / M_{\odot} \sim 14 $) did not suffer strong RPS during the first pericentric passage. In Coma-like halos ($\log M_h / M_{\odot} \sim 15$), almost all galaxies have suffered strong RPS during the first pericentric passage, which can remove all gas from low-mass galaxies but is insufficient to significantly reduce the gas content of more massive galaxies. In general, results from TNG and GAEA are consistent, with the RPS being only slightly stronger in TNG than in GAEA. Our findings suggest that most cluster galaxies will maintain a notable fraction of their gas and continue forming stars after the first pericentric passage, except for those with low stellar mass ($\log M_{\star}/{\rm M}_{\odot} <9.5$) in very massive halos ($\log M_{h}/{\rm M}_{\odot} > 15$).

We report the discovery of TOI-3493 b, a sub-Neptune-sized planet on an 8.15-d orbit transiting the bright (V=9.3) G0 star HD 119355 (aka TIC 203377303) initially identified by NASA's TESS space mission. With the aim of confirming the planetary nature of the transit signal detected by TESS and determining the mass of the planet, we performed an intensive Doppler campaign with the HARPS spectrograph, collecting radial velocity measurements. We found that TOI-3493 b lies in a nearly circular orbit and has a mass of 9.0+/-1.2 M_earth and a radius of 3.22+/-0.08 R_earth, implying a bulk density of 1.47+/-0.23 g/cm^3, consistent with a composition comprising a small solid core surrounded by a thick H/He dominated atmosphere.

The astronomy communities are widely recognised as mature communities for their open science practices. However, while their data ecosystems are rather advanced and permit efficient data interoperability, there are still gaps between these ecosystems. Semantic artefacts (e.g., ontologies, thesauri, vocabularies or metadata schemas) are a means to bridge that gap as they allow to semantically described the data and map the underlying concepts. The increasing use of semantic artefacts in astronomy presents challenges in description, selection, evaluation, trust, and mappings. The landscape remains fragmented, with semantic artefacts scattered across various registries in diverse formats and structures -- not yet fully developed or encoded with rich semantic web standards like OWL or SKOS -- and often with overlapping scopes. Enhancing data semantic interoperability requires common platforms to catalog, align, and facilitate the sharing of FAIR semantic artefacts. In the frame of the FAIR-IMPACT project, we prototyped a semantic artefact catalogue for astronomy, heliophysics and planetary sciences. This exercise resulted in improved vocabulary and ontology management in the communities, and is now paving the way for better interdisciplinary data discovery and reuse. This article presents current practices in our discipline, reviews candidate SAs for such a catalogue, presents driving use cases and the perspective of a real production service for the astronomy community based on the OntoPortal technology, that will be called OntoPortal-Astro.

The European Space Agency's Ariel mission, scheduled for launch in 2029, aims to conduct the first large-scale survey of atmospheric spectra of transiting exoplanets. Ariel achieves the high photometric stability on transit timescales required to detect the spectroscopic signatures of chemical elements with a payload design optimized for transit photometry that either eliminates known systematics or allows for their removal during data processing without significantly degrading or biasing the detection. Jitter in the spacecraft's line of sight is a source of disturbance when measuring the spectra of exoplanet atmospheres. We describe an improved algorithm for de-jittering Ariel observations simulated in the time domain. We opt for an approach based on the spatial information on the Point Spread Function (PSF) distortion from jitter to detrend the optical signals. The jitter model is based on representative simulations from Airbus Defence and Space, the prime contractor for the Ariel service module. We investigate the precision and biases of the retrieved atmospheric spectra from the jitter-detrended observations. At long wavelengths, the photometric stability of the Ariel spectrometer is already dominated by photon noise. Our algorithm effectively de-jitters both photometric and spectroscopic data, ensuring that the performance remains photon noise-limited across the entire Ariel spectrum, fully compliant with mission requirements. This work contributes to the development of the data reduction pipeline for Ariel, aligning with its scientific goals, and may also benefit other astronomical telescopes and instrumentation.

The High Energy Cosmic-Radiation Detection (HERD) facility has been proposed as a leading experiment on China's Space Station (CSS). Scheduled for installation around 2027, HERD is expected to operate for at least a decade. The main scientific objectives include indirect detection of dark matter with unprecedented sensitivity, studying the cosmic-ray spectrum and composition up to the knee, and observing all-sky gamma rays with energies above 100 MeV. HERD is designed as a large-acceptance telescope with a unique design aimed at maximizing its efficiency. It comprises a central 3D imaging calorimeter (CALO) made of LYSO crystals, encircled by four complementary subdetectors on its top and four lateral faces: the scintillating fiber tracking detector (FIT), the plastic scintillator detector (PSD), the silicon charge detector (SCD), and a transition radiation detector (TRD) on one lateral side. To fully harness HERD gamma-ray detection capabilities down to 100 MeV, an advanced ultra-low-energy gamma-ray (ULEG) trigger system has been developed. We present an extensive overview of the design, performance, and optimization of the gamma-ray trigger system supported by software simulations and preliminary results from the successful implementation of the HERD prototype at CERN's PS and SPS beam test campaigns in Fall 2023.

We study the AGB population of the galaxy M31, based on available HST and Spitzer data, to characterize the individual sources in terms of mass, metallicity and formation epoch of the progenitors. Particular attention is dedicated to the derivation of the dust production rate of the stars, in the attempt of determining the global current dust production rate of the galaxy, divided between the silicates and the carbonaceous dust contributions. We use results from stellar evolution modelling complemented by the description of the dust formation process in the wind, to be used in a population synthesis approach, based on the star formation history and age-metallicity relationship obtained in previous investigations. The comparison between the results from synthetic modelling and the data available are used for the characterization of AGB stars in M31. We find that the bulk of the AGB population of M31 is composed by low-mass stars of different metallicity formed between 6 Gyr and 14 Gyr ago, with an additional, significant contribution from the progeny of 1.7-2.5Msun stars formed during the secondary peak in the star formation, which occurred between 1 and 2 Gyr ago. The dust production rate of the galaxy is mostly provided by carbon stars, whose contribution is of the order of 4x10^{-4} Msun/yr, completed by silicates production from massive AGB stars, occurring at a rate of 6x10^{-5} Msun/yr. The implications of the present results on the reliability of AGB modelling are also commented.

The search for signs of life in the Universe has entered a new phase with the advent of the James Webb Space Telescope (JWST). Detecting biosignature gases via exoplanet atmosphere transmission spectroscopy is in principle within JWST's reach. We reflect on JWST's early results in the context of the potential search for biological activity on exoplanets. The results confront us with a complex reality. Established inverse methods to interpret observed spectra-already known to be highly averaged representations of intricate 3D atmospheric processes-can lead to disparate interpretations even with JWST's quality of data. Characterizing rocky or sub-Neptune-size exoplanets with JWST is an intricate task, and moves us away from the notion of finding a definitive "silver bullet" biosignature gas. Indeed, JWST results necessitate us to allow "parallel interpretations" that will perhaps not be resolved until the next generation of observatories. Nonetheless, with a handful of habitable-zone planet atmospheres accessible given the anticipated noise floor, JWST may continue to contribute to this journey by designating a planet as biosignature gas candidate. To do this we will need to sufficiently refine our inverse methods and physical models for confidently quantifying specific gas abundances and constraining the atmosphere context. Looking ahead, future telescopes and innovative observational strategies will be essential for the reliable detection of biosignature gases.

This paper introduces two astronomical methods developed through computational simulation to evaluate the historical dating of ancient astronomical sources. The first identifies a 1151-year planetary cycle based on the recurrence of visible configurations of Mercury to Saturn, including the Sun and Moon, from a geocentric perspective. The second, called SESCC (Speed-Error Signals Cross Correlation), statistically estimates the epoch of star catalogs by analyzing the correlation between positional error and proper motion in ecliptic latitude. Both methods are reproducible, data-driven, and yield results that contradict key tenets of the New Chronology proposed by Fomenko and Nosovsky, most notably the claim that the Anno Domini began in 1152 CE. Open-source code and analysis tools are provided for independent verification.

The study of molecular line emission is crucial to unveil the kinematics and the physical conditions of gas in star-forming regions. Our aim is to quantify the reliability of using individual molecular transitions to derive physical properties of the bulk of the H2 gas, looking at morphological correlations in their overall integrated molecular line emission with the cold dust. For this study we selected transitions of H2CO, CH$_3$OH, DCN, HC$_3$N, CH$_3$CN, CH$_3$OCHO, SO, and SiO and compared them with the 1.38 mm dust continuum emission at different spatial scales in the ALMAGAL sample, that observed a total of 1013 targets covering all evolutionary stages of the high-mass star-formation process and different conditions of clump fragmentation. We used the method of the histogram of oriented gradients (HOG) implemented in the tool astroHOG to compare the morphology of integrated line emission with maps of the 1.38 mm dust continuum emission. Moreover, we calculated the Spearman's correlation coefficient, and compared it with our astroHOG results. Only H$_2$CO, CH$_3$OH, and SO show emission on spatial scales comparable with the diffuse continuum emission. However, from the HOG method, the median correlation of the emission of each of these species with the continuum is only $\sim$24-29%. In comparison with the dense fragments these molecular species still have low values of correlation. On the other hand DCN, HC$_3$N, CH$_3$CN, and CH$_3$OCHO show a good correlation with the dense dust fragments, above 60%. The worst correlation is seen with SiO, both with the extended continuum emission and with compact sources. From the comparison of the results of the HOG method and the Spearman's correlation coefficient, the HOG method gives much more reliable results than the intensity-based coefficient in estimating the level of similarity of the emission morphology.

Accurately estimating the statistical properties of noise is important in space-based gravitational wave data analysis. Traditional methods often assume uncorrelated noise or impose restrictive parametric forms on cross-channel correlations, which could lead to biased estimation in complex instrumental noise. This paper introduces a spline-based framework with trans-dimensional Bayesian inference to reconstruct the full noise covariance matrix, including frequency-dependent auto- and cross-power spectral densities, without prior assumptions on noise shapes. The developed software $\mathtt{NOISAR}$ can recover the features of the noise power spectrum curves with a relative error $\leq 10\%$ for both auto- and cross-one.

Self-interacting dark matter (SIDM) theories predict that dark matter halos experience core-collapse in late-stage evolution, a process where the halo's inner region rapidly increases in density and decreases in size. This process can be modeled by treating the dark matter as a gravothermal fluid, and solving the fluid equations to predict the density profile evolution. This model is incomplete without calibration to N-body simulations, through a constant factor $\beta$ included in the thermal conductivity for the long-mean-free-path limit. The value of $\beta$ employed in the gravothermal fluid formalism has varied between studies, with no clear universal value in the literature. In this work, we use the N-body code Arepo to conduct a series of isolated core-collapse simulations across a range of scattering cross-sections, halo concentrations, and halo masses to calibrate the heat transfer parameter $\beta$. We find that $\beta$ is independent of cross-section, halo concentration, and halo mass for velocity independent elastic scattering cross-sections. We present a model for an effective $\beta$ as a function of a dimensionless cross-section, to describe halo evolution in the long mean free path limit, and show that it accurately captures halo evolution as long as the cross section is not too large. This effective model facilitates comparisons between simulations and the gravothermal model, and enables fast predictions of the dark matter density profile at any given time without running N-body simulations.

Nanoflares, which are consequences of braids in tangled magnetic fields, are an important candidate to heat the solar corona to million degrees. However, their observational evidence is sparse and many of their observational characteristics are yet to be discovered. With the high-resolution observations taken by the Extreme Ultraviolet Imager onboard the Solar Orbiter, here we study a series of ejections of plasma blobs resulted from a braided magnetic loops in the upper transition region and reveal some critical characteristics of such processes. The cores of these ejections have a size of about 700\,km, a duration less than 1 minute and a speed of about 90\,\kms. An important characteristic is that these plasma blobs are apparently constrained by the post-reconnection magnetic loops, along which they show an extension of up to about 2\,000\,km. The propagation of unbraiding nodes along the main axis of the tangled loops has a speed of about 45\,\kms. The separation angles between the post-reconnection loops and the main axis of the tangled loops are about 30\degree. The observations from the Atmospheric Imaging Assembly reveal that the braiding loops are upper transition region structures. Based on these observations, the typical magnetic free energy producing a blob is estimated to be about $3.4\times10^{23}$\,erg, well in the nano-flare regime, while the kinematic energy of a blob is about $2.3\times10^{23}$\,erg, suggesting that a majority of magnetic free energy in a magnetic braid is likely transferred into kinematic energy.

Context: Recent JWST measurements allow access to the near-infrared spectrum of the sub-Neptune TOI-270 d, for which two different interpretations, a high-metallicity miscible envelope and a lower metallicity hycean world, are currently in conflict. Aims: Here, we reanalyze the published data and reproduce previously retrieved molecular abundances based on an independent data reduction and a different retrieval framework. The aim of this study is to refine the understanding of TOI-270 d and highlight considerations for JWST data analysis. Additionally, we test the impact of data resolution on atmospheric retrieval calculations. Methods: We reduce one JWST NIRSpec G395H and one NIRISS SOSS GR700XD transit dataset using the Eureka! pipeline and a custom MCMC-based light curve fitting algorithm at the instruments' native resolutions. The atmospheric composition is estimated with the updated BeAR retrieval code across a grid of retrieval setups and spectral resolutions. Results: Our transit spectrum is consistent with previous studies, except at the red end of the NIRISS data. Our retrievals support a higher mean molecular weight atmosphere for TOI-270 d. We provide refined abundance constraints and find statistically favored model extensions indicating either sulfur-rich chemistry with species such as CS2, CS, and H2CS, or the possible presence of CH3Cl or CH3F. However, Bayesian inference cannot distinguish between these scenarios due to similar opacities below 4 microns. Conclusions: Our analysis reinforces TOI-270 d as a highly interesting warm sub-Neptune for atmospheric studies, with a complex chemistry in a cloud-free upper atmosphere. However, its exact nature remains uncertain and warrants further detailed photochemical modeling and observations.

In single dish neutral hydrogen (HI) intensity mapping, signal separation methods such as Principal Component Analysis (PCA) are used to clean the astrophysical foregrounds. PCA induces a signal loss in the estimated power spectrum, which can be corrected by a transfer function (TF). By injecting mock signals of HI into the data and performing the PCA cleaning, we can use the cleaned mock HI signal to cross-correlate with the original mock, and estimate the signal loss as a TF, $T(\vec{k})$. As expected, a correction of ${T}(\vec{k})^{-1}$ restores the cross-power between the HI and optical galaxies. However, contrary to intuition, the HI auto-power also requires a ${T}(\vec{k})^{-1}$ correction, not ${T}(\vec{k})^{-2}$. The ${T}(\vec{k})^{-1}$ correction is only known empirically through simulations. In this Letter, we show that the ${T}(\vec{k})^{-1}$ correction in auto-power is universal, and can be analytically proven using the quadratic estimator formalism through window function normalisation. The normalisation can also be used to determine the TF correction for any type of linear process. Using the window function, we demonstrate that PCA induces mode-mixing in the power spectrum estimation, which may lead to biases in the model inference.

Adding to the RESOLVE and ECO Gas in Galaxy Groups (G3) initiative, we examine possible drivers of group-integrated HI-to-halo mass ratios ($M_{\rm HI,grp}/M_{\rm halo}$) and group X-ray emission, including group halo mass ($M_{\rm halo}$), virialization as probed by crossing time ($t_{\rm cross}$), presence of active galactic nuclei (AGN), and group-integrated fractional stellar mass growth rate (FSMGR$_{\rm grp}$). G3 groups span $M_{\rm halo}=10^{11-14.5}\,M_\odot$ with comprehensive HI and AGN information, which we combine with X-ray stacking of ROSAT All-Sky data. We detect hot gas emission exceeding AGN and X-ray binary backgrounds confidently for $M_{\rm halo}=10^{12.6-14}\,M_\odot$ and unambiguously for $M_{\rm halo}>10^{14}\,M_\odot$, reflecting an inverse dependence of $M_{\rm\,HI,grp}/M_{\rm halo}$ and hot gas emission on halo mass. At fixed halo mass, $M_{\rm\,HI,grp}/M_{\rm halo}$ transitions to greater spread below $t_{\rm cross}\sim2$ Gyr. Dividing groups across this transition, lower-$t_{\rm cross}$ groups show elevated X-ray emission compared to higher-$t_{\rm cross}$ groups for $M_{\rm halo}>10^{13.3}\,M_\odot$, but this trend reverses for $M_{\rm halo}=10^{12.6-13.3}\,M_\odot$. Additionally, AGN-hosting halos below $M_{\rm halo}\sim10^{12.1}\,M_\odot$ exhibit a broad, $\sim$0.25 dex deep valley in $M_{\rm HI,grp}/M_{\rm halo}$ compared to non-AGN-hosting halos with correspondingly reduced FSMGR$_{\rm grp}$. When diluted by non-AGN-hosting halos, this valley becomes shallower and narrower, falling roughly between $M_{\rm halo}=10^{11.5}\,M_\odot$ and $M_{\rm halo}=10^{12.1}\,M_\odot$ in the overall $M_{\rm\,HI,grp}/M_{\rm\,halo}$ vs. $M_{\rm halo}$ relation. We may also detect a second, less easily interpreted valley at $M_{\rm halo}\sim10^{13}\,M_\odot$. Neither valley matches theoretical predictions of a deeper valley at or above $M_{\rm halo}=10^{12.1}\,M_\odot$.

Electromagnetic waves (EMWs) can be generated by gravitational waves (GWs) within magnetic field via the Gertsenshtein effect. The conversion probability between GWs and EMWs can be enhanced by inhomogeneities in the electron density and magnetic field within the magnetized plasma of both the Milky Way (MW) and the intergalactic medium (IGM) in the expanding universe. Polarized GWs can induce polarized EMWs, and the polarization properties of these EMWs can be altered by Faraday rotation as they propagate through magnetized plasma. Additionally, the polarization intensity of the EMWs may be weakened due to depolarization effects. In this study, we calculate the enhanced GW-EMW conversion in inhomogeneous magnetized plasma during the propagation of GWs through the universe and our galaxy. We analyze the polarization states of the EMWs generated by polarized GWs and discuss the depolarization effects induced by the medium's irregularities, as well as the differential Faraday rotation occurring in multi-layer polarized radiation. Our work provides alternative methods for detecting GWs and exploring their polarization states, and potentially constrain the parameters of the possible GW sources, especially the primordial black hole (PBH), contributing to the advancement of very-high-frequency GW detection and research.

Aims: This paper investigates the coevolution of metals and dust for 173 galaxies at $4.0<z<11.4$ observed with JWST/NIRSpec in the CEERS project. We focus on galaxies with extremely low dust attenuation to understand the physical mechanisms at play. Methods: We developed a new version of the \texttt{CIGALE} code that integrates spectroscopic and photometric data. By statistically comparing observations with modeled spectra, we derive physical parameters to constrain these mechanisms. Results: Our analysis reveals a population of 49 extremely low dust attenuation galaxies (GELDAs), consistent with $A_{FUV} = 0.0$ within $2\sigma$ and $M_{\star} < 10^9 M_\odot$. The stacked spectrum of GELDAs shows a very blue UV slope $\beta_{FUV} = -2.451 \pm 0.066$ and a Balmer decrement H$\alpha$/H$\beta = 2.932 \pm 0.660$, consistent with no dust and Case B recombination with minimal underlying absorption. Notably, GELDAs are more prevalent at $z > 8.8$ (83.3\%) than at lower redshifts (26.3\%), suggesting they could dominate in the early Universe. Using a far-infrared dust spectrum from the ALPINE sample, we study $M_{dust}$ vs. $M_{\star}$ trends. These exhibit upper and lower sequences connected by transitional galaxies. Our comparison with models indicates a critical transition around $M_{\star} \approx 10^{8.5}\,M_\odot$, from dust dominated by stellar sources (SNe and AGB stars) to dust growth via gas accretion. This corresponds to a metallicity of $12 + \log_{10}(O/H) = 7.60$ ($Z/Z_\odot \approx 0.1$), aligning with the point where ISM dust growth matches stellar dust production. The sample has a high gas fraction ($f_{\mathrm{gas}} \gtrsim 0.9$), with no significant gas expulsion, and high surface gas densities. This leads to low star formation efficiencies compared to sub-millimeter galaxies. GELDAs may help explain the observed excess of bright galaxies at $z \gtrsim 9$.

Forming giant planets are accompanied by circumplanetary disks, as indicated by considerations of angular momentum conservation, observations of candidate protoplanets, and the satellite systems of planets in our Solar System. This paper derives surface density distributions for circumplanetary disks during the final stage of evolution when most of the mass is accreted. This approach generalizes previous treatments to include the angular momentum bias for the infalling material, more accurate solutions for the incoming trajectories, corrections to the outer boundary condition of the circumplanetary disk, and the adjustment of newly added material as it becomes incorporated into the Keplerian flow of the pre-existing disk. These generalizations lead to smaller centrifugal radii, higher column density for the surrounding envelopes, and higher disk accretion efficiency. In addition, we explore the consequences of different angular distributions for the incoming material at the outer boundary, with the concentration of the incoming flow varying from polar to isotropic to equatorial. These geometric variations modestly affect the disk surface density, but also lead to substantial modification to the location in the disk where the mass accretion rate changes sign. This paper finds analytic solutions for the orbits, source functions, surface density distributions, and the corresponding disk temperature profiles over the expanded parameter space outlined above.

We introduce a novel method to constrain the Hubble constant ($H_0$) by combining fast radio bursts (FRBs) and their persistent radio sources (PRSs) through the observationally validated Yang relation, $ L_{\nu} \propto | \mathrm{RM} | $, which links PRS luminosity to the rotation measure (RM) of the associated FRB. Using a mock sample of PRSs, we demonstrate that the Yang relation can help to unravel the degeneracies among $H_0$, baryon density parameter $\Omega_b$, and baryon fraction in the intergalactic medium $f_{\mathrm{IGM}}$ in the traditional approach of using dispersion measure only to perform cosmological analyses. Our method employs a two-stage Markov Chain Monte Carlo (MCMC) analysis to constrain $H_0$. Using the available data of six observed PRS systems, we obtain a preliminary constraint of $H_0 = 75 \pm 30~\mathrm{km\,s^{-1}\,Mpc^{-1}}$. We briefly discuss possible refinements of the method by reducing residual degeneracies and systematic uncertainties using future data and physical modeling. Our results indicate that the Yang relation can potentially become a new probe for performing FRB cosmology.

By consistently using the effective field theory of inflationary fluctuations in the decoupling limit, we explicitly prove that the renormalized one-loop power spectrum of the primordial curvature perturbation freezes exactly on scales larger than its sound horizon.

We consider an appealing scenario for the production of purely gravitational dark matter in the background of warm inflation, a mechanism that maintains stable thermal bath during inflation. Through systematic investigation of various gravitational production channels, we reveal distinctive features compared to the standard inflation scenario. Notably, the inflaton annihilation channel in warm inflation exhibits markedly different thermodynamics from the standard inflation paradigm, leading to a suppression on the production of sub-inflaton-mass dark matter. For the production channel of inflationary vacuum fluctuations, we find a correlation of $\rho_\chi\propto m_\chi^{5/2}$ for the conformally coupled dark matter, which expands the feasible range of dark matter mass. Our results also indicates that a minimum temperature threshold of $10^{-6}M_P$ is necessary for warm inflation, which allows adequate dark matter production. With observational constraints, our results provide stringent limits on the mass range of purely gravitational dark matter with sufficient density: $10^{-8}-10^{-2}M_P$ for minimal coupling and $10^{-14}-10^{-2}M_P$ for conformal coupling.

Cosmic birefringence$-$the rotation of the polarization plane of light as it traverses the universe$-$offers a direct observational window into parity-violating physics beyond the Standard Model. In this work, we revisit the anisotropic component of cosmic birefringence, which leads to the generation of $B$-mode polarization in the cosmic microwave background (CMB). Using an exact theoretical treatment beyond the thin last-scattering surface approximation, we constrain the amplitude of anisotropic birefringence with combined polarization data from SPTpol, ACT, POLARBEAR, and BICEP. The joint analysis yields a best-fit amplitude of $A_{\rm CB} = 0.42^{+0.40}_{-0.34} \times 10^{-4}$, consistent with zero within $2\sigma$, and we place a 95\% confidence-level upper bound of $A_{\rm CB} < 1 \times 10^{-4}$. The constraint is not dominated by any single experiment and remains robust under the inclusion of a possible isotropic rotation angle. These results provide leading constraints on anisotropic cosmic birefringence from CMB $B$-mode polarization and illustrate the potential of upcoming experiments to improve sensitivity to parity-violating effects in the early universe.

An important prediction of inflation is the production of a primordial stochastic gravitational wave background. Observing this background is challenging due to the weakness of the signal and the simultaneous presence of an astrophysical background generated by many unresolved late-time sources. One possible way to distinguish between the two is to examine their anisotropies. In this paper we calculate the primordial correlation function of gravitational wave anisotropies in the cosmological background generated by axion inflation, where the inflaton is a pseudo-Nambu-Goldstone boson coupled to gauge fields. In this scenario, tensor modes arise not only from the standard amplification of vacuum fluctuations present in any inflationary model, but also from the inverse decay process of the produced gauge fields. The correlator of gravitational wave anisotropies consists therefore of two main components: the contribution from vacuum tensor modes and the contribution from tensor modes sourced by the gauge fields. Our analysis shows that, while the former, previously studied in the literature, is negligible, the one arising from the sourced tensor modes, normalized by the fractional energy density at interferometer frequencies, can reach values as large as $\mathcal{O}(10^{-1})$. This result shows that axion inflation can generate large anisotropies with the potential to be observed by gravitational wave detectors within a reasonable time frame.

Differential SAR interferometry (DInSAR), by providing displacement time series over coherent objects on the Earth's surface (persistent scatterers), allows to analyze wide areas, identify ground displacements, and study their evolution at large times. In this work we implement an innovative approach that relies exclusively on line-of-sight displacement time series, applicable to cases of correlated persistent-scatterer displacements. We identify the locus of the final positions of the persistent scatterers and automatically calculate the lower bound of the magnitude of the potential three-dimensional displacements. We present the results obtained by using Sentinel-1 data for investigating the ground stability of the hilly village Cazzaso located in the Italian Alps (Friuli Venezia Giulia region) in an area affected by an active landslide. SAR datasets acquired by Sentinel-1 from both ascending and descending orbits were processed using the SPINUA algorithm. Displacement time series were analysed in order to solve phase unwrapping issues and displacement field calculation.

We present a mechanism that allows thermal relic dark matter to annihilate efficiently in the Galactic Halo and in galaxy clusters, but not in the lower-velocity environments of dwarf spheroidal (dSph) galaxies. We realize this within a complete model in which the dark matter consists of two distinct states separated by a small mass splitting. An indirect detection signal is generated only through the coannihilations of these two states, requiring both to be present. In the halo of the Milky Way, the dark matter particles in the lighter state can be excited into the long-lived heavier state through scattering. Once excited, these heavier particles can coannihilate with those in the lighter state, yielding a gamma-ray signal with little or no suppression. By contrast, the dark matter particles in dwarf galaxies do not possess enough kinetic energy to be excited, thereby suppressing the coannihilation rate and corresponding indirect detection signals from those systems. This framework breaks the predictive relationship that ordinarily exists between these respective gamma-ray signals and complicates our ability to interpret the results of indirect detection searches.

We present the most sensitive search to date for light axion-like particles with masses below a micro-eV, using spectropolarimetric data collected from the Lick and Keck Observatories. The conversion of optical photons emitted from the surface of a magnetic white dwarf (MWD) into axions in the strong magnetic field around the star induces a nearly wavelength-independent linear polarization in the observed starlight. We analyze the Stokes parameters $(U, Q, I)$ measured with the Kast spectrograph at the Lick Observatory toward the MWDs SDSS J033320+000720 and ZTF J190132+145807, and with the LRISp-ADC instrument at the Keck Observatory toward ZTF J190132+145807, SDSS J002129+150223, and SDSS J100356+053825 to search for this effect. The data show no evidence of axion-induced linear polarization, and we set world-leading constraints on the axion-photon coupling $|g_{a\gamma\gamma}| \lesssim 1.7 \times 10^{-12} \,\mathrm{GeV}^{-1}$ at the $95\%$ confidence level for masses $m_a \lesssim 2 \times 10^{-7}\,\mathrm{eV}$.

Gravitational-wave (GW) observations of compact binaries have the potential to unlock several remarkable applications in astrophysics, cosmology, and nuclear physics through accurate measurements of the source luminosity distance and inclination. However, these parameters are strongly correlated when performing parameter estimation, which may hamper the enormous potential of GW astronomy. We comprehensively explore this problem by performing Bayesian inference on synthetic data for a network of current and planned second-generation GW detectors, and for the third-generation interferometer Einstein Telescope~(ET). We quantify the role of the network alignment factor, detector sensitivity, and waveform higher-order modes in breaking this degeneracy. We discuss the crucial role of the binary mass ratio: in particular, we find that ET can efficiently remove the error in the distance as long as the compact binary is asymmetric in mass.

We report on the mass measurement of the rapid proton-capture process key nuclide ${}^{84}$Mo and its vicinity, such as ${}^{78}$Y${}^{\rm m}$, ${}^{79}$Y, ${}^{83}$Nb, and ${}^{88}$Ru, using the multi-reflection time-of-flight spectrograph at RIKEN RIBF. For ${}^{78}$Y${}^{\rm m}$, ${}^{84}$Mo, and ${}^{88}$Ru, their masses are experimentally determined for the first time with uncertainties of $\delta m \approx 20~{\rm keV}$. The mass precision of ${}^{79}$Y and ${}^{83}$Nb is improved to 13 keV and 9.6 keV, respectively. The new $\alpha$-separation energy of ${}^{84}$Mo, 1.434(83) MeV, unambiguously rules out the possibility of forming the ZrNb cycle. The X-ray burst simulation with the new masses shows that our measurements effectively remove the large final abundance uncertainties in the $A=80-90$ mass region. The new mass values improve the prediction power for the composition of the nuclear ashes in X-ray bursts, including the production of light $p$-nuclei.

We study a dark energy model composed of a bare negative cosmological constant and a single ultra-light axion, motivated by the string axiverse. Assuming that intelligent observers can exist and observe an accelerating universe, we derive nontrivial constraints on both the axion mass and the bare cosmological constant. The axion mass is bounded from above to avoid fine-tuning of the initial misalignment angle near the hilltop, and from below because extremely light axions would require the bare cosmological constant to be unnaturally close to zero to achieve accelerated expansion. As a result, the anthropically allowed axion mass range typically lies around $m = \mathcal{O}(10)\, H_0$ for a decay constant close to the Planck scale, where $H_0$ is the observed value of the Hubble constant. In this framework, the dark energy equation of state parameter $w_0$ generically deviates from $-1$ by $\mathcal{O}(0.1)$, providing a natural explanation for why $w \ne -1$ may be expected. This outcome is intriguingly consistent with recent DESI hints of time-varying dark energy, and offers a compelling anthropic explanation within the $\Lambda$ + axion framework.

Based on an extended nuclear statistical equilibrium model, we investigate the properties of non-accreted crusts of young and warm neo-neutron stars, i.e., of finite-temperature inhomogeneous dense matter in beta equilibrium. We present two novel results and one known, but frequently ignored property of such matter. The first new feature is the appearance, in the deep inner crust, of an extensive and almost pure $^{14}$He layer that extends up to the density of the transition to homogeneous matter. This layer may challenge the idea of nuclear pasta phases, significantly impact the transport properties and the crust crystallization process. Second, we raise the question of the (in)stability of the inner crust with respect to diffusion of ions (buoyancy) and demonstrate that our crust is stable, in contrast with the predictions of some other models. Finally, we show that subsaturated stellar matter is thermodynamically stable with respect to density fluctuations, which rules out a first-order phase transition between inhomogeneous and homogeneous phases.

Recent results from the Atacama Cosmology Telescope (ACT) indicate a scalar spectral index $n_s \simeq 0.9743$, in excellent agreement with the prediction of linear inflation. However, the corresponding tensor-to-scalar ratio $r \simeq 0.0667$ is in tension with current observational bounds. In this work, we investigate how this tension can be alleviated in the Palatini formulation of gravity. We consider two classes of models based on simple monomial potentials: (i) models with a non-minimal coupling between the inflaton and gravity, and (ii) models including an $\alpha R^2$ term. In the first case, we find that a quadratic potential with a linear non-minimal coupling leads to the linear inflation attractor, with $r$ suppressed as $\xi$ increases. In the second case, we show that a linear potential can yield values of $r$ consistent with observations for sufficiently large $\alpha$. Our results demonstrate that simple monomial models can remain compatible with current observational constraints when embedded in the Palatini framework.

We demonstrate the prompt-delayed signals induced by knockout neutrons from the quasi-elastic scattering in neutrino experiments provides a new avenue for detecting light dark matter. As an illustration, we consider the detection of atmospheric dark matter in the liquid scintillator detectors. The results show that the constraint on the DM-nucleon interaction from KamLAND is approximately one order of magnitude more stringent than those obtained from the elastic nuclear recoil signals in dark matter direct detection experiments. Furthermore, a larger volume neutrino experiment, such as JUNO, is expected to significantly enhance the light dark matter detection sensitivity through the quasi-elastic scattering.