We study the effect of temperature on secular, compositional $g$-modes in the core of inviscid neutron stars. Using a chiral $SU(2)_f$ sigma model, we construct isentropic temperature profiles for hot and dense matter and find that the frequency of the global core $g$-mode's dependence on temperature is governed by the nuclear symmetry energy slope parameter $L$. As a result, the $g$-mode frequency of a warm neutron star can be either higher or lower than that of its cold counterpart, depending on $L$. Our results highlight the interplay of thermal effects and composition gradients, and demonstrate the potential of neutron star $g$-mode observations to constrain the density dependence of the symmetry energy.

The advent of high-fidelity imaging of supermassive black holes calls for efficient and robust data-analysis methods. In this work, we use $\texttt{Jipole}$, a differentiable, $\texttt{ipole}$-based radiative transfer code, to enable gradient-based analyses of images generated from state-of-the-art general relativistic magnetohydrodynamic (GRMHD) simulations. We compute image sensitivities, i.e., pixel-wise derivatives of the intensity with respect to model parameters, which form the Jacobian of the forward model and define a local map from parameter space to image space. Using these sensitivities in a mock data analysis, we find that GRMHD-based images generate a structured error landscape for parameter fitting, with anisotropies and local minima, making parameter exploration nontrivial but still tractable when guided by gradient information. We characterize this landscape through the Jacobian and assess the feasibility of gradient-based recovery under idealized, blurred, and noisy conditions. Our results show that automatic differentiation-computed image gradients can guide parameter exploration effectively even in the presence of noise. These findings establish a basis for efficient, high-precision model--data comparisons in black hole imaging and motivate the integration of these sensitivities into advanced inference frameworks.

We present JWST MIRI and NIRCam imaging of the inner ~5 kpc of the M82 superwind at 0.05-0.375'' (~0.9-6.5 pc) resolution. Targeted filters probe emission from polycyclic aromatic hydrocarbons (PAHs; F335M, F360M, F770W, F1130W) and continuum (F250M, F360M). The images reveal a network of cool wind filaments traced by PAHs. PAH surface brightness declines with the inverse square of distance to the midplane, suggesting that the incident radiation field from the starburst drives the observed PAH intensity out to 2.5 kpc. The 3.3/11.3 and 3.3/7.7 band ratios show uniformity with distance from the starburst, though comparisons with mid-IR dust emission models indicate a modest shift toward larger PAHs. Outside the disk, 11.3/7.7 increases moderately, reflecting that PAHs become more neutral with distance from the starburst as they are exposed to a declining radiation field and ionization parameter. Overall, PAHs in the wind are consistent with standard-to-large sizes and standard-to-high ionization states. Including Spitzer and Herschel data, PAH abundance (qPAH) is set at ~1% in the starburst and remains unchanging out to 5 kpc off the disk. This flat qPAH profile suggests that PAHs are shielded from the hot wind, perhaps residing in the surface layers of cool clouds, with possible replenishment from cloud interiors and enrichment of the halo from previous bursts. In this picture, clouds are not dense enough to promote PAH growth, and they likely undergo radiative cooling and mixing with the hot phase to survive the gauntlet for at least ~20 Myr.

We investigate the possibility of a binary supermassive black hole system at the center of MCG+11--11--032, a local (z = 0.036) Seyfert 2 galaxy. Prior work with stacked Swift/XRT spectra suggested the presence of two Fe K$\alpha$ lines (at 6.16 keV and 6.56 keV) with 2$\sigma$ confidence. This could be consistent with a prediction of several hydrodynamical models, in which each black hole hosts a mini-disk and contributes one Doppler-shifted Fe K$\alpha$ line to the total spectrum. Another study using a single exposure from Chandra/ACIS did not find evidence for a double line. Here, we conduct follow-up with two epochs of XMM-Newton/EPIC data spaced $\sim$6 months apart. After fitting our spectra with models from the previous two studies, we do not find evidence for a double iron line in either observation. Our best-fit model yields $\Gamma = 1.63^{+0.20}_{-0.21}$ and $N_\text{H}/10^{22} \text{ cm}^{-2} = 17.9^{+2.7}_{-2.4}$ for the first epoch, and $\Gamma = 1.46^{+0.22}_{-0.24}$ and $N_\text{H}/10^{22} \text{ cm}^{-2} = 17.1^{+2.7}_{-2.4}$ for the second. We compare our spectral parameters with those derived in past work on this source, finding broad agreement with prior datasets. Lastly, we discuss the properties of MCG+11--11--032 alongside samples of Seyfert 2 galaxies from the literature, finding that it is consistent with this population and the single AGN scenario.

Relativistic jets launched by active galactic nuclei are fundamental for understanding the physics of accreting supermassive black holes and their immediate environment, yet the mechanisms driving jet launching remain uncertain. In this study, we investigate the sub-parsec jet of 3C 84 using multi-epoch, multi-frequency, very long baseline interferometry (VLBI) observations with the European VLBI Network and the Very Long Baseline Array at 22 and 43 GHz. We analyse the evolution of the spectral index gradient in the core region to relate the observed structure to physical interpretations and to discriminate between competing jet launching models. Furthermore, we examine the impact of the ambient medium and magnetic field configuration on jet morphology and dynamics over time, and explore their connection to a coinciding $\gamma$-ray flare. Our spectral analysis reveals significant changes across three epochs, indicating dynamic activity between filamentary structures on sub-parsec scales, evolving magnetic fields, and a complex interaction with the surrounding medium, all of which shape the innermost jet and may influence its high-energy emission.

We present an Atacama Large Millimeter/submillimeter Array (ALMA) spectroscopic follow-up survey of an 870 $\mu$m-selected sample of dusty star-forming galaxies (DSFGs) in the GOODS-S field. We use these linescans to identify or confirm spectroscopic redshifts (spec-zs) for 20 sources. Including spec-zs from the literature, there are now secure or tentative spec-zs for 52 out of 75 DSFGs (69%). At $f_{870}>2.5$ mJy, the sample is 97% spectroscopically complete, allowing us to model the full DSFG redshift distribution down to nearly the confusion limit for a 15-m telescope at 850 $\mu$m. This is the highest completeness for an unbiased sample at this flux limit to date. We find that nearly all of the DSFGs in our sample that were targeted with JWST/NIRSpec were spectroscopically identified, without much dependence on near-infrared or submillimeter flux or redshift. However, only 29% of our sample have JWST spectroscopic coverage. We use the spec-zs to evaluate various photometric redshift (photo-z) estimates, finding that all methods exhibit an outlier fraction of at least $>20$%. Nearly all of the photo-z methods tend to overshoot the redshifts, leading to overestimates of the number of DSFGs at high redshift ($z>4$). Our results suggest that $\lesssim10$% of $f_{870} \gtrsim 2$ mJy DSFGs lie at $z>4$ and $\lesssim2$% lie at $z>5$, reflecting a steep decline in the abundance of massive dusty galaxies in the first 1.5 Gyr.

Galactic radio transients (GRTs) are mysterious short-lived (~days to months) radio transients that are quiet at all other wavelengths. Until now, roughly half a dozen such sources have been reported, predominantly towards the Galactic center. However, no unifying properties have been identified among these, leaving their nature, emission mechanism, and even classification poorly understood. Due to the lack of periodic and uniform radio observations over wide areas of the Galactic plane until now, the sample size of such transients remained limited. Here, we use radio observations from the Australian SKA Pathfinder's Variables and Slow Transients survey to discover six new radio transients along the Galactic plane that resemble GRTs. Detailed investigation of archival data suggests that these sources may be divided into two classes: sources that exhibit sporadic, pulse-like (minutes) radio emission, and sources that exhibit long-term (weeks) flaring-type radio emission. For the short-time variable sources, we draw similarities between optically bright long-period radio transients and our sample to propose wide-orbit (~days) white dwarf binaries as underlying sources. For sources that show long-term outbursts, we draw comparisons between dust-obscured outbursts from WD binaries and our sample. These results could imply that the ongoing wide-field radio surveys are uncovering radio emission from sub-populations of WD binaries that were previously unexplored.

The multiply imaged SN 2022riv was discovered through a search of galaxy cluster fields as part of a Hubble Space Telescope (HST) SNAP program to find highly magnified stars. The supernova (SN) was detected in the last-to-arrive image of a galaxy at redshift $z=1.522$ strongly lensed by the foreground galaxy cluster RX J2129.7+0005. Follow up James Webb Space Telescope (JWST) NIRSpec G140M and PRISM spectroscopy yields a Type Ia SN classification. Using the SALT3-NIR light-curve fitter, we obtain a cosmology-independent measurement of the magnification of $5.35\pm1.01$ for the last-to-arrive image of the SN, with multiple SALT SN spectral time-series models yielding consistent constraints. The last-to-arrive image of SN 2022riv we detect appeared adjacent to the brightest cluster galaxy (BCG) at a location with an exceptionally high stellar mass density ($\sim 1-2$ dex higher than that of SN Refsdal), where microlensing is expected to introduce a 20-50% modulation of the magnification. Analyzing six independent lens models of the cluster, we find that four predict the magnification with much greater precision ($p < 0.05$) than would be expected by random chance, given the large effect anticipated from microlensing. Five models yield magnifications of roughly $4-7$ (within $1\sigma$) prior to accounting for microlensing, whereas HoliGRALE favors a significantly higher value of $15.39 \pm 0.85$. After incorporating nominal microlensing, the HoliGRALE prediction is within $1\sigma$ tension with our measurement. A companion paper (Dalrymple et al.) will present constraints on the relative time delay of the image that arrived earlier.

Recent analyses combining DESI DR2 BAO with CMB and SNe Ia data have reported $2.8$--$4.2\sigma$ evidence for dynamical dark energy, but the significance depends strongly on the supernova sample, raising the question of whether this signal reflects new physics, dataset-specific systematics, or the choice of dark energy parameterization. We investigate this question by analyzing four SNe Ia compilations (Pantheon, Pantheon+, DES-Dovekie, and Union3) with DESI DR2 BAO and Planck CMB distance priors, using flux averaging, model-independent expansion rate extraction, parametric ($w_0 w_a$CDM) fits, and a non-parametric reconstruction of the dark energy density ratio $X(z) \equiv \rho_{\rm DE}(z)/\rho_{\rm DE}(0)$. Flux averaging reduces the $\Omega_m$ difference between SNe and DESI from ${\sim}2\sigma$ to ${\sim}1\sigma$ for Pantheon+ and DES-Dovekie. The reconstructed $X(z)$ for DESI DR2 + CMB + SNe is consistent with $\Lambda$CDM for Pantheon, Pantheon+, and DES-Dovekie except at $0.5

Stellar abundances reveal non-monotonic [Y/Eu] and [Sr/Eu] evolution, a systematic decline with increasing [Eu/H] at low metallicity, a minimum at $[\rm{Eu/H}] \sim -0.3$ and then a rise at high metallicity. This behavior requires at least three distinct neutron-capture sources operating on different timescales. We develop a one-zone chemical-evolution model constraining their typical delay-times, rates, and yield ratios. Reproducing the observed $\rm{[Y/Eu]}$ and $\rm{[Sr/Eu]}$ sequences requires, a delayed $r$-process channel (most likely binary neutron-star mergers) dominating Eu production ($\gtrsim 95\%$ of total Eu). A prompt channel preferentially producing first-peak elements with minimal Eu, explaining the increasing [Y/Eu] at decreasing [Eu/H] below $[\rm{Eu/H}] \lesssim -2.5$; and delayed AGB $s$-process enrichment with delays greater than $t_{min} = 0.3-0.6$\,Gyr reproducing the late-time upturn in Y (Sr). Our model quantitatively reproduces all constraints, including the large $\Delta[\rm{Y/Eu}] \approx 0.6$ dex variation between the late-time rise [Eu/H] and the minimum value, the location of the minimum at [Eu/H] $\sim -0.3$ and late-time rise. The first-to-second peak yield ratios correspond to $[\rm{Y/Eu}] \approx -0.3$ (prompt) and $\approx -0.8$ (BNS mergers). The observed $\Delta[\rm{Y/Eu}]$ amplitude establishes a model-independent lower limit on the first to second peak yield ratio $\gtrsim 3.4$ between the prompt and delayed channels, ruling out models with similar prompt and delayed yield ratios. These results demonstrate that explaining the observed heavy-element abundance patterns requires multiple channels with distinct nucleosynthetic signatures and operational timescales, providing constraints on the relative rates, delay times, and yield patterns of candidate production sites.

The lopsidedness of galaxies is a commonly observed phenomenon, and through different studies, it has been observed that nearly 30% of galaxies show this phenomenon. In this work, we study morphological lopsidedness in both stellar and gas disks in the inner and outer regions using Fourier analysis techniques and compare the results for a sample of nearby galaxies with different morphologies and environments. Although lopsidedness can result from diverse factors like tidal interactions, gas accretion, and internal instability, recent studies suggest it is a common feature that is not solely reliant on rare events, and moderate lopsidedness most likely results from the disk's response to a lopsided dark matter halo potential. Assuming lopsidedness originates due to a lopsided halo, we find the morphological and kinematic halo perturbation parameters in the same radial range. Unlike previous studies, we use 3D kinematic modelled rotation curves for finding kinematic lopsidedness and, hence, kinematic halo perturbation parameter. Although the detected linear correlation between them is not statistically significant for our small sample of eleven galaxies, this approach provides a more uniform and physically consistent framework to test the theoretically expected similarity between morphological and kinematic halo perturbation parameters. Further, within this framework, the discrepancy between them does not appear to depend on the nature of the rotation-curve asymmetry of the two sides of the galaxy, in contrast to trends seen in earlier studies. In future work, we plan to extend this analysis to a substantially larger sample in order to robustly assess these findings.

Polarized emission from dust and synchrotron radiation from the ISM are the dominant foregrounds for CMB polarization and are a major challenge for extracting the primordial signal on large angular scales. A key characteristic of the galactic foreground emission is its $EE/BB$ ratio. We argue that MHD shocks play an important role in setting the observed $EE/BB$ ratio. To support this, we first analyze quasi-linear magnetohydrodynamics (MHD) simulations to obtain an $EE/BB$ ratio that increases as $\sim k^2$, then show that with increasing energy injection rates, the $EE/BB$ ratio flattens to a value $\gtrsim 1$, approaching observational results. Looking at the distribution of the velocity divergence, a tail with power law $-7/2$ develops around the same injection rates where the $EE/BB$ ratio flattens. While the system becomes more isotropic, MHD shocks are intrinsically anisotropic and lead to the $E/B$ power asymmetry. We also observe total pressure balance among all our simulations, indicating slow wave dominance. Therefore, in the regime we consider, it is important to go beyond linear MHD equations to understand the foreground radiation.

Understanding how galaxies ignite and extinguish their star formation remains a cornerstone question in modern astrophysics. Recent JWST surveys have revealed an overabundance of massive quiescent galaxies in the first billion years of the Universe, challenging current models of galaxy evolution. In the nearby Universe, ram pressure stripping (RPS) is a major environmental mechanism capable of rapidly shutting down star formation, yet direct observation remains scarce at redshift $z\gtrsim1$, and its role at $z>2$ is even poorly constrained by simulations. Here, we utilize JWST and ALMA observations to present direct evidence of RPS in the post-starburst galaxy A2744-JF-z3, residing in a galaxy group at redshift 3.06, the earliest such detection to date. Spectroscopic diagnostics and spectral energy distribution modeling reveal the ongoing removal of cold gas and dust, coincident with the abrupt cessation of star formation. Contrary to hydrodynamical simulations that predict a reduced incidence of RPS at high redshift, our results instead imply that RPS can operate at $z>3$, suggesting a highly stochastic and impulsive stripping within a clumpy, filamentary intra-group and circumgalactic medium. These observations extend environmental quenching well into the epoch of galaxy assembly, highlighting RPS as a previously overlooked decisive pathway to rapid quenching in nascent groups and protoclusters in the early Universe.

Cosmic ray (CR) transport is usually modeled with a single diffusion coefficient, but this description captures only the growth of the variance and not the full transport process. Distinct transport mechanisms can share the same effective diffusion coefficient while producing different particle distributions and approaches to the diffusive limit. This limitation is especially relevant in realistic multiphase, structured, and time-dependent media, and is also reflected in observed environmental variations in CR transport near pulsar wind nebulae, supernova remnants, and molecular clouds. Particle-tracing studies also show clear departures from standard diffusion, including both superdiffusion and subdiffusion. We therefore develop a propagator-based framework centered on $P(x,t)$, the probability distribution of particle positions, or equivalently its Fourier-Laplace transform $P(k,s)$. This object is compact and statistically complete, and naturally exposes memory: the CR flux can depend on earlier gradients when unresolved trapping or phase changes are coarse-grained away. Using the Montroll-Weiss formalism, we show how to measure $P(k,s)$ directly from trajectories, how to recover the associated memory kernel, and how to represent broad kernels efficiently with a Prony expansion. Applied to a multiphase medium, the framework shows that slow regions can regulate escape without dominating the total residence-time budget. We also introduce an accelerated Monte Carlo method for coarse-grained transport, and show that if trapping structures evolve while particles are still sampling them, the static long-time limit need not be reached. This paper provides the foundation for future observational applications, particle-tracing measurements, and CR-MHD closures.

JWST has revealed an overabundance of very bright, blue galaxies at z>10, raising fundamental questions about how star formation and feedback operate at Cosmic Dawn. We present new JWST/NIRSpec MSA PRISM/CLEAR spectroscopy of JADES-GS-z14-0 (z=14.18) obtained with the JADES and OASIS programmes. While the rest-frame UV continuum flux level and shape are consistent between the two datasets, the OASIS spectrum shows a 10$\sigma$ detection of the CIII]$\lambda\lambda1907,1909$ emission line, with a luminosity three times higher than that measured in the JADES data. This difference is naturally explained by the offset in shutter placement between OASIS and JADES, implying that the CIII] emission is spatially displaced by $\sim400$ pc from the stellar continuum. The non-detection of CIII] in NIRCam medium-band imaging indicates that the emitting region is extended on scales $\gtrsim165$ pc, with a surface brightness below the detection threshold. Interpreting this diffuse, carbon-enriched gas as the result of ongoing or past outflows, we infer a mass outflow rate of $\dot{M}_{\rm out}\sim160~{\rm M_\odot\,yr^{-1}}$. We compare it with the star-formation rate (SFR) and derive a mass-loading factor of $\eta = \dot{M}_{\rm out}/{\rm SFR} = 4-15$, suggesting highly efficient feedback at very early times. Finally, we show that, if outflows are one of the mechanisms regulating star formation in JADES-GS-z14-0, the instantaneous star-formation efficiency in massive haloes is constrained to $\epsilon_\star\lesssim0.08$. These results support a scenario in which outflows play a crucial role during the earliest phases of galaxy formation. Comparing our results with the current theoretical galaxy formation model, we conclude that a combination of moderate star-formation efficiency and reduced dust attenuation can account for the emergence of luminous galaxies at the highest redshifts.

In this second paper in a series dedicated to characterizing shear layers via 2D hybrid (kinetic ions -- fluid electrons) simulations, we study the dynamical role of nonthermal particles (cosmic rays, CRs), either spontaneously generated or pre-existing. We initialize Kolmogorov-type sinusoidal velocity shear flows unstable to the Kelvin--Helmholtz instability, which evolve nonlinearly into turbulence. Particles with large gyroradii act as long-range messengers that promote momentum exchange between layers, hence introducing a form of CR viscosity. Even when not energetically dominant, increasing the CR energy density generally enhances momentum transfer, provided that their gyroradii are smaller than the shear lengthscale. We consider flows ranging from subsonic to supersonic and assess the rate of shear dissipation, the partition of the initial kinetic energy among heating, thermal ion acceleration, CR reacceleration, and magnetic-field amplification, and the maximum energy attained by accelerated particles.

Protostellar winds can theoretically lift solids from the planet-forming disks, but direct evidence for launched dust has been scarce so far. Numerous atomic lines that are unique to mid-infrared (IR) wavelengths reveal refractories eroded from dust grains and provide information on wind properties in the earliest stages of the star formation process. We present JWST/MIRI-MRS spectral imaging of the inner 2000 au of the BHR71-IRS1 blueshifted side of the outflow. Atomic line intensities are compared to shock models to constrain the physical conditions and elemental abundances of the outflowing gas. Dust continuum maps are constructed from PSF-subtracted cubes, and the dust spectral energy distribution is analyzed. The ionized central jet of BHR71-IRS1 is spatially resolved and imaged for the first time, revealing a unique inventory of refractory, volatile, and noble-gas fine-structure lines (Fe, Ni, Co, Cl, S, Ne, Ar). The emission is concentrated along four bright knots that wiggle along the jet axis. PSF-subtracted continuum maps reveal extended mid-IR continuum emission co-spatial with the jet bullets and within the H$_2$-traced outflow cone. Spectral energy distributions along the jet are fit together with the extinction, revealing a warm (200-400 K) and a cold (70-90 K) dust component. Shock modeling constrained by the mid-IR lines indicates a decline in shock velocity from 70 to 35 km s$^{-1}$ and pre-shock density from $>$10$^5$ to $ 4\times 10^4$ cm$^{-3}$ with distance from the protostar. Gas-phase Fe and Ni are measurably depleted relative to Solar abundances, consistent with a substantial fraction of refractories remaining locked in grains in spite of the shocks. These JWST observations provide direct evidence that dust is launched in a Class 0 jet and at least partly survives shock processing.

