Papers for Friday, Dec 05 2025

Recent analyses combining DESI DR2 BAO, Planck CMB, and Pantheon+ SNe have reported mild but intriguing deviations from the LambdaCDM model. A central challenge is to determine whether these deviations reflect genuine dynamical dark energy behavior or instead arise from cross-probe inconsistencies, prior choices, or mismatches in likelihood construction. Previous work demonstrated that imposing a biased supernova-motivated prior on Omega_{m0} can artificially displace the BAO-inferred w_0,w_a values from the LambdaCDM expectation. A complementary pedagogic study further showed that the differing degeneracy geometries of BAO, CMB, and SNe can generate apparent dark energy evolution even when the underlying cosmology is exactly LambdaCDM. In this manuscript, we develop a controlled tension injection framework that provides a systematic means of quantifying how probe-level tensions influence inferred dark energy parameters. Self-consistent BAO, CMB, and SNe mock datasets are augmented with parameterized shifts in (Omega_{m0}, H_0), supernova absolute calibration, and the BAO sound-horizon scale $r_d$. The resulting datasets are analyzed through a unified MCMC pipeline, enabling a direct assessment of how these controlled tensions propagate into biases in (w_0, w_a) and the pivot equation-of-state parameter w_p. This forecasting framework provides practical guidance for identifying probe combinations that are most susceptible to spurious signatures of dynamical dark energy, and helps ensure robust multi-probe inference in forthcoming precision surveys.

One of the nearest and best studied open clusters, Pleiades is an important cornerstone of stellar astrophysics. Despite its role as reference coeval stellar population, its multiplicity properties remain vaguely determined. The combined use of Gaia DR3 multiband photometry, astrometric parameter RUWE, non-single star solutions along with available ground-based spectroscopic, high angular resolution, and polarimetric observations enable more robust constraints on the binary star population in the cluster. Several conclusions may have broader implications for other stellar populations. Twin binaries, with mass ratio close to $q\sim 1$, tend to have lower RUWE, increasing their membership selection probability, relative to $q\sim 0.5$ systems that are disfavored. The frequently observed peak in mass ratio distribution for $q\sim 1$ binaries may be partially attributed to this bias. Photometrically fitted mass ratio is underestimated for double-lined spectroscopic binaries in agreement with other authors. Differential extinction photometrically mimics stellar binarity. An area of enlarged absorption is traced by increased polarization south of the Merope star and excluded from the analysis to avoid this bias. The fraction of systems with $q>0.6$ companions is measured to be $f=16.4\%^{+2.6}_{-0.6}$ for $m>0.5~M_\odot$ stars, which is larger than recent Gaia-based estimates, but compatible with the pre-Gaia values for Pleiades and the field population. Binary fraction shows no steady increase with stellar mass in the 0.5 $-$ 1.2 $M_\odot$ range, while mass ratio has a bimodal distribution with a minimum near $q\sim 0.7$.

The Vera C. Rubin Observatory is slated to observe nearly 20 billion galaxies during its decade-long Legacy Survey of Space and Time. The rich imaging data it collects will be an invaluable resource for probing galaxy evolution across cosmic time, characterizing the host galaxies of transient phenomena, and identifying novel populations of anomalous systems. While machine learning models have shown promise for extracting galaxy features from multi-band astronomical imaging, the large dimensionality of the learned latent space presents a challenge for mechanistic interpretability studies. In this work, we present Minuet, a low-dimensional diffusion autoencoder for multi-band galaxy imaging. Minuet is trained to reconstruct 72x72-pixel $grz$ image cutouts of 6M galaxies within $z<1$ from the Dark Energy Camera Legacy Survey using only five latent dimensions. By using a diffusion model conditioned on the transformer-based autoencoder's output for image reconstruction, we achieve semantically-meaningful latent representations of galaxy images while still allowing for high-fidelity, probabilistic reconstructions. We train a series of binary classifiers on Minuet's latent features to quantify their connection to morphological labels from Galaxy Zoo, and a conditional flow to produce posterior distributions of SED-derived redshifts, stellar masses, and star-formation rates. We further show the value of Minuet for nearest neighbor searches in the learned latent space. Minuet provides strong evidence for the low intrinsic dimensionality of galaxy imaging, and introduces a class of astrophysical models that produce highly compact representations for diverse science goals.

AM Canum Venaticorum (AM CVn) stars are ultra-compact binary systems composed of a white dwarf (WD) primary accreting from a H-deficient donor. They are important as potential progenitors of Type Ia supernovae and laboratories for gravitational-wave studies, yet their evolutionary history remains unsolved. Three formation channels have been proposed: the WD channel, the He-star channel, and the cataclysmic variable (CV) channel. We aim to provide the first accurate measurements of the fundamental parameters of the accretor in ZTFJ225237.05-051917.4, including the abundances of key elements such as C, N, and Si, by analysing UV spectra obtained with the Hubble Space Telescope. These measurements provide new insight into the system's evolutionary history and establish it as a benchmark to develop our pipeline for application to a larger sample of AM CVns. We determine the binary parameters from photometric modelling and constrain the atmospheric parameters of the WD accretor, including Teff, logg, and chemical abundances, by fitting the UV spectrum with synthetic spectral models. We then infer the system's formation channel by comparing our results with theoretical evolutionary models. We measure a Teff=23300$\pm$600K and a surface gravity of logg=8.4$\pm$0.3, which implies an accretor mass of 0.86$\pm$0.16 solar masses. We find a high N/C abundance ratio by mass of >153. The accretor is significantly hotter than previous estimates based on simplified blackbody fits to the spectral energy distribution, underscoring the importance of detailed spectral modelling for determining accurate system parameters. Our results show that UV spectroscopy is well-suited to constraining the formation channels of AM CVn systems. We conclude that the He-star channel can be excluded based on the high N/C ratio, while the WD and CV channels remain consistent with the observations.

Giant impacts between planetary embryos are a natural step in the terrestrial planet formation process and are expected to create disks of warm debris in the terrestrial regions of their stars. Understanding the gas and dust debris produced in giant impacts is vital for comprehending and constraining models of planetary collisions. We reveal the distribution of millimeter grains in the giant impact debris disk of HD 172555 for the first time, using new ALMA 0.87 mm observations at $\sim$80 mas (2.3 au) resolution. We modeled the interferometric visibilities to obtain basic spatial properties of the disk, and compared it to the disk's dust and gas distributions at other wavelengths. We detect the star and dust emission from an inclined disk out to $\sim$9 au and down to 2.3 au (on-sky) from the central star, with no significant asymmetry in the dust distribution. Radiative transfer modeling of the visibilities indicates the disk surface density distribution of millimeter grains most likely peaks around $\sim$5 au, while the width inferred remains model-dependent at the S/N of the data. We highlight an outward radial offset of the small grains traced by scattered light observations compared to the millimeter grains, which could be explained by the combined effect of gas drag and radiation pressure in the presence of large enough gas densities. Furthermore, SED modeling implies a size distribution slope for the millimeter grains consistent with the expectation of collisional evolution and flatter than inferred for the micron-sized grains, implying a break in the grain size distribution and confirming an overabundance of small grains.

Every dark matter halo forms with a $\rho\propto r^{-1.5}$ density cusp at its center. For warm dark matter (WDM), these prompt cusps can be massive enough to influence the kinematics of dwarf galaxies. By implementing prompt cusps in the Galacticus galaxy formation model, we show that the measured velocity dispersions of Tucana V and Triangulum II are serious outliers for dwarf galaxies arising in WDM models. For thermal-relic dark matter, the three faintest Milky Way satellites together constrain the particle mass to be $m_\chi>5.8$ keV at 95 percent confidence or $m_\chi>9.4$ keV at 90 percent confidence. Improved velocity dispersion measurements for these systems could greatly refine this constraint, as could identification and kinematic characterization of more such galaxies.

We aim to study the atmospheric properties of the warm Neptune GJ 436 b by combining a set of five transit events observed with the CARMENES spectrograph with one transit from CRIRES$^+$ so as to provide the most constrained results possible at high resolution. We removed telluric and stellar signals from the data using SysRem and potential planetary signals were investigated using the cross-correlation technique. Following standard procedures for undetected species, we performed injection recovery tests and Bayesian retrievals to place constraints on the detectability of the main near-infrared absorbers. In addition, we simulated ELT/ANDES observations by computing end-to-end in silico datasets with EXoPLORE. No molecular signals were detected in the atmosphere of GJ 436 b, which is consistent with previous studies. Combined CARMENES-CRIRES$^+$ injection-recovery and Bayesian retrieval analyses show that the atmosphere is likely covered by high-altitude clouds ($\sim$ $1$ mbar) at low and intermediate metallicities or, alternatively, is very metal-rich ($\gtrsim$ $900\times$ solar), which would suppress spectral features without invoking clouds. Simulations of ELT/ANDES observations suggest a boost by nearly an order of magnitude to the upper limit in the photon-limited regime, reaching $0.1$ mbar at $10$-$300\times$ solar metallicities. The joint analysis of all useful transit observations from CARMENES and CRIRES$^+$ provides the most stringent constraints to date on the atmospheric properties of GJ 436 b. Complementary CCF-based and retrieval approaches consistently indicate that the atmosphere is either cloudy or highly metal enriched. Any weak near-infrared absorption lines, if present, are likely to be below current detection limits. However, according to our simulations, these features may be revealed with ELT/ANDES even in single-transit observations.

In this work, we explore the nature of $z>1$ galactic bars. Once thought to be highly transient, our results demonstrate otherwise. Our sample consists of nine massive ($>10^{10.5}\,\rm M_{\odot}$) star-forming barred-spiral galaxies at $z_{\rm spec} \sim 1.5$. Using rest-frame near-IR (F444W) JWST/NIRCam imaging, we apply ellipse fitting along with 1D and 2D morphological modeling to directly measure bar properties. We find that five galaxies host flat surface brightness profiles (bar Sérsic index $<0.4$), indicative of highly evolved, "mature" bars. By contrast, only two galaxies show exponential profiles, characteristic of young bars, and these are also shorter in absolute length than the flat bars. We therefore conclude that a large fraction of bars at this epoch have already matured, thereby indicating the presence of well-settled disks required to facilitate bar formation and sustained evolution well before $z\sim1.5$. To assess the gravitational impact of the bars, we calculate the maximum transverse-to-radial force ratio ($Q_{b}$). We find that $Q_{b}$ values are comparable to, or weaker than, those of bars in the local Universe, Seven of the nine bars show only a marginal increase in strength with maturity (from exponential to flat bars). Contrarily however, the remaining two bars are flat, but have the lowest $Q_{b}$ values in our sample. We hence propose that the mature bars at $z\sim 1.5$ may experience phases of weakening due to rapid gas inflows and/or minor mergers. In conclusion, our work sheds light on the rapidly evolving nature of high-z bars and paves the way for larger statistical studies.

We present JWST/NIRSpec IFU observations of a candidate runaway supermassive black hole at the tip of a 62 kpc-long linear feature at z=0.96. The JWST data show a sharp kinematic discontinuity at the tip, with a radial velocity change of $\approx 600$ km/s across 0.1'' (1 kpc). The velocity gradient, together with the projected post-shock flow velocity of $\approx 300$ km/s, is well described by a simple shock-compression model of a supersonic object, with a velocity of $v_{BH} = 954^{+110}_{-126}$ km/s and an inclination $i=29^{+6}_{-3}$ deg. The previously puzzling kinematics along the linear feature, with the observed radial velocity decreasing from $\approx 300$ km/s near the tip to $\approx 100$ km/s closer to the former host galaxy, are naturally explained as gradual downstream mixing of shocked gas with the circumgalactic medium through turbulent entrainment. The runaway black hole interpretation is further supported by the morphology of the gas at the tip of the wake and an analysis of the [OIII]/H$\alpha$, [NII]/H$\alpha$, [SII]/H$\alpha$, and [SIII]/[SII] line ratios. The line ratios are consistent with fast radiative shocks and rapid cooling, with best-fit shock velocities that are in good agreement with expectations from the black hole velocity and the shock geometry. Energy conservation over the lifetime of the wake suggests a SMBH mass of $M_{BH} \gtrsim 10^7$ M$_{\odot}$. These results confirm that the wake is powered by a supersonic runaway supermassive black hole, a long-predicted consequence of gravitational-wave recoil or multi-body ejection from galactic nuclei.

We analyse black-hole scaling relations at the high-mass end, focusing in particular on the regime of ultra-massive black holes, $\mathrm{M}_\mathrm{BH} > 10^{10}\,\mathrm{M}_\odot$ (UMBHs). In a sample of 16 Brightest Cluster Galaxies (BCGs) without previous black-hole mass measurements we discover 8 UMBHs based on direct dynamical detections with triaxial Schwarzschild models. This first sample of triaxial black-hole mass determinations increases the number of known UMBHs by a factor of two and dramatically increases the constraints for BH mass scaling relations at the high-mass end. We find that BCGs are outliers in the canonical BH - $\sigma$ relation, while the size of their depleted cores - the central light-deficient region - is a much better unbiased predictor of the black hole mass and should be used as a proxy at the high-mass end. BCGs smoothly join the trend already established for massive core galaxies in previous studies. This also holds for tight correlations between core size and sphere-of-influence radius and core size and core density. All these relations strongly support the black-hole binary model for the formation of the centers of the most massive galaxies.

The observational efficiency of pulsars, defined as the ratio of the observationally derived isotropic-equivalent luminosity, $4\pi d_{obs}^2 F_{obs}$, where $F_{obs}$ is the average pulsed energy flux of a pulsar and $d_{obs}$ is its estimated distance, to its energy budget, shows a wide range of values. This dispersion is believed to be a combination of beaming effects, different geometries, and case-by-case variability of the emission mechanism efficiency, but it is not clear in what proportion. In this work we focused on the gamma-ray range and analysed the four main ingredients that likely contribute to this dispersion: the geometrical term arising from the anisotropic emission (beaming), viewing and inclination angles, the uncertainty on the pulsar distance, the uncertainty on the moment of inertia, and the intrinsic efficiency of the mechanism producing the gamma-ray emission. Estimating the expected ranges of the moment of inertia and the distance errors, and considering a geometrical and spectral model that we have recently used to fit the light curves and spectra of the entire gamma-ray pulsar population, we estimate the a priori distribution of the first three ingredients in order to obtain the a posteriori distribution of the intrinsic efficiency of the mechanism. We found the latter to peak at $\sim 5-15 \%$ and to have a dispersion of around one order of magnitude. That is, we found the intrinsic efficiency of the mechanism to be the leading factor in the observed dispersion. In addition, we found little sensitivity of these results on different distributions of the estimated pulsar distance errors, and saw that the weak, alleged correlation with the spin-down power can only explain part of the observed dispersion. This methodology can be easily applied to other geometrical models of the emission, to test the sensitivity of these results on the beaming distribution.

