Papers for Wednesday, Apr 23 2025

Motivated by the recent Year-2 data release of the DESI collaboration, we update our results on time-varying dark energy models driven by the Cohen-Kaplan-Nelson bound. The previously found preference of time-dependent dark energy models compared to $\Lambda$CDM is further strengthend by the new data release. For our particular models, we find that this preference increases up to $\approx 2.6\,\sigma$ depending on the used supernova dataset.

In the presence of a weak gravitational wave (GW) background, astrophysical binary systems act as high-quality resonators, with efficient transfer of energy and momentum between the orbit and a harmonic GW leading to potentially detectable orbital perturbations. In this work, we develop and apply a novel modeling and analysis framework that describes the imprints of GWs on binary systems in a fully time-resolved manner to study the sensitivity of lunar laser ranging, satellite laser ranging, and pulsar timing to both resonant and nonresonant GW backgrounds. We demonstrate that optimal data collection, modeling, and analysis lead to projected sensitivities which are orders of magnitude better than previously appreciated possible, opening up a new possibility for probing the physics-rich but notoriously challenging to access $\mu\mathrm{Hz}$ frequency GWs. We also discuss improved prospects for the detection of the stochastic fluctuations of ultra-light dark matter, which may analogously perturb the binary orbits.

The light curves of radioactive transients, such as supernovae and kilonovae, are powered by the decay of radioisotopes, which release high-energy leptons through $\beta^+$ and $\beta^-$ decays. These leptons deposit energy into the expanding ejecta. As the ejecta density decreases during expansion, the plasma becomes collisionless, with particle motion governed by electromagnetic forces. In such environments, strong or turbulent magnetic fields are thought to confine particles, though the origin of these fields and the confinement mechanism have remained unclear. Using fully kinetic particle-in-cell (PIC) simulations, we demonstrate that plasma instabilities can naturally confine high-energy leptons. These leptons generate magnetic fields through plasma streaming instabilities, even in the absence of pre-existing fields. The self-generated magnetic fields slow lepton diffusion, enabling confinement and transferring energy to thermal electrons and ions. Our results naturally explain the positron trapping inferred from late-time observations of thermonuclear and core-collapse supernovae. Furthermore, they suggest potential implications for electron dynamics in the ejecta of kilonovae. We also estimate synchrotron radio luminosities from positrons for Type Ia supernovae and find that such emission could only be detectable with next-generation radio observatories from a Galactic or local-group supernova in an environment without any circumstellar material.

The DESI collaboration, combining their Baryon Acoustic Oscillation (BAO) data with cosmic microwave background (CMB) anisotropy and supernovae data, have found significant indication against $\Lambda$CDM cosmology. The significance of the exclusion of $\Lambda$CDM can be interpreted entirely as significance of the detection of the $\boldsymbol{w_a}$ parameter that measures variation of the dark energy equation of state. We emphasize that DESI's DR2 exclusion of $\Lambda$CDM is quoted in the articles for a combination of BAO and CMB data with each of three different and overlapping supernovae datasets (at 2.8-sigma for Pantheon+, 3.8-sigma for Union3, and 4.2-sigma for DESY5). We show that one can neither choose amongst nor average over these three different significances. We demonstrate how a principled statistical combination yields a combined exclusion significance of 3.1-sigma. Further we argue that, based on available knowledge, and faced with these competing significances, the most secure inference from the DESI DR2 results is the 3.1-sigma level exclusion of $\Lambda$CDM, obtained from combining DESI+CMB alone, while omitting supernovae.

Several enigmatic dusty sources have been detected in the central parsec of the Galactic Center. Among them is X7, located at only $\sim$0.02 pc from the central super-massive black hole, Sagittarius A* (Sgr A*). Recent observations have shown that it is becoming elongated due to the tidal forces of Sgr A*. X7 is expected to be fully disrupted during its pericenter passage around 2035 which might impact the accretion rate of Sgr A*. However, its origin and nature are still unknown. We investigated the tidal interaction of X7 with Sgr A* in order to constrain its origin. We tested the hypothesis that X7 was produced by one of the observed stars with constrained dynamical properties in the vicinity of Sgr A*. We employed a set of test-particle simulations to reproduce the observed structure and dynamics of X7. The initial conditions of the models were obtained by extrapolating the observationally constrained orbits of X7 and the known stars into the past, making it possible to find the time and source of origin by minimizing the three-dimensional separation and velocity difference between them. Our results show that ejecta from the star S33/S0-30, launched in $\sim$1950, can to a large extent, replicate the observed dynamics and structure of X7, provided that it is initially elongated with a velocity gradient across it, and with an initial maximum speed of $\sim$600~km~s$^{-1}$. Our results show that a grazing collision between the star S33/S0-30 and a field object such as a stellar mass black hole or a Jupiter-mass object is a viable scenario to explain the origin of X7. Nevertheless, such encounters are rare based on the observed stellar dynamics within the central parsec.

We obtain constraints in a 12 parameter cosmological model using the recent DESI Data Release (DR) 2 Baryon Acoustic Oscillations (BAO) data, combined with Cosmic Microwave Background (CMB) power spectra (Planck Public Release (PR) 4) and lensing (Planck PR4 + Atacama Cosmology Telescope (ACT) Data Release (DR) 6) data, uncalibrated type Ia Supernovae (SNe) data from Pantheon+ and Dark Energy Survey (DES) Year 5 (DESY5) samples, and Weak Lensing (WL: DES Year 1) data. The cosmological model consists of six $\Lambda$CDM parameters, and additionally, the dynamical dark energy parameters ($w_0$, $w_a$), the sum of neutrino masses ($\sum m_{\nu}$), the effective number of non-photon radiation species ($N_{\textrm{eff}}$), the scaling of the lensing amplitude ($A_{\textrm{lens}}$), and the running of the scalar spectral index ($\alpha_s$). Our major findings are the following: i) With CMB+BAO+DESY5+WL, we obtain the first 2$\sigma$+ detection of a non-zero $\sum m_{\nu} = 0.19^{+0.15}_{-0.18}$ eV (95%). Replacing DESY5 with Pantheon+ still yields a $\sim$1.9$\sigma$ detection. ii) The cosmological constant lies at the edge of the 95% contour with CMB+BAO+Pantheon+, but is excluded at 2$\sigma$+ with DESY5, leaving evidence for dynamical dark energy inconclusive, contrary to claims by DESI collaboration. iii) With CMB+BAO+SNe+WL, $A_{\textrm{lens}} = 1$ is excluded at $>2\sigma$, while it remains consistent with unity without WL data - suggesting for the first time that the existence of lensing anomaly may depend on non-CMB datasets. iv) The Hubble tension persists at 3.6-4.2$\sigma$ with CMB+BAO+SNe; WL data has minimal impact.

This study investigates the open clusters SAI 72 and SAI 75 using Gaia DR3 data, employing the Automated Stellar Cluster Analysis (ASteCA) tool to determine their structural and fundamental properties, including center coordinates, size, age, distance, mass, luminosity, and kinematics. Based on membership probabilities (P >= 50%), we identified 112 and 115 stars as probable members of SAI 72 and SAI 75, respectively. Radial density profile (RDP) analysis yielded cluster radii of 2.35 arcmin for SAI 72 and 2.19 arcmin for SAI 75. The spectral energy distribution (SED) fitting was performed to refine metallicity, distance, and color excess parameters, ensuring consistency within 1 sigma of isochrone-based estimates. Isochrone fitting of the color-magnitude diagram (CMD) suggests ages of 316 Myr and 302 Myr, with corresponding distances of 3160 +/- 80 pc and 3200 +/- 200 pc. We derived their Galactic positions, projected distances (X_sun, Y_sun), and vertical displacements (Z_sun). Mass function analysis estimates cluster masses of 612 +/- 174 solar masses for SAI 72 and 465 +/- 90 solar masses for SAI 75. Kinematic studies indicate that both clusters have reached dynamical equilibrium. The AD diagram method provided convergent point coordinates of (A, D)_o = (97.016 +/- 0.09, 4.573 +/- 0.05) for SAI 72 and (99.677 +/- 0.10, 1.243 +/- 0.09) for SAI 75. Orbital analysis confirms that both clusters follow nearly circular trajectories with low eccentricities and minor variations in apogalactic and perigalactic distances. Furthermore, we determine that SAI 72 and SAI 75 originated beyond the solar circle at R_Birth = 10.825 +/- 0.068 kpc and R_Birth = 9.583 +/- 0.231 kpc, respectively. Their maximum heights above the Galactic plane, Z_max, are 109 +/- 9 pc for SAI 72 and 232 +/- 24 pc for SAI 75, reinforcing their classification as part of the young stellar disc population.

We present an analysis of AGN activity within recently quenched massive galaxies at cosmic noon ($z\sim 2$), using deep Chandra X-ray observations of the Ultra-Deep Survey (UDS) field. Our sample includes over 4000 massive galaxies ($M_\ast > 10^{10.5}$ M$_{\odot}$) in the redshift range $1 < z < 3$, including more than 200 transitionary post-starburst (PSB) systems. We find that X-ray emitting AGN are detected in $6.2 \pm 1.5$ per cent of massive PSBs at these redshifts, a detection rate that lies between those of star-forming and passive galaxies ($8.2 \pm 0.5$ per cent and $5.7 \pm 0.8$ per cent, respectively). A stacking analysis shows that the average X-ray luminosity for PSBs is comparable to older passive galaxies, but a factor of $2.6 \pm 0.3$ below star-forming galaxies of similar redshift and stellar mass. The average X-ray luminosity in all populations appears to trace the star-formation rate, with PSBs showing low levels of AGN activity consistent with their reduced levels of star formation. We conclude that, on average, we see no evidence for excess AGN activity in the post-starburst phase. However, the low levels of AGN activity can be reconciled with the high-velocity outflows observed in many PSBs, assuming the rare X-ray detections represent short-lived bursts of black hole activity, visible $\sim$5 per cent of the time. Thus, X-ray AGN may help to maintain quiescence in massive galaxies at cosmic noon, but the evidence for a direct link to the primary quenching event remains elusive.

This study comprehensively analyzes three open star clusters: SAI 16, SAI 81, and SAI 86 using Gaia DR3 data. Based on the ASteCA code, we determined the most probable star candidates (P >= 50%) and estimated the number of star members of each cluster as 125, 158, and 138, respectively. We estimated the internal structural parameters by fitting the King model to the observed RDPs, including the core, limited, and tidal radii. The isochrone fitting to the color-magnitude diagram provided log(age) values of 9.13 +/- 0.04, 8.10 +/- 0.04, and 8.65 +/- 0.04 and distances of 3790 +/- 94 pc, 3900 +/- 200 pc, and 3120 +/- 30 pc for SAI 16, SAI 81, and SAI 86, respectively. We also calculated their projected distances from the Galactic plane (X_sun, Y_sun) as well as their vertical distances (Z_sun), Galactocentric distances (R_gc), and total masses (M_C) in solar units, which are about 142 +/- 12, 302 +/- 17, and 192 +/- 14 for SAI 16, SAI 81, and SAI 86, respectively. Examining the dynamical relaxation state indicates that all three clusters are dynamically relaxed. By undertaking a kinematic analysis of the cluster data, the space velocity was determined. We calculated the coordinates of the apex point (A_o, D_o) using the AD diagram method along with the derivation of the solar motion parameters (S_sun, l_A, b_A). Through our detailed dynamic orbit analysis, we determined that the three SAI clusters belong to the young stellar disc, confirming their membership within this component of the Galactic structure.

