We investigate whether low-frequency pulsar observations can provide LoS magnetic field estimates and whether these are consistent with synthetic LoS signatures extracted from a three-dimensional CME reconstruction constrained by Solar Orbiter data. We analyze a CME occultation of the LoS to PSR J1022+1001 on 20 August 2021, observed simultaneously with LOFAR and NenuFAR. From LOFAR, we derive time-resolved dispersion measure (DM) and rotation measure (RM) and isolate the CME contributions using background estimates for interstellar, solar wind and ionospheric components. We then infer the density-weighted LoS-averaged magnetic field component <B||>_PSR from the ratio delta-RM/delta-DM. In parallel, we reconstruct the CME using a semi-empirical 3DCORE model fitted to Solar Orbiter in-situ magnetic field observations at 0.65 au. We sample the modeled magnetic field along the pulsar LoS using fixed spatial sampling points and compute synthetic LoS-averaged signatures <B||>_3D for different flux rope configurations. The derived <B||>_PSR increases from approximately -9 nT to a peak near 63 nT during the observed interval. Comparison with synthetic signatures shows that the polarity and temporal evolution of the LoS signal are strongly dependent on the flux rope configuration and only a South-West-North (SWN) configuration (confirmed by Solar Orbiter in-situ data) reproduces the observed sign and overall evolution, whereas alternative configurations are incompatible. The modeled amplitudes, however, are systematically larger than the pulsar-derived values by roughly a factor of five. We show that simultaneous low-frequency pulsar DM and RM measurements can provide LoS magnetic field estimates for a CME and can be used to test CME magnetic structure against data-constrained three-dimensional reconstructions.

We simulate Venus' evolution with a coupled one-dimensional solar-atmosphere-lithosphere-mantle-core model to predict currently unobservable features and its eruptive mass flux. We identified four distinct evolutionary pathways that simultaneously match the atmospheric abundances of water and carbon dioxide as well as the lack of a core dynamo. These scenarios are characterized by I) generally monotonic cooling, II) a low mantle melt fraction in which Venus' volcanically active phase is ending, III) a small inner core, and IV) oscillations of internal properties. Through random forest classification we determined that the key parameters that distinguish these types are the initial mantle water abundance, the mantle viscosity, the dehydration stiffening strength, the eruption efficiency, and the melting point of the core. In each of the plausible histories, Venus retains at least one Earth ocean's worth of water in its mantle and remains volcanically active today. Venus' lack of a current geodynamo allows thermal histories with an initially large inner core in our parameter sweep. In 88% of plausible histories we found that Venus possessed a past magnetic field. The results strongly disfavor recent high eruption rate estimates, but are consistent with lower estimates. Current resurfacing estimates also strongly disfavor the low melt scenario, implying that Venus is not nearly volcanically ``dead.'' These predictions are testable with anticipated data and the model can be applied to exoplanets to predict their properties.

Extreme mass-ratio inspirals (EMRIs), comprising a stellar-mass compact object (CO) orbiting a supermassive black hole (BH), are key targets for future space-based gravitational-wave (GW) observatories. Incorporating the spin of the secondary body into waveform models not only enhances measurement precision but also offers insight into the spin distribution of stellar-mass objects. In this work, we construct the flux and waveform for an EMRI with a spinning secondary in a Kerr background under the linear-spin approximation. Using the radiative prescription (half-retarded minus half-advanced field), we derive orbit-averaged evolution equations for the fundamental constants of motion, including the energy, angular momentum, Carter-like constant, and the parallel spin component. This framework provides a tractable route to generating waveforms that incorporate the secondary spin, with the potential for further simplification in future work.

We report the discovery of bimodal structure in the drift rate distribution of upward-drifting burst clusters from the hyperactive repeating fast radio burst FRB 20240114A. Using unsupervised machine learning (UMAP dimensionality reduction combined with HDBSCAN density-based clustering) applied to 233 upward-drifting burst clusters from the FAST telescope dataset, we identify a distinct subpopulation of 45 burst clusters (Cluster C1) with mean drift rates 2.5x higher than typical upward-drifting burst clusters (245.6 vs 98.1 MHz/ms). Gaussian mixture modeling reveals strong evidence for bimodality (delta-BIC = 296.6), with clearly separated modes (Ashman's D = 2.70 > 2) and a statistically significant gap in the distribution (11.3 sigma). Crucially, we demonstrate that this bimodality persists when restricting the analysis to single-component (U1) burst clusters only (delta-BIC = 19.9, Ashman's D = 2.71), confirming that the result is not an artifact of combining single- and multi-component burst clusters with different drift rate definitions. The extreme-drift subpopulation also exhibits systematically lower peak frequencies (-7%), shorter durations (-29%), and distinct clustering in multi-dimensional feature space. These findings are suggestive of two spatially separated emission regions in the magnetosphere, each producing upward-drifting burst clusters with distinct physical characteristics, although confirmation requires observations from additional epochs and sources.

In this study, the behavior of high-frequency quasi-periodic oscillations (QPOs) is investigated around a Kerr black hole immersed in a uniform Bertotti-Robinson magnetic field. The motion of the test particle is analyzed by determining the geodesic equations and evaluating the corresponding orbital, radial, and vertical epicyclic frequencies. These fundamental frequencies are used to construct the theoretical framework of QPO models based on parametric and forced resonance mechanisms. Observational data obtained from several black hole X-ray binaries (GRO J1655-40, XTE J1550-564, XTE J1859+226, GRS 1915+105, H1743-322, M82~X-1, and Sgr~A$^{*}$) are used to constrain the black hole parameters through Bayesian inference and Markov Chain Monte Carlo (MCMC) analyses. For the X-ray binaries GRO J1655-40, GRS 1915+105, H1743-322, and M82~X-1, nonzero values of the dimensionless parameter $b=Bm$ are obtained at the $68\%$ confidence level within the framework of the parametric resonance model, while only upper bounds at the $90\%$ confidence level are obtained for the remaining sources. In contrast, in the case of the forced resonance model, only an upper bound at the $90\%$ confidence interval is obtained for the magnetic field parameter for all considered X-ray binary sources. The analysis indicates that the value of the magnetic field parameter is small but not negligible, producing minor modifications to particle dynamics and epicyclic frequencies. The influence of the magnetic field is further examined through the properties of the innermost stable circular orbit and the radiative properties of the thin accretion disk, including the energy flux and temperature profiles, within the allowed parameter range inferred from the MCMC analysis.

We present a simplified and general description of the high-redshift information in acoustic scale measurements from the cosmic microwave background and large-scale structure. The transverse distance interval between photon--baryon decoupling and a late epoch in the matter era provides an analytically tractable summary statistic thereof and a general diagnostic of the current tension between the Dark Energy Spectroscopic Instrument and the CMB. We show that this "matter-era distance excess" is unlikely to be explained by modified dynamics at low redshift. We then analytically derive the matter-era distance interval's sensitivity to new physics at high redshift, including nonstandard recombination, nonminimal dark matter dynamics, and spatial curvature; in particular, we explain how this observable represents a direct geometric measurement of (and underlies the current incompatibility with) neutrino masses. Finally, we demonstrate that phenomenological models of dynamical dark energy mediate the matter-era distance excess in a manner reliant on their unphysical, extrapolated behavior at high redshift. Invoking alternative explanations of the excess removes the CMB's contribution to the evidence for these models; the residual preference of around $1.7\sigma$ mostly derives from DESI's two lowest-redshift measurements of the Alcock--Paczynski distortion, without which it drops to $0.5 \sigma$.

Mass discrepancies in galaxies are empirically known to appear only below a characteristic acceleration scale a0. Here we show that this behaviour is not limited to galaxies: it extends continuously across the full hierarchy of self-gravitating stellar systems, from gas-rich dwarfs and spirals to massive early-type galaxies, and further down to compact stellar clusters. We introduce the Milgromian dynamics (MOND) depth index DM, together with dynamical maturity index T = tcross/tH, dynamical collisionality index T1 = tcross/trelax, with tcross being the crossing time, tH the Hubble time and trelax the median two-body relaxation time, and the MOND acceleration index A = abar/a0. We uncover a well-defined two-dimensional dividing surface in dynamical space. The "dark matter phenomenon" is found only in systems that are both in the deep-MOND regime (abar < a0) and collisionless (trelax > tH), while high-acceleration, collisional systems (abar > a0, trelax << tH), including globular clusters and UCDs, show no evidence for a mass discrepancy. This clean dynamical separation defines a new, physically motivated classification scheme for stellar systems, unifying galaxies and clusters under one framework. The observed division emerges naturally within the MOND framework and provides a useful diagnostic for examining how different gravitational paradigms account for the origin of the mass discrepancy.

We present the first spectropolarimetric time-series analysis of the Maunder Minimum analog HD 166620, using 12 nights of data from CFHT/SPIRou and a single epoch from CFHT/ESPaDOnS. While individual Stokes $V$ profiles exhibit no significant polarization signatures, we leverage the rotational coverage of the SPIRou dataset to compute a grand average LSD profile. Forward modeling of the cumulative Stokes $V$ signal, assuming a purely axisymmetric dipole, yields a best-fit dipole field strength of $B_{\rm dip} = 1.10^{+0.95}_{-0.90}$G ($3\sigma$). This field strength matches simulations of the solar dipole during the Maunder Minimum phase. Our results are consistent with independent constraints on the dipole field strength from an LBT/PEPSI snapshot and exclude the presence of strong non-axisymmetric fields potentially missed by this single-epoch observation. These findings provide direct empirical evidence that the transition to weakened magnetic braking involves a weakening of the large-scale magnetic field and suggest that HD 166620 represents a state comparable to the Sun near the peak activity of a grand minimum.

A recent development in the study of magnetic cataclysmic variable stars (mCVs) has been the identification of asynchronously spinning mCVs with orbital periods <2 h that have significantly higher white dwarf spin-to-orbital period ratios than their longer-period counterparts. We report the discovery of two additional mCVs in this class. The first, Gaia21akb, is a candidate asynchronous polar at the period minimum. While TESS photometry cannot, in isolation, lead to a conclusive identification of the orbital period, the probable orbital period of 1.29 h would be the second-shortest of any known polar and would result in a spin-to-orbit ratio of 0.9879. The second system in our study, ZTF18aazmehw, is an eclipsing mCV with a 1.50 h orbital period and a spin-to-orbit ratio of 0.867. Contrary to expectations for an asynchronous polar, ZTF18aazmehw does not show discernible evidence of pole switching and might possess a disk-like structure. The increasing number of short-period asynchronous mCVs with large spin-to-orbit ratios lends credence to theoretical predictions that asynchronously rotating mCVs with sufficiently strong white dwarf magnetic fields can achieve synchronization when their orbital separations have shrunk sufficiently.

The evolution of starspots of the giant primaries of RS CVn systems is typically detected indirectly with photometric and spectroscopic monitoring. These observations suggest slowly-evolving stellar surfaces and can constrain differential rotation as starspots move with respect to one another. However, starspot latitudes are difficult to constrain without resolved images of the stellar surfaces from which the unambiguous locations of starspots are determined. We imaged the active RS CVn primary $\zeta$ And with the 330-m-baseline Center for High Angular Resolution Astronomy Array for three epochs over approximately six rotations of the star. The resultant images show a more complicated picture of stellar activity than expected from the contemporaneous photometry and earlier Doppler images. The spot structures change on the timescale of rotation, making differential rotation difficult to study. Our observations show changes in the polar spot, growing over time. We do not detect the secondary star in the interferometric data, though the observations are sensitive to the predicted 0.75 $M_\odot$ main-sequence star, and we suggest the companion may be a white dwarf.

Scientists across all disciplines share a common challenge: the divide between their theoretical knowledge and the specialized skills and time needed to build interactive tools to communicate this expertise. While large language models (LLMs) offer unparalleled acceleration in code generation, they frequently prioritize functional syntax over scientific accuracy, risking visually convincing but scientifically invalid results. This work advocates the Scientist-AI-Loop (SAIL), a framework designed to harness this speed without compromising rigor. By separating domain logic from code syntax, SAIL enables researchers to maintain strict oversight of scientific concepts and constraints while delegating code implementation to AI. We illustrate this approach through two open-source, browser-based astrophysics tools: an interactive gravitational lensing visualization and a large-scale structure formation sandbox, both publicly available. Our methodology condensed development to mere days while maintaining scientific integrity. We specifically address failure modes where AI-generated code neglects phenomenological boundaries or scientific validity. While cautioning that research-grade code requires stringent protocols, we demonstrate through two examples that SAIL provides an effective code generation workflow for outreach, teaching, professional presentations, and early-stage research prototyping. This framework contributes to a foundation for the further development of AI-assisted scientific software.

In this Letter, we establish a robust, physically motivated classification method using a Support Vector Machine (SVM) trained on a "gold-standard" sample of 24 GRBs with spectroscopically confirmed progenitors (associated SNe or KNe). By isolating the prompt main spike to excise contamination from extended emission, we derive a quantitative classification index, I_SVM = 5.01 log_10 E_p,i - 1.25 log_10 E_iso - 0.34 log_10 T_90,z - 12.90 (units: keV, 10^52 erg, s). Events with I_SVM > 0 are classified as mergers. Analysis of the standardized classification weights reveals that the discriminative power of E_p,i is approximately 5 times that of T_90,z, while E_iso contributes a weight comparable to E_p,i. This quantitatively demonstrates that spectral hardness and energetics, rather than duration, are the primary physical signatures distinguishing mergers from collapsars. The derived boundary implies a stringent hardness ceiling for collapsars, while mergers are identified as outliers with excessive hardness relative to their energy budget. The classifier successfully identifies the nature of historic test cases, including the ultra-long GRB 111209A (collapsar) and the short GRB 050709 (merger), independent of instrumental eras. This tool paves the way for cleaning archival and future high-redshift GRB samples for precision cosmology.

Halo star clusters serve as vital tracers for the formation and evolution of the Andromeda galaxy. In this work, we present physical parameters for 29 M31 halo star clusters, derived from a combination of spectroscopic and photometric data. Low-resolution spectra were acquired using the BFOSC spectrograph on the NAOC Xinglong 2.16-m telescope. For the photometric analysis, we utilized uSC and vSAGE bands from the SAGE survey, complemented by archival data from GALEX(NUV, FUV), PAN-STARRS(grizy) and the 2MASS(JHK). Ages and metallicities were determined via ULySS (Vazdekis et al. and pegase-hr) SSP model and the Bruzual & Charlot (2003) (BC03) stellar population synthesis models. The derived parameters show good agreement with literature values. Notably, for three of these clusters, this study represents the first combined photometric and spectroscopic analysis.