Using a high-resolution simulation of a dwarf galaxy, we quantify the energetic importance of magnetic fields within the different phases of its interstellar medium (ISM) on parsec scales. We show that, whilst overall the magnetic field is only energetically dominant for a small fraction of the ISM, it becomes important in the thermally unstable regime (45.2% of the mass is magnetically dominated), and in the majority of the cold neutral medium (66.1% of the mass). In the molecular gas, the magnetic field dominates more of the total mass budget (39.8%) than thermal energy, turbulent kinetic energy, or gas self-gravitating potential energy. However, much of this gas will be CO-dark. This suggests that magnetic forces are non-negligible during the formation of cold dense gas, which will slow its collapse and lead to an increase in the fraction of cold atomic, and molecular gas in the ISM. Consequently, star-forming clouds may be surrounded by a larger reservoir of cold gas than would otherwise be expected.

Aerosols and clouds are expected to be ubiquitous in exoplanet and brown dwarf atmospheres, where they can have a significant impact on transmission and emission spectra. The cloud code Virga is capable of quickly modeling cloud particle sizes as a function of altitude, and has recently been updated to include functionality for aggregates (ranging from very fluffy chains to compact fractals). We analyze the effect that these aggregates have on transmission spectra for typical warm Neptune and hot Jupiter environments, as well as their effect on emission spectra for an L-type brown dwarf, over the wavelength range 0.3 - 15 um. We find significant, measurable differences in spectra when particle shape is changed (particularly the shortest wavelengths where particle morphology strongly affects the scattering slope). We provide some intuitive rules for how non-absorbing aggregates impact spectra: when particle sizes are small compared to the wavelength of light, the most elongated and chain-like particles have the highest opacities. When particles are large, the inverse is true (the most compact shapes have the highest opacities). We present an explanation for these effects in terms of the dynamics of how the particles form and move through the atmosphere, as well as in terms of fundamental optics theory. Given the significant impact that particle shape can have on spectra, we strongly encourage the community to include shape as a free parameter in future case studies, atmospheric models, and retrievals.

Recent observational studies suggest that the interplay between internal and environmental mechanisms, in particular, the combined action of stellar bars and ram pressure stripping (RPS) may influence central star formation activity in jellyfish galaxies. However, current evidence relies on small samples, leaving open whether bars play a significant role during stripping. In this study, we analyse a sample about five times larger than those used in previous works, comprising 176 galaxies identified as RPS candidates based on optical morphological indicators such as asymmetries, debris tails, and displaced star-forming regions. To assess the impact of these processes, we examine radial $u-r$ colour profiles from SDSS imaging as tracers of the specific star formation rate (sSFR). We classify galaxies by bar presence and RPS signatures, and construct comparison samples through stepwise matching in stellar mass and environment to disentangle the individual and combined effects of bars and RPS on stellar population gradients. Our results show that central rejuvenation signals emerge in RPS candidate galaxies, becoming most evident when bars and RPS act together. Barred RPS galaxies are systematically bluer at all radii than their non-RPS counterparts, while unbarred systems display only mild or no central differences, suggesting that the observable outcome of RPS depends on the stripping stage. Furthermore, barred galaxies exhibit flatter central colour profiles than unbarred ones -- a robust signature across all matched configurations. These findings highlight the key role of bars in amplifying environmental effects on the stellar populations of jellyfish galaxies, underscoring how internal structures can modulate the observable signatures of environmental processes in galaxies.

How is the production of primordial black holes (PBHs) in single-field models of inflation impacted by the presence of additional scalar fields? We consider the effect of a spectator field - a free scalar field with sub-Hubble mass, no direct coupling to the inflaton, and which makes a subdominant contribution to the total energy density - in the context of single-field models of inflation featuring a transient phase of ultra-slow roll (USR) evolution. Despite the modest title, a spectator field can have a dramatic impact: the slow-roll evolution of the spectator prevents the combined inflaton-and-spectator system from entering into USR, which naively might be expected to preclude the production of PBHs. However, we demonstrate that the growth of perturbations is maintained or enhanced by the spectator, through the rich interplay of curvature and isocurvature perturbations. We show in a model-independent way that the single-field phase of ultra-slow-roll is replaced by two turns in field space encompassing a phase of tachyonic instability for the isocurvature perturbations and a transfer of power from isocurvature to curvature modes. Furthermore, we highlight a degeneracy between the fine-tuning of the feature in the inflaton potential and the parameters of the spectator, leading to an overall resilience of model predictions to parameter variations. This makes it easier for the underlying PBH model to accommodate both high-precision CMB constraints and production of PBHs in the asteroid-mass range.

The Large and Small Magellanic Clouds (LMC and SMC) are the Milky Way's nearest interacting galaxy pair, offering a unique laboratory for studying tidal effects on galactic disks. Despite extensive survey efforts, the three-dimensional geometry of the Clouds, particularly the putative warp of the LMC, remains poorly constrained due to incompleteness in their crowded centers and the low stellar density of their peripheries, which demand wide-field coverage. Using red-clump (RC) stars as standard candles, corrected for age- and metallicity-dependent population effects with empirically calibrated color-magnitude relations and spatially resolved star-formation histories, we construct the most detailed distance map of the Magellanic System to date. Based on $\sim$2.3 million RC stars from \textit{Gaia} DR3 combined with modern reddening maps, we measure median heliocentric distances of $50.62 \pm 2.32$~kpc for the LMC (to $\sim$23$^{\circ}$) and $60.75 \pm 2.85$~kpc for the SMC (to $\sim$12$^{\circ}$). The maps reveal substructures including the LMC Northern Arm, southern hooks, the Magellanic Bridge, and SMC peripheral over-densities, with refreshed distance estimates. Fitting the LMC disk within $7^{\circ}$ yields a global inclination of $i = 25.32^{\circ} \pm 0.10^{\circ}$ and line-of-nodes position angle of $\theta = 142.34^{\circ} \pm 0.21^{\circ}$. Most strikingly, we find the LMC periphery is warped azimuthally into a U-shaped structure reaching vertical amplitudes of $\sim$7~kpc at radii of $\sim$15~kpc. We interpret this warp as a tidal response to the most recent close passage of the SMC $\sim$300 Myr ago, providing new constraints on the orbital history and dynamical evolution of the Magellanic Clouds.

Feedback from active supermassive black holes alters the distribution of matter in the Universe by injecting energy in the neighbouring hot gaseous medium, which leads to ejection of gas from the halos of galaxy groups and massive galaxies. Recent cosmological simulations such as FLAMINGO calibrate their feedback model on the baryon fractions of galaxy groups to tune the efficiency of gas ejection. However, recent observational constraints from optically selected groups and the kinetic Sunyaev-Zel'dovich effect yield lower baryon fractions than previous studies, which indicates that feedback may be more ejective than previously thought. Here we show that models involving highly ejective feedback are inconsistent with the scaling relations of local galaxy groups in the mass range $10^{13}-10^{14}M_\odot$. We study the X-ray luminosity-temperature relation in a sample of 44 galaxy groups with high-quality XMM-Newton observations. We show that highly ejective models under-predict the luminosity of galaxy groups at fixed mass at high significance ($5.7\sigma$). This conclusion is robust against selection effects and is obtained from directly measurable and minimally correlated quantities. We point out that turning observable quantities into gas fraction estimates is challenging, especially in the context of stacking large samples of heterogeneous systems. We argue that calibrating feedback models on baryon fractions is prone to systematic uncertainties and that observable scaling relations are better suited for this task.

This work shows how human physical reasoning can guide machine-driven symbolic regression toward discovering empirical laws from observations. As an example, we derive a simple equation that classifies fast radio bursts (FRBs) into two distinct Gaussian distributions, indicating the existence of two physical classes. This human-AI workflow integrates feature selection, dimensional analysis, and symbolic regression: deep learning first analyzes CHIME Catalog 1 and identifies six independent parameters that collectively provide a complete description of FRBs; guided by Buckingham-$\pi$ analysis and correlation analysis, humans then construct dimensionless groups; finally, symbolic regression performed by the machine discovers the governing equation. When applied to the newer CHIME Catalog, the equation produces consistent results, demonstrating that it captures the underlying physics. This framework is applicable to a broad range of scientific domains.

Little red dots (LRDs) newly discovered with JWST are active galactic nuclei (AGN) that may represent black hole (BH) growth at the earliest cosmic epochs. These sources show puzzling features unlike typical AGNs, including red optical continua, weak hot-dust emission, and a lack of detectable X-rays. Previously, LRDs have often been interpreted as dust-reddened AGNs, leading to severe inconsistencies with the luminosity and BH mass densities inferred for previously known AGNs over $0<z<5$. The BH-envelope (BHE) model has been proposed to explain these characteristics, in which an accreting BH is enshrouded by a dense, optically thick gaseous envelope. In this Letter, we reanalyze the SEDs of $\sim 400$ photometric LRDs in the COSMOS-Web survey using the BHE model and reassess their implications for cosmological BH evolution. We find that the optical-NIR spectra of LRDs are well reproduced by blackbody emission with an effective temperatures of $4000-6000~\K$. Within the BHE framework, the inferred bolometric luminosities decrease by $\gtrsim1-2$ orders of magnitude compared to dust-reddened AGN assumptions. As a result, the revised luminosity function, BH accretion density, and BH mass function become consistent with those of AGNs at $z<5$. The stellar masses of LRD hosts are estimated by attributing the UV excesses to star formation. Although the resulting $M_{\rm BH}/M_\star$ ratio remains higher than the local empirical value, the excess is modest. Overall, the BHE model not only resolves the spectral features of LRDs but also brings their statistical properties into agreement with the broader cosmological BH population.

The LEBS Toolkit provides a comprehensive framework for astronomers and natural scientists to address development challenges through their field. The toolkit aligns with the Sustainable Development Goals (SDGs) and offers practical actions to achieve development outcomes. It emphasises the importance of effective communication between astronomy and development, highlighting four key components: Leadership, Externalities, Break Barriers to Entry and Support Functioning Markets & Institutions (LEBS). Using these components, scientists can promote accurate information dissemination, encourage positive community impacts, overcome monopolistic barriers, and improve market efficiency and institutional settings. The LEBS Toolkit serves as a guide for the scientific community in bridging the gap between scientific knowledge and practical development actions, ultimately contributing to global development goals.

Although TRAPPIST-1's temperate planets have the highest transmission signals of any known system, flares contaminate 50-70% of transits at the 1000 ppm level, far above 100 ppm secondary atmospheres. Efforts to mitigate flare contamination and assess impacts on radiation environments are each hampered by a lack of empirical spectral analysis and physics-based modeling. We present spectrotemporal analysis and radiative-hydrodynamic modeling of 5.5 hr of NIRISS and NIRSpec observations of 6 TRAPPIST-1 flares of 2.2-8.7x10^30 erg. Flare lines and continua are characterized using grid searches of RADYN beam-heating models spanning 10$^4\times$ in electron beam parameters. Best-fit models indicate these flares result from moderate-intensity beams with emergent electron fluxes of 10^12 erg s^-1 cm^-2 and energies $\leq$37 keV, although all models over-predict the Paschen jump. These models predict XUV, FUV, and NUV counterparts to the infrared peak fluxes of 8.9-28.9x10^27, 4.3-13.9x10^26, and 3.4-11.4x10^27 erg s^-1, respectively. Scaling the flare rate into the XUV suggests flaring contributes 1.35$_{-0.15}^{+2.0}\times$ quiescence yr$^{-1}$. We bin integrations of similar flare effective temperature to construct fiducial flare spectra from 2000-4500 K in order to develop separate empirical and RADYN-based mitigation pipelines. Both pipelines are applied to all 5.5 hr of R=10 data, resulting in maximum residuals from 1-2.8$\mu$m of 100-140 ppm and typical residuals of 54$\pm$14 and 65$\pm$17 ppm for the empirical and RADYN-based pipelines, respectively. Injection testing supports 3$\sigma$ detection capability for CO2 atmospheres with features of 150-250 ppm, with weak evidence (BF$\approx$3) still obtained at 130 ppm. Our results motivate multi-wavelength observations to improve model fidelity and test high-energy predictions.

Models of cosmic inflation generically predict a weak but potentially detectable amount of primordial non-Gaussianity (PNG), which can be used to obtain insights into the degrees of freedom during inflation and their interactions. The simplest types of PNG are the local and non-local (equilateral and orthogonal) shapes of the primordial three-point correlators, which are predicted by models with multiple light fields and derivative interactions in single-field inflation, respectively. In this paper we place constraints on local, equilateral, and orthogonal non-Gaussianities using the power spectrum and bispectrum extracted from first public release of the Dark Energy Spectroscopic Instrument (DESI). Our analysis makes use of higher-order clustering information through a consistent effective field theory (EFT) model for both the power spectrum and bispectrum at one-loop order. Using robust scale cuts where the EFT description is valid, we find the following constraints on PNG amplitudes: $f^{\rm loc}_{\rm NL}=-0.1\pm 7.4$, $f^{\rm equil}_{\rm NL}=719\pm 390$, $f^{\rm orth}_{\rm NL}=-200\pm 100$ (at $68\%$ CL). Non-local PNG constraints can be further improved by combining high-redshift DESI with legacy BOSS data and using simulation-based priors on bias parameters, yielding the strongest large-scale structure constraints to date $f^{\rm equil}_{\rm NL}=200\pm 230$, $f^{\rm orth}_{\rm NL}=-24\pm 86$. Our constraint on $f^{\rm loc}_{\rm NL}$ is competitive with the cosmic microwave background (CMB) limit; the combination gives $f^{\rm loc}_{\rm NL}=-0.0\pm 4.1$, $18\%$ stronger than the CMB only result, which represents the strongest bound on multi-field inflation yet obtained.

We present 4,110 strong gravitational lens candidates, 3,887 of which are new discoveries, selected from a sample of 5,837,154 luminous red galaxies (LRGs) observed with the Dark Energy Spectroscopic Instrument (DESI). Candidates are identified via the presence of background ionized oxygen [O II] nebular emission lines in the foreground LRG spectra which may originate from the lensing of higher redshift star-forming galaxies. Using the measured foreground redshift, background redshift, and integrated flux of the background [O II] doublet, we integrate over impact parameters to compute the probability that each candidate is a lens. We expect 53% of candidates to be true lenses with Einstein radii ranging from 0.1'' to 4'', which can be confirmed with high-resolution imaging. Confirmed strong lenses from this sample will form a valuable cosmological dataset, as strong gravitational lensing is the only method to directly measure dark matter halo substructure at cosmological distances. We independently recover the host of the multiply imaged gravitationally lensed type Ia supernova iPTF16geu (Goobar et al. 2017). Monitoring these lenses for future multiply lensed transients will enable a) H0 measurements via time-delay cosmography and b) substructure measurements via flux ratios.