Multiphase gas can be found in many astrophysical environments, such as galactic outflows, stellar wind bubbles, and the circumgalactic medium, where the interplay between turbulence, cooling, and viscosity can significantly influence gas dynamics and star formation processes. We investigate the role of viscosity in modulating turbulence and radiative cooling in turbulent radiative mixing layers (TRMLs). In particular, we aim to determine how different amounts of viscosity affect the Kelvin-Helmholtz instability (KHI), turbulence evolution, and the efficiency of gas mixing and cooling. Using idealized 2D numerical setups, we compute the critical viscosity required to suppress the KHI in shear flows characterized by different density contrasts and Mach numbers. These results are then used in a 3D shear layer setup to explore the impact of viscosity on cooling efficiency and turbulence across different cooling regimes. We find that the critical viscosity follows the expected dependence on overdensity and Mach number. Our viscous TRMLs simulations show different behaviors in the weak and strong cooling regimes. In the weak cooling regime, viscosity has a strong impact, resulting in laminar flows and breaking previously established inviscid relations between cooling and turbulence (albeit leaving the total luminosity unaffected). However, in the strong cooling regime, when cooling timescales are shorter than viscous timescales, key scaling relations in TRMLs remain largely intact. In this regime -- which must hold for gas to remain multiphase -- radiative losses dominate, and the system effectively behaves as non-viscous regardless of the actual level of viscosity. Our findings have direct implications for both the interpretation of observational diagnostics and the development of subgrid models in large-scale simulations.

The dust attenuation of galaxies is highly diverse and closely linked to stellar population properties and the star dust geometry, yet its relationship to galaxy morphology remains poorly understood. We present a study of 141 galaxies ($9<\log(\rm M_{\star}/\rm M_{\odot})<11.5$) at $1.7<z<3.5$ from the Blue Jay survey combining deep JWST/NIRCam imaging and $R\sim1000$ JWST/NIRSpec spectra. Using \texttt{Prospector} to perform a joint analysis of these data with non-parametric star-formation histories and a two-component dust model with flexible attenuation laws, we constrain stellar and nebular properties. We find that the shape and strength of the attenuation law vary systematically with optical dust attenuation ($A_V$), stellar mass, and star formation rate (SFR). $A_V$ correlates strongly with stellar mass for starbursts, star-forming galaxies and quiescent galaxies. The inclusion of morphological information tightens these correlations: attenuation correlates more strongly with stellar mass and SFR surface densities than with the global quantities. The Balmer decrement-derived nebular attenuation for 67 of these galaxies shows consistent trends with the stellar continuum attenuation. We detect a wavelength-dependent size gradient: massive galaxies ($\rm M_{\star}\gtrsim 10^{10}~M_{\odot}$) appear $\sim30\%$ larger in the rest-optical than in the rest-NIR, driven by central dust attenuation that flattens optical light profiles. Lower-mass systems exhibit more diverse size ratios, consistent with either inside-out growth or central starbursts. These results demonstrate that dust attenuation significantly alters observed galaxy structure and highlight the necessity of flexible attenuation models for accurate physical and morphological inference at cosmic noon.

Self-interacting neutrinos provide an intriguing extension to the Standard Model, motivated by both particle physics and cosmology. Recent cosmological analyses suggest a bimodal posterior for the coupling strength $G_{\rm eff}$, favoring either strong or moderate interactions. These interactions modify the scale-dependence of the growth of cosmic structures, leaving distinct imprints on the matter power spectrum at small scales, $k\,>\,0.1\,{\rm Mpc}^{-1}$. For the first time, we explore how the 21-cm power spectrum from the cosmic dawn and the dark ages can constrain the properties of self-interacting, massive neutrinos. The effects of small-scale suppression and enhancement in the matter power spectrum caused by self-interacting neutrinos propagate to the halo mass function, shaping the abundance of small- and intermediate-mass halos. It is precisely these halos that host the galaxies responsible for driving the evolution of the 21-cm signal during the cosmic dawn. We find that HERA at its design sensitivity can improve upon existing constraints on $G_{\rm eff}$ and be sensitive to small values of the coupling, beyond the reach of current and future CMB experiments. Crucially, we find that the combination of HERA and CMB-S4 can break parameter degeneracies, significantly improving the sensitivity to $G_{\rm eff}$ over either experiment alone. Finally, we investigate the prospects of probing neutrino properties with futuristic Lunar interferometers, accessing the astrophysics-free 21-cm power spectrum during the dark ages. The capability of probing small scales of these instruments will allow us to reach a percent-level constraint on the neutrino self-coupling.

Bayesian field-level inference of galaxy clustering guarantees optimal extraction of all cosmological information, provided that the data are correctly described by the forward model employed. The latter is unfortunately never strictly the case. A key question for field-level inference approaches then is where the cosmological information is coming from, and how to ensure that it is robust. In the context of perturbative approaches such as effective field theory, some progress on this question can be made analytically. We derive the parameter posterior given the data for the field-level likelihood given in the effective field theory, marginalized over initial conditions in the zero-noise limit. Particular attention is paid to cutoffs in the theory, the generalization to higher orders, and the error made by an incomplete forward model at a given order. The main finding is that, broadly speaking, an $m$-th order forward model captures the information in $n$-point correlation functions with $n \leqslant m+1$. Thus, by adding more terms to the forward model, field-level inference is made to automatically incorporate higher-order $n$-point functions. Also shown is how the effect of an incomplete forward model (at a given order) on the parameter inference can be estimated.

Context. To date, more than a hundred debris disks have been spatially resolved. Among them, the young system HD 120326 stands out, displaying different disk substructures on both intermediate (30-150 au) and large (150-1000 au) scales. Aims. We present new VLT/SPHERE (1.0-1.8 $\mu$m) and ALMA (1.3 mm) data of the debris disk around HD 120326. By combining them with archival HST/STIS (0.2-1.0 $\mu$m) and archival SPHERE data, we have been able to examine the morphology and photometry of the debris disk, along with its dust properties. Methods. We present the open-access code MoDiSc (Modeling Disks in Scattered light) to model the inner belt jointly using the SPHERE polarized and total intensity observations. Separately, we modeled the ALMA data and the spectral energy distribution (SED). We combined the results of both these analyses with the STIS data to determine the global architecture of HD 120326. Results. For the inner belt, identified as a planetesimal belt, we derived a semi-major axis of 43 au, fractional luminosity of 1.8 x 10-3 , and maximum degree of polarization of 45-57 % at 1.6 $\mu$m. The spectral slope of its reflectance spectrum is red between 1.0 and 1.3 $\mu$m and gray between 1.3 and 1.8 $\mu$m. Additionally, the SPHERE data show that there could be a halo of small particles or a second belt at distances <150 au. Using ALMA, we derived in the continuum (1.3 mm) an integrated flux of 541-581 $\mu$Jy. We did not detect any 12CO emission. At larger separations (>150 au), we highlight a spiral-like feature spanning hundreds of astronomical units in the STIS data. Conclusions. Further data are needed to confirm and better constrain the dust properties and global morphology of HD 120326.

The Milky Way has a large population of dwarf galaxy satellites. Their properties are sensitive to both cosmology and the physical processes underlying galaxy formation, but these properties are still not properly characterized for the entire satellite population. We aim to provide the most accurate systemic dynamical and metallicity properties of the dwarf galaxy Boötes II (Boo II). We use a new spectroscopic sample of 39 stars in the field of Boo II with data from the Fiber Large Array Multi Element Spectrograph (FLAMES) mounted on the Very Large Telescope (VLT). The target selection is based on a combination of broadband photometry, proper motions from Gaia, and the metallicity-sensitive narrow-band photometry from the Pristine survey that is ideal for removing obvious Milky Way contaminants. We found 9 new members, including 5 also in the recent work of Bruce et al. (2023), and the farthest member to date (5.7 half-light radii from Boo II centroid), extending the spectroscopic spatial coverage of this system. Our metallicity measurements based on the Calcium triplet lines leads to the detection of the two first extremely metal-poor stars (EMPS, [Fe/H] < -3.0) in Boo II. Combining this new dataset with literature data refines Boo II's velocity dispersion (5.6km/s), systemic velocity(-126.8 km/s) and shows that it does not show any sign of a significant velocity gradient. We are thus able to confirm the kinematic and metallicity properties of the satellite as well as identify new members for future high-resolution analyses.

Understanding gas flows between galaxies and their surrounding circum-galactic medium (CGM) is crucial to unveil the mechanisms regulating galaxy evolution, especially in the early Universe. However, observations of the CGM around massive galaxies at z>6 remain limited, particularly in the cold gas phase. In this work, we present multi-configuration ALMA observations of [CII]$\lambda158,\mu$m and millimetre continuum emission in the z~6.4 quasar PSOJ183+05, to trace the cold CGM and investigate the presence of outflows. We find clumpy [CII] emission, tracing gas up to a ~6 kpc radius, consistent with the interface region between the interstellar medium (ISM) and CGM. The [CII] kinematics shows a rotating disk and a high-velocity, biconical outflow extending up to 5 kpc. The inferred mass outflow rate is ~930 MSun/yr, among the highest at z>6, and comparable to the star-formation rate. These findings suggest that quasar-driven outflows can rapidly transfer energy and momentum to the CGM, without immediately quenching star formation in the host galaxy ISM. This supports a delayed feedback scenario, in which outflows reshape CGM conditions and regulate future gas accretion over longer timescales. Combining high-resolution and sensitive ALMA data with observations from JWST and MUSE will be crucial to map the CGM across its different phases and build a comprehensive picture of the baryon cycle in the first massive galaxies.

Galaxy mergers are pivotal events in the evolutionary history of galaxies, with their impact believed to be particularly significant in dwarf galaxies. We report the serendipitous identification of an isolated merging dwarf system with a total stellar mass of M$_{\rm \star}$$\sim$10$^{9.7}$M$_{\rm \odot}$, located in the centre of a cosmic void. This system is one of the rare examples, and possibly the first, of merging dwarf galaxy pairs studied within the central region of a cosmic void. Using CAVITY PPAK-IFU data combined with deep optical broadband imaging from the Isaac Newton Telescope, we analysed the kinematics and ionized gas properties of each dwarf galaxy in the system by employing a full spectral fitting technique. The orientation of this merging pair relative to the line of sight allowed us to determine the dynamical mass of each component, showing that both had similar dynamical masses within galactocentric distances of up to 2.9 kpc. While the gas-phase metallicity of both components is consistent with that of star-forming dwarf galaxies, the star formation rates observed in both components exceed those typically reported for equally massive star-forming dwarf galaxies. This indicates that the merger has presumably contributed to enhancing star formation. Furthermore, we found no significant difference in the optical g-r colour of this merging pair compared to other merging dwarf pairs across different environments. While most merging events occur in group-like environments with high galaxy density and the tidal influence of a host halo, and isolated mergers typically involve galaxies with significant mass differences, the identified merging pair does not follow these patterns. We speculate that the global dynamics of the void or past three-body encounters involving components of this pair and a nearby dwarf galaxy might have triggered this merging event.

The trans-neptunian object (58534) 1997 CQ$_{29}$ (a.k.a. Logos) is a resolved wide binary in the dynamically Cold Classical population. With Hubble Space Telescope resolved observations where the primary Logos is well separated from its secondary Zoe it can be established that Logos has a time-variable brightness. Logos' brightness varied by several tenths of a magnitude over a short timescale of hours while the brightness variability of Zoe was on a longer timescale. New unresolved ground-based observations obtained with the Lowell Discovery Telescope and the Magellan-Baade telescope confirm at least one highly variable component in this system. With our ground-based observations and photometric constraints from space-based observations, we suggest that the primary Logos is likely a close/contact binary whose rotational period is 17.43$\pm$0.06 h for a lightcurve amplitude of 0.70$\pm$0.07 mag while Zoe is potentially a (very) slow rotator with an unknown shape. Using the Candela software, we model the Logos-Zoe system and predict its upcoming mutual events season using rotational, physical, and mutual orbit parameters derived in this work or already published. Zoe's shape and rotational period are still uncertain, so we consider various options to better understand Zoe. The upcoming mutual event season for Logos-Zoe starts in 2026 and will last for four years with up to two events per year. Observations of these mutual events will allow us to significantly improve the physical and rotational properties of both Logos and Zoe.