The ~ 200,000 stars observed by the Kepler mission have provided unprecedented constraints across astrophysics. With the advent of modern spectroscopic and photometric surveys, new limits in stellar characterizations are within reach. In this work, we report a compilation of atmospheric parameters (Teff, logg, and [M/H]) for the Kepler stars by crossmatching with several spectroscopic and spectro-photometric surveys. We use these to calculate bolometric corrections, which combined with color-magnitude diagram (CMD) information from Gaia yield self-consistent luminosities on a survey-by-survey basis. These properties will aid in future explorations of Kepler data towards new astrophysical insights. We make our catalog publicly available online in Zenodo (doi:https://doi.org/10.5281/zenodo.18620911).

We present the X-ray analysis of coronal properties in a statistically representative sample of 23 mostly radio-quiet AGN from the SUBWAYS campaign (SUpermassive Black holes Winds in XrAYs), focusing on quasars at redshifts $0.1 < z < 0.4 $ and bolometric luminosities $2 \times 10^{44} <L_{bol}(erg/s) < 2 \times 10^{46}$. The main aim of this work is to investigate the properties of the hot corona through the study of the hard X-ray band emission, including a proper treatment of the soft X-ray band. High-quality X-ray spectra from XMM-Newton, complemented by NuSTAR data extending up to 30-40 keV in the rest frame, are available for this sample. The soft X-ray band (0.3-2 keV) spectrum is best fitted by a warm corona model with a median temperature of 0.40 keV, and an optical depth in the range $\tau$=20 - 40, consistent with previous results on lower luminosity sources. The hard X-ray band is well described using a hot corona model, with a median high-energy cut-off of 87 keV, at the lower end of the distribution of typical values found in Seyfert galaxies (100 - 200 keV). The derived median value of the optical depth ($\tau$ = 1 - 5) suggests the presence of a moderately optically thick corona. Combining the SUBWAYS results with literature samples at low and high redshift, we assemble the largest sample to date of AGN with E$_{cut}$ and accretion parameter measurements, finding a significant anticorrelation of E$_{cut}$ with both $\lambda_{Edd}$ and $L_{bol}$ with the median E$_{cut}$ decreasing from 250 - 300 keV at low accretion rates and luminosities to 90 - 100 keV at high accretion rates and luminosities - consistent with enhanced coronal cooling, possibly driven by pair-production. These results favor cooler, optically thicker coronae in luminous AGN compared to those in lower-luminosity Seyfert galaxies.

The icy material within protoplanetary disks plays a central role in planet formation, yet remains poorly characterized by observations. We present 1.6-28$\mu$m spectra of five disks obtained as part of the JWST Edge-on Disk Ice (JEDIce) program, representing the largest survey of disk ices to date. The major ice species H$_2$O, CO$_2$, and CO are detected towards all disks, and exhibit a wide range of absolute optical depths and optical depth ratios across the sample. This is suggestive of a range of ice abundances and compositions, but quantitative constraints will require radiative transfer modeling. All disks exhibit ice features across the entire spatial region where the IR continuum is detected; vertically elevated ice grains therefore seem to be ubiquitous in disks. The CO ice is consistently dominated by apolar CO:CO$_2$ mixtures, implying that the disk ice compositions are neither completely reset nor pristinely inherited from the protostellar stage. The presence of these mixtures also suggests that entrapment may be important in shaping the spatial distribution of CO within the disks. Small molecules commonly seen in protostellar ices (CH$_4$, CH$_3$OH, NH$_3$) are generally not detected in our sample, though tracers of ammonium salts (OCN$^-$ and the 6.85 $\mu$m band) are common, potentially reflecting an evolution towards comet-like ice compositions. The spectra also contain a wealth of information about the micron-sized dust, atomic and molecular gas, and PAH content, which together with the ice constraints will provide a comprehensive picture of the chemical, physical, and dynamical state of these systems.

In this work, we introduce TUNeS (Temporal UNet emulator for Structure formation), a neural network framework for accelerating N-body simulations by predicting the nonlinear evolution of the matter density field from an initial particle distribution. TUNeS employs a two-stage modeling strategy, combining particle-based inference with a density-field refinement on a regular grid, enabling accurate reconstruction of both large- and small-scale structures. The model is designed to operate across redshift, taking particle snapshots at arbitrary input redshifts and predicting density fields at arbitrary target redshifts. In this work, we evaluate its performance using simulations initialized at $z=100$, with predictions generated at multiple lower redshifts. Trained on only eight N-body simulations, TUNeS reproduces reference results with good agreement in both Gaussian and non-Gaussian statistics, including two-point correlations, one-point distributions, peak counts, and three-dimensional Minkowski functionals. In particular, at $k \simeq 1\,h\,\mathrm{Mpc}^{-1}$, the power spectrum error remains at the few-percent level. End-to-end inference from $256^3$ particles to a $256^3$ density grid can be completed in $\sim25\,\mathrm{second}$ on a single GPU. Thanks to its architectural design, the model naturally scales to larger particle numbers and larger volumes through particle batching and window-based refinement.

The emergence of habitable conditions on the early Earth and on rocky exoplanets requires persistent energy sources that can drive both prebiotic chemistry and climate warming under magnetically active young G to M stars. To quantify the contribution of stellar energetic particle (StEP) events associated with superflares to the atmospheric chemistry of young planets with primitive atmospheres, we carried out a suite of laboratory proton irradiation experiments on mildly reduced gas mixtures. We present first proton irradiation experiments of $N_{\mathrm{2}}O$/$CO_{\mathrm{2}}$ rich gas mixtures that yield abundant nitrous oxide ($N_{\mathrm{2}}O$) at mixing ratios up to 1000 ppmv, together with amino acid precursors including glycine, corresponding to global production rates of order $2\times10^{10}$ kg/yr on the early Earth. Our photochemical modeling of StEP driven proton irradiation reproduces the experimentally inferred $N_{\mathrm{2}}O$ production rates and provides self-consistent atmospheric $N_{\mathrm{2}}O$ profiles. We then use these profiles of $N_{\mathrm{2}}O$ as input to a 3D global climate model to evaluate the radiative and climatic impact of StEP generated $N_{\mathrm{2}}O$ in primitive atmospheres representative of the early Earth and young rocky exoplanets. Our results show that frequent StEP events can help alleviate the faint young Sun paradox on the early Earth and can maintain temperate surface conditions on young rocky exoplanets beyond the outer edges of habitable zone, while simultaneously enhancing the buildup of prebiotic molecules. Together, these processes may constitute a robust pathway toward early planetary habitability.

The TeV neutrino signal observed by IceCube from the active galactic nucleus (AGN) NGC 1068 can probe its innermost coronal regions. If these neutrinos originate from hadrons accelerated within a magnetized turbulent corona, their intensity and spectrum depend on the turbulent magnetic field strength and turbulence coherence scale. The gamma rays accompanying neutrino production are absorbed in this optically thick environment, in a way that depends sensitively on the size of the corona. By a joint fit of the IceCube and Fermi-LAT observations, we translate the multimessenger signal from NGC 1068 and the tentative signal from NGC 7469 into quantitative constraints on coronal properties. NGC 1068, with a significant TeV neutrino excess, favors a compact, strongly magnetized corona with a large turbulence coherence length relative to the coronal size. NGC 7469, with two $\sim 100$ TeV neutrino events, points instead to a somewhat larger corona with much smaller coherence length and high magnetization, but a very small fraction of energy in non-thermal protons. We obtain the diffuse flux from a population of Seyfert galaxies identical to either NGC 1068 or NGC 7469. Finally, we consider a third scenario, motivated by the spectral break observed in the diffuse neutrino flux at tens of TeV, with coronal properties intermediate between the two point-source-inspired models. To enable detailed comparisons with the IceCube and electromagnetic observations, we release our model predictions in a GitHub repository.

Testing the $\Lambda$CDM model requires cosmological probes spanning the wide redshift interval between Type Ia Supernovae (SNe Ia, $z\lesssim2.9$) and the Cosmic Microwave Background (CMB, $z\approx1100$). Gamma-Ray Bursts (GRBs), observed up to redshift $z=9.2$, offer the opportunity to explore this regime. Here, we investigate how many GRBs are needed to become a useful cosmological probe capable of independently testing deviations from $\Lambda$CDM suggested by the recent DESI BAO observations. We develop forecasts based on the two-dimensional X-ray and optical Dainotti relations, between the luminosity at the end of the plateau phase and its rest-frame duration. Using simulated GRB samples constructed from the observed population, we evaluate the constraining power of GRBs on cosmological parameters within the $w$CDM and $w_0w_a$CDM models, both independently and in combination with CMB observations. Our results show that GRB samples containing several tens to hundreds of well-characterized plateau can already approach the precision currently achieved by CMB measurements on the Dark Energy (DE) equation-of-state parameter $w$. Particularly, a sample of $\sim66$ optical GRBs can reach a precision $\sigma_w \approx 0.47$, comparable to that obtained from Planck within the $w$CDM framework. Such sample sizes are already attainable through Machine Learning techniques that double the number of GRBs using inferred redshifts. These forecasts indicate that future GRB observations, when combined with next-generation transient missions and improved statistical techniques, will provide an independent high-redshift probe of cosmic expansion and will play an important role in testing the robustness of potential Dynamical DE signals suggested by other cosmological datasets.

A particularly intriguing subclass of rocky exoplanets are the ultra-short period (USP) worlds that orbit their host stars in less than a day. These planets are particularly rare around M dwarf stars, with so far only ten that have a constrained mass and radius. We present the validation and characterization of the ultra-short period (0.3-days), Earth-sized planet TOI-4552b orbiting a nearby (27.26-pc away) M4.5V dwarf. Complementing the TESS photometry, ground-based transit observations from LCO, ExTrA and SPECULOOS validated the planetary radius and cleared the field of any contaminants. Speckle imaging with Zorro (Gemini-S) rules out false positive scenarios caused by eclipsing binary sources. Spectroscopic observations with NIRPS and HARPS were used to obtain stellar abundances, constrain the planetary mass, and, in conjunction with the transit observations, estimate the orbital parameters. TOI-4552 is a quiet star exhibiting no short-term stellar variations seen in photometric or radial velocity data that can be associated to stellar rotation. TOI-4552b ($M_p=1.83\pm0.47\,M_e$, $R_p=1.11\pm0.04\,R_e$) lies between the Earth-like and iron-rich composition tracks on the Mass-Radius diagram. The EXOPIE interior structure model, without constraints from refractory abundance ratio, yields a core mass fraction (CMF) of 0.54 and a bulk density of 7.74g/cm$^3$. Since the CMF spans a wide range due to the large uncertainty on the mass, the definitive interior composition cannot be determined with the current dataset. TOI-4552b hints as being marginally more iron-rich compared to the Earth but confirmation of its status requires additional, precise radial velocity measurements. Combined with its high emission spectroscopic metric (ESM=19.5), negligible stellar activity and short orbital period, TOI-4552b emerges as a compelling target for atmospheric and surface composition studies with JWST.

We study the non-Gaussianity of the large-scale clustering of high-redshift halos, seeking to assess which terms of standard bias expansions are needed to understand these highly biased populations. We find that the clustering can be well modeled with only linear and quadratic bias parameters while assuming a Gaussian underlying matter field. Our analysis focuses on AbacusSummit halos at redshift $z=8$. We work with halos of mass at least $1\times10^{11}h^{-1}M_\odot$ in boxes of side length $2h^{-1}$Gpc. Measurements of bias coefficients are made by fitting bias expansions to the halo power spectrum and bispectrum. Tidal bias is not detected with only a ~$0.1\sigma$ deviation from $0$, but we see a $17\sigma$ level detection for a bias term of the form $\delta^2$. A bias term of the form $\delta^3$ is weakly detected at the $1.3\sigma$ level. Nonlinear matter is also detected at a $1.3\sigma$ level. To test how bias evolves, we run one test at $z=5$. We use a mass threshold for halos that gives the same variance in the halo field as our $z=8$ sample. Bias is smaller at $z=5$ and a tidal bias is detected at the $1\sigma$ level. Bias coefficients at $z=5$ match a linear evolution of the $z=8$ bias coefficients to within $10\%$.

Nuclear Star Clusters (NSCs) are dense stellar systems located at the centers of galaxies. Employing Enzo-Abyss, which integrates hydrodynamics with a direct N-body solver, we introduce a simulation capable of resolving the evolution of NSCs within a live galaxy. This includes live dark matter, gaseous dynamics, star formation and feedback, collisional dynamics for star clusters. The evolution of NSCs is typically shaped by two main processes: mergers of star clusters and in-situ star formation. Our simulation enables investigation of the contributions of these mechanisms to the growth of NSCs. This work focuses on the impact of stellar physics and gas content on the growth of NSCs within a dwarf galaxy. To this end, we carry out four simulations, a fiducial simulation, one without supernova feedback, one with low star formation efficiency, and one with higher galactic gas content. This study shows a likelihood that both mergers and in-situ star formation contribute to NSC evolution comparably. In addition, mergers result in disruption of dense gas clumps within star clusters, indicating that in-situ star formation is suppressed when mergers occur. However, the limitations -- such as the lack of individual star physics and limited spatial/particle mass resolution -- hinder drawing a definite conclusion. Nevertheless, with further development, our simulations will serve as a cornerstone that untangles the complex interplay between mergers and in-situ star formation in shaping the structure and mass of NSCs, thereby providing insights into their formation and evolution.

High-energy astrophysical events, particularly Gamma Ray Bursts (GRBs), have been proposed as significant contributors to mass extinction events on Earth-like planets in most of the galaxy, internal to our radius in it. This paper examines the extent to which GRBs may reset the evolutionary progress of complex life through repeated extinction-level disruptions. While resilient extremophiles may survive even the most intense GRBs, more complex surface-dwelling organisms are vulnerable to indirect atmospheric effects, primarily UV exposure following ozone depletion. By identifying evolutionary milestones and estimating how frequently GRBs would need to occur to prevent recovery between such milestones, this work proposes that GRBs could act as evolutionary filters, limiting the emergence of advanced life, but only much closer to the galactic center. We consider the implications for searches of various biosignatures versus technosignatures.

Changing-look active galactic nuclei (CL-AGN) exhibit spectroscopic and photometric changes on timescales of months to years, making them powerful laboratories for studying accretion variability onto supermassive black holes. Motivated by the growing relevance of large spectro-photometric time-domain surveys, especially the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), we compiled a master catalog of known CL-AGN from the literature and evaluated its spatial overlap with the Rubin survey footprint. Using a geometric cross-match based on sky coordinates, we identify 79 sources located in high-cadence regions of the main survey footprint (Wide-Fast-Deep, or WFD), including 5 particularly favorable targets lying within the Deep Drilling Fields (DDFs) of COSMOS and XMM-LSS. These sources represent especially promising candidates for future variability studies in the Rubin era. This Research Note presents a first proof of concept for connecting known CL-AGN with Rubin observing fields, while the full catalog and a more comprehensive analysis will be presented in a forthcoming paper.

We present Active Learning for Accelerated Bayesian Inference (\texttt{alabi}): an open-source Python package for performing Bayesian inference with computationally expensive models. Given a forward model and observational data to construct a likelihood and priors, \texttt{alabi}\ uses a Gaussian Process (GP) surrogate model trained to predict posterior probability as a function of input parameters, and employs active learning to iteratively improve GP predictive performance in high-likelihood regions where the GP is most uncertain. \texttt{alabi}\ provides a uniform interface for using Markov chain Monte Carlo (MCMC) with different packages, including the affine-invariant sampler \texttt{emcee}, and nested samplers \texttt{dynesty}, \texttt{multinest}, and \texttt{ultranest}. This approach facilitates accurate estimation of the desired posterior distribution, while reducing the number of computationally expensive model evaluations required by factors of thousands. We demonstrate the performance of \texttt{alabi}\ on a variety of test cases, including where inference is challenging due to complex posterior structure or high dimensionality. We show that \texttt{alabi}\ offers a substantial improvement for likelihood functions with evaluation times $\gtrsim 1$\,s, speeding up MCMC computations by a factor of $10-1000\times$ when tested on problems with up to 64 dimensions.

We identify different dust features in our compilation of infrared spectra for 267 planetary nebulae (PNe) from the Spitzer, ISO, and IRAS telescopes. We classify 206 objects according to their dust type: mixed dust (MD), oxygen-rich dust (ORD), carbon-rich dust (CRD), PNe with only polycyclic aromatic hydrocarbons (PAHs) in their spectra (oPAH), and featureless (F) PNe. We study statistically the distributions of surface brightness and diameter of PNe with different types of dust as well as their distributions in the this http URL find that both MD and ORD PNe are closer to the Galactic centre than CRD and oPAH PNe, and that the Galactic distributions of each pair of groups are statistically compatible, suggesting that they have similar progenitors. Since oPAH PNe have, on average, larger diameters and lower surface brightness than CRD PNe, we suggest that oPAH PNe are evolved CRD PNe. On the other hand, F PNe have the lowest surface brightness and the largest diameters, suggesting they could contain evolved PNe from any initial type of dust. Among the PNe with silicates, we find that only a few ORD PNe have just amorphous silicates in their spectra, and their distributions of Galactocentric distances and Galactic heights suggest that they had low-mass progenitors. MD PNe with both amorphous and crystalline silicates have the largest surface brightness and the smallest diameters and might be the earliest stages of PNe with the most massive and metal-rich progenitors.

We investigate the evolutionary pathways of green valley (GV) galaxies drawn from the SDSS-IV/MaNGA survey. The GV sample is divided into fast- and slow-rotating galaxies based on stellar spin, and their stellar and gas-phase metallicities are compared. Fast-rotating galaxies exhibit systematically higher metallicities than slow-rotating galaxies in both gas and stars. However, the gas-phase difference is significant only at low stellar masses, while the stellar metallicity offset persists across the full mass range. Using a simple yet physically motivated chemical evolution model, optimised to jointly fit gas-phase metallicities and integrated stellar spectra, we reconstruct the star formation and chemical enrichment histories of individual galaxies and constrain gas inflow and outflow parameters. At low stellar masses, fast- and slow-rotating galaxies show similar gas-infall and star formation timescales, but the the slower population experienced stronger outflows which reduce their chemical content in both gas and stars. At high masses, the combination of reduced pristine gas inflow and more efficient gas removal in slow-rotating galaxies produce gas-phase metallicities comparable to fast-rotating galaxies but systematically lower stellar metallicities. These differences suggest distinct evolutionary pathways for GV galaxies. Slow-rotating galaxies likely experienced more mergers, usually associated with strong gas removal processes, leading to their systematically lower metallicities. At low masses, stronger supernova-driven outflows reduce their chemical content while leaving star-formation timescales similar to fast-rotating galaxies. At high masses, merger-triggered AGN feedback may rapidly deplete and suppress gas infall, producing the shorter star-formation timescales seen in slow-rotating galaxies. Alternative environmental and assembly-driven scenarios are also discussed.

We present the first atmospheric observation and characterisation of the aligned, 1468 K hot Jupiter, NGTS-2b, with one JWST NIRSpec/G395H transit. These observations complete the GO 3838 observing campaign of the BOWIE-ALIGN program, which aims to investigate the link between hot Jupiter atmospheric composition and formation history through the atmospheric analysis of planets orbiting F stars that are aligned and misaligned with the host stellar spin axis. The 2.84-5.18 micron spectrum shows weak absorption features attributed to H$_2$O and CO$_2$ absorption, which our free chemistry retrievals fit with posteriors that converge on high mean molecular weight solutions attained through significant H$_2$O mixing ratios. By comparing our results to interior modelling, we show that some of these solutions exceed the 43.5x solar upper limit we obtained from our interior structure models. Such solutions are likely due to cloud-metallicity degeneracies and insufficient wavelength coverage to resolve them. We show that, in the case of our observations, the likelihood distribution of H$_2$O abundances is flat and uninformative, such that our retrievals are biased by the prior. Additionally, our statistically favoured atmospheric solution contains absorption from SO. The chemical abundances retrieved with this model are likely not astrophysically feasible and we demonstrate that the presence of SO is driven by only two data points. Our equilibrium chemistry retrievals hint at a subsolar C/O ratio and supersolar metallicity; however, we find wide posterior distributions that extend to solar values.

The release of heat by a planetary embryo modifies the local density perturbations, forming thermal lobes in its vicinity, and thereby altering the torque exerted by the disk on the embryo. In laminar disks, these thermal torques can dominate the disk-embryo interaction, rendering the classical Lindblad and corotation torques largely subdominant. The aim of this work is to investigate how turbulence driven by the MRI instability affects the thermal lobes formed around a planetary embryo, and to analyze the resulting torque acting on the embryo. We evaluate the thermal torques exerted on a planetary embryo of mass $M_p=0.33M_{Mars}$ and on a planetary core with mass $M_{p}=1M_{\oplus}$, each embedded in a turbulent gaseous protoplanetary disk, by means of high-resolution 3D magnetohydrodynamics simulations that include thermal diffusion and an initially toroidal magnetic field. The magnetic field strength is characterized by the $\beta$-plasma parameter with $\beta\in\{50,1000\}$. We consider two values for the luminosity of the planetary embryo: $L=0$ (cold thermal lobes) and $L=L_c$ (hot thermal lobes), where $L_c$ represents the critical luminosity. We find that, even in the presence of a weak magnetic field and irrespective of the luminosity, for both planetary masses, the development of turbulence in the disk (which takes between 1.5 to 3 orbital periods) completely disrupts the thermal lobes. As a result, the torque acting on both the planetary embryo and the Earth-mass core displays a strongly oscillatory behavior. This suggests that planets with masses in the range $0.03M_{\oplus}\lesssim M_{p}\lesssim 1M_{\oplus}$ experience stochastic migration, as expected in turbulent disks. Thermal torques become inefficient in turbulent regions of protoplanetary disks, such as outside the dead zone, in both the inner and outer disk regions where the magnetorotational instability operates.

The information recoverable from galaxy spectra depends fundamentally on spectral resolution, yet assembling large samples at high resolution remains observationally expensive. We present a deep-learning framework for spectral super-resolution that enhances low-resolution galaxy spectra by a factor of $\sim$10 in resolving power ($R\sim100$ to $R\sim1000$). The model is trained on 1,187 paired JWST/NIRSpec observations from the JADES program, where low-resolution prism spectra are matched with medium-resolution grating spectra (G140M, G235M, G395M) combined into a unified reference covering 1-5 $\mu$m. Our three-stage architecture performs an initial super-resolution, infers the redshift from the coarse reconstruction, and then applies a physics-informed residual refinement that uses attention across emission-line tokens to learn inter-line relationships and predict parametric line profiles, alongside a convolutional branch for continuum corrections. Evaluated on a 20% held-out sample, the model achieves noise-limited residuals over most of the spectral range and systematically improves the signal-to-noise ratio of key diagnostic lines including [OII], H$\beta$, [OIII], and H$\alpha$, often by factors of several. The super-resolved spectra successfully deblend features that are entirely unresolved at prism resolution, such as the [OIII] $\lambda\lambda4959,5007$ doublet and H$\beta$. As a proof of concept using JWST data, this approach is readily extensible to the low-resolution grism spectroscopy that will be delivered by Euclid and the Roman Space Telescope, potentially enabling population-level diagnostics across millions of galaxy spectra that would otherwise be inaccessible at grism resolution.

We present the weak gravitational lensing dataset from the Super-pressure Balloon-Borne Imaging Telescope (SuperBIT), which imaged 30 galaxy clusters during its 45 night flight in April to May 2023. SuperBIT is a first-of-its-kind balloon-borne imaging telescope that achieved near diffraction-limited observations in near-space conditions above 98% of the Earth's atmosphere. We use the metacalibration algorithm to obtain calibrated galaxy shapes for our target clusters and several calibration fields, enabling unbiased reconstruction of the weak-lensing signal. We employ several diagnostics throughout the pipeline, including assessments of point-spread function (PSF) modeling residuals and their impact on weak-lensing measurements, as well as tests for correlations between galaxy shapes and measured galaxy and PSF properties. To assess the multiplicative shear bias of the pipeline, we analyze a parallel set of simulated images that incorporate the real observing conditions from the flight, including measured SuperBIT PSFs, observed sky backgrounds, and detector noise, yielding a bias of $(1.1 \pm 7.8)$~per~cent.

The galaxy-halo connection traces processes by which galaxies form and evolve. The halo occupation distribution (HOD) describes the relationship between galaxies and their host dark matter haloes. Measurements of the galaxy two-point correlation function (2PCF) allow us to extract information about the HODs of observed galaxy samples. Several parametric HOD models have been proposed in the literature, but the choice of parameterization restricts the space of possible HODs. To resolve this issue, we introduce a non-parametric HOD fitting method in which we train an emulator to learn the mappings among the galaxy 2PCF, physical properties used to select galaxy samples, and the HOD, all obtained from simulated past lightcones constructed with the Santa Cruz semi-analytic models. Implementing this emulator within a likelihood analysis framework, we derive constraints on the HOD of a galaxy sample when provided with a measurement of its 2PCF. Using the emulator to accelerate likelihood evaluations, we test the non-parametric HOD approach on a set of 2PCFs for mock galaxy samples drawn from the TNG100-1 simulation and selected above threshold values of stellar mass and star formation rate. Our framework is able to recover TNG100-1 HODs within 0.2 dex. We use the TNG100-1 mocks to tune the reported uncertainties to estimate those expected in the analysis of observations. Comparing to parametric HOD modeling routines applied to the same mock galaxy samples, our approach consistently infers the HOD with comparable or greater precision and accuracy.