We introduce a new kinematic disturbance parameter, $\Delta V_{\star-g}$ (pronounced `DVSG'), which takes advantage of integral field spectroscopy (IFS) to quantify differences between a galaxy's stellar and gas velocity maps. The motivation behind $\Delta V_{\star-g}$ is to capture disturbances in the kinematics of a galaxy that might be missed by alternative methods, while also attempting to minimize bias towards galaxy properties or features of the IFS data. We first detail the reasons for introducing this parameter, and explain how the $\Delta V_{\star-g}$ value of a galaxy can be calculated. We then present initial results using $\Delta V_{\star-g}$ to quantify the kinematic disturbance of obscured active galactic nuclei (AGN) found in the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. We find that there is no statistically significant difference between the $\Delta V_{\star-g}$ distributions of AGN and a control sample (matched in mass and redshift) of inactive galaxies. This suggests that AGN triggering may not be preferentially caused by any distinct kinematic disturbance process, or combination of processes, beyond those observed in inactive galaxies.

We summarise the properties and nature of a peculiar group of B-type stars called FS CMa stars. These stars show the B[e] phenomenon, i.e., their spectra exhibit both forbidden emission lines and infrared excess. Such properties point to an extended circumstellar gas and dust component. Although the phenomenon has been explained in most B[e] stars, the origin and nature of FS CMa stars is disputed. Here, we focus on the merger hypothesis, for which evidence has recently been discovered.

We anticipate that hundreds of thousands of distant, strongly gravitationally lensed sources will be detectable with the European Space Agency's (ESA) Euclid mission and the Rubin Observatory Legacy Survey of Space and Time. We consider the virtues and shortcomings of the Singular Isothermal Elliptical Potential (SIEP) with Parallel External Shear (XS_||) for these systems. Its principal virtue is that it admits an analytic forward model that gives image positions and magnifications as functions of the source position (and shape for extended sources). Preliminary experiments suggest a speed-up of a factor in excess of 10,000 compared with conventional models that instead map from the image plane to the source plane and require iteration to converge upon a unique source. A second virtue is that the Witt--Wynne geometric representation of SIEP+XS_|| permits the quick visual verification of the model's adequacy for a particular lensed system. Unfortunately, the model's strictly elliptical lens equipotential is inconsistent with strictly elliptical surface mass density contours. The Witt--Wynne construction might nonetheless yield a sufficiently good first approximation to accelerate convergence to one's preferred lens model.

'Lava worlds'-Earth-sized planets hot enough (Teq >~ 1100 K) to melt their dayside silicate surfaces-have emerged as promising candidates for atmospheric detection and characterization. Thermal emission observations show an apparent dichotomy: the hottest lava worlds have colder daysides than the temperature of a maximally emitting bare rock, indicating the likely presence of thick and/or reflective atmospheres while the coldest ones do not. However, where in instellation flux this potential bifurcation occurs is uncertain. We present a JWST MIRI LRS eclipse of the ultra-short period (USP) lava world HD 3167 b (Teq = 1786 K, R = 1.6 Rearth, P = 0.96 d) that helps bridge this gap. We measure the white light eclipse depth to be 38 +/- 11 ppm, more than 5 sigma lower than the expected eclipse depth of a dark, maximally hot bare rock. We use this to derive a dayside brightness temperature that is best explained by the presence of an atmosphere that cools the dayside by reflecting incoming starlight and/or efficiently redistributing heat to the planet's nightside. An atmosphere is further compatible with the planet's slight under-density compared to an Earth-like composition. The corresponding dayside emission spectrum is not precise enough to constrain atmospheric composition, motivating follow-up spectroscopic observations with JWST NIRSpec. Lastly, we use our observation and existing data to refine key planetary parameters of the HD 3167 system. HD 3167 b is currently the least irradiated USP super-Earth with evidence for an atmosphere.

Recent population-level studies of sub-Neptune atmospheres have identified a parabolic trend in transmission spectrum amplitude for planets with Teq ~ 500-800 K. While the trend has been commonly attributed to hydrocarbon aerosols, we lack a first-principles explanation of its underlying chemical mechanism. Previous work has focused on the role of methane photolysis and subsequent polymerization, but with limited reaction networks that truncated at C2-species and couldn't reproduce the observed parabolic trend. In this work, enabled by a computer-automated, rate-based chemical network generator, we construct the most comprehensive carbon reaction network for exoplanet atmospheres to date. We explicitly model the formation of polycyclic aromatic hydrocarbons (PAHs), which are well established as soot precursors in combustion chemistry. We calculate the chemical timescales of hydrocarbon species through an eigenvalue timescale method and model their quenched abundances across a range of C/O, metallicities, and Teq. In this framework, the deep atmosphere acts as a "soot factory" analogous to a combustion engine, transporting the ingredients for hydrocarbon aerosol formation to the JWST-observable region of the atmosphere, where it may be further augmented by photochemistry. We find that the predicted abundances of PAHs peak near 600 K, and fall off toward higher and lower Teq, in agreement with the observed parabolic trend and existing JWST and HST observations. We also show that PAH abundances are expected to vary with C/O and metallicity, thus providing a natural explanation for observed diversity among planets with similar Teq.

In protoplanetary discs, the presence of dust traps can significantly alter the transport of solids from the outer to the inner regions, and hence they are often invoked as an explanation for the chemical diversity of inner discs observed with JWST (e.g., varying oxygen abundances and C/O ratios). As a detailed treatment of dust transport around dust traps is computationally expensive, earlier works investigating the impact of outer traps on the inner disc composition have often used simplified dust models representing the size distribution with a single effective size and drift speed. In this paper, we revisit the impact of outer traps on dust transport using the state-of-the-art one-dimensional dust evolution code \texttt{DustPy}, which simulates the transport and evolution of dust particles including detailed coagulation and fragmentation. We quantify and map the leakiness of dust traps across a broad parameter space, performing over 300 simulations while varying the disc viscosity, turbulence strength, planet mass and location, and dust fragmentation velocity. We find that dust traps are leakier than previously thought, on a broader parameter space, such that most outer traps (r > 5 au) will result in a long-lived O-rich inner disc with gas-phase C/O < 1. In similar conditions (e.g., carved by the same planet mass), we find inner traps are much leakier than outer traps, though their relative efficiency in reducing the pebble flux is time-dependent. Highly blocking traps altering the inner disc composition dramatically (leading, e.g., to C/O > 1) are possible to set up but necessitate low viscosity and weak turbulence, along with efficient planetesimal formation by the streaming instability. In that case, we find that is the formation of planetesimals, rather than the dust traps themselves, that is capable of significantly altering the inner disc composition.

Massive quiescent galaxies (QGs) at high redshifts are likely progenitors of massive elliptical galaxies in the local Universe. Recent observations, such as the discovery of QGs in overdensity (galaxy groups and proto-clusters) at high redshifts, have highlighted the importance of the relationship between star formation activity in galaxies and the surrounding environment. We spectroscopically confirm a galaxy group associated with a massive QG at $z_\mathrm{spec}=4.53$ from the Lyman break feature using Subaru/FOCAS. This group consists of at least three star-forming galaxies within 150 pkpc of the QG, which suggests the importance of physical association with other galaxies for galaxy quenching. In order to understand the role of the surrounding environment, we also perform a statistical analysis to characterize the typical environment of QGs at high redshifts. By selecting QGs using the SFR-based selection in the COSMOS field, we find that the fraction of QGs is higher in group or cluster-like environment at $12$ demonstrates that environmental quenching (primarily driven by galaxy mergers and interactions) plays a major role in the formation and evolution of massive QGs at high redshifts.

RR Lyrae stars are widely considered tracers of ancient (greater than 10 Gyr), metal-poor stellar populations. However, recent kinematic and photometric studies suggest the existence of a metal-rich RRL sub-population associated with the thin disc and intermediate ages (approximately 2-5 Gyr), challenging canonical evolutionary models. We aim to provide the first spectroscopic confirmation of a member of this elusive population. Specifically, we target a metal-rich RRL candidate recently identified photometrically as a member of the intermediate-age open cluster Trumpler 5. We obtained high-resolution spectroscopy using PEPSI at the LBT and GHOST at Gemini-South Telescope. We measured radial velocities from multiple epochs to constrain cluster membership and derived detailed chemical abundances (Mg, Ca, Sc, Ti, Mn, Fe, Cu, Zn, Y, and Ba) to compare the RRL's composition with that of red clump stars in the cluster. The RRL's systemic velocity Vgamma = 50.57 +0.78/-0.36 km/s is in excellent agreement with the cluster mean V = 50.76 +/- 0.49 km/s. Combining RVs, proper motions, and parallax, the probability of the star being a background interloper is negligible (approximately 0.002%, better than 4 sigma). We derived a metallicity of [Fe/H] = -0.40 +/- 0.05, matching the cluster value. While most abundance ratios (Mg, Ti, Mn, Cu, and Zn) align with cluster members, the RRL exhibits significant depletion in Ca, Sc, Y, and Ba. Notably, [Sc/Fe] is under-abundant by approximately 0.6 dex relative to the cluster stars, following trends seen in field metal-rich RRLs. We provide strong constraints on the membership status between an RRL variable and an intermediate-age open cluster [...]

Phylogenetic methods, traditionally used in biology to trace the evolutionary relationships among species, are emerging as a powerful framework to reconstruct evolutionary processes in galaxies from chemical information. We apply galactic phylogenetics to study the chemical evolution of stellar populations in distinct regions of a simulated disc galaxy, assessing its capability to unveil assembly histories. We used a high-resolution simulation that follows the chemical enrichment of an isolated disc galaxy, by different stellar progenitors. We track gas particles as they turn into stars and inherit their parent gas chemical composition. Target particles are selected to store the chemical history of each chemical element considered in the simulation. Two regions were analysed: an inner ring, influenced by early bar-driven inflows, and an outer ring, shaped by spiral arms. We built phylogenetic trees for stellar populations in each region and quantified their structure using the Corrected Colless index, a standard metric of tree balance used in biology. The inner ring tree reveals a compact clade of old stars enriched by rapid SNII feedback, followed by a hierarchical sequence with increasing SNIa and AGB contributions. In contrast, the outer ring exhibits more symmetric, caterpillar-like trees with smoother abundance gradients, consistent with more prolonged star formation and efficient local mixing. Chemical enrichment rates corroborate these trends, showing fast early enrichment in the inner ring and gradual, spatially extended enrichment in the outer disc. The structural indices differ significantly between the two regions and converge robustly even for modest stellar samples (NSSP = 100). Galactic phylogenetics provides a novel and complementary tool to decode the fossil record of galaxies.

With n-body simulations and analytic approximations we study the dynamics and stability of low eccentricity misaligned test particles around binary systems with varying mass fraction and eccentricity. General relativity (GR) plays a primary role in determining the motion of an outer particle since it drives apsidal precession of the binary orbit. The effects of GR can drive particle instability close to the binary orbit, depending upon the binary parameters and the initial inclination of the particle. For the binary parameters we consider, we find instability up to a semimajor axis of about 8 ab, where ab is the binary semimajor axis. In particular, we identify and analyse three different regions of instability that are driven by GR in the phase plane of the initial semimajor axis and the initial inclination of the particle. The results have implications for circumbinary orbits and circumbinary disks on all scales, but are particularly important around supermassive black hole binaries where the effects of GR can be strong.

OB runaway stars are massive stars moving through interstellar space at high velocities (up to 200 km/s), produced by dynamical ejections in young massive clusters or supernova explosions in massive binaries. They can travel several hundred parsec before exploding as supernovae, affecting the dynamical and chemical evolution of the Galaxy. The Vel OB1 association, one of the largest OB associations, hosts about 20 O-type and more than 50 B-type stars. We aimed to identify OB runaways in this region, quantify their number, identify their parent clusters, and understand their production channels and impact on the surrounding medium. Using Gaia DR3 coordinates, parallaxes, and proper motions, we identified OB runaways by measuring their peculiar velocity. We inspected infrared WISE images to identify wind bow shocks and reconstructed runaway trajectories to locate parent clusters and estimate travel times. We identified six young stellar clusters hosting most of the massive-star population in Vel OB1 (distance 1.6-2.1 kpc; age 1-10 Myr) and derived a threshold velocity of 15 km/s to classify runaways. We identified 25 OB runaways (including HMXB VelaX-1) and one F-type runaway. We detected 16 arc-like features, six associated with runaways selected by peculiar velocity, and ten bow shocks aligned with runaway proper motions. Parent clusters are identified for seven runaways, most likely ejected dynamically. The runaway fraction is about 30%. Wind bow shocks from OB runaways reveal valuable information on local ISM conditions.

We discuss the effective field theory of large scale structure in terms of a single scalar degree of freedom, corresponding to the velocity potential of the matter fluid in a $\Lambda$CDM universe. This cosmic ``pion'' field is nonlinearly related to the overdensity and the gravitational potential, and corresponds to the Goldstone boson of spontaneously broken spacetime symmetry, allowing us to organize perturbation theory in a systematic way while keeping the symmetries manifest. We develop the effective field theory of the pion field to next-to-leading order, and we use it to calculate the corrections to the power spectrum and to check that these are consistent with the consistency relations of spontaneously broken spacetime symmetry. We compare our results against computer simulations for the evolution of large scale structure in the pion field picture, and we make use of N-body simulations to measure EFT coefficients and analyze the growth of additional degrees of freedom in the deep nonlinear regime. We conclude with a discussion of how the pion field picture may help suggest new variables for analyzing simulations and experimental surveys of large scale structure.

Cosmology is living through fascinating times, where new observations from ground and space telescopes are questioning the established paradigm, the so-called Lambda Cold Dark Matter model. The particle nature of Dark Matter is severely constrained by underground experiments, while recent observations by galaxy surveys indicate that the cosmological constant (Lambda) may not be constant after all. Furthermore, observations at high redshift of fully-formed galaxies with massive black holes at their centers by the James Webb Space Telescope, as well as black holes with unexpected properties observed by gravitational wave detectors LIGO-Virgo, are driving an in-depth revision of our assumptions in models of structure formation and the evolution of the universe. I propose to explore two new paradigms to account for Dark Matter and Dark Energy, based on known physics without the need for new particles nor new degrees of freedom. I will extend the primordial spectrum of fluctuations to small scales with new statistical properties to provide a viable Primordial Black Hole scenario for Dark Matter, and will include non-equilibrium thermodynamics in the expanding universe, in the form of General Relativistic Entropic Acceleration, to explain Dark Energy. My proposal could provide a unified explanation for a plethora of interrelated multi-epoch, multi-scale and multi-probe observations from present and future Gravitational Wave detectors, Large Scale Structure observatories and Cosmic Microwave Background experiments. It emphasizes the need to develop new theoretical ideas hand-in-hand with observations to acquire a deeper understanding of our universe. If these ideas are correct, they will open a new window into the early universe and a new fundamental understanding of gravity in the late universe.

We investigate the current viability of a well-known coupled dark energy scenario in which fermionic cold dark matter (DM) interacts with a spin-0 dark energy component through a non-trivial field dependence of the DM mass. This ultra-light scalar mediates a fifth force between DM particles, which can leave signatures on cosmological scales. We use state-of-the-art data on the cosmic microwave background from Planck's CamSpec likelihood, baryon acoustic oscillations from the second DESI data release as well as the supernovae of Type Ia (SNIa) from Pantheon+ and DES-Dovekie. We perform the analysis considering both a flat potential and a Peebles-Ratra (PR) potential for the scalar field in order to assess the impact of the potential slope on the fitting performance of the model. While for a constant potential the scalar field dynamics is insensitive to the sign of the coupling parameter $\beta$, the PR potential breaks the existing symmetry in the solutions at late times and could induce a difference at the phenomenological level between positive and negative values. We study for the first time if it is actually the case, finding no important asymmetry in the fitting results. In the light of the aforesaid datasets, we find in all cases a peak at $|\beta|\sim 0.03$ - less pronounced than reported in some recent works -, excluding the no-coupling scenario at $\sim 95\%$ CL. The model is able to explain an effective crossing of the phantom divide, with the equation-of-state parameter lying within the $2\sigma$ bands of model-agnostic reconstructions. Our results are very robust under changes in the SNIa sample used in the analysis and is not significantly altered when we replace a constant potential with the PR one, although the latter is crucial to produce the aforesaid crossing. In passing, we also provide constraints obtained with the PR potential in the uncoupled case.

Meteorites are classified as either non-carbonaceous- (NC) or carbonaceous (CC), representing bodies that likely formed in the inner- or outer solar system, respectively. Despite its location in the inner solar system, the Earth is thought to contain either minor- (~6 %) or substantial amounts (~40 %) of outer solar system material. However, because neither interpretation leverages variations among multiple isotopic systems simultaneously, Earth's provenance remains equivocal. Here, we examine variations in 10 nucleosynthetic isotope anomalies among planetary- and meteorite parent bodies to show that the linear extension of an array defined by NC bodies in any two isotopic anomalies always intersects the observed isotopic composition of the bulk silicate Earth to within 1 standard deviation. The Earth therefore formed exclusively from inner solar system material whose composition did not vary over the course of accretion and was, on average, unlike that of any chondrite. Extension of the NC array yields isotopic compositions for Mercury and Venus that are more extreme than for Earth, implying a spatial or temporal gradient during the formation of the terrestrial planets.

The jets and lobes of nearby radio galaxies are ideal laboratories to explore the spectral and dynamical evolution of the radio-emitting plasma that emanates from active galactic nuclei. Here, we present a high-resolution radio continuum study of the low-redshift (z=0.01713), Fanaroff-Riley I (FRI) radio galaxy 3C 449 performed by combining radio data at 145 MHz acquired with the LOw Frequency ARray (LOFAR) and archival Very Large Array (VLA) data at 1365, 1485, 4985, and 8485 MHz. Our LOFAR maps have angular resolutions of 20"x20" and 6.0"x6.0", and show the full extent of the known radio emission at the highest angular resolution to date. Our spectral index maps show the distribution of the spectrum in the 145-8485 MHz frequency range over a region that extends beyond 2.5'. The average 145-8485 MHz spectrum is consistent with a single power law and stays approximately constant over the inner ~50" of both jets. Beyond ~50", on both source sides, the higher-frequency spectrum steepens, indicating no significant downstream sites of particle acceleration. By modelling the spectrum under the assumption of equipartition and of a constant magnetic field across the source, we derive maps of the highest radiative age of the particles. At the outer edges of both the northern tail and southern lobe, the spectral age is ~150 Myr. If the latter age were representative of the dynamical source age, the average expansion speed of both jets during the source lifetime would be supersonic, with Mach numbers of M~4.1 and M~2.8 for the northern and southern jets, respectively. As numerical simulations of FRI jets suggest that the source's current expansion is subsonic, the high Mach numbers might arise either from the source being highly supersonic for a small fraction of its lifetime or from a severe underestimation of the spectral age due to particle acceleration on scales of hundreds of kpc.

Self-gravitating systems of $N$ particles are chaotic. We wonder how chaotic the Galaxy is, and what the consequences are. We therefore simulate the dynamical evolution of a galaxy-scale distribution of point masses in order to measure the degree of chaos in such a system. These calculations were performed using the softened gravitational $N$-body tree-code Bonsai, with up to 40 million equal-mass particles. Smaller simulations were performed to establish the scaling of the Lyapunov time $t_L$ with $N$. We establish the relations between the degree of chaos, the number of particles, and the softening length in the gravitational force calculation of large-scale $N$-body simulations. The moment the bar forms appears insensitive to infinitesimal perturbations to the initial realisation. In contrast, the bar strength and its further evolution sensitively depend on such perturbations. Interestingly enough, the run-to-run variation in the bar strength has its maximum around the maximum bar strength, and drops to the moment the bar buckles. The galaxies we simulated are highly chaotic, but the softening in the simulations suppresses chaos. Still, our models show considerable variations in the macroscopic behaviour due to infinitesimal perturbations to the initial conditions. Real galaxies, however, should be orders of magnitude more chaotic than our simulations, and we are unable to quantify their consequences. Smooth galactic potentials to study individual stellar orbits should be handled with caution on timescales longer than the Lyapunov time. Extrapolating to the number of stars in the Galaxy, ignoring planets and other minor bodies, we conclude that the Milky Way-size galaxies are chaotic on a timescale $\lesssim 0.1$ Myr.

Be/X-ray binaries, the most common class of high-mass X-ray binaries, are characterised by OBe companions, but display a rich variety of X-ray behaviours. One of the most atypical systems is X0726-260, which also has the earliest optical counterpart among the whole Milky Way and Magellanic Cloud sample. We intend to improve the characterisation of the optical counterpart, LS 437, and to constrain the physical mechanisms responsible for the anomalous properties of X0726-260. We analyse high-quality, high-resolution optical spectroscopy of LS 437 with standard quantitative methodology to derive stellar parameters. We also make use of archival X-ray monitoring. We derive a moderate projected rotational velocity v sin i $\approx$ 155 km/s and a spectral type O7.5 Ve (Teff = 36 000 K), which makes LS 437 substantially earlier than any other Oe star in an X-ray binary. At this spectral type, the stellar wind likely contributes significantly to mass accretion, and the X-ray light curve is strongly suggestive of an orbitally modulated wind accretor. The source shows marked carbon depletion, whereas nitrogen is only slightly above solar abundance. LS 437 is the earliest Oe star known in the Galaxy, alongside HD 155806. Long-term X-ray lightcurves of X0726-260 strengthen the identification of a persistent 34.5 d periodicity as the orbital period, demonstrating that the X-ray emission is orbitally modulated and no X-ray outbursts have occurred over the past 30 years. Likewise, emission features in the optical spectrum indicate a remarkably stable circumstellar disk, with no sign of major structural changes over the past 40 years. All these characteristics set X0726-260 clearly apart from typical Be/X-ray binaries.

Previous theoretical studies indicate that the inner disk in dwarf novae evaporates into a high-temperature, optically thin, and geometrically thick accretion flow during quiescence, with the inner edge moving toward the white dwarf at the onset of an outburst. We incorporate this process into the numerical model developed by Kimura & Osaki (2023) and test the code on two representative dwarf novae, SS Cyg and U Gem. By modeling the inner accretion flow, we calculate the optical, ultraviolet (UV), and X-ray luminosities. Our results show that evaporation suppresses the inside-out outbursts without requiring a radially dependent viscosity parameter in the cold state. The observed time delay between the rise in UV luminosity and the onset of the optical outburst is more than one day, which is successfully reproduced when the inner disk is truncated at several $\times 10^{9}$ cm in the standard evaporation model. However, while the modeled accretion rate at the inner disk edge in U Gem accounts for its quiescent X-ray luminosity, the rate in SS Cyg remains insufficient. This discrepancy in SS Cyg suggests that SS Cyg may require either more efficient evaporation or an additional mass supply into the coronal cavity via gas-stream overflow. By accounting for disk evaporation, our simulations offer a refined version of the disk instability model for dwarf nova outbursts that naturally explains the observed multiwavelength light curves.

The outgassing signatures of Io, Europa, Enceladus, Triton, and Io-like exomoons are the focus of this review chapter. The rocky volcanic world of Io is unique in our Solar System, with plumes reaching to hundreds of kilometres in altitude. Io-like exomoons could leave signatures strong enough to be detected with ground-based telescopes. The icy moons Europa and Enceladus, with their subsurface oceans, are currently the best candidates for life. Triton is different in many ways and raises unexplored questions. Our knowledge of these active moons is derived from space- and ground-based observations. To understand their origin, we discuss moon formation in general, before examining evidence and signatures of plumes on these moons. Given the accessibility of subsurface oceanic material through the occurrence of plumes, we expand on possibilities to investigate biosignatures.

Changing-Look AGN (CL-AGN) are AGN which transition between Seyfert types, challenging AGN unification models. Most CL-AGN have been identified via repeat spectroscopy, making it difficult to determine the duration and magnitude of the CL-AGN transition. As such, the physical mechanisms behind this transition are still unknown. We use synthetic photometry in combination with ZTF light curve data to develop a new criterion to identify photometric CL-AGN transitions based on changes in g-band magnitude and g-r color. We find that a CL-AGN criterion of $| \Delta g| > 0.4$ mag and $| \Delta (g-r)| > 0.2$ mag recovers a photometric transition in $9.6^{+4.9}_{-3.4}\%$ of CL-AGN hosts over the six-year ZTF survey, including a candidate repeating changing-look event in SDSS J084957.78+274728.9. Using simulated AGN light curves, we estimate the false positive rate among the simulated Seyferts to be $1.6^{+0.19}_{-0.17}\%$. We find that the rate of similar flares among Type 1 Seyferts is $1.2^{+0.87}_{-0.50}\%$ , and among Type 2 Seyferts is $\leq 0.39\%$ over six years. Photometric CL-AGN transitions last between 21 and 560 days, with a median duration of 360 days, consistent with the thermal or orbital timescales for AGN disks. We do not detect a correlation between black hole mass and transition duration, likely due to the small sample of detected photometric transitions. This method can be applied to the upcoming Legacy Survey of Space and Time to identify CL-AGN candidates and test theories of their origins

Turbulent mixing remains one of the primary uncertainties in the modeling of stellar interiors. In stellar evolution simulations, regions where mixing occurs are typically identified using instability criteria. A particularly interesting situation arises when nuclear reactions produce inversions in the mean molecular weight within stellar interiors. Under these conditions, the material can become unstable to either thermohaline or a Rayleigh-Taylor instabilities. We demonstrate that the standard criterion adopted in stellar evolution calculations does not accurately distinguish between these two regimes. We derive an alternative criterion and show that chemically driven convection in stellar interiors might be viable under much smaller mean molecular weight inversions than it is normally assumed. We investigate whether inversions in the mean molecular weight can trigger chemically driven convection above the red giant branch bump (RGBB) or during the helium core flash. We find that the inversion at the base of the convective envelope above the RGBB is too weak and short-lived to sustain steady-state convection. In contrast, rapid carbon production at the base of the He-flash-driven convective zone can maintain a steady chemically driven convective region. This process could significantly alter our understanding of the He-core flash and warrants further study.

We examine the positions of substellar objects in mid-infrared color-magnitude and color-color diagrams to distinguish between cloudy and cloud-free atmospheres. Using Spitzer mid-infrared spectra of 113 M5-T9 ultracool dwarfs, we derive synthetic photometry for the JWST MIRI F560W, E'770W, F1000W, and F 1280W filters, which cover key absorption features including the ~9 um silicate signa-ture. We find that diagrams involving F770W and F1000W best separate L-type objects with silicate clouds in their photospheres. L dwarfs with mE77ow - mF1000w < 0.03 mag are seven times more likely to host cloudy atmospheres. Diagrams using F1000W and F1280W are less informative due to the lower signal of the spectra at long wavelengths. Current model predictions struggle to reproduce the positions of cloudy, warm brown dwarfs, likely because atmospheric models underestimate the ~9 um silicate feature. Cloudy Sonora Diamondback models better match the observed trends, although this may reflect improvements capturing indirect effects of clouds on the 6.25 um water absorption feature rather than accurately modeling the silicate feature itself. Our analysis indicates that JWST MIRI photometry can efficiently identify new cloudy extrasolar atmospheres for targeted spectroscopic follow-up, optimizing the use of telescope time.

The 100,000-year problem concerns the dominant period of glacial-interglacial cycles over the past 800,000 years and their correlation with Earth's orbital eccentricity, despite eccentricity's weak influence on solar radiation. Two theories compete: the astronomical theory, in which orbital forcing drives the cycles with amplification from Earth system feedbacks, and the geochemical theory, in which internal dynamics dominate with orbital forcing synchronising oscillations. We investigate these theories using conceptual models. Augmentations to the Budyko energy balance model fail to reproduce the 100,000-year period, revealing formulation limitations. Linearised versions of existing non-linear ice volume models perform comparably to their full counterparts, indicating the data does not necessitate non-linear dynamics. We develop two simple linear models: a feedforward model aligned with the astronomical theory and a feedback model aligned with the geochemical theory. The feedforward model reproduces the ice volume record well and offers a novel explanation for the absence of eccentricity's 400,000-year period, arising from oceanic heat storage and tropospheric energy responding with differing phase lags. Conservative estimates show bulk ocean temperature variation can be explained by eccentricity alone, challenging the geochemical theory's core assumption. We also show that widespread use of Q65 may bias models towards geochemical explanations by underrepresenting eccentricity. The feedback model's improvement is concentrated around Marine Isotope Stage 11, suggesting this anomalous interglacial reflects Earth-based events rather than a general requirement for feedback mechanisms. We conclude that 800,000 years of glacial cycles can be largely reproduced by a linear astronomical model, emphasising the importance of parsimony when interpreting palaeoclimate data.

We model the dynamical mass density profiles of 14 strong gravitational lens galaxies from the Sloan Lens ACS (SLACS) sample using spatially resolved kinematics obtained from Keck KCWI integral-field spectroscopy. We use the Jeans Anisotropic Modeling (JAM) method, combining 2D kinematic maps with joint constraints from lens models from Hubble Space Telescope imaging. We use informative priors on the anisotropy and intrinsic shape from local galaxies to help break the residual mass-anisotropy degeneracy (MAD). We find nearly isothermal power-law total mass density slopes ($\rho_{\rm tot}\propto r^{-\gamma}$) for the sample with a mean of $\gamma = 2.04\pm0.02$ with intrinsic scatter of $0.08^{+0.03}_{-0.02}$. We fit explicitly for deviations from the pure power-law form that are fully sensitive to the mass-sheet degeneracy (MSD) and constrain the value of the mass-sheet parameter $\rm \lambda_{int}$ for each individual galaxy to an average precision of 5.8%. The mean value of $\rm \lambda_{int}$ for the sample is $1.01\pm0.03$, with intrinsic scatter of $0.11\pm0.03$. Values of $\rm \lambda_{int}$ for individual objects and the scatter in the sample are consistent to $1\sigma$ uncertainty with those found by the Time-Delay COSMOgraphy collaboration's 2025 milestone analysis, which used a spherical analysis of the same dataset, but azimuthally averaged. We thus conclude that on average power-law mass profiles are a good first-order description of the SLACS sample and do not introduce measureable bias in time-delay cosmography. However, our analysis indicates that more flexible mass models should be able to reproduce the highly detailed kinematic datasets more accurately.

Observations of ultra-dense substructures in strong lensing systems challenge the standard cosmological model at small scales. Self-interacting dark matter (SIDM), as an alternative to the cold and collisionless dark matter (CDM) of the standard cosmological model, provides a natural mechanism for forming such structures via gravothermal core collapse. We show that strong gravitational lensing of fast radio bursts (FRBs) provides an effective approach to detecting these substructures and probing dark matter self-interactions. Core-collapsed SIDM halos exhibit steeper central density profiles than CDM halos, enhancing the lensing cross section and producing longer time delays between FRB images. We compute lensing properties of core-collapsed subhalos and host halos, including maximal impact parameters and time-delay distributions. We demonstrate that future all-sky monitors, such as BURSTT, SKA2-Low, and SKA2-Mid, which are expected to detect $10^{5}$--$10^{7}$ FRBs over a decade, can measure time-delay distributions with high statistical significance. Modeling collapsed halos with a cored power-law density profile with inner slope $\gamma=3$ and assuming no excess beyond the singular isothermal sphere lens model, we show that our strategy can probe self-interaction cross section strengths of $\sigma_{\text{SI}}/m \gtrsim \min\{18,\, 40\lambda_{\text{sub}}\}\,\text{cm}^2/\text{g}$, where $\lambda_{\text{sub}}$ parameterizes the collapse time of a subhalo relative to that of the isolated case.

The Tianlai Dish Pathfinder Array (TDPA) is a radio interferometer designed to test techniques for 21 cm intensity mapping in the post-reionization universe as a means of measuring large-scale cosmic structure. Using 9 nights of observations targeting the North Celestial Pole (NCP) field, totaling approximately 107 hours of integration time, we analyze data in the frequency range 700-800 MHz (corresponding to redshift $z \sim 0.9$). We do the data format conversion, radio frequency interference (RFI) flagging, calibration, imaging and point source subtraction, and foreground removal via Singular Value Decomposition (SVD). The spherically averaged power spectrum $\Delta^2(k)$ is obtained. This work successfully establishes and validates a comprehensive data analysis framework for the TDPA. We identify key improvements including sky model refinement, increased integration time, and pipeline optimization that will enable future detection of the 21 cm signal through auto-correlation and cross-correlation with optical galaxy surveys.

We present JWST/MIRI-MRS observations of ISO-Oph 37, a highly inclined flat-spectrum ($\lesssim$1 Myr old) source, to investigate the chemical composition and dynamical origin of its inner-disk gas. The spectrum reveals a rich combination of molecular emission and absorption: H$_2$O, CO, and OH are detected in emission, while strong absorption is observed from CO, H$_2$O, CO$_2$, HCN, C$_2$H$_2$, and CH$_4$, with no detectable ice absorption features. LTE slab modeling of the absorption yields excitation temperatures of $T_{\rm ex}\sim400-600$ K and column densities of $\log N/{\rm cm}^{2}\sim16-19$, characteristic of warm gas located within the inner few au. The absorption lines are significantly blueshifted relative to the systemic velocity, with mid-IR lines exhibiting larger shifts than near-IR CO absorption. This velocity structure points to a velocity- and temperature-stratified molecular disk wind. In this framework, the absorption directly samples disk material lifted from the inner disk surface, preserving the chemical imprint of the wind-launching region. Along the line of sight, ISO-Oph 37 is unusually hydrocarbon-rich compared to other known absorption systems (GV Tau N and IRS 46), exhibiting high (C$_2$H$_2$+CH$_4$)/HCN, (C$_2$H$_2$+CH$_4$)/CO and H$_2$O/CO column density ratios, while the CO and HCN columns remain broadly typical. We find that these molecular ratios are best explained by enhancement of both hydrocarbons and water, driven by inward drift and sublimation of icy pebbles and by thermal processing of carbonaceous grains at the soot line. ISO-Oph 37 thus demonstrates that carbon-rich inner-disk chemistry can be established early in disk evolution and that it can be directly probed through molecular absorption in disk winds.

We present a method of spatially mapping microwave kinetic inductance detector (MKID) arrays, in a dark setup. MKIDs are superconducting natively multiplexed resonators which enable kilopixel arrays, such as for the proposed Probe far-Infrared Mission for Astrophysics (PRIMA). In such telescope applications one must map the spatial location of each MKID with their individual resonance frequencies. Traditional LED arrays or beam-mapping methods become increasingly difficult as pixel spacing decreases, e.g., 900 {\mu}m separated MKIDs in the spectrometer module of PRIMA. Our new mapping technique uses a cryogenic interferometer in reflection mode. As on-resonance signals reflect from an MKID, they accrue a phase proportional to the path-length, exactly corresponding to their physical distance on the feedline. Specifically, we use a superconducting transmission line that has nonlinear kinetic inductance. The slow-wave structure of this nonlinear device is designed to have a signal speed of 0.64% the speed of light, enabling a compact system. Current biasing this line allows for varying the wave speed and ensuring that the phase measured is periodic within a nulling interferometric mode. Using this setup, we measure a length ordering that reflects the bimodal MKID distribution of a 44 pixel array of MKIDs designed for PRIMA which contains the same spacing as the final kilopixel array design.

The exponential growth of data from modern radio telescopes presents a significant challenge to traditional single-pulse search algorithms, which are computationally intensive and prone to high false-positive rates due to Radio Frequency Interference (RFI). In this work, we introduce FRTSearch, an end-to-end framework unifying the detection and physical characterization of Fast Radio Transients (FRTs). Leveraging the morphological universality of dispersive trajectories in time-frequency dynamic spectra, we reframe FRT detection as a pattern recognition problem governed by the cold plasma dispersion relation. To facilitate this, we constructed CRAFTS-FRT, a pixel-level annotated dataset derived from the Commensal Radio Astronomy FAST Survey (CRAFTS), comprising 2{,}392 instances across diverse source classes. This dataset enables the training of a Mask R-CNN model for precise trajectory segmentation. Coupled with our physics-driven IMPIC algorithm, the framework maps the geometric coordinates of segmented trajectories to directly infer the Dispersion Measure (DM) and Time of Arrival (ToA). Benchmarking on the FAST-FREX dataset shows that FRTSearch achieves a 98.0\% recall, competitive with exhaustive search methods, while reducing false positives by over 99.9\% compared to PRESTO and delivering a processing speedup of up to $13.9\times$. Furthermore, the framework demonstrates robust cross-facility generalization, detecting all 19 tested FRBs from the ASKAP survey without retraining. By shifting the paradigm from ``search-then-identify'' to ``detect-and-infer,'' FRTSearch provides a scalable, high-precision solution for real-time discovery in the era of petabyte-scale radio astronomy.

White-light flares (WLFs) are crucial for understanding the energy transport and heating processes in the lower solar atmosphere. Systematic studies are highly necessary. However, most WLFs are very weak and difficult to detect. To address this, we propose a new method of detecting WLFs. Through the observations of SDO/HMI, the light curve of each pixel in the flaring region can be obtained. By subtracting the slowly varying background, we obtained a series of rapidly varying radiative pulses. Pixels for which radiative pulses during flares significantly exceed those occurring before and after the flare were identified as WL emission regions. We applied our method to the detection of the X2.2 flare on September 6, 2017 and validated the method. We found that the WL emission in this flare exhibits two phases, and that different regions show distinct WL emission properties. We also detected the WL emission in all the flares (1 X-class, 2 M-class, and 20 C-class) occurred in active region NOAA 12887. It was found that 15 of the 23 flares are WLFs (1 X-class, 2 M-class, and 12 C-class). The occurrence rate of WLFs in this active region is $\sim65\%$. Surprisingly, the occurrence rate of WLFs in C-class flares even reaches up to $60\%$. It should be noted that most of these C-class WLFs are below C5.0. In addition, a C1.0 WLF was identified; this is the lowest GOES-class event with confirmed WL emission to date. These results demonstrate that WL emission is ubiquitous in most flares, even down to C-class events.

We report a photometric and asteroseismic analysis of KIC 3868420, a newly identified high-amplitude $\delta$ Scuti-$\gamma$ Doradus hybrid star located in the Hertzsprung gap - a short-lived and rarely observed post-main-sequence phase. Using four years of Kepler long-cadence photometry, we detect 36 significant frequencies, including 11 independent modes spanning both low- and high-frequency regimes. Grid-based modeling with MESA and GYRE, including rotation, shows that five independent frequencies match a combination of radial p-modes and non-radial g-modes, supporting its hybrid nature. The best-fit models yield an evolved post-main-sequence star ($M \sim 2.26 - 2.30 M_\odot$, $R \sim 4.41 - 4.43 R_\odot$, $\tau \sim 5.4 \times 10^{8}$ yr), although degeneracies from rotation and unknown inclination preclude a unique solution. KIC 3868420 thus represents a rare example of a high-amplitude hybrid pulsator in the Hertzsprung gap, demonstrating that high-amplitude pulsation can coexist with mixed p/g-mode behavior in a rapidly evolving star. This finding highlights the value of space-based photometry for identifying such rare objects.

We present a deep learning framework to enhance the identification of Ly$\alpha$ emitters (LAEs) in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), an untargeted spectroscopic survey of LAEs at $1.9 < z < 3.5$ without imaging pre-selection. We primarily address the low signal-to-noise ratio (S/N) regime ($4.8 \leq \mathrm{S/N} \leq 5.5$), where LAE candidates suffer from substantial noise contamination. To distinguish LAE candidates from artifacts and sky residuals, we employ a convolutional neural network (CNN) trained on two-dimensional spectral images of single emission lines. The training sample is constructed from the HETDEX COSMOS catalog, with external validation from ancillary observations and our participatory science project, \textit{Dark Energy Explorers}. For small-format, low-resolution spectroscopic data, the model achieves a balanced accuracy, precision, and recall of $94.1\%$, $97.5\%$, and $97.5\%$, respectively, in the high-S/N regime ($\mathrm{S/N}>5.5$), and $85.1\%$, $78.2\%$, and $84.4\%$ in the low-S/N regime. Using HETDEX LAEs independently identified by DESI spectroscopy, the model recovers $99\%$ and $93\%$ of the high- and low-S/N LAEs, respectively. Visual attribution indicates that the CNN attends to smooth, spatially extended central emission in true positives and to irregular or noisy features in true negatives. Applied to the full HETDEX catalog, the CNN enables an S/N threshold down to 4.8 by suppressing spurious spikes across $z\sim 1.9$--$2.5$ in the redshift distribution. Our approach facilitates HETDEX cosmological analyses by mitigating false positives in galaxy clustering and highlights the value of domain-specific deep learning for refining low-S/N spectroscopic identification in untargeted surveys.

The physical origin of optical changing-look AGNs (CLAGNs), characterized by the appearance or disappearance of broad emission lines, is thought to be mainly driven by the variation of the black-hole (BH) accretion rate. In this work, we explore this issue based on a sample of {224} CLAGNs with UV-to-optical continua, where the UV radiation is more sensitive to the accretion state near the BH horizon. We find that the luminosity correlation of $L_{3000}$--$L_{5100}$ at 3000$\rm Å$ and 5100$\rm Å$ becomes steeper at low luminosities (e.g., $L_{3000}\lesssim10^{44}\rm erg/s$), where the sources with high luminosities are roughly consistent with the prediction of a standard accretion disk. At lower luminosities, the observations are more consistent with the prediction of a truncated disk. The whole sample has a median bolometric Eddington ratio of $\sim$2.2\%, which is consistent with the critical value for state transition in X-ray binaries. Such transitions can significantly alter the UV-to-optical continuum, largely due to variations in the truncation radius, even when the change in the overall accretion rate is minimal. The deficit of ionization photons resulting from an increase in the truncation radius will lead to the weakening or disappearance of broad lines, which triggers the AGN changing-look.

$\mu$ Herculis is a bright, nearby quadruple system. Its brightest member, $\mu$ Her Aa, displays solar-like oscillations, establishing the system as a crucial benchmark for asteroseismology, provided that its mass can be determined independently of stellar models. We aim to resolve the full hierarchical architecture of the system and determine precise, model-independent dynamical masses for all four components (Aa, Ab, B, and C), along with a consistent astrometric solution for the system's centre of mass. We performed a joint fit of radial velocities, relative astrometry and absolute astrometry from \textit{Hipparcos}, \textit{Gaia} DR3, and ground-based catalogues, spanning nearly three centuries. Our forward-modelling framework simultaneously constrains the Keplerian orbits of the inner Aa--Ab and B--C subsystems, the wide A--BC orbit, and the sky motion and parallax of the total centre of mass. Leveraging several complementary datasets and the decisive 2023 periastron passage of the Aa--Ab pair, we precisely determine all orbital parameters and obtain sub-percent precision on the component masses: $M_{\rm Aa} = 1.134 \pm 0.007\,M_{\odot}$, $M_{\rm Ab} = 0.2286 \pm 0.0006\,M_{\odot}$, $M_{\rm C} = 0.445 \pm 0.005\,M_{\odot}$, and $M_{\rm B} = 0.417 \pm 0.005\,M_{\odot}$. We derive a system parallax of $\varpi_{\rm CM} = 120.069 \pm 0.089\,\mathrm{mas}$ that reconciles and improves upon the individual \textit{Hipparcos} and \textit{Gaia} DR3 values.

We present a focused multiwavelength study of blue straggler stars (BSSs) in the intermediate-age open cluster Tombaugh 2, located in the outer Galactic disk, to constrain the dominant formation pathways of BSSs in a low-density environment. Cluster members are identified using Gaia DR3 astrometry through a Gaussian Mixture Model, yielding a clean sample of high-probability members. Color-magnitude diagram analysis indicates an age of 1.74 Gyr. The radial surface density profile is well described by a King model, indicating a centrally concentrated overall structure, while the cluster exhibits only weak or no clear evidence of mass segregation among its stellar populations. We identify 26 BSS candidates and 2 YSS candidates. Spectral energy distributions constructed from ultraviolet, optical, and infrared photometry reveal that 9 BSSs (32%) exhibit significant ultraviolet excess, indicating an additional hot component. Binary SED decomposition identifies stripped companions with effective temperatures Teff $\sim$ (1.5-8) $\times$ 10$^4$ K and radii R $\sim$ 0.04-0.28 R_$\odot$, consistent with proto-white dwarfs, extremely low-mass pre-helium white dwarfs, and young hot remnants formed through recent mass transfer. A slight central concentration of BSSs, together with stripped companions, suggests that binary mass transfer is an important formation channel, with no evidence for merger-driven formation. Multi-epoch VLT/FLAMES spectroscopy reveals radial-velocity variability in several systems, providing independent evidence for binarity. Our results highlight that optical-infrared photometric analyses alone may fail to detect hot compact companions, while spectroscopy and ultraviolet observations provide complementary constraints, with ultraviolet data offering a direct probe of such companions in intermediate-age open clusters.

We investigate whether the Type Ibn supernova SN 2023uqf, reported close in time and direction to the $\sim$442 TeV IceCube alert IC-231004A, is physically consistent with a shock--circumstellar medium (CSM) interaction scenario. One-dimensional radiation-hydrodynamics calculations with {\tt STELLA} reproduce the ZTF optical light curves with a dense helium-rich CSM following $\rho_\mathrm{CSM} \propto r^{-3}$ and a CSM density parameter $D'\approx 50$. Using the shock evolution and CSM conditions inferred from the optical data, we model time-dependent cosmic-ray acceleration and hadronic neutrino production during the interaction phase. The inferred shock and CSM properties open a short-lived window in which multi-PeV hadron acceleration and efficient hadronic interactions can coexist, making SN 2023uqf a plausible transient PeVatron candidate. After folding the predicted neutrino emission through the IceCube effective area, we obtain an expected number of $\sim10^{-5}-10^{-4}$ track-like events at $d = 723$ Mpc, depending on the alert selection. In the low-count regime, the model predicts a detection-time weighting for a rare event, and the detection time of IC-231004A falls within the high-weight interval while its energy scale is compatible with the modeled spectrum. Although a single event cannot establish a definitive association, our results show that the optically inferred environment of SN 2023uqf is consistent with a transient PeVatron window and illustrate how interaction-powered Type Ibn supernovae can be tested as high-energy neutrino sources.

Magnetic field amplification is an integral part of the process of particle acceleration at non-relativistic shocks. It is necessary to reach the maximum energies required by observations, especially in supernova remnants, thought to be sources of the bulk of Galactic cosmic rays. Such amplification can be caused by the acoustic instability that develops when small density perturbations interact with the cosmic-ray pressure gradient in the upstream of a cosmic-ray-modified shock. The vorticity induced by the nonlinear development of the instability may lead to turbulence, which amplifies the pre-existing magnetic fields. To study this phenomenon, we use the PLUTO code to carry out 2D (and some 3D) magnetohydrodynamical simulations of the evolution of small density perturbations in the presence of an assigned cosmic-ray pressure gradient. Adopting more realistic values of Mach number and cosmic-ray acceleration efficiency than previously assumed in the literature, we show that the acoustic instability can transform small density perturbations into large nonlinear structures while the fluid crosses the precursor region of a cosmic-ray-modified shock. We study the power spectrum of turbulent magnetic fluctuations that may be important to scatter particles. We comment on the possible constructive interference between acoustic and non-resonant streaming instabilities. We discuss limitations of previous and current numerical investigations in accessing spatial scales where turbulence is expected to turn nonlinear, and outline perspectives for future investigations.

The cubic-kilometer neutrino telescopes have opened neutrino astronomy as an observational discipline. The recent detection of KM3-230213A, the highest-energy neutrino ever observed at ~220 PeV, as a near-horizontal muon track underscores that the ultra-high-energy frontier is accessed through horizontal directions where the Earth's opacity above ~100 TeV confines the observable sky to a narrow band around and above the horizon. Yet extending general-purpose detector architectures into this regime requires disproportionate increases in instrumentation, cost, and logistical complexity. A compelling alternative is to deploy specialized detectors that target this natural geometry. POLARIS (Pacific Ocean Large Area Radial Instrumented Sparse array) is a sparse planar deep-water Cherenkov array optimized for neutrino-induced muon tracks from horizontal directions in the multi-TeV to PeV regime. By rotating the conventional vertical string layout into a radial planar configuration, the detector presents maximal cross-section to horizontal tracks while naturally suppressing the down-going atmospheric background. With only 1100 optical modules, the five-arm design reaches point source and diffuse flux sensitivities at PeV energies competitive with detectors deploying several times more instrumentation. As a dedicated $\nu_\mu$ track detector, POLARIS provides the muon-flavor channel that tau-optimized experiments such as TAMBO and Trinity do not cover, enabling full flavor composition measurements from astrophysical sources. Using the Prometheus simulation framework, this study demonstrates that targeted sparse geometries can open new discovery space at the high-energy frontier at a fraction of the cost of general-purpose arrays.

Variability is a well known phenomenon in low-mass young stellar objects, but in recent years the monitoring of methanol masers and infrared continuum emission has permitted the detection of both burst-like episodes and periodic variations also in high-mass (proto)stars. Multi-epoch studies on large samples of these objects have become possible thanks to the NEOWISE database, which surveyed the sky in the mid-IR for about a decade. Our goal is to analyse the mid-IR emission from the well studied massive protostar IRAS20126+4104 and confirm the hypothesis that such emission is periodic, as proposed in previous studies. We take advantage of the NEOWISE, ALLWISE, and Spitzer databases to obtain 24 images of the 3.4 $\mu$m emission from IRAS20126+4104 spanning 19 years, with $\sim$6 months sampling over a decade. With these data we create a light curve for each lobe of the bipolar nebulosity/outflow associated with the protostar. Our results confirm that the IR emission from IRAS20126+4104 varies regularly with a period of $\sim$6.8 yr. The period is the same for both lobes, but their emissions are anticorrelated with a phase difference of $\sim$2.5 yr. The variation is consistent with that found in previous studies for the 6 GHz CH$_3$OH masers and the near-IR emission from the lobes. After discussing four possible ``clocks'' that could determine the observed periodicity, we rule out all but a model involving rotation of the star with a spot obscuring $\sim$20% of the stellar surface. The long rotation period implies that the 12 $M_\odot$ protostar is bloated, with a radius of $\sim$200 $R_\odot$.

The shaping of planetary nebulae on their evolution from asymptotic giant branch circumstellar envelopes to their final, most often axisymmetrical, form is still a process with many unknown details. The key to understanding the whole shaping process is the study of the transition objects called pre-planetary nebulae (pPNe). In this context, modelling tools must be kept to the standard of radio telescope capabilities, so we can make the most of the data they collect. In this work we first present the newest update of the SHAPE and SHAPEMOL modelling tools, adding ten new molecular species to be reproduced together with other general improvements. Later, we put this new update into practice to study M1-92, a pPN with a rich chemistry that can provide valuable information on its origin and shaping. We created a 3D morpho-kinematical model of the nebula in SHAPE that is able to reproduce 23 line profiles from the IRAM 30m telescope and HIFI/HSO and five maps from IRAM NOEMA. The observational dataset is reproduced simultaneously under the same physical conditions, adjusting only the relative abundance of the different species. We obtained a full description of the nebula's physical and chemical properties, and we provide the total estimates for mass (0.79 $M_\odot$), linear momentum (4.10$\times10^{39}$ g cm s$^{-1}$), and kinetic energy (6.48$\times10^{45}$ erg) as well as their detailed distribution across the nebula. We also analysed the isotopic ratios, finding robust discrepancies (values of 10 versus 30) in the $^{12}$C/$^{13}$C ratio across structures depending on their age.

Accurate characterization of filamentary structures in star-forming clouds is essential for understanding star formation. Traditional methods fit observed surface density profiles $\Sigma(r)$ with slope $\gamma$ and width $H$ using the Plummer function, assuming $\beta=\gamma+1$ and $h\approx H$ for the volume density slope and width. These assumptions are inconsistent with the finite nature of filaments. We present a new fitting method that explicitly accounts for finite cylindrical geometry and establishes self-consistent empirical relationships between the parameters of $\Sigma(r)$ and those of the volume density profile $\rho(r)$ with slope $\beta$ and width $h$. The method was validated on model profiles and applied to California filaments. The slope difference $\delta\equiv\beta-\gamma$ falls below unity for shallow ($\beta\lesssim 2$) and compact profiles; $h$ and $H$ can differ by over an order of magnitude for extended filaments with shallow slopes. Accurate parameter recovery requires high resolvedness $R\equiv H/O\gtrsim 10$ (where $O$ is the beam width); at lower resolvedness, slopes are severely overestimated and filaments remain unresolved even when $H\gg O$. The traditional Plummer function yields systematically overestimated slopes. Accurate deconvolution requires a priori knowledge of the true parameters, creating a fundamental circular problem whose only robust solution is obtaining sufficiently high angular resolution. Current far-infrared observations typically lack sufficient resolution, and some previously reported filament properties may require reinterpretation.

We review the role of primordial black holes (PBHs) for illuminating the dark ages of the cosmological evolution and as dark matter (DM) candidates. We elucidate the role of phase transitions for primordial black hole formation in the early Universe and focus our attention on the cosmological QCD phase transition within a recent microscopical model. We explore the impact of physics beyond the Standard Model (SM) on the cosmic equation of state and the probability distribution for the formation of PBHs which serve as candidates for DM and contribute to present-day binary black-hole merger events.

We investigate the impact of a late-time transition in the standardized absolute magnitude $M$ on the best-fit values of cosmological parameters using the Pantheon+ dataset. Extending previous analyses which focused on flat $\Lambda$CDM, we examine this transition within flat $\Lambda$CDM, wCDM, and CPL cosmologies, as well as a model-independent cosmographic expansion, employing both frequentist ($\chi^2$ minimization with \textit{AIC}/\textit{BIC}) and Bayesian (MCMC and Nested Sampling) inference frameworks. We confirm that the data consistently favor a step in absolute magnitude of $\Delta M \simeq 0.19~\mathrm{mag}$ at a characteristic distance of $d_{\mathrm{crit}} \approx 20~\mathrm{Mpc}$. The inclusion of this transition leads to a statistically significant improvement in the quality of fit and has a distinct impact on parameter estimation: it induces a systematic increase in the inferred Hubble constant of approximately $2\%$ across all tested models. In contrast, we find that the dynamical parameters governing the background expansion, including the matter density $\Omega_m$ and the dark energy equation of state ($w_0, w_a$), remain stable and largely unaffected. These results indicate that the $20~\mathrm{Mpc}$ feature acts primarily as a low-redshift calibration shift rather than a modification of the late-time expansion history.

Magnetic fields play a fundamental role in the formation of protostellar winds. In the magneto-centrifugal models, poloidal magnetic fields launch winds from accretion disks, and fast-rotating gas twists the fields into toroidal geometry that collimates and accelerates winds through magnetic hoop stress. However, toroidal fields in protostellar winds remain observationally unresolved. Here we report polarization observations of carbon monoxide emission toward the NGC1333 IRAS 4A protostellar outflow. The inferred magnetic fields are perpendicular to the outflow axis and aligned with the rotational structure of the outflow, indicating toroidal fields with strengths of a few milligauss, sufficient to collimate and accelerate the outflow at several hundred astronomical units from the protostar. A linear correlation is found between the curl of plane-of-the-sky magnetic field and the line-of-sight electric current density. Our analysis provides better constraints on ion-electron drift velocity in protostellar outflows and supports rotating outflows driven by the magneto-centrifugal mechanism.

Context. Be stars are presently viewed as B-type stars surrounded by a disc fueled by the star itself during episodicexcretion events. The origin of these events are poorly this http URL. This article aims to determine whether or not surface equatorial Kelvin waves can be unstable and therefore canplay a role in the triggering of the Be this http URL. We first derive an analytical expression for gravito-inertial modes in the shallow-water framework. Then, weinvestigate numerically the evolution of equatorial Kelvin modes as system parameters vary. The study is extended tothick-layer configurations with a constant density fluid. We then analyze the stability of these modes under differentialrotation and viscous this http URL. We show that equatorial Kelvin waves still exist in a spherical shell of finite thickness, but that their equatorialconfinement is weaker. At low azimuthal wavenumbers, Kelvin waves are in the inertial waves frequency band and thusget specificities of inertial waves like shear layers associated with singularities of the Poincaré equation. These shearlayers are new dissipative structures for Kelvin waves. When a radial (shellular) differential rotation is imposed, we showthat equatorial Kelvin waves can be destabilised provided that differential rotation and viscosity are in an appropriaterange. The non-monotonic behaviour of the growth rate of the instability is traced back to the rise of a critical layerwhere the fluid azimuthal velocity equals the phase speed of the surface this http URL. This study provides new insights into the behavior of equatorial Kelvin waves in astrophysics, particularlyin rapidly rotating stars. The results reinforce the idea that gravito-inertial waves, and more specifically the equatorialKelvin waves, can be unstable and thus be key parts in the mechanisms leading to the Be phenomenon.

We present a new metal enrichment analysis of a unique galaxy group IC 1262 using archival Chandra and GMRT observations, focusing on metal transport via radio jet, sloshing cold fronts, and shock front. This group shows two sloshing cold fronts along the east and north-west direction which is nearly orthogonal to the north - south orientated radio jet. We report discontinuities in the metallicity at the location of previously detected cold fronts, a more prominent one towards the eastern direction. In addition, the gas inside the cold fronts is 45$\pm$8 per cent more enriched than the gas outside the cold front, suggesting the role of sloshing in transporting metals through the IGrM. We also confirm the presence of a previously reported shock front with higher significance and with greater details. Across this shock, we detect a significant metallicity drop from 0.45$\pm$0.05 $Z_{\odot}$ to 0.22$\pm$0.04 $Z_{\odot}$, located at a projected distance of 78$\pm$2 kpc in the southern direction. The shock could potentially account for the region of gas enrichment seen in the abundance map and profile, which could be the result of a non-Maxwellian electron distribution in its vicinity. This should be considered a contributing factor rather than the sole cause of the observed discontinuity in the abundance. Furthermore, our spectral analysis reveals two temperature X-ray gas preferentially aligned with the radio-jet axis, indicating a possible influence of radio AGN activity on the surrounding gas.

In this study, we investigate the dynamical evolution of a massive binary black hole (MBHB) in the Leo I dwarf spheroidal galaxy model and examine how dark matter along with stellar matter's gravitational interactions influence its long-term behavior. Using high-resolution direct N-body simulations, we follow the orbital evolution of the binary within a realistic model of the Leo I stellar and dark matter distribution. We found that the binary separation decreases from an initial 300-parsec orbit to roughly 1 parsec over a period of about 2 Gyr, primarily driven by dynamical friction and stellar hardening. The orbital evolution then stalls at this scale, illustrating the well-known final parsec problem. During this phase, the binary also develops increasing orbital eccentricity and produces a modest redistribution of the inner mass profiles in some cases. We then further estimate the final stage of the system's evolution using gravitational-wave emission models and find that the binary is unlikely to merge within a Hubble time. The prolonged dynamical friction phase appears to be related to the low stellar and dark matter densities in Leo I. These results suggest that massive binary black holes in dwarf spheroidal galaxies such as Leo I will not contribute to the gravitational-waves detectable from LISA even if dark matter is considered.

We present PHLEGETHON, a fully compressible, Eulerian magnetohydrodynamic (MHD) code designed for multidimensional simulations in stellar astrophysics. The code uses a time-explicit, second-order, finite-volume method optimized to model a wide range of dynamical processes in stars, from very low-Mach-number turbulent convection in the cores of massive stars to supersonic flows in subsurface convection zones. PHLEGETHON employs low-dissipation Riemann solvers and a well-balanced method to accurately capture slow flows arising from strongly stratified media. The induction equation is solved using a staggered constrained-transport method to ensure divergence-free evolution of the magnetic field. The MHD equations are coupled to arbitrary nuclear reaction networks solved in a time-implicit approach, together with super-time-stepping for efficient treatment of thermal diffusion. Equations of state appropriate for stellar plasmas are available, accounting for partial ionization, electron degeneracy, and electron-positron pair production. The code is implemented in a compact and user-friendly manner, and it scales to tens of thousands of CPU cores using MPI-based domain decomposition. We perform several verification tests to demonstrate the accuracy and versatility of the code, and present simulations of magnetoconvection in a core-collapse supernova progenitor star. The rich variety of physical effects and numerical methods implemented in PHLEGETHON enables the code to model diverse multidimensional processes that play a crucial role in stellar-interior dynamics, such as reactive convection, convective boundary mixing, internal-wave excitation, and magnetic-field amplification mechanisms. Within a single framework, these phenomena can be investigated across a wide range of stellar evolutionary stages, from main-sequence stars to supernova progenitors. PHLEGETHON is publicly accessible online.

Core-collapse supernovae (CCSNe) in luminous infrared galaxies (LIRGs) can have extreme line-of-sight host galaxy dust extinctions, which leads to a large fraction of the events remaining undetected by optical and infrared surveys. This population of undetected CCSNe is important to constrain in order to determine the cosmic CCSN rates. Our aim is to confirm and refine our estimates for the undetectable fraction of CCSNe in LIRGs in the local Universe. Our study is based on the near-infrared K-band multi-epoch SUNBIRD survey monitoring dataset of a sample of nine LIRGs using the Gemini-South telescope with the multi-conjugate GSAOI/GeMS laser guide star adaptive optics system. We determined the limiting magnitudes for CCSN detection for each epoch in our dataset with artificial supernova injection and image subtraction methods. Subsequently, we used a Monte Carlo method to determine the combined effects of limiting magnitudes, survey cadence, CCSN subtype distribution, and their light curve evolution diversity. The intrinsic CCSN rates of the sample galaxies were estimated based on detailed modelling of their spectral energy distribution. Finally, we combined the resulting CCSN detection probabilities with the intrinsic CCSN rates for the dataset, and compared that against the real CCSN detections over the survey period. Based on our GSAOI/GeMS dataset, assuming optical or near-infrared example surveys with capabilities to detect CCSNe in local LIRGs with host extinctions of $A_V =$ 3 or 16 mag, respectively, the resulting total undetectable fractions are $86.0^{+4.7}_{-5.9}$ % and $53.6^{+15.6}_{-19.6}$ %. When folding in the results from our previous near-infrared adaptive optics assisted LIRG monitoring dataset, the corresponding total undetectable fractions are $88.3^{+2.6}_{-3.2}$ % and $61.4^{+8.5}_{-10.6}$ %, respectively.

Light curves from binary systems containing white dwarfs with neutron star or stellar-mass black hole companions (WD+NS and WD+BH) with edge-on orbital planes potentially show self-lensing/eclipsing signals. Here, we evaluate the properties and detectability of these signals in the NASA's Transiting Exoplanet Survey Satellite (TESS), and the Nancy Grace Roman Space Telescope (Roman) observations. WD+NS systems with orbital periods $T\lesssim25~$days mostly have considerable finite-source sizes with the normalized source radii $\rho_{\star}\gtrsim1$. WD+BH systems with $T\gtrsim3$ days have $\rho_{\star}\lesssim1$, and $\rho_{\star}\sim0.01$ for BHs with a few tens solar-mass. Our analytical calculations show the probabilities of occurring self-lensing signals in WD+NS and WD+BH systems are $\sim10^{-3},~10^{-2}$, and maximize for systems with low-mass WDs revolving massive NSs/BHs. We simulate their light curves and generate synthetic data for them by applying the observing protocols of these two satellites. We assume self-lensing signals are detectable if (i) $1\leq T\leq T_{\rm{obs}}$ (where $T_{\rm{obs}}=62~\rm{and}~27.4$ days are the Roman and TESS continuous observing windows), (ii) $\rm{SNR}\ge3,~6$, their signals are (iii) deeper than twice the photometric error, and (iv) covered by at least one datum. Systems with detectable self-lensing signals in the TESS and Roman observations on average have small inclination angles $i\lesssim0.2^{\circ}$, with the orbital periods $\sim6,~19~$days, and their signals last $\sim[6,~30]~\rm{minutes}$. The TESS and Roman efficiencies for detecting these signals are $\sim2-6\times10^{-4}$ and $\sim2-12\times10^{-10}$. Although detecting these self-lensing signals by Roman is impossible, the TESS telescope potentially manifests at least one self-lensing signal due to these binary systems, if $8\%,~\rm{and}~3\%$ of WDs have NS and BH companions.

Within the volume-limited subsample at $z<0.06$ of the Zwicky Transient Facility (ZTF) DR2 sample, we confirm a statistically significant excess of Type Ia supernovae (SNe Ia) at $z \simeq 0.02$-$0.04$, previously reported but not explained by survey selection effects. Forward simulations assuming a uniform volumetric SN Ia rate and realistic ZTF detection efficiencies fail to reproduce the feature at the $5$-$7\sigma$ level. We also detect an excess in the rates compared to our survey simulations at $z \simeq 0.08$ and $0.14$, albeit at smaller significance. To investigate the origin of these inhomogeneities, we compare the observed SN distribution to constrained reconstructions of the local matter density field from the Manticore project, based on Bayesian forward modelling of the 2M++ galaxy catalogue. While SN overdensities are spatially associated with prominent nearby structures such as the Perseus, Coma, and Hercules superclusters, the amplitude of the SN excesses significantly exceeds that expected from matter overdensities alone. By reconstructing a redshift-dependent volumetric SN Ia rate, we find that local enhancements can reach factors of two to five within specific clusters, while the sample-averaged rate remains consistent with previous low-redshift measurements. These results indicate that the SN Ia rate is not a linear tracer of the underlying matter density and suggest a strong environmental dependence in dense structures. We discuss possible physical origins and highlight the implications for low-redshift SN cosmology, including correlated peculiar velocities and additional covariance beyond standard linear corrections.

Here we present a survey aimed at detecting hydroacoustic signals from fireballs using the six hydrophone stations operated as part of the Comprehensive Test Ban Treaty Organisation (CTBTO) International Monitoring System. We identified 30 fireballs where propagation paths to stations exist. These included high energy fireballs (E $\geq$ 5 kT), those which occurred over favorable locations for coupling into the deep ocean as well as a selection of bolides close to CTBTO hydrophone stations. The largest of these impactors were $>$ 5 meters in diameter. From theoretical and empirical considerations we show that direct hydroacoustic shock transmission is the most likely source mechanism, though large meteorites impacting the ocean surface from a fireball might be detectable in extreme cases. We find one possible instance of a fireball occurring on Sep 2, 2003 off the coast of Alaska, where a linked hydroacoustic signal with the expected timing and backazimuth is detected. However, given the size of our survey and the random background rate of signals, this detection is statistically weak. We conclude that hydroacoustic detection in the SOFAR channel of fireballs is very rare. Using our chosen set of signal processing parameters, assuming direct path H-phase signals, adopting a signal celerity range of 1.42-1.55 km/s we find no unambigous detections in 53 station-fireball pairs. Based on SOFAR-equivalent yields derived assuming the minimum detectable amplitude signal family association is representative of the noise background in our survey we estimate a conditional upper limit for fireball coupling efficiency of order 10$^{-10}$. A single well recorded airplane impact provides an empirical estimate for the energy coupling of surface ocean impacts to the SOFAR channel of 10$^{-4}$ for high velocity surface impacts.

While photoionization rates of Ganymede's O$_2$ dominated atmosphere are well constrained, the contribution of electron-impact ionization is rather uncertain. Previous quantitative estimates have relied on assumptions about densities and energy distributions of precipitating electrons, or on rare spacecraft measurements that cannot be unambiguously mapped to the regions of ionization. In this study, we present a novel approach to quantify electron-impact ionization rates directly through OI 1356 Å emission brightness observations. The analysis of measured cross sections reveals that the ionization-to-excitation ratio is limited to 10-60 over all electron energies, reducing the uncertainty of estimating ionization rates to a factor less than 6. We apply this method to Juno UVS observations of Ganymede's aurora. We find that the OI 1356 Å brightness of the auroral ovals is well described by 3-5° latitude wide Gaussian distributions centered on the open-closed field line boundary, with an average peak of 120 R. The average brightness outside the ovals in the polar and equatorial background regions is ~8 R. From these observations, we derive a global map of electron-impact ionization rates, which are at least an order of magnitude higher than photoionization rates. The estimated total global ionization rate is 1.3-7.6$\times$10$^{26}$ s$^{-1}$, with average column rates of ~5$\times$10$^{9}$ cm$^{-2}$s$^{-1}$ in the ovals and ~3$\times$10$^{8}$ cm$^{-2}$s$^{-1}$ in the background regions. Comparison of radio occultation measurements with predicted electron densities indicates that transport processes are the dominant loss mechanism in Ganymede's ionosphere. The rate of ionospheric outflow of O$_2^+$ is 0.1-2$\times$10$^{26}$ s$^{-1}$ or 0.5-11 kg s$^{-1}$, indicating 0.03-0.5 cm Myr$^{-1}$ erosion of Ganymede's surface ice.

Generally, large solar flares accelerate electrons to high energies more efficiently than microflares. However, some microflares, known as hard microflares (HMFs), also produce high-energy electrons, as indicated by their flat hard X-ray (HXR) spectra. These events are typically associated with footpoints located in or at the edge of sunspots. The mechanisms behind this efficient acceleration, and their connection to solar energetic electrons (SEEs), remain unclear. We compare, for the first time, HXR spectra of HMFs with in-situ electron spectra of associated SEEs using Solar Orbiter STIX and EPD observations. This provides insight into acceleration processes and the transport of high-energy electrons into interplanetary space. We identify eight HMFs observed jointly by Solar Orbiter and Earth-based instruments that are associated with SEEs, confirmed through timing and magnetic connectivity analysis. Each event is studied using HXR spectroscopy, SEE velocity-dispersion analysis, and in-situ electron spectral analysis. Seven of eight events show consistent timing between flare HXR emission and inferred electron injection, as well as good agreement with magnetic connectivity estimates. The known correlation between HXR photon and in-situ electron spectral indices extends to HMFs, which occupy the hard end of the distribution, even compared to larger flares. We conclude that HMFs produce prompt SEEs with hard spectra, demonstrating efficient electron acceleration without requiring large flare energy release. Their magnetic configuration, involving open field lines from the sunspot, suggests they may be an important contributor to filling the heliosphere with energetic particles.

The PLAnetary Transits and Oscillations of stars (PLATO) mission is expected to launch in January 2027. A total of 8\% of its data rate will be dedicated to complementary science targets selected from approved Guest Observer proposals. We seek to provide an open-source catalogue of variable stars in PLATO's first long-term observing field, LOPS2. We want to use existing observations from the Transiting Exoplanet Survey Satellite (TESS), which has observed many stars in LOPS2. We classified 38 million calibrated aperture light curves from the TESS-Gaia Light Curve pipeline (TGLC, $G\lesssim17$) for 6 million unique sources in LOPS2 with two machine learning frameworks -- a deep neural network and a feature-based gradient-boosted decision-tree ensemble. We combined their predictions to create this first version of the LOPS2 variability catalogue, performed manual vetting of a sub-sample classified light curves, and a statistical analysis of the results to validate our methodology and to assess the variability properties and parameters of the stars in the catalogue. Our classification resulted in the identification of approximately 72% of the light curves having dominant instrument- or pipeline-induced signal, with the remaining 28% representing 3.6 million individual candidate variable stars, including pulsating, rotating, and eclipsing stars. Candidate pulsators exhibit varied behaviour in terms of their frequencies, amplitudes, rotation, and fundamental parameters. To ensure purity of the samples, filtering on colour, luminosity, the dominant frequency and its amplitude, and presence of close neighbours is helpful. We provide the first version of our PLATO LOPS2 variability catalogue to the community for further study and scrutiny. It is to date one of the largest catalogues of variable stars from an automated classification pipeline.

The temporal and spatial coincidence between the gravitational wave (GW) event GW190425 and the fast radio burst (FRB) event FRB 20190425A raises the intriguing possibility of a physical connection between the two. The widely discussed possibility invoking the collapse of a supermassive neutron star as the merger product suffers the inconsistency between the model prediction and the measured inclination angle of the system. Here, we propose a novel physical mechanism to account for the association. We envisage a magnetar located at about 2.5 light hours away from the binary neutron star merger site. The kiloherz GWs generated by the merger are converted into kiloherz electromagnetic (EM) radiation via the Gertsenshtein-Zeldovich (GZ) effect near the magnetar. Subsequent inverse Compton scattering off the kilohertz EM waves by relativistic particles generates the observed gigahertz FRB emission. Our calculation reveals that, with appropriate parameter choices, the properties of FRB 20190425A can be reproduced.

Persistent radio sources (PRSs) are (sub-)parsec-scale compact non-thermal continuum sources associated with some repeating fast radio bursts (FRBs). Their nature is debated, but their properties provide insights into the FRB environment and progenitors. We measure the spectrum of the recently confirmed PRS associated with FRB 20190417A. Spectral features such as the self-absorption and cooling break can be used to constrain the age and size of PRSs and test theoretical models. We present observations made with the 1.26 GHz upgraded Giant Metrewave Radio Telescope (uGMRT) and observations from the 6 GHz Karl Jansky Very Large Array (VLA). With complementary archival data and the LOw Frequency ARray Two Meter Sky Survey (LoTSS), we characterise the spectrum of the PRS between 144 MHz and 6 GHz. The spectrum follows a power-law behaviour at gigahertz frequencies. The source is not detected at 144 MHz down to a $2\sigma=170 \; {\rm \mu Jy}$ sensitivity. We modelled the spectrum with a broken power law, obtaining a spectral index $\alpha = 0.20 \pm 0.05$ between 1-6 GHz. We placed a lower limit on the turn-over frequency of $> 370$ MHz ($95\%$ confidence). The flat spectrum and low-frequency turn-over of the target are consistent with the spectral properties predicted for magneto-ionic nebulae, inflated behind the supernova ejecta by a flaring young magnetar. Considering the multi-zone magnetar wind nebula scenario, we estimate an age of $t< 250$ yr and a radius of $R< 0.4$ pc for the target, which would thus be slightly older than the PRSs associated with FRB 20121102A and FRB 20190520B.

Context. Studying the kinematic and dynamic evolution of fast eruptive events from the middle to high solar corona is a primary objective of the Metis coronagraph on Solar Orbiter. During perihelion, Metis acquires visible light images at a 20s cadence, reaching a spatial resolution of around 2000 km at 0.28 au. This enables capturing coronal mass ejections (CMEs) and transient structures with unprecedented spatial and temporal resolution. Aims. On October 8-9, 2022, an extensive CME was observed by Metis at 0.3 au (resolution: 4.4 $10^{3}$ km/pixel). We aim to exploit this high resolution to resolve multiple substructures within the CME front, revealing plasma elements with distinct trajectories and speeds to provide a detailed kinematic characterization of the eruption. Methods. A normalization-based running difference algorithm was applied to enhance the complex morphology. Height-time diagrams were used to estimate propagation speeds and frequency variations. A 3D flux rope reconstruction, combined with multi-spacecraft coronagraphs and disk imagers, enabled tracking the CME from its initiation in the lower corona to approximately 5 solar radii. Joint observations with Solar Orbiter EUI-FSI provided insights into the eruption's onset, while Metis captured its development into the middle corona. Results. Metis observations resolve the CME's fine structure and internal plasma motions. The detection of circular, fast-propagating wavefronts (500 km/s, 3 minute period) at the western flank suggests wave excitation and magnetic reconfiguration processes. Multiple interpretations are proposed for these coronal wave trains, including quasi periodic propagating fast modes, offering new insights into wave generation and energy transport in the solar corona.

Gaia has the potential to deliver several tens of new dormant black holes (BHs) with low-mass stellar companions (hereafter, Gaia BHs) in the upcoming fourth data release. Three Gaia BHs are already known, but their formation pathways remain uncertain. Here, we perform a large parametric study to explore the formation of Gaia BHs from isolated binary systems with the population-synthesis code SEVN and compare our models with the properties of the three already reported Gaia BHs. Specifically, we explore the impact of accretion efficiency, mass transfer stability, natal kicks, angular momentum transport, and core-collapse supernova prescriptions. We find that models in which stable mass transfer is highly non-conservative and angular momentum is lost as a wind from the donor surface (Jeans mode) maximize the probability of forming dormant systems that match the properties of the observed Gaia BHs in terms of both orbital period and eccentricity, because such assumptions prevent the initial orbit from shrinking too much when the BH progenitor fills its Roche lobe. If we allow for common-envelope evolution, we find that models with common-envelope ejection efficiency $\alpha{} < 1$ predict dormant systems with orbital periods that are too short compared to the observed Gaia BHs. The eccentricity of the observed Gaia BHs, when combined with information about orbital period and BH mass, favors relatively large natal kicks, similar to those inferred from Galactic neutron stars. Finally, models in which the natal kicks are low - e.g. because they are modulated by fallback - result in the formation of a large population of dormant BHs with long orbital periods ($P_{\rm orb}>10^4$ days) and low eccentricity, which will be tested soon by the fourth Gaia data release.

We analyze JWST spectroscopy of the northwest filament of NGC7023, where the relatively soft radiation field results in a photodissociation region with an extended atomic hydrogen region, and strongly pronounced variations of the carbonaceous emission band profiles. We focus on the 16.4 and 17.4 um bands and their relation to the main bands at 3.3, 3.4, 5.2, 5.7, 6.2, 7.7, 8.6, 11.3, and 12.7 um, and aim to identify which bands and sub-features originate from co-spatial emission carriers. We apply a PAHFIT spectral decomposition to measure the emission bands and their individual sub-components, and produce maps that spatially resolve the main dissociation front (DF1). Nearly all emission maps peak at DF1, while the relative intensity in the atomic hydrogen region (ATM) varies strongly. We classify the features into spatial distribution types based on the intensity ratio in ATM relative to DF1. Most bands are of type I (low ATM/DF1; 3.3, 3.4, 5.2, 5.7, 11.3 um) or II (medium ATM/DF1; 16.2, 7.7, 8.6, 12.7, 16.4 um), while only few are of type III (high ATM/DF1; 11.0, 17.4 um). A breakdown of the 5.7, 7.7, 11.3 and 12.7 um bands into blue and red sub-components reveals that contributions on the blue side are of type III, while those on the red side are of type I or II. These strongly differing spatial distributions reveal that at least two different populations contribute to the 16-18 um range, and that these populations are also connected to the profiles of the 5.7, 7.7, 11.3, and 12.7 um bands. The maps further indicate a continued evolution of these profiles toward the central cavity of NGC7023, where fullerene emission (C60) was previously detected. We speculate that the population of emission carriers could be in an intermediate photochemical evolution stage that precedes fullerene formation.

While the connection between massive stars and supernova explosions is well established observationally, the link between massive stars and black hole formation remains elusive. Some massive stars may collapse directly to black holes without a successful supernova, and may therefore appear as disappearing stars. We investigate the expected photometric properties of such black hole progenitors by combining detailed single and binary stellar evolution models with physically motivated prescriptions linking pre-collapse core structure to explosion or direct collapse outcome, together with stellar atmosphere calculations, producing synthetic photometry across standard ultraviolet to infrared filters. Weighting by an initial mass function and empirical binary distributions, we predict both the observable distribution of progenitor brightness and colour and the rate of direct-collapse events, which we estimate to be about 0.4 per century for a galaxy forming stars at 1 Msun/yr. We find that black hole progenitors are predominantly hot and blue at the pre-collapse stage, with many in Wolf-Rayet phases and luminous in the ultraviolet, while only a minority are red supergiants. Consequently, searches that focus primarily on red supergiants are likely to miss a substantial fraction of direct-collapse events. Monitoring campaigns that include ultraviolet-sensitive observations of nearby star-forming galaxies therefore provide a promising route to detecting disappearing massive stars, offering a direct observational probe of black hole formation. Our results provide predictions to interpret such surveys and constrain the channels that lead to black hole formation.

The actinides, including thorium (Th), are the heaviest observable elements synthesized in the universe, holding clues to the extremes of the astrophysical and nuclear conditions of $r$-process sites. We present Th abundances based on high-resolution spectroscopy for 47 metal-poor stars, the largest homogeneously analyzed sample to date. The chemical evolution of Th exhibits a decrease in dispersion in [Th/H] and [Th/Fe] from $\sim$0.6 dex at the lowest metallicities to $\sim$0.2 dex at higher metallicities. We also find that Th and the lanthanides Eu and Dy are co-produced remarkably well, with average [Th/Eu]$\sim0.0$ across $-3.0 \lesssim$ [Fe/H] $\lesssim -1.5$, as well as across stars with $0.0\lesssim$ [Eu/Fe] $\lesssim2.5$. Even so, the absolute range of $\log\epsilon$(Th/Eu) is 1.02 dex, with an observed standard deviation of $\pm0.20$ dex and an intrinsic standard deviation of $\pm0.11$ dex at the lowest metallicities. We infer that $68\%$ of $r$-process events have $\log\epsilon$(Th/Eu) yields that only vary within a factor of $\pm1.3$ or $\pm30\%$, while $5\%$ of $r$-process events have $\log\epsilon$(Th/Eu) yields that vary by factors $>3.3$ approaching $\sim$10. This serves as a strong constraint for the nuclear and astrophysical models of $r$-process sites, and suggests that achieving an $r$-process site that is both prompt and produces a robust $\log\epsilon$(Th/Eu) ratio is a challenge for current models.

VLBI measurements of the sizes of compact extragalactic radio sources, jetted active galactic nuclei, provide data for probing the angular size--redshift relation, offering a complementary test to other distance--redshift methods. We analyse a significantly expanded dataset to reassess the angular size--redshift relation and its potential for constraining cosmological model parameters, focusing on the matter density parameter $\Omega_{\mathrm{m}}$ in a flat $\Lambda$ Cold Dark Matter Universe. This is the first major update of the compact-source angular size test in the past quarter of a century, using a dataset an order of magnitude larger than in previous studies. MCMC analysis on real data and on multiple mock catalogues to evaluate parameter constraints in the presence of observational scatter. In addition, we conducted a test with 100 randomized catalogues created by shuffling redshifts while preserving other observables. We also explored how astrophysical parameters depend on fixed cosmological models with different $\Omega_{\mathrm{m}}$ values. The randomization test showed that the posterior distributions from randomized data do not overlap with those from real observations, confirming that the measured angular size--redshift relation is physically meaningful. The astrophysical model parameter that describes the redshift dependence of the source angular size exhibits degeneracy with $\Omega_{\mathrm{m}}$. Simulated mock catalogues indicate that the method is able to constrain $\Omega_{\mathrm{m}}$ if the data scatter is below $\sim20\%$. Scaling estimates suggest that high-quality data of samples of several thousands to $\sim 100\,000$ sources, a standardisation calibration approach, and/or refining sample selection criteria are needed to fully exploit the potential of the angular size--redshift test with this type of objects (abridged).

We investigate how strongly late-time inferences about DE dynamics depend on the functional prior used to represent the expansion history. Using identical late-time combinations of CC, DESI BAO measurements, the Pantheon+ SN1a sample, and the H0DN prior, we compare a node-based reconstruction of the reduced Hubble function $E(z)$ with a representative family of smooth low-dimensional DE EoS parametrizations, including CPL. Over the redshift range constrained by the data, both approaches yield consistent $H(z)$, and, in the absence of H0DN, compatible values of $H_0$. However, a clear method dependence emerges at intermediate redshift ($z\sim1.7$): the reconstruction favors stronger deceleration, $q_{\rm Rec}(1.7)\simeq0.56-0.61$, whereas the smooth parametrizations cluster at $q(1.7)\simeq0.32-0.40$, implying a persistent $\sim2-3\sigma$ discrepancy across dataset combinations and parametrizations. For the EoS-based parametrizations, whose effective DE densities remain positive by construction, the preferred $w_{\rm DE}(1.7)<-1$ values correspond to NECB-violating (phantom-like) behaviour, but this is a less robust discriminator as $w_{\rm DE}$ becomes ill-conditioned as $\rho_{\rm DE}\to0$. In the effective-fluid mapping, the reconstruction accommodates the same late-time kinematical preference through a rapid descent of $\rho_{\rm DE}(z)$ toward very small values and a sign change, whereas the EoS-based parametrizations absorb it through smoother, and in several cases NECB-violating, evolution over $z\sim1-2$. Although the reconstruction improves the best-fit likelihood, especially with H0DN, Bayesian evidence continues to favor the simpler parametric descriptions. Our results isolate $z\sim1.5-2$ as the key window in which EoS-based DE parametrizations can compress localized kinematic structure and associated features of DE that are still permitted by current late-time data.

Supersonic winds from massive stars carry great amounts of kinetic power and modify the surrounding interstellar medium. Through this interaction a stellar bubble is formed. Theoretical studies and recent observations suggest that the winds of massive stars could be sources of Galactic cosmic rays. The first detection of synchrotron emission from the bubble of a single star was reported, indicating the presence of relativistic electrons. Studying the non-thermal emission from a single massive star can help to better understand the acceleration of particles taking place in massive star clusters. WR 102 is the perfect case of study. In this work, we present the first high-energy model for the bubble of WR 102: G2.4+1.4. We aim at fitting the radio data and predicting gamma-ray emission. We assume that both electrons and protons are accelerated at the wind shock. We applied a classical model for the stellar bubble and adopted a one-zone model for estimating the radiation produced by the relativistic particles near the acceleration region. Additionally, we computed the expected emission from the protons that diffuse to the outer regions of the bubble. Also, we estimated the leptonic and hadronic contributions expected from cosmic rays. We fitted the observations considering that 3% of the wind kinetic power goes into relativistic electrons, and a magnetic field of 250 $\mu$G. The dominant component at high energies is produced by locally accelerated protons reaching the shell. Protons might reach PeV energies in the wind bubble, but the predicted gamma-ray flux is too low to be detectable.

Little Red Dots (LRDs) are a population of compact, red sources that have emerged as one of the most puzzling findings of JWST. Variability provides a direct probe of their central engines. Here we present the first joint spectroscopic and photometric time-domain study of LRDs undertaken with the JWST TWINKLE slitless spectroscopy program. Surveying the FRESCO GOODS-North legacy field, TWINKLE monitors a complete, H$\alpha$-flux-limited sample of 18 LRDs at z = 3.9-6.8, achieving a rest-frame baseline of $\sim$140-220 days. We detect no variability in photometry, H$\alpha$ line flux, or line shape across the sample. If LRDs resembled AGN in reverberation mapping samples -- the foundation for black hole mass calibrations and luminosity scaling relations -- we would expect >10 sources to show measurable fluctuations. Observing none implies a 5.9$\sigma$ deficit. The non-detections hold across all broad H$\alpha$ emitters within TWINKLE's field of view -- the 18 V-shaped LRDs as well as 9 non-LRDs. Comparison with simulated light curves disfavors sub-Eddington accretion and is instead consistent with super-Eddington accretion, other mechanisms that suppress variability, or perhaps no AGN whatsoever. If LRDs do harbor black holes, calibrations derived from sub-Eddington systems may not apply, thereby explaining JWST's apparently "overmassive" black holes. These observations provide unique constraints on the physics of one of the most enigmatic populations discovered by JWST.

Scalar-induced gravitational waves (SIGWs) have recently attracted considerable interest, both as a possible explanation for the nanohertz signal reported by the Pulsar Timing Array (PTA) collaboration and for their connection with primordial black hole (PBH) physics. In addition to SIGWs, scalar-tensor-induced gravitational waves (STGWs) have emerged as a promising cosmological source of the stochastic gravitational wave background (SGWB). In this paper, we compute the STGWs generated during a generic matter-dominated (MD) era, as well as during an early matter-dominated (eMD) epoch followed by a sudden transition to the standard radiation-dominated (RD) stage, working in the Poisson gauge. We find that, in a purely MD age, the corresponding energy density rapidly dilutes, whereas in the presence of an eMD phase it remains non-vanishing due to the short duration of the eMD period. We then investigate whether the STGW signal could provide a dominant contribution to the $\texttt{NANOGrav 15-year}$ dataset and we forecast the prospects for its detection with future observations by the Square Kilometre Array (SKA). In particular, we consider STGWs generated during both eMD and RD eras, including their linear-order contributions. Our results show that the GWs induced by scalar-tensor mixing constitute a viable target for future, more sensitive detections of the SGWB.

Although Uranus and Neptune are commonly classified as ice giants, their exact compositions remain poorly constrained. Recent studies of outer Solar System bodies challenge the traditional view that these planets are primarily ice-dominated, suggesting that refractory material plays a more significant role. Determining the proportions of ice and rock within Uranus and Neptune is essential for understanding their formation and the evolutionary history of the Solar System. In this work we computed interior structure models for both planets and explored, within a Bayesian framework, the range of compositions that satisfy the available observational constraints. We quantified the resulting ice and rock fractions and analyzed their impact on the inferred internal structure. Our results suggest that the envelopes of both Uranus and Neptune are systematically enriched in refractory material, with median rock fractions of approximately 60% within the heavy-element component, similar to Pluto, Kuiper belt objects, and comets. In contrast, the deep interiors of the two planets exhibit distinct compositions: Neptune is best fit by relatively rock-rich mantles (median rock fraction of ~ 55%), whereas Uranus is inferred to have more ice-rich mantles (median rock fraction of ~ 41%), consistent with a more strongly stratified structure. These results point to compositional differences between Uranus and Neptune that may reflect divergent formation and evolutionary pathways.

Reddened Type-1 quasars offer a unique window into the structure and evolution of active galactic nuclei (AGN), yet their physical origin and the source of their reddening remain uncertain. Optical surveys often miss these dust-obscured objects, resulting in an incomplete view of the quasar population. In this work, we construct a sample of 6,600 Type-1 quasars at redshifts $0.5 \leq z \leq 1.2$ by combining deep optical imaging from HSC with mid-infrared photometry from WISE, enabling a more complete selection that is not biased against reddened objects. We perform detailed SED modeling using the CIGALE code, enhanced by synthetic photometry derived from SDSS spectra to better constrain the optical continuum. We classify quasars into blue and reddened Type-1 populations based on their continuum slopes and compare their SEDs and emission line properties. As expected from this definition, reddened Type-1 AGN show higher dust extinction, with a median $A_V = 0.60^{+0.32}_{-0.19}$ mag, compared to $A_V = 0.06^{+0.10}_{-0.03}$ mag for blue objects. But they also exhibit smaller torus half-opening angles, with a median of $25.7^{+10.1}_{-8.7}$ deg, compared to $33.3^{+11.1}_{-5.9}$ deg for blue objects. While such extinction could arise on either galaxy or nuclear scales, the systematically stronger narrow-line equivalent widths and weaker Balmer broad lines in reddened Type-1s indicate that the obscuration acts on nuclear scales, likely from dust concentrated near the polar axis. We discuss the possibility that these structural differences may be linked to a sub-pc outflow, that carries dusty gas into the polar region and evacuates the torus region.

The issue and proof of Gurzadyan theorem are presented concisely, avoiding tedious and unnecessary calculations that would mask what is essential. The goal is to provide a good mathematical and physical understanding of the theorem, making you want to learn more about its use in cosmology.

This paper investigates a phased array ground station capable of tracking multiple sources, multi-beamforming, electronic steering, easy scaling, and low cost. The project will develop a 20-pair dual-polarized yagi-UHF/VHF phased array with a pseudo-random layout, comparing parameters of random and uniform distributions. We will present several analyses: general analysis for side lobes across both elevation and azimuth, analysis of scaling with number of elements ("element sweep"), electronic beam steering analysis, mechanical beam steering analysis, electro-mechanical beam steering analysis, array density analysis, and reception/transmission spectra analysis.

During the voyages that led him to discover the new continent bearing his name, Amerigo Vespucci made interesting astronomical observations of the southern sky. In the past, his data have been interpreted with criteria that do not follow Vespucci's indications, resulting in identifications that are not credible or even leading to the assertion that the data themselves are incomprehensible. However, it is possible to construct a coherent picture of all the information, arriving at an identification that is in some cases very probable, in other cases almost certain, of the stars described by Vespucci. Analysis of documents shows that he made good-quality measurements, but his incomplete knowledge of ancient texts prevented him from distinguishing the new stars from the already known ones, giving rise to a period of confusion in 16th century celestial cartography.

We study the geodesic motion and epicyclic oscillations of massive test particles around a static, spherically symmetric black hole (BH) solution of the Einstein--Skyrme (ES) theory in Anti-de Sitter (AdS) spacetime. The lapse function of this BH depends on the Skyrme coupling $\eta$, a charge-like parameter $Q$ inherited from the Skyrme term, and the cosmological constant $\Lambda<0$. We first map out the horizon structure and identify three regimes-non-extremal BH (NEBH), extremal BH (EBH), and naked BH (NBH)-showing that the NEBH $\to$ EBH $\to$ NBH transition is governed by $Q$ rather than $\eta$, which enters $f(r)$ only as a constant shift. We then derive the effective potential (EP), locate the innermost stable circular orbit (ISCO), and compute the radiative efficiency, finding that $\mathcal{E}_{\rm ISCO}>1$ in AdS renders the standard Novikov-Thorne formula negative. The corrected radial epicyclic frequency $\Omega_r$ reveals a distinctive AdS signature: $\nu_r$ grows at large $r$ and overtakes the orbital frequency $\nu_\phi$, causing the periastron precession frequency $\nu_p = \nu_\varphi - \nu_r$ to change sign-a feature absent in asymptotically flat geometries. Adopting the relativistic precession (RP) model for quasi-periodic oscillations (QPOs), we perform a Markov chain Monte Carlo (MCMC) analysis using twin-peak QPO data from XTE~J1550-564, GRO~J1655-40, Sgr~A$^*$, and M82~X-1. The posteriors converge to $Q\approx 0.6$ across all sources, with orbital radii near $r\approx 4.2\,M$ and masses consistent with independent estimates, demonstrating that the ES-AdS BH accommodates the observed frequency pairs within physically motivated parameter ranges.

A new radiation solver, RAPRAL (RAdiation Prediction based on RAy tracing and Line-by-line) implemented in C++, is developed for simulating high-temperature thermochemical nonequilibrium radiative processes. RAPRAL integrates detailed line-by-line spectral modeling with a ray-tracing solution of the radiative transfer equation, enabling accurate resolution of both spectral features and spatial radiation transport. The adopted methods and their implementation are described in detail. To assess the overall capability and accuracy of RAPRAL, we first focus on the computation of atomic and molecular bulk spectral coefficients. Through comparison with the established code in the literature, RAPRAL demonstrates its ability to accurately capture key spectral features across a wide range of conditions. Moreover, RAPRAL is applied to predict afterbody radiative heating in the Fire II flight experiment, based on a two-temperature, 11-species air flowfield. The results demonstrate that the present approach provides reliable predictions of radiative heat flux and effectively captures the dominant radiation mechanisms. Overall, the presented results demonstrate that RAPRAL is a robust tool for simulating radiative processes in hypersonic air flows, and future versions will extend its capabilities to include species relevant to planetary atmospheres.

Next-generation (XG) gravitational-wave detectors, such as Cosmic Explorer (CE) and the Einstein Telescope (ET), will observe compact binary coalescences at unprecedented rates and signal-to-noise ratios (SNRs). Accurate sky localization of these sources is crucial for several aspects of the science case of CE and ET. The localization of most binary black hole (BBH) signals, which will spend at most a few minutes within the XG detector's effective sensitivity band, will continue to rely primarily on timing triangulation across a network of detectors. A key design choice for triangulation is the baseline between instruments. We investigate how the baseline affects the localization capabilities of a two-detector CE network, analyzing both fixed-parameter injections and a realistic BBH population consistent with the latest GWTC-4 results. For detector-frame total masses up to $\sim\!100\,{\rm M}_\odot$, we find that baselines corresponding to light travel times of $8-11$ ms ($\sim\!2300-3300$ km) offer a reasonable compromise, producing predominantly unimodal or bimodal sky localizations suitable for electromagnetic follow-up and statistical host galaxy identification and galaxy cross-correlation studies. Shorter baselines significantly degrade localization, particularly for high SNR events. Crucially, we find that adding a third detector to the network eliminates localization multimodality for a substantial fraction of sources. A network with two CEs and LIGO-India provides unimodal posteriors for a good fraction of events, whereas two CEs plus ET would provide unimodal posteriors for essentially all of them. These considerations should be useful to inform the development of the XG detector network.

We present GreyRing, a new model for the post-merger signal in black-hole binary coalescences based on the greybody factor of the remnant. The model accurately reproduces the full frequency-domain ringdown signal of a large set of comparable-mass, aligned-spin numerical relativity waveforms, achieving mismatches of order ${\cal O}(10^{-6})$ for the dominant $(\ell,m)=(2,2)$ mode, and typically outperforming state-of-the-art time-domain models. Building on this model, we introduce a novel consistency test of strong gravity based on the greybody factor: the remnant mass and spin inferred from GreyRing can be compared with those obtained through standard black hole spectroscopy. This agnostic test relies exclusively on the post-merger signal and does not require the inclusion of overtones or the choice of very early ringdown starting times, combining the advantages of inspiral-merger-ringdown consistency tests and traditional black hole spectroscopy. We apply the test to GW250114 and find that the remnant mass and spin inferred from GreyRing are consistent with those measured from the full signal. Remarkably, the inferred parameters can be measured with a precision comparable to, or slightly better than, that achieved with standard black-hole spectroscopy. Our greybody-factor waveform model allows for new precision tests of strong gravity using the ringdown signal.

We introduce pyEFPEHM, a post-Newtonian (PN) inspiral waveform model for eccentric and spin-precessing compact binaries that includes higher-order modes and matter effects. Accurate and efficient waveform models capturing these effects are essential for probing compact-binary formation channels and exploiting current and future gravitational-wave (GW) observations. pyEFPEHM extends pyEFPE, significantly improving its physical content and accuracy. In particular, we show that above 2.5PN order the quasi-circular contributions to the orbital phasing dominate at each PN order, and incorporate all available higher-order quasi-circular PN corrections to the phasing, including adiabatic tidal effects. We generalize the multiple-scale analysis solution of the spin-precession equations, extending it to higher PN orders and including all available quasi-circular corrections. Finally, we add eccentric corrections up to 1PN order in the waveform amplitudes, including the GW multipoles $(l,|m|)=(2,2),(2,1),(2,0),(3,3),(3,2),(3,1),(3,0),(4,4),(4,2),(4,0)$. We validate pyEFPEHM against analytical waveform models and numerical relativity simulations, showing that it provides a robust and computationally efficient description of the inspiral, with good agreement across a broad region of parameter space and up to close to merger. The accuracy degrades in the late inspiral for systems with very unequal masses ($m_2/m_1 \lesssim 0.1$), significant spins aligned with the orbital angular momentum ($|\chi_\mathrm{eff}| \gtrsim 0.5$), and high eccentricities ($e \gtrsim 0.6$), where the PN expansion is expected to break down. pyEFPEHM represents a significant step toward physically complete and efficient waveform modeling of eccentric and precessing binaries, providing a foundation for future extensions including higher-order corrections, calibration to numerical relativity, and merger ringdown modeling.

We report the properties of electron distributions and turbulence during a Coronal Mass Ejection (CME) in April 2023 observed by Magnetospheric Multiscale (MMS). The CME exhibits a clear sheath and magnetic cloud (MC), and within the MC, the solar wind becomes sub-Alfvénic for two hours. We investigate plasma and turbulence properties of the sub-Alfvénic CME wind and compare them with those in the super-Alfvénic solar wind in the MC and CME sheath. Electrons within the sub-Alfvénic MC show significantly higher temperatures than those in the CME sheath and the super-Alfvénic MC, with their one-dimensional distributions revealing super-thermal tail and a depletion in electron populations between 15-50 eV. Within the CME sheath, isolated regions of electron heating are observed, where parallel energy flux is enhanced up to ~1 keV. Magnetic field fluctuations within the sub-Alfvénic MC interval exhibit negligible cross helicity and steeper-than-Kolmogorov scaling in the inertial range, with no clear spectral break. These fluctuations also show reduced intermittency at ion and sub-ion scales, emerging intermittency at electron scales, and weak magnetic compressibility. Together, these observations point to the presence of weak magnetohydrodynamic (MHD) turbulence within the sub-Alfvénic MC, resembling conditions commonly observed in planetary magnetospheres such as Jupiter's.

It has been assumed that hydrogen is negligibly incorporated into core-forming metals below $\sim$3 GPa, and therefore the presence of hydrogen in iron cores of small terrestrial bodies including the moon has not been considered. Here we performed high-pressure melting experiments on the Fe-H system under H$_2$-saturated conditions, combined with synchrotron X-ray diffraction (XRD) measurements. Results demonstrate substantial depression of the Fe-H melting curve compared to that for Fe at 1.0-3.3 GPa, indicating that hydrogen is incorporated into liquid iron even at low pressures less than 1 GPa and the solubility is enhanced with increasing pressure. Based on the density of liquid Fe-H derived from diffuse scattering signal in XRD data, we found that the solubility of hydrogen in liquid iron is about 0.9 wt% at 3.6 GPa and likely enhanced to 1.2 wt% at 5 GPa corresponding to lunar core conditions. The 1.2 wt% H causes 9 % density reduction, which might fully explain the observed density deficit of the lunar core with respect to iron, depending on the density estimate from seismological data.

We study graviton production from an oscillating inflaton condensate during reheating by systematically comparing Boltzmann and Bogoliubov descriptions for inflaton potentials of the form $V(\phi)\propto\phi^n$ around the minimum. The Bogoliubov framework provides a unified description of graviton production, capturing both perturbative and non-perturbative effects across short and long wavelengths, whereas the Boltzmann approach is restricted to perturbative production at short wavelengths. For the quadratic case ($n=2$), we find that the two approaches yield identical graviton spectra at short wavelengths, indicating that the Boltzmann treatments fully captures perturbative gravitational production in this regime. For steeper potentials ($n>2$), however, we identify a sizable contribution arising from the non-adiabatic transition between inflation and reheating. This component is naturally incorporated in the Bogoliubov formalism but absent in the Boltzmann description, and we show that it is important over a broad range of momenta. We derive analytic approximations within both frameworks that clarify the physical origin and scaling behavior of the spectrum. Our results delineate the regime of validity of Boltzmann approaches and show that, for steeper inflaton potentials, graviton production is governed by non-adiabatic transition dynamics for which the Bogoliubov formalism provides the most appropriate description.

The observed unusual behaviors of the orbits of Trans-Neptunian objects as well as the gravitational anomalies detected by the Optical Gravitational Lensing Experiment can be explained by assuming the existence of a ninth planet in the Solar System, having a mass of the order of $5-10M_{\oplus}$, and located at the distance of 300-1000 AU from the Sun. Since no optical counterpart of Planet 9 was observed, it is reasonable to assume that it has a very low luminosity. Various proposals on the nature of Planet 9 have been advanced, including the possibility that it is a black hole, an axion or a dark matter star. We propose that dark matter heating of Planet 9 could generate a thermal radio flux that could allow its observational detection, even if Planet 9 is a very dark object. We estimate the dark matter impact parameter, the mass and the kinetic energy deposition rates, as well as the surface temperature of Planet 9. By adopting a specific model for the time evolution of the planet, and assuming a long time capture of dark matter, the surface temperature of Planet 9, and the spectral features of the emitted radiation are obtained. Our results indicate that dark matter capture may provide an efficient mechanism for the heating of Planet 9, and also provide a specific observational signature of the planet. The numerical evaluations depend on the unknown value of the dark matter-ordinary matter interaction cross-section, with the estimates obtained as a function of its ratio and the saturation cross section for dark matter to deposit its entire energy. For a value of this ratio of $10^{-10}$, and for a dark matter density of the order of $1.32\times 10^{-17}$ g/cm$^3$, in a few Gyrs the surface temperature of Planet 9 can reach values of the order of 200 K, or even higher, with a maximum wavelength of around $\lambda_{max}=1.44\times 10^{-3}$ cm, situated in the infrared domain.

The radial acceleration relation reveals a nearly universal acceleration scale of order $10^{-10}\,\mathrm{ms^{-2}}$ in galactic dynamics, whose origin remains unexplained within conventional cold dark matter scenarios. We propose that this scale arises from an intrinsic infrared spectral property of the dark sector. Specifically, we hypothesize that a long-lived, color-neutral gluonic vacuum component survives the post-inflationary expansion era and, at large distances, develops a spectrally rigid lowest-weight structure. The microscopic seed for this infrared organization is provided by the QCD trace anomaly, which breaks classical scale invariance and, through dimensional transmutation, generates an intrinsic infrared scale in the gluonic sector. Requiring Lorentz covariance together with a positive-energy lowest-weight unitary realization then selects the Anti de Sitter algebra $\mathfrak{so}(2,3)$ as the simplest symmetry admitting a discrete tower of states with a representation-theoretically protected gap. The associated gap introduces a finite correlation length $r_{\texttt{c}}$ that controls the large-scale coherence of the dark sector. A self-gravitating condensate dominated by the lowest-weight mode generates a characteristic acceleration $g^{}_\star = G M_h / r_{\texttt{c}}^2$, naturally of the same order as the observed galactic acceleration scale, within standard Newtonian gravity. In this framework, the galactic acceleration scale appears as the gravitational imprint of a trace-anomaly-seeded infrared spectral gap in a coherent gluonic dark sector, rather than as a consequence of modified gravity or of galaxy-by-galaxy formation histories.

We investigate the cosmological implications of non-polynomial quasi-topological gravity (NPQTG), a novel class of modified gravitational theories in which the background dynamics is encoded in a single function of the Hubble parameter. This framework provides a minimal and theoretically consistent extension of general relativity, incorporating higher-curvature effects while preserving second-order field equations and avoiding higher-derivative instabilities. We first establish the general conditions for cosmological viability and construct explicit realizations, including polynomial, quartic, power-law and non-polynomial models, demonstrating how different functional forms lead to distinct expansion histories. Focusing on the quartic and power-law cases, we show that the resulting cosmological evolution reproduces the standard thermal history of the Universe and gives rise to an effective dark-energy sector of geometric origin, with dynamical equation-of-state behavior that can lie in the quintessence or phantom regime. We then confront the models with observational data from Type Ia Supernovae, Cosmic Chronometers, and Baryon Acoustic Oscillations, using a Bayesian MCMC analysis. We find that both models provide an excellent fit to the data, remaining fully compatible with current constraints and statistically competitive with $\Lambda$CDM. Our results demonstrate that NPQTG offers a simple and efficient framework for describing late-time cosmic acceleration with dynamical dark energy, while maintaining theoretical consistency and observational viability.

It is generally believed that high energy radiation (power-law components) can be mostly produced by a hot corona gas in the accreting black holes. There is a very popular hybrid disk radial coupling model that the inner part of cool Keplerian disk (or Shakura-Sunyaev disk) can produce advection-dominated accretion flow or corona-like structure, which can also generate outflows/jets. Here we argue that this simple coupling model cannot explain the whole hardness-intensity diagram of the spectral states and their transitions, and associated jets of a $X-$ray binary. Based on recent theoretical works on advective disk structures around a black hole, as well as many observational behaviors of a source, we conclude that there should be a third component of hot accretion flow with the radial coupling model, which can successfully explain all the spectral states and transitions. Interestingly, this model also provides a new scenario for the jet generation, launching, and evolution during the states with help of created barrier of the inner flow. We have also find out the jet kinetic power expression with our new jet generation scenario.

Bridging theory and observations is a key task to understand galaxy formation and evolution. With the advent of state-of-the-art observational facilities, an accurate modelling of galaxy observables through radiative transfer simulations coupled to hydrodynamic simulations of galaxy formation must be performed. We present a novel pipeline, dubbed RTGen, based on the Monte Carlo radiative transfer code RADMC-3D , and explore the impact of the physical assumptions and modelling of dust and gas phases on the resulting galaxy observables. In particular, we address the impact of the dust abundance, composition, and grain size, as well as model the atomic-to-molecular transition and study the resulting emission from molecular gas. We apply Monte Carlo radiative transfer a posteriori to determine the dust temperature in six different hydrodynamic simulations of isolated galaxies. Afterwards, we apply ray tracing to compute the spectral energy distribution, as well as continuum images and spectral line profiles. We find our pipeline to predict accurate spectral energy distribution distributions of the studied galaxies, as well as continuum and CO luminosity images, in good agreement with literature results from both observations and theoretical studies. In particular, we find the dust modelling to have an important impact on the convergence of the resulting predicted galaxy observables, and that an adequate modelling of dust grains composition and size is required. We conclude that our novel framework is ready to perform high-accuracy studies of the observables of the ISM, reaching few tens percent convergence under the studied baseline configuration. This will enable robust studies of galaxy formation, and in particular of the nature of massive clumps in high-redshift galaxies, through the generation of mock images mimicking observations from state-of-the-art facilities such as JWST and ALMA.

The analysis of biorelevant molecules in returned mission samples such as from the carbonaceous asteroid (162173) Ryugu is key to unravelling the role of extraterrestrial organics in the evolution of life. Coordinated analyses using chemically non-destructive techniques at the finest length-scales on pristine samples are particularly important. Here, we identify the chemical signature of uncommon globular and nitrogen-containing diffuse composite organic matter in asteroid Ryugu and map the distribution of biorelevant molecules therein with unprecedented detail. Using a novel electron-microscopy-based combination of vibrational and core-level spectroscopy, we disentangle the chemistry and nanoscale petrography of these organics. We show that some of these organics contain soluble and highly aliphatic components as well as NHx functional groups, that have formed in outer solar nebula environments before parent body incorporation or were synthesized by subtle fluid reactions on the final Ryugu asteroid. These novel coordinated analyses will open up new avenues of research on these types of precious and rare asteroidal dust samples.

The IceCube South Pole Neutrino Observatory has discovered the presence of a diffuse astrophysical neutrino flux at energies of TeV and beyond using neutrino induced muon tracks and cascade events from neutrino interactions. We present two analyses sensitive to neutrino events in the energy range \SI{1}{TeV} to \SI{10}{PeV}, using more than 10 years of IceCube data. Both analyses consistently reject a neutrino spectrum following a single power-law with significance $>4\,\sigma$ in favor of a broken power law. We describe the methods implemented in the two analyses, the spectral constraints obtained, and the validation of the robustness of the results. Additionally, we report the detection of a muon neutrino in the MESE sample with an energy of $11.4^{+2.46}_{-2.53} $\,\si{PeV}, the highest energy neutrino observed by IceCube to date. The results presented here show insights into the spectral shape of astrophysical neutrinos, which has important implications for inferring their production processes in a multi-messenger picture.

The Vera C. Rubin Observatory LSST is expected to discover tens of millions of new Active Galactic Nuclei (AGNs). The survey's exceptional cadence and sensitivity will enable UV/optical/NIR monitoring of a significant fraction of these objects. The unprecedented number of sources makes spectroscopic follow-up for the vast majority of them unfeasible in the near future, so most studies will have to rely on photometric redshifts estimates which are traditionally much less reliable for AGN than for inactive galaxies. This work presents a novel methodology to constrain the photometric redshift of AGNs that leverages the effects of cosmological time dilation, and of the luminosity and wavelength dependence of AGN variability. Specifically, we assume that the variability can be modeled as a damped random walk (DRW) process, and adopt a parametric model to characterize the DRW timescale ($\tau$) and asymptotic amplitude of the variability (SF$_\infty$) based on the redshift, the rest-frame wavelength, and the AGN luminosity. We construct variability-based photo-$z$ priors by modeling the observed variability using the expected DRW parameters at a given redshift. These variability-based photometric redshift (VAR-PZ) priors are then combined with traditional SED fitting to improve the redshift estimates from SED fitting. Validation is performed using observational data from the SDSS, demonstrating significant reduction in catastrophic outliers by more than 10% in comparison with SED fitting techniques and improvements in redshift precision. The simulated light curves with both SDSS and LSST-like cadences and baselines confirm that, VAR-PZ will be able to constrain the photometric redshifts of SDSS-like AGNs by bringing the outlier fractions down to below 7% from 32% (SED-alone) at the end of the survey.

The polarized light of the cosmic microwave background is sensitive to new physics that violates parity symmetry. For example, the interaction of photons with the fields of elusive dark matter and dark energy could cause a uniform rotation of the plane of linear polarization across the sky, an effect known as cosmic birefringence. We extract the cosmological rotation angle, $\beta$, using Bayesian analysis of parity-violating correlations, $EB$ and $TB$, of polarization data from the Atacama Cosmology Telescope (ACT) Data Release 6. We use prior probabilities for instrumental miscalibration angles derived from the optics model for the ACT telescope and instruments, and marginalize over a residual intensity-to-polarization leakage. We measure $\beta = 0.215^\circ\pm 0.074^\circ$ (68\% confidence level), which excludes $\beta=0$ with a statistical significance of $2.9\sigma$. Although there remain systematics in the ACT data that are not understood and do not allow us to draw strong cosmological conclusions, this result is consistent with previous independent results from the WMAP and Planck missions. It is suggestive that independent data sets and analyses using different methodologies have yielded the same sign and comparable magnitudes for $\beta$.

We present $100$ full-sky quasar spectrophotometric mock catalogs with smooth redshift evolution from $z=0$ to $z\sim 4$, tailored to analyze the Gaia-unWISE Quasar Catalog (Quaia). In particular, we apply a novel hierarchical nonlocal nonlinear bias scheme (Hicobian) to dark matter fields generated through Augmented Lagrangian Perturbation Theory on the lightcone (WebON code), calibrating the free parameters of the bias model on Abacus quasar HOD mock catalogs tuned to reproduce DESI Early Data Release observations in real and redshift space. After having obtained such accurate spectroscopic catalogs, we inject in the mocks the observational effects characterizing the Quaia catalog: (i) spectrophotometric redshift uncertainties, (ii) the angular selection function, and (iii) the redshift number counts distribution. We assess the accuracy of our catalogs by validating a number of summary statistics: the full-sky QSO maps, the redshift uncertainty distributions as a function of redshift, the redshift $n(z)$ distribution, the angular power spectra and their normalized covariance matrices, and the angular two-point correlation functions. We find excellent agreement between these metrics from the mocks and from the Quaia catalog. We publicly release the mock catalogs to the community.

Using a two-fluid approach, we consider the properties of relativistically nonlinear (arbitrary $a_0$), circularly polarized \EM\ waves propagating along magnetic field in electron-ion and pair plasmas. Dispersion relations depend on how wave intensity scales with frequency, $a_0 (\omega)$. For superluminal branches, the nonlinear effects reduce the cut-off frequency, while the general form of the dispersion relations $\omega(k)$ remains similar to the linear case. For subluminal waves, whistlers and Alfven, a new effect appears: dispersion curves effectively terminate at finite $\omega^\ast - k^\ast$, where the group velocity becomes zero. Qualitatively, subluminal modes with fluctuating electric field larger than the guide field, $E_w (\omega) \geq B_0$, cannot propagate. In extended systems, e.g., within magnetospheres of neutron stars, this leads to opening of the magnetosphere by a strong wave.

Cosmological N-body simulations have long suggested that the pseudo phase-space density (PPSD), $\rho/\sigma^3$, of cold dark matter haloes follows the universal relation $\rho/\sigma^3 \propto r^{\chi}$, with $\chi \approx -1.875$, as predicted by spherical secondary-infall similarity solutions. This power law appears to hold despite the fact that neither the density $\rho(r)$ nor velocity dispersion $\sigma(r)$ follow universal power law relations individually, even at fixed mass. We analyze 246 host haloes from the \textit{Symphony} suite of high-resolution cosmological zoom-in simulations, to consistently measure PPSD profiles across host masses from $10^{11}$ to $10^{15} M_\odot$. We find that the PPSD systematically deviates from a power law, and that haloes with larger deviations from Jeans equilibrium systematically develop steeper average PPSD slopes. This result suggests that the PPSD is not universal; instead, it is linked to a halo's degree of dynamical equilibrium, which is ultimately set by halo formation history. As a result, we show that secondary halo properties such as concentration and accretion rate inherit significant correlations with the PPSD slope. Moreover, our hosts' PPSD profiles are remarkably consistent with predictions from one-dimensional self-similar fluid collapse models, indicating that three-dimensional structure, velocity anisotropy, and filamentary accretion all play negligible roles in shaping the PPSD. Thus, we argue that the PPSD is shaped by mass assembly alone, and that its non-universality reflects the diversity of halo growth histories.

We carry out an independent re-analysis of the Dark Energy Spectroscopic Instrument (DESI) public dataset, focusing on extensions to the standard cosmological model, $\Lambda$CDM. Utilizing the dataset and Effective Field Theory (EFT)-based pipeline described in Paper 1, we constrain cosmological models with massive neutrinos ($\Lambda$CDM+$M_\nu$), spatial curvature ($o\Lambda$CDM), dynamical dark energy ($w_0w_a$CDM), and their combinations using the power spectrum and bispectrum of DESI galaxies and quasars. Our work also presents the first measurements of relevant non-minimal cosmological parameters from the combination of cosmic microwave background (CMB) and DESI full-shape (FS) data, which are made possible thanks to carefully chosen priors on EFT parameters. We find that the addition the FS likelihood to DESI's baryon acoustic oscillation (BAO) data improves the limits on the spatial curvature by a factor of two over the BAO only results, though the improvements are less significant with the CMB data. The dark energy equation of state figure-of-merit increases both with and without the supernovae data (SNe), by $\approx30\%$ and $\approx20\%$ relative to the CMB+BAO and CMB+BAO+SNe results, respectively. Our FS likelihood also yields the strongest CMB-independent constraint on the total neutrino mass $M_\nu<0.32\,{\rm eV}$, with the $30\%$ improvement due to the bispectrum. In combination with the CMB, we find a $14\%$ improvement assuming the $\Lambda$CDM+$M_\nu$ model (yielding $M_\nu<0.059\,{\rm eV}$), but this increases to $22\%$ when using non-minimal backgrounds: $M_\nu<0.097\,{\rm eV}$ in $o\Lambda$CDM+$M_\nu$ and $M_\nu<0.13\,{\rm eV}$ in $w_0w_a$CDM+$M_\nu$. Overall, our work illustrates that robust and substantial gains in constraining power can be obtained by incorporating the FS power spectrum and bispectrum measurements in analyses of non-minimal cosmological models.

We present the quasar catalog from Data Releases 10 to 12 of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Quasar Survey, comprising quasars observed between September 2021 and June 2024. We robustly identified $11,346$ quasars, of which $5,386$ are newly discovered objects not present in the Million Quasars catalog. This release brings the total number of quasars identified by the 12-year LAMOST survey to $67,521$, of which $29,513$ are newly discovered. While the absolute flux calibration for LAMOST quasar spectra from Data Releases 6 to 9 was previously performed using the SDSS/PanSTARRS1 multi-band photometric data, the inherent variability of quasars can affect the flux accuracy. To address this limitation, we recalibrated the LAMOST spectra using (quasi-)simultaneous photometric data from Zwicky Transient Facility (ZTF), which has conducted high-cadence sky monitoring since March 2018. Based on the recalibrated single-epoch spectra, we estimated the emission line fluxes, continuum fluxes, and virial black hole masses. These improved spectra facilitate direct comparison with the spectra of common quasars from the Sloan Digital Sky Survey (SDSS), enabling searches for rare quasars, such as changing-look quasars exhibiting the appearance or disappearance of broad emission lines and broad absorption line quasars. The combined dataset of photometry and multi-epoch spectra will enhance the detections of AGN-related transients, such as Bowen fluorescence flares and extreme variability quasars, thereby improving our understanding of quasar variability.

The effective field theory (EFT) of dark energy provides a model-independent framework for studying cosmology within scalar-tensor theories. In this work, we explore how the time evolution of the cosmological background, inferred from cosmic chronometer measurements of the Hubble parameter, can be used to reconstruct the relevant EFT functions. Our approach enables the direct determination of these EFT functions from observational data without assuming any specific cosmological model. This makes it possible to test the background evolution of a wide range of dark energy models, including the $\Lambda$CDM model. We further demonstrate how the reconstructed EFT functions can be applied to constrain concrete theories, such as the quintessence model.

The observational link between long gamma-ray bursts (GRBs) and broad-lined stripped-envelope core-collapse supernovae (SNe Ic-BL) is well established. Significant progress has been made in constraining what fraction of SNe Ic-BL may power high- or low-luminosity GRBs when viewed at small off-axis angles. However, the GRB-SN connection still lacks a complete understanding in the broader context of massive-star evolution and explosion physics. Models predict a continuum of outcomes for the fastest ejecta, from choked to ultra-relativistic jets, and observations from radio to X-rays are key to probing these scenarios across a range of viewing angles and velocities. Here, we present results from a coordinated radio-to-X-ray campaign targeting nearby (z<=0.1) SNe Ic-BL designed to explore this diversity. With eight new radio-monitored events and updated data for one previously observed SN, we further tighten constraints on the fraction of SNe Ic-BL as relativistic as SN 1998bw/GRB 980425. We identify SN 2024rjw as a new radio-loud event likely powered by strong interaction with circumstellar material (CSM), and add evidence supporting a similar interpretation for SN 2020jqm. We also establish new limits on the properties of radio-emitting ejecta with velocities consistent with cocoons from choked jets, highlighting SN 2022xxf as a promising cocoon-dominated candidate. These results refine our understanding of the continuum linking ordinary SNe Ic-BL, engine-driven explosions, and GRBs, and contribute to building a sample that will inform future multi-messenger searches for electromagnetic counterparts to high-energy neutrinos.

Odd Radio Circles (ORCs) are a new class of extragalactic object, with large rings of faint radio continuum emission typically spanning 100s of kpc; their origins are unknown. Previous optical spectroscopy of the central galaxy in ORC4, a classic isolated ORC, revealed spatially-extended ionized gas with strong [OII] emission and line ratios consistent with LINER emission. We present new Keck/KCWI+KCRM integral field spectroscopy covering multiple strong optical emission lines to measure the extent, morphology, and spatially-resolved kinematics and line ratios of the ionized and neutral gas in the ORC4 central galaxy. We find that [OII] is the strongest optical emission line in this massive, old galaxy, and the [OII] emission is detected to larger radial extent than the other optical lines. The gas kinematics show strong spatial asymmetries, high velocity gradients (>100 km/s), and high velocity dispersion (~200 km/s). The emission line ratios are most consistent with shock models with shock velocities of ~200-300 km/s and are not fit well by stellar or AGN photoionization models. These findings are consistent with a model in which the gas in the ORC4 central galaxy is the result of shock ionization in and around the central galaxy, likely due to mixing and cooling of gas associated with the event that created the large-scale radio ring of emission that identified this source as an ORC.

The rotation-activity relationship of G-type stars results from surface magnetic fields emerging from the interior. How the magnetic flux and its emergence rate scale with rotation rate are not well understood, both observationally and theoretically. We aim at constraining the emerging magnetic flux as a function of the rotation rate in solar-type stars by numerical simulations compared to empirical constraints set by direct measurements of stellar magnetic fields. We use our Flux Emergence And Transport (FEAT) model for stars with a range of power-law slopes for the dependence of emerging flux on rotation. Complementing this with a heuristic account of the main flux components, we model the resulting mean unsigned field strength as a function of the rotation rate. We compare the results with the Zeeman-intensification measurements and spectropolarimetric data of solar-type stars. Deviations of the model from observations of G stars correlate strongly with stellar metallicity ($r=0.83$) and effective temperature ($r=-0.76$), with a combined coefficient of 0.90, reflecting the dependence of magnetic activity on these two parameters. Correcting for these effects with multilinear regression, we find that magnetic flux emergence rates must scale steeply with rotation power-law exponent of about 1.9) to reproduce observed field strengths, significantly exceeding the estimates in the literature. We also provide correction factors for metallicity and temperature for measurements of early-G-type stellar magnetic fields. Stellar magnetic flux emergence rates scale steeply with rotation, requiring active-region fields to dominate the total surface flux on rapid rotators, whereas small-scale-dynamo fields dominate for slow rotators like the Sun. Metallicity significantly influences the rotation-magnetism relationship, necessitating sample-dependent corrections for accurate stellar dynamo modelling.

There is growing evidence from gravitational-wave observations that some merging black holes are created from previous mergers. Using the prediction that these hierarchically merged black holes have dimensionless spin magnitudes of $\chi \approx 0.69$, we identify a subpopulation in the gravitational-wave data consistent with a hierarchical-merger origin in dense star clusters. This subpopulation's primary mass distribution peaks at $17.0^{+18.3}_{-4.4},\mathrm{M}_{\odot}$, which is approximately twice as large as its secondary mass distribution's mode ($10.5^{+29.7}_{-4.7},\mathrm{M}_{\odot}$), and its spin tilt distribution is consistent with isotropy. Our inferred secondary mass distributions imply that isolated binary evolution may still be needed to explain the entirety of the $9\,\mathrm{M}_{\odot}$ peak. Surprisingly, we find that the rate of hierarchical mergers may evolve more steeply with redshift than the rest of the population ($98.0\%$ credibility): the fraction of all binary black holes that are hierarchically formed at $z=0.1$ is $0.03^{+0.05}_{-0.02}$, compared to $0.09^{+0.11}_{-0.07}$ at $z=1$. This provides an explanation for the previously discovered broadening of the effective spin distribution with redshift. Our results have implications for star cluster formation histories, as they suggest the potential existence of a high-redshift population of massive, compact clusters.

Spectroscopic observations by the James Webb Space Telescope (JWST) have revealed young, compact, high-redshift ($z$) galaxies with high nitrogen-to-oxygen (N/O) ratios. GN-z11 at z=10.6 is one of these galaxies. One possible scenario for such a high N/O ratio is pollution from supermassive stars (SMSs), from which stellar winds are expected to be nitrogen-rich. The abundance pattern is determined by both galaxy evolution and SMS pollution, but so far, simple one-zone models have been used. Using a galaxy formation simulation, we tested the SMS scenario. We used a cosmological zoom-in simulation that includes chemical evolution driven by rotating massive stars (Wolf-Rayet stars), supernovae, and asymptotic giant branch stars. As a post-process, we assumed the formation of an SMS with a mass between $10^3$ and $10^5$ $M_\odot$ and investigated the contribution of its ejecta to the abundance pattern. The N/O ratio was enhanced by the SMS ejecta, and the abundance pattern of GN-z11, including carbon-to-oxygen and oxygen-to-hydrogen ratios, was reproduced by our SMS pollution model if the pollution mass fraction ranges within 10-30 per cent. Such a pollution fraction can be realized when the gas ionized by the SMS is polluted, and the gas density is $10^4$-$10^5$ cm$^{-3}$ assuming a Strömgren sphere. We also compared the abundance pattern with those of other N/O-enhanced high-$z$ galaxies. Some of these galaxies can also be explained by SMS pollution.

Type IIn supernovae (SNe IIn) are hydrogen-rich explosions embedded in dense circumstellar medium (CSM), which gives rise to their characteristic narrow hydrogen emission lines. The nature of their progenitors and pre-explosion mass loss remains, however, poorly understood. Using high-resolution Hubble Space Telescope (HST) imaging, we analyze the local stellar environments of a volume-limited sample (z < 0.02) of 31 SNe IIn. The environments of SNe IIn are found to be very diverse; the SN could reside within a star-forming region (Class 1), outside a star-forming region (Class 2), or in much older environments without any obvious signs of star formation (Class 3). The bright SNe IIn (Mpeak < -19.5 mag) predominantly occur in Class 1 environments, indicative of very massive progenitors, while the faint SNe IIn (Mpeak < -15.5 mag) are associated with Classes 2 and 3 environments, suggesting the least massive progenitors. Meanwhile, normal SNe IIn with -19.5 < Mpeak < -15.5 mag occur in all three types of environments, suggesting a diversity in their progenitor mass, lifetime, and evolutionary pathways. Moreover, the directly detected SN IIn progenitors are systematically brighter and/or bluer than the youngest stellar populations in their environments, suggesting that they were either in a non-quiescent state when observed or had experienced binary interactions. These results point to a significantly diverse origin for progenitors of SNe IIn, spanning a wide range of masses, evolutionary stages, and potential binary interaction histories.

A review of the current status of the field of Ultra-High-Energy Cosmic Ray (UHECR) including a summary of remaining open questions was presented in the white paper "Ultra-High Energy Cosmic Rays: at the Intersection of the Cosmic and Energy Frontiers" (Astropart. Phys. 147 (2023) 102794; arXiv:2205.05845). The authors concluded that two types of next-generation detectors are needed to answer these questions: high-accuracy instruments and detectors that maximize exposure at the highest energies. The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA), a proposed dual-satellite observatory, exemplifies the latter class and is designed to increase statistics of the highest-energy cosmic rays and to detect very-high-energy neutrinos following multi-messenger alerts. POEMMA-Balloon with Radio (PBR) implements a compact, balloon-borne version of the POEMMA concept, adapted for a Super-Pressure Balloon flight from Wanaka, New Zealand, with an expected campaign exceeding 20 days. PBR couples a wide field-of-view Schmidt telescope and a hybrid optical focal surface with a dedicated radio instrument to deliver simultaneous, complementary measurements of extensive air showers. The mission will validate the fluorescence detection strategy from space and raise technology readiness for a POEMMA-like space mission by observing UHECR-induced fluorescence light from suborbital altitudes, obtaining the first simultaneous optical Cherenkov and radio observations of high-altitude horizontal air showers above the cosmic-ray knee (E>3PeV), enabling energy-spectrum and composition studies at the PeV scale, and performing follow-ups of multi-messenger alerts to search for very-high-energy neutrinos via upward-going air showers. This paper summarizes the PBR payload and its expected performance.

The CASPER (Collaboration for Astronomy Signal Processing and Electronic Research) toolflow is a widely used framework for designing and implementing digital signal processing systems, particularly in the field of radio astronomy. It provides a set of tools and libraries that enable researchers to create custom hardware and software solutions for processing astronomical data. The CASPER toolflow has been instrumental in the development of Field-Programmable Gate Array (FPGA) based digital instruments for various radio telescopes, enabling for real-time data processing and analysis. However, the current frontend tool that CASPER uses for high-level FPGA design is based on Model Composer integrated into MATLAB/Simulink, which is a proprietary software. In this paper, we introduce Scilab as a new frontend tool for the CASPER toolflow. Scilab is an open-source software platform for numerical computation and data visualization, which offers a similar environment to MATLAB/Simulink for designing CASPER blocks, generating FPGA Intellectual Property (IP) cores, and simulating Digital Signal Processing (DSP) systems. We present our implementation of Scilab in the CASPER toolflow and demonstrate its capabilities by developing an FPGA based spectrometer on a RFSoC4x2, a commonly used CASPER platform well suited to radio astronomy applications. We have also developed Scilab support for other CASPER compatible platforms. Our results show that Scilab can successfully be used as an alternate frontend for CASPER-based designs.

It is well known that, in potentials that are quadratic near the minimum but shallower away, such as small $\alpha$ ($\ll M_P^2$) attractors, the inflaton condensate fragments into localized compact objects known as oscillons during self-resonance preheating. In this work we investigate the self-resonance in deformed $\alpha$-attractor T-model with a Gaussian feature near the minimum, distant from inflation's end. Linear analysis reveals altered resonance bands and deformed Floquet charts dependent on feature parameters. In fully nonlinear lattice simulations, we find that the gradient energy transfer is largely independent of the potential feature parameter $h$. In contrast, after resonance terminates, the subsequent evolution of gradient energy becomes strongly dependent on $h$. Statistical analysis reveals that models with the potential feature produce larger number of smaller oscillons, with a reduced energy stored in these objects, increasingly suppressed as the magnitude of $h$ grows. By tracking the total energy and the gradient energy contained in oscillons, we find that in models with nonzero $h$ oscillons are systematically shorter-lived, with this effect strengthening for larger $h$. The gravitational wave emission is dominated by the resonance stage and is strongly suppressed once oscillons form. Potential features leave the low-frequency spectrum largely unchanged but significantly modify the high-frequency tail. Although a complete reheating description requires external couplings and higher-resolution simulations, clear qualitative differences of cosmic expansion history already emerge within our simulated time window. These results highlight the important role of potential features in shaping reheating dynamics and their cosmological implications, and provide a deeper understanding of preheating dynamics and the properties of oscillons.

Pair-instability supernovae leave behind no compact remnants, resulting in a predicted gap in the distribution of stellar black-hole masses. Gravitational waves from binary black-hole mergers probe the relevant mass range and analyses of the LIGO-Virgo-KAGRA catalog (GWTC-4) indicate a possible mass cutoff at $40$-$50\,M_\odot$. However, the robustness of this result remains unclear. To this end, using full Bayesian parameter estimation, we simulate gravitational-wave catalogs with and without such a mass cutoff, then test whether its presence or absence is correctly inferred with parametric population models. For catalogs similar to GWTC-4, confident identification of a cutoff is not guaranteed, but the best constraints among our simulations are compatible with results from GWTC-4 when the model includes a cutoff. Conversely, spurious identification of a cutoff is unlikely. For catalogs expected by the end of the O4 observing run, uncertainty in the cutoff mass is reduced by $\gtrsim20\,\%$, but a cutoff at 40-50$M_\odot$ yields only a lower bound on the $^{12}\mathrm{C}(\alpha,\gamma)^{16}\mathrm{O}$ reaction rate, our most stringent constraints on the S-factor at $300\,\mathrm{keV}$ being $S_{300}\gtrsim125\,\mathrm{keV}\,\mathrm{b}$ at $90\,\%$ credibility. Relative uncertainties on the Hubble parameter $H_0$ from gravitational-wave data alone can still be up to $100\,\%$. We also analyze GWTC-4 with the nonparametric PixelPop population model, finding that some mass features are more prominent than in parametric models but a sharp cutoff is not required. However, the parametric model passes a likelihood-based predictive test in GWTC-4 and the PixelPop results are consistent with those from our simulated catalogs with a cutoff. Such tests are necessary to make astrophysical claims from gravitational-wave catalogs.

We introduce $\texttt{GPUmonty}$, a CUDA/C-based Monte Carlo radiative transfer code accelerated using graphics processing units (GPUs). $\texttt{GPUmonty}$ derives from the CPU-based code $\texttt{grmonty}$ and offloads the most computationally expensive stages of the calculation -- superphoton generation, sampling, tracking, and scattering -- to the GPU. Whereas $\texttt{grmonty}$ handles photons sequentially, $\texttt{GPUmonty}$ processes large numbers of superphotons concurrently, leveraging the single-instruction, multiple-thread (SIMT) execution model of modern GPUs. Benchmarks demonstrate a speedup of about $12\times$ relative to the original CPU implementation on a single GPU, with runtime limited primarily by register pressure rather than compute or memory bandwidth saturation. We validate the implementation through analytic tests for a optically thin synchrotron sphere, as well as comparisons with $\texttt{igrmonty}$ for scattering synchrotron sphere and GRMHD simulation data. Relative errors remain below a percent level and convergence is consistent with the expected $N_{\rm s}^{-1/2}$ Monte Carlo scaling. By significantly reducing computational costs, GPUmonty enables the extensive parameter space surveys and faster spectra modeling required to interpret horizon-scale observations of supermassive black holes. $\texttt{GPUmonty}$ is publicly available under the GNU General Public License.

When low-mass ($\lesssim 2$ $M_{\odot}$) red giants lose their envelopes to a companion just before the helium flash, the resulting mass transfer can produce binaries hosting hot subdwarfs, horizontal branch stars, and undermassive red clump stars. Recent work predicts a continuum of such products, from fully stripped hot subdwarfs to partially stripped horizontal branch and red clump stars, and suggests that young, metal-rich RR Lyrae can form when partial stripping leaves a helium-burning star in the instability strip. To enable direct comparison with observations, we model these binaries in a simulated Milky Way-like galaxy with a realistic metallicity-dependent star formation history and 3D dust map, generate epoch astrometry using Gaia's scanning law, and fit it with the cascade of astrometric models applied in Gaia DR3. We compare the simulated population to DR3 observations of hot subdwarfs, RR Lyrae, and red giants with high astrometric mass functions. The model significantly overpredicts the number of hot subdwarfs with astrometric binary solutions, partly because the predicted flux ratios are more unequal than observed. It also predicts $\gtrsim 100$ RR Lyrae with DR3 astrometric orbital solutions, while none are observed. We conclude that RR Lyrae in au-scale binaries may be substantially rarer than predicted. In contrast, the model plausibly explains the population of red clump stars with high astrometric mass functions, which we interpret as potential black hole impostors. We predict that $\sim 10 \times$ more stripped-star binaries will be detectable in DR4, whose sensitivity to longer periods will more strongly test wide-orbit systems.

We present an example cluster follow-up study with Euclid. Our target, a $z\sim 1.74$ candidate cluster nicknamed the 'Puddle', was initially discovered by the Massive and Distant Clusters of WISE Survey 2 as a $z_\mathrm{phot}\sim 1.65$ candidate cluster. It was also detected independently as a $z_\mathrm{phot}\sim 1.5$ candidate with the two cluster-finding algorithms in Euclid Quick Release 1 (Q1). A Keck MOSFIRE spectrum shows the brightest nucleus is at $z=1.74$ and is dominated by an active galactic nucleus. Our analysis focused on the galaxy population and the brightest cluster galaxy (BCG), and is based on Euclid and ancillary photometry. Compared to similar fields, we measured an overdensity of $110\pm 14$ galaxies with $H_\mathrm{E}\leq 22.25$ in a 2' radius around the BCG. About $18\pm 4$% of the completeness-corrected galaxy population is red, which is consistent with some clusters at $z>1.5$ but lower than others. Euclid imaging revealed that six or seven galaxies appear to be assembling to form the future BCG. Spectral energy distribution fitting suggests that the merging BCG has a stellar mass of $5.7\pm 0.3\times 10^{11}\,M_\odot$ and that it experienced a short burst of star formation $\sim 300\,$Myr ago. Its morphology and star-formation history suggest that the proto-BCG is a more evolved version of the merging core of SPT2349$-$56. These systems indicate that multiobject mergers might be a common BCG formation process. Assuming a similar density of mergers in the Euclid Wide Survey, we expect that Euclid will discover approximately 400 assembling BCGs by the end of its mission.

We point out that several terms in the third-order corrections to the slow-roll power spectra presented by Ballardini et al. [1] are incorrect. The authors of that work claim that their result differ from the ones originally presented by Auclair & Ringeval [2] due to some different approximation schemes. However, in our original work, all terms at all orders have been derived exactly and any difference between two expansions performed at the same pivot wavenumber signals a problem. As we show in this comment, Ballardini et al. [1] have misevaluated some definite three-dimensional integrals by integrating a truncated Taylor expansion instead of Taylor expanding an integral. Our claim is backed-up with a Monte-Carlo numerical integration of the incriminated three-dimensional integrals, which, unsurprisingly, matches the analytical value derived in Auclair & Ringeval [2].

As the oldest building blocks of our Galaxy, globular clusters retain the archaeological footprint of the early stellar environments. Accurate absolute ages of globular clusters are required to interpret this ancient record. Existing dating techniques often produce precise but discordant ages, suggestive of systematic errors in excess of 1 Gyr. The James Webb Space Telescope (JWST) has unlocked a new dating method that leverages the cooling behavior of previously unobservable brown dwarf members. With a largely independent set of systematic errors, this new method provides a new consistency test for more established methodologies. I present a likelihood-based histogram-free method to derive globular cluster ages from multi-band JWST photometry of cluster members near and below the hydrogen-burning limit. By applying the method to a large set of simulated observations, I establish that formal age errors (i.e. errors based on measurement uncertainties alone) under 0.2 Gyr are attainable for nearby globular clusters. I also evaluate the significance of associated systematic effects, including the chemical heterogeneity of globular clusters (multiple populations), unresolved binary systems and uncertainties in brown dwarf cooling rates. As with other methods of age determination, systematic effects dominate the error budget (in selected cases, by over an order of magnitude), but may be reduced with more sophisticated analysis. Finally, I provide a lookup table for determining the number of observations, exposure times and temporal baselines required to estimate the age of a given cluster to a prescribed precision.

Wide (soft) binaries are expected to be rapidly disrupted in dense stellar environments, yet they are observed in both the Galactic field and open clusters (OCs). In this paper, we investigate the formation and disruption of wide binaries in star clusters using direct N-body simulations. We perform simulations containing 10,000 objects with varying binary fractions and initial bulk rotation to give an in-depth look into the dynamical evolution of wide binaries in star clusters. We find that wide binaries dominate early disruption and formation processes during the initial high-density phase of cluster evolution. We propose two semi-analytical models to reproduce the evolution of the wide-binary population in simulations. The exponential model consists of an early, rapid-disruption phase with a time less than 10 Myr, driven by frequent encounters at high density, and a longer, relaxation-driven phase between 200 and 300 Myr. The broken power-law model provides break timescales when the decrease of wide binaries slows down during the early and long-term disruption. All timescales from both models agree with each other and decrease with increasing stellar density induced by high primordial binary fraction and cluster rotation. Wide binary disruption is mostly responsible for the early decline in the total binary fraction of the cluster. Such disruption leads to the decrease of radial binary fraction toward the cluster center until 500 Myr. Our results suggest low-density OCs or stellar groups younger than 10 Myr as the optimal environments for detecting wide binaries and provide a physical framework for understanding their contribution to the Galactic field population.

The James Webb Space Telescope provides an exciting opportunity to investigate stellar systems located in heavily obscured regions like the Galactic bulge. Possibly, the most enigmatic among them is Terzan 5: long classified as a globular cluster, it is now known to host distinct stellar populations with different iron abundances (ranging approximately from [Fe/H]=-$0.8$ to [Fe/H]=$+0.3$ dex). Indeed the chemical and structural properties collected so far suggest that it is the remnant of one of the primordial clumps that contributed to the early assembly of the bulge, a so-called "Bulge Fossil Fragment". Here we present a new photometric analysis of Terzan 5 based on JWST/NIRCam observations in the F115W and F200W filters, as well as archival HST/ACS optical (F606W and F814W) data. The dataset overcomes the severe and spatially variable extinction along the line of sight and yields the deepest color-magnitude diagram ever obtained for Terzan 5. Proper motion selections and high-resolution differential reddening corrections allow us to isolate bona fide cluster members and to provide an unprecedented view of the main-sequence turn-off region. We clearly identify two main components and determine their respective ages: the old, sub-solar component has an age of 12.5 $\pm$ 0.5 Gyr, while the super-solar component is significantly younger with an age of 4.7 $\pm$ 0.5 Gyr. Interestingly, we also find hints of an even younger main sequence turn-off and sub-giant branch, consistent with the presence of a further stellar component with an age of 3.8 $\pm$ 0.5 Gyr. There is also evidence of a blue plume populated by stars as bright as $m_{\rm F115W}\sim 17.4$, suggesting a prolonged period of star formation extending up to 2.5 Gyr ago.

Spectroscopic redshift determination from R ~ 450 slitless grism data is vulnerable to line-identification degeneracy when only a small number of features are detected. This paper presents a custom extraction pipeline for the 473 published Euclid Q1 lens candidates that fall within the NISP slitless-spectroscopy footprint, combining extended wavelength coverage (9200-18500 Angstroms), VIS-guided spectral deblending, and an optional LangGraph/Claude morphology and quality-control module. All quantitative results reported here use the deterministic heuristic fallback rather than the AI-assisted branch. The extended range yields new features below 12500 Angstroms in 312/473 source spectra (66%) and 324/473 deflector spectra (69%). A VIS-J color prior calibrated on 19711 DESI galaxies is used to anchor deflector redshifts. Against 29 DESI spectroscopic redshifts, blind deflector assignments agree in 2/29 systems at |dz/(1+z)| < 0.03; with the color prior, agreement is 3/29 at < 0.03 and 6/29 at < 0.05. The anchored pipeline assigns deflector redshifts to 403 systems; of these, 322 also have a source redshift satisfying z_src > z_def, forming complete pairs with valid lensing geometry. In an isolated-galaxy control sample, two internal redshifting methods agree for 20% of spectra at |dz/(1+z)| < 0.01 and 44% at < 0.05, indicating limited internal self-consistency even without lens blending. These results show that extended-coverage Euclid NISP spectra recover many additional candidate features and can support anchored lens redshifts, but they do not support reliable blind deflector redshifts from the grism data alone.

The first Euclid Quick Data Release (Q1) provides extensive imaging and spectroscopic data for hundreds of millions of photometric objects across several deep fields. Accurate classifications and photometric redshifts (photo-z) for these sources are crucial to maximizing the value of these data. In this work, we perform source classification and photo-z estimation for the Euclid Deep Field North (EDF-N) around the North Ecliptic Pole, using a deep learning framework (DeepDISC) that learns and infers using 9-band images simultaneously. We train three dedicated models for (1) source detection and classification, (2) galaxy photo-z, and (3) quasar photo-z. The Euclid Q1 input source catalog, and classifications and spectroscopic redshifts (spec-z) from the Dark Energy Spectroscopic Instrument Data Release 1 are adopted as our training data. DeepDISC source detection achieves overall completeness of ~93% and purity of ~80% if using the Euclid source catalog as the ground truth. Using a JWST source catalog within EDF-N as the reference, we estimate a true purity of ~ 90% for DeepDISC sources. About 99.2%, 99.0%, and 84.8% of stars, galaxies, and quasars, respectively, are correctly recovered with their spectroscopic classifications. The DeepDISC photo-zs show good agreement with spectroscopic redshifts, for both galaxies and quasars. Comparisons with other Euclid Q1 products demonstrate that DeepDISC provides comparable or improved performance in source detection/deblending, classification and photo-z, especially for quasars. These results demonstrate the potential of pixel-level deep learning approaches for large-scale sky surveys such as Euclid and Roman, which will continue to improve with better training labels. We release the full DeepDISC source catalog (~13 million objects) for EDF-N with classifications and photo-zs, including photo-z probability distributions.

Quasars exhibit stochastic variability across wavelengths, typically well described by a Damped Random Walk (DRW). Occasionally, however, they undergo extreme luminosity changes--known as flares--that represent significant departures from this baseline behavior and provide valuable probes of accretion disc dynamics and the physics of supermassive black hole fueling. Although modern transient surveys have spurred growing interest in flare detection, no systematic search has yet been conducted within the legacy SDSS Stripe 82 dataset, which contains 9,258 spectroscopically confirmed quasars observed over a ~10-year baseline. The principal statistical challenge is distinguishing these rare events from the ever-present stochastic variability. To address this, we present FLARE (Flare detection via physics-informed Learning, Anomaly scoring, and Recognition Engine), a modular three-stage framework for detecting transient flares in quasar light curves. FLARE models baseline DRW behavior, applies statistical anomaly scoring to flag candidate events, and employs a recognition engine to verify detections. For the Stripe 82 implementation, we deploy two complementary baselines--a physics-informed probabilistic Gated Recurrent Unit (GRU) trained on simulated DRW light curves, and an iterative Ornstein-Uhlenbeck (OU) process fitted directly to observed data with outlier masking--followed by Extreme Value Theory (EVT) for anomaly scoring. We benchmark twelve open-weight and proprietary Vision Language Models (VLMs) as recognition engines for final candidate verification. Detection is performed on r-band light curves, with candidates cross-checked against g-band data to rule out instrumental artifacts. Applying this framework, we identify 51 quasars exhibiting distinct flaring activity.

Hot Neptunes in the sub-Jovian desert offer unique insights into planetary system evolution, retaining signatures of dynamical processes that shaped their present-day architectures. Many of these planets exhibit polar orbits, yet the mechanisms responsible for these misalignments between the stellar spin axis and planet orbit normal remain under debate. GJ 436 b stands among the very few hot Neptunes with both a polar and an eccentric orbit, thereby preserving dynamical signatures that may have otherwise been erased by tidal circularization. We investigate the unusual orbital architecture of GJ 436, exploring von Zeipel-Lidov-Kozai migration induced by a distant companion as a mechanism to explain the present-day orbit of GJ 436 b. Using $\sim$20 years of archival radial velocity measurements and astrometric data from the Hipparcos-Gaia Catalog of Accelerations, we constrain a potential companion to $a_{c}<5.4$ AU for $m_{c}>0.05$ $M_{Jup}$ and $a_{c}<64$ AU for $m_{c}>24$ $M_{Jup}$ in the GJ 436 system at the $2\sigma$ confidence level, providing the most stringent constraints to date. We further perform three-body hierarchical secular simulations to determine which companion configurations could reproduce GJ 436 b's present-day orbit within the observationally constrained parameter space. Our dynamical modeling favors sub-Jovian masses on orbits with $a_\mathrm{c} \gtrsim 6.8$ AU, suggesting a substellar perturber. These observational and dynamical constraints can guide future companion searches and illuminate formation mechanisms for hot Neptune desert planets on polar orbits.

Protostellar streamers are elongated structures that channel material from larger scale onto disks, influencing their physical and chemical evolution. The M512 protostar in Orion/Lynds 1641 hosts one of the most massive and extended streamer discovered so far, offering a unique opportunity to study these processes. We investigate the morphology, chemistry, and origin of this streamer,and its potential impact on the protostellar disk. Using archival ALMA observations of C18O, DCO+, N2D+, and HCO+, we compare their spatial distributions through moment maps and spatial profiles. The streamer shows clear chemical stratification: C18O lies on the western side of the protostar, N2D+ is farther out to the east, and DCO+ is in the middle. This suggests that the structure has been shaped by environmental effects rather than tracing a single coherent infalling flow, with only the densest gas near the protostar likely to accrete onto the disk. Overall, the bulk of the streamer reflects the physical and chemical imprint of the surrounding cloud, highlighting the importance of environmental shaping in interpreting streamer-disk connections and their role in disk growth.

The scalar-tensor theory of gravity, a modified gravity theory, introduces a fundamental scalar field that can serve as dynamical dark energy, driving the late-time accelerated expansion of the Universe. In this work, we analyze electromagnetic (EM) radiations arising from scalar fields and compare these features with those induced by axion-like particles (ALPs). Scalar and ALP fields couple differently to the EM field due to their distinct parity properties, $\phi F_{\mu\nu} F^{\mu\nu}$ for scalar fields and $\phi F_{\mu\nu} \tilde{F}^{\mu\nu}$ for ALPs. Building on analytical methods developed for ALPs, this work presents a theoretical feasibility analysis that demonstrates how the scalar field could produce observable EM signatures from oscillating field configurations. We also show that resonance effects can amplify the EM radiation for the scalar field under specific conditions, and that the enhancement mechanisms depend on the coupling structure and the configuration of the background magnetic field. Resonance phenomena can accentuate the differences in signal strength and spectral features, potentially aiding future observations in distinguishing scalar fields from ALPs. Although our studies apply to general scalar fields, we embed them within the framework of scalar-tensor theory and discuss the mass and coupling parameter in the context of testing modified gravity. This work provides a theoretical framework for studying generic pure and pseudo-scalar fields on an equal footing and suggests new avenues for observational tests of modified gravity scenarios alongside ALP models.

Superconducting resonator parametric amplifiers are potentially important components for a wide variety of fundamental physics experiments and utilitarian applications. We propose and realise an operating scheme that achieves amplification through the use of non-degenerate pumps, which addresses two key challenges in the design of parametric amplifiers: non-continuous gain across the amplification band and pump tone removal. We have experimentally demonstrated the non-degenerate pumping scheme using a half-wave resonator amplifier based on NbN thin-film, and measured a peak gain of 26 dB and 3-dB bandwidth of 0.5 MHz. The two non-degenerate pump tones were positioned ~10 bandwidths above and below the frequency at which peak gain occurs. We have found the non-degenerate pumping scheme to be more stable compared to the usual degenerate pumping scheme in terms of gain drift over time, by a factor of 4. This scheme also retains the usual flexibility of NbN resonator parametric amplifiers in terms of reliable amplification in a ~4 K environment, and is suitable for cross-harmonic amplification. The use of pump tones at different frequencies allows phase-sensitive amplification when the signal tone is degenerate with the idler tone. A gain of 23 dB and squeezing ratio of 6 dB were measured.

Apparent horizon (AH) finders are essential for characterizing black holes and excising their interiors in numerical relativity (NR) simulations. However, open-source AH finders to date are tightly coupled to individual NR codes. We introduce BHaHAHA, the BlackHoles@Home Apparent Horizon Algorithm, the first open-source, infrastructure-agnostic library for AH finding in NR. BHaHAHA implements the first-ever hyperbolic flow-based approach, recasting the elliptic partial differential equation for a marginally outer trapped surface as a damped nonlinear wave equation. To enhance performance, BHaHAHA incorporates a multigrid-inspired refinement strategy, an over-relaxation technique, and OpenMP parallelization. When compared to a naïve hyperbolic relaxation implementation, these enhancements result in 64x speedups for difficult common-horizon finds on a single spacetime slice, enabling BHaHAHA to achieve runtimes within 10% of the widely used (single-core) AHFinderDirect and outperform it on multiple cores. For dynamic horizon tracking with typical core counts on a high-performance-computing cluster, BHaHAHA is approximately 2.1 times faster than AHFinderDirect at accuracies limited by interpolation of metric data from the host NR code. Implemented and tested in both the Einstein Toolkit and BlackHoles@Home, BHaHAHA demonstrates that hyperbolic relaxation can be a robust, versatile, and performant approach for AH finding.

Cosmological correlation functions of inflaton and graviton perturbations are the fundamental observables of early universe cosmology and remain a primary target for observations. In this work, we ask the following question: are these observables independent of one another? We find that in the parity-odd sector of inflationary perturbation theory, the answer is a resounding no! In earlier work we derived a correlator-to-correlator factorisation formula which states that parity-odd correlators factorise into lower-point correlators under some mild assumptions on the underlying theory. In this work, we show that these assumptions are satisfied in dynamical Chern-Simons gravity where the action of minimal inflation is augmented by a coupling between the inflaton and the gravitational Chern-Simons term. Such a theory gives rise to a parity-odd trispectrum of curvature perturbations, and we show that such a trispectrum can be expressed solely in terms of the bispectrum that arises due to the minimal coupling between the inflaton and graviton, and the graviton power spectrum which receives a parity-odd correction in this theory. The trispectrum is quadratic in this mixed inflaton-graviton bispectrum and can therefore be interpreted as a ``double copy". Our final expression for the parity-odd trispectrum is a relatively simple function of the external momenta that is rational and factorised.

We perform a Bayesian analysis of the neutron star (NS) equation of state (EoS) based on a wide set of Skyrme functionals, derived from previous nuclear physics inferences. The novelty of this approach lies in starting from the full multidimensional posterior distribution of nuclear matter parameters, consistent with a comprehensive set of static and dynamic nuclear structure observables. We construct unified EoSs for $npe\mu$ matter, where the inner crust of the NS is treated using an extended Thomas-Fermi method, providing for the first time a fully consistent Bayesian treatment of the correlation of bulk with surface as well as with spin-orbit and effective mass parameters. We then employ a standard Bayesian framework to identify those EoSs that satisfy astrophysical constraints from NS mass measurements, the tidal deformability from GW170817, and NICER mass-radius observations. We also examine NS observables, such as the crustal moment of inertia, which is crucial in understanding pulsar glitches. Compared to previous works, we observe an increase in both the NS surface thickness and the crustal moment of inertia.

We review the cosmic generalisation of the virial theorem known as the Layzer-Irvine equation, also independently derived by Dmitriev and Zeldovich. This equation has been studied in the literature for dark matter-dark energy interaction models, as well as in the context of alternative theories of gravity. We discuss results from the previous scenarios and point out future directions.

We derive relativistic Maxwell-Bloch equations for potential applications in astronomical environments, where various radiative processes are known to occur, including the maser action and Dicke's superradiance. We show that for both phenomena a radiating system's response is preserved at different relative velocities between the system's rest frame and the observer, while the relevant timescales and the radiation intensity transform as expected from relativistic considerations. We verify that the level of coherence between groups of emitters travelling at different speeds is unchanged in all reference frames. We also derive relativistic versions of the maser equations applicable in the steady-state regime.

The nature of the cosmological constant is a mystery. We don't understand its quantum origin but we associate it with the actual acceleration of the universe because it is the simplest description we had until recently of the present cosmological observations. However, this may change with the next generation of experiments. If we can convince ourselves that the cosmic acceleration is not due to a constant, this would open up new fascinating avenues. By exploring the simplest cosmological model in the bulk, that of an empty and flat space with a cosmological constant $\Lambda$, we find that its holographic correspondence makes sense as a theory of fundamental quantum degrees of freedom at the boundary. Moreover, we find that an observer in the bulk, making long-range gravitational observations, cannot distinguish the acceleration induced by the cosmological constant $\Lambda$ from that induced by the thermodynamic properties of the boundary, the de Sitter horizon, strictly at the level of the background cosmology. By including matter in the bulk we extend this holographic correspondence to GREA, where the quantum d.o.f. associated with the evolving boundary of the causal horizon induces an entropic acceleration that varies in time. Upcoming surveys such as DESI, Euclid, and the Vera Rubin Observatory (LSST) will test this framework through the growth of large-scale structures in the late universe, where GREA and $\Lambda$CDM differ quantitatively.

The multi-messenger detection of the binary neutron star (NS) merger GW170817 has revolutionized the field of gravitational wave (GW) astronomy. However, several important questions remain to be answered. One of these is the nature of the compact remnant leftover by GW170817 (short- or long-lived NS versus black hole). A key goal going forward is to understand the diversity of NS-NS merger remnants, and how such diversity maps onto their viability as gamma-ray burst (GRB) central engines. Here, we present a study aimed at assessing the sensitivity of triggered searches for intermediate-duration, post-merger GWs powered by long-lived GRB remnants using networks of current and future ground-based GW detectors and the Cross-Correlation Algorithm (CoCoA). We develop a Python-based framework to efficiently estimate CoCoA distance horizons for a broad range of post merger secular bar-mode waveforms and for different GW detector networks. This framework can be used to identify the most promising regions of parameter space in which to concentrate search efforts, helping design future search strategies to optimally balance search sensitivity and related parameter space gridding schema against computational cost.

The local solar gravitational potential forms a basin for ultralight dark matter (ULDM), with discrete energy levels. Even if barely populated, it introduces a new characteristic timescale in DM dynamics. This necessitates a generalization of the notion of coherence time. We find that, at long times, the phenomenon of recoherence emerges, whereby a subcomponent of ULDM exhibits a formally divergent coherence time. The fact that this generalized coherence time can significantly exceed the naive estimate implies an enhanced sensitivity for dark matter searches that accumulate data over extended observation periods.

This paper proposes a new avenue for understanding the cosmological singularity. The standard cosmological model contains a generic initial singularity usually referred to as the {\em big bang}. Herein, we present a novel idea to extend the description of our Universe beyond this limit. The proposal relies on rewriting physics in a purely Riemannian, {\em i.e.} locally Euclidean, four-dimensional space and the emergence of Lorentzian patches owing to the interaction of all matter fields to a clock field that is responsible for a signature change. If our Universe is contained within one of these patches, the initial singularity is replaced by a smooth boundary on which the signature of the physical metric flips. In this paper, we first define the model and draw the necessary conditions on its arbitrary functions for solutions to exist. Next, we prove the existence of solutions that lead to an emergent universe with a primordial (almost) de Sitter phase. To finish, we discuss the consequences of this construction for the universe on scales much larger than our observable Universe: a large ``Euclidean sea'' in which Lorentzian islands locally emerge and host an expanding universe potentially similar to ours. While speculative, this scenario has specific features that can be tested, and the present paper sets the basis for further phenomenological investigations.