Papers for Monday, May 04 2026

Merging galaxy clusters are a promising laboratory for measuring the self-interaction cross-section (SICS) of dark matter. However, previous studies have focused on galaxy-mass offsets, which numerical simulations have shown to be intrinsically small because galaxies remain tightly coupled to the dominant dark matter potential even with significant self-interaction. Their interpretation is further complicated by unknowns of the merger phase, geometry, and initial conditions. In this paper, we overcome these obstacles by introducing the shock-to-shock distance, traced by double radio relics, as a merger chronometer that time-stamps the post-pericenter dynamical phase. Because the propagation speed of merger shocks is nearly independent of the SICS, while the halo-to-halo distance is depressed by SIDM-induced drag, the ratio of the two distances translates directly into a constraint on sigma/m. Applying this method to a gold sample of eleven cluster mergers hosting symmetric double radio relics, we determine a 68% upper limit on the SICS of sigma/m < 0.22 cm^2/g. This is the first constraint from cluster collisions that fully marginalizes over mass uncertainty, viewing angle, collision speed, merger phase, impact parameter, and gas profile slope.

JWST has opened a new era in the study of rocky exoplanets, enabling direct characterization of their surfaces with mid-infrared spectroscopy. Different types of rock have distinct spectral features that are diagnostic of the chemical composition and other physical properties like surface texture. Measurements of these features can provide valuable clues about a planet's geologic history and interior processes. Here we report a JWST 5-12 micron thermal emission spectrum for the rocky exoplanet LHS 3844 b. It is best matched by a dark, low-silica surface, such as basalt or other olivine-rich materials. The spectrum rules out fresh powder surfaces; however, space weathering can darken the powders and make them more consistent with the data. The data also disfavor trace concentrations of CO$_2$ or SO$_2$ gas (with 5-sigma and 3-sigma upper limits of 100 mbar and 10 microbar, respectively). Taken together, these results are well fit by an old, space-weathered surface with no evidence of accumulated volcanic gases.

Solar prominences are cool and dense plasma structures floating in the hot solar corona. They are ubiquitous features in the solar atmosphere, but their formation mechanism is still unclear. Here we perform comprehensive fully three-dimensional numerical simulations of prominence formation including the physics necessary to describe all atmospheric layers of the sun. With appropriate initial conditions for the magnetic field, solar prominences form self-consistently in the simulations. The formation starts by the random ejection of a dense plasma seed from the chromosphere into the corona. Subsequently, the prominence is built up by a combination of plasma injections from the chromosphere and condensation of inflowing coronal plasma. The prominence properties qualitatively match those of observed prominences. Our findings demonstrate the importance of the dynamics at and below the solar surface in the formation and evolution of solar prominences. This suggests that subsurface dynamics should also be considered in the study of prominence eruptions, which can be associated with coronal mass ejections.

Galaxy clusters are the most massive gravitationally bound structures in the Universe and key probes of cosmic evolution. The large data volume expected from upcoming surveys requires efficient automated analysis methods for tens of thousands of clusters. We present a study using Convolutional Neural Networks (CNNs) to infer cluster structural parameters from weak gravitational lensing observations. Three architectures (VGG-Net, Inception-v4, Inception-ResNet-v2) were implemented in PyTorch and trained on 75,000 synthetic reduced shear maps generated with MOKA, simulating galaxy clusters at $z = 0.25$. The networks simultaneously predict five parameters: virial mass, NFW concentrations, substructure count, and substructure mass fraction. Tests on 5000 clusters show high accuracy for primary properties. With realistic noise ($n_{\rm gal}=30$, $\sigma_{\epsilon}=0.3$), mass predictions remain robust (RMS $\sim 1.02 \times 10^{14}$ M$_\odot$/h, $\sim20$\% deviation). Concentration estimates are stable, with VGG-22 achieving the lowest RMS. Substructure properties are more challenging, with systematic underestimation across models. Comparison with traditional shear profile fitting shows improved CNN performance. VGG-22 achieves near-unbiased mass estimates and significantly better concentration recovery, reducing systematic errors. These results demonstrate that CNNs provide an effective and scalable alternative to traditional methods, particularly suited for large survey datasets.[Abridged]

We report the discovery of two metastable neutral helium (He I*) absorbers in the Milky Way, and use the upgraded CRyogenic InfraRed Echelle Spectrograph on the Very Large Telescope to determine the helium isotope ratio, $^{3}$He/$^{4}$He, along these sightlines. We have also obtained deeper observations of a third sightline to report a $\lesssim4\%$ precision measure of $^{3}$He/$^{4}$He in the Orion Nebula. These data have allowed us to place a $2\sigma$ limit on the time-variability of He I* absorption in the Orion nebula, ${\rm d}\log_{10} [N({\rm He\,I}^{*})/{\rm cm}^{-2}]/{\rm d}t\leq7.2\times10^{-4}~{\rm dex~yr}^{-1}$ ($<0.17\%~{\rm yr}^{-1}$), suggesting that these absorbers are in radiative equilibrium. We compute new galactic chemical evolution models of the Milky Way, and use our observations to infer the primordial helium isotope ratio and a scaling factor for the yields reported by nucleosynthesis calculations. Based on the data and models that we report here, we infer a best-fit value ($^{3}$He/$^{4}$He)$_{\rm P}=(1.15^{+0.24}_{-0.21})\times10^{-4}$, which agrees with Big Bang nucleosynthesis calculations that assume the Standard Model of particle physics in combination with the baryon density inferred from the cosmic microwave background temperature fluctuations. We infer the stellar yield scale relative to the solar metallicity, $y/Z_{\odot}=2.12^{+0.31}_{-0.29}$, which is somewhat higher than previously found. Finally, we note that the forthcoming extremely large telescopes are poised to determine $^{3}$He/$^{4}$He in more metal-poor environments, to secure a model-independent determination of the primordial value.

The Probe far-Infrared Mission for Astrophysics (PRIMA) is a far-infrared (24-261 micron wavelengths) probe-class space observatory currently under Phase A study, which promises orders-of-magnitude improvement in mapping speed over its predecessors. PRIMA will field exquisitely sensitive kilopixel arrays of kinetic inductance detectors (KIDs) for the Far-Infrared Enhanced Survey Spectrometer (FIRESS) instrument. PRIMA will orbit in space at the Sun-Earth L2 point, where Planck found the energetic particle flux to be about 300/min/cm2. Thus, the possible effect of a high fluence of energetic particles on the detector sensitivity must be characterized. Previous work has suggested that bombardment of KIDs by ions can reduce the quasiparticle lifetime (Barends et. al. 2009), but the conditions of the experiment were not representative of a detector which is continuously held at sub-Kelvin temperatures in the energetic particle environment of L2 orbit. To better replicate the damage which would be produced by energetic particles in this environment, we developed a fully cryogenic irradiation experiment in which a stepper motor controls a screen which can block or reveal an alpha particle emitter. This setup can be used to irradiate aluminum KID arrays fabricated for FIRESS to well-controlled dose levels. In this work, we calculate the damage dose expected for a 5-year mission in L2 orbit, and we irradiate an array to approximately 62 percent of this level. Before and after irradiation, we measure the quasiparticle lifetimes, resonant frequencies, and quality factors of the detectors.