TOI-1232 is a G-dwarf star with a mass of $1.06_{-0.06}^{+0.07} M_\odot$, a radius of $1.07\pm 0.05 R_\odot$, and slightly higher metallicity than solar of Fe/H = $0.18 \pm 0.05$. The star hosts a transiting warm Jovian-mass planet, TOI-1232 b, with an orbital period of $P_{b} = 14.256_{-0.001}^{+0.001}$ days, identified with data from multiple sectors of the $\textit{TESS}$ space telescope. The $\textit{TESS}$ light curve of TOI-1232 is complex, as it is contaminated by a background eclipsing binary with a period of $1.37$ days. The TOI-1232 b was firmly confirmed by ground-based transit follow-up campaigns from Las Cumbres, Hazelwood, Brierfield, and ASTEP this http URL, the $\textit{TESS}$ transits of TOI-1232 b exhibit strong transit timing variations (TTVs) with a super-period of $235.5 \pm 0.7$ days and a semi-amplitude of 27 minutes. Radial velocity (RV) follow-up with the FEROS spectrograph confirms the planetary nature of the transiting candidate, while a self-consistent $N$-body analysis of RVs and TTVs pinpoints the presence of a second outer Saturn-mass companion, TOI-1232 c with a period of $P_{c} = 30.356_{-0.012}^{+0.010}$ days. The TOI-1232 warm-giant system is particularly important due to the evidence of two massive planets that reside near the 2:1 commensurability but are not locked in a mean motion resonance (MMR). Thanks to $\textit{TESS}$, we have revealed a handful of these rare systems. Hence, TOI-1232 is an important addition to understanding the formation and dynamical evolution of such compact, massive, warm giant planets.

Weak gravitational lensing is a widely used probe in cosmological analysis. It allows astrophysists to understand the content and evolution of the Universe. We are entering an era where we are not limited by the data volume but by systematic uncertainties. It is in this context that we present here a simple python-based software package to help in the computation of E-/B-mode decomposition, which can be use for systematic checks or science analysis. As we demonstrate, our implementation has both the high precision and speed required to perform this kind of analysis while avoiding a scenario wherein either numerical precision or computational time is a significant limiting factor.

More than 200 moons exist in our Solar System, yet no exomoon has been confirmed to date. While the innermost two planets of the Solar System lack natural satellites and most studies favour the existence of exomoons around long-period planets, some theoretical studies that take tidal dissipation, orbital decay, and migration processes into account suggest that exomoons may survive around short-period exoplanets. We investigated the impact of exomoons on planetary thermal phase curves and assessed their detectability within a theoretical framework. We simulated the thermal phase curves of exomoon-exoplanet systems, including mutual transits and occultations, and explored their dependence on planetary orbital periods across a wide range of systems. Close-in airless exomoons maintain large day-night temperature contrasts, amplifying the thermal phase-curve signal of the system. When the exomoon transits or is occulted by the exoplanet, the transit depth varies with the planetary phase, and the occultation depth varies with the exomoon's phase. The maximum occultation depth can reach $\sim$ 20 ppm for long-period systems. For short-period planets, the signal can reach up to $\sim$100 ppm, although such configurations may not be dynamically stable over long timescales. If exomoons are not accounted for, the planetary temperature distribution retrieved from observed thermal phase curves may overestimate the planetary day-night temperature contrast and underestimate the planetary horizontal heat transport. In principle, the periodic exomoon-exoplanet mutual occultation signal could be extracted using methods such as box-fitting least squares, providing a framework for future observational studies and instrument planning.

The physical origin of fast radio bursts (FRBs) remains an unsolved mystery in astrophysics, with the magnetar central engine model as the leading framework. Systematically searching for physical associations between FRBs and the energetic astrophysical transients (ATs) that form magnetars provides a critical test of this scenario, and key clues to FRB progenitors. We perform a systematic search for FRB-AT associations using a sample of 3765 unique FRBs, combining the second CHIME/FRB catalog with 124 additional localized FRBs with measured redshifts. We develop a 3D Bayesian inference framework that jointly incorporates angular separation, positional uncertainty, and redshift constraints to quantify the association probability of candidate pairs. Through spatial cross-matching, we identify 14 FRB-optical transient and 15 FRB-gamma-ray burst (GRB) candidate pairs. Our framework recovers the previously reported high-significance association between FRB 20180916B and AT 2020hur, with an association probability of 0.9998. For the proposed candidate FRB 20190309A and short GRB 060502B, our analysis yields an association probability of 0.83, which is insufficient to claim statistically significant association. No new statistically significant FRB-AT associations are found for all remaining candidates. Our work demonstrates that small angular separation alone is insufficient to confirm FRB-AT associations, and high-precision FRB localization is essential for definitive identification.

We analyze a clean sample of 1,118 late-type, face-on galaxies without AGN contamination from the MaNGA survey. Their photometric structures are quantified via two-component (bulge+disk) decompositions on deep $g$-band images from the DESI Legacy Survey. Using a disk central surface brightness of $\mu_{\rm 0,d,cor}$(g) = 22 $\pm$ 0.3 mag arcsec$^{-2}$ (corrected for inclination and cosmic dimming) as the classification threshold, we identify 159 low surface brightness (LSB) galaxies, 388 LSB candidates, and 571 high surface brightness (HSB) galaxies. LSB galaxies are predominantly low-mass ($M_\ast < 3 \times 10^{10}$ M$_\odot$), exhibiting 29\% larger effective radii, 15\% lower star formation rates (SFRs), and 12\% reduced gas-phase metallicities than HSB counterparts at comparable masses. These differences cause systematic offsets from standard scaling relations. Despite comparable gas content, LSB galaxies host older stellar populations, longer gas depletion times, and less efficient star formation. Spatially resolved analyses further reveal that LSB galaxies display centrally suppressed $\Sigma_{\rm SFR}$, flatter SFR gradients, and rising specific SFR profiles toward their outskirts. Together with steeper negative metallicity gradients, these trends suggest ongoing gas accretion fueling outer-disk star formation. Consistently, the outer regions of LSB galaxies exhibit stronger H$\delta_A$ absorption and lower D$_n$4000 indices, indicating fading A-star populations. Moreover, LSB galaxies show lower $\Sigma_{\ast}$ across all $R/R_e$ and more centrally depleted stellar mass profiles on an absolute radial scale, compared with HSB and large-size star-forming galaxies. Collectively, LSB galaxies represent a distinct population with slow evolution, inefficient star formation, and continued susceptibility to late-time gas accretion and peripheral star formation.

We present high spatial resolution ($\lesssim$1.0''), multi-wavelength observations of UGC 2369S, a nearby luminous infrared galaxy showing three distinct cores separated on kpc scales in near-infrared (NIR) imaging with significant X-ray emission. Utilizing optical/NIR adaptive optics (AO), radio, \chandra X-ray, as well as archival HST imaging, we perform a comprehensive study of AGN activity, obscuration, and host properties. As one of the clearest cases of a triple-nucleus merger at $\simeq$3 kpc separations, UGC 2369S is the first to be studied with high-resolution observations at multiple wavelength. We find that the northern core, having possibly the most massive black hole in the system ($\rm M_{BH}\simeq10^{8}\,M_{\odot}$) is consistent with a heavily obscured AGN. However, its high dust extinction ($\rm A_v>5$), hydrogen column density ($N_\mathrm{H}\gtrsim 10^{25}\,\rm cm^{-2}$) and non-detection of optical coronal lines and coronal X-ray emission leave the identification inconclusive. The other two cores show no evidence for black-hole activity and instead exhibit signatures of tidal disruption. From stellar mass surface density and stellar velocity dispersion maps, we infer that the strongly varying gravitational potential in this three-body system may have cannibalized the stellar bulge of the southwestern core, leaving a metal enriched remnant. An ongoing survey focusing on similar triple systems could help us understand how they evolve and help benchmark numerical simulations, providing insight into gravitational wave predictions and the formation of the most massive black holes.

We report our optical observations of the Crab pulsar using the Imager of MPPC-based Optical photoN counter from Yamagata (IMONY), a high-time-resolution photon-counting imager with 100 ns timing resolution, mounted on the 3.8 m Seimei telescope in Japan (f/D~6). The detector format was upgraded from a $4\times4$ to an $8\times8$ GAPD array with larger pixels ($100$ to $200~{\mu m}$), resulting in a 14".5 field of view on the Seimei telescope. We conducted nightly optical observations for one week, including two nights of simultaneous optical and radio observations with the 64 m Usuda radio telescope. Thanks to the large diameter of the Seimei telescope and the high time resolution of IMONY, we successfully detected optical Single Pulses (SPs) emitted in each rotation. Moreover, we found an optical peak timing drift of $30\pm7.9~\mathrm{\mu s}$ over three days, with a significance of $3.9{\sigma}$. The corresponding emission region size is 9.1 km, which is equivalent to 0.006 times the light cylinder radius of the Crab pulsar. We ruled out the possibility of a pulsar glitch and suggested that the optical pulsed emission region of the Crab pulsar may fluctuate due to the spatial drift and variations in the magnetospheric caustics.

We present the first statistical census of emission-line variable active galactic nuclei (EVA) at cosmic noon by combining untargeted and deep HETDEX spectroscopy with multi-epoch spectra from SDSS, DESI, and LAMOST. Anchoring all candidates to a HETDEX spectroscopic epoch and requiring AGN classification in either the HETDEX or the external epoch(s), we identify a homogeneous sample of 100 EVA at z~1.5, including 98 newly identified. Emission-line variability is selected primarily through statistically significant line-flux changes, supplemented by extensive visual inspections using contemporaneous photometric light curves. The resulting incidence fraction is $f_{\rm EVA} \approx 0.9\%$. The rest-frame intervals between spectroscopic epochs span $\sim$1--10 yr, with brightening and dimming events exhibiting statistically indistinguishable characteristic timescales ($\Delta T\sim2.2$ and $\sim2.6$ yr, respectively). A key result is the characterization of the Baldwin effect in the time domain: while many EVA follow the ensemble Baldwin effect (eBeff) between two epochs, a substantial fraction exhibit apparent anti-eBeff responses. Time-resolved spectroscopy of an individual source reveals that the intrinsic EW--luminosity relation is non-stationary, with the line-to-continuum responsivity systematically evolving from stronger to weaker across successive variability cycles; sparse two-epoch sampling of this evolving intrinsic Baldwin evolution (iBeff) naturally produces both eBeff-like and anti-eBeff behaviors. Finally, EVA show no strong preference for extreme Eddington ratios but exhibit a mild tendency toward lower $\lambda_{\rm Edd}$ values relative to matched control samples, driven primarily by sources observed in their dim states. Together, these results establish a coherent framework for interpreting emission-line variability in AGN at the peak epoch of cosmic black hole growth.

As stars evolve, they undergo significant changes in their physical properties, which can have a profound impact on the planets orbiting them. In particular, the mass lost through stellar wind may be partially accreted by orbiting planets. We present the results of 18 simulations of one-planet systems with planetary masses of 0.5, 1, 2.5, 5, 10, and 13~$\mathrm{M}_\mathrm{J}$, each at initial orbital distances of 5, 10, and 20~AU, orbiting a 2~M$_\odot$ star through its red giant branch and thermally pulsating asymptotic giant branch phases. Our results show that planets with smaller orbits and higher masses accrete more stellar wind material than their wider-orbit and lower-mass counterparts, although the total mass accreted across all simulations remains small compared to their initial planetary mass. Even for the most massive planet, 13 $\mathrm{M}_\mathrm{J}$ at 5 AU, the total mass accreted was $\sim0.56$\% of the planet's initial mass; nevertheless, we find that the accretion luminosities of the simulated planets, with the exception of one planet, exceed their expected equilibrium luminosities, suggesting that such emission could be potentially detected. This result is key for the detection of planets around AGB stars, which have no confirmed detections as of yet. We also estimated the accretion and luminosities of two detected two-planet systems over a few orbits, obtaining results consistent with the one-planet simulated systems. Additional tests without wind accretion and with stellar wind drag force showed that, while both have a negligible effect on the orbital evolution, wind accretion remains relevant for the planetary luminosity.

Using a sample of 116 wide binary systems as coeval and chemically homogeneous stellar pairs, we investigate the factors governing lithium depletion in main-sequence stars. We recover the well-established morphology of the lithium--effective temperature ($T_{\mathrm{eff}}$) relation, including the Li dip (6200--6600\,K), the Li plateau (6000--6200\,K), and a linear trend for cooler stars ($T_{\mathrm{eff}}$ $<$ 6000\,K), where lithium abundance increases by $\sim$0.15\,dex per 100\,K. We demonstrate that the apparent correlation between projected rotational velocity ($v\sin i$) and lithium abundance is secondary to the underlying $T_{\mathrm{eff}}$ dependence; $v\sin i$ is not an independent driver of lithium depletion in our sample. Notably, we identify an anomalous system within the Li dip where the primary star exhibits a $\sim$1.4\,dex lithium excess compared to its secondary companion at nearly identical $T_{\mathrm{eff}}$. We discuss two plausible origins for this anomaly: external enrichment via planetesimal accretion or planetary engulfment, and binary interactions with an unresolved tertiary companion. Our results confirm $T_{\mathrm{eff}}$ as the dominant parameter controlling lithium depletion, while highlighting that additional, non-standard processes can occasionally produce significant lithium enrichment.

Comet 12P/Pons-Brooks is a periodic comet with an orbital period of approximately 71 years. Because of the period duration, aphelion of 35.3 au, and highly inclined orbit (74.2 deg), it is classified as a Halley-type comet. In this work, we present the results of the analysis of photometric and spectral observations of comet 12P/Pons-Brooks while it was at a heliocentric distance beyond 1.1 au, before perihelion passage. Quasisynchronous observations were carried out from July 2023 to March 2024 using 0.61-m and 1.3-m telescopes at the Skalnaté Pleso Observatory, 1-m telescope Zeiss-1000 Sanglokh International Astronomical Observatory, and 0.7-m telescope AZT8 at the observation station of Taras Shevchenko National University of Kyiv. Photometric observations were conducted using B, V, and R filters of the Johnson-Cousins and Bessel photometric systems. We completed our data collection with broadband photometric and spectral observations on the 2-m HCT telescope at the Indian Astronomical Observatory. During this period, the dust coma morphology and photometric data measured in the R filter, such as the apparent, absolute magnitudes and dust activity level, indicate several repeated outbursts in brightness and dust activity occurred around July 22, September 25, October 17-23, and a significant outburst was recorded on November 1-7. We detected strong gas emission features in the cometary spectrum belonging to CN, C2, and C3 molecules. Furthermore, gas production rates of these molecules were estimated using Haser's model. A comparison of dust activity level and CN gas production rate shows that dust activity increased significantly as the comet approached the Sun, while gas activity exhibited more moderate changes.

We present a feasibility study of the Antarctic 12m Terahertz Telescope (ATT12), a next-generation facility to be constructed at New Dome Fuji in Antarctica, designed to open up the FIR and THz windows for extragalactic astronomy. While ATT12 will enable a wide range of Galactic and extragalactic science, this paper focuses on its potential for studies of dusty star-forming galaxies (DSFGs) across cosmic time. Using realistic atmospheric transmission models and the planned instrumental specifications of heterodyne spectrometers and wide-field multi-color continuum cameras, we assess the expected sensitivity and scientific capabilities. We show that spectroscopic observations will enable detections of [CII]158um from galaxies with log(LIR/Lsun)>12 out to z~7, while [OIII]88um will remain observable for HyLIRG-class systems up to z~10. Line ratios including [OIII]52/88um, [NII]122/205um, and [OIII]/[NIII] will provide unique diagnostics of electron density and O/N abundance at z~4-8. Wide-field continuum surveys with the continuum cameras (KIDS-1/2; 300-850 GHz) will reach confusion-limited depths of ~1-2 mJy over ~10,000 deg$^2$, detecting of order $10^{6}$-$10^{7}$ DSFGs with log(LIR/Lsun)>12 at z<5 and $\lesssim10^{3}$--$10^{4}$ HyLIRGs up to z~7 or higher. Higher-frequency cameras (KIDS-3/4; >850 GHz) are designed for targeted follow-up observations and to extend coverage toward the THz regime. Taken together, ATT12 will provide the first statistically representative samples of DSFGs across cosmic time and, through synergy with ALMA, JWST, and the proposed FIR Probe PRIMA, will establish a multi-wavelength framework in which ATT12 discovers large samples through wide-area surveys, ALMA provides high-resolution follow-up of gas and ISM structure, JWST probes stellar populations and metallicity in the rest-frame optical/NIR, and PRIMA delivers ultra-sensitive FIR spectroscopy.

We leveraged reliable age and distance estimates from LAMOST-DR3 and APOGEE-DR17+AstroNN combined with \gaia\ data to perform a detailed analysis of the stellar age distribution in the Milky Way's (MW) outer disc using giant stars. Selecting stars near the midplane ($|z|<0.3$ kpc) on near-circular orbits ($\lambda_c > 0.9$), we analysed these independent datasets that employed different age-estimation methods. Our stringent kinematic selection criteria effectively exclude halo stars, ensuring that the observed age trends reflect genuine disc properties rather than contamination from older halo populations. Our results reveal a 'U-shaped' stellar age profile, where a negative gradient in the inner disc transitions to a positive gradient in the outer disc region. We identify the minimum in the stellar age profile at $R_{\rm min}=11.28 \pm 0.58$ kpc and $R_{\rm min}=12.15\pm 0.62$ kpc for the APOGEE-DR17 and LAMOST-DR3 samples, respectively. Using N-body+SPH simulations, we demonstrate that $R_{\rm min}$ corresponds to the break radius in the stellar density profile ($R_{\rm br}$), marking the edge of the Galaxy's star-forming disc. This break arises from a sharp decline in the star formation rate, with the outer positive age gradient produced by the radial migration of stars born inside $R_{\rm br}$. The cessation of star formation in the outer disc might be due to several mechanisms, including the dynamical influence of the bar's outer Lindblad resonance, the onset of the Galactic warp, or thermally regulated star formation. Overall, our results support the picture that the MW has a Type II (down-bending) stellar disc with a break at $R_{\rm br} \approx 11.28-12.15$ kpc, where the combination of star-formation cut-off and radial migration produces the observed U-shaped age profile.

There is now strong evidence that the Milky Way (MW) hosts a nuclear stellar disc (NSD). However, whether the NSD is purely axisymmetric or contains a nuclear bar remains unresolved. Since approximately $50\%$ of barred galaxies with MW-like mass in the local Universe host a nuclear bar, investigating whether the MW hosts one is of interest. We conduct a systematic analysis to identify robust kinematic diagnostics capable of determining whether the MW hosts a nuclear bar. Using N-body simulations, we explore the kinematic signatures indicative of a nuclear bar. Using the phase-space coordinates longitude $(\ell)$, latitude $(b)$, proper motions ($\mu_\ell$ and $\mu_{\rm b})$ and line-of-sight velocity $(v_{\rm los})$, we test various diagnostics assuming different nuclear bar orientations. We also evaluate how sample size, dust extinction and bar amplitude influence the efficacy of the diagnostics. We identify two independent kinematic diagnostics capable of revealing a nuclear bar in the MW: (1) the vertex deviation, $l_{\rm v}$, of the ($v_{\ell}-v_{\rm los}$) velocity ellipse; and (2) The asymmetry in the $\mu_{\ell}$ vs $\ell$ distribution. While both are impacted by the sample size and extinction, the vertex deviation proves more robust, especially when combining stars from multiple observational fields. We also assess the correlation between the line-of-sight velocity and the $h_3$ Gauss-Hermite moment ("skewness") of the line-of-sight velocity but find no clear distinction between an NSD and a nuclear bar based on this metric. Our results suggest that data from the current KMOS survey may allow a marginal detection of a nuclear bar using the vertex deviation method. A companion paper provides further validation and detailed analysis of this approach. Nonetheless, future surveys will provide the high quality data necessary to fully exploit the diagnostics outlined in this study.

Context. The High Resolution Telescope of the Polarimetric and Helioseismic Imager on Solar Orbiter (SO/PHI-HRT) operates in an extreme observational environment, observing the Sun as close as $0.28$ au. The high thermal load and large illuminating field puts high demands on the instrument in terms of both imaging performance and false light control. Aims. To characterise the amount of stray light (false light) within SO/PHI-HRT, apply a correction, and re-compare the data products with the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory (SDO/HMI). Methods. We analyse solar limb profiles and a Mercury transit to quantify the amount of stray light and add a correction term when partially reconstructing the SO/PHI-HRT images. For the comparison with SDO/HMI we use data from the 2023 March Solar Orbiter inferior conjunction and compare the magnetic fields on a pixel-by-pixel basis. Results. Increased continuum intensity contrast in the quiet Sun, and darker intensity levels are found in strong magnetic features. Consequently, much stronger fields are inferred in these features. Comparing the stray light corrected data with that from the standard SDO/HMI data products results in a much closer agreement across all vector magnetic field components, particularly when the cadence and noise levels are identical. In most solar features, SO/PHI-HRT infers stronger fields than the SDO/HMI line-of-sight magnetograms. Compared to the vector magnetic field from SDO/HMI the two are very well aligned, with only slight differences in the strongest field regions (where $|\mathbf{B}|>1600$ G or $|\mathbf{B_{\text{LOS}}}|>1300$ G).