Context. Extreme-ultraviolet (EUV) observations have revealed small-scale transient brightenings that may share common physical mechanisms with larger-scale solar flares. A notable feature of solar and stellar flares is the presence of quasi-periodic pulsations (QPPs), which are considered a common and potentially intrinsic characteristic. Aims. We investigate the properties of QPPs detected in EUV brightenings, which are considered small-scale flares, and compare their statistical properties with those observed in solar and stellar flares. Methods. We extracted integrated light curves of 22,623 EUV brightenings in two quiet Sun regions observed by the Solar Orbiter/Extreme Ultraviolet Imager and identified QPPs in their light curves using Fourier analysis. Results. Approximately 2.7 % of the EUV brightenings exhibited stationary QPPs. The QPP occurrence rate increased with the surface area, lifetime, and peak brightness of the EUV brightenings. The detected QPP periods ranged from approximately 15 to 260 seconds, which is comparable to the periods observed in solar and stellar flares. Consistent with observations of QPPs in solar and stellar flares, no correlation was found between the QPP period and peak brightness. However, unlike the trend observed in solar flares, no correlation was found between the QPP period and lifetime/length scale. Conclusions. The presence of QPPs in EUV brightenings supports the interpretation that these events may be small-scale manifestations of flares, and the absence of period scaling with loop length further suggests that standing waves may not be the primary driver of QPPs in these events.

The cluster environment has been shown to affect the molecular gas content of cluster members, yet a complete understanding of this often subtle effect has been hindered due to a lack of detections over the full parameter space of galaxy star formation rates and stellar masses. Here we stack CO(2-1) spectra of z~1.6 cluster galaxies to explore the average molecular gas fractions of galaxies both at lower mass (log(M/solar mass)~9.6) and further below the Star Forming Main Sequence (SFMS; DeltaMS~ -0.9) than other literature studies; this translates to a 3sigma gas mass limit of ~7x10^9 solar masses for stacked galaxies below the SFMS. We divide our sample of 54 z~1.6 cluster galaxies, derived from the Spitzer Adaptation of the Red-Sequence Cluster Survey, into 9 groupings, for which we recover detections in 8. The average gas content of the full cluster galaxy population is similar to coeval field galaxies matched in stellar mass and star formation rate. However, when further split by CO-undetected and CO-detected, we find that galaxies below the SFMS have statistically different gas fractions from the field scaling relations, spanning deficiencies to enhancements from 2sigma below to 3sigma above the expected field gas fractions, respectively. These differences between z=1.6 cluster and field galaxies below the SFMS are likely due to environmental processes, though further investigation of spatially-resolved properties and more robust field scaling relation calibration in this parameter space are required.

Axion quark nuggets (AQNs) are hypothetical objects with nuclear density that would have formed during the quark-hadron transition and could make up most of the dark matter today. These objects have a mass greater than a few grams and are sub-micrometer in size. They would also help explain the matter-antimatter asymmetry and the similarity between visible and dark components of the universe, i.e. $\Omega_{\text{DM}} \sim \Omega_{\text{visible}}$. These composite objects behave as cold dark matter, interacting with ordinary matter and producing pervasive electromagnetic radiation. This work aims to calculate the FUV electromagnetic signature in a 1 kpc region surrounding the solar system, resulting from the interaction between antimatter AQNs and baryons. To this end, we use the high-resolution hydrodynamic simulation of the Milky Way, FIRE-2 Latter suite, to select solar system-like regions. From the simulated gas and dark matter distributions in these regions, we calculate the FUV background radiation generated by the AQN model. We find that the results are consistent with the FUV excess recently confirmed by the Alice spectrograph aboard New Horizons, which corroborated the FUV excess initially discovered by GALEX a decade ago. We also discuss the potential cosmological implications of our work, which suggest the existence of a new source of FUV radiation in galaxies, linked to the interaction between dark matter and baryons.

Spotted stars in eclipsing binary systems allow us to gather significant information about the stellar surface inhomogeneities that is otherwise impossible from only photometric data. Starspots can be scanned using the eclipse (or transit) mapping technique, which takes advantage of the passage of a companion star (or planet) in front of a spotted giant star in a binary system. Based on the characteristics of their ultra-precise space photometric light curves, we compile a list of eclipsing binaries whose primary component is a spotted subgiant or giant star, with the aim of applying the eclipse mapping technique to them. Eclipsing binaries with giant primaries were selected from Transiting Exoplanet Survey Satellite (TESS) light curves by visual inspection. Spots showing up as bumps during eclipses are modeled with an eclipse mapping technique specialized for two stars, and the number of spots are found with the help of Bayes factors. The full light curves themselves were analyzed with time series spot modeling, and the results of the two approaches were compared. We present a catalog of 29 eclipsing close binaries with active giant components and analyze TIC 235934420, TIC 271892852 and TIC 326257590 from the Continuous Viewing Zones (CVZ) of TESS. Remarkable agreement is found between the starspot temperatures, sizes, and longitudes from the eclipse mapping results and the corresponding full light curve solutions. Spots are always present at the substellar points of the tidally locked binaries. Data from the TESS CVZ allow us to follow the changes of spot patterns on yearly timescales.

Observations of auroras on exoplanets would provide numerous insights into planet-star systems, including potential detections of the planetary magnetic fields, constraints on host-star wind properties, and information on the thermal structures of planets. However, there have not yet been any discoveries of auroras on exoplanets. In this paper, we focus on the search for infrared auroral emission from the molecular ion H$_3^+$, which is common in the atmospheres of solar system planets Jupiter, Saturn, and Uranus. Using Keck/NIRSPEC high-resolution spectroscopy, we search for H$_3^+$ emission from two hot Jupiters, WASP-80b and WASP-69b. We do not see any evidence of emission in the observed spectra when cross-correlating with an H$_3^+$ spectral model or when using an auto-correlation approach to search for any significant features. We therefore place upper limits on the total emission of $5.32 \times 10^{18}$ W for WASP-80b and $1.64 \times 10^{19}$ W for WASP-69b. These upper limits represent the most stringent limits to date and approach the regime of emission suspected from theoretical models.

Any attempt to understand the ubiquitous nature of the magnetic field in the present universe seems to lead us towards its primordial origin. For large-scale magnetic fields, however, their strength and length scale may not necessarily originate from a singular primordial mechanism, namely inflationary magnetogenesis, which has been a popular consideration in the literature. In this paper, we propose a minimal scenario wherein a large-scale magnetic field is generated from the inflationary perturbation without any non-conformal coupling. Due to their origin in the inflationary scalar spectrum, these primordial fields are inherently weak, with their strength suppressed by the small amplitude of scalar fluctuations. We then consider the coupling between this large-scale weak primordial magnetic field and a light axion of mass $<10^{-28}$ eV, which is assumed to be frozen in a misaligned state until the photon decoupling. After the decoupling, when the universe enters into a dark age, the light axion coherently oscillates. By appropriately tuning the axion-photon coupling parameter $\alpha$, we demonstrate that a large-scale magnetic field of sufficient strength can indeed be generated through tachyonic resonance. We further show that the produced magnetic field induces a unique spectrum with multiple peaks of secondary gravitational waves, which the upcoming CMB-S4 can probe through B-mode polarization. The strength can be sufficient enough to violate the PLANCK bound on tensor-to-scalar ratio $r \lesssim 0.036$. Such a violation leads to a constraint on $\alpha \lesssim 80$. With this limiting value of the coupling, we find that present-day magnetic field strength could be as high as $10^{-10}$ Gauss at $\Mpc$ scale, consistent with observation.

A large fraction of the baryon budget at $z<1$ resides in large-scale filaments in the form of diffuse intergalactic gas, and numerous studies have reported a significant correlation between the strength of the absorptions produced by this gas in the spectra of bright background sources, and impact parameter to cosmic filaments intersected by these sightlines. However, a similar relation is harder to determine for the warm-hot phase of the intergalactic gas, since its higher Doppler parameter and significantly lower neutral gas fraction makes this gas difficult to detect in absorption. We use a sample of 13 broad Ly$\alpha$ absorbers (BLAs) detected in the HST/COS spectrum of a single QSO ($z\sim0.27$), whose sightline intersects several inter-cluster axes, to study the relation between BLAs and the large-scale structure of the Universe. Given their Doppler parameters of $b>40$ km s$^{-1}$, BLAs are good tracers of warm-hot intergalactic gas. We use VLT/MUSE and VLT/VIMOS data to infer local overdensities of galaxies at the redshifts of the BLAs, and to assess the potential association of the BLAs with nearby galaxies. We find that out of the 13 BLAs in our sample, four are associated with a strong overdensity of galaxies, and four with tentative overdensities. The remaining five are located at redshifts where we do not identify any excess of galaxies. We find that these overdensities of galaxies at the redshift of BLAs are local, and they vanish when larger cosmic volumes are considered, in terms of a larger velocity offset to the BLA or larger impact parameter to the QSO sightline. Finally, we find a positive correlation between the total hydrogen column densities inferred from the BLAs, and the relative excess of galaxies at the same redshifts, consistent with the picture where warm-hot gas resides deep within the gravitational potential well of cosmic filaments.

Strong gravitational lensing is a powerful tool for probing the nature of dark matter, as lensing signals are sensitive to the dark matter substructure within the lensing galaxy. We present a comparative analysis of strong gravitational lensing signatures generated by dark matter subhalo populations using two different approaches. The first approach models subhalos using an empirical model, while the second employs the Galacticus semi-analytic model of subhalo evolution. To date, only empirical approaches have been practical in the analysis of lensing systems, as incorporating fully physical models was computationally infeasible. To circumvent this, we utilize a generative machine learning algorithm, known as a normalizing flow, to learn and reproduce the subhalo populations generated by Galacticus. We demonstrate that the normalizing flow algorithm accurately reproduces the Galacticus subhalo distribution while significantly reducing computation time compared to direct simulation. Moreover, we find that subhalo populations from Galacticus produce comparable results to the empirical model in replicating observed lensing signals under the fiducial dark matter model. This work highlights the potential of machine learning techniques in accelerating astrophysical simulations and improving model comparisons of dark matter properties.

Inverse Compton scattering by the thermal motions of electrons is believed to produce polarized hard X-rays in active galactic nuclei and black-hole binaries. Meanwhile, plasma within the plunging region of the black hole free falls into the event horizon with a bulk relativistic speed, which could also imprint polarization on up-scattered photons but has not been discussed in detail. To examine this, we computed polarimetric signatures via general relativistic ray-tracing of a toy model consisting of an accreting, geometrically thin plasma with moderate optical depth, falling onto the black hole with a bulk relativistic speed within the plunging region. We show that the maximum spatially unresolved linear polarization could be as large as approximately $7 - 8$ percent when the black hole is viewed near edge-on, while the corresponding resolved linear polarization could be roughly $50$ percent. The large discrepancy between the two is due to 1) dilution from the radiation outside the plunging region and 2) substantial cancellations of the Stokes $Q$ and $U$ fluxes. The resultant polarization contributed by bulk Comptonization could nevertheless exceed that of thermal electron scattering in a Novikov-Thorne disk. Our results thus suggest a new model for imprinting considerable polarization on the electromagnetic observables of accreting black holes. Measurements of X-ray polarization from black-hole binaries and the central black hole of active galactic nuclei could provide direct detection of the plunging region and help constrain plasma properties in the immediate vicinity of the event horizon.