MASCARA-5~b/TOI-1431~b is an ultra-hot Jupiter \citep[$P_{\rm orb}=2.650237\pm0.000003\,{\rm d}$, $T_{\rm eq}=2370\pm70\,{\rm K}$, $M_{\rm p}=3.12\pm0.18\,M_{\rm Jup}$, $R_{\rm p}=1.49\pm0.05\,R_{\rm Jup}$;][]{addison2021} orbiting a bright Am star ($V=8.0\,{\rm mag}$). Recent time-series observations obtained with PEPSI@LBT during the planet's post-eclipse phases have revealed Fe~{\sc i} emission lines indicative of a thermally inverted atmosphere. These observations demonstrate that MASCARA-5~b is well-suited to atmospheric characterization via emission spectroscopy, thereby motivating further follow-up observations covering additional orbital phases to constrain the planet's atmospheric chemistry, thermal structure, and dynamics. Here we present pre-eclipse time-series observations obtained with the high-resolution optical spectrograph EXPRES@LDT. Our analysis confirms the previous detection of gas-phase Fe~{\sc i} on MASCARA-5~b's dayside (with a $5.5\sigma$ significance obtained from two nights of observations) and the fact that the thermal profile is inverted with lower and upper temperatures $\sim2000\,{\rm K}$ and $\sim4500\,{\rm K}$, respectively. A search for Fe~{\sc ii} and Cr~{\sc i} did not yield any plausible detections. We also find that the pre-eclipse signal exhibits a non-negligible blueshift of $-3.2\pm1.4\,{\rm km/s}$ potentially caused by winds.

Modulated ion acoustic waves (IAWs) -- including triggered ion acoustic waves (TIAWs) and frequency-dispersed ion acoustic waves (FDIAWs) -- are increasingly recognized as efficient drivers of electron heating in the solar wind through nonlinear wave-particle interactions. Identification of these events in the Parker Solar Probe (PSP) FIELDS burst-mode archive has so far relied on expert visual inspection and does not scale to the full mission. We present EMBER (Electron heating from Modulated Burst-mode Event Recognition), an open-source pipeline that converts PSP FIELDS Digital Burst Memory (DBM) voltage bursts into log-scaled Fourier spectrograms and applies a multi-detector, background-only anomaly detection suite. The suite combines physics-motivated detectors, classical outlier detectors, and deep learning detectors. The EMBER ensemble recovers 93% of the anomalous events at 1% FAR (1 false positive per 100 held-out backgrounds). Coincident SWEAP/SPAN diagnostics show that flagged intervals exhibit core perpendicular electron temperatures above the adiabatic cooling expectation and elevated Te/Ti, reproducing the preferential-heating phenomenology established by prior manual studies without any use of electron temperatures in the detection step.

With recent advances in exoplanet observational techniques enabling the discovery of increasingly smaller planets, a crucial question emerges in the search for habitable planets: how small can a planet be and still maintain an atmosphere? We present results from the Smaller Than Earth Habitability Model (STEHM) which examines how small a planet can be and still maintain a long-term (multi-gigayear) atmosphere for planets from 1.0$R_\oplus$ down to 0.5$R_\oplus$. The model is based on a stagnant lid planet orbiting within the habitable zone of a sun-like star. Our model demonstrates that planets $\geq$0.8$R_\oplus$ can maintain their atmospheres under our Earth-like default conditions for a solar analog star, while smaller planets lose their atmospheres. Variations from the default Earth-like values cause mostly minor variations to the planet size boundary results, with some changes allowing $\geq$0.7$R_\oplus$ planets to maintain their atmosphere. Initial carbon inventory emerges as the most influential parameter for atmospheric retention, though orders of magnitude difference to Earth values are required to make a significant difference to longevity of atmospheric retention. Planets with substantial initial carbon content, large amounts of heat producing elements, cool initial mantle temperatures and low core radius fractions show the best atmospheric retention capabilities. Our results indicate that atmospheric retention on small planets depends strongly on their formation conditions and early evolution, providing important constraints for future observations of rocky exoplanets and their potential habitability.

We analyze our new 144 ks deep Chandra observation of the sloshing cold front cluster RXJ2014.8-2430. Previous observations of RXJ2014.8-2430 with XMM-Newton shows evidence of a large scale, sloshing cold front around 800 kpc away from the cluster core. Previous shallow Chandra data also shows evidence of two younger cold fronts closer to the core. Our new deeper Chandra data allow us to analyze the fine, small scale structure of these three cold fronts. Using both beta model subtraction and Gaussian Gradient Magnitude filtering, we confirm the locations of the three cold fronts, as well as discover a large concave structure southeast of the cluster core near the outermost cold front, which could be a large Kelvin-Helmholtz instability or a gas cavity from AGN activity. Analyzing the three cold fronts, we measure the widths of the cold fronts and find them to be consistent with or lower than the Coulomb mean free paths within error, signifying that diffusion is suppressed across the cold fronts. If the concave feature is the inner rim of a cavity, we find that it has a radius in the range 200-330kpc, and would have $PV$ values in the range of $5.7 \times 10^{60}$ - $2.7 \times 10^{61}$ erg. These values would make it consistent with the some of the most powerful bubbles observed.

We explore the theoretical viability of modeling a decaying dark matter sector through a unified scalar field approach. Using exact analytical solutions of the Friedmann constraints, we map the fluid phenomenology onto a scalar field potential. Our analysis reveals that physical viability, specifically the existence of a well-defined potential minimum; inevitably forces the dark energy equation of state into the phantom domain. To resolve the kinetic pathologies at late times, we propose reinterpreting the framework within a complex scenario, mapping the imaginary transition to the angular dynamics of a $U(1)$ phase. This mapping naturally yields an ultra-light mass scale of $m_\phi \sim 10^{-33} \ \text{eV}$, classifying the model as a unified dark fluid. Finally, we employ a dynamical approach to study the effects of non-minimal coupling, proving that the phantom-dominated epoch acts as a stable, late-time cosmic attractor in this kind of cosmological scenario.

The Nuclear Bulge of the Milky Way harbors stellar populations that provide crucial insights into galaxy formation processes and serve as a nearby analog for understanding bulge formation in external galaxies. However, detailed studies of this region are severely hampered by extreme and highly variable interstellar extinction, which obscures the intrinsic stellar properties and impedes accurate stellar mass determinations. Our goal is to measure the extinction law towards the Nuclear Bulge and to estimate its stellar density. We developed a method to determine the extinction law towards the Nuclear Bulge by kinematically selecting red clump stars belonging to this region. We created a high-spatial resolution reddening map, and computed stellar mass with completeness-corrected red clump star counts, scaled from empirical measurements. We find a total-to-selective extinction ratio of $\mathrm{A_K/{E_{H-K}} = 1.259 \pm 0.074}$, and an extinction ratio of $\mathrm{A_H/A_K = 1.794 \pm 0.046}$, which are consistent with previous works. The high-spatial resolution reddening map shows clear filamentary structures, and a gradient in the extinction over the giant molecular cloud G0.253+0.016 (i.e., the Brick). From the star counts, we measured a stellar mass of $\mathrm{12.2~\pm2.6\times10^8~M_{\odot}}$ for the Nuclear Bulge, in agreement with other mass estimates.