Hydrogen cyanide (HCN) is crucial for the RNA World hypothesis, forming biomolecules essential for early life. Life likely emerged around 4 billion years ago during the early Archean Eon, a period on Earth with a fainter sun, frequent impacts, and a weakly reducing atmosphere. Warm little ponds (WLPs) are hypothetical protective aqueous environments that help explain the emergence and evolution of fragile prebiotic chemistry in such a hostile environment. WLPs need to undergo cycles of evaporation and rehydration, concentrating prebiotic molecules that increase the likelihood of (de-)polymerisation and forming early RNA molecules. We use a 1-D model of atmospheric chemistry to compare atmospheric HCN delivery to WLPs with exogenous sources. Using early Archean Earth as our baseline, we examine the sensitivity of atmospheric HCN delivery to the atmospheric C/O ratio, semi-major axis, assumed stellar host type, and methane budget, exploring conditions across rocky exoplanets. We find that atmospheric HCN delivery is sensitive to these parameters but its values generally exceed that of meteoritic delivery and our baseline Archean Earth. Planetary atmospheres with higher C/O ratios within the habitable zones of G stars and those closely orbiting M-dwarfs deliver the most atmospheric HCN. We find that atmospheric HCN delivery is remarkably robust, so this molecule is likely not the rate limiting step for the emergence of prebiotic chemistry on rocky exoplanets. This finding, with important caveats, potentially increases the probability of life emerging on other worlds.

$\require{mediawiki-texvc}$Galaxy assembly was already well underway in the first 400 Myr of cosmic time, as recently revealed by JWST. However, the contribution of these early galaxies to cosmic reionisation remains uncertain. Here we present new JWST/NIRSpec observations of GS-z13-1-LA obtained as part of the OASIS and JADES programmes, whose combined deep (56 h) NIRSpec/PRISM spectrum confirms the Lyman-$\mathrm{\alpha}$ line detection and blue UV continuum at redshift $z = 13.1$ presented in a previous work. The measured Lyman-$\mathrm{\alpha}$ emission (rest-frame equivalent width of $66_{-9}^{+10}\,Å$) and steep continuum slope ($\beta_\text{UV} \approx -3$) point towards GS-z13-1-LA hosting a remarkably hot and powerful ionising source, and allow at most a modest contribution from the nebular continuum. The steep turnover of the continuum is still present, but less pronounced in the new OASIS spectrum. Combined, this implies that ionising photons may escape GS-z13-1-LA at a sufficient rate to weaken the other, still undetected UV lines, and to lead the formation of a small ionised bubble ($R_\text{ion} \approx 0.2\,\mathrm{pMpc}$). A yet larger bubble could alleviate the required ionising production efficiency of GS-z13-1-LA from $\xi_\mathrm{ion} \approx 10^{26.4}\,\mathrm{Hz\,erg^{-1}}$ down to $\approx 10^{25.9}\,\mathrm{Hz\,erg^{-1}}$, still extremely high but more readily reconcilable with stellar models. In turn, this would require a notable overdensity of galaxies with highly efficient ionising capabilities, a scenario for which tentative evidence is found in the form of 16 nearby photometric candidates and one spectroscopically confirmed source, JADES-GS-z13-0. The new OASIS observations therefore confirm the overall picture of GS-z13-1-LA as an early beacon of reionisation, providing compelling evidence for its start only 330 Myr after the Big Bang.

Using $JWST$/NIRCam data over a 0.28 deg$^{2}$ area from COSMOS-Web survey, together with $HST$/ACS data, we investigate early-type fraction of massive galaxies with $M_{star}>10^{10.5}M_{\odot}$ at $0.2<z<2.0$, and explore the formation of their early-type morphology. We measure concentration index $C$ ($=R_{80}/R_{20}$) and asymmetry index $A$, and select early-type galaxies with $C>C_{n=2.5}$ and $A_{cor}<0.2$. Here $C_{n=2.5}$ is the concentration expected for a Sersic profile with $n=2.5$ under the spatial resolution and depth of the data, and $A_{cor}$ is the asymmetry corrected for resolution effects. The fraction of early-type galaxies with $M_{star}>10^{11}M_{\odot}$ ($=10^{10.5}$-$10^{11}M_{\odot}$) decreases with increasing redshift from ~70% (~40-60%) at $z$ ~ 0.3 to ~20-25% (~15-25%) at $z$ ~ 1.8. We also examine the evolution of their $R_{20}$ and $R_{80}$, which enclose 20% and 80% of the total flux of the galaxy, respectively. The median $R_{80}$ shows strong mass dependence and significant redshift evolution, whereas the median $R_{20}$ shows little dependence on either stellar mass or redshift. In contrast, morphological differences are more pronounced in $R_{20}$ than in $R_{80}$: the median $R_{20}$ of early-type galaxies is smaller than that of late-type and irregular galaxies by 0.25-0.45 and 0.3-0.6 dex, respectively. The median SSFR of sample galaxies strongly correlates with $R_{20}$, and early-type galaxies have lower SSFRs by ~1 dex. We further find that early-type galaxies at $z>1.3$ have younger mass-weighted stellar ages of $t_{mw}<2$ Gyr than late-type and irregular ones. Their SSFRs, $t_{mw}$, and morphological properties suggest that these high-$z$ early-type galaxies experienced rapid formation of a dense stellar core through starburst, followed by quenching of star formation, and subsequently resumed star formation ~1-2 Gyr later.

Proto-stellar clusters, likely progenitors of globular clusters, are compact with typical mass $\sim 10^6\,{\rm M}_\odot$ and size $\sim 1\,{\rm pc}$, as revealed recently by JWST observations at $z\sim 10$. Sufficiently high compactness can provide a time window for early-formed stellar black holes (BHs) to accrete primordial gas. We develop a model to determine the final spin distribution of stellar BHs which grow in mass via gas accretion within compact gaseous proto-stellar clusters. The velocity shear within a BH's sphere of influence induces the formation of an accretion disk which is repeatedly disrupted by stochastic perturbations to the BH motion. We assume low initial BH spins $a_{*,{\rm ini}} = 0.01$, and restrict initial BH masses below the upper BH mass gap, $m_{\rm BH,ini} < 55\,{\rm M}_\odot$. Our analysis shows a strong BH spin-mass correlation, obtained within $\sim 10 \,{\rm Myr}$ when gas is depleted. Low-spin BHs, $a_{*} \leq 0.3$, are predominantly low-mass, $m_{\rm BH} \lesssim 25\,{\rm M}_\odot$, in contrast to high-spin black holes, $a_{*} \geq 0.7$, which are predominantly high-mass, $m_{\rm BH} \gtrsim 65\,{\rm M}_\odot$. Notably, there exist also low-spin, high-mass outliers with $\sim 1$ mass-gap BH per cluster expected to have $a_{*} \sim 0.1$. The general trend, however, expressed by the median spin as a function of final BH mass is well fit by a high-spin saturating exponential with transition mass $\sim 50\,{\rm M}_{\odot}$. For $m_{\rm BH} \geq 100\,{\rm M}_\odot$ the median spin is $\bar{a}_{*} \sim 0.90$ with the central $68\%$ of the distribution spanning $a_{*} \sim 0.70 - 0.96$, in striking agreement with the estimated spins of the gravitational-wave signal GW231123. These spin values persist up to the highest masses generated by our mechanism, $m_{\rm BH} \sim 10^3\,{\rm M}_\odot$.

Large high-quality samples of HII regions and their parent Giant Molecular Clouds (GMC) are now available for local galaxies. It is therefore possible to investigate links between the CO and H$\alpha$ luminosity functions and whether massive stars form in GMCs of all masses. The CO luminosity functions (LF), representing the distribution of GMC masses, are consistently steeper than the H$\alpha$ luminosity functions. The CO LF invariably steepens in the outer disk where fewer massive GMCs are present beyond the median cloud galactocentric distance. The H$\alpha$ LF also steepens in the outer disk for most of the galaxies examined. Using Salpeter, Kroupa, and Chabrier Initial Mass Functions (IMF) along with stellar mass-luminosity-radius relations, we compute numerically the bolometric luminosity and H$\alpha$ emission from young star clusters. The cluster masses are linked to the GMC mass by assuming that the cluster mass is a constant fraction (3\%) of the parent cloud mass. In particular, results for a fully stochastic IMF are compared to suggestions that very massive stars only form in massive clusters or clouds. Within the limits of the observations -- no small molecular clouds or low-luminosity HII regions can be detected at the typical $\sim 10$~Mpc distance of the sample galaxies -- we find no evidence for a maximum stellar mass which varies with cloud or cluster mass.

Astrophysical luminous objects such as the first stars have not yet formed in the Dark Ages. However, primordial black holes (PBHs) always exist throughout cosmic history since the inflation epoch. During the Dark Ages, PBHs may accrete the ambient gas and release radiation like astrophysical luminous objects, change the cosmic radiation field, the thermal status of the intergalactic medium (IGM), and the hydrogen spin temperature. The accretion rate is modulated by the relic supersonic relative streaming velocities between dark matter (DM) and baryons, imprinting Velocity Acoustic Oscillations (VAOs) features in the 21 cm power spectrum. Such VAOs features could be a promising probe for detecting the PBHs in Dark Ages. We find that even if PBHs comprise only a small fraction of DM, they can generate VAOs wiggles with a relative amplitude up to about 30% in Dark Ages. For example, for PBHs with a mass at recombination of 200 solar masses and mass fraction in the total DM f_PBH,rec around 1e-13 at the recombination era, VAOs features appear at redshift around 20; if f_PBH,rec is around 3e-10, then VAOs features could appear as early as redshift around 40. Moreover, the redshift evolution of the VAOs features exhibits clearly separated stages dominated by inhomogeneous Ly-alpha scattering, and inhomogeneous X-ray heating, respectively. It reflects the characteristics of PBHs (mass and fraction in total DM) and their interactions with the IGM. We also estimate that, the VAOs wiggles at redshift around 20 are detectable for the upcoming SKA-low AA*, while wiggles at redshift around 40 are detectable for an hypothetic lunar surface-based interferometer array in the future.

Polycyclic aromatic hydrocarbons (PAHs) are fundamental to understanding the interstellar medium (ISM) of several astrophysical objects. Normally present in Starburst (SB) galaxies, they have also been more frequently detected in active galaxy nuclei (AGNs), suggesting an inner dusty torus that can shield the radiation from the central black role. In this work, we analyze the 6.2 $\mu$m PAH band of SB-, AGN- and mixed-dominated spectra from 175 IDEOS database galaxies. After fitting of the band, the sources were distributed into the Peeters' A, B and C classes according to their profile peak positions. Class A objects are predominant in 80% of the entire sample, which could indicate the presence of PAHs with nitrogen incorporation. The water ice absorption at 6.0 $\mu$m was also studied in eleven objects, and it affected the PAH band poorly. A prominent second spectral feature after 6.3 $\mu$m is present in ten galaxies. Fitting both PAH profiles at 6.2 $\mu$m changes all the fit results: the first profile is consistently blue-shifted and classified as class A due to the presence of the second component. Further studies are needed to better comprehend these PAH trends in galactic environments.

The organosulfur biosignature gases dimethylsulfide (DMS) and dimethlydisulfide (DMDS) have recently been claimed to be present in the atmosphere of sub-Neptune exoplanet K2-18b, leading to the suggestion of possible extraterrestrial life. Abiotic formation pathways for DMS and DMDS in reducing atmospheres have also been proposed, raising concern over the use of DMS and DMDS as biosignature gases more generally. In this paper we independently test and contrast the proposed abiotic formation pathways for DMS and DMDS using K2-18b as a case study, and explore the wider implications for the atmospheric carbon and sulfur chemistry of hydrogen-rich sub-Neptunes. We demonstrate that one proposed formation pathway is capable of producing observable abundances of abiotic DMS and DMDS, however it depends sensitively on the energy barrier of the limiting step, which remains unmeasured experimentally. In contrast, hydrocarbons including C2H6 are formed abundantly in such atmospheres and offer a plausible alternative explanation to the reported suggestions of organosulfur compounds on K2-18b, having previously been shown to be degenerate observationally. Finally, we demonstrate that sulfur hazes form via the photochemistry of H2S and condense in the atmosphere of K2-18b even at trace abundances. We propose that variation in atmospheric sulfur abundance can explain the diversity of haziness observed across the sub-Neptune population so far with JWST.

We present a statistical analysis for the interaction-induced central concentration of HI gas distributions and its connection with interaction-induced central star formation enhancement, using a large sample of $\sim 10^4$ galaxies from the ALFALFA and SDSS surveys. By adopting the HI profile parameter $K$, an indicator of gas concentration inferred from the integrated 21 cm emission line, we find that galaxies with more centrally concentrated HI (higher $K$ values) or enhanced specific star foramtion rate (sSFR) exhibit significantly stronger clustering and higher probability of hosting a nearby neighbor on scales below $100h^{-1}\mathrm{kpc}$, which is more pronounced in low-mass galaxies. Furthermore, by utilizing the enhancement functions for a sample of galaxy pairs, we directly trace the evolution of HI concentration and sSFR enhancement as a function of projected separation. Our findings indicate that tidal interactions drive a statistical synchrony between the central concentration of atomic gas and the enhancement of central star formation. Gas concentration appears to be a necessary condition for central star formation enhancement in interacting systems at all but the smallest separations. Compared to satellite galaxies, central galaxies exhibit stronger enhancement of gas fraction, gas concentration and sSFR, suggesting the role of environmental regulation.

Gamma-Ray Burst GRB 241030A (z = 1.411) exhibited a bright afterglow (similar to GRB 221009A), detected across gamma-ray, X-ray, UV, and optical bands, providing a probe of GRB afterglow physics. We compiled multi-wavelength observations spanning from a minute to a week after the prompt emission, processing the data through a unified photometry pipeline. We analysed the observations both analytically and using Bayesian inference with two independent models. Our models assume that the afterglow emission arises from the strong forward shock of a laterally structured jet, with possible contributions from synchrotron self-Compton (SSC) scatterings. Our models reproduce X-ray to optical data, favouring a jet propagating into a constant-density interstellar medium, with a viewing angle within the jet core. However, both analyses require parameter values that are extreme compared to expectations from standard theory. In particular, our results imply extremely energetic jets despite regular prompt energy, leading to a very inefficient prompt emission. Furthermore, the jets are inefficient at accelerating particles, with low electron and magnetic energy fractions, leading to significant SSC emission. Our analyses indicate that the jets have large opening angles and propagate in high-density media. If the afterglow is indeed powered by radiation emitted behind a strong forward shock, our results place GRB 241030A within a sub-class of GRBs characterised by extreme kinetic energies, large jet opening angles, and very low prompt emission efficiencies, with strong SSC radiation. These predictions are difficult to reconcile with typical expectations from other GRBs. We therefore suggest that the afterglow of GRB 241030A is not solely powered by forward shock emission.

We present HECATEv2, the second release of the Heraklion Extragalactic Catalogue (HECATE), an all-sky, value-added galaxy catalogue comprising 204733 galaxies from the HyperLEDA database with recession velocity <14000 km/s (D~200 Mpc). This release focuses on qualitative upgrades of the provided information while maintaining the same parent galaxy sample as HECATEv1. Improvements include a new cosmology-based distance framework, expanded and homogenised optical and mid-infrared photometry from SDSS-DR17/NSA, PS1-DR2, and AllWISE, and new quality-control flags for stellar contamination, incorrect photometry, and coordinate inconsistencies. We also extend the galaxy-size coverage and derive stellar population parameters for a substantially larger fraction of the sample. Star-formation rates (SFR) and stellar masses (Mstar) are now available for >70% of galaxies using updated mid-IR/optical calibrations that account for stellar population age and dust attenuation, while gas-phase metallicities are derived for ~90%. Activity classifications are provided for >50% of galaxies based on spectroscopic and/or photometric diagnostics, and supermassive black hole masses for ~86%. In terms of L$_{B}$,L$_{Ks}$,SFR, and Mstar, HECATEv2 is among the most complete local-Universe catalogues with spectroscopic redshifts. We also provide spatial completeness maps as a function of distance and luminosity, highlighting variations across the sky. Compared to other catalogues (e.g. GLADE+, NED-LVS), HECATEv2 offers broader (optical, near- and far-IR photometry, metallicity, activity classifications) or comparable (mid-IR photometry, SFR, Mstar) coverage, making it a robust reference for studies of SMBH-host galaxy connections, gravitational-wave and high-energy transient hosts, population analyses, and rare galaxy subpopulations.

We present results from our ongoing campaign to follow up the satellite candidates from the Identifying Dwarfs of MC Analog GalaxiEs (ID-MAGE) survey. Previously, we published a list of 355 unresolved satellite candidates identified around 36~nearby LMC- and SMC-mass hosts (D$=$4$-$10~Mpc). We present the velocities of 83 satellite candidates from new Green Bank Telescope \hi\ observations, optical long-slit spectra, and the Dark Energy Survey Instrument Data Release 1. Based on their velocities, we identify six candidates as probable satellite galaxies ($6.5\times10^5\leq M_\star/M_\odot\leq1.5\times10^7$) and 77 as background galaxies. Our results underscore the ability of spectroscopic follow-up to effectively separate satellites from background galaxies. Using the refined sample, we update our previously derived estimates for the average satellite population per host and find 1.7$\pm$0.7 (1.0$\pm$0.3) satellites per LMC-mass (SMC-mass) host. Our current satellite sample includes 25 galaxies confirmed by distances or velocities. This set includes the complete satellite populations of three hosts (UGC~04422: zero satellites, UGC~08201: zero satellites, NGC~3432: four satellites), which we compare to simulations and known satellite systems from the literature. Our sample is nearly complete for the most massive satellites (M$_\star > 10^7~M_\odot$). We find these massive satellites have a quenched fraction of 10--25\%, placing them between the $<$5\% quenched fraction of isolated galaxies and the 40--70\% quenched fraction of MW-analog satellites with $10^7~M_\odot < $ M$_\star < 10^8~M_\odot$. This demonstrates the impact that low-mass galaxies have on the evolution of their satellites.

Baryon acoustic oscillation (BAO) analyses usually report the anisotropic shift parameters $\alpha_\perp(z)$ and $\alpha_\parallel(z)$ relative to a fiducial cosmology, and these quantities are primarily used for cosmological parameter inference. Here we show that they can also be used to construct a direct internal consistency test of the background geometry. In particular, we derive a new null test of flat Friedmann-Lemaître-Robertson-Walker (FLRW) geometry written entirely in terms of the reported BAO shift parameters. The test is calibration free: the sound-horizon ratio $r_{\rm d}/r^{\rm fid}_{\rm d}$ cancels identically, so the relation is independent of the absolute BAO scale. We also derive a calibration-free reconstruction of the deceleration parameter $q(z)$ from the radial BAO sector. Applying these results to anisotropic DESI DR2 BAO measurements, we find no evidence for a breakdown of the flat-FLRW distance relation within current uncertainties. Our results show that anisotropic BAO measurements already provide a nontrivial internal geometric consistency test before performing any model fit.