In certain scenarios, the accreted angular momentum of plasma onto a black hole could be low; however, how the accretion dynamics depend on the angular momentum content of the plasma is still not fully understood. We present three-dimensional, general relativistic magnetohydrodynamic simulations of low angular momentum accretion flows around rapidly spinning black holes (with spin $a = +0.9$). The initial condition is a Fishbone-Moncrief (FM) torus threaded by a large amount of poloidal magnetic flux, where the angular velocity is a fraction $f$ of the standard value. For $f = 0$, the accretion flow becomes magnetically arrested and launches relativistic jets but only for a very short duration. After that, free-falling plasma breaks through the magnetic barrier, loading the jet with mass and destroying the jet-disk structure. Meanwhile, magnetic flux is lost via giant, asymmetrical magnetic bubbles that float away from the black hole. The accretion then exits the magnetically arrested state. For $f = 0.1$, the dimensionless magnetic flux threading the black hole oscillates quasi-periodically. The jet-disk structure shows concurrent revival and destruction while the gas outflow efficiency at the event horizon changes accordingly. For $f \geq 0.3$, we find that the dynamical behavior of the system starts to approach that of a standard accreting FM torus. Our results thus suggest that the accreted angular momentum is an important parameter that governs the maintenance of a magnetically arrested flow and launching of relativistic jets around black holes.

The transfer of material between planetary bodies due to impact events is important for understanding planetary evolution, meteoroid impact fluxes, the formation of near-Earth objects (NEOs), and even the provenance of volatile and organic materials at Earth. This study investigates the dynamics and fate of lunar ejecta reaching Earth. We employ the high-accuracy IAS15 integrator withing the REBOUND package to track for 100,000 years the trajectories of 6,000 test particles launched from various lunar latitudes and longitudes. Our model incorporates a realistic velocity distribution for ejecta fragments (tens of meters in size), derived form large lunar cratering events. Our results show that 22.6% of lunar ejecta collide with Earth, following a power-law $C(t) \propto t^{0.315}$ with half of the impacts occurring within ~10,000 years. We also confirm that impact events on the Moon's trailing hemisphere serve as a dominant source of Earth-bound ejecta, consistent with previous studies. Additionally, a small fraction of ejecta remains transiently in near-Earth space, providing evidence that lunar ejecta may contribute to the NEO population. This aligns with recent discoveries of Earth co-orbitals such as Kamo'oalewa (469219, 2016 HO3) and 2024 PT5, both exhibiting spectral properties consistent with lunar material. These findings enhance our understanding of the lunar ejecta flux to Earth, providing insights into the spatial and temporal patterns of this flux and its broader influence on the near-Earth environment.

In this work we examine the baryon acoustic oscillations (BAO) in 2D angular and redshift space $\{\theta, \Delta z\}$, with $\Delta z$ denoting the redshift difference between two given angular shells. We thus work in the context of tomographic analyses of the large scale structure (LSS) where data are sliced in different redshift shells and constraints on Cosmology are extracted from the auto and cross-angular spectra of two different probes, namely the standard galaxy angular density fluctuations (ADF, or 2D clustering), and the galaxy angular redshift fluctuations (ARF). For these two observables we study by first time how the BAO peak arises in the $\{\theta, \Delta z\}$ plane. Despite being a weak feature (particularly for $\Delta z \neq 0$), a Fisher forecast analysis shows that, a priori, most of the information on cosmological and galaxy bias parameters is carried by the BAO features in shell auto- and cross-angular power spectra. The same study shows that a joint probe analysis (ADF+ARF) increases the Fisher determinant associated to cosmological parameters such as $H_0$ or the Dark Energy Chevallier-Polarski-Linder (CPL) parameters $\{w_0,w_a\}$ by at least an order of magnitude. We also study how the Fisher information on cosmological and galaxy bias-related parameters behaves under different redshift shell configurations: including cross-correlations to neighbour shells extending up to $(\Delta z)^{\rm tot}\sim 0.6$ ($(\Delta z)^{\rm tot}\sim 0.4$) for ADF (ARF) is required for Fisher information to converge. At the same time, configurations using narrow shell widths ($\sigma_z \leq 0.02$) preserve the cosmological information associated to peculiar velocities and typically yield Fisher determinants that are about two orders of magnitudes larger than for wider shell ($\sigma_z>0.02$) configurations.

Using JWST near-infrared data of the inner Orion Nebula, \citet{Pearson_McCaughrean_2023} detected 40 binary systems they proposed to be Jupiter-Mass Binary Objects (JuMBOs) -- although their actual nature is still in debate. Only one of the objects, JuMBO\,24, was detected in the radio continuum. Here, we report on new radio continuum (10 GHz) Karl G. Jansky Very Large Array (VLA) detections of the radio counterpart to JuMBO\,24, and on an unsuccessful search for 5 GHz continuum emission with the High Sensitivity Array (HSA). From our new VLA detections and adopting a distance to the region, we set an upper limit of $\simeq 6$~km~s$^{-1}$ to the velocity of the radio source in the plane of the sky. This upper limit favors an origin for this source similar to that of stars, that is, from a stationary contracting core. The nature of the radio emission remains uncertain but the lack of strong variability (all VLA observations are consistent with a steady flux of $\sim$50 $\mu$Jy), of detection on long HSA baseline, and of detectable circular polarization in VLA data do not favor a non-thermal origin.

We present new rotational period estimates for 216 Jupiter Trojans using photometric data from the Zwicky Transient Facility (ZTF), including 80 Trojans with previously unknown periods. Our analysis reveals rotation periods ranging from 4.6 hours to 447.8 hours. These results support the existence of a spin barrier for Trojans larger than 10 km, with periods clustering between 4 and 4.8 hours. This spin barrier is roughly twice as long as that observed for main-belt asteroids, suggesting that Jupiter Trojans have significantly lower bulk densities, likely due to a higher fraction of ices and volatile materials in their composition. We identify three new Trojans with reliable rotation periods near the spin barrier, doubling the number of known Trojans in this critical period range. Using these results, we estimate a mean density of approximately 0.52 g/cm^3 for rubble-pile Trojans. Our findings support the growing evidence that many Trojans are rubble-pile bodies with distinct physical properties compared to main-belt asteroids. Looking forward, we anticipate that data from the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will provide rotational period estimates for several hundred thousand Trojans, down to objects as small as 1 km, enabling a more detailed investigation of their rotational properties and internal structure.

The Black Hole Explorer (BHEX) will be the first sub-mm wavelength Space Very-Long-Baseline Interferometry (VLBI) mission. It targets astronomical imaging with the highest ever spatial resolution to enable detection of the photon ring of a supermassive black hole. BHEX is being proposed for launch in 2031 as a NASA Small Explorers mission. BHEX science goals and mission opportunity require a high precision lightweight spaceborne antenna. A survey of the technology landscape for realizing such an antenna is presented. Technology readiness (TRL) for the antenna is discussed and assessed to be at TRL 5. An update on our technology maturation efforts is provided. Design studies leading to the conceptual design of a metallized carbon fiber reinforced plastic (CFRP) technology based antenna with a mass of only $\dim 50$ kg, incorporating a 3.4 m primary reflector with a surface precision of < 40 $\mu$m to allow efficient operation up to 320 GHz are outlined. Current plans anticipate attaining TRL6 in 2026 for the BHEX antenna. Completed design studies point to a large margin in surface precision which opens up opportunities for applications beyond BHEX, at significantly higher THz frequencies.

Fast radio bursts (FRBs) are luminous radio transients with millisecond duration. For some active repeaters, such as FRBs 20121102A and 20201124A, more than a thousand bursts have been detected by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The waiting time (WT) distributions of both repeaters, defined as the time intervals between adjacent (detected) bursts, exhibit a bimodal structure well-fitted by two log-normal functions. Notably, the time scales of the long-duration WT peaks for both repeaters show a decreasing trend over time. These similar burst features suggest that there may be a common physical mechanism for FRBs~20121102A and 20201124A. In this paper, we {revisit} the neutron star (NS)--white dwarf (WD) binary model with an eccentric orbit to account for the observed changes in the long-duration WT peaks. According to our model, the shortening of the WT peaks corresponds to the orbital period decay of the NS-WD binary. We consider two mass transfer modes, namely, stable and unstable mass transfer, to examine how the orbital period evolves. Our findings reveal distinct evolutionary pathways for the two repeaters: for FRB~20121102A, the NS-WD binary likely undergoes a combination of common envelope (CE) ejection and Roche lobe overflow, whereas for FRB~20201124A the system may experience multiple CE ejections. These findings warrant further validation through follow-up observations.

This study presents a detailed timing analyses of two cataclysmic variables (CVs), [PK2008] HalphaJ115927 and IGR J14091-610, utilizing the optical data from the Transiting Exoplanet Survey Satellite (TESS). Periods of 7.20$\pm$0.02 h, 1161.49$\pm$0.14 s, and 1215.99$\pm$0.15 s are presented for [PK2008] HalphaJ115927, and are interpreted as the probable orbital, spin, and beat periods of the system, respectively. The presence of multiple periodic variations suggests that it likely belongs to the intermediate polar (IP) category of magnetic CVs. Interestingly, [PK2008] HalphaJ115927 exhibits a unique and strong periodic modulation at 5.66$\pm$0.29 d, which may result from the precession of an accretion disc, similar to the IP TV Col. The detection of a spin signal of 576.63$\pm$0.03 s and inferred orbital signal of $\sim$ 15.84 h supports the classification of IGR J14091-610 as an IP. The identification of such a long orbital period adds a new example to the limited population of long-period IPs. The observed dominant signal at the second harmonic of the orbital frequency also suggests ellipsoidal modulation of the secondary in this system. The observed double-peaked spin pulse profile in [PK2008] HalphaJ115927 likely results from two-pole accretion, where both poles contribute to the spin modulation, and their geometry allows equal visibility of both accreting poles. In contrast, IGR J14091-610 exhibits a single-peaked sinusoidal like spin pulse, attributed to the changing visibility of the accretion curtains due to a relatively low dipole inclination. The present observations indicate that accretion in both systems occurs predominantly through a disc.

Baryonic effects created by feedback processes associated with galaxy formation are an important, poorly constrained systematic effect for models of large-scale structure as probed by weak gravitational lensing. Upcoming surveys require fast methods to predict and marginalize over the potential impact of baryons on the total matter power spectrum. Here we use the FLAMINGO cosmological hydrodynamical simulations to test a recent proposal to approximate the matter power spectrum as the sum of the linear matter power spectrum and a constant multiple, $A_{\rm mod}$, of the difference between the linear and non-linear gravity-only power spectra. We show that replacing this constant multiple with a one-parameter family of sigmoid functions of the wavenumber $k$ allows to us match the predictions of simulations with different feedback strengths for $z \leq 1, k < 3~h\cdot{\rm Mpc}^{-1}$, and the different cosmological models in the FLAMINGO suite. The baryonic response predicted by FLAMINGO models that use jet-like AGN feedback instead of the fiducial thermally-driven AGN feedback can also be reproduced, but at the cost of increasing the number of parameters in the sigmoid function from one to three. The assumption that $A_{\rm mod}$ depends only on $k$ breaks down for decaying dark matter models, highlighting the need for more advanced baryon response models when studying cosmological models that deviate strongly from $\Lambda$CDM.

One of the most important discoveries in modern cosmology is cosmic acceleration. However, we find that today's universe could decelerate in the statistically preferred Chevallier-Polarski-Linder (CPL) scenario over the $\Lambda$CDM model by cosmic microwave background, type Ia supernova and DESI's new measurements of baryon acoustic oscillations. Using various datasets, at a beyond $5\,\sigma$ confidence level, we demonstrate that the universe experiences a triple deceleration during its evolution and finally reaches the state of the ``Big Stall", which predicts that: (i) the universe suddenly comes to a halt in the distant future; (ii) its eventual destiny is dominated by dark matter rather than dark energy ; (iii) it ultimately retains an extremely small fraction of dark energy but exerts an extremely large pressure. Our findings profoundly challenge the established understanding of cosmic acceleration and enrich our comprehension of cosmic evolution.