We have examined the stellar structure of 10 nearby, low stellar mass (10^6 to 6 x 10^7 Msolar) dwarf irregular galaxies by fitting ellipses as a function of surface brightness on ultra-deep V images. These are compared to far ultraviolet images as tracers of the star formation. We find that the often asymmetrical distribution of large patches of star formation activity in dwarfs, even out to low disk surface brightness levels, skews the broad-band optical isophotes in these galaxies. We also looked for evidence of the presence of a stellar halo. Possible hints of such are found in several galaxies from irregularities in the ellipses, but a stack of seven of the galaxies shows a pure exponential out to a V surface brightness of 32.3 mag/arcsec^2 where the stellar surface density is 0.0013 +/- 0.0011 Msolar/pc^2. The extended stellar component, most likely a disk structure, is probably due to internal evolutionary processes rather than external accretion. The UBVI colors of the annuli are consistent with ages of 1-6 Gyr for the far outer stellar disk.

High-resolution spectroscopic observations of helium emission lines provide a powerful probe of accretion geometry in classical T Tauri stars, revealing regions not well traced by hydrogen lines. Parallel studies in the planetary-mass regime are lacking. In this work, we investigate helium emission from the nearby (47 pc), wide-orbit (~84 au), ~13 $M_{Jup}$ accreting circumbinary companion Delorme 1 (AB)b and use resolved line profiles to constrain their origin. We analyse 33 high-S/N VLT/UVES spectra spanning near-ultraviolet to optical wavelengths at R~50,000. We detect seven He I lines at >5$\sigma$ confidence - 3890, 4027, 4473, 4923, 5017, 5877, and 6680 Ȧ - with significant epoch-to-epoch variability. The He I 5877, 4923, 4473, and 4027 Ȧ lines are asymmetric, showing a narrow component near 0 km/s and a broad component redshifted by ~15 km/s. The accretion luminosity ($1.3^{+1.6}_{-0.7}\times 10^{-5} L_{\odot}$) and mass accretion rate ($0.7^{+0.9}_{-0.4} \times 10^{-8} M_{Jup} yr^{-1}$) inferred from the median He I line luminosities are broadly consistent with, but slightly higher than, estimates from the ultraviolet excess. We conclude that protoplanet Delorme 1 (AB)b shows asymmetric He I profiles analogous to those of classical T Tauri stars, but with much smaller narrow- and broad-component widths. The triplet-singlet line ratio, a strong correlation with ultraviolet excess and the near-zero, redshifted velocities obtained for the narrow component suggest that it originates within the post-shock region, close to the planet surface. The persistent redshift of the broad component, its line width, and velocity correlation with the narrow component imply an origin within the shock structure, closer to the shock front. Emission seems to be dominated by accretion based on the obtained accretion luminosities, but a contribution from chromospheric activity may be present.

We investigate the properties of known RR Lyrae in the Vera C. Rubin Observatory Data Preview 1 (DP1) fields and compare those with the predictions based on stellar pulsation models tailored to the Legacy Survey of Space and Time (LSST) filters. The cross-match of the DP1 data with two public variable star catalogs resulted in $\sim 600$ RR Lyrae with adequate light curve sampling in five (out of seven) DP1 fields. The majority of RR Lyrae are in the 47 Tucanae and Fornax fields. We estimated photometric metallicities for these RR Lyrae using the theoretical metallicity-color relation based on $gri$-band data, and find a good agreement with literature values where the light curve sampling is sufficient for fitting template light curves accurately. The distance modulus to all RR Lyrae in DP1 fields were determined using the theoretical period-luminosity-metallicity (PLZ) relations and the $W_{gr}$ period-Wesenheit-metallicity (PWZ) relation which has the smallest metallicity term. The distances based on PWZ relations are in good agreement with the literature values with a mean offset of $0.01\pm0.36$~mag. However, the PLZ-based distance moduli are systematically large which could be due to the theoretical calibration uncertainties that include evolved horizontal branch models. The predicted period-amplitude relations based on evolved models are also inconsistent with the amplitudes based on DP1 light curves. We conclude that the metallicity and distance estimates are sensitive to the template fitting to sparsely sampled light curves in DP1 data and future data release will significantly improve these determinations for RR Lyrae stars.

Precession is a very common phenomenon for small-scale astronomical objects. However, the precession of galactic disks, occurring on a scale larger than kilo-parsec, has barely been studied in the literature. Quantifying this precession in observations remains challenging due to the lack of high-resolution dynamical data. Cosmological simulations, where gravitational interactions are self-consistently modeled, offer a unique avenue for investigating disk precession. Leveraging the IllustrisTNG simulations, we trace the evolution of spin orientation in Milky Way-like galaxies over cosmic time. We find that disk precession is ubiquitous in galaxies and significantly affects galaxy evolution. The precession is driven by the external tidal torque originating from the anisotropic matter distribution within $30\ \mathrm{kpc}$, and is violent at $\mathrm{z} > 1$ and becomes gentler but significant at $\mathrm{z} \sim 0$, when the disks are considered dynamically settled. Disk precession can induce significant cold gas warp, which is often observed in the Milky Way and nearby galaxies. We predict that the Milky Way is precessing at a rate of $\simeq3-10$ degrees per billion years at current epoch based on its observed warp. Violent precession can heat the orbits of stars, which may eventually produce prolate elliptical galaxies. The tidal torque from central galaxies can cause the precession of nearby satellite galaxies and causes their disks to point towards the centrals, which explains the observational radial alignment. We also find that the precession of accreted cold gas stream, regulated by the galaxies' torque, is crucial for the evolution of disk galaxies.

The physics governing the boundary between the most massive neutron stars (NSs) and the least massive black holes (BHs) is currently uncertain, but could potentially be constrained with new observations. While NSs have been observed with masses up to $\sim2~M_{\odot}$, there is a dearth of electromagnetic observations of compact objects in the $\sim2-5~M_{\odot}$ range, known as the lower mass gap. Recent observations of gravitational-wave (GW) signals from binary mergers detected by the LIGO-Virgo-KAGRA (LVK) collaboration indicate that this gap is likely not empty. Rapidly distinguishing whether a candidate GW event has components in this purported mass gap can indicate the likelihood of a detectable electromagnetic counterpart, and thus inform decisions for follow-up observations. In this work we train a neural network model, GWSkyNet-MassGap, that simultaneously predicts the probability that a candidate merger has a component in the lower mass gap ($P_{\mathrm{MassGap}}$) and the probability that it involves a NS ($P_{\mathrm{NS}}$). We find that the model is able to infer information about the source chirp mass to predict $P_{\mathrm{MassGap}}$ and $P_{\mathrm{NS}}$, leading to correct predictions for high-mass mergers with $\mathcal{M}_c\gtrsim15~M_{\odot}$, but less accurate predictions for lower-mass systems which require knowledge of the binary mass ratio to break the mass degeneracy. For candidate events in the first part of LVK's fourth observing run (O4a), the model has a mean prediction error of 9% for $P_{\mathrm{MassGap}}$ and 6% for $P_{\mathrm{NS}}$. The model could be further developed to rapidly predict the source chirp mass for candidate events in future observing runs.