Precise astrometric and photometric measurements of celestial point sources are fundamental to modern astronomy. These measurements, used to determine object positions, motions, and fluxes, are based on observational models that have evolved from empirical centroiding rules to rigorous probabilistic formulations at the pixel level. This review summarizes key contributions that formalized this transition and analyzes seminal works addressing both the theoretical limits and the empirical performance of estimators. Central to these developments is the derivation of fundamental bounds, such as the Cramér-Rao Lower Bound (CRLB), and the assessment of widely used estimators, including Maximum Likelihood (ML), Least Squares (LS), and Weighted Least Squares (WLS). These studies show that, while the CRLB sets a theoretical benchmark, practical estimators achieve it only under specific signal-to-noise ratio (SNR) regimes, with notable discrepancies in high-SNR conditions. Moreover, recent results demonstrate that jointly estimating source flux and background significantly improves photometric precision compared to sequential approaches. Looking ahead, the increasing complexity of astronomical surveys, driven by massive data volumes, dynamic observational conditions, and the integration of machine learning, poses new challenges to reliable inference. In this context, tools from statistical theory, including performance bounds and theoretically grounded estimators, remain critical to guide algorithm design and ensure robust astrometric and photometric pipelines.

We present the results of constraining the angular diameter distance to blazar TXS 0506+056. We used data obtained with the 15 GHz VLBA in MJD 54838-60262 and data from the 15 GHz OVRO 40 m single dish telescope in MJD 54474-59023. We used a variability timescale and a causality argument of a linear size to measure the angular diameter distance to the source. To constrain the Doppler factor, we applied the relation between the rest-frame brightness temperature of the emission region and the observed brightness temperature. To calculate the observed brightness temperature, the angular size and flux density variation of the emission region are required. The angular size of the emission region (i.e., the VLBA core) was obtained from a FWHM, which is a circular Gaussian model-fitting parameter that ranges from 0.048-0.228 mas, and its uncertainty is determined to be 1.8-13 %. Using the OVRO SD light curve, we obtained a variability timescale of $128.0_{-0.3}^{+0.2}$ days and a peak flux density of $1.750_{-0.104}^{+0.015}$ Jy for the largest flare that peaked on MJD $58921.7_{-5.5}^{+2.6}$. We assumed a disk brightness geometry, equipartition brightness temperature ($5\times10^{10}$ K), and perfect radius. Using the VLBA core sizes obtained near the flare peaks, we found consistent distance measurement results with the $\Lambda$CDM model within 1$\sigma$ uncertainties. We suggest that the best distance from the source is $941_{-64}^{+59}$ Mpc, which is comparable with the $\Lambda$CDM distance of $948.2\pm13.5$ Mpc. The distance measurement should indeed be taken at the peak of a flare. We found that the decomposed timescale allowed us to obtain consistent distances with the $\Lambda$CDM. We strongly suggest to decompose light curves when the variability timescales are to be obtained properly.