The variability mechanisms from jetted AGNs are still under debate. Here the damped random walk (DRW) model, implemented through Gaussian Processe (GPs), is used to fit the $ZTF$ long-term optical light curves of 1684 $\gamma$-ray emission jetted AGNs. This analysis yields one of the largest samples with characteristic optical variability timescales for jetted AGNs. A single DRW model from GPs can fit the optical light curve of most jetted AGNs well/potentially well, while there are still some jetted AGNs whose light curve can not be fitted well by a single DRW model. After the jet power, proxied by gamma-ray luminosity, is introduced as a new parameter, new relationships among intrinsic variability time scales, black hole mass and jet power are discovered for efficient accretion AGNs ($\tau^{\rm in} \propto M_{\rm BH}^{0.29^{+0.06}_{-0.06}}P_{\rm jet}^{-0.3^{+0.03}_{-0.03}}$ with scatter of approximately 0.09~dex) and for inefficient accretion AGNs ($\tau^{\rm in} \propto M_{\rm BH}^{0.06^{+0.07}_{-0.07}}P_{\rm jet}^{0.37^{+0.11}_{-0.11}}$ with scatter of approximately 0.14~dex), respectively. Our results support that the optical variability of jetted AGNs with efficient accretion may originate within the standard accretion disk at UV emitting radii similar to non-jetted AGNs, and is directly related to the acceleration of shock in the jet and then enhanced through the beaming effect in beamed AGNs. For the jetted AGNs with inefficient accretion, the intrinsic timescale is consistent with the escape timescale of electrons.

X-ray observations can be used to effectively probe the galactic ecosystem, particularly its hot and energetic components. However, existing X-ray studies of nearby star-forming galaxies are limited by insufficient data statistics and a lack of suitable spectral modeling to account for X-ray emission and absorption geometry. We present results from an X-ray spectral study of M51 using 1.3-Ms Chandra data, the most extensive for such a galaxy. This allows the extraction of diffuse X-ray emission spectra from spiral arm phase-dependent regions using a logarithmic spiral coordinate system. A hierarchical Bayesian approach analyzes these spectra, testing models from simple 1-T hot plasma to those including distributed hot plasma and X-ray-absorbing cool gas. We recommend a model fitting the spectra well, featuring a galactic corona with a lognormal temperature distribution and a disk with mixed X-ray emissions and absorption. In this model, only half of the coronal emission is subject to internal absorption. The best-fit absorbing gas column density is roughly twice that inferred from optical extinction of stellar light. The temperature distribution shows a mean temperature of $\sim 0.1$ keV and an average one-dex dispersion that is enhanced on the spiral arms. The corona's radiative cooling might balance the mechanical energy input from stellar feedback. These results highlight the effectiveness of X-ray mapping of the corona and cool gas in spiral galaxies.

We present a novel reinforcement learning (RL) approach for solving the classical 2-level atom non-LTE radiative transfer problem by framing it as a control task in which an RL agent learns a depth-dependent source function $S(\tau)$ that self-consistently satisfies the equation of statistical equilibrium (SE). The agent's policy is optimized entirely via reward-based interactions with a radiative transfer engine, without explicit knowledge of the ground truth. This method bypasses the need for constructing approximate lambda operators ($\Lambda^*$) common in accelerated iterative schemes. Additionally, it requires no extensive precomputed labeled datasets to extract a supervisory signal, and avoids backpropagating gradients through the complex RT solver itself. Finally, we show through experiment that a simple feedforward neural network trained greedily cannot solve for SE, possibly due to the moving target nature of the problem. Our $\Lambda^*-\text{Free}$ method offers potential advantages for complex scenarios (e.g., atmospheres with enhanced velocity fields, multi-dimensional geometries, or complex microphysics) where $\Lambda^*$ construction or solver differentiability is challenging. Additionally, the agent can be incentivized to find more efficient policies by manipulating the discount factor, leading to a reprioritization of immediate rewards. If demonstrated to generalize past its training data, this RL framework could serve as an alternative or accelerated formalism to achieve SE. To the best of our knowledge, this study represents the first application of reinforcement learning in solar physics that directly solves for a fundamental physical constraint.

Recent observations and statistical studies have revealed that a significant fraction of hydrogen-poor superluminous supernovae (SLSNe-I) exhibit light curves that deviate from the smooth evolution predicted by the magnetar-powered model, instead showing one or more bumps after the primary peak. However, the formation mechanisms of these post-peak bumps remain a matter of debate. Furthermore, previous studies employing the magnetar-powered model have typically assumed a fixed magnetic inclination angle and neglected the effects of magnetar precession. However, recent research has shown that the precession of newborn magnetars forming during the collapse of massive stars causes the magnetic inclination angle to evolve over time, thereby influencing magnetic dipole radiation. In this paper, therefore, we incorporate the effects of magnetar precession into the magnetar-powered model to develop the precessing magnetar-powered model. Using this model, we successfully reproduce the multi-band light curves of 6 selected representative SLSNe-I with post-peak bumps. Moreover, the derived model parameters fall within the typical parameter range for SLSNe-I. By combining the precessing magnetars in SLSNe-I and long GRBs, we find that the ellipticity of magnetars is related to the dipole magnetic field strength, which may suggest a common origin for the two phenomena. Our work provides a potential explanation for the origin of post-peak bumps in SLSNe-I and offers evidence for the early precession of newborn magnetars formed in supernova explosions.

In this work, we present EDRIS (French for Distance Estimator for Incomplete Supernova Surveys), a cosmological inference framework tailored to reconstruct unbiased cosmological distances from type Ia supernovae light-curve parameters. This goal is achieved by including data truncation directly in the statistical model which takes care of the standardization of luminosity distances. It allows us to build a single-step distance estimate by maximizing the corresponding likelihood, free from the biases the survey detection limits would introduce otherwise. Moreover, we expect the current worldwide statistics to be multiplied by O(10) in the upcoming years. This provides a new challenge to handle as the cosmological analysis must stay computationally towable. We show that the optimization methods used in EDRIS allow for a reasonable time complexity of O($N^2$) resulting in a very fast inference process (O(10s) for 1500 supernovae).

This paper reports on the detection of a likely explosive outflow in the high-mass star-forming complex G34.26+0.15, adding to the small number (six) of explosive outflows detected so far. ALMA CO(2-1) and SiO(5-4) archival observations reveal multiple outflow streamers from G34.26+0.15, which correlate well with H2 jets identified from Spitzer-IRAC 4.5 um and [4.5]/[3.6] flux ratio maps. These nearly linear outflow streamers originate from a common center within an ultracompact HII region located in the complex. The velocity spread of the outflow streamers ranges from 0 to 120 km/s. The radial velocities of these streamers follow the Hubble-Lemaître velocity law, indicating an explosive nature. From the CO emission, the total outflow mass, momentum, and outflow energy are estimated to be ~264 M_sun, 4.3*10^3 M_sun km/s, and 10^48 erg, respectively. The event triggering the outflow may have occurred about 19,000 years ago and could also be responsible for powering the expanding UC HII region, given the similar dynamical ages and positional coincidence of the UC HII region with the origin of the outflow. The magnetic field lines in the region associated with G34.26+0.15 also appear to align with the direction of the outflow streamers and jets, possibly being dragged by the explosive outflow.

SiO masers from AGB stars exhibit variability in intensity and polarization during a pulsation period. This variability is explained by radiative transfer and magnetic properties of the molecule. To investigate this phenomenon, a 3D maser simulation is employed to study the SiO masers based on Zeeman splitting. We demonstrate that the magnetic field direction affects maser polarization within small tubular domains with isotropic pumping, and yields results that are similar to those obtained from 1D modelling. This work also studies larger clouds with different shapes. We use finite-element domains with internal node distributions to represent the maser-supporting clouds. We calculate solutions for the population inversions in all transitions and at every node. These solutions show that saturation begins near the middle of a domain, moving towards the edges and particularly the ends of long axes, as saturation progresses, influencing polarization. When the observer's view of the domain changes, the plane of linear polarization responds to the projected shape and the projected magnetic field axis. The angle between the observer's line of sight and the magnetic field may cause jumps in the plane of polarization. Therefore, we can conclude that polarization is influenced by both the cloud's major axis orientation and magnetic field direction. We have investigated the possibility of explaining observed polarization plane rotations, apparently within a single cloud, by the mechanism of line-of-sight overlap of two magnetized maser clouds.

The ultra-hot Jupiter (UHJ) TOI-2109b marks the lower edge of the equilibrium temperature gap between 3500 K and 4500 K, an unexplored thermal regime that separates KELT-9b, the hottest planet yet discovered, from all other currently known gas giants. To study the structure of TOI-2109b's atmosphere, we obtained high-resolution emission spectra of both the planetary day- and nightsides with CARMENES and CRIRES$^+$. By applying the cross-correlation technique, we identified the emission signatures of Fe I and CO, as well as a thermal inversion layer in the dayside atmosphere; no significant H$_2$O signal was detected from the dayside. None of the analyzed species were detectable from the nightside atmosphere. We applied a Bayesian retrieval framework that combines high-resolution spectroscopy with photometric measurements to constrain the dayside atmospheric parameters and derive upper limits for the nightside hemisphere. The dayside thermal inversion extends from 3200 K to 4600 K, with an atmospheric metallicity consistent with that of the host star (0.36 dex). Only weak constraints could be placed on the C/O ratio ($>$ 0.15). The retrieved spectral line broadening is consistent with tidally locked rotation, indicating the absence of strong dynamical processes. An upper temperature limit of 2400 K and a maximum atmospheric temperature gradient of 700 K/log bar could be derived for the nightside. Comparison of the retrieved dayside T-p profile with theoretical models, the absence of strong atmospheric dynamics, and significant differences in the thermal constraints between the day- and nightside hemispheres suggest a limited heat transport efficiency across the planetary atmosphere. Overall, our results place TOI-2109b in a transitional regime between the UHJs below the thermal gap, which show both CO and H$_2$O emission lines, and KELT-9b, where molecular features are largely absent.

Cosmic rays are often modeled as charged particles. This allows their non-ballistic propagation in magnetized structures to be captured. In certain situations, a neutral cosmic ray component can arise. For example, cosmic ray neutrons are produced in considerable numbers through hadronic pp and p$\gamma$ interactions. At ultrahigh energies, the decay timescales of these neutrons is dilated, allowing them to traverse distances on the scale of galactic and cosmological structures. Unlike charged cosmic rays, neutrons are not deflected by magnetic fields. They propagate ballistically at the speed of light in straight lines. The presence of a neutral baryonic cosmic ray component formed in galaxies, clusters and cosmological filaments can facilitate the escape and leakage of cosmic rays from magnetic structures that would otherwise confine them. We show that, by allowing confinement breaking, the formation of cosmic-ray neutrons by high-energy hadronic interactions in large scale astrophysical structures can modify the exchange of ultra high-energy particles across magnetic interfaces between galaxies, clusters, cosmological filaments and voids.