We investigate constraints on the high-density equation of state (EOS) of neutron star matter by analyzing the probability distributions of the endpoints of mass-radius M(R) sequences within a Bayesian weighting framework. Starting from two representative hadronic baseline EOSs, SFHo and DD2, matched at higher densities to an extended linear sigma model description and constrained to approach perturbative QCD (pQCD) results, we construct families of causal hybrid EOSs spanning a broad range of stiffness at supranuclear densities. Observational constraints from the binary neutron-star merger GW170817, mass-radius measurements from the Neutron Star Interior Composition Explorer (NICER), and candidate low-mass and mass-gap compact objects are incorporated through Bayesian likelihood weighting. This approach allows us to determine probability distributions for the maximum neutron-star mass M$_{\rm TOV}$ and the corresponding radius R$_{\rm TOV}$, i.e., the endpoints of the M(R) sequences. We find that the maximum-mass distributions are largely determined by observational constraints and show only weak sensitivity to the choice of baseline EOS, favoring values around 2.2-2.3 M$_\odot$ when the most robust constraints are applied. In contrast, the corresponding radius distributions exhibit a stronger dependence on the underlying hadronic EOS, with typical preferred values near $12\pm 1$ km. Additional tidal-deformability constraints further restrict the allowed parameter space and disfavor very stiff EOS realizations when interpreted together with the possible mass-gap neutron-star candidate. Our results demonstrate that endpoint distributions of M(R) sequences provide a sensitive and complementary diagnostic for constraining the high-density behavior of the neutron-star EOS within a multimessenger Bayesian framework.

The large-scale clustering of galaxies encodes both geometric and dynamical information about the Universe. The Baryon Acoustic Oscillations (BAO) phenomenon provides a standard ruler that constrains the cosmic expansion history, while Redshift Space Distortions (RSD) probe the growth of structure through the peculiar velocity field. In this work, we present a joint analysis of BAO and growth rate parameter, $f\sigma_{8}$, in the Local Universe out to $z = 0.1$, using the $65,331$ galaxy distances of CosmicFlows-4++ database. A distinctive property of this catalogue is the availability of real space galaxy positions in addition to the redshift space coordinates. Fitting an empirical model to the measurements we obtain $r_{\rm{BAO}}^{\rm{real}} = 132\pm 8\,h^{-1}\,{\rm Mpc}$ in real space, and $r_{\rm{BAO}}^{z} = 139 \pm 7\,h^{-1}\,{\rm Mpc}$ in redshift space, at redshift $z = 0.07$. Modeling the enhancement of the correlation function within the Kaiser formalism, we derive a constraint on the growth rate parameter $f\sigma_8 = 0.344 \pm 0.105$. This analysis demonstrates how the combination of real and redshift space clustering measurements enables a simultaneous probe of important observables of the large-scale structure. Their joint detection in the same dataset, therefore, provides a self consistent view of the structure and evolution of the Local Universe. This study may be used for consistency analyses of upcoming surveys, as DESI and 4MOST, that will also provide data in both real and redshift space.

We investigate the origin of non-Markovianity in stochastic inflation and its implications for nonlinear perturbation theory. In the Schwinger--Keldysh formulation, the noise terms sourcing the infrared (IR) Langevin equations are determined by ultraviolet (UV) modes evolving on top of the stochastic IR background. Since the UV-mode evolution generally depends on the past history of the IR sector, the resulting stochastic dynamics is intrinsically non-Markovian. Working perturbatively, we derive the UV-mode solutions up to second order and decompose the corresponding noise contributions into two parts. The first is a ``deterministic'' contribution, generated by the functional Taylor expansion of the first-order UV solution around the background trajectory. The second is a genuinely ``stochastic'' contribution, originating from terms in the UV-mode equations that are quadratic in the noise variables and are usually neglected in the standard formulation of stochastic inflation. Under this conventional truncation, the deterministic contribution reduces to a Markovian correction in attractor backgrounds, whereas it could become history dependent in non-attractor phases and gives rise to non-Markovian terms involving integrals over first-order IR perturbations. We finally show that the stochastic contribution is of the same perturbative order as the deterministic one, which indicate that the conventional truncation is generically inconsistent and quadratic-noise terms may be required for a consistent treatment of nonlinear perturbations in stochastic inflation. Our analysis clarifies the perturbative structure of non-Markovianity and provides the basis for a systematic treatment of quadratic-noise effects beyond the standard formulation.

We aim to characterize the thermodynamic and dynamical conditions leading to condensation in cluster cores, and to assess the role of magnetic fields. We implement a Monte-Carlo tracer particle algorithm in the GPU-accelerated code AthenaPK, and run a purely hydrodynamical and a magnetohydrodynamical (MHD) simulations of an idealized cool-core cluster. We identify the subset of hot ICM tracers that undergo a transition to the cold phase and reconstruct their histories over a lookback time of $300\,\mathrm{Myr}$ prior to condensation. In both runs, the large majority of tracers transitioning to the cold phase follow a thermodynamic pathway driven by mixing, whereby hot ambient gas is entrained onto low-entropy seed clumps that subsequently grow into larger clouds and filaments. In the hydrodynamical run, these seeds form mainly via in-situ cooling at the edges of AGN cavities. In the MHD run, the cold gas cycle is more complex: AGN outflows occasionally shred portions of existing filaments into fragments which are then uplifted, seeding further condensation. In the MHD run, the properties of condensing tracers begin to diverge from the background ICM significantly earlier than in the hydrodynamical run (${\sim}150\,\rm Myr$ before the cooling transition versus ${\sim}30\,\rm Myr$), with vorticity and magnetic energy growing together. The turbulent Mach number at condensation is also systematically lower than in the hydrodynamical run. We examine the post-condensation evolution of individual cold structures in the MHD run, namely a massive core filament and two isolated clouds in quiescent regions. We find that magnetic tension dominates over ram pressure as the primary drag force, significantly reducing the clouds' terminal velocity. Our results demonstrate that magnetic fields substantially impact the assembly history and kinematic properties of the cold phase in cool-core clusters.

Black holes are expected to end their lifetime in a burst of Hawking radiation, emitting all Standard Model particles at ultra-high energies. The evaporation of a nearby primordial black hole (PBH) has been proposed as an explanation for the high-energy neutrino-like event reported by KM3NeT. Such a scenario requires the source to be extremely close to Earth, implying detectable gamma-ray and cosmic-ray emission. Accounting for the time-dependent field of view of gamma-ray observatories, we find that current experiments should have observed a pre-burst signal, while neutrino telescopes would also detect lower-energy events before the burst. The absence of such multimessenger signals strongly disfavors a PBH origin of the KM3-230213A event in the minimal 4D Schwarzschild scenario.