A Coronal Mass Ejection (CME) was detected crossing the radio signals transmitted by the Mars Express (MEX) and Tianwen-1 (TIW) spacecraft at a solar elongation of $4.4^{o}$. The impact of the CME was clearly identifiable in the spacecraft signal SNR, Doppler noise and phase residuals observed at the University of Tasmania's Very Long Baseline Interferometry (VLBI) antenna in Ceduna, South Australia. The residual phases observed from the spacecraft were highly correlated with each other during the transit of the CME across the radio ray-path despite the spacecraft signals having substantially different Doppler trends. We analyse the auto- and cross-correlations between the spacecraft phase residuals, finding time-lags ranging between 3.18-14.43 seconds depending on whether the imprinted fluctuations were stronger on the uplink or the downlink radio ray-paths. We also examine the temporal evolution of the phase fluctuations to probe the finer structure of the CME and demonstrate that there was a clear difference in the turbulence regime of the CME leading edge and the background solar wind conditions several hours prior to the CME radio occultation. Finally, autocorrelation of the MEX two-way radio Doppler noise data from Ceduna and closed-loop Doppler data from ESA's New Norcia ground station antenna were used to constrain the location of the CME impact along the radio ray-path \add{to a region 0.2 AU from the Sun, at a heliospheric longitude consistent with CME origin at the Sun. The results presented demonstrate the potential of the multi-spacecraft-in-beam technique for studying CME structures in great detail, and providing measurements that complement the capabilities of future solar monitoring instruments.

We derive a grey linear diffusion equation for photons with respect to inertial (or lab-frame) space and time, using asymptotic analysis in 1D planar geometry. The solution of the equation is the comoving radiation energy density. Our analysis does not make use of assumptions about the magnitude of velocity; instead we derive an asymptotic scaling in the lab frame such that we avoid apparent non-physical pathologies that are encountered with the standard static-matter scaling. We permit the photon direction to be continuous (as opposed to constraining the analysis to discrete ordinates). The result is a drift-diffusion equation in the lab frame for comoving radiation energy density, with an adiabatic term that matches the standard semi-relativistic diffusion equation. Following a recent study for discrete directions, this equation reduces to a pure advection equation as the velocity approaches the speed of light. We perform preliminary numerical experiments comparing solutions to relativistic lab-frame Monte Carlo transport and to the well-known semi-relativistic diffusion equation.

We investigated the occurrence rate of the sustained gamma ray emission (SGRE) events from the Sun using Fermi Large Area Telescope (LAT) data covering solar cycle (SC) 24 and the rising and maximum phases of SC 25. Due to a solar array drive assembly malfunction starting in March 2018, only small number (15) of SGRE events were observed during SC 25. Over the first 61 months, the average sunspot number (SSN) increased from 56.9 in SC 24 to 79.0 in SC 25. Fast and wide (FW) CMEs and decameter-hectometric (DH) type II bursts also increased significantly in SC 25 by 29% and 33%, respectively when normalized to SSN. Therefore, the increase in solar activity should result in a higher number of SGREs in SC 25. We estimated the number of SGREs in SC 25 using three methods. (i) If the SGRE number varies commensurate FW CMEs and DH type II bursts, SC 25 should have 45-47 SGRE events. (ii) In SC 24, ~17% of FW CMEs and 25% of DH type II bursts were associated with SGRE events. At this rate, SC 25 should have 46 SGRE events. (iii) Since SGRE events are invariably associated with >100 keV hard X-ray (HXR) bursts, we identified DH type II bursts associated with >100 keV HXR bursts from Fermi gamma ray burst monitor (GBM) during LAT data gaps. Based on our finding that SGRE events in SCs 24 and 25 were all associated with HXR bursts of duration > ~5 min, we found only 27 of the 79 LAT-gap type II bursts had >100 keV HXR bursts with duration > ~5 min. These DH type II bursts are likely to indicate SGRE, bringing the total number of SGRE events to 42 (15 observed + 27 inferred). Thus, the three methods provide similar estimates of the number of SGRE events in SC 25. We, therefore, conclude that that SC 25 is stronger than SC 24 based on the estimated number SGRE events. Halo CMEs, ground level enhancement (GLE) events, and intense geomagnetic storms are also consistent with a stronger SC 25.

We investigate the clustering and halo occupation distribution (HOD) of DESI Y1 emission-line (ELGs) and luminous red (LRGs) galaxies at $0.8<z<1.1$, including their cross-correlation (ELG$\times$LRG), using the AbacusSummit suite and a new Halo Occupation Model (HOMe) for galaxy multi-tracers. This integrates intra-halo dynamics, halo exclusion, and quenching, bridging insights from hydrodynamical, HOD, abundance-matching, and semi-analytic studies. Leveraging full phase-space information from the Uchuu N-body simulation, and sampling satellites from dark-matter particle positions via physically motivated prescriptions, HOMe reproduces the anisotropic clustering down to $s=200\,h^{-1}$kpc with unprecedented accuracy. Model parameters are inferred solely from two-point statistics using a two-level Bayesian framework, yielding high-fidelity ELG, LRG and cross-reference catalogs. We find that satellite ELGs behave as incoherent flows within their parent halos, dominating the clustering below $4\,h^{-1}$Mpc. The HOD from the best-fit HOMe has the following properties: (i) 90.50% (85.91%) of ELGs (LRGs) are central galaxies without satellites, residing in halos of $M_{\rm vir}\sim6.6\times10^{11}\,(1.2\times10^{13})\,h^{-1}{\rm M}_\odot$; (ii) the ELG$\times$LRG cross-correlation is governed by central-central pairs and shaped by halo exclusion on $2-5\,h^{-1}$Mpc scales; (iii) 9.50% (14.09%) of ELGs (LRGs) are satellites, of which 1.09% (3.52%) inhabit halos with a central galaxy of the same species in a maximally conformal configuration, 7.02% (0.005%) orbit complementary hosts in a minimally conformal state, and 0.58% (10.57%) are orphans. HOMe high sensitivity precisely captures the dynamics of satellites in different host environments, opening a promising avenue for understanding systematics, the dynamical nature of dark matter, potentially distinguishing gravity models.

In order to realize multi-spectral imaging of space coronagraph, a compact filter wheel mechanism is designed. The filters with different spectral transmittance can be cut into the optical path at different times by this mechanism. Because the image contrast of space coronagraph is very high, the high stability requirement for the optical unit of coronagraph is put forward. Small modulus worm gear and worm are taken by the mechanism, in order to realize compact structure, high stiffness, the unidirectional 360 degrees rotation and reverse self-locking function. The precision, stiffness, mechanical properties and reliability of the mechanism are analyzed in the paper. The results show that the position accuracy of the filter wheel can meet the requirement of 0.5mm. The first order modal of the mechanism is 313Hz. The results of vibration test indicate that stiffness, dynamic performance and reliability of the mechanism can be meet. Therefore, the design of filter wheel in this paper can ensure the multi-spectral imaging requirements under complex spatial conditions.

Unstable mass transfer may occur during white dwarf-neutron star (WD-NS) mergers, in which the WD can be tidally disrupted and form an accretion disk around the NS. Such an accretion disk can produce unbound wind ejecta, with synthesized $^{56}\mathrm{Ni}$ mixed in. Numerical simulations reveal that this unbound ejecta should be strongly polar-dominated, which may cause the following radioactive-powered thermal transient to be viewing-angle-dependent. This issue has so far received limited investigation. We investigate how the intrinsically non-spherical geometry of WD-NS wind ejecta affects the viewing-angle dependence of the thermal transients. Using a two-dimensional axisymmetric ejecta configuration and incorporating heating from the radioactive decay of $^{56}\mathrm{Ni}$, we employ a semi-analytical discretization scheme to simulate the observed viewing-angle-dependent photospheric evolution, as well as the resulting spectra and lightcurves. The observed photosphere evolves over time and depends strongly on the viewing angle: off-axis observers can see deeper, hotter inner layers of the ejecta and larger projected photospheric areas compared to on-axis observers. For a fiducial WD-NS merger producing 0.3 solar mass of ejecta and 0.01 solar mass of synthesized $^{56}\mathrm{Ni}$, the resulting peak optical absolute magnitudes of the transient span from ~ -12 mag along the polar direction to ~ -16 mag along the equatorial direction, corresponding to luminosities of $10^{40}$-$10^{42}$ erg s$^{-1}$. The typical peak timescales are expected to be 3-10 d. We for the first time explore the viewing-angle effect on WD-NS merger transients. Since their ejecta composition and energy sources resemble those of supernovae, yet WD-NS merger transients are dimmer and evolve more rapidly, we propose using "mini-supernovae" to describe the thermal emission following WD-NS mergers.

We present high-precision chemical abundances for 20 FGK stars hosting planets observed in JWST Cycle 2 GO programs. Using high-resolution, high-signal-to-noise ratio spectra from the ESO and Keck archives, we perform a strict line-by-line differential analysis relative to the Sun to derive stellar parameters and abundances of 19 elements from C to Zn. The stars span effective temperatures of 4500-6500 K and metallicities from -0.57 to +0.50 dex. The sample includes hosts of both gas giants and terrestrial planets, allowing direct comparison between stellar composition and planetary properties. Several of the giant planets orbit metal-rich stars. The detailed abundance patterns show clear chemical diversity, including carbon-enhanced but mildly metal-poor stars (TOI-824, TOI-1130, GJ 9827) and $\alpha$-enhanced metal-poor stars (TOI-561, GJ 9827, TOI-824). These variations trace differences in protoplanetary disk composition and may influence planetary interiors and atmospheric chemistry. The planet-hosts show a range of [C/O] ratios, and the diverse [Mg/Si] ratios may suggest varied interior compositions for their rocky planets. This homogeneous stellar abundance, together with future uniform JWST planetary atmosphere measurements, provides a foundation for exploring the planet mass-metallicity relation and the connection between stellar chemistry and planetary formation pathways. These results constitute the first step in a larger survey spanning multiple JWST cycles to systematically examine how host star composition shapes exoplanetary systems.

We utilize multi-epoch line-of-sight velocity measurements from the Milky Way Survey of the Dark Energy Spectroscopic Instrument to estimate the binary fraction for member stars in the dwarf spheroidal galaxy Ursa Minor. Our dataset comprises 670 distinct member stars, with a total of more than 2,000 observations collected over approximately one year. We constrain the binary fraction for UMi to be $0.61^{+0.16}_{-0.20}$ and $0.69^{+0.19}_{-0.17}$, with the binary orbital parameter distributions based on solar neighborhood observation from Duquennoy \& Mayor (1991) and Moe \& Di Stefano (2017), respectively. Furthermore, by dividing our data into two subsamples at the median metallicity, we identify that the binary fraction for the metal-rich ([Fe/H]>-2.14) population is slightly higher than that of the metal-poor ([Fe/H]<-2.14) population. Based on the Moe \& Di Stefano model, the best-constrained binary fractions for metal-rich and metal-poor populations in UMi are $0.86^{+0.14}_{-0.24}$ and $0.48^{+0.26}_{-0.19}$, respectively. After a thorough examination, we find that this offset cannot be attributed to sample selection effects. We also divide our data into two subsamples according to their projected radius to the center of UMi, and find that the more centrally concentrated population in a denser environment has a lower binary fraction of $0.33^{+0.30}_{-0.20}$, compared with $1.00^{+0.00}_{-0.32}$ for the subsample in more outskirts.

The elements in the universe are synthesized primarily in stars and supernovae, where nuclear fusion favors the production of even-Z elements. In contrast, odd-Z elements are less abundant and their yields are highly dependent on detailed stellar physics, making theoretical predictions of their cosmic abundance uncertain. In particular, the origin of odd-Z elements such as phosphorus (P), chlorine (Cl), and potassium (K), which are important for planet formation and life, is poorly understood. While the abundances of these elements in Milky Way stars are close to solar values, supernova explosion models systematically underestimate their production by up to an order of magnitude, indicating that key mechanisms for odd-Z nucleosynthesis are currently missing from theoretical models. Here we report the observation of P, Cl, and K in the supernova remnant Cassiopeia A using high-resolution X-ray spectroscopy with XRISM, with the detection of K at above the 6$\sigma$ level being the most significant finding. Supernova explosion models of normal massive stars cannot explain the element abundance pattern, especially the high abundances of Cl and K, while models that include stellar rotation, binary interactions or shell mergers agree closely with the observations. Our observations suggest that such stellar activity plays a significant role in supplying these elements to the universe.

We report detailed analyses of deep, high-resolution spectra of three Galactic planetary nebulae (PNe) with high abundance discrepancy factors (ADFs), Hf2-2, M1-42 and NGC6153, obtained with the Ultraviolet and Visual Echelle Spectrograph (UVES) on the 8.2m Very Large Telescope (VLT). These spectra were carefully reduced, including rigorous absolute flux calibration, yielding detections of ~410-800 emission lines in each PN. Plasma diagnostics and abundance calculations were critically performed using nebular lines. In all three PNe, the electron temperatures derived using the collisionally excited lines (CELs) are higher than those yielded by the HI Balmer and Paschen jumps, while the temperatures yielded by the OII and NII optical recombination lines (ORLs) are very low, <2000 K, indicating that the heavy-element ORLs probe cold nebular regions. The ORL abundances of N, O and Ne are systematically higher than the corresponding CEL values, confirming high ADFs in the three objects. Position-velocity (PV) diagrams were created, and spatio-kinematical studies show that CELs come from the outer nebular regions, while the ORL-emitting regions are close to nebular center. Additionally, the velocity indicated by CEL line-splitting decreases with ionization potential, which was not obvious in ORLs. These spatial and kinematic differences support two distinct components of ionized gas: a cold, metal-rich component and a warmer component with normal metallicity. Heavy elements are strongly enriched in the cold gas, while its H^+ fraction is low but still produces significant HI emission, affecting CEL abundance estimates.

We find the constraint on the population of ultra-light regular primordial black holes (RPBHs) by using isocurvature gravitational waves (GW). If ultra-light RPBHs dominated the early Universe, the initial isocurvature perturbation is converted into curvature perturbation that induce second-order GW background upon evaporation of RPBHs. The upper limit of extra relativistic degrees of freedom $\Delta N_{\rm eff}$, which could be inferred from Big Bang Nucleosynthesis or Cosmic Microwave Background observations, places a constraint on the maximum energy density of GW, which in turn can be used to constrain the RPBH population. As RPBHs have different lifetime from their singular counterparts, the constraint must be modified accordingly. While the formalism that we provide is generic, we work out explicitly the case of Simpson-Visser metric for demonstration. The RPBHs associated with this metric have lower Hawking temperature and smaller horizon size, leading to a longer lifetime than the singular Schwarzschild black holes. This implies a stronger constraint on RPBH population as they dominate the Universe for a longer period of time and generate stronger GW.

We study the possibility of registering high-energy extensive air showers (EAS) by the onshore detector facility simultaneously with the trigger of the Baikal-GVD neutrino telescope. The location of the surface detector array on the shore of Lake Baikal is motivated by the fact that permanent placement of detectors on the surface of the lake is challenging. Within the given geometry, simultaneous registration is possible for EAS with a zenith angle of about 76 degrees within the limited solid angle. The corresponding inclined EAS are dominated by muons and significantly attenuated. The installation will make it possible to obtain an estimate of the number of high-energy muons in EAS. This, subsequently, would make it possible to verify EAS modeling and calculations of the atmospheric neutrino flux. The detector may also be used for cross-calibration of energy and direction measurements by the neutrino telescope. We use the CORSIKA program to simulate the registration of EAS on the shore close to the Baikal-GVD. The propagation of ultra-high energy muons produced by EAS through 3.5 km of water and their registration by Baikal-GVD is simulated using the PROPOSAL package. We calculate the minimal total area of the onshore detectors suitable for several cosmic ray energy thresholds, starting from 1 PeV. The report presents estimates of the number of events jointly registered by the onshore installation and the Baikal-GVD the EAS registration detector systems with different areas.

Observational studies have established a remarkably tight power-law relationship between the H I masses and sizes of late-type galaxies, known as the H I mass-size relation. This relation has been shown to persist across various models of a galaxy's H I surface density profile. Using the Simba cosmological simulations, we investigate the robustness of this relation under different feedback prescriptions, including cases where specific feedback mechanisms are absent. While the global properties of galaxies are significantly affected by changes in feedback, the H I mass-size relation remains intact. Moreover, its parameters consistently align with the best available empirical measurements. We analyze the H I mass distributions of galaxies and demonstrate that, regardless of the feedback scenario, galaxies within a given H I mass bin exhibit outer H I radial profiles well approximated by an exponential function. Furthermore, the exponential decline rate remains remarkably similar across different physical prescriptions. We attribute the persistence of the H I mass-size relation to this inherent self-similarity in the H I mass distributions.

One of the imminent science goals of pulsar timing arrays (PTAs) is the detection of a continuous gravitational wave (CGW) emitted by an individual supermassive black hole binary (SMBHB). SMBHBs that cause CGWs with GW frequencies $f_\mathrm{GW} > 10 \text{nHz}$ have undergone significant orbital evolution, hence a change of $f_\mathrm{GW}$ over time. In PTA data sets with sufficiently long observational time span, this means that the Earth and pulsar terms' contributions to the CGW signal signature can eventually become resolvable. Since the pulsar term is accumulated incoherently and thus often treated as an additional source of noise, this separation can prove to be beneficial for the detection of the CGW signal in the PTA data set. We aim to investigate to what extent resolvable Earth and pulsar terms affect currently used techniques for CGW searches with PTA data sets, that treat the pulsar term as an additional source noise. We focus on the dependency of the pulsar term frequencies on the pulsar's distance. We aim to answer the question of whether adding more distant pulsars to a PTA data set can mitigate biases and improve the detection of CGWs. We show that under ideal conditions, more distant pulsars can facilitate the CGW search with PTA data sets. Bayesian parameter estimation is yielding better parameter constraints and the frequentist per-frequentist optimal statistic search becomes more stable. However, using the realistic data set simulations, it was found that other configuration parameters of a PTA, such as the anisotropic distribution of pulsars and the effective number of pulsars in a PTA, can play a crucial role to the importance of this effect.

Stellar streams offer one of the most sensitive probes of the Milky Way`s gravitational potential, as their phase-space morphology encodes both the tidal field of the host galaxy and the internal structure of their progenitors. In this work, we introduce a framework that leverages Flow Matching and Simulation-Based Inference (SBI) to jointly infer the parameters of the GD-1 progenitor and the global properties of the Milky Way potential. Our aim is to move beyond traditional techniques (e.g. orbit-fitting and action-angle methods) by constructing a fully Bayesian, likelihood-free posterior over both host-galaxy parameters and progenitor properties, thereby capturing the intrinsic coupling between tidal stripping dynamics and the underlying potential. To achieve this, we generate a large suite of mock GD-1-like streams using our differentiable N-body code \textsc{\texttt{Odisseo}}, sampling self-consistent initial conditions from a Plummer sphere and evolving them in a flexible Milky Way potential model. We then apply conditional Flow Matching to learn the vector field that transports a base Gaussian distribution into the posterior, enabling efficient, amortized inference directly from stream phase-space data. We demonstrate that our method successfully recovers the true parameters of a fiducial GD-1 simulation, producing well-calibrated posteriors and accurately reproducing parameter degeneracies arising from progenitor-host interactions. Flow Matching provides a powerful, flexible framework for Galactic Archaeology. Our approach enables joint inference on progenitor and Galactic parameters, capturing complex dependencies that are difficult to model with classical likelihood-based methods.

The transition regime between giant planets (GPs) and brown dwarfs (BDs) is still an open subject of study in exoplanetary science. A complete understanding of the population of long-period GPs and BDs would be pivotal in understanding this topic, but the number of such objects with precisely measured orbital and physical parameters is still small. Moreover, their dynamical influence on smaller companions in inner orbits is still unclear. Within the GAPS programme, we aim to confirm and characterize sub-stellar companion candidates from Gaia DR3, and to study the potential presence of additional lower-mass planets in their systems. We present the results of an intensive high-precision RV monitoring of HD 128717, which hosts the astrometric candidate Gaia-ASOI-009. We used the HARPS-N spectrograph at TNG to collect a high-cadence RV time series of the target. We used MCMC analyses to refine the Gaia DR3 orbital solution of the companion and, finally, performed a combined model of RV and proper motion anomaly (PMa) to derive the complete 3-D orbit of the companion. We also ran a suite of numerical simulations to confirm our results. We confirm the sub-stellar nature of Gaia-ASOI-009, i.e. Gaia-6 B: from the combined RV+PMa fit, we confirm that it is a high-eccentricity low-mass brown dwarf with $P_\text{B} = 9.37^{+0.06}_{-0.05}$ yr, $M_\text{B} = 19.8 \pm 0.5$ $M_\text{J}$, $e_\text{B}=0.85$, $i_\text{B} = 130^{\circ}$. The derived orbital solution differs significantly from the one published in Gaia DR3. Through a series of dedicated simulations, we demonstrate that this discrepancy arises from a degeneracy in the Gaia DR3 astrometric solution. Specifically, the combination of Gaia-6 B long orbital period and high eccentricity, both poorly constrained by the limited timespan of DR3, led to an incorrect solution characterized by a shorter period and lower eccentricity.

Recent studies have shown that the irregular O-C variations observed in many non-Blazhko RR Lyrae stars may result from random, cycle-to-cycle (C2C) variations in their light curves. However, centuries-long data series reveal that the O-C diagrams of Blazhko stars exhibit particularly large-amplitude, irregular variations. In this Letter, we extend the previous investigation of non-Blazhko stars to Kepler Blazhko stars to explore the role of C2C variations in the O-C diagrams. We derived the O-C diagrams from Kepler space telescope light curves using a precise template-fitting method. Based on their Fourier analyses, we also constructed residual O-C diagrams that were pre-whitened for frequencies associated with the Blazhko effect. We then fitted the same statistical models to both types of O-Cs that we had previously applied to non-Blazhko stars. The optimal statistical model includes the C2C variation for 74% of the O-C curves in our Blazhko sample, and the parameter describing the strength of the C2C variation is significantly larger than that obtained for non-Blazhko stars. This may explain the strong irregular O-C variations previously observed in Blazhko stars. Furthermore, we found a strong positive correlation between the C2C variation strength and the amplitude of the frequency-modulation component of the Blazhko effect, indicating a connection between the two phenomena.

The metallicity of the star-forming environment is a fundamental parameter shaping the evolution of protoplanetary disks and the formation of planetary systems, yet its influence remains poorly constrained. We present a spectroscopic study of low-mass pre-main sequence (PMS) stars ($M < 1 \, M_\odot$) in the exceptionally metal-poor cluster Dolidze~25 ($Z \approx 0.2 \, Z_\odot$), using VLT/MUSE observations to probe accretion processes and disk evolution in a subsolar environment. We identify 132 cluster members using a combination of \textit{Gaia} astrometry and spectroscopic youth indicators, including lithium absorption and Balmer emission. The stellar parameters are derived using low-metallicity BT-Settl models yielding effective temperatures, extinctions, luminosities enabling robust estimates of stellar masses and ages. Mass accretion rates ($\dot{M}_\mathrm{acc}$) derived from H$\alpha$ emission span $10^{-10}$--$10^{-8} \, M_\odot\,\mathrm{yr}^{-1}$ with a median value of \(8 \times 10^{-10}\,M_\odot\,\mathrm{yr}^{-1}\). These rates are comparable to those in solar-metallicity regions of similar age, such as Lupus and Orion, indicating minimal metallicity dependence in accretion processes. Our analysis shows that using solar-metallicity templates to fit low-metallicity stars leads to systematic overestimations of \(T_\mathrm{eff}\) (by approximately \(300\,\mathrm{K}\)) and \(A_V\) (by around \(0.5\,\mathrm{mag}\)), underscoring the importance of employing metallicity-matched models for reliable characterization in low-\(Z\) environments. We present flux-calibrated, extinction-corrected spectra of these metal-poor PMS stars as a valuable resource for future investigations of disk evolution in subsolar regimes.

We investigate the stellar velocity dispersion functions (VDFs) of pressure-supported galaxies in the EAGLE cosmological simulations. The central stellar velocity dispersion is one of the fundamental dynamical tracers of the total mass of galaxy subhalos, alongside luminosity and stellar mass. Because it reflects the gravitational potential, the stellar velocity dispersion is expected to be relatively insensitive to feedback from Active Galactic Nuclei (AGN), a critical process that regulates the connection between other galaxy observables and subhalo mass. To examine the impact of AGN feedback, we analyze the VDFs from five EAGLE simulation runs, each adopting a different AGN feedback model: one "standard", two "enhanced", one "reduced", and one with no AGN feedback. We compute the stellar velocity dispersions of pressure-supported galaxies using member stellar particles, mimicking fiber spectroscopy. The VDFs from the standard and enhanced AGN feedback models show little difference. However, contrary to our initial expectation that the VDF shape would be largely insensitive to AGN feedback, the simulations with reduced and no AGN feedback show a significant excess of high velocity dispersion galaxies ($\sigma_{*}$ > 200 km s$^{-1}$) and a deficit of low dispersion galaxies (100 < $\sigma_{*}$ (km s$^{-1}$) < 200), compared to those with standard or enhanced AGN feedback. The presence of high velocity dispersion galaxies in the no-AGN model arises from enhanced central star formation, due to the absence of AGN-driven gas heating or expulsion. Our results demonstrate that the shape of the theoretical VDF is sensitive to the strength of AGN feedback. These predictions offer a theoretical benchmark for future observational studies of the galaxy velocity dispersion function using large-scale spectroscopic surveys.

Microchannel plates (MCPs) were the driving detector technology for ultraviolet (UV) astronomy over many years, and still today MCP-based detectors are the baseline for several planned UV instruments. The development of advanced MCP detectors is ongoing and pursues the major goals of maximizing sensitivity, resolution, and lifetime, while at the same time decreasing weight, volume, and power consumption. Development efforts for an MCP-based detector system for the UV are running at IAAT at the University of Tübingen. In this publication, we present our latest results towards coating aluminum gallium nitride (AlGaN) photocathodes directly on MCPs, to improve quantum detection efficiency in the far- and extreme-UV. Furthermore, we report on the implementation of a non-iterative centroiding algorithm for our coplanar cross-strip anode directly in an FPGA.

The luminosity of embedded protostars is commonly measured via observations of the dust continuum spectral energy distribution from millimetre to infrared wavelengths. However, this method cannot be applied to embedded protostars in binary or multiple systems, where their components are usually unresolved over this extended wavelength range. We propose a new method, based on the idea that a molecule formed (mainly) on the grain surfaces only emits lines in the region where it thermally sublimates from the grain mantles, heated by the photons emitted by the embedded source. In this respect, carbonyl sulfide (OCS) is an optimal molecule, because of its low binding energy and rotational lines in the millimetre. We apply the method to the protobinary system NGC1333 IRAS4A, using ALMA high-spatial resolution ($\sim$50 au) observations of the OCS(19-18) line as part of the ALMA Large Program FAUST. We also present new quantum mechanics calculations of the OCS binding energy distribution, essential for the application of the method. We found that the two binary components, A1 and A2, have a comparable luminosity within the error bars, 7.5$\pm$2.5 and 7$\pm$1 L$_\odot$, respectively. We discuss the reliability of the estimated luminosities and the potential of this new method for measuring the luminosity of embedded protostars in binary and multiple systems.

Galaxy clusters trace the densest regions of the cosmic web and are crucial laboratories for studying the thermodynamic and chemical evolution of the intracluster medium (ICM). We present a Chandra study of the massive galaxy cluster SPT-CL J0217-5014 ($z \sim 0.53$; $M_{\rm 500} \sim 3 \times 10^{14}~\rm M_{\odot}$), previously reported as a Swift serendipitous clusters with the highest Fe abundance ($\sim 1.3\pm 0.4$ $\rm Z_{\odot}$ within $\sim 1'.7$) and a potentially disturbed morphology. The X-ray morphology reveals a disturbed ICM with a surface brightness edge at $\sim 0'.26$ ($\sim 100$ kpc) to the west and a tail-like feature extending towards the east. The best-fit metal abundance within 1'.5 ($\sim 0.7\rm R_{500}$) is $0.61_{-0.23}^{+0.26}~\rm Z_{\odot}$. The derived central electron number density, entropy, and cooling time classify this system as a non-cool-core cluster, suggesting that merger activity has likely disrupted the possible pre-existing cool core. At larger radii ($\sim 1' - 2'$), we detect excess X-ray emission to the south, spatially aligned with a filamentary distribution of red galaxies, indicating ongoing accretion along an intracluster filament. Based on the DESI DR9 cross-matched optical clusters and photometric redshifts, we identify three nearby, lower-mass clusters that likely trace the large-scale structures, suggesting that SPT-CL~J0217-5014 is the primary node of a dynamically active environment where past mergers and anisotropic accretion along cosmic filaments have shaped the present-day ICM.

Interstellar objects (ISOs) ejected from planetary systems carry kinematic signatures of their formation environments. The properties of these velocity distributions govern the ISOs' propagation and dynamical evolution in the Galactic potential. We investigate how planetary system architecture influences ISO production during post-gas-disc dynamical instabilities using N-body simulations. We explore the ISO production outcomes of 2461 randomly generated systems spanning total system masses of 300-800 Earth masses and multiplicities of 3-7 planets. By integrating planets embedded in a disc of test particles for 10 Myr, we find that evolving systems can be broadly divided into two distinct classes based on their initial architectures. Catastrophic systems are characterized by high multiplicities and orbitally compact configurations, or by high-mass planets in systems with large mass asymmetries. These systems eject a large fraction of their planetesimals (median 59 percent) and, depending on the ejection pathway, produce high-speed ISOs (median 2.9 km/s). In contrast, quiet systems have lower masses and multiplicities and do not undergo significant orbital rearrangement, yet still eject a median of 28 percent of planetesimals at lower velocities (median 1.6 km/s). This dichotomy points to distinct ejection pathways, involving either violent global instabilities or more gradual, diffusive processes. Overall, we find that ISO ejection velocities are typically low, on the order of a few km/s. Although ISOs subsequently experience dynamical heating as they orbit the Galaxy, their velocity distributions retain signatures of their progenitor systems' architectures and histories, underscoring the potential use of ISOs in Galactic archaeology.

A central question for understanding interplanetary coronal mass ejection (ICME) physics and improving space weather forecasting is how ICMEs evolve in interplanetary space. We have updated one of the most comprehensive in situ ICME catalogs to date, which now includes 1976 events from 11 space missions covering over 34 years, from December 1990 to August 2025. We have combined existing catalogs including magnetic obstacles and identified and added boundaries of an additional 807 (40.8%) events ourselves. With this catalog, we demonstrate the most extensive analysis to date of total ICME magnetic field values as a function of heliocentric distance. Parker Solar Probe has observed 6 ICMEs at $< 0.23$ au (until April 2025), and Solar Orbiter and BepiColombo have added more events near 0.3~au, bridging the major observational gap towards the solar corona. Our main result is that a single power law can describe the evolution of the mean total magnetic field (exponent value of $k=-1.57$) and maximum field ($k=-1.53$) in the magnetic obstacle (MO), from 0.07 to 5.4 au. Extending the power law to the solar photosphere reveals a strong inconsistency with magnetic field magnitudes observed in the quiet Sun and active regions by 2 and 4 orders of magnitude, respectively. We introduce a multipole-type power law with two exponents, $k_1=-1.57$, and $k_2=-6$, relating the ICME magnetic field magnitude to an average solar active region field strength. These results present important observational constraints for the evolution of ICMEs from the Sun to the heliosphere.

The kinematic evolution of hypercompact H II (HC H II) regions around young high-mass stars remains poorly understood due to complex interactions with parental environs. We present ALMA QUARKS/ATOMS 1.3 mm/3 mm observations (the highest resolution $\sim0.01$ pc) of a deeply embedded HC H II region (diameter $\sim0.015$ pc, electron density $\sim2\times10^{5}$ cm$^{-3}$) exhibiting a striking $\gtrsim20$ km s$^{-1}$ global redshift seen in optically thin H30$\alpha$/H40$\alpha$ recombination lines relative to its parental hot molecular core within a hub-filament system. The 1.3 mm continuum data reveal a distinct 0.1-pc arc and a perpendicular 0.04-pc tail. We propose that this morphology arises from a dynamic champagne flow: the slow expansion of HC H II region into a pre-existing filament forms the arc and associated low-velocity (few km s$^{-1}$) SiO shocks. Meanwhile, in the opposite direction ionized gas escapes along a steep density gradient traced by the tail and high-velocity (20 km s$^{-1}$) SiO emission. We reject the bow shock scenario in which ionized gas co-moves with a runaway high-mass star because shocked gas in the arc aligns with the hub velocity, contradicting the bow shock prediction. Non-LTE radiative transfer modeling further rules out infall of ionized gas as the velocity shift origin. We conclude that this exceptional HC H II region is undergoing a few-thousand-year transition phase of "hatching out of the egg": the ionized gas of HC H II region has just broken out of its parental hot core and now is flowing outward supersonically. This work highlights how anisotropic density distributions induce supersonically anisotropic ionized flows that govern HC H II region evolution.

Recently, a class of dormant black hole binaries with luminous companions (dBH-LC) has been observed, such as $Gaia$ BH1, BH2, and BH3. Unlike previously discovered X-ray BH binaries, this type of dBH-LC has relatively long orbital periods (typically more than several tens to a few hundred days) and shows very weak X-ray emission. Therefore, studying the formation and evolution of the whole dBH-LC population is also a very interesting problem. Our aim is to study the contribution of massive stars to the dBH-LC population under different evolutionary models (isolated binary evolution (IBE) and hierarchical triple evolution), and different formation channels (such as mass transfer, common envelope evolution). Using the Massive Objects in Binary Stellar Evolution (MOBSE) code, the Triple Stellar Evolution (TSE) code, and the latest initial multiple-star distributions, we model the populations of massive stars. Finally, we calculate the orbital properties, mass distributions, and birthrates of the BH-LC populations formed under these different conditions. In the Milky Way, we calculate that the birthrate of dBH-LC formed through IBE is about 4.35$\times$$10^{-5}$ ${\rm yr}^{-1}$, while the birthrate through triple evolution is about 1.47$\times$$10^{-3}$ ${\rm yr}^{-1}$. This means that the birthrate from triple evolution is one to two orders of magnitude higher than that from IBE. We find that in triple evolution, the main formation channel of dBH-LC is post-merger binaries formed from inner binary mergers triggered by von Zeipel$-$Lidov$-$Kozai oscillations.

We present a population-scale catalogue of 287,872 supermassive black hole masses with high accuracy. Using a deep encoder-decoder network trained on optical spectra with reverberation-mapping (RM) based labels of 849 quasars and applied to all SDSS quasars up to $z=4$, our method achieves a root-mean-square error of $0.058$\,dex, a relative uncertainty of $\approx 14\%$, and coefficient of determination $R^{2}\approx0.91$ with respect to RM-based masses, far surpassing traditional single-line virial estimators. Notably, the high accuracy is maintained for both low ($<10^{7.5}\,M_\odot$) and high ($>10^{9}\,M_\odot$) mass quasars, where empirical relations are unreliable.

The core collapse of a massive star at the end of its life can give rise to one of the most powerful phenomena in the Universe. Because of violent mass motions that take place during the explosion, core-collapse supernovae have been considered a potential source of detectable gravitational waveforms for decades. However, their intrinsic stochasticity makes ineffective the use of modelled techniques such as matched filtering, forcing us to develop model independent technique to unveil their nature. In this work we present MUSE pipeline, which is based on a classification procedure of the time-frequency images using a Convolutional Neural Network. The network is trained on phenomenological waveforms that are built to mimic the main common features observed in numerical simulation. The method is finally tested on a representative 3D simulation catalog in the context of Einstein Telescope, a third generation GW telescope. Among the three detector geometries considered here, the 2L with a relative inclination of $45^\circ$ is the one achieving the best results, thus being able to detect a Kuroda2016-like waveform with an efficiency above $90\%$ at 50 kpc.

We present the design and testing of a prototype multiplexing kinetic inductance detector (KID) readout electronics for the PRobe far-Infrared Mission for Astrophysics (PRIMA) space mission. PRIMA is a Probe-class astrophysics mission concept that will answer fundamental questions about the formation of planetary systems, the co-evolution of stars and supermassive black holes in galaxies, and the rise of heavy elements and dust over cosmic time. The readout electronics for PRIMA must be compatible with operation at Earth-Sun L2 and capable of multiplexing more than 1000 detectors over 2.5 GHz bandwidth while consuming around 30 W per readout chain. The electronics must also be capable of switching between the two instruments, which have different readout bands: the hyperspectral imager (PRIMAger, 2.6-4.9 GHz) and the spectrometer (FIRESS, 0.4-2.4 GHz). The PRIMA readout electronics use high-heritage SpaceCube digital electronics with a build-to-print SpaceCube Mini v3.0 board using a radiation-tolerant Kintex KU060 field programmable gate array (FPGA) and a custom high-speed digitizer board, along with RF electronics that provide filtering and power conditioning. We present the driving requirements for the system, as well as the hardware, firmware, software, and system-level design that meets those requirements.

High-energy irradiation in X-rays and UV (XUV) can transform the planetary atmospheres through photoevaporation and photochemistry. This is more crucial for M stars, whose habitable zones for Earth-like planets are located within a few percent of an AU. Transiting exoplanets around M dwarfs offer the opportunity to study their characteristics and habitability conditions. L 98-59 is an M3 dwarf hosting six Earth-like planets, with two of them in the habitable zone of the star. X-ray observations made in 2020 and 2021 detected significant flares above a quiescent luminosity of 4-10 x 10^26 erg/s. We present the results from two short XMM-Newton observations of L 98-59, which are part of a monitoring survey to detect long-term X-ray variability and activity cycles. In October 2024 the X-ray quiescent luminosity of the star was about 5.9 x 10^25 erg/s, and it was about 6.3 x 10^26 erg/s in February 2025. We speculate that in late 2024 the star had a minimum of activity; in 2021 the star was near a maximum of an activity cycle, and in 2025 it was at the middle of the cycle. We suggest a coarse estimate of the period of about 2 years and a peak-to-peak amplitude of about 10, which is the highest among the stars with a known X-ray cycle other than the Sun. We also infer that even the outer planet in the habitable zone, L 98-59f, is exposed to an X-ray dose between 100 and 1600 times the X-ray irradiation of the Earth in the XMM band.

Modern radio interferometers deliver large volumes of data containing high-sensitivity sky maps over wide fields-of-view. These large area observations can contain various and superposed structures such as point sources, extended objects, and large-scale diffuse emission. To fully realize the potential of these observations, it is crucial to build appropriate sky emission models which separate and reconstruct the underlying astrophysical components. We introduce aim-resolve, an automatic and iterative method that combines the Bayesian imaging algorithm resolve with deep learning and clustering algorithms in order to jointly solve the reconstruction and source extraction problem. The method identifies and models different astrophysical components in radio observations while providing uncertainty quantification of the results. By using different model descriptions for point sources, extended objects, and diffuse background emission, the method efficiently separates the individual components and improves the overall reconstruction. We demonstrate the effectiveness of this method on synthetic image data containing multiple different sources. We further show the application of aim-resolve to an L-band (856 - 1712 MHz) MeerKAT observation of the radio galaxy ESO 137-006 and other radio galaxies in that environment. We observe a reasonable object identification for both applications, yielding a clean separation of the individual components and precise reconstructions of point sources and extended objects along with detailed uncertainty quantification. In particular, the method enables the creation of catalogs containing source positions and brightnesses and the corresponding uncertainties. The full decoupling of sky emission model and instrument response makes the method applicable to a wide variety of instruments or wavelength bands.

We aim to investigate the interplay between dense molecular gas, star formation, and active galactic nucleus (AGN) feedback in the luminous infrared galaxy (LIRG) NGC 1068 at sub-kiloparsec scales. We present the HCN (4-3) and HCO$^+$ (4-3) maps of NGC 1068, obtained with JCMT as part of the Mapping the dense molecular gas in the strongest star-forming galaxies (MALATANG) project, and perform spatially resolved analyses of their correlations with infrared luminosity and soft X-ray emission. Spatially resolved relations between the luminosities of infrared dust emission and dense molecular gas tracers ($L_{\rm IR}-L'_{\rm dense}$) are found to be nearly linear, without clear evidence of excess contributions from AGN activity. The spatially resolved X-ray emission ($L^{\rm gas}_{0.5-2\,\mathrm{keV}}$) displays a radially-dependent twofold correlation with the star formation rate (SFR), suggesting distinct gas-heating mechanisms between the galaxy center and the outer regions. A super-linear scaling is obtained in galactic center regions with SFR surface density ($\Sigma_{\rm SFR}$) $>$ 8.2 $\times$ 10$^{-6}$ $M_\odot$ yr$^{-1}$ kpc$^{-2}$: log($L^{\rm gas}_{0.5-2\,\mathrm{keV}}$/erg s$^{-1}$) = 2.2 log(SFR/$M_\odot$ yr$^{-1}$) + 39.1. We further found a statistically significant super-linear correlation ($\beta = 1.34$ $\pm$ 0.86) between $L^{\rm gas}_{0.5-2\,\mathrm{keV}}$/SFR and HCN(4-3)/CO(1-0) intensity ratio, whereas no such trend is seen for HCO$^+$(4-3)/CO(1-0) or CO(3-2)/CO(1-0). These findings indicate that AGN feedback does not dominate star formation regulation on sub-kiloparsec scales, and that the excitation of dense gas traced by HCN (4-3) may be more directly influenced by high-energy feedback processes compared to HCO$^+$ (4-3) and CO (3-2).

Recent detections of very high energy (VHE; GeV-TeV) photons from gamma-ray burst (GRB) afterglows, most notably the extreme event GRB 221009A, require refined models that include realistic jet structures and complex circumburst environments. The jet's angular structure is crucial for shaping afterglow emission. Our recent work demonstrates that Gaussian jets, with their smooth angular decline, naturally produce early bright peaks for on-axis observers and delayed, softer, dimmer peaks at higher inclinations. The gradual decline suppresses excessive lateral expansion, unlike the sharp edge in top-hat jets, making Gaussian jets a compelling alternative to both top-hat and other structured-jet models. Here we implement a Gaussian structured-jet model to explain TeV afterglows from adiabatic forward shocks propagating in a wind-driven medium. We show that the TeV peak time and flux depend sensitively on jet geometry, kinetic energy, wind density, and on microphysical parameter ratios that scale the SSC component. We identify the afterglow parameter space that is favourable for detecting sub-TeV photons with the Cherenkov Telescope Array (CTA), finding that only about ten per cent of simulated TeV events exceed CTA sensitivity in a wind medium. These detections arise from near core-aligned views, with high kinetic energy and wind density, moderate initial Lorentz factor and downstream magnetic field, and a relatively large fraction of energy in nonthermal electrons. Applying this model to GRB 221009A, we perform multi-band fits including wind-modified dynamics, Klein-Nishina effects, and EBL attenuation, and find that a mildly off-axis geometry reproduces the observed X-ray and GeV-TeV light curves.

During the evolution of stars on the asymptotic giant branch (AGB), thermal pulses lead to the formation of strongly stratified layers in the outer regions of the CO core, which might lead to inversions in the chemical gradient. Such inversions would produce instabilities beyond the ones predicted by the Schwarzschild criterion and the standard use of mixing length theory (MLT). We used a set of MLT equations that consider the impact of the background chemical gradients. This extension of MLT is referred to in this work as MLT$\sharp$, to make a distinction between both prescriptions. We applied MLT$\sharp$ in tandem with the more general Ledoux instability criterion. We computed the evolution in the AGB phase and compared the chemical profiles resulting from MLT, MLT$\sharp$ and the double diffusive GNA theory. We continued the evolution through a post-AGB thermal pulse and performed a pulsational analysis of the resultant GW Vir models to asses $g$-mode pulsation periods. Finally, we tested our results with pulsation properties of known GW Vir stars derived from recent observations. We find that the much simpler MLT$\sharp$ set of equations closely reproduces the results from the GNA theory. As such, MLT$\sharp$ offers a simple way to include chemically driven convection in stellar evolution computations. Stellar evolution simulations show that Rayleigh-Taylor and thermohaline instabilities can play an important role during the TP-AGB. We obtained significantly different chemical profiles using a standard MLT approach compared to those resulting from our MLT$\sharp$ and GNA computations. Our adiabatic pulsational analysis shows that these differences in the chemical stratification leave clear mode-trapping signatures in the pulsation spectrum of the GW Vir models.

We present a systematic investigation of long-term mid-infrared (MIR) color variability in 1,718 Narrow-Line Seyfert 1 galaxies (NLSy1s) using 14-year \textit{WISE}/NEOWISE monitoring data. Through Pearson correlation analysis between photometric magnitude and color, we identify: (1) a radio-quiet NLSy1 (RQ-NLSy1) population comprising 230 bluer-when-brighter (BWB) sources, 131 redder-when-brighter (RWB) sources, and 1,323 objects showing weak or statistically insignificant color variations; and (2) a radio-loud NLSy1 (RL-NLSy1) population containing 5 BWBs, 2 RWBs, and 27 sources with weak/no color variations. Our analysis reveals that the BWB tendency strengthens significantly in galaxies with redder mean MIR colors $\rm \left<W1-W2\right>$ and lower starlight contamination. Furthermore, this color-change pattern demonstrates that the most bolometric luminous sources exhibit the most pronounced BWB behavior. While similar trends exist for black hole mass and Eddington ratio, bolometric luminosity appears to be the primary physical driver. Potential origins of these variations (e.g., host galaxy contribution, accretion disk variability, and dust reprocessing) are discussed. We conclude that temperature-dependent dust reprocessing dominates the observed BWB, RWB, and no/weak variation patterns. This interpretation may also apply to similar MIR color variations observed in other extragalactic MIR transients, such as tidal disruption events, ambiguous nuclear transients, and changing-look AGNs. In addition, we find no significant difference in long-term MIR color variations between RL-NLSy1s and RQ-NLSy1s, however, RL-NLSy1s show significantly greater dispersion in intrinsic variability amplitude compared to RQ-NLSy1s due to jet-induced complexity, where non-thermal synchrotron emission from relativistic jets obscures thermal dust signatures.

Observations have revealed evidence of photochemical processing in protoplanetary disks. This processing occurs in the photon dominated layer, the optically thin regions of the disk high above the disk midplane. It remains unclear, however, how much this photochemical processing impacts the compositions of the planets and their building blocks within the disk. Here we use the oxygen isotopic compositions of Solar System solids, which has been attributed to photochemistry in the solar nebula, to quantitatively evaluate whether this processing could have produced the conditions needed to provide the diversity of compositions seen in the Solar System. We do this by modeling the chemical evolution while fine dust grows into the building blocks of the planets. We find that the oxygen isotopic evolution cannot be attributed to processing in the solar nebula and must instead be inherited from the parent molecular cloud. Further, our results indicate that the observed photochemical processing in protoplanetary disks does not significantly impact the compositions of planets that form within.

The Vera C. Rubin Observatory's 10-Year Legacy Survey of Space and Time (LSST) is expected to produce a hundredfold increase in the number of transients we observe. However, there are insufficient spectroscopic resources to follow up on all of the wealth of targets that LSST will provide. As such it is necessary to be able to prioritise objects for followup observations or inclusion in sample studies based purely on their LSST photometry. We are particularly keen to identify tidal disruption events (TDEs) with LSST. TDEs are immensely useful for determining black hole parameters and probing our understanding of accretion physics. To assist in these efforts, we present the Many Artificial LSST Lightcurves based on the Observations of Real Nuclear transients (MALLORN) data set and the corresponding classifier challenge for identifying TDEs. MALLORN comprises 10178 simulated LSST light curves, constructed from real Zwicky Transient Facility (ZTF) observations of 64 TDEs, 727 nuclear supernovae and 1407 AGN with spectroscopic labels using Gaussian process fitting, empirically-motivated spectral energy distributions from SNCosmo and the baseline from the Rubin Survey Simulator. Our novel approach can be easily adapted to simulate transients for any photometric survey using observations from another, requiring only the limiting magnitudes and an estimate of the cadence of observations. The MALLORN Astronomical Classification Challenge, launched on Kaggle on 15/10/2025, will allow competitors to test their photometric classifiers on simulated LSST data to find TDEs and improve upon their capabilities prior to the start of LSST.

The distance to the Vela Junior supernova remnant (RX J0852.0-4622 or G266.2-1.2) has long remained uncertain, limiting our understanding of its physical properties. Using VLT/MUSE integral field spectroscopy, we uncover chemical and kinematic connections between the nebula surrounding its Central Compact Object (CXOU J085201.4-461753) and the nearby Herbig-Haro outflow of Ve 7-27 (Wray 16-30), indicating a shared nitrogen-rich, Fe-peak-enhanced environment. This link ties stellar birth and death, with the young star Ve 7-27 embedded in material expelled by Vela Junior's massive progenitor, and the remnant's blast wave is expanding through the same medium. Adopting the Gaia-based distance to Ve 7-27, we revise Vela Junior's distance to $1.41\pm0.14$ kpc. At this distance, the remnant's physical radius is $23.3\pm2.3$ pc, and X-ray proper motions of the northwestern rim correspond to shock speeds of $(2.8\pm0.7)\times10^3$ to $(5.6\pm1.5)\times10^3$ km s$^{-1}$. These imply an age of $\sim$1.6-3.3 kyr and a very low ambient density, indicating that Vela Junior is expanding within a highly rarefied wind-blown cavity carved by a massive progenitor -- consistent with the non-detection of strong thermal X-ray emission. This distance update also resolves long-standing inconsistencies, with major implications for its energy budget, particle acceleration efficiency, and compact object evolution.

We investigate the post-glitch rotational evolution of pulsars by analyzing the 2007 giant glitch of PSR J1718$-$3718 using a vortex creep model that incorporates both inward and outward nonlinear vortex motion, along with a time-varying external torque. A comprehensive fitting framework is developed, constrained by prior knowledge of moment of inertia participation from previous glitch studies. We apply a Markov Chain Monte Carlo approach to quantify uncertainties and parameter correlations. The model reproduces the observed timing data and yields physically consistent values for moment of inertia fractions and creep timescales. Our results indicate that inward creep and a long-term change in external torque dominate the observed increase in spin-down rate, pointing to structural changes within the star-likely triggered by a crustquake that initiated both vortex motion and a change in the moment of inertia. We estimate that the glitch involved approximately $2.4 \times 10^{12}$ inward-moving vortices and $\sim 142$ crustal plates with a typical size of $\sim 0.03$ km. This study demonstrates that detailed post-glitch modeling of sparse timing data can simultaneously constrain internal superfluid dynamics and external torque evolution, providing a quantitative framework to probe the structural properties of neutron star interiors.

We investigate the formation of primordial black holes (PBHs) in inflationary scenarios featuring an ultra-slow-roll (USR) plateau, focusing on two coexisting production channels: PBHs originating from relic vacuum bubbles where the inflaton got trapped on the plateau, and PBHs arising from standard adiabatic density perturbations. From detailed numerical simulations we find that the bubbles are generically surrounded by type-II curvature fluctuations. Special attention is given to the distribution of initial conditions, including the relevant mean profiles and shape dispersion around them. For the adiabatic channel, we extend the logarithmic template formula $\zeta[\zeta_G]$, which maps the Gaussian curvature perturbation to the fully non-Gaussian one while incorporating mode evolution, and we compare this with numerical results obtained using the $\delta N$ formalism. While the template departs from numerical results near its logarithmic divergence, it still provides accurate threshold values for PBH formation in the parameter range relevant to our analysis. Finally, we compute the PBH mass functions for both channels. We find that the adiabatic channel dominates over the bubble-induced channel by a factor $\sim \mathcal{O}(10-10^{2})$, and that both contributions are largely dominated by the mean profiles.

Radiation pressure on spectral lines is a promising mechanism for powering disc winds from accreting white dwarfs (AWDs) and active galactic nuclei (AGN). However, in radiation-hydrodynamic simulations, overionization reduces line opacity and quenches the line force, which suppresses outflows. Here, we show that small-scale clumping can resolve this problem. Adopting the microclumping approximation, our new simulations demonstrate that even modest volume filling factors ($f_V \sim 0.1-0.01$) can dramatically increase the wind mass-loss rate by lowering its ionization state -- raising $\dot{M}_{\rm wind}$ and yielding $\dot{M}_{\rm wind}/\dot{M}_{\rm acc}\!\gtrsim\!10^{-4}$ for such modest filling factors. Clumpy wind models produce the UV resonance lines that are absent from smooth wind models. They can also reprocess a significant fraction of the disc luminosity and thus dramatically modify the broad-band optical/UV SED. Given that theory and observations indicate that disc winds are intrinsically inhomogeneous, clumping offers a physically motivated solution. Together, these results provide the first robust, self-consistent demonstration that clumping can reconcile line-driven wind theory with observations across AWDs and AGNs.

Much has been learned about Galactic cosmic rays in the past decade: On the observational side, the spectra of cosmic ray nuclei have been directly measured with high precision, resolving chemical composition up to TV rigidities. At even higher rigidities, direct detection is making contact with indirect observations from air shower arrays. A number of breaks have been found in the nuclear spectrum, which was previously thought to be a pure power law up to the knee. Data from air shower arrays also show interesting features in the arrival directions of cosmic-ray nuclei. On the theoretical side, more sophisticated models are able to explain the various spectral breaks either with transitions between different classes of sources or with changes in the transport regime. Yet, it has become clear that our ignorance of the structure of the Galactic magnetic fields, both on large and small scales, is limiting precision predictions. Turning this problem into an opportunity though, we can use Galactic cosmic rays as a laboratory for the study of Galactic magnetic fields. In this review talk, delivered at the 39th International Cosmic Ray Conference (ICRC2025), I have summarised what is known about the spectrum and anisotropies of Galactic cosmic rays, what is not known yet and what can be learnt in the future.

We present the first spatially resolved investigation of near-infrared coronal lines in dwarf galaxies hosting active galactic nuclei (AGN), using JWST/NIRSpec integral field spectroscopy. Coronal lines (CLs), which are forbidden transitions from highly ionized species with ionization potentials up to 450 eV, act as sensitive tracers of the AGN ionizing continuum and feedback processes. Across four dwarf galaxies with ionized gas outflows traced by the optical [O III] lines, we report the detection of 16 unique species of near-infrared CLs. Line ratio diagnostics indicate that photoionization from the AGN dominates the excitation of CLs. We find that the coronal line region in dwarf galaxies, traced by the various CLs, extends up to 0.5 kpc, and can constitute up to 10% of their host galaxy size. Correlations between CL luminosities and [O III] ionized gas outflow properties are consistent with a scenario in which AGN-driven outflows likely facilitate the detection of CLs and contribute to their extent. Several CLs, including [Si VI], [Si VII], and [Mg VIII], exhibit a secondary broad component (W$_{80}$ > 300 km/s). If we interpret this spatially compact gas as part of an outflow, this would indicate that the outflowing gas includes a wide range of ionizations. The estimated energetics imply this highly ionized component is compact yet powerful enough to perturb gas in the central regions of the host dwarfs. These results indicate that AGN in low-mass galaxies may produce outflows capable of influencing their structure and evolution.

Context. While supernova remnants (SNRs) are observed to produce up to 1 M$_\odot$ of dust, the amount of dust destroyed by the forward shock (FS) is poorly constrained, raising the question whether they are net dust producers or destroyers. Aims. We aim to estimate the dust destruction efficiency of SNR FSs in a realistically turbulent interstellar medium (ISM) during their most destructive phase, and assess dust shielding by high density filaments during this period. Methods. We run 3D turbulence simulations for different turbulent Mach numbers (0-3) and average ISM densities (1-100 cm$^{-3}$) to resemble observations of the turbulent ISM. We then set off a supernova to trace its 3D magnetohydrodynamical evolution for 10 kyr. Finally, we run post-processing simulations to study the dust transport and destruction by the SNR FS, considering gas and plasma drag, kinetic and thermal sputtering, and grain-grain collisions, and either silicate or carbonaceous dust. Results. The dust destruction rate of the FS strongly depends on the average ISM density and turbulence strength, varying between 27-92% (0.85-11.0 M$_\odot$) in the studied 10 kyr. Overall, dust is less efficiently destroyed in a low density medium (1 cm$^{-3}$, 27-57%) than in intermediate (10 cm$^{-3}$, 46-92%) and high densities (100 cm$^{-3}$, 73-87%). The FS destroys 8-34% less dust in high Mach turbulence compared to a homogeneous medium. Furthermore, carbonaceous grains are more robust (up to 21% more) than silicates. Conclusions. Filaments can partly shield dust from destruction in the first 10 kyr, however, always more than 0.85 M$_\odot$ of dust is destroyed, making most SNRs dust sinks under the conditions explored in this work. The destruction efficiency of the SNRs with less than 1 M$_\odot$ of destroyed dust has not yet plateaued so that they are most likely also net dust destroyers by the end of their lifetime.

In this work we derive analytic expressions and numerical recipes for finding the effective observed position of sources close enough on sky that their Point Spread Functions (PSF), modelled as Gaussian profiles, overlap. In particularly we derive these for an elongated PSF, with a long and short axis, such as we would see from an instrument with a rectangular or elliptical mirror (relevant, for example, for the Gaia mission). We show that in this case the problem can be reduced to a one dimensional brightness profile with extrema along the line connecting the two sources, with an effective PSF width that depends on the relative orientation of the PSF and its degree of elongation. The problem can then be expressed in units of this effective width to be a function of the relative separation and light ratio alone (thus reducing to a rescaling of the un-elongated case). We derive the minimum light ratio, for a given separation and effective width, above which two sources will be resolved. We map out numerical procedures for finding the positions of these extrema across all possible cases. Finally we derive the positional offset and deviance associated with observing a fixed pair of blended sources from a variety of orientations, showing that this can be a significant source of excess noise.

Strong gravitational lensing images are subject to shape distortions due to inhomogeneities along the line of sight. The leading order shape distortion is shear, which, if measurable, will be a complementary cosmological probe to traditional cosmic shear. In Hogg et al. (2025a), we modelled 23 of the SLACS strong lenses, studying the line-of-sight (LOS) shear under a variety of shear and mass model parametrisations. In this work, we model 27 additional lenses, extending our sample of LOS shear constraints to 45 in total. We find a mean shear magnitude of $0.11\pm 0.024$, showing that a significant fraction of the lenses modelled in this work possess LOS shears with unexpectedly large magnitudes, $|\gamma_{\rm LOS}| > 0.1$, even when an octupolar distortion is included in the lens mass. We further investigate if factors such as lens and source redshift, filter and PSF, or flux and signal-to-noise ratio in the lensed arcs correlate with shear. We find that none of these features play a statistically significant role in the production of unusually large shear magnitudes.

Scattering between dark matter and baryon can cool intergalactic medium temperature and deepen the 21cm signal. Such interactions have been proposed to explain the unusually deep 21cm absorption signal reported by EDGES in 2018. This paper explores the potential to detect dark matter - baryon scattering with the Hongmeng project, an upcoming lunar orbiting satellite experiment dedicated to measuring global 21cm signal between redshifts $11-46$. We self-consistently forward-model the simulated sky temperature data, jointly varying both the astrophysical and foreground models. We show that even with a very conservative observational strategy in which the experiment only takes data when both the Earth and the Sun are shielded by the Moon, Hongmeng can tighten the current constraints on dark matter - baryon scattering cross-section $\sigma_0$ by a factor of 21 after the full mission which lasts for five years. The prospective upper limit on $\sigma_0$ can reach $4 \times 10^{-43} {\rm cm^2}$ for dark matter masses between 0.1 MeV and 0.4 GeV. Even after only one month of observation, a factor of 3 improvement relative to current $\sigma_0$ limits can be expected.

Estimating physical properties of galaxies from wide-field surveys remains a central challenge in astrophysics. While spectroscopy provides precise measurements, it is observationally expensive, and photometry discards morphological information that correlates with mass, star formation history, metallicity, and dust. We present a conditional flow matching (CFM) framework that leverages pixel-level imaging alongside photometry to improve posterior inference of galaxy properties. Using $\sim10^5$ SDSS galaxies, we compare models trained on photometry alone versus photometry plus images. The image+photometry model outperforms the photometry-only model in posterior inference and more reliably recovers known scaling relations. Morphological information also helps mitigate the dust--age degeneracy. Our results highlight the potential of integrating morphology into photometric SED fitting pipelines, opening a pathway towards more accurate and physically informed constraints on galaxy properties.

The comoving cosmic star formation rate density peaks at z~2-3, with dusty star-forming galaxies being significant contributors to this peak. These galaxies are characterized by their high star formation rates and substantial infrared luminosities. The formation mechanisms remain an open question for these galaxies, particularly with respect to how such intense levels of star formation are triggered and maintained. We aim to resolve CO(3-2) emission toward two strongly lensed galaxies, SPT0125-47 and SPT2134-50, at z~2.5-2.8 to determine their morphology and physical properties. We used high-resolution ALMA band 3 observations of CO(3-2) emission toward both sources to investigate their properties. We performed parametric and nonparametric lens modeling using the publicly available lens modeling software PyAutoLens. We divided the CO(3-2) emission line into two bins corresponding to the red and blue portions of the emission line and nonparametrically modeled the source plane emission for both bins. We performed a basic analysis of the morphology and kinematics in the source plane using nonparametric lens modeling of the red and blue bins. We found tentative evidence of a velocity gradient across both sources and no evidence of any clumpy structure, companions, or ongoing mergers. The previously calculated high star formation rates and low depletion times of both SPT0125-47 and SPT2134-50 suggest that these galaxies are undergoing a dramatic phase in their evolution. Given the lack of evidence of ongoing interactions or mergers in our source plane models, we suggest that the intense star formation was triggered by a recent interaction and/or merger. We also consider the possibility that these galaxies might be in the process of settling into disks.

Determining the physical properties of galaxies during the first billion years after the big bang is key to understanding both early galaxy evolution and how galaxies contributed to the epoch of reionization. We present deep ALMA observations of the redshifted [OIII] 88um line for the gravitationally lensed ($\mu = 11.4\pm1.9$) galaxy A383-5.1 (z=6.027) that has previously been detected in [CII] 158um. Recent James Webb Space Telescope (JWST) imaging identified this sub-L* galaxy as a ''Little Red Dot'' (LRD). With a line luminosity of $L_{\rm [OIII]} = (1.29\pm0.24)\times10^8$ L$_\odot$ (corrected for lensing magnification) A383-5.1 is one of the faintest galaxies with combined [CII] and [OIII] detections. The ALMA data reveal no dust continuum emission, consistent with previous observations. The high line luminosity ratio of [OIII]/[CII] $\sim 14\pm5$ is consistent with A383-5.1 being low-metallicity and dust-poor. The non-detection of dust continuum in bands 6 and 8 is consistent with the high [OIII]/[CII] ratio and suggests a presence of a strong ultraviolet radiation field, which would be less affect by dust attenuation, implying that galaxies of this type could contribute significantly to the ionization of the intergalactic medium. The presence of strong ionizing field could provide an important piece of information for understanding the nature of LRDs and their role in cosmic reionization.

We report on new orbit-minimum photometry and revised radial-velocity fitting that provide an improved measurement of the mass of the neutron star (NS) in pulsar PSR~J0952$-$0607 at $M_NS = 2.35\pm 0.11 M_\odot$. With its fast spin and unusually low magnetic field, this NS has evidently experienced unusual evolution, likely connected with its high mass, which is now $2.5\sigma$ above that of the heaviest pulsar with a white dwarf companion, as measured by Shapiro delay techniques. By tightening the mass measurement, we also raise the maximum (commonly called Tolman-Oppenheimer-Volkoff) NS mass to $M_{\rm TOV} > 2.27\,M_\odot$$(2.12\,M_\odot)$ at $1\sigma$$(3\sigma)$ confidence, which improves bounds on the dense-matter equation of state. While the statistical error decreases and systematic issues should be modest, uncertainties remain; we comment briefly on these factors and prospects for further improvement.

We present the spectrum of $n_s$ constraints from current CMB data (Planck, ACT, SPT-3G) combined with DESI BAO data, and highlight the interplay of $n_s$ with the optical depth to reionization $\tau$. The spectral index $n_s$ of the primordial power spectrum provides a window into early universe, and constraints on $n_s$ play an important role in discriminating early universe models such as models of cosmic inflation. Historically constrained by cosmic microwave background (CMB) experiments, the constraints on $n_s$ shift upward when CMB data is combined with the latest baryon acoustic oscillation (BAO) data from the Dark Energy Spectroscopic Instrument (DESI). Recent work explained the origin of this and the relation to the BAO-CMB tension between CMB experiments and DESI BAO, and as a case study presented constraints on $n_s$ from the combination of Atacama Cosmology Telescope (ACT) DR6 data and DESI DR2 data. Here we present constraints from Planck (PR3 and PR4), ACT, the South Pole Telescope (SPT), and the combination of all three CMB experiments, CMB-SPA, with and without DESI DR2 BAO data, and with and without CMB lensing data. In all cases the constraint on $n_s$ is shifted upwards when DESI is included, with the largest shift exhibited by ACT. This is accompanied by a commensurate shift in the constraint on the optical depth to reionization $\tau$, which is again greatest for ACT. When CMB data are combined into CMB-SPA and combined with DESI the $n_s$ constraint disfavors at more than $2\sigma$ the inflation models preferred by Planck alone, such as Higgs, Starobinsky, and exponential $\alpha$-attractors, in favor of other models, such as polynomial $\alpha$-attractors. This work motivates the further study of the tension between CMB and DESI BAO data, and of the rich interplay between $n_s$ and $\tau$.

This paper examines four case studies of life-detection claims in astrobiology, covering both biosignatures and technosignatures: the 1877 "canals" on Mars, the 1976 Mars Viking landers experiments, the 2020 phosphine detection on Venus, and the 2020 Breakthrough Listen Candidate 1 (BLC1) signal. We analyse the process of discovery for each case, including how they were detected, the media reception, the ensuing scientific debate, the correction processes, and the time it took until an expert consensus was reached. We identify lessons learned while providing scientists, the scientific community, and science communicators with recommendations for approaching future claims of astrobiological discoveries. To avoid potential cognitive biases and mitigate premature conclusions, we stress the need for clear communication of uncertainties, as well as thorough debate and verification processes among the scientific community. These responsible approaches can strengthen the credibility of scientists, cultivate a supportive scientific community, and help astrobiology flourish as a field.

Before we discuss the deflection of light in a gravitational field, we give a brief overview of some basic physical formulas on photon properties, generation and propagation. The much debated problems of the redshift and the photon propagation in a gravitational field is then considered and applied to the calculation of the speed of light. Many citations are given in direct quotations to avoid any misunderstandings. If the quotations are in German, an English translation is provided. Based on this speed, calculated and measured results are recalled on the deflection of light, with emphasis on the deflection near the Sun. We conclude that the speed of light and the deflection angle can be determined by energy and momentum conservation principles.

This work proposes a direct link between the hierarchy problem and Weakly Interacting Massive Particles (WIMPs): we suggest that the small mass of the Higgs boson arises from being dynamically driven to the scale of the WIMP. Such a special electroweak vacuum is singled out by lying close to the critical boundary of a phase transition, as recently explored in a new class of cosmological solutions to the hierarchy problem. They generically predict the Higgs potential to be destabilised just above the weak scale. Intriguingly, the requirement for new physics to achieve this coincides with two independently well-motivated expectations: a split spectrum of light fermions and heavy bosons, as anticipated from naturalness, and the so-called "WIMP miracle". A WIMP with mass around the weak scale not only happens to have the correct thermal relic abundance to be the dark matter (DM), it can also give rise to the necessary critical boundary at the TeV scale through its Yukawa couplings to the Higgs. We use a higgsino-like singlet-doublet model to illustrate our Higgs-DM criticality scenario and show that if this WIMP DM mass is observed to be greater than ~1.2 TeV then it necessarily implies a strong bound on the Higgs mass and an upper bound on the scale of heavy new physics that restores vacuum stability. It can be thoroughly probed in direct detection experiments, astrophysical signals and future collider searches, further motivating a comprehensive exploration of the remaining heavy WIMP parameter space.

We investigate extreme-mass ratio inspirals (EMRIs) with orbital periods exceeding the observational timescale of mHz gravitational wave observatories. In their early, highly eccentric phases, these systems generate transient gravitational wave bursts during pericentre passages, separated by long quiescent intervals; we designate these signals ``popcorn EMRIs.'' We utilize a steady-state analytical model based on the continuity equation in phase space to estimate the population in a Milky Way-like galaxy. The normalization of this model is linked to the solution of the Fokker-Planck equation describing stellar relaxation. Adopting a conservative one-year observation baseline ($P>1$ year), we estimate the steady-state population of popcorn EMRIs. We forecast an observable burst rate of 5 to 44 events per year. The low duty cycle ($\sim 10^{-4}$) confirms their manifestation as isolated transients. Individual bursts from the Galactic Centre exhibit high detectability. Analyzing these intrinsically transient signals demands a rigorous mathematical framework, as standard windowing techniques distort burst morphology. We establish an analytical foundation using standard smoothing techniques commonly used in real analysis. This yields the mathematically correct definition for the Fourier transform of transient signals, justifying the use of the direct Fourier transform without ad hoc windowing and ensuring the integrity of spectral analysis.

We study the electroweak first-order electroweak phase transition (FO-EWPT) within the Standard Model Effective Field Theory (SMEFT) framework induced by dimension-6 operators. Such phenomena can be probed independently via \textit{di}-Higgs production at the collider experiments as well as via the detection of gravitational waves (GW). There are three dimension-6 SMEFT operators that simultaneously modify the Higgs potential at tree level and contribute to the \textit{di}-Higgs production at the hadron colliders. With \textit{di}-Higgs production being suppressed at current LHC runs, we aim to probe this production at high luminosity (HL) and high energy (HE) runs of the LHC to achieve better sensitivity of dimension-6 SMEFT operators. The correlations among these operators are analyzed in the context of probing FO-EWPT, emphasizing the complementarity between future GW observations and upgraded LHC searches.

Inflation generally assumes a field with nonzero potential that leads to inflationary expansion happening at arbitrarily early times. We demonstrate potentially observable consequences of inflation with a finite initial time in a model in five-dimensional warped anti-de Sitter space, with both a UV and an IR brane present during inflation. Considering an inflaton with an approximately flat potential localized on the UV brane, we derive the resulting brane motion in the bulk and the 4D effective action describing the dynamics. A concrete model allows us to evaluate possible consequences of a starting point of inflation. The background evolution is driven by the fast roll of the radion at early times and the slow roll of the inflaton at late times. We find that the action has the form of a two-field hyperbolic inflation model, the two fields being the radion and the inflaton, both of which have a time-dependent background solution. This setup is holographically dual to an inflaton coupled to a strongly coupled confined sector in which the confinement scale is larger than the Hubble scale, and a confinement scale whose ratio to the 4D Planck scale evolves cosmologically. Focusing on the period when the equation of state becomes that of inflation, we find that the presence of the IR brane leads to deviations from the approximate de Sitter background in addition to those from the slow-roll parameters of the inflaton potential. We quantify the effect of the presence of the IR brane on the two point function of the adiabatic scalar perturbations and tensor perturbations. The dominant deviations occur at large scales: the adiabatic power spectrum has a blue tilt, while the tensor power spectrum shows oscillatory features. We present numerical fits to the shape of the adiabatic power spectrum, and discuss the implications for cosmic microwave background (CMB) analysis.

Nucleon short-range correlations (SRCs) and their high-momentum tails (HMTs) encode key short-range dynamics in nuclei and dense matter. This review provides a concise overview of SRC features relevant to the Equation of State (EOS) of isospin-asymmetric nuclear matter. We summarize empirical and theoretical properties of the single-nucleon momentum distribution $n(k)$, emphasizing the role of the neutron--proton tensor force, the dominance of correlated np pairs, and the enhancement of minority-species HMTs. Links to nucleon effective E-masses, quasi-deuteron components, and orbital entanglement are briefly noted. We examine how SRC-induced HMTs modify kinetic and potential contributions to the EOS in both non-relativistic and relativistic frameworks, including the softening of the kinetic symmetry energy and departures from the isospin parabolic approximation of asymmetric nuclear EOS. Sensitivity to high-momentum components and generalizations to arbitrary dimensions are also highlighted. Implications for heavy-ion reactions are summarized, including effects on particle yields, collective flows, deeply sub-threshold particle production and hard photon emission, driven by modified initial nucleon momentum distributions and abundant high relative-momentum np pairs during the reaction. Finally, we outline SRC-HMT consequences for neutron-star matter, covering proton fractions, tidal deformabilities, $Z$-factors, cooling, and the core--crust transition, as well as possible connections to dark-matter interactions in dense environments.

We investigate the production of dark matter (DM) via a heavy $Z'$ mediator during the post-inflationary reheating epoch. In particular, we study production from three mechanisms which are smoothly connected to one another: WIMP-like freeze-out, FIMP-like freeze-in, and ultra-relativistic freeze-out (UFO). This is the first systematic study of $Z'$ portal DM which includes UFO. We find that much of the available parameter space for keV to TeV DM lies in the UFO regime for $ 1 \text{ TeV}\lesssim M_{Z'} \lesssim 1 \text{ PeV}$. When the mediator mass $M_{Z'}$ is greater than both the DM mass and the reheating temperature, UFO is a robust mechanism for producing cold DM. Although UFO DM is initially "hot" after freeze-out, it can easily become cold before structure formation if freeze-out occurs during post-inflationary reheating. Compared to standard freeze-in, UFO can accommodate significantly stronger interaction strengths (stronger couplings and/or smaller mediator masses).

We present a mechanism for gravitational wave (GW) production from isocurvature perturbations in spectator scalar fields during inflation. These energetically subdominant fields develop blue-tilted power spectra through inflationary dynamics, generating second-order scalar perturbations that source a stochastic GW background. The mechanism naturally satisfies CMB constraints at large scales while producing enhanced signals at smaller scales across a broad frequency range $10^{-20} - 1$ Hz. We perform comprehensive numerical and analytical calculations of the complete isocurvature spectrum evolution, including gravitational particle production, reheating dynamics, and scalar-induced GW generation. For spectator fields with effective masses $0.5 \lesssim m_{\chi,\mathrm{eff}}/H_I$, the resulting GW energy density reaches $\Omega_{\text{GW}} h^2 \sim 10^{-20}$-$10^{-12}$, accessible to pulsar timing arrays, space-based interferometers, and next-generation CMB experiments. Our analysis reveals that GW-induced constraints exceed current isocurvature bounds. We examine both unstable (curvaton-like) and stable (dark matter) spectator fields, demonstrating strong sensitivity to reheating temperature, inflaton-spectator coupling, and decay dynamics. This framework establishes isocurvature-sourced GWs as a powerful probe of early universe physics, enabling simultaneous constraints on inflationary dynamics, dark matter production, and reheating through coordinated multi-frequency GW observations.

For kinetic modeling of plasma processes in space, a rejection-sampling procedure for generating a Kappa distribution in particle-in-cell (PIC) simulation is proposed. A Pareto distribution is employed as an envelope distribution. The procedure only requires uniform variates, and its acceptance efficiency is $\approx 0.73$--$0.8$.

We study the linear stability of nonrelativistic Proca stars under generic perturbations. Using a combination of analytic and numerical methods, we demonstrate that, as expected, the ground state is always mode-stable. Additionally, we identify several mode-stable spherically symmetric excited states, including stationary states of constant and radial polarization, as well as multi-frequency states in case that the spin-spin selfinteraction vanishes. The existence of stable excited states may have implications for spin-$1$ ultralight dark matter models.

Space-based gravitational wave (GW) detectors are designed for wave sources in the millihertz band with different locations and orientations. Time-delay interferometry (TDI) technique is an indispensable ingredient in space-borne GW detection that effectively suppresses the laser phase noise. The abundant TDI solutions derived in the literature also feature distinct angular-dependent sensitivities. Because a GW source's angular location is unknown prior to the signals' detection, a solid-angle average is often performed when analyzing the sensitivity function of a given TDI combination. The present study explores the angular dependence of the detector's sensitivity. This detail is relevant, because once the initial detection is achieved, the source's location can be extracted and used to provide information on a refined TDI combination tailored for the specific GW source. As the TDI technique is a post-processing algorithm, such a procedure can be implemented in practice. We evaluate the angular dependence of the detector's response function to the GW signals for different TDI combinations as a function of the orientation angles. Moreover, we classify the response functions into seven categories at the low-frequency limit, leveraging the characteristics of the underlying geometrical TDI combinations. By further averaging out the azimuthal angle $\phi_D$ in the detector's plane, the main features of the resulting response functions and their zenithal dependence with respect to the GW source are scrutinized. The findings presented in this work provide pertinent insights for ongoing space-borne detector programs.

We unite the dual quarkyonic model with the quark-meson coupling (QMC) model to construct a novel nuclear model based on the quark degrees of freedom, which can cover a wide range of nuclear density from low density to the crossover region. In the model, the relativistic, gaussian quark wavefunction is used to describe the nucleon structure. We first evaluate the energy density, chemical potential, pressure and sound velocity within the ideal Fermi gas picture. In this case, those physical quantities are discontinuous or divergent at the quark saturation density, where quarkyonic phase emerges. To remove such singular behavior, we next introduce an infrared regulator, and combine the dual quarkyonic model and the QMC model to include the nuclear interaction -- we call it the quarkyonic quark-meson coupling (QQMC) model. We then find that the quark saturation density depends strongly on the nucleon size. For example, when $r_p = 0.6\, (0.8)$ fm, where $r_p$ is the root-mean-square radius of proton, the quark saturation density is about $3.6\,(1.5) \times \rho_0$ in symmetric nuclear matter, where $\rho_0$ is the nuclear saturation density. It is notable that the nuclear interaction is quite important to consider physical quantities quantitatively. In fact, the QQMC model can produce the sound velocity which is consistent with that inferred from the observed data of several neutron stars. Furthermore, pressure in symmetric or pure neutron matter deduced from the experiments of heavy-ion collisions at high energy can be explained by the QQMC model as well. We discuss in detail the formulation for the QQMC model and the physical quantities calculated by the model.

The detection of gravitational waves from compact binary coalescences has provided significant insights into our Universe, and the discovery of new and unique gravitational wave candidates from independent searches remains an ongoing field of research. In this work, we built a hybrid search pipeline that combines matched filtering and deep learning to identify stellar-mass binary black hole candidates from detector strain data. We first present results from a targeted injection study to benchmark the sensitivity of our method and compare it with existing search pipelines. We demonstrate that our hybrid approach has comparable sensitivity for injections with a source-frame chirp mass greater than 25$\,$M$_{\odot}$, and below this threshold our sensitivity drops off for signals with a network SNR less than 15. We also observe that our search method can identify a significant population of unique candidates. Furthermore, we conduct an offline search for gravitational wave candidates in the third observing run of the LIGO-Virgo-KAGRA Collaboration (LVK), yielding 31 candidates previously reported by the LVK with a probability of astrophysical origin $p_{\rm astro}\geq0.5$. We identify two other candidates: one previously reported only in a search conducted by the Institute for Advanced Study, and one previously unreported promising new candidate with a $p_{\rm astro}$ of 0.63. This unique candidate has a high chirp mass and a high probability that the primary black hole is an intermediate-mass black hole.

We investigate the imprints of chaos in gravitational waves from extreme-mass-ratio inspirals configuration, where a stellar massive object, confined in a harmonic potential, orbits a supermassive Schwarzschild-like black hole embedded in a Dehnen-type dark matter halo. In our first paper [1], we demonstrated the system's transition from non-chaotic to chaotic dynamics by analyzing Poincaré sections, orbital evolution, and Lyapunov exponents across different energies and dark matter halo parameters. In this work, we compute the gravitational waveforms of the small celestial object along different chaotic and non-chaotic orbits by implementing the numerical kludge scheme. We further perform a spectral analysis of the gravitational waveforms from such orbits. In particular, we show that when the system is in a chaotic state, the gravitational wave signals are characterized by broader frequency spectra with finite widths, enhanced amplitude and energy emission rate, distinctly differentiating them from the signals generated during the system's non-chaotic state. Through recurrence analysis we also show that the time series of gravitational waveforms strain carry unique information on the motion of chaotic dynamics, which can be used to distinctly differentiate from non-chaotic to chaotic motion of the source. Furthermore, we discuss the potential detectability of these orbits for upcoming observatories like LISA, TianQin, and Taiji, emphasizing the significant potential for detecting chaotic imprints in gravitational waves to substantially enhance our understanding of chaotic dynamics in black hole physics and the dark matter environments of galactic nuclei.

Relativistic black hole collisions are one of the most dramatic astrophysical events that can be imagined. They could provide the ideal condition for searching for possible new physics beyond general relativity. However, such events are presumably rare and difficult to occur under normal conditions. Black holes in a triple system can be accelerated to the relativistic limit and may harbor the chance for a relativistic collision. In this paper, we study the relativistic black hole collisions in a massive black hole background and the capabilities of several current and future gravitational wave detectors in detecting such signals.

Adaptive optics enables the deployment of interferometer-based spectroscopy without the need for moving parts necessary for scanning the interferometer arms. Here, we employ a Michelson Interferometer in conjunction with a Spatial Light Modulator (SLM) for determining the spectral profile of a narrow-band light source. Interestingly, we observe that the fringes across the interferometer output beam are inherently shifted in wavelength even when a constant phase profile is provided to the SLM. We calibrate the spectral shifts as a function of fringe spatial location by measuring the incident light spectrum at various points across the fringe pattern, and observe that the spectral peak traces out a `teardrop' shape, whose width is dependent on the spectral bandwidth of the source, the relative tilt and path difference between the two arms of the interferometer, and the divergence of the beam. Next, we demonstrate that this inherent spectral variation of the fringes can be used to perform fast single-snapshot spectroscopy of narrow-band light sources, while a time-varied phase profile provided to the SLM leads to multi-step spectroscopy with lower noise, higher resolution, and better contrast. Our findings establish that the Michelson Interferometer can be used to perform spectroscopy of any source within a certain spectral range from simple images of the fringe pattern, so as to facilitate exciting applications towards hyperspectral imaging.

The violation of the null energy condition (NEC) is essential for constructing nonsingular cosmological scenarios, such as Genesis cosmology, which avoids the initial singularity by initiating cosmic evolution from an asymptotically Minkowski state. To address theoretical concerns regarding the accumulation of negative energy, the smeared null energy condition (SNEC) has been proposed as a quantum-motivated, semi-local bound on NEC violation. In this work, we examine the implications of the SNEC conjecture for Genesis models, typically constructed within generalized Galileon theories. Our results demonstrate that SNEC imposes nontrivial restrictions on the viability of Genesis models, highlighting the SNEC conjecture as a powerful tool for constraining nonsingular cosmological scenarios.

We show that stationary black holes and the evolving universe belong to fundamentally different thermodynamic regimes: black holes obey ordinary Bekenstein Hawking thermodynamics, whereas cosmology necessarily follows memory-bearing teleodynamics. We show that teleodynamics is not valid for black holes, but is unavoidable in an expanding cosmology. This provides a dynamical, semi-classical realization of the thermodynamic split conjecture and identifies memory accumulation as the natural source of deviations from the area law in cosmology. Our results suggest that quantum gravity should not seek to extrapolate black hole thermodynamics to the universe, but instead must incorporate horizon memory as a fundamental microscopic ingredient and consider cosmological constructions consistent with that.

Considering the earth or any other celestial body the main sources of the gravitational as well as of the magnetic field lie inside the body. Above the surface both fields are in good approximation harmonic vector fields determined by their values at the body's surface or any other surface enclosing the body. The intensity problem seeks to determine harmonic vector fields vanishing at infinity and with prescribed intensity of the field at the surface. This problem constitutes a nonlinear boundary value problem, whose general solvability is not yet established. In this paper {\em axisymmetric} harmonic fields ${\bf H}$ outside the unit sphere $S^2$ are studied and, given an axisymmetric Hölder continuous intensity function $I\neq 0$ on $S^2$, the existence of infinitely many solutions of the intensity problem is proved. These solutions can more precisely be characterized as follows: fix a number $\de \in \nat\setminus \{1 \}$ and a meridional plane $M$ through the symmetry axis $S\!A$, and in $M$ a unit circle $S^1$ (symmetric with respect to $S\!A$) and, furthermore, $2\, N$, $N \in \nat_0$, points $z_n \in M$ (symmetric with respect to $S\!A$, avoiding $S\!A$, and outside $S^1$), then the existence of an (up to a sign) unique harmonic field ${\bf H}$ is established that vanishes at (the axisymmetric circles piercing $M$ at) $z_n$ and nowhere else, that has intensity $I$ at $S^2$ and (exact) decay order $\de$ at infinity. The proof is based on the solution of a nonlinear elliptic equation with discontinuous coefficients, which are, moreover, singular at the symmetry axis. Its combination with fixed boundary conditions was the basis of a recent treatment of the ``geomagnetic direction problem'' \cite{KR22}. Here we have instead natural boundary conditions, which provide less information, and which require, therefore, in part new solution techniques and sharper estimates.

We present results from long-term, numerical-relativity simulations of binary neutron star mergers modeled using both, fully tabulated, finite-temperature, equations of state and their corresponding hybrid representations. The simulations extend up to 150 ms which allows us to assess the role of the treatment of finite-temperature effects on the dynamics of the hypermassive neutron star remnant. Our study focuses on the analysis of the spectra of the post-merger gravitational-wave signals and on how these are affected by the treatment of thermal effects in the two EOS representations. Our simulations highlight distinct differences in the GW frequency evolution related to the thermal modeling of the EOS, demonstrating that deviations from established quasi-universal relations become significant at late post-merger phases. Furthermore, we investigate the stability of the HMNS against convection. Employing both the Ledoux criterion, necessary condition for the development of convective instabilities, and the Solberg-Høiland criterion, a generalized criterion for axisymmetric perturbations based on a combined analysis of the Brunt-Väisälä frequency and of the epicyclic frequency, we show that differential rotation and thermal stratification in the HMNS give rise to local (yet sustained) convective patterns that persist beyond 100 ms after merger. Those convective patterns, while substantially different between tabulated and hybrid EOS treatments, trigger the the excitation of inertial modes with frequencies smaller than those attained by the fundamental quadrupolar mode, and are potentially within reach of third-generation GW detectors. The late-time excitation of inertial modes, previously reported in studies based on hybrid EOS, is fully supported by the tabulated, finite-temperature EOS simulations presented here, which account for thermal effects in a more consistent way.