We provide the community with a homogeneous catalogue of small, significant substructures (henceforth NESTs) extracted from the spatial distribution of Young Stellar Objects (YSOs) in a large, consistent sample of star-forming regions. The catalog allows us to explore the relevance of small scale spatial substructure and discuss the interpretation of NESTs as tracers of star formation activity and remnants of the star formation process. We apply our procedure to consistent catalogues of YSOs to obtain NESTs in a sample of star-forming regions. We apply a photometric classification scheme to obtain the evolutionary stage of YSOs and statistically explore the distribution of class 0/I objects as a proxy of recent star formation activity. The region sample is diverse (in distance, size, structure, and global evolutionary stage), and we consequently find different structural properties and star formation histories. Most NESTs in regions with high recent star formation activity show even higher levels of activity. Moreover, the proportion of NESTs with higher activity than the region average increases with the global level of activity of the region. In approximately half of the regions we also find significant spans in the evolutionary stages of the NESTs, consistent with gradients and episodes of star formation. The combination of NESTs with a statistical exploration of the star formation history within each region provides robust and powerful insights into the star formation process. Our results support the role of NESTs as pristine remnants of star formation in highly active regions,stressing the role of fragmentation. The combination of small structures with large scale spatio-evolutionary patterns suggests hyerarchical, prolonged, dynamic, and complex star formation scenarios.

3I/ATLAS is the third confirmed interstellar object passing through the Solar System. In this work, we conduct narrowband radio technosignature search toward 3I/ATLAS using the Five-hundred-meter Aperture Spherical Telescope (FAST) L-band multibeam receiver from October 2025 to January 2026 on 4 separate dates (i.e. Mars closest, perihelion, Earth closest and flew away from Earth, respectively). We carry out frequency-drifting signal searching with signal-to-noise ratio (SNR) over 10 within 1.05-1.45 GHz via \texttt{bliss} pipeline. These signal hits are grouping into event by beam, frequency and drift rate matching, the events are then filtered by cluster analysis and drift rate cut-off. We also characterized the events by their significance in SNR, structure tensor as well as principal component analysis (PCA). No credible narrowband radio technosignature are detected from 3I/ATLAS after visual inspections. The null results place constraints on the presence of transmitters above $2.862\times 10^{-3}$ W. We further introduce a Bayesian inference framework to constraints on the existence probability and characteristic power of hypothetical transmitters using physically motivated priors to bracket plausible transmitter scenarios.

Neptune's present axial tilt of approximately 28 deg. with respect to its orbital plane can be explained by collisions that its primordial core may have experienced with surrounding planetary embryos during the final stages of its formation. Alternatively, Neptune could have attained its present mass solely through pebble accretion, without the formation of nearby planetary embryos. The embryo-collision scenario has the advantage of naturally explaining the large axial tilts observed in the ice giants. To account for these tilts without invoking late-stage catastrophic collisions, an alternative mechanism must be considered. In this work, I propose that Neptune's current axial tilt could result from the interaction between Triton's tidally evolving orbit and Neptune's spin axis, causing it to resonate with the solar system eigenfrequency s_8. Starting from a Triton-mass satellite captured via the binary planetesimal disruption mechanism, I show that orbital evolutions bringing the satellite near Triton's present orbit can induce a spin-s_8 resonance capable of producing a significant axial tilt of Neptune's spin axis. I develop a model for planetary spin-axis evolution using Euler's equations for a rigid body, which is incorporated into classical numerical integrations of the Newtonian equations of motion. I also include a tidal model to account for the satellite's semimajor-axis decay and orbital circularization. Several numerical simulations are performed with this model, including Neptune as the central body, the newly captured satellite, the Sun and the three other giant planets. Increases in Neptune's obliquity are observed, exceeding 50 deg. in some cases. An obliquity above 20 deg. is obtained for roughly 1/3 of the cases. If Neptune initially had a near-zero obliquity, its current value could therefore have been naturally engendered by the tidal evolution of Triton.

Cosmic-ray pitch-angle transport in magnetohydrodynamic (MHD) turbulence is governed by the interplay between magnetic mirroring and gyroresonant scattering. We develop a guiding-center (GC) Langevin model with explicit mirror drift and gyroresonant diffusion to describe the pitch angle evolution. This model accurately captures our test-particle simulation results in three-dimensional MHD turbulence, driven both solenoidally and compressively. We find that magnetic mirroring can drive anomalous pitch-angle diffusion at large pitch angles (including $90^\circ$) with non-Markovian memory effects, which arises from trapping of particles in magnetic wells. Gyroresonant scattering controls the escape rate from these wells. Across $M_{\rm A}$, large-pitch-angle particles are jointly regulated by mirror trapping and gyroresonant escape, exhibiting a transition from anomalous to normal diffusive pitch-angle transport as scattering strengthens, whereas small-pitch-angle particles remain gyroresonance-dominated and diffusive throughout. The pitch angle transport is found to be dominated by the compressible perturbations with marginal influence from Alfvén modes. In compressible turbulence with realistic damping accounted for, transit time damping (TTD) treatment fully recovers mirror interactions.

Traditional weak gravitational lensing shear estimators are carefully calibrated but struggle to fully capture realistic galaxy morphologies, point-spread-function (PSF) effects, blending, and noise in deep surveys, while blindly trained machine learning (ML) models can introduce significant calibration biases. Here we construct a fully D$_4$-equivariant deep neural network for galaxy shape measurement whose architecture enforces symmetry under 90$^{\circ}$ rotations and mirror transformations, and adopt the Analytical Calibration framework (AnaCal) to calibrate the model using its backpropagated gradients. For isolated galaxies in LSST-like single-band simulations, we demonstrate that our approach achieves $\sim$10% lower shape noise than the traditional moment-based Fourier Power Function Shapelets estimator in the high-noise regime, equivalent to a $\sim$20% gain in effective galaxy number density, while simultaneously achieving multiplicative biases consistent with zero across a wide range of noise levels, PSF sizes and ellipticities, and magnitude selection cuts, with all measurements satisfying $|m| {<} 10^{-3}$ (i.e., within the 0.2% LSST requirement) and most at the ${\sim}10^{-4}$ level. We demonstrate this framework on isolated single-band galaxy images with Gaussian noise and known PSF, establishing a rigorous, physics-informed foundation for future extensions of ML-based shear estimation to blended sources and multi-band observations in Stage-IV surveys. All codes and data products will be made publicly available upon acceptance.

Here we present broadband radio modeling of GRB 240205B, using observations with the Australia Telescope Compact Array (ATCA) and the South African MeerKAT radio telescope. Our observations include an automatically triggered early-time ATCA observation that began approximately 13 minutes after the gamma-ray signal and continued for 12 hours, resulting in the earliest detected GRB radio afterglow to date at about 35 minutes post-burst. Following this initial detection, we conducted an extensive radio follow-up campaign for more than 5 months. Although the observations beyond one day post-burst are well described by a standard forward shock model, the observation before one day post-bust reveals an additional synchrotron component, which can be explained as the reverse shock. This component would have been missed without the automated ATCA rapid-response trigger. We find that a combined reverse and forward shock model in a stellar wind medium best describes the radio afterglow. We constrain the spectral breaks due to synchrotron self-absorption and the minimum electron energy, and we use the light-curve peaks to constrain the microphysical parameters. We put GRB 240205B in the context of the growing sample of GRBs with radio detections in the first hours after the gamma-ray trigger. Using our rapid response observation, we estimate the highest model independent constraint on a GRB minimum bulk Lorentz factor of around 100 at about 35 minutes post burst. We also discuss future prospects of detecting similar long GRBs at centimeter wavelengths, as well as potential improvements to future strategies for targeting their radio afterglows.

The nature of the so-called G objects orbiting the Galactic Center remains unresolved. These sources exhibit compact Br$\gamma$ emission, extreme infrared colors, and remarkable dynamical stability through close passages to the central supermassive black hole, challenging conventional interpretations as stars or unbound gas clouds. We investigate the hypothesis that G objects are the remnants of neutron stars that have been converted into low-mass black holes through the capture of primordial black holes, a viable dark-matter candidate. We construct a population-level framework linking the abundance and spatial distribution of these remnants to the neutron-star population, the inner dark-matter density profile, and the primordial black-hole mass and abundance. Within this framework, the observed G-object population and the long-standing deficit of ordinary radio pulsars in the Galactic Center emerge as complementary consequences of the same conversion process. We further identify a suite of observational signatures-across infrared, radio, X-ray, and microlensing channels-that render this scenario empirically testable and distinguishable from stellar-envelope models. Our results show that G objects can act as sensitive probes of compact-object capture physics and of dark matter on sub-galactic scales.

Io's internal heat is generated by Jupiter-driven tidal dissipation and Laplace resonance. This energy partially melts the mantle, but the melt fraction, depth, and spatial distribution of dissipation remain poorly constrained. Tidal deformation is linked to the mantle's physical state via a parametric approach accounting for melting onset depth and latent heat of fusion. Io is modeled as a three-layer body comprising a fluid core, a viscoelastic mantle, and an elastic lithosphere. The degree-2 potential Love number k2 is computed by solving spheroidal oscillation equations with an adapted California Planetary Geophysics Code (CPGC). Mantle properties (viscosity, shear modulus, Andrade parameter beta) are iteratively updated based on the local melt fraction, and mantle compressibility is explicitly incorporated. Results show a self-consistent shallow-mantle enhancement of tidal heating. Melt presence decreases effective viscosity and increases anelasticity, amplifying upper-mantle dissipation. Reproducing the observed real part of k2 in 1D models requires melt fractions below the rheologically critical threshold. Compressibility yields higher values for the real part of k2, reinforcing this upper bound. Mass flux analysis confirms melt percolation capacity exceeds thermodynamic production, implying efficient drainage. The reference Andrade parameter beta strongly influences the imaginary Love numbers (k2, h2, l2) and the predicted libration amplitude. These constraints support a heterogeneous mantle characterized by a magmatic sponge structure rather than a global magma ocean, establishing a quantitative framework linking Io's interior, partial melting, and tidal dissipation to Juno observations.

Spin-orbit obliquity measurements of hot-Jupiter systems constrain giant planet migration and tidal evolution. In binary systems, combining stellar obliquities with the orbit-orbit angle ($\gamma$) between the planetary and stellar companion orbits provides further insight into the dynamical influence of stellar companions. Here we aim to determine the projected obliquity ($\lambda$) of the hot Jupiter TOI-1333Ab ($P\approx4.72$ d, $M_{\rm p}\approx2.4$ M$_{\rm J}$) and place the system in the context of hot-Jupiter migration and tidal realignment in binary systems. We analysed spectroscopic observations obtained during planetary transit to model the Rossiter-McLaughlin effect and derive the projected obliquity. We combined this measurement with published system parameters and constraints on the wide stellar companion orbit to assess plausible migration scenarios. We measure a projected obliquity of $\lambda=-5 \pm 10^\circ$, showing that TOI-1333Ab is well aligned with the stellar spin axis of its F-type host star. The low obliquity and its modest eccentricity ($e=0.073^{+0.092}_{-0.052}$) are consistent with either disc-driven migration or high-eccentricity migration followed by efficient tidal circularisation and realignment. With an effective temperature of $6274\pm94$ K, the host star lies above the canonical Kraft break where the systems are frequently misaligned. Despite this, we find the system to be well aligned. In comparison with other planetary systems in binaries, TOI-1333 occupies a relatively isolated region in projected obliquity-orbit-orbit angle ($\gamma=81.5\pm1.1^\circ$) space, making it a valuable system for studying the interplay between migration, tides, and stellar companions.

We present a high-efficiency, high resolution on-chip filterbank spectrometer designed for line intensity mapping and broadband wave-like dark matter searches. This spectrometer maximizes sensitivity to the faint, aggregate cosmic signal for line intensity mapping while providing redshift precision necessary to resolve the clustering of large-scale structure. For broadband dark matter searches, it partitions broadband signals into narrow frequency bins that directly correspond to dark matter masses, allowing signal frequency to be mapped with high precision to the particle's mass. Existing superconducting filterbank architectures used by the mm-wave community are limited by a 50\% inherent efficiency limit and are highly sensitive to thin-film dielectric loss. The design presented in this paper addresses these bottlenecks by eliminating the termination resistor and employing a niobium-on-silicon coplanar waveguide resonant structures for the filterbanks. Sonnet electromagnetic simulations of a 10-channel device around 90 GHz demonstrates >74\% per channel efficiency and a resolving power of $R=1211\pm105$. Sensitivity analyses confirm that the design is robust against typical fabrication uncertainties with the exception of dielectric thickness, providing a scalable technology solution for the next generation of millimeter-wave spectroscopic experiments.

We assess the impact of non-ideal, continuously rotating half-wave plates (HWPs) on cosmic microwave background (CMB) polarization measurements targeting large angular scale signal. Such hardware solutions are used in or planned for multiple modern CMB efforts, both ground-based, for instance, small aperture telescopes of Simons Observatory or satellite borne, such as LiteBIRD. Using a frequency-dependent parametric model based on the Mueller matrix formalism, we characterize the induced mixing of Stokes parameters. Through end-to-end simulations, we propagate these effects from time-ordered data to cosmology via map-making and component-separation stages, quantifying their impact on the $B$-modes power spectrum and the tensor-to-scalar ratio, $r$. Our analysis shows that neglecting the frequency dependence of a three-layer HWP gives rise to significant polarization leakage, biases foreground spectral parameters, and leads to residual contamination in the recovered CMB maps. To mitigate these effects, we investigate multiple analysis strategies progressively incorporating a more complete description of the instrumental response. At the map-making level, this requires generalizing the standard pointing matrix to account for the full time- and frequency-dependent instrumental response. We find that standard HWP models, reduce the biases only down to $r \sim 10^{-2}$, while a more advanced approach based on a generalization of both map-making and component separation, implemented using JAX, can suppress it down to $r \sim 7 \times 10^{-4}$. Finally, we extend this approach to a time-domain component-separation, enabling a statistically consistent treatment of instrumental response in the presence of time-domain features. We demonstrate its feasibility and validate it by performing a full end-to-end analysis, recovering results in good agreement with the map-based ones.

We present the first semi-empirical strong-line calibrations to determine metallicity in Active Galactic Nuclei (AGNs) that use the directly observable X-ray luminosity ($Ł_{\rm X}$) instead of the dimensionless ionization parameter ($U$). The calibrations are derived from an extensive grid of photoionization models computed with the {\sc Cloudy} code, which are compared with observational data of Seyfert nuclei from the Burst Alert Telescope (BAT) AGN Spectroscopic Survey (BASS). In this first paper, we develop new calibrations for two key optical metallicity diagnostics based on the $N2$ and $O3N2$ indices, which are valid in a metallicity range of $8.0 \lesssim \logOH \lesssim 9.1\, {\rm or}\, 0.2 \lesssim (Z/Z_{\odot}) \lesssim 2.6$, with precision of $1\sigma \approx 0.22$ dex ($N2$) and $\approx 0.20$ dex ($O3N2$). We systematically investigate the influence of the AGN spectral index $(\aox)$, narrow-line region (NLR) gas density (\Ne), the characteristic peak temperature of the Big Blue Bump $(T_{\rm BB})$, and $Ł_{\rm X}$. We find a strong, opposing secondary dependence on $Ł_{\rm X}$ for both indices. We demonstrate that neglecting this parameter overlooks systematic offsets intrinsic to the diagnostics, leading to metallicity errors of up to $\sim 1.0$ dex, particularly for the least and most luminous sources. This framework offers a more precise characterization of chemical enrichment in the NLRs of AGNs by leveraging their intrinsic X-ray emission to mitigate these systematic biases.

We present a comprehensive chemical and age orbital cartography of the Galactic high-$\alpha$ disc using subgiant stars with precise ages, element abundances, and full phase-space information from the \textsl{LAMOST--Gaia} data set. Specifically, we map how average [Fe/H], [$\alpha$/Fe], and age vary across present-day kinematic and orbital coordinates. We analyse the data in full and across mono-abundance populations to measure element abundance-orbital and age-orbital gradients across orbital actions and angular-momenta. Our results show that the high-$\alpha$ disc exhibits clear and coherent gradients in [Fe/H], [$\alpha$/Fe], and age with orbits; these gradients are much stronger and sharper in orbital space than in present-day kinematics, showing that orbital diagnostics recover the intrinsic disc structure of old disc populations more effectively than instantaneous kinematic coordinates. We find that older high-$\alpha$ populations display qualitatively similar element abundance--orbital and age--orbital trends to stars in the low-$\alpha$ disc, although the high-$\alpha$ gradients are generally shallower. The presence of these ordered correlations indicates that the old high-$\alpha$ disc is structured, and preserved a strong fossil record of its early assembly despite the Milky Way's subsequent accretion history. This result implies that later mergers did not fully erase the chemical-orbital and age-orbital structure imprinted during the high-$\alpha$ disc's earliest formation epoch. All together, these findings indicate that the Galactic high-$\alpha$ disc formed mainly through inside-out and upside-down growth.

Quantum field theories provide fundamental models of complex interacting systems, from high-energy physics and cosmology to condensed matter. However, solving these models in non-perturbative and dynamical regimes is often extremely challenging, particularly in more than one spatial dimension. Analog simulation using tunable synthetic quantum systems can both verify existing theoretical predictions and lead to new physical insights. Here, we realize quantum simulation of massive relativistic fields in $2+1$ dimensions (two spatial dimensions and time), using two coherently coupled spin components in a uniform two-dimensional Bose-Einstein condensate. Specifically, we encode the paradigmatic sine-Gordon model in the field describing the relative phase, $\phi$, of the two components. We show that, in the perturbative regime, collective field excitations exhibit a relativistic dispersion with a tuneable mass gap. We also observe explicitly non-perturbative phenomena, including the existence of topological domain walls across which $\phi$ rapidly winds by $2\pi$. Our work opens possibilities for studies of cosmologically relevant phenomena including preheating, dynamics of topological defects, and relativistic false-vacuum decay.

Thermal misalignment provides an alternative to the standard misalignment mechanism for the cosmological production of scalar dark matter. In this framework, feeble couplings to particles in the thermal bath generate a finite-temperature potential that drives the scalar towards large field values early in the radiation era, dynamically inducing the misalignment before the onset of scalar oscillations. As a result, the relic abundance is controlled primarily by particle masses and couplings rather than the initial field value. As a light spectator field, the scalar acquires inflationary fluctuations that are uncorrelated with the adiabatic curvature mode, generically sourcing isocurvature perturbations. We show that, unlike standard misalignment, where light scalars are strongly constrained by cosmic microwave background bounds on dark matter isocurvature for high-scale inflation, thermal misalignment can naturally suppress the isocurvature signal. This occurs through a novel late-time phase offset between the background zero mode and the superhorizon perturbations, which reduces the final dark matter density contrast. Thermal misalignment therefore provides a new and generic route to isocurvature-safe scalar dark matter.

We propose the idea that our Universe is a realization among different possible branches, which can be observationally tested through the modified dispersion relation of the gravitational waves. We achieve this through a framework of spatially constrained vector fields. We show that the simplest realizations of such theories in flat and cosmological spacetimes do not introduce new propagating modes, but they give rise to tensor perturbations that differ from those of standard general relativity. We further show that such theories admit stealth black hole solutions, and we recover weak gravitational potentials, thus passing the solar system experiments. Finally, we discuss the implications of such theories and propose further generalizations.

Tentative evidence suggests that the cores of massive neutron stars consist of deconfined quark matter. We argue that the formation of such a quark matter core during a galactic supernova could be accompanied by the emission of gravitational waves in the MHz band. These signals constitute a new target for high-frequency gravitational wave detectors, demonstrating that such detectors may offer unique opportunities for testing quantum chromodynamics in an otherwise inaccessible regime.

The meV mass range has emerged as a focal point in axion physics, where advances in theory, cosmology, astrophysics, and experimental techniques converge. Axions in this mass range are theoretically well motivated, can arise in ultraviolet-complete models, and can have significant cosmological impacts as dark matter or dark radiation. In parallel, their efficient production in stellar and supernova environments provides powerful astrophysical probes. Here, we provide a comprehensive overview of meV axions across these domains, highlighting both established results and open questions. We discuss the theoretical underpinnings of meV axions, their cosmological and astrophysical signatures, and the diverse experimental strategies -- ranging from helioscopes and haloscopes to quasiparticle systems and large-volume Cherenkov detectors -- that aim to explore this regime. The convergence of these approaches emphasizes the pivotal role of the meV mass range for axion discovery in the coming years, identifying meV axions as a key probe for testing beyond-Standard-Model physics. This review document is the direct outcome of the discussions at the dedicated workshop "The meV Mass Axion Frontier: Challenges and Opportunities", held at Laboratori Nazionali di Frascati (IT) on 27--28 October 2025, and organized by the EU funded COST Action "Cosmic WISPers in the Dark Universe: Theory, astrophysics, and experiments" (CA21106, this https URL). Its aim is to provide an overview of current efforts in meV axion research, their motivations, and the research goals that animate the community involved in this search.