Using ''color cut'' method we obtained from SDSS DR16 catalog 733 red quasar candidates which amounted to approximately 4% of the objects from the initial sample. Then we estimated the radiative efficiency, spins, inclination angles, and corresponding new SMBH masses for all 733 objects using three theoretical models. Obtained spin distributions contain a large percentage of objects with retrograde rotation. It may indicate that these are either very young objects or objects that formed as a result of mergers. The dependencies of the estimated spin values on SMBH masses show strong correlation with linear fit slope 0.9-1.0 which allows us to assume that red quasars are likely to contain both Seyferts and NLS1, and that the main mechanism of SMBH mass growth in these objects is disk accretion.

We present the final 1.2mm and 2mm source catalogs and the redshift distribution of the mm-selected population from the NIKA2 Cosmological Legacy Survey (N2CLS) in the COSMOS field. Our aim is to provide a comprehensive dataset for studying the physical properties and evolution of high-redshift DSFGs. N2CLS covers ~1070 arcmin2 with median noise levels of 315$\mu$Jy and 91$\mu$Jy at 1.2 and 2mm, respectively. Sources are extracted with a S/N threshold of 3.9, ensuring >80% purity. Multi-wavelength counterparts are identified using high-resolution interferometric (sub-)mm data (NOEMA, ALMA) and radio observations (VLA, MeerKAT). Redshifts are compiled from spectroscopic and photometric catalogs (e.g., COSMOSWeb). The N2CLS master catalog includes 323 sources detected at >80% purity, with 104 sources detected in both bands, 197 only at 1.2mm, and 22 only at 2mm. Multi-wavelength identifications are secured for ~89% of the sample. The redshift distribution of 1.2mm sources peaks at 2.8$\pm$0.1, consistent with the epoch of peak cosmic star formation. In the total sample, we lack redshift for ~2% of the identified galaxies, plus 34 sources for which no accurate positional proxy is available, preventing the identification of a multi-wavelength counterpart. We identify 66 galaxies at z>4. The observed redshift distribution agrees well with the SIDES simulations, while four other galaxy evolution models are statistically inconsistent with the data. N2CLS is the largest contiguous deep survey to date with this depth and homogeneous coverage. This homogeneous coverage is important, as 25% of N2CLS sources lack a SCUBA2 850 mic counterpart, likely because the strongly non-uniform noise distribution of the SCUBA2 map results in lower sensitivity in parts of the this http URL released data products provide a legacy dataset for studying dust-obscured galaxy evolution.

A kinematic analysis of the young stellar association TWHya has been performed. The components of the displacement matrix in the Ogorodnikov-Milne linear model have been estimated both graphically and by solving the basic kinematic equations. The association's volume expansion with a coefficient of $K_{xyz}=103\pm9$ km s$^{-1}$ kpc$^{-1}$ was confirmed, which yields a dynamical age estimate of $t=9.7 \pm0.8$ Myr. Using the graphical method, estimates of the association's proper rigid-body rotation parameters $\omega$ around the galactic axes x and y have been obtained for the first time, with velocity values in the range of 50-70 km s$^{-1}$ kpc$^{-1}$ and errors in their determination of 14-19 km s$^{-1}$ kpc$^{-1}$. However, these values are not confirmed by another method. For example, when solving kinematic equations only using proper motions, all three components of rigid body rotation do not differ significantly from zero, $(\omega_x,\omega_y,\omega_z)=(4,7,11)\pm(5,5,5)$ km s$^{-1}$ kpc$^{-1}$.

The surface characterization of rocky exoplanets via emission spectroscopy represents a frontier of current (JWST) and future (HWO) observational efforts. Here, we implement new features in the open-source retrieval code $\texttt{POSEIDON (v1.4)}$ to fully account for an emitting and reflecting planetary surface and an overlying absorbing and scattering atmosphere. We show that realistic rocky surfaces (with wavelength-dependent albedos derived from laboratory measurements) affect emission spectra by imparting mid-infrared diagnostic absorption features, imprinting pseudo-features due to atmospheric transparency windows, and flipping absorption features to emission via surface-atmosphere interface pseudo-temperature inversions. We demonstrate that current JWST spectral data can distinguish between tenuous (low surface pressure, $\leq$ 1 bar) and thick (high surface pressures, $\geq$ 0.1 bar) atmospheres by performing atmosphere + surface retrievals on published JWST emission data of the rocky worlds TOI-1685b and 55 Cancri e. We then explore JWST MIRI LRS's capability to constrain surface geology of rocky worlds, finding that with sufficient SNR retrievals can distinguish between granite-like and basaltic surfaces for synthetic datasets. Finally, we provide an open-source database of lab-derived surface albedos (in the form of directional-hemispherical reflectances), organized by geologic classification and include supplemental tables developed to foster future collaboration between geology and exoplanet science. Our atmosphere + surface retrieval technique provides a pathway to probe geologic processes on rocky exoplanets, showing that upcoming JWST data for terrestrial worlds will enable a deeper exploration of rocky surfaces beyond our Solar System.

Coating thermal noise in high-reflectivity test-mass mirrors is a major limitation for future gravitational-wave detectors, especially in the 10--300 Hz band. ET-Pathfinder therefore requires mirror coatings that combine very high reflectance at 1.55 micrometer with low thermal noise under cryogenic conditions. Conventional dielectric Bragg mirrors provide high reflectance but require thick coatings, whereas metasurface mirrors can reduce coating-related noise but are limited by fabrication tolerances and line-edge roughness. We present a hybrid metasurface--Bragg mirror concept tailored to ET-Pathfinder. The design combines a fabrication-tolerant one-layer metasurface, an anti-resonant Fabry--Perot spacer, and a reduced dielectric Bragg stack. Optical performance is evaluated using full-wave electromagnetic simulations, while fabrication robustness is assessed with a truncated-Gaussian Monte Carlo analysis. Line-edge roughness is included as a systematic edge-smoothing effect. The resulting reflectance distributions are used to determine the minimum Bragg-stack support required to meet system-level specifications. The ideal metasurface exceeds 99.999% reflectance. When fabrication uncertainties and line-edge roughness are included, the metasurface reflectance is limited to about 99.9% at the 95% yield level. The remaining transmission can be compensated by a supporting Bragg stack with as few as seven layer pairs. For this configuration, the hybrid mirror achieves a total thermal displacement noise about one order of magnitude below the projected ET-Pathfinder coating-noise budget. These results show that fabrication-limited metasurface reflectance can be compensated within a hybrid architecture, enabling reduced coating thickness and thermal noise for next-generation gravitational-wave detectors.