In recent years, high-precision high-cadence space photometry has revealed that stochastic low frequency (SLF) variability is common in the light curves of massive stars. We use the data from the Transiting Exoplanet Survey Satellite (TESS) to study and characterize the SLF variability found in a sample of 49 O- and B-type main-sequence stars across six Cygnus OB~associations and one low-metallicity SMC star AV~232. We compare these results to 53 previously studied SLF variables. We adopt two different methods for characterizing the signal. In the first, we follow earlier work and fit a Lorentzian-like profile to the power density spectrum of the residual light curve to derive the amplitude $\alpha_0$, characteristic frequency $\nu_{\rm char}$, and slope $\gamma$ of the variability. In our second model-independent method, we calculate the root-mean-square (RMS) of the photometric variability as well as the frequency at 50\% of the accumulated power spectral density, $\nu_{50\%}$, and the width of the cumulative integrated power density, $w$. For the full sample of 103 SLF variables, we find that $\alpha_0$, $\gamma$, RMS, $\nu_{50\%}$, and $w$ correlate with the spectroscopic luminosity of the stars. Both $\alpha_0$ and RMS appear to increase for more evolved stars whereas $\nu_{\rm char}$ and $\nu_{50\%}$ both decrease. Finally, we compare our results to 2-D and 3-D simulations of subsurface convection, core-generated internal gravity waves, and surface stellar winds, and find good agreement between the observed $\nu_{\rm char}$ of our sample and predictions from sub-surface convection.

We search for possible GeV-TeV gamma-ray imprints of ultrahigh-energy (UHE; $\gtrsim 0.1$ EeV) cosmic ray (CR) acceleration in the large-scale structures surrounding the brightest gamma-ray burst (GRB) explosion, GRB 221009A. Using 1.25 years of post-event Fermi Large Area Telescope (LAT) data, we construct a 1 GeV - 1 TeV test-statistic (TS) map within 15 Mpc of the burst. We identify two peaks in the TS map with TS $\geq 9$. The most significant peak, J1911.8+2044, exhibits gamma-ray emission in pre-burst LAT data. The other peak, J1913.2+1901, coincides with a 664.6 GeV photon recorded $\sim191.9$ days after the GRB trigger and located at about $0.75^{\circ}$ from the GRB localization. The per-photon 95% containment angle for the LAT is about $0.25^{\circ}$ in the 100 GeV - 1 TeV energy range. We explore two possible origins for the $\gamma$-ray emission: (1) UHECRs from GRB 221009A propagating through a magnetized cosmological volume in its vicinity, and (2) UHE or very-high-energy (VHE; $\gtrsim 100$ GeV) $\gamma$-ray emission from GRB 221009A, propagating in the same volume. In both cases, electromagnetic cascade emission is induced in the structured region embedding the burst. If any TS features are related to large-scale imprints induced by cosmic rays, it might be further evidence that GRB 221009A accelerated UHECRs. However, our results show that alternative scenarios without invoking UHECRs cannot be ruled out, and the observed high-energy photon could be unrelated to GRB 221009A.

Neutron stars (NSs) probe the high-density regime of the nuclear equation of state (EOS). However, inferring the EOS from observations of NSs is a computationally challenging task. In this work, we efficiently solve this inverse problem by leveraging differential programming in two ways. First, we enable full Bayesian inference in under one hour of wall time on a GPU by using gradient-based samplers, without requiring pre-trained machine learning emulators. Moreover, we demonstrate efficient scaling to high-dimensional parameter spaces. Second, we introduce a novel gradient-based optimization scheme that recovers the EOS of a given NS mass-radius curve. We demonstrate how our framework can reveal consistencies or tensions between nuclear physics and astrophysics. First, we show how the breakdown density of a metamodel description of the EOS can be determined from NS observations. Second, we demonstrate how degeneracies in EOS modeling using nuclear empirical parameters can influence the inverse problem during gradient-based optimization. Looking ahead, our approach opens up new theoretical studies of the relation between NS properties and the EOS, while effectively tackling the data analysis challenges brought by future detectors.

We present radial density profiles, as traced by luminous galaxies and dark matter particles, for voids in eleven snapshots of the \texttt{TNG300} simulation. The snapshots span 11.65~Gyr of cosmic time, corresponding to the redshift range $0 \le z \le 3$. Using the comoving galaxy fields, voids were identified via a well-tested, watershed transformation-based algorithm. Voids were defined to be underdense regions that are unlikely to have arisen from Poisson noise, resulting in the selection of $\sim100-200$ of the largest underdense regions in each snapshot. At all redshifts, the radial density profiles as traced by both the galaxies and the dark matter resemble inverse top-hat functions. However, details of the functions (particularly the underdensities of the innermost regions and the overdensities of the ridges) evolve considerably more for the dark matter density profiles than for the galaxy density profiles. At all redshifts, a linear relationship between the galaxy and dark matter density profiles exists, and the slope of the relationship is similar to the bias estimates for \texttt{TNG300} snapshots. Lastly, we identify distinct environments in which voids can exist, defining ``void-in-void" and ``void-in-cloud" populations (i.e., voids that reside in larger underdense or overdense regions, respectively) and we investigate ways in which the relative densities of dark matter and galaxies in the interiors and ridges of these structures vary as a function of void environment.

(Aims) Hot Dust Obscured Galaxies (Hot DOGs) are a population of hyper-luminous, heavily obscured quasars. Although nuclear obscurations close to Compton-thick are typical, a fraction show blue UV spectral energy distributions consistent with unobscured quasar activity, albeit two orders of magnitude fainter than expected from their mid-IR luminosity. The origin of the UV emission in these Blue excess Hot DOGs (BHDs) has been linked to scattered light from the central engine. Here we study the properties of the UV emission in the BHD WISE J020446.13-050640.8 (W0204-0506). (Methods) We use imaging polarization observations in the $R_{\rm Special}$ band obtained with the FORS2 instrument at VLT. We compare these data with radiative transfer simulations to constrain the characteristics of the scattering material. (Results) We find a spatially integrated polarization fraction of $24.7\pm 0.7$%, confirming the scattered-light nature of the UV emission of W0204-0506. The source is spatially resolved in the observations and we find a gradient in polarization fraction and angle that is aligned with the extended morphology of the source found in HST/WFC3 imaging. A dusty, conical polar outflow starting at the AGN sublimation radius with a half-opening angle $\lesssim 50~\rm deg$ viewed at an inclination $\gtrsim 45~\rm deg$ can reproduce the observed polarization fraction if the dust is graphite-rich. We find that the gas mass and outflow velocity are consistent with the range of values found for [OIII] outflows through spectroscopy in other Hot DOGs, though it is unclear whether the outflow is energetic enough to affect the long-term evolution of the host galaxy. Our study highlights the unique potential for polarization imaging to study dusty quasar outflows, providing complementary constraints to those obtained through traditional spectroscopic studies.

Determining the composition of an exoplanet atmosphere relies on the presence of detectable spectral features. The strongest spectral features, including DMS, look approximately Gaussian. Here, I perform a suite of Gaussian feature analyses to find any statistically significant spectral features in the recently published MIRI/LRS spectrum of K2-18b (N. Madhusudhan et al. 2025). In N. Madhusudhan et al. 2025, they claim a 3.4-$\sigma$ detection of spectral features compared to a flat line. In 5 out of 6 tests, I find the data preferred a flat line over a Gaussian model, with a $\chi^{2}_{\nu}$ of 1.06. When centering the Gaussian where the absorptions for DMS and DMDS peak, I find ln(B) = 1.21 in favour of the Gaussian model, with a $\chi^{2}_{\nu}$ of 0.99. With only $\sim$2-$\sigma$ in favour of Gaussian features, I conclude no strong statistical evidence for spectral features.

The infrared (IR) to X-ray luminosity ratio (IRX) is an indicator of the role of the dust plays in cooling of hot gas in supernova remnants (SNR). Using the 3D dynamics of gas and interstellar polydisperse dust grains we analyze the evolution of SNR in the inhomogeneous medium. We obtain spatial distributions of the surface brigthness both of the X-ray emission from hot gas inside SNR and the IR emission from the SNR swept-up shell, as well as, the average gas temperature in the SNR, $T_X$. We find that the IRX changes significantly (by a factor of $\sim 3-30$) as a function of impact distance within the SNR and its age. In a low inhomogeneous medium the IRX drops rapidly during the SNR evolution. On the other hand, if large inhomogeneities are present in the medium, the IRX is maintained at higher levels during the late SNR evolution at radiative phase due to replenishment of dust in the hot gas by incompletely destroyed fragments behind the shock front. We show that the onset of the radiative phase determines the evolution of the $T_X - {\rm IRX}$ diagram. We illustrate that decreasing gas metallicity or density leads to high values of temperature and IRX ratio. We discuss how our results can be applied to the observational data to analyse the SNR older than 10 kyr (i.e. when the mass of the swept-up dust in the shell is expected to exceed that produced in the SNR) in the Galaxy and Large Magellanic Cloud.

Long Gamma Ray Bursts (lGRBs) are associated with jets in Type Ic broadline supernovae. The Collapsar model provides a theoretical framework for the jet formation from the core collapse of a massive star in such supernovae. The GRB can only be produced after a successful jet break out from the star. Under this formalism the GRB duration ($t_{\rm{90}}$) has been hypothesized to be the difference between the central engine activity duration ($t_{\rm{eng}}$) and the jet breakout time ($t_{\rm{bo}}$), that is $t_{\rm{90}} = t_{\rm{eng}} - t_{\rm{bo}}$. This disallows $t_{\rm{90}} > t_{\rm{eng}}$ and puts a lower bound on successful lGRB jet central engine duration ($t_{\rm{eng}} > t_{\rm{bo}}$), various numerical simulations have shown otherwise. This study considers a photospheric GRB emission from a relativistic jet punching out of a Wolf-Rayet-like star. We use the bolometric lightcurve generated to calculate the lGRB duration ($t_{\rm{90}}$) for varying engine duration. We find for longer engine duration the lGRB lightcurve reflects the jet profile and $t_{\rm{90}} \approx t_{\rm{eng}}$. While for shorter engine duration, the $t_{\rm{90}}$ has photospheric radius ($R_{\rm{ph}}$) dependence. This can be modeled by a relation, $t_{\rm{90}} = t^{\rm{90}}_{\rm{eng}} + 0.03\left(\frac{R_{\rm{ph}}}{c}\right)$, where c is the speed of light, with a lower bound on $t_{\rm{90}}$ for a successful lGRB. This relation should be most relevant for possible low-luminous lGRBs originating from a collapsar with central engine duration comparable to the jet breakout time.

The CHNOS elemental budgets of rocky planets are crucial for their structure, evolution and potential chemical habitability. It is unclear how the nonlocal disk processes affecting dust in planet-forming disks affect the CHNOS elemental budgets of nascent planets both inside and outside the Solar System. We aim to quantify the coupled effect of dynamical and collisional processes on the initial refractory CHNOS budgets of planetesimals, forming interior to the water ice line for a Solar and non-Solar composition consistent with the star HIP 43393. Methods. We utilize the SHAMPOO code to track the effects of dynamical and collisional processes on 16000 individual dust monomers. Each monomer is here assigned a refractory chemical composition and mineralogy informed by the equilibrium condensation code GGCHEM given the P-T conditions at the initial position of the monomer. Monomers travel embedded in aggregates through a young class I disk, whose structure is calculated with the ProDiMo code. Furthermore, monomers are allowed to undergo dehydration and desulfurization. We find that solid material becomes well-mixed both radially and vertically. For both the Solar and HIP43393 compositions, the solid phase in the disk midplane regions interior to r~0.7AU can become enriched in hydrogen and sulfur by up to 10at% relative to predictions from purely local calculations. This originates from the inward radial transport of hydrated and sulfur-bearing minerals such as lizardite and iron sulfide. Nonlocal disk processing in a young turbulent, massive disk can lead to significant compositional homogenization of the midplane dust and by extension of the initial composition of planetesimals. Planetesimals forming at r<0.7AU may become enriched in hydrated minerals and sulfur, which could result in more widespread aqueous alteration interior to the water iceline compared to planetesimals that emerge...