In this paper, we pursue an investigation of the consequences of a mixing between supersymmetric partners - the photon and photino - analogous to the so-called Primakoff effect, but induced by a Lorentz-symmetry violating (LSV) fermionic-condensate background. In our framework, the LSV parameters are introduced as members of a non-dynamical superfield. As a consequence, we show that naturally there appears a mixing term between the gauge boson and the gaugino, which can be readily seen in the superspace/superfield approach. We inspect the kinetic photon-photino mixing matrix in the scenario of stellar physics which we apply our results to. Bounds on the strength of the fermionic LSV background are can be set by invoking the energy loss argument and the solar data we adopt.

Spacetime is foamy due to quantum fluctuations. Various gedanken experiments show that distances fluctuate by amounts consistent with the holographic principle, hence the name "holographic quantum foam" (HQF). One important prediction of HQF is that necessarily there exists a dark sector in the universe. The resulting cosmology is found (at least qualitatively) to be consistent with observations. Interestingly the quanta of the dark sector are found not to obey the familiar (fermionic or bosonic) statistics, but the exotic statistics known as infinite statistics (or quantum Boltzmann statistics). The most important challenge now is to check if HQF is consistent with experiments/observations. One way is to look for observational evidence of blurred distant point-sources due to physics at the Planck scale. For over two decades it has been debated whether those tiny inherent uncertainties in time and path-length can accumulate in transiting electromagnetic wavefronts from quasars and Gamma-Ray Bursts (GRBs). But a recent event is special: GRB221009A was extremely bright and energetic. That allowed follow-up across the whole spectrum from the optical/near-infrared through to X-rays, and including the highest-ever-recorded energy gamma-rays; all consistent with blurring by HQF. Those data, and a calculation of the HQF-widened point-spread function (PSF) for real telescopes viewing a GRB are presented.

We study the cosmological implications of the minimal non-linear realisation of scale invariance within the Standard Model (SM). This framework provides a technically natural explanation for the hierarchy between the Planck scale and the electroweak scale and introduces only a light, feebly coupled dilaton field beyond the SM particles. Although the model is almost indistinguishable from the minimal SM at low energies, its cosmological consequences differ dramatically. In particular, the electroweak Higgs field remains trapped in the symmetric phase until the Universe cools to very low temperatures, $T_c^{(\chi)}\sim 28$ MeV, where the first-order QCD chiral symmetry-breaking phase transition triggers the electroweak phase transition. This scenario offers intriguing possibilities for the production of primordial black holes, low-frequency gravitational waves, and multi-quark and lepton nuggets, which we explore in some detail using simplified approximations.

We review two types of fractional Kaehler potentials $K$ which reduce, along the inflationary path, to the form $N/(1-\phi^{q_{\rm M}})^{p}$ with $q_{\rm M}=1$ or $2$ and $0.1\leq p\leq10$. Their coexistence, within a non-linear sigma model, with chaotic inflationary potentials of the form $\phi^n$ (where $n=2$ or $4$) determines, independently from $q_{\rm M}$ and $n$, a class of $(p,N)$-inflationary attractors which leads to observables compatible with the ACT DR6. An implementation of these models in the context of supergravity can be also achieved by introducing two chiral superfields and a monomial superpotential, linear with respect to the inflaton-accompanying field, and supplementing the $K$'s above with a shift symmetry. Although inflation is attained for subplanckian inflaton values, the tensor-to-scalar ratio obtained for certain $N$ values can be possibly observable in the near future.

We investigate turbulent spin-down dynamics in self-gravitating Bose-Einstein condensates, comparing purely bosonic and axionic (higher-order interacting) systems. Through simulations of the Gross-Pitaevskii-Poisson system, we study condensates pinned to a crust potential undergoing rapid rotation slowdown. We find that axionic condensates exhibit more uniform density profiles and smaller sizes compared to their bosonic counterparts for similar interaction strengths, which facilitates earlier vortex entry. The sudden spin-down triggers vortex depinning and a turbulent cascade. For comparable sizes, both systems exhibit a short-lived Kolmogorov energy cascade ($k^{-5/3}$ scaling) followed by a transition to Vinen turbulence ($k^{-1}$ scaling). Crucially, their responses diverge with increasing interaction strength (and thus condensate size): the axionic system increasingly deviates from Kolmogorov scaling because of enhanced vortex retention, a trend quantitatively confirmed by analyzing the vortex fraction and its dependence on the final rotation frequency. Spectral analysis reveals that the growth of incompressible energy is primarily driven by quantum pressure during vortex detachment, rather than by compressible flows. The compressible spectrum shows thermalization ($k$ scaling). Our results demonstrate how distinct nonlinearities govern vortex dynamics and turbulent dissipation in self-gravitating quantum fluids.

The worldwide LIGO-Virgo-KAGRA network of gravitational-wave (GW) detectors continues to increase in sensitivity, thus increasing the quantity and quality of the detected GW signals from compact binary coalescences. These signals allow us to perform ever-more sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. This paper is the first of three, where we present the results of a suite of tests of GR using the binary signals included in the fourth GW Transient Catalog (GWTC-4.0), i.e., up to and including the first part of the fourth observing run of the detectors (O4a). We restrict our analysis to the 91 confident signals, henceforth called events, that were measured by at least two detectors, and have false alarm rates $\le 10^{-3} \mathrm{yr}^{-1}$. These include 42 events from O4a. This first paper presents an overview of the methods, selection of events and GR tests, and serves as a guidemap for all three papers. Here we focus on the four general tests of consistency, where we find no evidence for deviations from our models. Specifically, for all the events considered, we find consistency of the residuals with noise. The final mass and final spin as inferred from the low- and high-frequency parts of the waveform are consistent with each other. We also find no evidence for deviations from the GR predictions for the amplitudes of subdominant GW multipole moments, or for non-GR modes of polarization. We thus find that GR, without new physics beyond it, is still consistent with these GW events. The results of the two additional papers in this trio also find overall consistency with vacuum GR, with more than 90% of the events being consistent with GR at the 90% credible level. While one of the ringdown analyses finds the GR value in the tails for its combined results, this may be due in part to catalog variance.

In this second of three papers on tests of general relativity (GR) applied to the compact binary coalescence signals in the fourth Gravitational-Wave Transient Catalog (GWTC-4.0), we present the results of the parameterized tests of GR and constraints on line-of-sight acceleration. We include events up to and including the first part of the fourth observing run (O4a) of the LIGO Virgo KAGRA detectors. As in the other two papers in this series, we restrict our analysis to the 42 confident signals, measured by at least two detectors, that have false alarm rates $\le 10^{-3} \mathrm{yr}^{-1}$ from O4a, in addition to the 49 such events from previous observing runs. This paper focuses on the eight tests that constrain parameterized deviations from the expected GR (or unaccelerated) values. These include modifications of post-Newtonian (PN) parameters, spin-induced quadrupole moments different from those of a binary black hole, and possible dispersive or birefringent propagation effects. Overall, we find no evidence for physics beyond GR, for spin-induced quadrupole moments different from those of a Kerr black hole in GR, or for line of sight acceleration, with more than 90% of the events including the null result (no deviation) within their 90% credible intervals. We discuss possible systematics affecting the other events and tests, even though they are statistically not surprising, given noise. We improve the bounds on deviations from the GR PN coefficients by factors of 1.2-5.5 and provide illustrative translations to constraints on some modified theories. Also, we update the bound on the mass of the graviton, at 90% credibility, to $m_g \leq 1.92\times 10^{-23} \mathrm{eV}/c^2$. Thus, we see that GR holds, and many of the bounds on possible deviations derived from our events are the best to date.

This is the third paper of the set recording the results of the suite of tests of general relativity (GR) performed on the signals from the fourth Gravitational-Wave Transient Catalog (GWTC-4.0), where we focus on the remnants of the binary mergers. We examine for the first time 42 events from the first part of the fourth observing run of the LIGO-Virgo-KAGRA detectors, alongside events from the previous observation runs, restricting our analysis to the confident signals, which were measured in at least two detectors and that have false alarm rates $\le 10^{-3} \mathrm{yr}^{-1}$. This paper focuses on seven tests of the coalescence remnants. Three of these are tests of the ringdown and its consistency with the expected quasinormal mode spectrum of a Kerr black hole. Specifically, two tests analyze just the ringdown in the time domain, and the third test analyzes the entire signal in the frequency domain. Four tests allow for the existence of possible echoes arriving after the end of the ringdown, which are not expected in GR. We find overall consistency of the remnants with GR. When combining events by multiplying likelihoods (hierarchically), one analysis finds that the GR prediction lies at the boundary of the $98.6^{+1.4}_{-9.4}\%$ ($99.3^{+0.7}_{-4.5}\%$) credible region, an increase from $93.8^{+6.1}_{-20.0}\%$ ($94.9^{+4.4}_{-18.2}\%$) for GWTC-3.0. Here the ranges of values comes from bootstrapping to account for the finite number of events analyzed and suggest that some of the apparently significant deviation could be attributed to variance due to the finite catalog. Since the significance also decreases to 92.2% (96.2%) when including the more recent very loud event GW250114, there is no strong evidence for a GR deviation. We find no evidence for post-merger echoes in the events that were analyzed. (Abridged)

Binary neutron star mergers and proto-neutron stars provide unique environments where dense matter is hot, lepton rich, and potentially undergoes a transition from hadronic to deconfined quark matter. We investigate the thermodynamics and stellar properties of hybrid matter under such conditions. The hadronic phase is described within a covariant density functional framework, while the quark phase is modeled using a Nambu-Jona-Lasinio (NJL) model that includes repulsive vector interactions, the axial $U_A(1)$-breaking 't Hooft determinant interaction, and two-flavor color-superconducting (2SC) pairing. The phase transition between hadronic and quark matter is constructed using a mixed-phase prescription that enforces baryon and lepton number conservation, allowing us to follow thermodynamic trajectories at fixed entropy per baryon and fixed lepton fraction. We analyze the phase structure of dense matter at finite temperature and study the composition of the hadronic, mixed, and quark phases in both neutrino-trapped and neutrino-free regimes. Our results show that neutrino trapping significantly modifies the particle composition and shifts the onset of deconfinement to higher densities. Using the resulting equations of state, we compute static stellar configurations and examine the influence of temperature and lepton content on the mass-radius relation of hybrid stars. Hot, neutrino-rich configurations are found to have larger radii and slightly higher maximum masses than their cold counterparts.

The production of very-high-energy (VHE, $E_{\gamma} \gtrsim 100$ GeV) gamma rays resulting from the scattering of high-energy cosmic-ray protons off axion-like particles (ALPs) populating the dark matter halo of the Milky Way is investigated. By employing the latest instrument response functions for current and future facilities, we demonstrate that ground-based VHE gamma-ray observatories, such as H.E.S.S., CTAO, and SWGO, provide a promising and complementary avenue to probe the yet uncharted ALP-proton coupling $g_{ap}$. Our results show that these experiments can reach sensitivity to couplings above $10^{-2}$ in the $1 - 10^{8}$ eV ALP mass range, a region that remains largely unexplored by supernova and neutron star cooling observations. Interestingly, we demonstrate that this search channel is capable of probing QCD axion dark matter models, assuming two benchmark models for it: the Kim-Shifman-Vainshtein-Zakharov (KSVZ) Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) models, specifically within the MeV mass range. These findings highlight the potential of VHE gamma-ray astronomy to provide unique constraints on the interaction between ALPs and the baryonic sector.

We perform a search for an X-ray monochromatic line arising from dark matter (DM) decay in the halo of the Large Magellanic Cloud. An emission line can be expected from two well-motivated DM candidates: sterile neturinos and axion-like particles (ALPs). We analyze the eROSITA-DE DR1 datasets in the energy range between 1 and 9 keV. No evidence for a DM line is found, and we set lower limits on the DM lifetime. We then recast these bounds into upper limits on the active-sterile neutrino mixing angle $\sin^2(2\theta)$ and on the ALP to photon coupling $g_{a\gamma}$, for DM masses between 2 and 18 keV. These results set new strong constraints for masses below 5 keV.

We investigate the phenomenon of black hole superradiance in the presence of torsion within the framework of Poincaré gauge theory. In particular, in contrast to the classical approach of General Relativity, we show that the inclusion of torsion in the space-time geometry enables the energy extraction from rotating black holes by Dirac fermions via chiral asymmetry, while preserving the Pauli exclusion principle.

We review recent progress in the understanding of the preheating stage of Higgs inflation formulated within the Einstein-Cartan framework of gravity. This setup smoothly interpolates between the metric and Palatini formulations of the theory, leading to a distinctive phenomenology in an intermediate regime. Following the end of inflation, the Higgs field undergoes a non-trivial out-of-equilibrium evolution driven by tachyonic instabilities and nonlinear self-interactions, which fragment the inflaton condensate and give rise to well-localized oscillon configurations. While early studies suggested the formation of long-lived oscillons and the possibility of an extended matter-dominated phase, more recent analyses show that self-interactions at small field values render these objects transient, eventually triggering their decay and the onset of radiation domination. We discuss the implications of this dynamics for the thermal history of the Universe, the inflationary observables, and the generation of stochastic gravitational waves.

The way in which the large-scale cosmic environment affects galactic properties is not yet understood. Dark matter halos, which embed galaxies, initially evolve following linear theory. Their subsequent evolution is driven by non-linear structure formation in the halo region and in its outer environment. In this work, we present the first study where we explicitly control the linear part of the evolution of the halo, thus revealing the role of non-linear effects on halo formation. We focus specifically on the effect of proximity to a large cosmological filament. We employ the splicing method to keep fixed the initial density, velocity, and potential fields where a halo will form while changing its outer environment, from an isolated state to one where the halo is near a large filament. In the regime of Milky Way-mass halos, we find that mass and virial radius of such halos are not affected by even drastic changes of environment, whereas halo spin and shape orientation with respect to a massive filament is largely impacted, with fluctuations of up to 80 % around the mean value. Our results suggest that halo orientation and shape cannot be predicted accurately from a local analysis in the initial conditions alone. This has direct consequences on the modeling of intrinsic alignment for cosmic shear surveys, like Euclid. Our results highlight that non-linear couplings to the large-scale environment may have an amplitude comparable to linear effects, and should thus be treated explicitly in analytical models of dark matter halo formation.

DESI has reported a dynamical dark energy (DE) signal based on the $w_0 w_a$CDM model that is in conflict with Hubble tension. Recalling that the combination of DESI DR1 BAO and DR1 full-shape (FS) modeling are consistent with $\Lambda$CDM, in this letter we comment on the status of fluctuations in DR1 BAO documented in \cite{DESI:2024mwx, Colgain:2024xqj} in the DR2 update. In particular, we note that neither DR1 BAO nor DR2 BAO nor DR2 BAO+CMB confronted to the $w_0 w_a$CDM model with relaxed model parameter priors confirm late-time accelerated expansion today. Translating DESI BAO constraints into flat $\Lambda$CDM constraints, we observe that the LRG1 constraint remains the most prominent outlier, a distinction now held jointly with ELG1, LRG2 switches from smaller to larger $\Omega_m$ values relative to Planck-$\Lambda$CDM, and ELG data drive the relatively low $\Omega_m$ in the full DR2 BAO. We observe that one cannot restore $w_0 = -1$ within one $1 \sigma$ by removing either LRG1 or ELG1 or LRG2, but LRG2 in DR2, in contrast to LRG1 in DR1, now has the greatest bearing on $w_0 > -1$. We conclude that BAO has yet to stabilise, but the general trend is towards greater consistency with DESI DR1 FS modeling results, where there may be no dynamical DE signal in DESI data alone.

The latest results from the Atacama Cosmology Telescope Collaboration have moved the preferred perturbation spectral index $n_{\rm s}$ closer towards one. We reanalyse constraints on the simplest version of the curvaton model, using $n_{\rm s}$ and the tensor-to-scalar ratio $r$. We show that in the massless curvaton case the model gives the same locus of $n_{\rm s}$-$r$ predictions as the power-law model $V(\phi) \propto \phi^p$, but with a different and more elegant physical interpretation. The model gives an excellent account of current observational data, both in the regime where the curvaton dominates the observed perturbations and with an admixture of inflaton-originated perturbations of up to about one-third of the total. Addition of the newest South Pole Telescope SPT-3G D1 dataset maintains this conclusion in favor of the curvaton model. Forthcoming large-scale structure surveys, such as that of the recently-launched SPHEREx probe, could discriminate between curvaton and single-field inflation models using the bispectrum.

The observation of neutron star mergers with gravitational waves (GWs) has provided a new method to constrain the dense-matter equation of state (EOS) and to better understand its nuclear physics. However, inferring nuclear microphysics from GW observations necessitates the sampling of EOS model parameters that serve as input for each EOS used during the GW data analysis. The sampling of the EOS parameters requires solving the Tolman-Oppenheimer-Volkoff (TOV) equations a large number of times -- a process that slows down each likelihood evaluation in the analysis on the order of a few seconds. Here, we employ emulators for the TOV equations built using multilayer perceptron neural networks to enable direct inference of nuclear EOS parameters from GW strain data. Our emulators allow us to rapidly solve the TOV equations, taking in EOS parameters and outputting the associated tidal deformability of a neutron star in only a few tens of milliseconds. We implement these emulators in \texttt{PyCBC} to directly infer the EOS parameters using the event GW170817, providing posteriors on these parameters informed solely by GWs. We benchmark these runs against analyses performed using the full TOV solver and find that the emulators achieve speed ups of nearly \emph{two orders of magnitude}, with negligible differences in the recovered posteriors. Additionally, we constrain the slope and curvature of the symmetry energy at the 90\% upper credible interval to be $L_{\rm sym}\lesssim106$ MeV and $K_{\rm sym}\lesssim26$ MeV.

Joint observations of extreme blazars with Fermi Large Area Telescope (LAT) and Imaging Atmospheric Cherenkov telescopes (IACT) have been previously used to derive lower bounds on intergalactic magnetic field (IGMF). We update these previous bounds using a set of extreme blazars that are detected in the Very-High-Energy (VHE, photon energies above 100 GeV) band by both Fermi/LAT and IACTs. We measure IGMF-dependent suppression of secondary delayed gamma-ray flux from electron-positron pairs deposited in the intergalactic medium by VHE gamma-rays interacting with Extragalactic Background Light. From overall 22 extreme blazars detected by Fermi/LAT and IACTs in the VHE band, seven have their spectral characteristics inconsistent with the possibility of zero magnetic field along their lines of sight, even under the most restrictive assumption that the sources have only switched on at the start of VHE band observations. Adopting this assumption, we derive a "conservative" lower bound on the IGMF strength at the level of 2e-17 G. The tightest bound is imposed by the signal of 1ES 0502+675, a source that has not been considered in the IGMF analysis before. Our bound is comparable to the bound derived by MAGIC collaboration, but is weaker than that previously derived from analysis of Fermi/LAT and HESS telescope data, even though our dataset includes that data. We clarify the origin of this discrepancy.

Various ``cosmic shorelines" have been proposed to delineate which planets have atmospheres. The fates of individual planet atmospheres may be set by a complex sea of growth and loss processes, driven by unmeasurable environmental factors or unknown historical events. Yet, defining population-level boundaries helps illuminate which processes matter and identify high-priority targets for future atmospheric searches. Here, we provide a statistical framework for inferring the position, shape, and fuzziness of an instellation-based cosmic shoreline, defined in the three-dimensional space of planet escape velocity, planet bolometric flux received, and host star luminosity. We circumvent the need to estimate individual host stars' historical X-ray and extreme ultraviolet fluences by including luminosity in the definition of the shoreline, explicitly modeling how sharply such drivers of atmospheric escape intensify toward lower-luminosity M dwarf stars and marginalizing over the associated uncertainties. Using Solar System and exoplanet atmospheric constraints, under the assumption that one planar boundary applies across a wide parameter space, we find the critical flux threshold for atmospheres scales with escape velocity with a power-law index of $p=5.9_{-0.43}^{+0.61}$, steeper than the canonical literature slope of $p=4$, and scales with stellar luminosity with a power-law index of $q=1.17_{-0.20}^{+0.28}$, steep enough to disfavor atmospheres on Earth-sized planets out to the habitable zone for stars less luminous than $\log_{10} (L_\star/L_\odot) = -2.23_{-0.21}^{+0.18}$ (roughly spectral type M4V). This model provides quantitative predictions for the probability any planet may have an atmosphere, which can be rigorously tested by upcoming JWST Rocky Worlds observations.