We present the discovery of a second coherent radio burst from the thermally emitting neutron star 2XMM J104608.7$-$594306 in our follow-up observations with the Murriyang Ultra-Wideband Low receiver. This burst shows complex morphology with multiple components and wideband emission spanning from 704 to 4032MHz. We measured a steep spectral index of $\alpha=-2.18\pm0.16$. Our polarimetric analysis demonstrates that the burst is highly polarised with a linear and circular polarisation fraction of 54% and 22%, respectively. We identified an orthogonal jump in the polarisation position angles of the burst, resembling those seen in radio pulsars. We compared this burst with the first radio burst detected from the source with MeerKAT. These two bursts detected in a total of 40 hours on source with MeerKAT and Murriyang, combined, show that 2XMM J104608.7$-$594306 can emit sporadic radio emission with luminosity jumps comparable to those seen in the bright bursts from SGR 1935+2154. This suggests that previously thought radio-quiet neutron stars such as X-ray dim isolated neutron stars and central compact objects could exhibit rare radio bursting activity.

Quintessential inflation models provide a framework that simultaneously describes inflation and dynamical dark energy, the latter of which has recently received growing support from DESI observations. A distinctive feature of these models is the kination phase after inflation, which enhances primordial gravitational waves at high frequencies. In this work, we study a class of alpha-attractor quintessential inflation models using a fully numerical approach that follows the scalar-field evolution from inflation to the dark-energy-dominated era, allowing us to compute with high precision both the dynamics of dark energy and the primordial gravitational wave spectrum. Using the latest observational data, including DESI and ACT, we constrain the model parameters and show that the model becomes disfavored once constraints from the gravitational-wave contribution to the effective number of relativistic degrees of freedom, {\Delta} Neff, are included. This is because the model predicts a scalar spectral index ns that becomes too small to remain consistent with observations when the gravitational-wave abundance is constrained to stay below the {\Delta} Neff bound. Finally, we present the resulting primordial gravitational wave power spectrum computed using our constrained parameter values, which highlights prospects for detection by future CMB B-mode experiments at low frequencies and by gravitational-wave interferometer experiments at high frequencies.

Supermassive black holes (SMBH, $M_{\rm{BH}} > 10^8 M_\odot$) powering luminous quasars already exist one billion years after the Big Bang, yet their connection to their star-forming host galaxies, their relation to the general galaxy population and their contribution to Reionisation remains deeply enigmatic. JWST is finding numerous Active Galactic Nuclei (AGN) in high-redshift galaxies with black hole masses that appear to be over-massive compared to their host's stellar mass, but rarely as massive as those found in luminous quasars. Here we report JWST/NIRSpec observations revealing overmassive SMBH in two ultra-luminous Lyman-$\alpha$ emitters at $z\sim6.6$ that exhibit rare double-peaked Lyman-alpha profiles. The broad Balmer lines indicate black hole masses $M_{\rm{BH}}\simeq 2\times10^8 M_\odot$, matching that found in faint $z\sim 6-7$ quasars, and very high BH-to-stellar-mass ratio ($\sim 0.1-0.2$) that exceed the local relation by a factor $\sim$400-800. Stellar population modelling favours young ages ($<50$ Myr), inconsistent with the sustained average Eddington-rate accretion required to reach the observed BH masses by $z=6.6$. The double-peak Lyman-$\alpha$ profiles require a large ionised bubble and high photoionisation rate that is consistent with the ionising output of quasars powered by black holes of similar mass, thus constraining the cessation of the last quasar episode to $<1$ Myr. We interpret both systems as post-quasar galaxies in which AGN feedback has delayed stellar mass assembly, and propose that episodic quasar activity partially explains the unexpected prevalence of large ionised bubbles deep into the Epoch of Reionisation.

Transition region (TR) loops are arcade-like features in the solar transition region, with temperatures roughly between $2\times10^4$ K and $6\times10^5$ K. They are a fundamental building block of TR, which are results of the coupling between the magnetic field and the TR plasma. Their dynamics is closely related to the transport of energy and mass through the TR. Studies on this class of loops since the launch of the Interface Region Imaging Spectrograph (IRIS) have revealed that they are distinct from coronal loops. Observations have revealed that they are associated with many small-scale dynamic phenomena in the TR, from which one can infer the physics behind the energy and mass transfer in a magnetically confined environment at TR temperature. This review summarises the observational results of TR loops, showing their morphology, dynamics, plasma parameters, their relationship with flux emergence, their heating properties, and their implication in the heating of the solar atmosphere. This class of magnetic loops is much less well understood than their coronal counterparts. This review also concludes with several critical questions that need to be answered in the coming era with more advanced observational techniques and more precise and realistic simulations.

We investigate the physical origin and location of high-energy neutrino emission in active galactic nuclei (AGN) using the blazar GB6 J1542+6129 as a case study, testing whether neutrinos are produced in compact regions near the black hole or in parsec-scale jets. This question is central to understanding the conditions under which hadronic processes become efficient in AGN environments. We perform a multimessenger analysis combining ~17 years of Fermi-LAT gamma-ray data, including a 5% adaptively binned light curve and Bayesian block decomposition, with ~14 years of VLBI/MOJAVE observations to derive the Doppler factor evolution of the radio core. These are compared to the temporal properties of a suspected IceCube neutrino flare with a duration of $147^{+110}_{-25}$ days. We find that the suspected neutrino flare precedes both a gamma-ray flare and a pronounced increase in the VLBI core Doppler factor by ~1 year. This delay is consistent with the propagation time of a disturbance from the central engine to the 15GHz radio core. The duration of the post-flare gamma-ray activity is similar to that of the neutrino flare. The broadband gamma-ray spectral energy distributions remain consistent in shape across the full, flare, and post-flare intervals, indicating stable particle acceleration conditions. The temporal ordering favors neutrino production upstream of the VLBI core. GB6 J1542+6129 provides evidence for spatially separated neutrino and gamma-ray/radio emission regions in AGN. The observations are consistent with a disturbance-driven, multi-zone scenario in which neutrinos are produced in a compact, photon-rich region near the central engine, while the same disturbance later enhances Doppler-boosted leptonic emission at the parsec-scale VLBI core. These results demonstrate the power of multimessenger observations in constraining the origin of astrophysical neutrinos.

This study performs a multi-wavelength astrometric and photometric examination of a high-confidence sample $(N = 120,418)$ derived from a parent population of 2.36 million unique WDSS-seeded systems. By establishing an empirical polynomial ridge line for the broader Gaia-2MASS-Pan-STARRS subset, we calculated magnitude residuals $(\Delta G)$ to probe the systemic limits of the Gaia single-star model. Results reveal a distinct "Detection Gap" manifested as a tri-modal distribution: 14,705 stars $(12\%)$ were identified as overt Astrometric Discordance failures $(\mathrm{RUWE} > 1.4)$, while a significant subset of candidates exhibits signs of Astrometric Suppression -- where dual-flux profiles are absorbed into a stable single-star solution $(\mathrm{RUWE} < 1.4)$ despite the physical presence of a companion. Crucially, while the raw failure rate reaches $12\%$ globally, we identify an asymptotic Intrinsic Binary Residual (IBR) of $\approx 7.0\%$ $(\approx 8,429$ sources) that persists independently of stellar density. Utilizing a "Triple Constraint" framework -- astrometric noise (RUWE), photometric excess $(\Delta G)$, and the absence of official Non-Single Star (NSS) classification -- we identify a $5.9\%$ Detection Gap subset within the Gaia-2MASS audit chain that is consistent with a population of "orphaned" binaries clustered near the theoretical -0.75 magnitude "Binary Ridge." This $7\%$ floor is interpreted here as an apparent global sensitivity limit in the Gaia pipeline, suggesting that local stellar mass density models -- which rely on single-star mass-to-light ratios -- may require a quantifiable correction to accurately reflect the local baryonic mass budget.