Type IIn supernovae (SNe) resembling SN 2009ip (09ip-like SNe) originate from the interaction between circumstellar material (CSM) and the ejecta. This subclass not only shares similar observational properties around the maximum, but is commonly characterized by a long duration precursor before its maximum. Investigating the observed properties of the precursor provides constraints on the mass-loss history of the this http URL present observational data of SN 2023vbg, a 09ip-like type IIn SN that displayed unique observational properties compared to other 09ip-like SNe. SN 2023vbg showed a long-duration precursor at $M_g\sim-14$ mag lasting for $\sim100$ days, followed by a bright bump at $M_g\sim-17$ mag at 12-25 days before the maximum. The luminosity of the precursor is similar to those of other 09ip-like SNe, but the bright bump has not been observed in other this http URL reaching the peak luminosity, the light curve exhibited a peculiar smooth this http URL the H$\alpha$ profile displays two velocity components ($\sim 500$ and $3000\ \mathrm{km\ s^{-1}}$), a broad component observed in other 09ip-like SNe was not detected. We suggest that these properties are explained by the difference in the CSM structure as compared to other 09ip-like SNe; SN 2023vbg had an inner denser CSM component, as well as generally smooth CSM density distribution in a more extended scale, than in the others. Such diversity of CSM likely reflects the diversity of pre-SN outbursts, which in turn may mirror the range of evolutionary pathways in the final stages of the progenitors.

QSEBs are small-scale magnetic reconnection events in lower solar atmosphere. Sometimes, they exhibit transition region counterparts, known as UV brightenings. Magnetic field extrapolations suggest that QSEBs can occur at various locations of a fan-spine topology, with UV brightening occurring at null point through a common reconnection process. We aim to understand how complex magnetic configurations like interacting fan-spine topologies can cause small-scale dynamic phenomena in lower atmosphere. QSEBs were detected using k-means clustering on Hbeta observations from Swedish 1-m Solar Telescope (SST). Further, chromospheric inverted-Y-shaped jets were identified in the Hbeta blue wing. Magnetic field topologies were determined through potential field extrapolations from photospheric magnetograms using the Fe I 6173 A line. UV brightenings were detected in IRIS 1400 A SJI. We identify two distinct magnetic configurations associated with QSEBs, UV brightenings, and chromospheric inverted-Y-shaped jets. The first involves a nested fan-spine structure where, due to flux emergence, an inner 3D null forms inside fan surface of an outer 3D null with some overlap. QSEBs occur at two footpoints along the shared fan surface, with UV brightening located near the outer 3D null point. The jet originates close to the two QSEBs and follows the path of high squashing factor Q. We discuss a comparable scenario using a numerical simulation. In second case, two adjacent fan-spine topologies share fan footpoints at a common positive polarity patch, with the QSEB, along with a chromospheric inverted-Y-shaped jet, occurring at the intersection having high Q values. This study demonstrates through observational and modelling support that associated QSEBs, UV brightenings, and chromospheric inverted-Y-shaped jets share a common origin driven by magnetic reconnection between interacting fan-spine topologies.

The detection of the 244 EeV Amaterasu event by the Telescope Array, one of the most energetic ultrahigh-energy cosmic rays (UHECRs; $E\gtrsim0.1$ EeV) observed to date, invites scrutiny of its potential source. We investigate whether the nearby blazar PKS 1717+177, located within $2.5^\circ$ of the reconstructed arrival direction, could explain the event under a proton-primary hypothesis. Using a one-zone jet model, we fit the multi-wavelength spectral energy distribution of the source, incorporating both leptonic and hadronic cascade emission from photohadronic interactions. Our model supports a cosmic-ray origin of the very-high-energy ($E\gtrsim 100$ GeV) $\gamma$-ray flux and predicts a subdominant neutrino flux, an order of magnitude lower than from TXS 0506+056. Under Lorentz invariance violation, protons above a specific energy can propagate over hundreds of Mpc without significant energy loss for certain parameter choices. In such a scenario, our analysis indicates negligible deflection in the Galactic magnetic field, implying a strong extragalactic magnetic field, placing a lower bound on the field strength. Our findings provide a compelling multi-messenger framework linking UHECRs, $\gamma$ rays, and neutrinos and motivate targeted searches by current and future high-energy neutrino telescopes during increased $\gamma$-ray or X-ray activity of this blazar.

Tidal disruption events (TDEs) occur when stars pass close enough to supermassive black holes to be torn apart by tidal forces. Traditionally, these events are studied with computationally intensive hydrodynamical simulations. In this paper, we present a fast, physically motivated two-stage model for TDEs. In the first stage, we model the star's tidal deformation using linear stellar perturbation theory, treating the star as a collection of driven harmonic oscillators. When the tidal energy exceeds a fraction $\gamma$ of the star's gravitational binding energy (with $\gamma \sim \mathcal{O}(1)$), we transition to the second stage, where we model the disrupted material as free particles. The parameter $\gamma$ is determined with a one-time calibration to hydrodynamical simulations. This method enables fast computation of the energy distribution $\mathrm{d}M/\mathrm{d}E$ and fallback rate $\mathrm{d}M/\mathrm{d}T$, while offering physical insight into the disruption process. We apply our model to MESA-generated profiles of middle-age main-sequence stars. Our code is available on GitHub.

Jets and disc winds play an important role in the evolution of protoplanetary discs and the formation of planetary systems. However, there is still a lack of observational data regarding the presence and parameters of outflows, especially for close young binaries. In this study, we aim to find the HH flow near the young sub-arcsecond binary DF Tau and explore its morphology. Narrow-band H$\alpha$ and H$_2$ 2.12 $\mu$m imaging and spectroscopic observations of DF Tau and its vicinity were performed. We have discovered several emission nebulae near the binary, which likely result from the interaction of gas outflow from the binary components with the surrounding medium. The outflow appears to occur both in the form of jets, generating numerous Herbig-Haro objects (HH 1266 flow), and as a weakly collimated wind responsible for the formation of the ring-like nebula around the binary and the rim of the cometary globule. We have found that the angle between the jet and the counter-jet is $168^\circ$ and discuss the complex morphology of the HH flow.

Deconvolution, imaging and calibration of data from radio interferometers is a challenging computational (inverse) problem. The upcoming generation of radio telescopes poses significant challenges to existing, and well proven data reduction pipelines due to the large data sizes expected from these experiments, and the high resolution and dynamic range. In this manuscript, we deal with the deconvolution problem. A variety of multiscalar variants to the classical CLEAN algorithm (the de-facto standard) have been proposed in the past, often outperforming CLEAN at the cost of significantly increasing numerical resources. In this work, we aim to combine some of these ideas for a new algorithm, Autocorr-CLEAN, to accelerate the deconvolution and prepare the data reduction pipelines for the data sizes expected by the upcoming generation of instruments. To this end, we propose to use a cluster of CLEAN components fitted to the autocorrelation function of the residual in a subminor loop, to derive continuously changing, and potentially non-radially symmetric, basis functions for CLEANing the residual. Autocorr-CLEAN allows for the superior reconstruction fidelity achieved by modern multiscalar approaches, and their superior convergence speed. It achieves this without utilizing any substep of super-linear complexity in the minor loops, keeping the single minor loop and subminor loop iterations at an execution time comparable to CLEAN. Combining these advantages, Autocorr-CLEAN is found to be up to a magnitude faster than the classical CLEAN procedure. Autocorr-CLEAN fits well in the algorithmic framework common for radio interferometry, making it relatively straightforward to include in future data reduction pipelines. With its accelerated convergence speed, and smaller residual, Autocorr-CLEAN may be an important asset for the data analysis in the future.

The enhanced primordial scalar power spectrum is a widely studied mechanism for generating primordial gravitational waves (PGWs), also referred to as scalar-induced gravitational waves (SIGWs). This process also plays a pivotal role in facilitating the formation of primordial black holes (PBHs). Traditionally, the ultra slow-roll (USR) mechanism has been the predominant approach used in the early universe. In this framework, the second slow-roll parameter $\epsilon_2$, is typically set to $-6$ or lower for a brief period -- marking a significant departure from the standard slow-roll condition where $\epsilon_2 \simeq 0$. Such conditions often emerge in models with inflection points or localized features, such as bumps in the potential. In this paper, we challenge the conventional assumption that $\epsilon_2 \lesssim -6$ is a prerequisite for substantial amplification of the scalar power spectrum. We demonstrate that any negative value of the second slow-roll parameter can indeed enhance the scalar power spectrum through sub-horizon growth, establishing this as a necessary and sufficient condition for amplification. Consequently, this mechanism facilitates the generation of both PGWs and PBHs. To illustrate this, we examine a standard scenario where a brief USR phase is embedded between two slow-roll (SR) phases. By systematically varying $\epsilon_{2}$ values from $-1$ to $-10$ in the USR region, we investigate the amplification of the power spectrum and its implications for PGWs and PBHs production, particularly in the context of ongoing and future cosmological missions.

We show that the regularization of the second order pole in the pole inflation can induce the increase of $n_s$, which may be important after the latest data release of cosmic microwave background (CMB) observation by Atacama Cosmology Telescope (ACT). Pole inflation is known to provide a unified description of attractor models that they can generate a flat plateau for inflation given a general potential. Recent ACT observation suggests that the constraint on the scalar spectral index $n_s$ at CMB scale may be shifted to a larger value than the predictions in the Starobinsky model, the Higgs inflation, and the $\alpha$-attractor model, which motivates us to consider the modification of the pole inflation. We find that if we regularize the second order pole in the kinetic term such that the kinetic term becomes regular for all field range, we can generally increase $n_s$ because the potential in the large field regime will be lifted. We have explicitly demonstrated that this type of regularized pole inflation can naturally arise from the Einstein-Cartan formalism, and the inflationary predictions are consistent with the latest ACT data without spoiling the success of the $\alpha$-attractor models.

We present estimators for quantifying intrinsic alignments in large spectroscopic surveys that efficiently capture line-of-sight (LOS) information while being relatively insensitive to redshift-space distortions (RSD). We demonstrate that changing the LOS integration range, {\Pi}max, as a function of transverse separation outperforms the conventional choice of a single {\Pi}max value. This is further improved by replacing the flat {\Pi}max cut with a LOS weighting based on shape projection and RSD. Although these estimators incorporate additional LOS information, they are projected correlations that exhibit signal-to-noise ratios comparable to 3D correlation functions, such as the IA quadrupole. Using simulations from Abacus Summit, we evaluate these estimators and provide recommended {\Pi}max values and weights for projected separations of 1 - 100 Mpc/h. These will improve measurements of intrinsic alignments in large cosmological surveys and the constraints they provide for both weak lensing and direct cosmological applications.

The natural environment of the Earth can act as a sensitive detector for dark matter in ultralight axions. When axions with masses between $1\times10^{-15}\,{\rm eV}$ and $1\times10^{-13}\,{\rm eV}$ pass through the Earth, they interact with the global geomagnetic field, generating electromagnetic (EM) waves in the extremely low-frequency range ($0.3$--$30\,{\rm Hz}$) through axion-photon coupling. This paper is one of a series of companion papers for~\cite{Taruya:2025zql}, focusing on the data analysis method and search results for an axion signal. Utilizing the theoretical predictions of axion-induced EM spectra from a companion study, we analyzed long-term observational data of terrestrial magnetic fields in this frequency band to search for axion-induced signals. Our analysis identified 65 persistent signal candidates with a signal-to-noise ratio (SNR) greater than 3. Aside from these candidates, we placed a new upper bound on the axion-photon coupling parameter, significantly refining the previous constraint from CAST by at most two orders of magnitude down to $g_{a\gamma} \lesssim 4\times10^{-13} \,{\rm GeV}^{-1}$ for the axion mass around $3 \times 10^{-14}\,{\rm eV}$.