We present a model-independent null test of the late-time cosmological response to a reduced sound horizon, as typically required by early-universe solutions to the Hubble tension. In this approach, we phenomenologically impose a shorter sound horizon without modeling early-universe physics to isolate its impact on late-time dark energy inference. Using baryon acoustic oscillations (BAO), supernovae (SN), big bang nucleosynthesis (BBN), and local $H_0$ data, while explicitly avoiding CMB anisotropies, we examine how this calibration shift propagates into constraints on the dark energy equation of state. We find that lowering $r_d$ systematically drives the $w_0$-$w_a$ posterior toward less dynamical, quintessence-like behavior, bringing it closer to $\Lambda$CDM. This result underscores that some of the apparent evidence for evolving or phantom-like dark energy may reflect early-universe assumptions rather than genuine late-time dynamics. More broadly, our analysis highlights the importance of carefully disentangling calibration effects from physical evolution in interpreting forthcoming results from DESI and future surveys.

The majority of massive stars are part of binary systems that may interact during their evolution. However, not many RSGs are known binaries, and only a few have constrained orbital parameters. We search the available TESS photometry for eclipsing companions of RSGs. We focus on the best candidate, VV Cephei type binary KQ Pup, which is made up of a RSG, KQ Pup A, and a B-type companion, KQ Pup B (orbital period of 26 yr). We use photometry, spectroscopy, and newly taken interferometric data with VLTI-GRAVITY. Using TESS, we discovered eclipses with a period of $17.2596 \: \rm d$, associated with KQ Pup B, making it a Ba+Bb binary. The detection of the hydrogen Br$\gamma$ line with VLTI-GRAVITY enabled us to track the orbital motion of the Ba+Bb pair relative to A and determine the astrometric orbit of A+B. The dynamical masses agree with independent estimates from asteroseismology and evolutionary models. The results give a mass of $ \sim 10 \: \rm M_{\odot} $ for the RSG KQ Pup A and $ \sim 14 \: \rm M_{\odot} $ for the sum of the hot components Ba+Bb. We determined an orbital parallax of $\pi = 1.24^{+0.05}_{-0.04}\, \rm mas $, which is the first such parallax measurement for a RSG. KQ Pup represents a unique demonstration of mass transfer mechanism in wide eccentric RSG systems. The variability of Balmer emission lines and the detection of Br$\gamma$ are a strong signature of accretion to Ba+Bb near periastron. With the RSG filling its Roche lobe only by $\sim 70\%$ at periastron, the mass transfer is instead driven by accretion from its extended atmosphere via the Wind Roche Lobe Overflow. The accretion disk dissipates by apastron. Overall, we discovered that several previously assumed RSG binaries host eclipsing inner systems, corresponding to $\sim 10 \%$ of all known Galactic RSG binaries. This suggests that many of the remaining RSG binaries may also be hierarchical triples.

We study the inflationary gravitational wave background induced by Abelian gauge fields generated by non-minimal kinetic and axial couplings to the inflaton. We show that, up to slow-roll corrections, for coupling functions that share the same dependence on conformal time, the gravitational wave spectrum is nearly scale invariant. We also derive its amplitude for generic gauge field coupling parameters, within the slow-roll approximation. The coupling values and the scale of inflation for which the induced gravitational wave background is observable, while ensuring that back-reaction on the inflationary dynamics remains negligible, are calculated. We find that a sizeable axial coupling can boost this secondary gravitational wave signal above the standard inflationary background. In the course of our analysis, we also show how to analytically match tensor perturbations across an arbitrary number of eras with different equations of state.

The appearance of highly anisotropic planes of satellites around the Milky Way and other galaxies was long considered a challenge to the standard cosmological model. Some recent simulations have found flattened satellite systems to be common, but these have been described as either "transient", short-lived alignments, or "persistent", long-lived structures. Here we analyse Milky Way analogue systems in the cosmological simulation TNG-50 to resolve this apparent contradiction. We show that, as the satellite populations of individual hosts rapidly change, the observed spatial anisotropies of their satellite systems are invariably short-lived, with lifetimes of no more than a few hundred million years. However, when the progenitors of the same satellites are traced backwards, we find examples where those identified to form a plane at the present day have retained spatial coherence over several billion years. The two ostensibly conflicting predictions for the lifetimes of satellite planes can be reconciled as two perspectives on the same phenomenon.

We present results of Pop. III formation in the MEGATRON suite of simulations, which self-consistently follows radiation and non-equilibrium chemistry, and resolves gas at near-pc resolution of a Milky Way-mass halo at Cosmic Dawn. While the very first Pop. III stars form in halos with masses well below the atomic cooling limit, whose cooling is dominated by molecular hydrogen, the majority of Pop. III stars form in more massive systems above the $10^4$~K atomic cooling threshold. The shift in cooling regime of halos hosting new Pop. III stars occurs within $100$ Myr of the first Pop. III star as the Lyman-Werner (LW) background rapidly increases to $10^{-21}\,\rm erg\,s^{-1}\,cm^{-2}\,Hz^{-1}\,sr^{-1}$. We find that the global Pop. III star formation rate stabilizes to a value of $10^{-3}\,\rm M_\odot\,yr^{-1}$ at $z=20$. Among the three processes that quench Pop. III star formation in mini-halos, the LW background, gas starvation, and external chemical enrichment, the LW background is most important. A small fraction of haloes undergo multiple episodes of Pop. III star formation when the earlier forming stars all directly collapse to black holes. If the halos become massive enough, they can form up to $\sim100$ Pop. III stars in a single burst, which may be observable by JWST with moderate gravitational lensing. Pop. III stars form at a wide range of distances from UV-bright galaxies, with only $0.06\%$ of Pop. III stars forming within the virial radius of galaxies with $M_{\rm UV} < -17$. Finally, by tracking Pop. III star remnants down to $z=0$, we find that $75-80\,$% reside in the stellar halo of our simulated Milky Way analogue, while the remainder are gravitationally bound to lower-mass systems, including satellite halos.

We present a new method for generating emission spectra from polycyclic aromatic hydrocarbons (PAHs) in arbitrary radiation fields. We utilize the single-photon limit for PAH heating and emission to treat individual photon absorptions as independent events. This allows the construction of a set of single-photon emission "basis spectra" that can be scaled to produce an output emission spectrum given any input heating spectrum. We find that this method produces agreement with PAH emission spectra computed accounting for multi-photon effects to within $\simeq10\%$ in the $3-20~{\rm \mu m}$ wavelength range for radiation fields with intensity $U<100$. We use this framework to explore the dependence of PAH band ratios on the radiation field spectrum across grain sizes, finding in particular a strong dependence of the 3.3 to $11.2~\mu$m band ratio on radiation field hardness. A Python-based tool and a set of basis spectra that can be used to generate these emission spectra are made publicly available.

We present pkdgrav3, a high-performance, fully parallel tree-SPH code designed for large-scale hydrodynamic simulations including self-gravity. Building upon the long development history of pkdgrav, the code combines an efficient hierarchical tree algorithm for gravity and neighbor finding with a modern implementation of Smoothed Particle Hydrodynamics (SPH) optimized for massively parallel hybrid CPU/GPU architectures. Its hybrid shared/distributed memory model, combined with an asynchronous communication scheme, allows pkdgrav3 to scale efficiently to thousands of CPU cores and GPUs. We validate the numerical accuracy of pkdgrav3 using a suite of standard tests, demonstrating excellent agreement with analytic or reference solutions. The code was already used in several peer-reviewed publications to model planetary-scale impacts, where SPH's Lagrangian nature allows accurate tracking of material origin and thermodynamic evolution. These examples highlight pkdgrav3's robustness and efficiency in simulating highly dynamical, self-gravitating systems. pkdgrav3 thus provides a powerful, flexible, and scalable platform for astrophysical and planetary applications, capable of exploiting the full potential of modern heterogeneous high-performance computing systems.

We investigate the origin of unidentified, extended TeV source 1LHAASO J0500$+$4454, considering three possible origins: cosmic rays interacting with a molecular cloud (MC), particles accelerated in a currently undetected supernova remnant (SNR), and an energetic outflow powered by a pulsar. Upper limits on the CO and X-ray emission from the $\gamma$-ray emitting region disfavor the MC and SNR scenarios, respectively. If a nebula of inverse Compton scattering $e^{\pm}$ powers 1LHAASO J0500$+$4454, then SED modeling indicates that the current particle energy in the nebula is $\sim 4 \times 10^{48}$ erg. If the coincident magnetar SGR 0501$+$4516's rotational energy powered 1LHAASO J0500$+$4454, then a conservative energy budget calculation requires an initial magnetar spin period $P_{0} \lesssim 5$ ms and a spin-down timescale $\tau_{\rm sd} \lesssim 30$ yr, which has implications for the origins of magnetars.

Bulk meteoritic data show isotopic variability of $slow$-neutron-capture ($s$-process) origin in a several elements heavier than Fe. One peculiar feature is that the lighter $s$-process elements (e.g., Zr and Mo) present larger anomalies than the heavier $s$-process elements (e.g., Nd and W). To address this observation, we compared Zr and Nd data to model predictions of the s-process abundances at the surface of low-mass asymptotic giant branch (AGB) stars of initial metallicity from solar to twice solar. We found that the relative magnitude of the isotopic variability between these two elements can be matched by models of AGB stars of super-solar metallicity. The match is favoured by stronger convective overshoot, leading to a deeper dredge-up of the H-rich envelope into the He-rich region, and/or a smaller (~ half than standard) mass of the region rich in the $^{13}$C nuclei that produce free neutrons via the $^{13}$C($\alpha$,n)$^{16}$O reaction. We conclude that nucleosynthesis in AGB stars can match the difference in the magnitude of the bulk meteoritic variations in Zr and Nd, provided that super-solar metallicity stars are the original site of these signatures. The AGB stars that produced such variations could have belonged to the current population of old, super-solar metallicity stars seen in the galactic solar neighbourhood.

We investigate how galactic disk structures connect to the detailed properties of their host dark-matter halos using the TNG50 simulation. From the hydrodynamic and matched dark-matter-only runs, we measure a comprehensive list of halo properties describing density structure, angular momentum, shape, assembly history, and environment. Using the morphological decomposition developed in Paper I, we quantify the sizes, scale heights, and mass fractions of the disk components for galaxies at $0 \le z \le 4$. Random Forest (RF) regression shows that halo properties alone predict disk size and thickness with high accuracy, while Symbolic Regression (SR) provides compact empirical relations with slightly lower accuracy. Disk height is consistently easier to predict than disk size, and lower-mass halos yield higher accuracy than massive halos. Predictions based on halo properties measured in the hydro simulations outperform those based on halos matched in the dark-matter-only simulation, reflecting the imprint of baryonic restructuring on the inner halo. SHAP analysis reveals the most informative halo parameters include concentration, Einasto shape, global and inner spin, and recent mass accretion, though their importance varies across disk properties. We show correlations between disk size and the density-profile shape arise primarily from disk-induced modification of the inner halo, rather than a primordial connection. Finally, we point out that disks become more extended with respect to their host halos at higher redshift in low-mass halos, while massive high-redshift halos show the opposite trend. We provide SR-based prescriptions that accurately map halo properties to disk structures, offering practical tools for galaxy-halo modeling.

The exteriors of stellar and galactic dynamos are usually modeled as a current-free potential field. A more realistic description might be that of a force-free magnetic field. Here, we suggest that, in the absence of outflows, neither of those reflect the actual behavior when the magnetic field spreads diffusively into a more poorly conducting turbulent exterior outside dynamo. In particular, we explain why the usual ordering of the dipole magnetic field being the most slowly decaying one is altered, and that the quadrupole can develop a toroidal component that decays even more slowly with radial distance. This is a robust feature that persists even for spatially nonuniform magnetic diffusivities. It is best seen for spherical dynamo volumes and becomes more complicated for oblate ones. In either case, however, those fields are confined within a magnetosphere beyond which the field drops exponentially. We demonstrate that the Faraday displacement current, which plays a role in a vacuum, can safely be neglected in all cases. The superposition of magnetic fields from galaxies in the outskirts of the voids between galaxy clusters can therefore not explain the void magnetization of the intergalactic medium, reinforcing the conventional expectation that those fields are of primordial origin. For quadrupolar configurations, the synchrotron emission from the magnetosphere is found to be constant along concentric rings. The dipolar and quadrupolar configurations display large-scale radial trends that are potentially distinguishable with existing radio telescopes.

We present the results of customised single-pulsar noise analysis of 27 millisecond pulsars from the second data release of the Indian Pulsar Timing Array (InPTA-DR2). We model various stochastic noise sources present in the dataset using stationary Gaussian processes and estimate the noise budget of the InPTA-DR2 using Bayesian inference, involving model selection, Fourier harmonics selection, and parameter estimation for each pulsar. We check the efficacy of our noise characterisation by performing the Anderson-Darling test for Gaussianity on the noise-subtracted residuals. We find that all 11 pulsars with time baseline $\lesssim2.5\,\text{yr}$ show Gaussian residuals and do not have evidence for any red noise process in the optimal model, except for PSR J1944$+$0907, which shows presence of DM noise. PSRs J0437$-$4715, J1909$-$3744 and J1939$+$2134 show preference for the most complicated noise model, having achromatic and chromatic red noise processes. Only 4 out of 15 pulsars with time baseline $\gtrsim2.5\,\text{yr}$ show significant non-Gaussianity in noise-subtracted residuals. We suspect that this may require more advanced methods to model noise processes properly. A comparative study of six pulsars with data removed near solar conjunctions showed deviations from the parameter estimates obtained with the original dataset, indicating potential bias in red noise processes due to unmodeled solar-wind effects. The results presented in this work remain broadly consistent with the InPTA-DR1 noise budget, with better constraints obtained on noise processes for several pulsars and support for achromatic red noise in PSR J1012$+$5307 due to the extended time baseline.

We explore the evolution of sub-Neptune (radii between $\sim$1.5 and 4 R$_\oplus$) exoplanet interior structures using our upgraded evolution code, \texttt{APPLE}, which self-consistently couples the thermal and compositional evolution of the whole structure. We incorporate stably stratified regions with convective mixing and, for the first time, ab initio results on the phase separation of silicate-hydrogen mixtures to model silicate rain in sub-Neptune envelopes. We demonstrate that inefficient mantle cooling can retain sufficient heat to Gyr ages: inefficient heat transport from mantle to envelope alone keeps radii $\sim$10\% larger than predicted by adiabatic models at late times. Silicate rain can contribute an additional $\sim$5\% to the radius, depending on envelope mass and initial metal abundance. The silicate-hydrogen immiscibility region may lie in the middle or even upper envelope, far above the envelope-mantle boundary layer, and bifurcates the envelope into two an upper, hydrogen-rich region and a lower, metal-rich region above the mantle. If silicate rain occurs, atmospheres should appear depleted of silicates while radii remain inflated at late ages. To demonstrate the effects of inefficient mantle cooling, we present interior evolution models for GJ 1214 b, K2-18 b, TOI-270 d, and TOI-1801 b, showing that hot, liquid silicate mantles with thin envelopes reproduce their radii and mean densities, providing an alternative to water-world interpretations. These results imply that bulk compositions inferred from mean density must account for the mantle thermal state and the envelope mixing/phase-separation history; such thermal ``memories'' may constrain formation entropies and temperatures when metallicities are more precisely measured.

To detect ultra-high-energy neutrinos, experiments such as the Askaryan Radio Array and the Radio Neutrino Observatory in Greenland target the radio emission induced by these particles as they cascade in the ice. This is done by, amongst others, using deep in-ice antennas at the South Pole or in Greenland. A crucial step toward this goal is the characterization of the in-ice radio emission from cosmic-ray-induced particle showers. These showers form a primary background for neutrino searches, but can also be used to validate the detection principle and provide calibration signals for in-ice radio detectors. In this work, we use the Monte-Carlo framework FAERIE to perform the first characterization of cosmic ray signals with simulations that incorporate both their in-air and in-ice emissions. We investigate cosmic ray signatures such as their radiation energy, timing, polarization and frequency spectrum and quantify how they depend on shower properties. These results provide key guidelines for cosmic-ray identification and cosmic-ray/neutrino discrimination in future in-ice radio experiments.

We study the local evolution of the Universe around galaxy formation sites in the EAGLE50 large-volume reference simulation. Using the reduced inertia tensor (r-TOI), we followed the anisotropic evolution of initially spherical Lagrangian volumes (LVs) centred at galaxy formation sites, both in dark matter (DM) and in cold baryons (CB), from very high redshift $z=15$ onward. We describe LV deformation in terms of the r-TOI eigen-directions, principal axes, their derived shape parameters, and the timescales for the freezing-out of these principal directions and axes. Of particular interest are the age of the Universe, $t_{\rm U}$, when the local Cosmic Web (CW) spine emerges, and that when anisotropic DM mass arrangements (i.e., migrant mass flows) cease. We find that the shapes LVs acquire along their evolution affect the halo and stellar mass of their central galaxy: prolate-shaped LVs show a tendency to host low-mass galaxies at $z=0$, while massive galaxies tend to form within triaxial or oblate LVs. Also, the local CW spine tends to set in earlier on in LVs that are to host massive galaxies than in those harbouring less massive galaxies. In addition, anisotropic DM-mass rearrangements stop late on average, at $t_{\rm U}\sim 10.5\,$Gyr, and even slightly later for CB. Interestingly, $z=0$ LVs with either flattened configurations in CB or those that are highly prolate in DM, are more likely to host rotation-dominated galaxies. This effect increases from $z=1$ to $z=0$. Finally, the CB spine of LVs that are more likely to host rotation-dominated galaxies emerges at later times.

The Galactic Be star binary MWC 656 was long considered the only known Be star+black hole (BH) system, making it a critical benchmark for models of massive binary evolution and for the expected X-ray emission of Be+BH binaries. However, recent dynamical measurements cast doubt on the presence of a BH companion. We present new multi-epoch ultraviolet spectroscopy from the Hubble Space Telescope, combined with high-resolution optical spectra, to reassess the nature of the companion. The far-ultraviolet spectra reveal high-ionisation features -- including prominent N v and He ii lines -- which are absent in the spectra of normal Be stars and are indicative of a hot, luminous companion. Spectral modelling shows that these features cannot originate from the Be star or from an accretion disc around a compact object. Instead, we find that the data are best explained by a hot ($T_\mathrm{eff} \approx 85$ kK), compact, hydrogen-deficient star with strong wind signatures, consistent with an intermediate-mass stripped star. Our revised orbital solution and composite spectroscopic modelling yield a companion mass of $M_2 = 1.48^{+0.55}_{-0.46}\,\mathrm{M}_\odot$, definitively ruling out a BH and disfavouring a white dwarf. MWC 656 thus joins the growing class of Be+stripped star binaries. The system's unusual properties -- including a high companion temperature and wind strength -- extend the known parameter space of such binaries. The continued absence of confirmed OBe+BH binaries in the Galaxy highlights a growing tension with population synthesis models.

We investigate the non-linear evolution of matter perturbations in $f(R)$ models with the Chameleon screening mechanism. The novel feature of our investigation is an iterative solution for the non-linear equation for the scalar field $\chi = \Phi - \Psi$, where $\Phi$ and $\Psi$ are the potentials that characterise scalar perturbations of the metric. We demonstrate the scheme on spherical perturbations - smooth, compensated top-hats of varying length scales. We find that the effect of the Chameleon mechanism is seen most prominently on scales where the size of the top-hat is comparable to the Compton scale of the background. There is a density enhancement near the outer edge of the top-hat and the top-hat does not retain its shape. We explain this well-known observation in the context of the spatio-temporal evolution of the Compton scale. Additionally, we find a slight enhancement of the density near the origin, a feature not reported previously in the literature. On scales much smaller or much larger than the background Compton length, including the Chameleon screening has no appreciable effect on the perturbations. In the former, the growth is enhanced as compared to GR and is almost the same as GR in the latter. Finally, we examine the non-linear density velocity divergence (DVDR) relation and find that for evolution affected by Chameleon screening, the DVDR is no longer one-to-one even for a single profile. The relation between density and velocity depends on the location within the perturbation.