We investigate the impact of baryonic feedback on two-point and three-point correlation functions (2PCFs and 3PCFs hereafter, respectively) involving galaxy density fields (g) and weak lensing shear fields (G), from simulated photometric catalogs of galaxies. Specifically, we baryonify high-resolution simulation using a baryonic correction model (BCM) and explore the consequences down to sub-arcminute (arcmin) scales, varying two model parameters with the largest impact on our probes: $M_{\rm c}$, which governs the amount of gas expelled beyond the halo boundary, and $\theta_{\rm ej}$, which encodes the maximal ejection radius relative to halo boundary. We create lensing maps and galaxy catalogs assuming survey properties of the upcoming Year-10 data for the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), and investigate the impact of baryonic feedback on the observed correlations, including the galaxy--galaxy--shear (ggG) and the galaxy--shear--shear (gGG) 3PCFs, which are measured, for the first time from simulations, with \textsc{TreeCorr}. Focusing on equilateral 3PCFs, we find that small scales are more heavily affected by baryonic effects than the corresponding 2PCFs, by up to 90 percent depending on the probe, redshift and BCM model. The galaxy--galaxy--galaxy (ggg) 3PCF is significantly affected at scales smaller than about 4 arcmin; a similar effect occurs at 10 arcmin for the ggG 3PCF, at 40 arcmin for the gGG 3PCF, and at about a degree for the shear--shear--shear (GGG) 3PCF. These four three-point statistics, which are collectively referred to as the $4\times3$PCFs, can be used at large scales to robustly constrain cosmological parameters. At smaller scales, their enhanced sensitivity to baryonic effects provides valuable leverage for constraining the BCM parameters and supplying informative priors. [Abridged]

We present a framework for generating lightcone simulations tailored to the analysis of Stage-IV cosmic shear data using Higher-Order Statistics (HOS). We revisit key design choices from previous simulation campaigns and re-optimize several internal parameters, benchmarking accuracy through changes in $\chi^2$ of cosmic shear statistics under survey conditions mimicking 10 years of observations from the Legacy Survey of Space and Time (LSST). We find that discretizing the lightcone uniformly in scale factor yields higher accuracy than commonly adopted schemes such as uniform spacing in redshift or comoving distance. While $N_{\rm part} = 1024^3$ simulation particles (corresponding to a mass resolution of $m_{\rm part} = 2.08\times10^{10}M_\odot$) is sufficient to model two-point statistics up to $\ell = 5000$, we observed significant instabilities on our full suite of HOS as the number of mass shells used in the lightcone construction, $N_{\rm shells}$, is varied. In contrast, simulations with $N_{\rm part} = 2048^3$ particles ($m_{\rm part} = 2.60\times10^{9}M_\odot$) robustly reproduce all statistics considered. In this higher-resolution configuration, $N_{\rm shells}$ can be reduced to $\sim50$ with only minor deviations, no larger than $0.1-0.3\sigma$ relative to our highest-resolution case ($N_{\rm shells}\sim100$). This has been explicitly verified through a comparison between our fiducial lightcone production mode based on slicing particle snapshots and an exact lightcone mode where individual particle trajectories are solved for at runtime. We further show that the particle density per pixel can be downsampled by a significant amount for $z>1.5$, saving large computational resources with no impact on the resulting statistics. These results guide the design of upcoming simulation campaigns geared towards forward-modeling and emulation-based analyses of Stage-IV data.

We present the results of AGN--host-galaxy decomposition using JWST/NIRCam, HST/ACS, and HST/WFC3 imaging of the North Ecliptic Pole Time Domain Field (NEP-TDF). The light-profiles of 36 NIRCam-selected AGN candidates are modeled for measurement of their point sources, and point source-subtracted host-galaxy emission is used in SED modeling for star formation rate (SFR) estimation. Offsets from the canonical star-forming main sequence (SFMS) show that the host galaxies form two distinct groups distinguished by their star formation: a ``bridge'' between the moderate SFRs of radio sources and low SFRs of X-ray sources, and a cleanly-separated ``branch'' above $\Delta \rm SFMS = -1$ whose SFR trends positively with AGN fraction. Branch galaxies include late-type galaxies with X-ray and radio detections and more dominant point sources that are most certainly AGN, while bridge galaxies have predominantly early-type morphologies with weaker point sources that may be due to compact stellar bulges. Both groups show evidence of recent transition between star formation and quiescence, but neither group shows preference for higher or lower stellar mass or redshift, suggesting that star formation in NIRCam-selected AGN-hosts is more strongly determined by AGN activity than by stellar mass.

In [26] it had been found that gravitational particle production (to be more specific, gravitational vacuum polarization) results in an effective increase in the directly measured value of the Hubble constant $H_0$ while it does not affect the value of the Hubble constant derived from energy densities $\bar{H}_0$. It had also been pointed out that this may explain why the Hubble constant $H_0$ determined from direct measurements (such as in SN Ia measurements) and the Hubble constant determined from indirect measurements (such as in CMB calculations in the framework of $\Lambda$CDM) are different. In the present study, first I correct a misidentification in \cite{Erdem-Universe}, namely, $\hat{H}_0=\left(\frac{\bar{H}_0}{H_0}\right)\bar{H}_0$ (rather than $\bar{H}_0$) is the value of the Hubble constant measured in CMB and BAO measurements. Then I extend the analysis to the $\sigma_8$ tension, and to determination of the Hubble constant through observations of fast radio bursts. It is observed that inclusion of the effect of gravitational vacuum polarization essentially does not neither mitigate nor exacerbate the $\sigma_8$ tension (while it mitigates or relieves the Hubble tension). This result is significant in the light of the studies in literature that question existence of a true $\sigma_8$ tension. Moreover, the present framework predicts that the value of the Hubble constant measured in fast radio bursts is $\hat{H}_0$ as in CMB and BAO measurements. This may be checked with observations in future after more precise and conclusive measurements of $\hat{H}_0$, $\bar{H}_0$, $H_0$.

After the Planck mission's launch in 2009, bolometers of its High Frequency Instrument (HFI) were considerably affected by cosmic rays, which necessitated several years of post-treatment to clean the data. To study the susceptibility of high sensitivity cryogenic detectors to particle impacts, IAS has developed the DRACuLA facility. We present the results of the latest two test campaigns performed on new generation of detectors.