Neutron stars offer powerful astrophysical laboratories to probe the properties of dark matter. Gradual accumulation of heavy, non-annihilating dark matter in neutron stars can lead to the formation of comparable-mass black holes, and non-detection of gravitational waves from mergers of such low-mass black holes can constrain such dark matter interactions with nucleons. These constraints, though dependent on the currently uncertain binary neutron star merger rate density, are significantly more stringent than those from direct detection experiments and provide some of the strongest limits on heavy, non-annihilating dark matter interactions. Additionally, dark matter with baryon number-violating interactions can induce excess heating in cold neutron stars and is thus significantly constrained by thermal observations of cold neutron stars.

We examine the impact of non-perturbative quantum corrections to the entropy of both charged and charged rotating quasi-topological black holes, with a focus on their thermodynamic properties. The negative-valued correction to the entropy for small black holes is found to be unphysical. Furthermore, we analyze the effect of these non-perturbative corrections on other thermodynamic quantities, including internal energy, Gibbs free energy, charge density, and mass density, for both types of black holes. Our findings indicate that the sign of the correction parameter plays a crucial role at small horizon radii. Additionally, we assess the stability and phase transitions of these black holes in the presence of non-perturbative corrections. Below the critical point, both the corrected and uncorrected specific heat per unit volume are in an unstable regime. This instability leads to a first-order phase transition, wherein the specific heat transitions from negative to positive values as the system reaches a stable state.

We analyze the damping of inflationary gravitational waves (GW) that re-enter the Hubble horizon before or during a post-inflationary era dominated by a meta-stable, right-handed neutrino (RHN), whose out-of-equilibrium decay releases entropy. Within a minimal type-I seesaw extension of the Standard Model (SM), we explore the conditions under which the population of thermally produced RHNs remain long-lived and cause a period of matter-domination. We find that the suppression of the GW spectrum occurs above a characteristic frequency determined by the RHN mass and active-sterile mixing. For RHN masses in the range $0.1$ - $10$ GeV and mixing $10^{-12} \lesssim |V_{eN}|^2 \lesssim 10^{-5}$, we estimate such characteristic frequencies and the signal-to-noise ratio to assess the detection prospects in GW observatories such as THEIA, $\mu$-ARES, LISA, BBO and ET. We find complementarity between GW signals and laboratory searches in SHiP, DUNE and LEGEND-1000. Notably, RHN masses of $0.2$ - $2$ GeV and mixing $10^{-10} \lesssim |V_{eN}|^2 \lesssim 10^{-7}$ are testable in both laboratory experiments and GW observations. Additionally, GW experiments can probe the canonical seesaw regime of light neutrino mass generation, a region largely inaccessible to laboratory searches.

Multi natural inflation is studied in the context of warm inflation. We study the warm multi natural inflation scenario with both linear and cubic dissipation coefficients. The model is motivated by axion-like inflation models with coupling to non-Abelian gauge fields through a dimension five coupling and dissipation originating from sphaleron decay in a thermal bath. Both cases of dissipation coefficients can be compatible with current observations. In the case of the cubic dissipation coefficient, we find that the curvature perturbation starts to grow suddenly when a transition from a weak dissipation to a strong dissipation regime occurs at the later stage of the inflation. We also show that such rapid growth of the curvature perturbation on small scales gives rise to abundant scalar induced gravitational waves, which may be detectable with future gravitational wave detectors such as DECIGO and ET. On the other hand, there are also other parameter regions of the model, in the warm inflation regime of weak to strong dissipation and with sub-Planckian axion decay constant, that can lead to overproduction of primordial black holes on small scales, which are constrained by nucleosynthesis bounds, thus ruling out the model in this region of parameters.

Lorentz-invariance violation (LV) at energy scales approaching the Planck regime serves as a critical probe for understanding quantum gravity phenomenology. Astrophysical observations of gamma-ray bursts (GRBs) present a promising avenue for testing LV-induced spectral lag phenomena; however, interpretations are complicated by degeneracies between LV effects and intrinsic emission delays. This study systematically investigates three competing time delay models: Model A (LV delay combined with a constant intrinsic delay), Model B (energy-dependent intrinsic delay without LV), and Model C (LV delay combined with energy-dependent intrinsic delay). We utilize mock GRB datasets generated under distinct delay mechanisms and employ Bayesian parameter estimation on simulated observations of 10 GRBs. Our findings demonstrate that Model C consistently recovers input parameters across all datasets. In contrast, Models A and B struggle to reconcile data generated under alternative mechanisms, particularly when confronted with high-energy TeV photons from GRB 190114C and GRB 221009A. Our analysis confirms that the incorporation of energy-dependent intrinsic delays in Model C is essential for establishing robust LV constraints, effectively resolving prior ambiguities in the interpretation of multi-GeV and TeV photon emissions. The results validate Model C as a generalized framework for future LV searches, yielding a subluminal LV scale of \(E_{\rm LV} \simeq 3 \times 10^{17}\) GeV based on realistic datasets. These findings are consistent with earlier constraints derived from Fermi-LAT datasets. This work underscores the necessity for joint modeling of LV and astrophysical emission processes in next-generation LV studies utilizing observatories such as LHAASO and CTA.

This study evaluates the possibility of efficient radio emission generation in the bow shock region of hot Jupiter-type exoplanets. As a source of energetic electrons, the shock drift acceleration mechanism at a quasi-perpendicular shock is proposed. Electrons reflected and accelerated by the shock propagate through the relatively dense stellar wind plasma and excite plasma waves; therefore, a plasma emission mechanism is considered as the source of the resulting radio waves. Using the bow shock of the hot Jupiter HD 189733b as a case study, the properties of the energetic electron beam, the excited plasma waves, and the resulting radio frequencies are estimated. An energy-based analysis is carried out to identify the range of stellar wind parameters for which radio emission from the bow shock of the exoplanet HD 189733b could be detectable by modern astronomical instruments.

We describe an exact solution representing a bouncing cosmology in the Minimal Exponential Measure (MEMe) model. Such a solution, obtained by means of the linearization around small values of the characteristic energy scale q of the theory, has the peculiarity of representing a complete bounce model that can be used to explore quantitative processes in non-singular cosmologies.

We present a reanalysis of 17 gravitational-wave events detected with Advanced LIGO and Advanced Virgo in their first three observing runs, using the new IMRPhenomTEHM model -- a phenomenological time-domain multipolar waveform model for aligned-spin black-hole binaries in elliptical orbits with two eccentric parameters: eccentricity and mean anomaly. We also analyze all events with the underlying quasi-circular model IMRPhenomTHM to study the impact of including eccentricity and compare the eccentric and quasi-circular binary hypotheses. The high computational efficiency of IMRPhenomTEHM enables us to explore the impact of two different eccentricity priors -- uniform and log-uniform -- as well as different sampler and data settings. We find evidence for eccentricity in two publicly available LVK events, GW200129 and GW200208_22, with Bayes factors favoring the eccentric hypothesis over the quasi-circular aligned-spin scenario: $\log_{10}\mathcal{B}_{\mathrm{E/QC}}\in\left[1.30^{+0.15}_{-0.15}, 5.14^{+0.15}_{-0.15}\right]$ and $\log_{10}\mathcal{B}_{\mathrm{E/QC}}\in\left[0.49^{+0.08}_{-0.08}, 1.14^{+0.08}_{-0.08}\right]$, respectively. Additionally, the two high-mass events GW190701 and GW190929 exhibit potential eccentric features. For all four events, we conduct further analyses to study the impact of different sampler settings. We also investigate waveform systematics by exploring the support for spin precession using IMRPhenomTPHM and NRSur7dq4, offering new insights into the formation channels of detected binaries. Our results highlight the importance of considering eccentric waveform models in future observing runs, alongside precessing models, as they can help mitigate potential biases in parameter estimation studies. This will be particularly relevant with the expected increase in the diversity of the binary black hole population with new detectors.

Space plasmas in various astrophysical setups can often be both very hot and dilute, making them highly susceptible to waves and fluctuations, which are generally self-generated and maintained by kinetic instabilities. In this sense, we have in-situ observational evidence from the solar wind and planetary environments, which reveal not only wave fluctuations at kinetic scales of electrons and protons, but also non-equilibrium distributions of particle velocities. This paper reports on the progress made in achieving a consistent modeling of the instabilities generated by temperature anisotropy, taking concrete example of those induced by anisotropic electrons, such as, electromagnetic electron-cyclotron (whistler) and firehose instabilities. The effects of the two main electron populations, the quasi-thermal core and the suprathermal halo indicated by the observations, are thus captured. The low-energy core is bi-Maxwellian, and the halo is described for the first time by a regularized (bi-)$\kappa$-distribution (RKD), which was recently introduced to fix inconsistencies of standard $\kappa$-distributions (SKD). In the absence of a analytical RKD dispersion kinetic formalism (involving tedious and laborious derivations), both the dispersion and (in)stability properties are directly solved numerically using the numerical Arbitrary Linear Plasma Solver (ALPS). The results have an increased degree of confidence, considering the successful testing of the ALPS on previous results with established distributions.

Gravitational waves from merging binary black holes present exciting opportunities for understanding fundamental aspects of gravity, including nonlinearities in the strong-field regime. One challenge in studying and interpreting the dynamics of binary black hole collisions is the intrinsically geometrical nature of spacetime, which in many ways is unlike that of other classical field theories. By exactly recasting Einstein's equations into a set of coupled nonlinear Maxwell equations closely resembling classical electrodynamics, we visualize the intricate dynamics of gravitational electric and magnetic fields during inspiral, merger and ring-down of a binary black hole collision.

The KM3NeT collaboration recently reported the detection of an ultra-high-energy neutrino event, dubbed KM3-230213A. This is the first observed neutrino event with energy of the order of $\mathcal{O}(100)$PeV, the origin of which remains unclear. We interprete this high energy neutrino event results from the Dirac fermion dark matter (DM) $\chi$ decay through the right-handed (RH) neutrino portal assuming the Type-I seesaw mechanism for neutrino masses and mixings. Fruthermore, dark matter $\chi$ is assumed to charged under $U(1)_X$ dark gauge symmetry, which is sponetaneously broken by the vacuum expectation value of the dark Higgs $\Phi$. In this scenario, the DM can decay into a pair of Standard Model (SM) particles for $v_\Phi \gg m_\chi$, which we assume is the case. If the DM mass is around $440$ PeV with a lifetime $5\times 10^{29}$ sec, it can account for the KM3-230213A event. However, such heavy DM cannot be produced through the thermal freeze-out mechanism due to overproduction and violation of unitarity bounds. We focus on the UV freeze-in production of DM through a dimension-5 operator, which helps in producing the DM dominantly in the early Universe. We have also found a set of allowed parameter values that can correctly account for the DM relic density and decay lifetime required to explain the KM3NeT signal. Moreover, we have generated the neutrino spectra from the two-body decay using the HDMSpectra package, which requires the dark Higgs vacuum expectations value (VEV) to be much larger than the DM mass. Finally, the large value of the dark Higgs field VEV opens up the possibility of generating GW spectra from cosmic strings. We have found a reasonable set of parameter values that can address the KM3NeT signal, yield the correct value of the DM relic density through freeze-in mechanism, and allow for possible detection of GW at future detectors.

A novel goodness-of-fit strategy is introduced for testing models for angular power spectra characterized by unknown parameters. Using this strategy, it is possible to assess the validity of such models without specifying the distribution of the estimators of the angular power spectrum being used. This holds under general conditions, ensuring the applicability of the method across diverse scenarios. Moreover, the proposed solution overcomes the need for case-by-case simulations when testing different models -- leading to notable computational advantages.