Cosmological inference from weak lensing (WL) surveys is increasingly limited by uncertainties in baryonic physics, which suppress the non-linear matter power spectrum on small scales. Multi-probe analyses that incorporate complementary tracers of the gas distribution around haloes offer a pathway to calibrate these effects and recover unbiased cosmological information. In this work, we forecast the constraining power of a joint analysis combining fiducial data from a Stage-IV WL survey with measurements of the dispersion measure from fast radio bursts (FRBs). We evaluate the ability of this approach to simultaneously constrain cosmological parameters and the astrophysical processes governing baryonic feedback, and we quantify the impact of key FRB systematics, including redshift uncertainties and source clustering. We find that, even after accounting for these effects, a 3$\times$2-point analysis of WL and FRBs significantly improves cosmological constraints, reducing the degradation factor on $S_8$ by $\sim 80\%$ compared to WL alone. We further show that FRBs alone are sensitive only to a degenerate combination of the key baryonic parameters, $\log_{10} M_{\rm c}$ and $\eta_{\rm b}$, and that the inclusion of WL measurements breaks this degeneracy. Finally, we extend our framework to incorporate galaxy clustering measurements using Luminous Red Galaxy and Emission Line Galaxy samples, performing a unified 6$\times$2-point analysis of WL, dispersion measures of FRBs, and galaxy clustering. While this combined approach tightens constraints on $\Omega_{\rm m}$ and $\log_{10} M_{\rm c}$, it does not lead to a significant improvement in $S_8$ constraints beyond those obtained from WL and FRBs alone.

In galaxy clusters, hierarchical assembly predicts the formation of stellar substructures and intracluster light (ICL), a diffuse stellar component tracing the global cluster potential. Because these features are extremely faint, alternative tracers such as globular clusters (GCs) provide a powerful tool to study cluster assembly. We use deep VLT/FORS $V$- and $I$-band imaging to investigate the GC population in the nearby Hydra I galaxy cluster ($\sim 45.7$ Mpc). GC candidates were selected from the $VI$ colour-magnitude diagram and divided into blue and red subpopulations. We find a clear spatial dichotomy: red GCs are concentrated around the massive central galaxies NGC 3311 and NGC 3309, while blue GCs are more extended and offset from the centre, coinciding with a secondary peak of X-ray-emitting gas. In the central regions, GC spatial distributions further depend on stellar population properties: young metal-rich GCs are more extended and may be linked to ram-pressure stripping, whereas old metal-poor GCs are more centrally concentrated, possibly originating from disrupted dwarf galaxies. Comparing the GC number density profiles to the surface brightness profile of NGC 3311, we find that the red GCs closely follow the galaxy light, while the blue population significantly deviates from it and traces the global gravitational potential of the cluster. This is also reflected in the specific frequency of blue GCs, which is approximately $\sim 5\times$ higher in the ICL-dominated outskirts than in the inner regions dominated by red GCs. Finally, we present a novel method to constrain the evolution of the galaxy luminosity function of the cluster using GC specific frequencies and colour distributions, yielding a past faint-end slope of $\alpha=-1.81^{+0.16}_{-0.16}$ compared to $\alpha=-1.41^{+0.08}_{-0.05}$ today, consistent with high-redshift observations and cosmological simulations.

NGC 4051 is a nearby (16.7 Mpc), Narrow Line Seyfert 1 galaxy (NLS1), which has a low black hole mass of $10^6$ M$_{\odot}$. It is also known for its rapid X-ray variability, on timescales of kilo-seconds and has a complex, multi component wind in both the soft X-ray and Fe K bands. Here we present the first high resolution XRISM Resolve spectrum of NGC 4051, which was captured in a historically bright state for a 150 ks exposure. XRISM resolves two blue-shifted Fe K shell absorption troughs in the mean spectrum, which can be ascribed to H-like iron and arises from two outflow components with outflow velocities of 0.025c and 0.04c. A time dependent spectral analysis shows that the iron K absorption is variable on timescales of less than a day, increasing in velocity over the duration of the observation. The velocity changes may be explained either by the passage of two separate transiting absorbers, of different velocities, or by a single accelerating outflow of approximately constant column density. In the latter case, the wind acceleration is likely to be too large to be caused by radiation pressure and instead magnetic driving is favored to accelerate the wind up to 0.04c. The outflow can originate from an accretion disk wind, whose kinetic power is sub-Eddington in contrast to recent examples of winds from powerful, luminous quasars observed by XRISM.

Analyses of baryon acoustic oscillations (BAO) commonly employ template-based methods to extract compressed parameters from the clustering of dark-matter tracers, which are then interpreted in terms of ratios of the sound-horizon scale and cosmological distances relative to a fiducial cosmology. A small mismatch between the sound-horizon scale derived from the standard analytic formulation (integral over the sound speed) and the effective scale imprinted in clustered matter can, however, introduce a systematic bias in cosmological inference. We extend previous work to a broader class of cosmological models, quantify this bias for surveys with DESI-like precision, and propose strategies to correct for the effect. We find that the induced bias becomes a significant fraction of the statistical uncertainty for deviations from the fiducial cosmology, at the level of $|\Delta \Omega_m| = 0.03$ and $|\Delta N_\mathrm{eff}| = 0.3$, and for very precise data corresponding to a forecasted Year-5 DESI survey (or other stage IV dark energy galaxy surveys). We present several ways to correct for this effect, suitable for a variety of applications. We therefore recommend that analyses exploring such parameter regimes either apply the proposed corrections or include an appropriate systematic error budget.

SDSS~J143016.05+230344.4 ($z=0.08105$) has been proposed as a candidate pre-coalescence supermassive black hole binary and shows remarkable multiwavelength variability. Its radio evolution provides a direct probe of the compact emitting region and of the physical origin of the late-time activity. We aim to localize the variable radio emission, characterize its spectral evolution, and constrain whether the radio brightening is produced by a newly emerging compact component, external absorption, or dissipation in a structured circumnuclear environment. At all epochs, the radio emission is dominated by a single unresolved milliarcsecond core with $T_{\rm B} \gtrsim 10^{7}$ K, constraining the variable emission to $\lesssim 0.3$ pc. The broadband spectra require two synchrotron self-absorbed components: a persistent low-frequency component with $\nu_{\rm p,steady} \approx 0.74$ GHz and $S_{\rm p,steady} \approx 1.22$ mJy, and a flare component whose turnover evolves from $(6.35 {\rm GHz}, 0.18 {\rm mJy})$ in 2022 February-May to $(8.61 {\rm GHz}, 0.38 {\rm mJy})$ in 2022 December, and then to $(5.83 {\rm GHz}, 0.25 {\rm mJy})$ in 2023 March-April. The flare contribution at 15 GHz reaches $\sim 80\%$ and matches the near-epoch VLBI recovery fraction, showing that the high-frequency brightening arises from a newly formed compact synchrotron component. A second brightening of the 15.2 GHz VLBI core is detected between 2023 September and 2024 February, while the source remains unresolved. Equipartition scalings imply characteristic radii of $\sim 5 \times 10^{-4}$ pc for the flare and $\sim 9 \times 10^{-3}$ pc for the steady component, and indicate a steep inner circumnuclear density profile, $n \propto R^{-1.7}$. The delayed radio flare is best explained by dissipation in an outflow or jet-base disturbance propagating through a structured circumnuclear medium.

Measurements of the kinematic cosmic dipole continue to show an intriguing tension between the value inferred from the CMB and that obtained from high-redshift source number counts. While the measured dipole direction appears consistent, the amplitude, set by the observer's peculiar velocity $v_{o}$, remains in significant disagreement. In this paper, we propose using strong gravitational lenses with well-measured Einstein radii to estimate the kinematic cosmic dipole, through the relativistic aberration of the Einstein angle induced by the observer's motion. We show that this effect could be detected solely from measurements of the Einstein radius in wide, high-resolution imaging surveys such as Euclid. However, the precision achievable using Einstein-radius measurements alone, without redshift or lens-galaxy mass information, appears insufficient to discriminate between the CMB value of $v_{o}$ and that derived from source number counts at high statistical significance. Nevertheless, we demonstrate that including a large sample of lenses with available kinematic information, either via the Fundamental Plane relation or, ideally, through spectroscopic velocity-dispersion measurements, drastically reduces the noise and substantially improves the constraining power of this method. We show that, for a realistic sample of strong lenses detected by Euclid and complemented with spectroscopic velocity dispersion measurements from 4MOST or DESI, it is possible to discriminate between the CMB- and source-number-counts-inferred values at the $\sim 5\sigma$ level using a new, fully independent method. We further demonstrate that this technique is only weakly sensitive to strong-lensing selection effects, with selection biases and threshold effects estimated to be well below the 1% level.

Very low-luminosity objects in nearby star-forming regions have been identified as promising proto-brown dwarf candidates. The study of their multiplicity can shed light on the dominant formation mechanism of these substellar objects. We aim at studying the multiplicity of the very low luminosity object IRAM 04191+1522. To do so, we have obtained 0.89mm ALMA observations with a very extended configuration, achieving an angular resolution of ~0.04 arcsec (6 au at 140 pc). We have complemented our data with new VLA observations, and ALMA archival data at 1.3mm. As a result, we resolve IRAM04191+1522 into a close binary candidate for the first time. The binary is detected in the ALMA continuum data with a projected separation of ~80 mas, or 11 au at a distance of 140 pc. The two sources are oriented in the East-West direction, with the eastern component being brighter and more extended than the western one, which is marginally resolved. The analysis of C18O(2-1) archival data reveals gaseous material in rotation around the binary, presumably from a circumbinary disk with ~27 au of radius centered on the faintest ALMA component. A fit of the position-velocity diagram allows us to estimate a total dynamical mass for the system of 50+-40 MJup. Therefore, we classify IRAM04191 as a tight proto-brown dwarf binary candidate. The VLA data reveals the detection of a single object closer to the western ALMA source, and with a spectral index consistent with a radio jet.

Interstellar comets arrive with key ingredients for deep-space platforms already in place: volatile inventories convertible to propellant, natural rotation providing continuous attitude variation, and hyperbolic trajectories that carry them through the inner Solar System and back out to interstellar space. Rather than constructing spacecraft from scratch, we ask what \emph{minimal modification} is required to steer such a body along a controlled trajectory. The answer is surprisingly modest. By relaxing full six-degree-of-freedom control to forward-cone steering -- sufficient for practical navigation -- we show that \emph{four thrusters suffice}: one primary jet and three secondary jets at $120^\circ$ intervals. The secondary jets synthesize continuous in-plane steering, while the primary jet provides low-bandwidth attitude shaping: as the body rotates, the primary-jet torque direction sweeps predictably over a cycle, enabling out-of-plane steering via phase-scheduled firing. We formalize reachability under bounded-curvature constraints, characterize the rotation-mediated steering envelope, discuss enabling requirements including non-solar power at large heliocentric distances, and identify operational regimes and observable signatures implied by active trajectory control. The setting of a nutating axis is briefly considered and conjectured to preserve core results. The findings contribute to the broader effort of understanding the dynamics and control of small-body missions and offer a reference architecture relevant to long-horizon deep-space exploration and to potential planetary-defense concepts.

Interstellar objects (ISOs) motivate a coupled mission-design and inference question relevant to spacecraft dynamics and control in extreme environments: if volatile-rich, rotating comet-like bodies were used for sustained deep-space navigation by exploiting pre-existing hyperbolic motion and in-situ propellant, what stability requirements arise under non-gravitational forcing, and what astrometric signatures might distinguish active stabilization from uncontrolled natural dynamics? We develop a stability-theoretic framework for trajectory tracking with jet-actuated correction, and show that high-speed transit geometry -- including debris-belt avoidance and encounter phasing -- tightly constrains feasible trajectories, making long-horizon tracking stability mission-critical. We model tracking residuals as the balance of disturbances and corrective action, and derive stability conditions across four levels: disturbance-energy stability, outer-loop contraction, actuator-memory stability, and rotation-mediated (Floquet) stability. The analysis implies residual diagnostics that can motivate empirical tests: under comparable forcing, effective stabilization is expected to strengthen short-horizon error correction, reduce event-conditioned persistence and variance clustering, regularize standardized innovations, and yield bounded post-shock recovery. More broadly, the framework provides a reference for deep-space guidance and control under nonlinear, multi-field disturbances and for planetary-defense concepts involving attitude shaping or impulsive kinetic impact.

We discuss in some details a novel algorithm for performing partial-sky spherical harmonic transforms (SHT), building on the Fourier-sphere method of Reinecke et al (2023) handling efficiently high numbers of arbitrary locations on the sphere. Our main motivations are Cosmic Microwave Background lensing from the South Pole Telescope, and the South Pole Observatory program targeting primordial gravitational waves from inflation, requiring high-resolution, numerically intensive work on small sky fractions. We achieve speed-up factors ranging from 3 to 10 on SPT-3G main field and BICEP3 deep footprint, and much more on smaller patches. More generally, the algorithm eliminates in our case study the usual disadvantages of arbitrary pixelisations in comparison to isolatitude pixelisations or flat-sky approximations, making it ideal for ambitious workflows that require repeated SHTs on limited sky regions.

We present a new careful and comprehensive analysis the observations of the satellites of Jupiter from the Sidereus Nuncius that extends and complements previous similar studies. Each observation is compared to the predictions obtained using a modern sky simulator, verifying and trying to understand them individually. The work considers both the information that can be extracted from the sketches and the angular measurements reported by Galilei. Angular measurements allow assessing the absolute accuracy in relation to modern ephemerides. We evaluate the performances of the telescope in terms of separation power of close-by satellites and the inefficiency in the detection connected to the proximity to the disk. A sinusoidal fit of the data, allows measuring the relative major semi-axes of the satellites' orbits and their periods with a statistical precision of 2-4\% and 0.1-0.3\% respectively. The posterior fit error is used to estimate the measurements precision. We show that with this data one can infer in a convincing way the third law of Kepler for the Jupiter system. The 1:2:4 orbital resonance between the periods of Io and Europa/Ganymede can be determined with \% precision. In order to obtain these results it is important to separate the four datasets. This operation was an extremely difficult task for Galilei. Nevertheless we show how some indication on the periods emerge from the using the modern Lomb-Scargle technique on the full data-set. We briefly extend the use of the simulator to verify the accuracy in the seven observations of the Moon and the performance in reproducing the Pleiades, the Orion belt, the Orion head and the Beehive cluster. Finally we present images obtained with a replica of the telescope that highlights the challenges of these observations thus confirming the excellence underlying this amazing set of early scientific data.

We consider the formation of metastable cosmic strings in a left-right unified theory. The produced monopoles are diluted by a stage of F-term hybrid inflation (FHI) which is realized consistently with the SUSY breaking and a global $U(1)$ $R$ symmetry in the context of a $U(1)_{\rm R}\times U(1)_{B-L}$ extension of MSSM. The hidden-sector Kaehler manifold enjoys an enhanced $SU(1,1)/U(1)$ symmetry with the scalar curvature determined by the achievement of a SUSY-breaking de Sitter vacuum without ugly tuning. FHI turns out to be compatible with data -- including the recent ACT results --, provided that the magnitude of the emergent soft tadpole term is confined in the range $(0.1-70)$ TeV, and it is accompanied with the production of cosmic strings. Their dimensionless tension $G\mu_{\rm cs}\simeq(1-11)\cdot10^{-8}$ interprets the present observations from PTA experiments on the stochastic background of gravitational waves. The $\mu$ parameter of MSSM arises by appropriately adapting the Giudice-Masiero mechanism and facilitates the out-of-equilibrium decay of the $R$ saxion at a reheat temperature lower than about $34$ GeV. The SUSY mass scale turns out to lie in the PeV region.

We study magnetic conversion of ultra-relativistic axion-like particles (ALPs) into photons in compact-star environments, focusing on the hot, transient conditions of core-collapse supernova (SN) remnants and neutron-star mergers (NSMs). We address previously overlooked uncertainties, particularly the suppression caused by ejected matter near the stellar surface, a region crucial to the conversion process. We derive analytical expressions for the transition rate; they reveal the influence of key parameters and their uncertainties. We update constraints using historical gamma-ray data from SN~1987A and find $g_{a\gamma}<5\times10^{-12}~{\rm GeV}^{-1}$ for $m_a\lesssim10^{-9}$ eV. We also forecast sensitivities for a future Galactic SN and for NSMs, assuming observations with Fermi-LAT or similar gamma-ray instruments. We distinguish ALPs -- defined as coupling only to photons and produced via Primakoff scattering -- from axions, which also couple to nucleons and emerge through nuclear bremsstrahlung. We omit pionic axion production due to its large uncertainties and inconsistencies, though it could contribute comparably to bremsstrahlung under optimistic assumptions. For the compact sources, we adopt time-averaged one-zone models, guided by numerical simulations, to enable clear and reproducible parametric studies.

We present a quantum quadratic gravity inflationary scenario that can accommodate the new cosmological constraints, which have disfavored Starobinsky inflation. The theory is asymptotically free in the ultraviolet, but 1-loop running is found to dynamically lead to slow-roll inflation toward the infrared. When a large number of matter fields contribute to the beta functions, the spectral index and the tensor-to-scalar ratio can be phenomenologically viable. We find that as inflation ends, the theory approaches its strong coupling regime and general relativity must emerge, as an effective field theory, as the universe must reheat and enter its standard radiation era. In order to avoid strong coupling, a minimum tensor-to-scalar ratio of 0.01 is predicted for this theory. Our framework offers a laboratory for connecting a concrete ultraviolet completion (quantum quadratic gravity) with inflationary dynamics, reheating, and precise cosmological observations.

The electroweak sphaleron rate in the high temperature phase of the Standard Model is inversely proportional to the weak-isospin conductivity. So far, only electroweak interactions were included in its computation. Here we take into account quark scattering through strong interactions at leading-log order. These reduce the quark contribution to the conductivity by up to 15 %, and the total conductivity by up to 6 %.

The non-conservation of the energy-momentum tensor in $f(R,T)$ gravity can be interpreted as an effective manifestation of dissipation. Motivated by this, we propose a new formulation of $f(R,T)$ gravity based on the Herglotz variational principle, which extends the usual {Hamilton} variational principle to dissipative systems by allowing the Lagrangian to depend explicitly on the action. The resulting gravitational field equations extend those of $f(R,T)$ gravity by including Herglotz contributions. In the Newtonian limit, these contributions modify the gravitational potential, allowing us to constrain the Herglotz vector through Mercury's perihelion precession and the relativistic light deflection. The Herglotz corrections lead to a scaling law consistent with observations from the Cassini spacecraft. Examining two representative cosmological models, the Herglotz vector effectively reduces to a single function that, under suitable conditions, can play the role of a cosmological constant, providing an alternative mechanism for the Universe's accelerated expansion. Within the Herglotz variational approach, the linear $f(R,T)=R+\alpha T$ model, previously ruled out in the standard formulation due to its fixed deceleration parameter, becomes consistent with observations {for a certain choice of the Herglotz vector.

We present new spherically symmetric solutions of the Einstein-Klein-Gordon equations in a quasi-stationary approximation that describe self-gravitating scalar field configurations around a black hole, including angular momentum number $\ell$. An approach analogous to the one which gives rise to $\ell$-boson stars is used here to construct self-gravitating ``gravitational atoms" with $\ell\ge0$. We refer to these new solutions as {\it noble gravitational atoms}, by analogy with noble atoms, which are characterized by closed electron shells. We show that, in the proper limit, noble gravitational atoms approach $\ell$-boson stars globally, displaying noticeable differences only in a region very close to the event horizon. Noble gravitational atoms with $\ell>0$ sometimes present density maxima located at relatively large radii, with small density close to the horizon for $\ell>1$. Furthermore, they do not always present the typical density spike at the event horizon if $\ell > 0$; on the contrary, they sometimes exhibit a small dip there. When $\ell=0$, a spike can appear, but its contribution to the total mass density is always negligible. The size, density, and lifetime of these objects vary significantly depending on the parameters, being in some cases as large as galaxies, as dilute as dark matter, and as long-lived as the Universe itself.