We propose a simple and efficient pathway for the formation of precursors to core nucleobases in DNA and RNA using a suite of computational chemistry methods. Benzene, which is thermochemically stable in N2- or CO2-dominated atmospheres, could have formed via upper-atmospheric photochemistry or surface lightning and accumulated on the early Earth or Mars. However, nitrogen insertion into the benzene ring to form pyrimidine and purine is widely considered to be challenging. We propose that nitrogen incorporation occurred through HCN 1,4-cycloaddition to benzene's {\pi}-system, followed by a C2H2 fragmentation mechanism, as confirmed by quantum chemistry calculations. This pathway, potentially facilitated by photochemistry at the ocean surface or episodic impact events on local reservoirs, can lead to pyrimidine formation, which can further react with NH3 and HCN to produce purine. Extending this pathway to early Mars, our photochemical model simulates heterocyclic compound formation under cold, dry surface conditions that favor high benzene and HCN concentrations but lack liquid water. We thus propose that organics formed during dry phases may have later dissolved into surface waters during wet phases and become concentrated as ocean sediments. This result supports Mars Sample Return efforts focused on ancient aqueous environments likely to retain prebiotic signatures.

We introduce a diagnostic -- the quantum noise fraction $\beta$ -- that determines the maximum sensitivity improvement achievable through quantum enhancement for any gravitational wave detector. Applied to the landscape of proposed high-frequency (kHz-GHz) detectors, this diagnostic reveals that resonant mass detectors operating through tidal coupling are thermally dominated ($\beta \approx 0$) at all frequencies below ~230 MHz at dilution temperatures, rendering squeezing and entanglement limited in effectiveness. Only above this thermal frontier, defined by $\hbar \omega = k_B T \ln 3$, does the quantum regime become accessible. We identify a single concrete realization: a bulk acoustic wave resonator at 1 GHz and 10 mK ($\beta = 0.98$), and propose a gravitational wave detector employing squeezed phononic states via circuit QED readout. An array of $10^4$ such resonators with 10 dB mechanical squeezing reaches $\sqrt{S_h} = 7.6 \times 10^{-26}/\sqrt{\rm Hz}$ -- still a factor ~$10^9$ above the BBN bound on stochastic backgrounds at 1 GHz, indicating that the sensitivity gap remains predominantly classical in origin and that concurrent advances in classical detector parameters will be required.

The growth of adiabatic density perturbations during an era of early matter domination induces $\mathcal{O}(1)$ fluctuations in pre-inflationary QCD axion dark matter across a broad, string-theory-motivated parameter space. Remarkably, at $\Lambda$CDM matter-radiation equality the scale of these perturbations coincides with the quantum Jeans scale, so they collapse to solitonic ``axion stars''. These axion stars have densities up to $10^4\,\mathrm{eV}^4$, and, including their surrounding halos, they contain as much as $50\%$ of dark matter. Direct searches for a smooth axion background can be suppressed, but transient enhancements or indirect astrophysical signals at axion masses $m_a\lesssim 10^{-5}\,{\rm eV}$ would point to a non-standard cosmological history.

Accurate models of merger remnants are increasingly important for gravitational-wave science, including precision tests of gravity with ringdown, inference of black-hole populations, and modeling hierarchical mergers. For eccentric binaries, remnant mass, spin, and recoil carry nontrivial imprints of eccentricity that are both physically informative and more challenging to model, yet remain less developed than in the quasi-circular case. We present two new models trained on numerical-relativity (NR) simulations of unequal-mass, non-spinning eccentric binary black holes: NRSurE_q4NoSpin_Remnant, which predicts remnant properties, and NRSurE_q4NoSpin_Dynamics, a time-domain surrogate for the evolution of eccentricity and mean anomaly. Both models are trained on NR simulations over a three-dimensional parameter space with mass ratios $q \leq 4$, eccentricity $e < 0.23$, and mean anomaly $\ell \in [0,2\pi)$ radians, where both $e$ and $\ell$ defined at $t=-1000M$ relative to peak amplitude and $M$ is the total mass. We highlight some applications, including the phenomenological impact of eccentricity on remnant properties and the enhancement or suppression of recoil. We also provide error estimates for all modeled quantities, supporting reliable use in current and future gravitational-wave parameter-estimation analyses. Both models will be made available through open-source codes.

During the propagation of interplanetary coronal mass ejections (ICMEs), evolution of the ICME-driven shock along with interactions with other solar wind structures, planetary bodies, and general changes to their morphology can alter particle acceleration efficiency and transport effects at their associated shocks. While the underlying mechanisms for these processes have been studied, the connection between the radial evolution of the ICME-driven shock during propagation and resulting gradual Solar Energetic Particle (SEP)and Energetic Storm Particle (ESP) intensities, composition, and acceleration has yet to be fully understood. The current distributed array of spacecraft at varying heliocentric distances provides a welcome opportunity to statistically analyze the radial dependency of particle populations and acceleration mechanisms present at ICME-driven shocks. We compile a database of 39 multipoint ICME events from 2016-2023, which are observed in situ by at least two of the following spacecraft: Parker Solar Probe (PSP), Solar Orbiter, ACE, Wind, and STEREO-A. Using the magnetic field, plasma, and ion compositional data provided by these spacecraft, we derive both local shock and ESP spectral shape parameters. By comparing the changes in these parameters at different stages of ICME propagation, we analyze the connection between the evolution of the local shock conditions and the spectral shape. We find evidence to suggest a consistent increase in shock acceleration efficiency with heliocentric distance while the parent ICME is within 0.7 au, followed by a reduction in shock efficiency at further distances.

In this work, we report on the high-resolution infrared spectrum of CCH$^+$ ($^3\Pi$) recorded in the range $3066-3184$~cm$^{-1}$ by means of leak-out spectroscopy. This spectral range covers the fundamental of the CH stretching mode and a highly excited bending vibrational mode. Based on this data (385 ro-vibrational lines), accurate spectroscopic descriptions of the ground and the two vibrationally excited states of CCH$^+$ were obtained. Besides the band origins, spin-orbit coupling constants, rotational constants, centrifugal distortion constants and $\Lambda$-doubling constants for the ground and excited vibrational states have been derived. This effective Hamiltonian analysis allowed a search for pure rotational lines of CCH$^+$ in its electronic and vibrational ground state using a two-color millimeterwave - infrared scheme. We observed all rotational transitions from $J^{\prime\prime} = 2$ up to $J^{\prime\prime} = 6$ within the $\Omega = 2$ lowest energy fine structure component with resolved hyperfine splittings. This data has already guided the first detection of CCH$^+$ in space toward the Orion Bar photo-dissociation region, and has the potential to support further astronomical searches for CCH$^+$ either through radio or infrared spectroscopy, for example with the James Webb Space Telescope.

We provide a detailed formulation for calculating the angular power spectrum of the cosmic background photons arising from the dark matter decay or annihilation in a comprehensive manner. We pay particular attention to the case of dark matter decaying or annihilating into line photons. It is pointed out that taking account of the energy resolution of the detector is essential to correctly evaluate the angular power spectrum. We apply our formulation to the observational data from infrared, optical, X-ray and gamma-ray telescopes.