Papers for Friday, Apr 24 2026

The surface and subsurface of worlds beyond Mars remain largely unexplored. Yet these worlds hold keys to fundamental questions in planetary science - from potentially habitable subsurface oceans on icy moons to ancient records preserved in Kuiper Belt objects. NASA's success in Mars exploration was achieved through incrementalism: 22 progressively sophisticated missions over decades. This paradigm, which we call Planetary Exploration 2.0 (PE 2.0), is untenable for the outer Solar System, where cruise times of a decade or more make iterative missions infeasible. We propose Planetary Exploration 3.0 (PE 3.0): a paradigm in which unvisited worlds are explored by a single or a few missions with radically adaptive space systems. A PE 3.0 mission conducts both initial exploratory science and follow-on hypothesis-driven science based on its own in situ data returns, evolving spacecraft capabilities to work resiliently in previously unseen environments. The key enabler of PE 3.0 is software-defined space systems (SDSSs) - systems that can adapt their functions at all levels through software updates. This paper presents findings from a Keck Institute for Space Studies (KISS) workshop on PE 3.0, covering: (1) PE 3.0 systems engineering including science definition, architecture, design methods, and verification & validation; (2) software-defined space system technologies including reconfigurable hardware, multi-functionality, and modularity; (3) onboard intelligence including autonomous science, navigation, controls, and embodied AI; and (4) three PE 3.0 mission concepts: a Neptune/Triton smart flyby, an ocean world explorer, and an Oort cloud reconnaissance mission.

Recent LIGO-Virgo-KAGRA (LVK) analyses have revealed complex structure in the binary black hole (BBH) population, including distinct features in the primary mass spectrum and nontrivial spin-mass correlations. However, the phenomenological models used to capture these features often lack analytic transparency, making it difficult to isolate robust physical laws from modeling artifacts. To address this, symbolic regression is applied to the posterior inference products of the GWTC-4 catalog, discovering compact, closed-form analytic expressions for four key population relationships: (i) the merger-rate evolution with redshift; (ii) the mass-ratio dependence of the effective-spin distribution; (iii) the redshift evolution of the effective-spin distribution; and (iv) the conditional mass-ratio distributions associated with the 10 solar mass and 35 solar mass primary mass peaks. This framework successfully compresses both rigid and highly flexible models into differentiable phenomenological laws, dynamically recovering a consistent low-redshift merger-rate slope without assuming an a priori power-law form. The exact analytic derivatives provided by symbolic regression show that the mass ratio--effective spin and redshift--effective spin correlations are robustly driven by broadening of the posterior widths rather than shifts in the mean. Furthermore, qualitatively distinct functional forms for the mass-ratio distributions conditioned on the 10 solar mass and 35 solar mass primary mass peaks are identified. These closed-form expressions enable exact analytic gradient diagnostics and compact surrogate summaries, particularly for flexible numerical posteriors that are not otherwise available in low-dimensional analytic form. They also facilitate rapid downstream calculations for rate forecasting, formation channel comparison, and stochastic background estimation.

Radio continuum observations are widely used to study cosmic ray (CR) electron populations and transport processes in star-forming galaxies, but their interpretation relies on several simplifying assumptions. Here, we revisit three common assumptions: that some vertical radio profiles can be explained by CR advection alone, that radio spectra directly trace the galaxy-wide CR electron spectrum, and that bremsstrahlung and Coulomb losses are negligible for radio-emitting electrons. We model radio emission using time-dependent CR electron evolution in a magnetohydrodynamical simulation of an isolated Milky Way-mass galaxy. CR electron spectra are evolved self-consistently along Lagrangian tracer particles with the CREST framework, including injection at supernova remnants, advection with the gas, and spatially and temporally varying radiative losses. We compare these results to commonly adopted steady-state models. We find that advection-only transport in self-consistently driven galactic winds fails to reproduce the extended vertical radio intensity profiles observed in edge-on galaxies, despite reproducing the observed steepening of spectral indices with height. This is because slowly accelerating winds keep electrons in strong cooling environments for too long. Matching observed radio haloes with advection alone requires unrealistically high midplane wind velocities, implying that additional transport or re-acceleration processes are required. Although galaxy-integrated CR electron spectra at radio-emitting energies are similar across models, the resulting synchrotron spectra differ systematically because radio emission is biased toward young electrons in dense, strongly magnetised regions. Finally, we show that bremsstrahlung and Coulomb losses significantly shape radio spectra even when their loss rate is subdominant and therefore cannot be neglected.

We investigate the connection between supermassive black holes (SMBHs), their host galaxies, and large-scale dark-matter halos using broad-line X-ray AGN from the XMM--XXL and Stripe\,82X surveys, together with galaxies from VIPERS and SDSS/Stripe\,82. Building on the homogeneous host-galaxy catalogue presented in Paper~I, we test whether AGN with a given black-hole mass, $M_{\rm BH}$, inhabit different large-scale environments from non-AGN galaxies with similar host properties. We first examine the empirical $M_{\rm BH}$--$M_{\star}$ relation of the AGN sample. We find a shallow trend with substantial scatter, likely driven by flux-limited selection effects and uncertainties in virial black-hole mass estimates. The ratio $M_{\rm BH}/M_{\star}$ decreases with increasing stellar mass, and AGN lying above and below the empirical relation show different median host properties, consistent with non-synchronous SMBH and stellar growth. We then divide the AGN into two black-hole mass bins, $8.0 \le \log(M_{\rm BH}/M_\odot) < 8.5$ and $8.5 \le \log(M_{\rm BH}/M_\odot) < 9.0$, and construct galaxy control samples matched in $M_{\star}$, SFR, and sSFR using a multivariate nearest-neighbour method. From AGN--galaxy cross-correlation functions, we infer the characteristic halo masses of AGN and matched galaxies. In the lower-$M_{\rm BH}$ bin, AGN occupy halos statistically indistinguishable from those of their controls. In the higher-$M_{\rm BH}$ bin, we find a mild indication that AGN may reside in somewhat more massive halos, with a difference of about 0.4 dex, although still consistent within the uncertainties. If confirmed with larger samples, this would suggest that halo-scale processes become important mainly at the highest $M_{\rm BH}$.

GD-1 is among the longest, coldest stellar streams in the Milky Way, making it an ideal target for probing dark matter substructure through dynamical heating. We present a catalog of 608 spectroscopically confirmed GD-1 members from the first three years of Dark Energy Spectroscopic Instrument (DESI) observations. This constitutes the largest homogeneous spectroscopic sample of GD-1, doubling the number of members previously available only through heterogeneous compilations combining multiple surveys with different systematics. Using these data, we derive updated stream tracks in sky position, proper motion, and radial velocity that extend over $100^\circ$ of the stream. We apply a Gaussian mixture model to decompose the stream into a dynamically cold thin component ($\sigma_V = 2.49\pm 0.28$ km s$^{-1}$, width $= 0.23\pm0.01^\circ$) and a kinematically hot cocoon ($\sigma_V = 6.13\pm0.75$ km s$^{-1}$, width $= 2.18\pm0.17^\circ$). The cocoon contains $\sim30\%$ of members and its velocity dispersion is consistent with $\sim11$ Gyr of heating by cold dark matter subhalos. We also detect a large proper motion dispersion ($41.36\pm4.98$ km s$^{-1}$) along the stream direction in the cocoon component. This feature indicates a significant line-of-sight distance spread in the cocoon, and its origin will be further explored in a forthcoming paper. These measurements demonstrate the power of DESI spectroscopy for characterizing the multi-component phase-space structure of stellar streams and constraining small-scale dark matter substructure.

The James Webb Space Telescope (JWST) has revealed a puzzling population of massive black holes in the first billion years, many of which are over-massive compared to their hosts (obese black holes), and reside in metal-poor hosts, posing a challenge for theoretical models at these early epochs. In this work, we compare the observational properties of astrophysically-seeded black holes using the DELPHI semi-analytic model and cosmologically-seeded primordial black holes (PBHs) using the PHANES analytic model. We explore the growth of light ($\sim 100 M_\odot$) and heavy ($\sim 10^{3-5}M_\odot$) seeds through mergers and accretion (both Eddington-limited and at super-Eddington rates) in the astrophysical scenario; PBHs (seeded between $10^{0.5-6}M_\odot$) only grow through accretion at sub-Eddington rates. Comparing to observables at $z \sim 5-10$, the only model that can be ruled out is the one where we allow Eddington-limited accretion onto light seeds. The observed high values of the black hole mass-stellar mass relation ($0.3-1$) can be reproduced by both PBHs and heavy seeds accreting at super-Eddington rates. However, only the PBH and Eddington-limited heavy seeding models can simultaneously reproduce the observed black hole masses (${\rm M_{bh}}$), stellar masses ($M_*$), and extremely low host metallicities ($Z \leq 0.01 Z_\odot$) inferred at $z \sim 7-10$. Crucially, we find PBHs show decrease in the black hole mass-stellar mass ratio with increasing halo mass at all redshifts, contrary to any astrophysical black hole model. Selecting systems at $z \sim 7$ with ${\rm M_{bh}}/M_* > 0.1$ and bolometric luminosities $\sim 10^{44-46} {\rm erg~s^{-1}}$ that show a negative black hole to stellar mass ratio and reside in $10^{9-11}M_\odot$ halos offer a promising clustering-based discriminant of PBH seeding models.

Gravitational wave (GW) detector LISA will observe near-coalescence extreme mass ratio inspirals (EMRIs), which typically form in galactic central accretion disks. Gas torques on EMRI will alter its GW-driven inspiral trajectory from the vacuum expectation, leading to potentially LISA-observable GW dephasing ($\Delta\psi_{\rm gas}$). Most studies compute $\Delta\psi_{\rm gas}$ for a thin, laminar disk, with negligible flow turbulence, where the disk exerts a fairly well-understood linear torque ($T_{\rm lin}$). However, these disks must be turbulent due to magneto-rotational instability in the inner regions. Hence, we present a proof-of-concept general, agnostic prescription for the turbulent torque ($T_{\rm turb}$) acting on an EMRI by modeling it as a Gaussian distribution around $T_{\rm lin}$, based on recent advances from a global hydrodynamical (HD) study. We compute $\Delta\psi_{\rm gas}$ for the ``golden'' circular EMRI with total source mass $M=10^6~{\rm M}_\odot$ and mass ratio $q=5\times10^{-5}$ in its final four-year evolution at redshift $z=0.276$ and signal-to-noise ratio (SNR) $=50$ by varying Eddington ratio ${\rm f}_{\rm Edd}$, turbulence normalization $C$ ($=~360$ in the aforementioned HD study), disk aspect ratio $h_0$, and turbo-viscous coefficient $\alpha$ in a reasonable parameters space. We find that for ${\rm f}_{\rm Edd}\gtrsim0.3$, $C\gtrsim300$, $h_0\gtrsim0.03$, and $\alpha\gtrsim0.1$, gas-induced dephasings are unobservable if only considering $T_{\rm lin}$ but could become detectable ($\Delta\psi_{\rm gas}>8/$SNR) if EMRIs exhibit chaotic migration due to turbulent gas flow. Hence, this work motivates running MHD simulations of accretion disks with embedded LISA EMRIs in the early in-spiral phase over long enough timescales to understand the evolution of their orbital elements and the imprint of the turbulent environment on their gravitational waveforms.

By monitoring the times of arrival of radio pulses from millisecond pulsars, Pulsar Timing Arrays (PTAs) serve as unique gravitational wave (GW) laboratories in the nanohertz band. To date, the primary astrophysical sources of GWs targeted in this frequency range have been inspiraling supermassive black hole binaries (SMBHBs) on circular and eccentric orbits. In this work, we demonstrate that, thanks to the so-called pulsar term in the timing residual waveform of GW signals, PTAs can probe individual SMBHBs that merged before timing observations began. We refer to the latter as \emph{zombie binaries}. Using SMBHB population models consistent with current PTA constraints, we find that while the probability of detecting such systems in existing PTA datasets remains low, the Square Kilometer Array observatory is expected to achieve sufficient sensitivity to have a few zombie binaries with a signal-to-noise ratio exceeding 3 in its data. Although their confident identification might be challenging, this new class of PTA sources opens a novel window for studying the most massive SMBHBs in our local universe.

Relativistic effects are sensitive to subtle changes in dark energy. These effects grow on very large scales and at high redshifts, which will be the reach of upcoming surveys. We investigate these effects in both the linear and the angular galaxy power spectra in a late-time universe dominated by cold dark matter and k-essence, focusing on three core models (dilaton, tachyon, and DBI scalar fields) and contrasting their predictions with those of the concordance model. By enforcing identical present-day cosmological parameters, we isolate the imprints of k-essence dynamics and perturbations on very large scales. We found that relativistic corrections dominate on very large scales and grow with redshift, but are largely insensitive to k-essence microphysics in Fourier space, leading to strong degeneracies among the models. However, in the angular power spectrum, where line-of-sight integrals are naturally included, relativistic effects are significantly amplified, yielding better sensitivity to clustering k-essence. In particular, the tachyon exhibits clear deviations across multipoles and redshifts, with distinct imprints in the Doppler and the combined (velocity and gravitational) potentials contributions. Furthermore, our results show that neglecting relativistic corrections can lead to systematic misestimation of deviations of k-essence from the cosmological constant. Our results show the relativistic angular galaxy power spectrum as a more consistent and robust probe of ultra-large-scale physics. These findings underscore the need for full relativistic modelling in next-generation surveys that are targeting horizon-scale modes, where the imprint of non-standard dark energy is most pronounced.

Pulsar timing arrays record gravitational waves from supermassive black hole binaries at two spacetime points: an Earth term, measured when the wave passes the Earth, and a pulsar term, measured when the wave passed each pulsar at an earlier epoch. We show that a future $\mu$Hz-band detection of a nearby massive binary by a mission such as $\mu$Ares would turn PTA pulsar terms into targeted probes of binary evolution. In analogy with supernova light echoes, each pulsar term acts as a gravity echo: a dated snapshot of the binary at an earlier stage of its inspiral. Together, the $\mu$Hz Earth-term measurement and the nHz pulsar-term echoes provide a temporal baseline that neither detector could access alone. For a fiducial equal-mass binary with total mass $10^9\,M_\odot$ at 80~Mpc, we find a combined pulsar timing array echo signal-to-noise ratio of 33, with up to 24 pulsars individually resolving the signal among pulsars with 50-year baselines. The angular dependence of the single-pulsar echo sensitivity alone enables independent sky localization of the source to $\sim$10--100~deg$^2$, and the resolved pulsar-term frequencies directly measure the binary inspiral rate hundreds to thousands of years ago. With sufficient pulsar distance precision, a small set of anchor pulsars could additionally phase-connect the array and trace the post-Newtonian evolution coherently over kpc baselines. The source population required for gravity echoes is drawn from the same massive-end census responsible for the observed nanoHertz stochastic background.

While massive black hole binaries (MBHBs) merge at gravitational-wave frequencies above the pulsar timing array (PTA) sensitivity band, we show that they leave orphaned low-frequency contributions in the PTA pulsar term. Due to the light-propagation time between each pulsar in the array and Earth, the pulsar term acts as a time-delayed probe of a chirping merger with a specific frequency response determined by the direction of origin and intrinsic properties of the MBHB. We provide a detailed consideration of how such a multiband signal would manifest in a full PTA, demonstrate an approach to stack these orphaned pulsar terms across the array, and discuss prospects for an archival, multiband search in conjunction with MBHB mergers observed in astrometric data or spaceborne interferometers like LISA.

We compute the abundance of cosmological gravitational waves produced during the evolution of an ultralight vector (spin-1) dark matter field. A homogeneous background vector field breaks spatial isotropy, requiring a Bianchi I geometry and inducing a mixing between the scalar, vector, and tensor perturbation sectors. We derive the perturbation equations in this background and show that, as a consequence of this mixing, scalar perturbations act as a source of tensor modes, generating a stochastic GW background. The production and cosmological evolution of these gravitational waves are implemented in a modified version of CLASS, from which we obtain the present-day spectrum.

In core-collapse supernovae (CCSNe), deleptonization normally favors $\nu_e$ over $\bar{\nu}_e$ emission. However, lepton-number emission self-sustained asymmetry (LESA) can make the energy-integrated emitted lepton-number flux negative along some directions. We derive a simple diagnostic for this transition and test it in 33 independent 3D CCSN simulations: 25 Princeton/Fornax models ($8.1$--$100\,M_\odot$) and 8 Garching models, including non-, slow-, and fast-rotating $15\,M_\odot$ cases. Of 23 non-black-hole-forming Princeton models, 22 cross the threshold, with median onset $t_c=225\,\mathrm{ms}$, IQR $162$--$264\,\mathrm{ms}$, and cross-model scatter $\mathrm{CV}=18.6\%$. Full-sky flux-sign searches show that the threshold identifies the anti-LESA-pole transition, distinguishing the global LESA-driven crossing from early localized turbulent crossings. The fast-rotating Garching $15\,M_\odot$ model, where rapid rotation suppresses the LESA dipole, is correctly classified as a non-crosser without using any rotation parameter. Both black-hole-forming Princeton models cross near $250\,\mathrm{ms}$ post-bounce and remain above threshold for $1807$ and $2463\,\mathrm{ms}$ before collapse. Thus, in the next nearby CCSN, the emitted $\bar{\nu}_e$ energy flux may exceed the $\nu_e$ flux along some lines of sight. Such directions may also correlate with sustained fast flavor instability, although testing this requires local phase-space distributions or dedicated linear stability analysis. The relevant quantity here is the energy-integrated emitted flux field, i.e. a luminosity difference per steradian, not a neutrino number flux.

The internal structure of dark matter halos on sub-galactic scales remains a key open question, particularly in the context of the core-cusp problem. Ultra-faint dwarf galaxies (UFDs), owing to their extreme dark matter dominance, provide a promising laboratory to probe these density profiles through stellar tracers. In this work, we assess the capability of the Chinese Space Station Telescope (CSST) to detect and characterize wide binary stars in the nearby UFD Segue 1, using mock observations. We generate mock binary populations based on our existing $N$-body simulations and incorporate realistic CSST observational conditions, including the expected deep-field limiting magnitude ($g \sim 27.5$ mag) and a photometric completeness of approximately $90\%$. The two-point correlation function (2PCF) of stellar pairs is used as a statistical tool to recover the binary fraction under these assumptions. We find that CSST can robustly detect wide binaries at the $3\sigma$ level for binary fractions as low as $f_b \gtrsim 0.01$, provided a stellar sample size of $N_{\mathrm{star}} \gtrsim 2300$. However, distinguishing between cusped and cored dark matter profiles is significantly more demanding, requiring $N_{\mathrm{star}} \gtrsim 6000$ and $f_b \gtrsim 0.1$ within $\sim 40\mathrm{kpc}$.

Various processes can induce long-lived overdense rings and arcs in protoplanetary and AGN accretion discs, such as the accumulation of gas at the outer edge of the dead zone, or the infall of material. Using the local approximation of dynamical friction, we investigate the orbital evolution of a low-mass highly-eccentric point-mass accretor (perturber) embedded in an isothermal disc hosting a density ring. We specifically consider the regime in which the eccentricity exceeds four times the disc aspect ratio. For prograde perturbers, orbits that cross the ring progressively circularize while their semi-major axes converge toward the ring radius. As a result, perturbers accumulate, forming a population ring superimposed on the gaseous ring. The ring therefore acts as a migration trap for these eccentric orbits. We also find that prograde orbits tangent to the ring, either at apocentre or pericentre, remain tangential throughout their evolution; perturbers confined to these trajectories experience the highest accretion rates. In contrast, retrograde perturbers always migrate inward. Once the semi-major axis becomes smaller than the ring radius, the eccentricity grows, but not enough for the orbit to intersect the ring again. We also discuss how feedback effects, such as jet launching and thermal torques, could modify the effective forces acting on the perturbers.

The Lunar Surface Electromagnetics Experiment (LuSEE-Night) is a joint NASA-DOE-ESA low-frequency radio telescope that will reach the lunar far side in 2027. The unknown dielectric properties of the subsurface at the LuSEE-Night landing site impose the most significant limitation for precision instrument calibration, as reflections from the lunar subsurface can change the primary beam at the 10-20% level. Simulations of these effects have provided insight and concern, showing that the lunar subsurface modeled as a lossy dielectric can absorb a large amount of the power of the sky signal. While this absorption may not strongly impact the signal-to-noise ratio in a sky-noise-dominated regime, it could complicate the beam pattern and make the signal more difficult to model and interpret. We have simulated the far-field properties of the LuSEE-Night beam for varying dielectric profiles of the lunar subsurface. We find that varying the properties of the lunar subsurface has the most significant impact around the antenna resonance, impacting its amplitude, position and width. Conversely, changing the properties of the foreground impacts the data across the band. We use a Bayesian inference pipeline to jointly estimate parameters of a galactic foreground model and dielectric properties of the lunar subsurface around the LuSEE-Night landing site and find that parameters of both the galaxy and subsurface properties can be estimated jointly. While the modeling is somewhat idealized, we believe that the results are largely robust owing to the fact that spectral variations for plausible subsurface and galaxy models have very different spectral signatures.

We present SDO2IRIS$^2$: a visual transformer model that translates a combination of images of the chromosphere and transition region (TR), observed by AIA, and a line-of-sight magnetogram, provided by HMI, into temperature, line-of-sight velocity (v$_{los}$), velocity of the turbulent motions (v$_{turb}$), and electron density (n$_{e}$) in the chromosphere. Using the thermodynamic variables obtained from the inversion of the chromospheric lines Mg II h&k, observed by IRIS, as the target of the model, and the intensity images in the chromosphere and TR, and the photospheric magnetogram as the input, the predicted T and n$_{e}$ show a strong correlation ($\approx 0.80$) for $\approx$80% of the test inverted data, a moderate-to-strong correlation ($\approx0.63$) for 70% of the v$_{turb}$ of the target test inverted data, while for the $v_{los}$, the correlation is weak. Therefore, the predicted values by SDO2IRIS$^2$ may be used as an estimation of the thermodynamics in the chromosphere, either as a stand-alone result or as complementary information to other chromospheric data observed simultaneously. The execution time employed by SDO2IRIS$^2$ to obtain the thermodynamic values in the chromosphere is of the order of a few minutes, being $\le10$ minutes when using a CPU, and $\le5$ minutes when using a GPU. SDO2IRIS$^2$ opens a new avenue for the use of SDO data thanks to the inversions provided by IRIS observables.

Models predict that chemical enrichment and gas redistribution should proceed rapidly once star formation begins, yet direct observational constraints at the earliest cosmic epochs have been scarce. Here we present evidence that metal-enriched gas in multiple ionic phases was already present around galaxies before the midpoint of cosmic reionization. Using JWST/NIRSpec rest-frame ultraviolet spectroscopy of three galaxies at redshifts $z=7.2-9.3$, we detect blueshifted metal absorption in all three systems; across the sample, the detected transitions span neutral, low-ionization, and high-ionization species, including O I, Si II, C II, Si IV, and C IV. These absorption features show velocity offsets of order $|\Delta v| \sim 50$--$250\,\mathrm{km\,s^{-1}}$, predominantly blueshifted relative to the systemic redshifts of the host galaxies derived from nebular emission lines. This ionic coexistence within a broadly shared velocity structure, together with the observed equivalent-width ratios, is consistent with outflowing or otherwise kinematically disturbed galaxy-associated gas, similar to that seen at lower redshift. The observations therefore indicate that metal-enriched gas associated with galaxies was already kinematically disturbed at very early times, requiring rapid metal production in the early generations of stars. These results show that key conditions for baryon cycling were established in at least a subset of luminous galaxies within the first several hundred million years of cosmic time, well before the completion of reionization.

We investigate whether the observed radio-infrared-$\gamma$-ray correlation in star-forming galaxies is a geometric effect rather than a signature of local cosmic-ray (CR) calorimetry. Using the GALPROP framework, we generate synthetic observations for external viewers from a grid of 3D Milky Way models with varied CR source, gas, interstellar radiation, and magnetic field distributions, all normalised to reproduce local CR data. We find that a tight, quasi-linear correlation arises naturally from line-of-sight integration through the extended, radially-structured disc, even when local calorimetry is absent. The correlation's properties depend strongly on viewing geometry, preserving its form under moderate inclination but breaking down in edge-on views where galactic components are stratified. We conclude that the correlation is primarily an emergent property of geometric projection, not local physics. This implies that its scatter is likely not random noise but a diagnostic of underlying galactic structure and viewing angle.

Calcium-aluminum-rich inclusions (CAIs) in carbonaceous chondritic meteorites are the oldest relics in the solar system. Notably, their radiogenic age feature a brief (100 kyr) condensation episode. In contrast, the reservoirs of the short-lived isotopes in CAIs, presumably supernovae or asymptotic giant stars, pollutes star-forming regions in giant molecular cloud complexes (GMC) over much longer (Myr) duration. Through a series of numerical simulations, we show here the possibility that, within an extended region (2$\sim$3 AU), nearly all ``pre-solar'' CAI-loaded grains in the infall clouds were sublimated and re-condensed during the early ($ \lesssim 10^5$ yr) infall and formation of class-0 disks. We adopt a set of initial conditions from a previous hydrodynamic simulation of the collapse of GMC and the formation of young stellar clusters. We analyze the evolution of the disk's thermal distribution and dynamical structure resulting from the interaction between circumstellar disks and infalling gas. Our follow-up simulations, with much higher resolution, show significant and rapid changes in the disk orientation and morphology due to the dynamic infall of external streamers. Warps and global spiral density waves commonly appear. They lead to intense dissipation which heats the gas to sufficiently high temperature to sublimate prior-generation CAIs. This solid-to-gas phase transition is followed by subsequent cooling and re-condensation. The CAI contained in the meteorites today could be the relics of the last episode of major infall onto class 0 disks.

The estimation of stellar atmospheric parameters for large-scale samples, particularly metal-poor stars, is a cornerstone of Galactic archaeology. In this work, we optimized a photometric filter design tailored to measuring stellar metallicities for very metal-poor stars with [Fe/H]$< -1$.The optimal configurations consist of a central wavelength $\lambda_{\rm c}$ = 3960 Angstrom with a bandwidth $\Delta\lambda$ = 80 Angstrom for giant stars, and $\lambda_{\rm c} $= 3920 Angstrom with $\Delta\lambda$ = 80 Angstrom for dwarf stars. By applying these optimized filters to synthetic photometry derived from Gaia XP spectra, we inferred metallicities for both populations. Both internal and external validations demonstrate high precision across a wide metallicity range: 0.18-0.19 dex for $-2 \le \rm [Fe/H] \le -1$, 0.23-0.33 dex for $-3 \le \rm [Fe/H] \le -2$, and approximately 0.39 dex for the most metal-poor regime, successfully extending down to $\rm [Fe/H] \approx -4$ for giant stars, $\rm [Fe/H] \approx -3.3$ for dwarf stars. Finally, we present a catalog of approximately 14.5 million metal-poor stars with robust $\rm [Fe/H]$ measurements, along with more than ten thousand red giant ultra metal-poor candidates with $\rm [Fe/H] < -4.0$, providing a valuable resource for exploring the early formation and chemical evolution of the Milky Way.

Black-hole accretion systems exhibit a characteristic coexistence of activities: broad-band X-ray variability, hot coronae, wide-angle winds, and both steady and discrete jets. This coexistence suggests a persistently time-dependent magnetic background in which noisy fluctuations and explosive release are both essential. In this paper, we connect them all to intermittent magnetic reconnection and propose a Synchronized Spin Model (SSM) in which multiple local dynamos in a rotating accretion flow are represented as interacting macro-spins. Their synchronization, partial synchronization, excursion, and reversal define a compact set of collective variables that organize both timing statistics and large-scale morphology. In this picture, multiscale magnetic reconnection sustains coronal heating, flares, intermittent outflows, and discrete jet activity, while the same synchronization dynamics produce amplitude modulation and demodulation, providing a route to $1/f$-like variability, rms--flux/Taylor-like scaling, and approximately log-normal statistics of the demodulated envelope. We further argue that, although the continuous flux distribution in black-hole systems is more naturally discussed in multiplicative or log-normal terms, broader event-catalog statistics remain useful for describing suitably defined burst hierarchies, particularly by analogy with solar and stellar flare systems. The hard/soft cycle of X-ray binaries is then interpreted as motion through magnetic state space.

Be stars are rapid rotators generally produced by binary interactions. The single Be stars in the observations pose challenges to the Be star formation theory. In this paper, we propose a new pathway for the formation of single Be stars, in which the Be star is taken as the ejected companion star from a Type Ia supernova (SN Ia) explosion. Recent numerical simulations suggest that explosive oxygen burning, initiated via the convective Urca process in certain helium (He) stars near the Chandrasekhar mass limit, can set off a SN Ia. Based on this proposition, we further demonstrate that about $0.4\%$ of He star + main-sequence (MS) star binaries may evolve into single Be stars, where the MS star is spun up due to the mass accretion from the He star, and then the He star explodes as a SN Ia. We employ detailed binary evolutionary simulations and find the parameter space that would produce single Be stars via the SN Ia channel. Around $22\%$ of Be stars from the SN Ia progenitor channel exhibit peculiar tangential velocities exceeding $24\ \rm km/s$, classifying them as runaway stars. This suggests that the SN Ia channel plays a meaningful role in forming single Be stars, particularly within the runaway star population.

Most Active Galactic Nuclei (AGN) are Radio Quiet, with radio emission that may arise from star-formation activity, AGN-driven winds, weak jets, and coronal activity. Disentangling these mechanisms is challenging and requires detailed multi-wavelength investigation, but it is crucial for quantifying AGN feedback in galaxy evolution. We present a detailed radio investigation of 21 X-ray selected AGN in the Supermassive Black Hole Winds in X-Rays (SUBWAYS) sample (log Lbol = 44.9-46.3 erg/s, z=0.1-0.5), selected to systematically search for Ultra-Fast Outflows (UFOs). UFOs are detected in 30% of the targets, making the sample particularly well-suited for investigating the role and signatures of multi-scale outflows at different frequencies. We build the radio SED of the sources complementing our proprietary data, collected with the JVLA at 1.5 and 6 GHz, with images from LoTSS and other publicly available radio surveys between 150 and 1400 MHz. We investigate the role and occurrence of the aforementioned mechanisms, with particular interest in outflows and their possible relation with UFOs. We combined information on spectral indices, luminosities, and morphologies of the radio emission with properties derived in other wavebands, such as Star Formation Rate, X-ray luminosity, Eddington ratio or the UFO kinetic luminosity. All the sources are detected and are mostly consistent with RQ AGN. For 80% of the sources the data suggest the presence of an outflow (wind or weak jet). Interestingly, our results indicate that AGN with UFOs tend to have larger radio extension and a steep radio spectrum consistent with outflows. Moreover, the radio emission of the 6 UFO hosts is consistent with predictions from wind-driven shock models, possibly indicating a direct connection between the two phases. Alternatively, this may reflect physical conditions favouring the rise of both phenomena.

We investigate electron--positron pair creation through the Schwinger process in newborn magnetars with millisecond spin periods and surface dipole fields close to or above the QED critical field, $B_{\rm Q} = 4.414\times10^{13}\,\mathrm{G}$. In the unscreened field scenario, we derive the analytical global pair creation flux and recast it into a compact form with accurate analytic approximations. For a fiducial model with $B_{\rm p} = 10^{14}\,\mathrm{G}$ and $P_0 = 1\,\mathrm{ms}$, the Schwinger channel exceeds the classical Goldreich--Julian particle supply by many orders of magnitude and becomes the dominant source of charges at the earliest stage of the magnetar. The associated discharge removes about $90\%$ of the initial rotational energy within 30 ms, suppresses the gravitational-wave loss channel, and implies that the observable millisecond phase is extremely short in this unscreened scenario. The rapid energy release over such a short timescale may also provide a viable power source for astrophysical transients. Extending the same fiducial model to $10^4\,\mathrm{yr}$ gives spin periods of order seconds, linking newborn millisecond magnetars to the mature magnetar population.

We report the closed-form expression for Hill's surfaces in the circular restricted three-body problem. The solution $\phi(r,\theta)$, derived in the primary-centric spherical coordinate system, is deduced from a cubic equation delivering at most two roots on each side of a separatrix. The famous patterns (tadpole, horseshoe and peanut shapes, Roche lobes and Hill's quasi-spheres) are exactly produced.

We present a new implementation of the SFUMATO code, called SFUMATO#, for solving self-gravitational radiation hydrodynamics problems using adaptive mesh refinement (AMR) with the CUDA/HIP programming frameworks. The code incorporates a multigrid solver for self-gravity, radiation transfer with M1 closure and reduced speed of light approximation, non-equilibrium chemistry, thermal evolution, and sink particle schemes. We develop new non-equilibrium chemistry and thermal solvers based on a linearized implicit method, whose accuracy is validated through a series of test problems by comparison with solutions obtained using the Newton-Raphson method. By incorporating the heat capacity of dust grains, the dust temperature can be evolved without iterative energy-balance calculations. From the perspective of computational cost, we demonstrate that adopting an increased pseudo dust heat capacity accelerates the chemistry solver while preserving accuracy, even when the value is increased by up to three orders of magnitude relative to the realistic value. In addition, we perform a suite of test problems to confirm the validity of the other components of our implementation. The code supports multi-GPU execution via MPI-based parallelization. We measure the strong-scaling performance of the hydrodynamics and self-gravity solvers on both uniform and AMR grids, as well as the overall code performance using a giant molecular cloud simulation. We find that the computational cost of the self-gravity solver increases with the number of MPI processes, indicating that efficient parallel performance is achieved only when the number of devices is chosen such that the cost of the self-gravity solver remains comparable to that of the other components.

White-light flares (WLFs) are usually tend to be those very large flares. Nevertheless, several small and compact WLFs have been reported and thought to be produced by low-height magnetic reconnection. However, whether low-height magnetic reconnection can efficiently produce WLFs remains unclear. For the first time, we conduct a statistical study of the WL emission in compact flares to address this question. Using over a decade observations from the \textit{Solar Dynamics Observatory} (SDO), we identify 28 compact flares, including 19 C-class and 9 B-class flares. We find these compact flares can be classified into three types based on the magnetic configuration of the flare, corresponding to the U-shape loop (type I), the flux emergence near sunspot (type II), and the fan-spine like structure (type III). For each type, the flares numbers are 9 (7 C-calss and 2 B-class), 9 (3 C-calss and 6 B-calss) and 10 (9 C-calss and 1 B-calss), respectively. We find the occurrence rate of WLFs in compact flares is $\sim60.7\%$ (17/28), and for the C-class the rate can be up to $\sim89.5\%$ (17/19). No WLF was found in B-class compact flares. The occurrence rates for three types are $\sim77.8\%$ (7/9), $\sim11.1\%$ (1/9) and 90\% (9/10), respectively. And for the C-class flares, the occurrence rates for three types are 100\% (7/7), $\sim33.3\%$ (1/3) and 100\% (9/9), respectively. Our results suggest type-I and type-III compact flares are more likely to produce WL emissions.

We present optical B and V band photometry of GP Com - an ultracompact binary consisting of an accreting white dwarf and helium secondary component. Our data set contains 7.7 hours observations in V band with the 2.0m telescope and 2.9 hours simultaneous observations in B and V bands with the 1.5m telescope of the Rozhen National Astronomical Obsevatory, Bulgaria. The observations cover of about 13 orbital periods. We find an orbital modulation with amplitude 0.04-0.05 mag in B and V bands. Adopting that it is due to a bright spot, we estimate its temperature 19700 \pm 3000 K. We estimate mass accretion rate onto the white dwarf of about 2.10^{-12} M_sun/yr, consistent with the predicted rate for a cool donor. The data are available on Zenodo: this http URL.

Context. For stellar evolution models we rely on mass-loss rate prescriptions that show features that lack empirical confirmation, such as the bi-stability jump. This jump is an increase in mass loss in the decreasing temperature regime Teff 28-21 kK. Although papers compared empirical results to prescriptions,a large observational sample of B stars for which the wind has been scrutinised over different metallicities is still lacking. Aims. By modelling of both ultraviolet (ULLYSES) and optical (XShootU) spectra, we determined the stellar and wind parameters, of 24 SMC B stars ranging in Teff from 13 to 29 kK. By combining this sample with LMC studies, we compared the wind behaviour of B stars in two different metallicity regimes. We compared our results to common mass-loss rate prescriptions to test features present in these models and their metallicity dependence. Methods. We have used the model atmosphere code fastwind and the genetic algorithm code Kiwi-GA to fit the UV and optical spectra simultaneously. This allows us to determine wind properties including clumping behaviour. Results. The metallicity trends present in the mass-loss prescriptions (Z^(0.41-1.4)) explored here overestimate the empirical metallicity dependence in the B-star regime, which appears very weak. We do not find an increase in mass-loss rate at approximately spectral type B1. We show that on average 40% of the wind mass is located in the wind medium between the clumps. We compiled a sample of 80+ O and B stars in the SMC and LMC. From a comparison we find a clear difference in O- and B-type metallicity dependence. Conclusions. The lack of a bi-stability jump in the B-star regime and a weak metallicity dependence for the same stars offers new empirical constraints to models of line-driven winds. As differences between these models are large (reaching orders of magnitude) such constraints are much needed.

Non-conservative mass transfer plays a central role in close-binary evolution, yet its effects on mass-transfer stability are uncertain. One widely adopted prescription, isotropic re-emission, is often assumed to promote stability compared to conservative mass transfer. We investigate the impact of isotropic re-emission on the stability of mass transfer in binaries with radiative envelopes that undergo delayed dynamical instability (DDI). We assess whether simplified criteria used in binary population synthesis codes accurately capture stability boundaries under varying mass-transfer efficiencies. We perform detailed stellar evolution calculations for a set of representative binaries undergoing DDI. Varying the mass-transfer efficiency beta, we track the onset of instability and quantify the corresponding critical mass ratio. We compare our results with predictions from the commonly used zeta-method, which is based on mass-radius exponents. We find that a lower mass-transfer efficiency destabilizes mass transfer in DDI systems, whereas the zeta-method predicts that isotropic re-emission should stabilize it. The discrepancy arises because the zeta-method fails to capture the full evolution of the orbit and mass ratio during pre-instability mass transfer. In some cases, the critical mass ratio is underestimated by nearly a factor of two. Our findings show that isotropic re-emission can reduce, rather than enhance, DDI stability, underscoring the limitations of using fixed critical mass ratios and zeta-based criteria. This highlights the need for calibrated prescriptions that capture the time-dependent evolution of mass ratio and orbital separation, with direct implications for modelling X-ray binaries, symbiotic stars, and double white dwarfs, including their transient rates and delay-time distributions.

We investigate the viability of Tsallis holographic dark energy (THDE) models, focusing on the role of the infrared (IR) cutoff in the growth of cosmic structures. Considering two commonly used choices of the cutoff, the particle horizon and the future event horizon, we analyze the evolution of linear matter perturbations and compute the growth factor, growth rate, and the observable $f\sigma_8(z)$. These predictions are compared with observational data from redshift-space distortion measurements. We find that the growth history is highly sensitive to the choice of IR cutoff. Models based on the future event horizon are consistent with observational data and can provide a fit comparable to, or slightly better than, the $\Lambda$CDM model for suitable values of the Tsallis parameter $\delta$. In contrast, models constructed using the particle horizon generally fail to reproduce the observed growth of structure. These results demonstrate that the viability of THDE models depends crucially on the choice of IR cutoff and highlight the importance of structure formation as a stringent test of generalized holographic dark energy scenarios.

We present ultra-deep ($\simeq 20-30$ hours), rest-frame UV spectroscopy with NIRSpec/JWST of two UV-bright galaxies at $z\sim 8.7$, CEERS-1019 and CEERS-1025 ($Z_{\rm neb} \simeq 0.1 Z_{\odot}$). The spectra reveal exceptionally strong P-Cygni profiles in wind lines (NV $\lambda$1240 and CIV $\lambda$1550) and significant broad and strong HeII $\lambda$1640 emission ($\rm EW\simeq 2-4$ A). We compare the observations with synthetic stellar population models at $Z_{\star} \simeq 0.1 Z_{\odot}$, both including and excluding very massive stars (VMS). Models including VMS provide a markedly improved fit to the data relative to non-VMS models ($\Delta$AIC and $\Delta$BIC $> 70$), which fail to reproduce the observed strengths of the wind features. A comparison with empirical spectra of VMS-dominated systems in the local Universe further supports this interpretation. The best-fit VMS models imply extremely young ages of the stellar populations ($\simeq 1.5-2.0$Myr) and high ionizing photon production efficiencies ($\log \xi_{\rm ion} [\rm Hz erg^{-1}] \gtrsim 25.8$), exceeding those inferred from non-VMS models by $\sim 0.1-0.2$ dex. These results provide evidence for an overabundance of VMS at high-$z$ with an IMF extending well beyond $100 M_{\odot}$, and highlight their potential role in shaping the rest-frame UV spectra, chemical enrichment, and ionizing output of galaxies in the early Universe.

We present the first high-resolution X-ray spectral analysis of Cygnus X-1 using XRISM. The observation was carried out from April 7 to 10, 2024, covering the orbital phase range 0.65--0.17 during its low/hard state. Taking advantage of the exceptional energy resolution of the Resolve instrument, we examined highly ionized iron absorption lines and characterized the ionization states, column densities, and line-of-sight velocities of the absorbing plasma. Spectral analysis revealed an ionization parameter of approximately 3, column densities of a few times 10^21 cm^-2, and a blueshifted velocity of approximately 100 km s^-1. The observation was divided into two phases: before and after orbital phase phi_orb = 0.9, corresponding to non-dipping and dipping intervals. While only weak absorption features were present before phi_orb = 0.9, strong absorption by He-like and H-like Fe appeared during the dipping phase. We measured equivalent widths of 2.3 eV, 0.4 eV, and 1.2 eV for He-like Fe K-alpha, and H-like Ly-alpha1 and Ly-alpha2, respectively, demonstrating the capability of XRISM Resolve to securely detect narrow absorption features of only a few eV. These measurements trace the motion of the absorbing material and offer insight into the kinematics and spatial distribution of the wind in the vicinity of the black hole. These findings enhance our understanding of wind-fed accretion in Cygnus X-1 and highlight the importance of continued high-resolution X-ray observations to further constrain the physical properties of winds and accretion flows in high-mass X-ray binaries.

JWST has revealed a population of super-luminous early galaxies with a volume density in excess of most expectations. The spectra reveal diverse properties: while some reveal strong emission lines characteristic of galaxies in the midst of strong bursts, others show weak emission lines that could reflect old stellar populations, large escape fractions, or post-burst star formation histories. Through the JWST Cycle 4 large program SPURS, we have obtained ultra-deep (29 hr) rest-frame UV spectroscopy of a z=9.3 super-luminous ($M_{\rm UV}=-21.66$) galaxy with large assembled stellar mass (1.6$\times$10$^9$ $M_\odot$) and extremely weak emission lines (H$\beta$ EW $\approx25$~Å). The strong stellar wind features and rest-optical line ratios suggest the galaxy is already significantly enriched, with a metallicity of 0.4--0.7~Z$_\odot$. The interstellar absorption lines reveal outflows ($v\simeq -161$~km~s$^{-1}$) with a large neutral gas covering fraction, suggesting that the weak emission lines are not due to large escape fractions. The combination of the Balmer break, weak emission lines, and stellar wind features constrains the star formation history, indicating a recent burst of star formation lasting 10--20 Myr followed by a downturn over the last 10~Myr. The observations suggest that $z\gtrsim 9$ weak emission line galaxies such as this source can be explained by stochastic star formation, provided that the downturns in star formation are recent (i.e., <10 Myr prior to observation). The ultra-deep grating spectrum enables the IGM damping wing to be characterized, decoupling the effects of local absorption. The smooth Ly$\alpha$ break indicates that this source, one of the most massive galaxies known at z>9, is likely situated in a small ionized bubble ($0.29_{-0.09}^{+0.11}$~pMpc), as is common at large neutral hydrogen fractions ($\bar{x}_{\rm HI}=0.81_{-0.21}^{+0.14}$).

We demonstrate that discrete self-similarity (DSS), originally discovered in the collapse of a massless scalar field in an asymptotically flat system, survives in primordial black hole (PBH) formation within an expanding cosmological background. Using fully relativistic numerical simulations of massless scalar-field collapse in an Friedmann-Lemaître-Robertson-Walker universe, we resolve the critical regime down to $|p-p_c|\sim 10^{-8}$, where $p$ and $p_c$ respectively are a parameter of the family of initial data and its threshold value, and find clear log-periodic oscillations in the PBH mass scaling relation. The detailed structure of these oscillations differs from that previously reported in the asymptotically flat case, exhibiting a more pronounced asymmetry between peaks and troughs. Analyzing two distinct families of initial data (Gaussian and piecewise rational curvature profiles), we find critical exponents and DSS periods that differ slightly but are broadly consistent within uncertainties. The presence of DSS implies characteristic log-periodic modulations in the PBH mass spectrum, with potential consequences for PBH abundances and the spectrum of induced gravitational waves.

We present the first high-resolution spectroscopy of the black hole high-mass X-ray binary Cygnus X-1 with XRISM, including orbital-phase-resolved analyses and tentative evidence of short-term variability in the Fe-K band on second timescales. Using data from the Performance Verification phase in April 2024, we analyzed spectral variability across orbital phases with the Resolve microcalorimeter and the Xtend CCD imager. The unprecedented resolution of Resolve reveals variability in highly ionized Fe absorption lines. The absorption features show orbital-phase-dependent variability in column density, ionization state, and blueshifted velocity, suggesting structural variations in the focused stellar wind along the line of sight. We also find indications of subtle broadening of the neutral Fe emission profile. In addition, intensity-sorted spectroscopy during dip phases suggests possible variability on timescales of a few seconds in the absorption features, consistent with cooler, denser and lower-ionized gas clumps. Although the statistical significance is limited, these results hint that the stellar wind and the X-rays from the accretion disk around the black hole may interact on timescales as short as a few seconds. These XRISM results constrain wind-fed accretion in Cyg X-1 and highlight Resolve's capability to probe plasma environments in high-mass X-ray binaries.

Context: ASPIICS is a giant-baseline visible light solar coronagraph, which relies on the millimetric positioning performance of the precision formation flying Proba-3 mission of the European Space Agency. Proba-3 was launched on 5 Dec 2024, and since then ASPIICS observes the solar corona with the field of view (1.1-3) R_sun. Aims: Diffraction, in particular diffraction of solar disk light on the external occulter, is known to provide a major source of straylight in coronagraphs. We aim to analyze diffracted light visible in ASPIICS images, compare it with the analytical-numerical diffraction model reported earlier, and fine-tune the model. Methods: We compare diffraction effects visible in ASPIICS data with simulated diffraction images; in particular, we compare the geometrical properties and the radiometric signal. The properties of the diffraction described in previous works suggest how to fine-tune the model in order to achieve a better correspondence with the observations. Results: Early ASPIICS observations, where diffraction is pronounced, fully confirm all the qualitative properties of diffracted light suggested by the model. After fine-tuning of the model we see quantitative correspondence of the level of 30\% -- 50\%, depending on the configuration. Conclusions: The performed analysis allows (a) to validate our analytical-numerical model and justify the assumptions, and (b) to estimate the contribution of the diffracted light in the ASPIICS images. In the majority of the field of view the diffracted light is two orders of magnitude below the coronal signal.

The mass-loss process of red supergiant (RSG) and asymptotic giant branch (AGB) stars and its relation to variability is poorly constrained. We study two evolved stars, the Mira-type AGB star R Car and the extreme RSG VX Sgr. Our sample comprises 54 VLTI-GRAVITY snapshots taken over 7 years, being the largest VLTI time-series dataset to date. We determine the angular diameter as a function of time. The radii of the photosphere ($R_{\star}$) and atomic atmospheric layers are variable and relate to the light curve with phase shifts, showing a maximum radius near visual brightness minima. The more extended CO layers show longer, irregular periods and maximum extensions of $\sim 1.3-1.7 \: R_{\star}$ for R Car, and of $\sim 1.5-2.2 \: R_{\star}$ for VX Sgr. Comparison with CO5BOLD simulations revealed a similar behavior. Furthermore, during 2020-2021, VX Sgr exhibited an extreme mass-loss event similar to that of Betelgeuse, preceded by two strong shocks and culminating with the extreme expansion of H$_2$O and CO layers, both up to $\sim 2.2 \: R_{\star}$. During this event, we detected Brackett $\gamma$ and Balmer emission lines, both of which are signatures of a shock propagating through the atmosphere. The Mira R Car showed a photospheric radius $R_{\star} = 280 \pm 25 \: \rm R_\odot$, with a fundamental mode (FM) pulsation amplitude $\sim13 \%$ of $R_{\star} $. During its active cycle, the RSG VX Sgr showed $R_{\star} = 1556 \pm 110 \: \rm R_\odot$ with FM amplitude $ \sim13 \%$ of $R_{\star} $, the same as R Car. During its quiescent cycle, it showed $R_{\star}= 1456 \pm 108 \: \rm R_\odot$ and low-amplitude pulsations near the first overtone, only $\sim4 \%$ of $R_{\star} $. This supports a steady mass loss for Miras related to stable, large-amplitude FM pulsation, whereas the mass-loss process for RSGs may be dominated by extreme events connected to changes in the pulsation mode.

We present multi-wavelength study of the $\gamma$/X-ray transient EP250416a (also designated GRB 250416C), triggered by the Einstein Probe (EP) Wide-field X-ray Telescope and also by SVOM and Konus-Wind. Observations spanning the gamma-ray, X-ray, and optical bands facilitated detailed analysis of the burst's prompt emission, afterglow evolution, and physical origin. EP250416a exhibits a burst duration of 30 s in X-ray and 17.7 s in gamma-rays, with joint spectral fitting of 0.5-5000 keV data gives $E\rm_{peak}=342_{-232}^{+90}$ keV. Optical spectroscopy of the afterglow, acquired with the Gemini Multi-Object Spectrograph (GMOS) on Gemini South, yielded a redshift of $z=0.963$. Accounting for the measured redshift, the isotropic energies are $E\rm_{X,iso}=2.7_{-0.5}^{+0.9}\times10^{50}$ erg and $E\rm_{\gamma,iso}=7.34_{-2.1}^{+5.1}\times10^{51}$ erg, aligning with the Amati relation for long GRBs. The fluence ratio $\rm S(25-50~keV)/S(50-100~keV)=0.78_{-0.15}^{+0.1}$ classifies EP250416a as an X-ray rich (XRR) GRB. The X-ray afterglow shows an initial shallow decay ($\alpha \approx -0.5$) transitioning to a canonical decay phase ($\alpha \approx -1$), with a very late jet break at $t\sim 1.5\times 10^6$ s, corresponding to a jet half-opening angle of $\theta _j=10.6_{-1.8}^{+1.9}$ degrees. EP250416a is optically dark, as it shows only a faint $r$-band detection ($r=24.16$ mag) from Gemini South-GMOS and a low optical-to-X-ray spectral index $\beta_{\rm OX} = 0.3$. This may be attributed to significant host-galaxy extinction, with a required $A_V^{\text{host}}=5.5\ \text{mag}$ derived from the extinction curve model.

The distribution of galaxies, halo abundance, and peculiar velocities are influenced by non-linear gravitational interactions, making the study of non-linear evolution crucial for accurate cosmological predictions. We explore these aspects using N-body simulations. Theoretical models of the halo mass function (HMF) can be formulated without referencing a cosmological model or input power spectrum. HMF obtained from N-body simulations show systematic deviations of 5-20\% from theoretical predictions. The physical origin of deviations may result from cosmology, the power spectrum, or both. We examine HMF deviations from universality for scale-free power spectra with an Einstein-de Sitter cosmology. We demonstrate that the mass function exhibits an explicit dependence on the slope of the input power spectrum. We find that an effective index of the $\Lambda$CDM model can correspond to the HMF from scale-free cosmologies as a first approximation. Furthermore, structure formation has led to deviations from homogeneity and isotropy on scales up to at least $100$ Mpc/h, expected to affect measurements of $H_0$. We revisit this issue of the concordance model. We find a correlation between errors in $H_0$ estimates and the density around the observer. Further, our mock observations reveal that deviations of up to 5\% can occur in Milky Way-sized halos. While this finding alone does not fully resolve the Hubble tension, it may account for part of it. It is essential to understand the limitations of N-body simulations to avoid misinterpreting data. We show that the missing power at small scales introduces errors in the root-mean-square fluctuations and in the simulated mass function. Our analytical calculation indicates that mode coupling between small and large scales depends on resolving collapsed halos. Therefore, accurate mode coupling estimates require sufficient halos in the simulation.

We apply the slip-back mapping method of Titov et al. 2009 and Lionello et al. 2020 to a thermodynamic MHD simulation to track topological changes in the magnetic field at a range of temporal cadences. The method constitutes the logical successor to a simple open-field map for a steady-state model, as it tracks changes in the open and closed fields for a time-dependent model by tracking individual magnetic elements as they advect across the map, rather than simply tracing field line connectivity from each cell. Through careful categorization of the slip-back mapping values and analysis of the flux changes, we not only effectively track the open flux but can recover the flux processed through interchange reconnection as well. The field lines involved in these processes are shown to follow lines of high squashing factor, as proposed by interchange reconnection-driven slow solar wind theory. The time-dependent model, which is scaled to solar minimum-like activity, projects that a median value of 3.5% of the total open flux in any given 24-hour interval has been processed through interchange reconnection. This corresponds to a relatively high proportion of the total open flux changes over time in the heliosphere. Our results show that not only is this method a useful tool for accurately tracking topological change in time-dependent simulations, but that its inherent complexity can be visually reduced into an intuitive 2D plot that simply and effectively communicates temporal changes.

We report the discovery of DESI-HVS1, a hypervelocity star (HVS) candidate identified from DESI DR1 spectroscopy and Gaia DR3 astrometry. DESI-HVS1 is an old, low-mass, metal-poor F-type star with a mass of $0.8\,M_\odot$, an age of $\sim14.1$~Gyr, and $\mathrm{[Fe/H]}=-1.6$. It is located at a heliocentric distance of $3.77^{+0.39}_{-0.36}$~kpc and has a Galactocentric total velocity of $523^{+46}_{-47}\,\mathrm{km\,s^{-1}}$, marginally exceeding the local escape speed, corresponding to an unbound probability of $P_{\rm ub} \sim 50\%$. Backward orbit integrations show that DESI-HVS1 had a closest approach to the Galactic Centre (GC) of $0.40^{+0.23}_{-0.11}\,\mathrm{kpc}$, with a velocity of $682^{+22}_{-35}\,\mathrm{km\,s^{-1}}$ and a flight time of $12.89^{+0.92}_{-0.74}\,\mathrm{Myr}$. The reconstructed orbit exhibits a clear perigalactic turning point and only a single crossing of the Galactic midplane ($P_{\rm cross} > 0.95$). These properties suggest that DESI-HVS1 is most naturally explained by the Hills mechanism, although alternative scenarios cannot be entirely ruled out. Its discovery provides the first strong evidence for an old, low-mass HVS candidate consistent with a GC origin, indicating that the apparent dominance of young, massive GC-origin HVSs is likely a consequence of observational selection effects.

Using the ssrAp star 33 Lib (HD137949) as an example, we show that magnetic braking caused by the interaction of a strong magnetic field with a stellar wind can play a key role in slowing the rotation of ssrAp stars. Numerical modeling of stellar rotation spin-down in the MESA package, taking into account the evolution of magnetic fields and the stellar wind, shows that stars with rotation periods of up to 80 years and longer can form. Moreover, braking by a magnetized wind makes it possible to estimate the mass-loss rate for stars of moderate mass ($1.25 M_\odot < M < 2 M_\odot$), which is difficult to do by other methods. We introduce the dimensionless parameter $\Xi$, which reflects the spin-down time and the stellar lifetime. Thus, when $\Xi \gg 1$, braking is important, whereas when $\Xi \ll 1$, it is negligible.

The 21-cm signal, one of the most promising probes of the high-redshift Universe, has traditionally been modelled without accounting for the effects of active galactic nuclei (AGN) in the pre-JWST era, primarily due to the lack of observational evidence for AGNs at z > 6. However, following the discovery of several AGNs at redshifts as high as z ~ 10 by JWST, it has become imperative to incorporate the impact of these early AGNs when predicting the 21-cm signal. Supposing that these AGNs are seeded by primordial black holes (PBHs), we study their impact with a semi-numerical model setup. Specifically, we extended the explicitly photon-conserving reionization framework, SCRIPT, including essential cosmic dawn physics and PBH contributions. This enables us to compute both the global signal and the power spectrum of the 21-cm line over the redshift range z ~ 30 - 5 within a self-consistent framework. Building on this setup, we then investigate the impact of different PBH mass functions (obeying current observational constraints) on the resulting signal. The X-ray heating from PBHs can substantially make the depth of the global 21-cm signal shallower and suppress the expected power amplitude during cosmic dawn. We also find that the choice of mass function plays a crucial role in shaping the 21-cm signal, and can, in fact, lead to significantly different predictions.

We introduce a family of phenomenological cosmological models featuring an interacting dark sector modulated by a sparseness scale parameter, in order to describe the late-time accelerated expansion of the universe. The sparseness scale, inspired by well-established saturation mechanisms in ecology and biology, is introduced in the interaction as a half-saturation constant that bounds the energy exchange between dark matter and dark energy, controls the dynamical behaviour of the physical variables and can prevent the phantom crossing. We consider three nonlinear interacting models, where two of them recover the linear interacting scenarios when the sparsity parameter vanishes. We examine the phase-space of the cosmological field equations by using the Hubble normalization approach. We determine the stationary points and their stability properties in order to reconstruct the asymptotics behaviour of the field equations. Such an analysis allows us to demonstrate the effects of the sparseness scale on the background dynamics. We test the interacting models with observational data. Specifically, we employ Supernovae catalogues, cosmic chronometers, Baryon Acoustic Oscillation measurements from DESI DR2, and redshift-space distortion measurements of the growth of large-scale structure through the $f$ and $f\sigma_8$ observables. The Bayesian analysis suggests that, for two of the three models, a vanishing sparsity parameter is disfavoured at more than the 95\% confidence interval, providing observational support for a nonzero sparseness scale in the dark sector interaction.

The Seyfert 1 galaxy NGC 3783 hosts a multiphase warm absorber (WA) that has been extensively studied in the X-ray band. High-resolution spectra from 2000-2001 revealed a complex outflow with multiple ionization and velocity components. Two decades later, new XMM-Newton and XRISM observations allow us to investigate the long-term evolution of these outflows. We perform joint spectral modeling of the XMM-Newton/RGS and XRISM/Resolve time-averaged spectra using the pion photoionization code within SPEX. We derive the ionization parameter, column density, turbulent velocity, and outflow velocity for each absorption component, and investigate their thermal stability and Absorption Measure Distribution (AMD) to characterize the physical and dynamical properties of the WA in NGC 3783 in 2024. We compare these results with the 2000-2001 epoch to assess long-term variability, stability, and possible changes in the absorber population. We identify eight WA components spanning log $\xi =$ 1.08-3.38 and outflow velocities of 480-1230 km s$^{-1}$. The ranges of column densities and turbulent velocities remain broadly consistent with the WAs from 2000-2001, but the earlier data contained more low-ionization, high-velocity components. The total column density in 2024 is 1.5 times larger than in 2000-2001, requiring replenishment by fresh material. The dominant Unresolved Transition Array (UTA) absorber (Comp. B3) has increased its column density by a factor of three while maintaining a similar ionization parameter. The WAs in NGC 3783 have undergone significant structural and dynamical evolution over the past 24 years.

We investigated the high-mass star formation activity in a subregion of the Sagittarius E star-forming complex, centered at (l,b) = (358.69 deg, 0.03 deg), where infrared and radio sources trace a prominent U-shaped structure that has not been identified in previous studies. We used radio continuum data from the Global View on Star Formation (GLOSTAR) survey, which is a wide-band radio (4-8 GHz) survey of the Milky Way that combines data from the Karl G. Jansky Very Large Array and the Effelsberg 100 m telescope. Using BLOBCAT source extraction software, we identified 49 compact radio sources. Based on multiwavelength associations and spectral index estimates, we identified GLOSTAR counterparts to 27 previously confirmed HII regions, detected radio emission from 3 WISE "radio-quiet" candidates, and report 5 new HII region candidates. The derived physical properties indicate that most are relatively evolved HII regions. We find around 50 cold dust clumps, predominantly toward the south and southeast. Mid-infrared flux-ratio maps ([4.5]/[3.6]) show localized shock enhancements along the arc and adjacent clumps, and 15 clumps exhibit SiO emission with broad components indicative of shocks. Together with CO data, the SiO velocity components delineate a continuous (>100 km/s) velocity bridge that links the far dust-lane inflow to the central molecular zone (CMZ) stream. The largest concentration of clumps and compact HII regions lies at this interface. These combined diagnostics favor a scenario in which bar-driven cloud-cloud collision at the far dust-lane-CMZ interface compressed the gas and triggered the observed high-mass star formation.

High radio frequencies observations with the Italian network of large single-dish radio telescopes resulted in ~450 solar images between 2018 and 2023 in K-band frequency range (18-26 GHz). Solar radio mapping at these frequencies allows the probing of the Active Regions (ARs) chromospheric magnetic field close to the Transition Region, where strong flares and coronal mass ejection events occur. Enhanced magnetic fields up to 1500-2000 G determine anomalous spectra in the ARs brightness compared to pure free-free emission, due to the addition of a steeper gyro-resonance component also associated with circular polarisation up to ~40%. When a significant AR spectral flattening is detected, the probability of a strong flare occurrence within ~30 hours is high (~89% in terms of statistical precision). Despite an approximate weekly cadence of our observations, only ~12% of strong flares are missed/unpredicted within this time interval. Through a correlation analysis, we assess the trade-off on the sensitivity and the robustness of this physics-based flare forecast method.

Cosmic rays (CRs) drive ionization and influence gas dynamics in molecular clouds (MCs), potentially impacting the resulting star formation outcomes. Although previous simulations of individual star formation have included methods for cosmic ray transport (CRT), none have been large enough to resolve the stellar initial mass function (IMF). We conduct numerical simulations following the collapse of a $20000 M_{\odot}$ MC and the subsequent star formation including CRT, both with and without CRs accelerated by winds from the young massive stars, and compare against a non-CRT simulation. We show that after the first massive stars form, the cavity produced by feedback is more pronounced in the CRT simulations because the external CRs are able to propagate inwards and compress the gas into higher density structures. This increases the subsequent star formation in the cloud; by the end of the simulation, the SFE in the CRT simulation including stellar wind CRs is 43 \% higher than the non-CRT simulation. The IMF is also top heavy in comparison, with a slope above 1 $M_{\odot}$ that is shallower by $\sim 20$ \%. These effects are also present in the simulation without wind-accelerated CRs, but they are not as pronounced; the SFE is only 16 \% higher than the non-CRT simulation, and the IMF high-mass slope is shallower by $\sim 10$ \%. These results may explain some of the observed top-heavy IMFs, which typically occur in high-CR environments such as the galactic center.

We model the multi-band lightcurves of 80 SNe Ic-BL, including 11 associated with lGRBs, using a magnetar engine model with $^{56}$Ni decay. We find that the data are all consistent with a magnetar central engine, and such a model yields high-quality fits across the sample. The medians with $1\sigma$ regions of the key parameters are $P_{\rm{i}}\sim2.04^{+1.84}_{-0.96}\,{\rm{ms}}$, $B_{\rm{p}}\sim3.96^{+3.28}_{-1.40}\times10^{15}\,{\rm{G}}$, $M_{\rm{ej}}\sim2.30^{+1.48}_{-1.02}\,M_\odot$, and $M_{\rm{Ni}}\sim0.18^{+0.14}_{-0.09}\,M_\odot$, with strong and statistically significant correlations observed for both $M_{\rm{ej}}-P_{\rm{i}}$ (anti-correlation) and $M_{\rm{Ni}}-M_{\rm{ej}}$ (correlation). Comparing the SN Ic-BL samples with and without lGRB association using fitting parameters, we find no significant difference between them, although the GRB-associated sample is slightly brighter, possibly due to an observational bias. Relative to ordinary SNe Ic, SNe Ic-BL have similar $^{56}$Ni and ejecta masses, suggesting comparable pre-SN progenitor properties, with differences possibly arising from the presence of a magnetar engine. In comparison with other possible magnetar-powered SESNe, including SLSNe Ic and FBOTs, we confirm a strong universal $M_{\rm{ej}}-P_{\rm{i}}$ correlation, indicating a common origin. SNe Ic-BL and SLSNe Ic have similar ejecta mass distributions, typically $M_{\rm ej}\gtrsim0.5\,M_\odot$, while FBOTs mostly lie below this value. Differences between SNe Ic-BL and SLSNe Ic may arise from magnetar properties, with SN Ic-BL magnetars rotating faster and having stronger fields. Moreover, the $P_{\rm{i}}-B_{\rm{p}}$ distribution of lGRB magnetars largely overlaps with that of SN Ic-BL magnetars. In connection with binary simulation results, we propose a unified physical classification and progenitor framework for magnetar-powered and ordinary SESNe.

(Abridged) We aim to investigate the inner regions of large and massive disks orbiting T Tauri stars, thought to be progenitors of systems with wide-orbit planets and possible cases of halted pebble drift. We analyze the MIRI spectra of three disks from the MINDS program: V1094 Sco, DL Tau, and IM Lup. The spectra reveal a striking diversity. V1094 Sco and DL Tau exhibit the highest C$_2$H$_2$/H$_2$O flux ratio in the MINDS sample of T Tauri disks. In V1094 Sco, even cold C$_4$H$_2$ is seen. In contrast, the IM Lup spectrum is dominated by O-bearing species. No one-to-one correspondence is found between the gas in the outer disk, as traced by the C$_2$H/C$^{18}$O flux ratio, and that of the inner disk as traced by the C$_2$H$_2$/H$_2$O flux ratio. To explain these results, we propose a scenario based on a toy model of halted pebble drift. We show that a volatile C/O ratio close to unity and low C and O abundances in inner disks arise only if: (1) ~95$\%$ of the icy grains are blocked in the outer disk, (2) the outer disk is chemically evolved, and (3) the gas in the outer disk has had time to reach the inner disk. DL Tau and perhaps V1094 Sco would be the rare examples for which all these conditions are met. Therefore, a high C$_2$H$_2$/H$_2$O flux ratio in pebble-rich disks would have a different origin than proposed for very-low mass stars, for which fast drift of O-rich pebbles would eventually leave a C-rich inner disk. We also show for the first time that the disks with high C$_2$H$_2$/H$_2$O flux ratio exhibit a prominent silica dust component, a result found in four disks published so far (V1094 Sco, DL Tau, CY Tau, DoAr 33). We propose that the reformation of dust at the sublimation front of silicates in a gas with super-solar (but below unity) C/O ratio leads to a silica stoichiometry (SiO$_2$). In turn, silica is a promising diagnostic of the C/O ratio in the inner disks.

We have searched for the 3.55 keV line from sterile neutrino decay using 3.1 megaseconds of RXTE cluster data. A 2.5$\sigma$ excess of emission over a thermal model is found over the energy span of the 3.55 keV line in the combined spectra of the eight clusters that individually have an excess. The residuals are added to increase the signal to noise ratio of the excess, which is then modeled with a Gaussian to simulate the instrumental spectral response. We find a significant correlation (r = 0.76) for a line centered at 3.6 keV with a model flux of 3.07 x 10$^{-5}$ ph cm$^{-2}$ s$^{-1}$. Mixing angle for detected clusters ranges from 0.35 to 6.2 x 10$^{-10}$. The decay rate inferred from the line flux is strongly correlated (r = 0.87) with cluster temperature, which is due to hotter, more massive clusters having a larger amount of dark matter. Approximately half of the decay line total flux comes from the Coma cluster. We fit the Coma cluster spectrum with two different three-component models. The first includes a Gaussian fixed at 3.55 keV to model soft emission. The second three-component model uses a second thermal component to model soft emission. The model fit parameters indicate that the second thermal component is modeling high-energy residuals rather than low ones, where the Gaussian is. Though our line fluxes exceed most reported detections and upper limits, they do not overproduce the dark matter. We conclude that some fraction of the marginally detected excess could be attributed to the decay line since low-temperature thermal emission and systematics fail to model it completely.

We have succeeded in operating a transition-edge sensor (TES) spectrometer and evaluating its performance at the SPring-8 synchrotron X-ray light source. The TES spectrometer consists of a 240 pixel National Institute of Standards and Technology (NIST) TES system, and 220 pixels are operated simultaneously with an energy resolution of $4$~eV at 6~keV at a rate of about 1~c/s/pixel. The tolerance for high count rates is evaluated in terms of energy resolution and live time fraction, leading to an empirical compromise of about 2 x 10^3 c/s/all pixels with an energy resolution of 5 eV at 6 keV. By utilizing the TES's wide-band spectroscopic capability, simultaneous multi-element analysis is demonstrated for a standard sample. We conducted X-ray absorption near-edge structure (XANES) analysis in fluorescence mode using the TES spectrometer. The excellent energy resolution of the TES enabled us to detect weak fluorescence lines from dilute samples and trace elements that have previously been difficult to resolve due to the nearly overlapping emission lines of other dominant elements. The neighboring lines of As K alpha and Pb L alpha2 of the standard sample were clearly resolved and the XANES of Pb L alpha2 was obtained. Moreover, the X-ray spectrum from the small amount of Fe in aerosols was distinguished from the spectrum of a blank target, which helps us to understand the targets and the environment. These results are the first important step for the application of high resolution TES-based spectroscopy at hard X-ray synchrotron facilities.

The James Webb Space Telescope and the future Habitable Worlds Observatory aim to discover exoplanet atmospheric spectra that detect life. Currently, most existing spectral "retrieval" algorithms focus on inferring the abundances of biogenic gases from these spectra. However, abundances are hard to interpret as signatures of life because they are modified by photochemistry, climate, and atmospheric escape. To address this problem, we develop a method for inferring the fluxes of gases at a planetary surface by inverting a coupled photochemical-climate model. As a proof-of-concept, we apply the approach to a synthetic 10-transit JWST NIRSpec Prism spectrum of TRAPPIST-1 e assuming it hosts a biosphere similar to the Archean Earth's. The retrieval confidently detects CO$_2$ and CH$_4$ and can constrain the flux of CH$_4$ into the atmosphere to within approximately 1.5 orders of magnitude (68$\%$ credible interval) provided that TRAPPIST-1's near-UV spectrum is accurately known. We demonstrate how inferred surface gas fluxes naturally fold into a probabilistic assessment of life, finding that ~ 80$\%$ of the surface gas flux posterior is consistent with a CH$_4$-producing metabolism for our nominal test case. As with any inverse problem, these results are conditional on a number of assumptions in our forward model. Overall, we argue that increasing the robustness of life detection on exoplanets requires moving beyond atmospheric abundances toward inference of the surface fluxes that sustain them.

Systematic errors in the parameter estimation (PE) of gravitational wave (GW) mergers can arise from various sources, including waveform systematics, noise mischaracterization, data analysis artifacts, and other unknown factors. In this study, we analyze selected events from the first three observing runs of the LIGO-Virgo-KAGRA (LVK) collaboration. We choose events that have been flagged in various studies as potentially affected by systematic errors. Here, we reanalyze these events using a couple of parametric models developed in previous work that incorporate uncertainties in both the phase and amplitude of the GW waveform. In this data-driven approach, we apply sufficiently broad priors on the uncertainty parameters to account for potential systematic errors. Our findings show that the proposed method effectively reduces systematic errors, even those arising from data artifacts, such as glitches occurring near a signal and the deglitching process in GW frame files. Similarly, inconsistent results from different waveform models become much more consistent in our framework. One noteworthy event we examine is GW191109\_010717, which is particularly interesting due to its anti-aligned spin properties. We report that, within our framework, the event still exhibits anti-aligned spin characteristics, but the inference results become consistent across raw and deglitched frame files, as well as across the waveform models used for this event (IMRPhenomXPHM, IMRPhenomXO4a, and NRSur7dq4). A similar trend is observed for the event GW200129\_065458, which previously yielded a high, but inconsistent precession parameter among different waveform models. In contrast, we observe a non-zero and consistent value of $\chi_{p}=0.60^{+0.31}_{-0.33}, 0.58^{+0.30}_{-0.29}$ and $0.56^{+0.31}_{-0.28}$ for the IMRPhenomXPHM, IMRPhenomXO4a, and NRSur7dq4 waveform models, respectively.

We use GRAVITY near-infrared (NIR) flare astrometry to test whether Sagittarius A* could be a solitonic boson star. We consider five spherically symmetric solitonic boson-star models with different effective radii, together with the Schwarzschild black hole. Treating the flares as hot spots on circular equatorial orbits, we analyze their centroid motions and images in these spacetimes and use them for parameter fitting. We perform the fitting using both $\chi^2$ analysis and Markov Chain Monte Carlo (MCMC) methods, which yield consistent results: the inferred masses of boson-star models are systematically larger than the established value of $4.3\times10^6M_\odot$. Notably, more diffusive boson stars exhibit imaging properties closer to those of a black hole, leading to mass estimates that are correspondingly closer to the established value. Overall, our results place stringent constraints on solitonic boson star interpretations of Sagittarius A*, although do not completely rule them out.

The Dyson Minds 2025 Workshop, held at the Center for Brains, Minds & Machines at MIT and organized by Penn State, MIT, and The Ultraintelligence Foundation, brought together researchers in astrophysics, engineering, artificial intelligence, computer science, and philosophy to examine "Dyson Minds" -- large-scale post-biological intelligences powered by energy harvested from supermassive black holes (SMBHs). Building on the ideas of F. J. Dyson (1960, 1966) and I. J. Good (1966), participants explored the physical, engineering, behavioral, and observational consequences of civilizations embodied as machinery operating near the universe's most powerful energy sources. The workshop aimed to develop new observational strategies capable of detecting signatures of such systems. Despite the highly cross-disciplinary scope, discussions centered on how a Dyson Mind might be constructed, how it might behave, and how those factors would shape strategies for the search for extraterrestrial intelligence. Key themes included the thermodynamic, mechanical, and stability limits of Dyson swarms; the trade-offs between power availability and communication latency in distributed minds; and how observability changes depending on whether Dyson Minds act as coherent entities or as loosely coordinated collectives. Across these topics, the consensus was that details of architecture and behavior strongly influence observational signatures. A major recommendation was to apply anomaly-detection methods to archival datasets, including those from WISE, JWST, and the Event Horizon Telescope, to identify unusual sources potentially overlooked by standard reduction pipelines. By integrating insights from multiple disciplines, the meeting advanced concrete, observation-focused strategies for future technosignature searches around SMBHs.

Hydrogen atom tunneling likely plays a substantial role in the gas-phase chemistry of astrochemical environments. To determine the potential effect that it has on the chemical modeling of aromatic molecules, we screened the kida.uva.2024 network, and our own expanded network to find reactions which could be significantly accelerated by hydrogen atom tunneling in the ISM. In total, 64 reactions were identified. The hydrogen abstraction reactions from H$_{2}$ to four key interstellar radicals (C$_{2}$H, OH, CN, and NH$_{2}$) were studied further using newly calculated potential energy surfaces and RRKM analyses to determine rate coefficients for a temperature of 10 K and a density of 2 $\times$ 10$^{4}$ cm$^{-3}$. Despite having low rate coefficients of 1.66 $\times$ 10$^{-15}$, 8.17 $\times$ 10$^{-16}$ and 3.15 $\times$ 10$^{-16}$ $\mathrm{cm^{3}\,s^{-1}}$ the C$_{2}$H, OH, and CN reactions are competitive in the ISM, due to large overall rates caused by the high abundance of molecular hydrogen. The calculated value for the NH$_{2}$ reaction, however, was much smaller and found to be inefficient at ISM conditions. The possible effects of all other considered reactions were studied with simulations using calculated collision limit rate coefficients. Upper and lower bounds were then placed on modeled aromatic abundances using the most significant reactions. Due to the dependence of calculated aromatic abundances on reactions involving c-C$_{6}$H$_{5}^{+}$ and the recent questions surrounding its reactivity, we also explored the abundance variations caused by reactions leading to or involving c-C$_{6}$H$_{5}^{+}$.

Massive stars with their strong ionizing radiation and strong stellar winds are the key feedback agents of the universe. Stellar winds of massive stars are often measured by fitting resonance lines in the UV using non-LTE stellar atmosphere models. So far, the line formation regions of these lines have not been measured empirically, preventing a comparison to the model's structures. We aim to conduct the first measurement of the resonance line formation regions in an early-type eclipsing binary in the SMC, namely AzV 75. We employ TESS and ASAS-SN photometry in combination with radial velocity measurements from multi-epoch HST UV spectra to derive the ephemeris. We examine the intensity changes in the C IV and N V resonance lines in the UV and combine them with a light-curve analysis to estimate the region in the wind where these lines are formed. AzV 75 has an orbital period P=165.66d, eccentricity e=0.42, mass ratio q=0.72, and inclination i=85.77°. With this orbital configuration, no secondary eclipse is expected. We report that the optically thick UV resonance lines exhibit flattening and shortening of the absorption trough, and weakening of their emission features, as they approach the phase of the expected secondary eclipse, while the continuum UV flux appears to remain unaffected. We illustrate that this can be explained by the primary's optically thick wind eclipsing the secondary star. The C IV and N V resonance line formation regions in the primary star extend up to 316 Rsol. The measured extend of the formation regions of resonance lines in a stellar wind are important benchmarks for 1D as well as 3D non-LTE stellar atmosphere models. A first comparison to 1D-stellar atmosphere models indicates that a classical beta-law with an exponent of beta=0.5 instead of beta=0.8 might be favoured for the primary star's velocity field.

The spectral kernel field equation R[k] = T[k] lacks a conservation-law analog. We prove (i) the fixed-point flow is strictly volume-expanding (tr DF > 0), precluding automatic conservation, and (ii) the conservation deficit per mode equals the Hessian stability margin exactly: D_m = -Delta'. Closing the deficit requires a scene-side compensating contribution, which we formalise as the distinction dynamics equation dc/dt = G[c, h_t], with MaxCal-optimal realisation G_opt. On fixed-topology 3D surface graphs we derive a conditional topology-preserving compression theorem: retaining k >= beta_0 + beta_1 modes (under a spectral-ordering assumption) preserves all Betti-number charges; we include a worked short-cycle counterexample (figure-eight) calibrating when the assumption fails. A triple necessary spectral diagnostic -- Fiedler-mode concentration, elevated curl energy, anomalous beta_1 -- is derived for planetary drainage networks at O(N) cost. Two internal real-data sequences serve as preliminary consistency checks; full benchmarks and adaptive-topology extensions are deferred.

Ringdown signals from perturbed black holes (BHs) offer a clean window into BH spacetime, strong-field gravity, and fundamental physics. Presently the quasi-normal modes of stellar-mass BH ringdowns have been successfully extracted in the ground-based gravitational wave (GW) observations. Looking ahead, the future space-borne observatories will listen to the ringdowns from massive BH binary coalescences more loudly and resolve multiple modes to unprecedented precision, which calls for efficient approaches to mitigate the sharply increasing computational burden. We develop a practical ringdown analysis pipeline for space-borne detectors by implementing FIREFLY, a novel acceleration algorithm validated in ground-based detectors, and for the first time demonstrate its compatibility and effectiveness with the time-delay interferometry (TDI) observables. With high fidelity, we achieve a $\sim 200$-fold speedup for a simulated ringdown signal including six modes, providing a viable and scalable route for multi-mode ringdown analysis in the space context. This new approach has sound statistical interpretation and is extensible to other GW sources in band.

We present a covariant and gauge-invariant formulation of the theory of radial adiabatic linear perturbations of self-gravitating, non-dissipative imperfect fluids within the theory of general relativity. By codifying the thermodynamical properties of the source into an equation of state and an ansatz on anisotropic pressure that involves both matter and kinematic variables, we obtain a set of equations that is directly applicable to a wide variety of thermodynamic theories for matter fields. As examples, we evaluate and compare the predictions of the Eckart theory, the Bemfica-Disconzi-Noronha-Kovtun theory, and the Truncated Israel-Stewart theory on the properties and evolution of radial adiabatic perturbations of stellar compact objects modeled by classical equilibrium solutions. Introducing a new solution of the Einstein field equations, and imposing causality, we propose an upper bound for the maximum compactness of dynamically stable stars with non-trivial radial and tangential pressures.

We use the 3D fully kinetic simulation to study different turbulence modes and turbulence anisotropy of relativistic turbulence in magnetically dominated collisionless plasmas. We extend the method developed by Cho & Lazarian (2002) for decomposing non-relativistic magnetohydrodynamic (MHD) turbulence into Alfvén, fast, and slow modes to the regime of collisionless plasmas. We find that Alfvén and slow modes are anisotropic, following the Goldreich & Sridhar (1995) scaling, while fast modes are isotropic. We observe a larger kinetic energy fraction of fast modes compared to that in the non-relativistic MHD turbulence, suggesting a stronger coupling of Alfvén and fast modes in relativistic magnetized turbulence in collisionless plasmas. We further examine the dynamic alignment and find a weaker scale dependence of the alignment angle than previously proposed. The dominant thermal fluctuations in the kinetic range can cause flattening of the turbulent velocity structure function and weakening of the turbulence anisotropy and dynamic alignment near the kinetic scales.

Accurate modeling of ion-molecule reaction networks is essential for understanding the chemical evolution of planetary ionospheres, particularly for giant planets where proton-transfer chains drive atmospheric composition. However, predicting reaction rates in these ultracold environments remains a challenge due to the non-trivial interplay between vibrational dynamics and quantum tunneling. In this work we present a chaos-diagnostic framework that integrates multireference electronic structure theory, Adiabatic Gauge Potentials (AGP), and Random Matrix Theory (RMT) to characterize the microscopic dynamics of proton transport. Using the formation of H+3 and the proton-bound cluster H+5 as representative model systems relevant to Jovian atmospheres, we demonstrate that the transition state acts as a dynamical bottleneck where quantum chaos is notably suppressed, effectively enhancing tunneling probabilities. We introduce a fragility index based on the AGP slope to quantify how specific vibrational modes reintroduce chaos and suppress reactivity. This diagnostic approach offers a generalizable, data-driven metric for identifying vibrationally gated pathways in complex astrochemical networks, providing a theoretical basis for refining kinetic models of planetary and interstellar plasmas

Astrophysics plays a pivotal role in the quest for axions and axion-like particles, offering guidance to experimental efforts and enabling the investigation of axion properties that cannot be probed otherwise. In this context, the extreme conditions in the interiors of compact stellar objects -- such as core-collapse supernovae, neutron stars, and binary neutron star mergers -- significantly enhance axion production, providing unparalleled sensitivity to extremely feeble couplings to Standard Model particles. In this context, the techniques of multimessenger astrophysics deepens the understanding of powerful transient events, maximizing the capabilities of current instruments to identify possible signatures of axion emission.

Ultra high energy cosmic rays have been see coming from the direction of the local cosmic void. We use this fact to argue that at least some of these these cosmic rays are relatively light magnetic monopoles and that their relative fraction above 1020 eV can be found from full sky observations.

We present a dynamical solution to the dark matter-baryon coincidence problem based on the neutron portal operator connecting the visible and dark sector asymmetries. This framework is motivated by the possibility that a strongly supercooled dark confinement phase transition accounts for the nano-Hz stochastic gravitational wave signal observed by pulsar timing arrays, while also generating the dark matter and baryon asymmetry in the Universe. We show that the GeV-scale mass of asymmetric dark matter can be naturally correlated with the (multi-)TeV scale cut-off for the neutron portal through its ultraviolet completion. The dark sector is governed by an approximate fixed point and confines once the heavy portal states are integrated out, dynamically generating a scale of $\mathcal{O} ({\rm GeV})$. We analyze both tree and loop-level ultraviolet completions and demonstrate how the resulting confinement scale is linked to the effective neutron portal scale. We also discuss cosmological constraints and experimental prospects in beam dump searches and colliders for probing the neutron portal.

Superradiant instability of ultralight bosons can produce clouds around rotating black holes, whose continuous gravitational wave (GW) emission is a promising observational target. Precise predictions of the signal frequency and its evolution are essential for detecting such continuous GWs. For axions, self-interactions can populate multiple superradiant modes via nonlinear couplings, and GW emission can occur through various channels. To calculate the frequency shifts of GWs emitted through these channels, we employ relativistic perturbation theory based on a bilinear form. We apply this framework to self-interaction effects for the first time, and also revisit the treatment of the self-gravity contribution. Our results provide a simple and unified framework for calculating frequency shifts, including cases in which multiple modes are excited, and are relevant for next-generation GW observations.

There has been a recent revival in understanding the spontaneous scalarization phenomenon in scalar-tensor gravity as a phase transition. Using the tools of the Landau theory, we now know that first-order transitions where scalarization occurs in a discontinuous manner is more prominent than what had been considered in the literature, and this might lead to novel observation channels. However, the examples so far have been restricted to specific quadratic scalar coupling terms and spherically symmetric stars. Here we explore the phase transition structure of scalarization for more general couplings, considering linear as well as quadratic terms in the conformal scaling factor of the theory. Moreover, we also investigate the effect of rotation on the scalarization phase transition. Both of these considerations are natural choices since the coupling in a scalar-tensor theory can appear at all orders, and astrophysical neutron stars commonly have angular momentum. The introduction of linear coupling leads to a complex solution space which is harder to explore. However, we demonstrate that the Landau model of scalarization enables us to systematically find the branches of scalarized solutions that are commonly overlooked in numerical searches, providing a novel tool. On the other hand, the main effect of stellar rotation is shifting the stellar masses at which the phase transition occurs to higher values, but the qualitative picture remains similar to what happens under spherical symmetry.

Since the dawn of the space age, observations of energetic particles in planetary radiation belts have been interpreted within a diffusive transport framework, even though the processes that populate and deplete these belts produce highly structured and spatially localized distributions. This exposes a fundamental problem: how can coherent phase-space structures evolving under collisionless dynamics give rise to observational signatures that appear consistent with diffusion-based transport? Here we show that diffusion-like behaviour can arise from an observational phase-mixing effect, independent of stochastic wave-particle transport. As spacecraft sample neighbouring drift shells while particles undergo electromagnetic drifts, spatially localized drift-phase structures are converted into rapidly decorrelating temporal signals, making them observationally indistinguishable from stochastic processes. We show that the effective lifetime of these structures is only a few drift periods, preventing the resolution of fine-scale structure. These results demonstrate that collisionless dynamics can mimic diffusive transport on short timescales, limiting the inference of particle acceleration processes and biasing transport estimates. This calls for a reassessment of diffusion-based interpretations of radiation belts at Earth, across the solar system, and in the radiation belts of ultra-cool brown dwarfs.

The recent detection of a stochastic gravitational wave background by pulsar timing arrays has opened a new window in understanding supermassive black hole binaries and in probing the universe at the early time. Recently, pulsar timing array (PTA) collaborations have been further paving the way to probe anisotropies in the stochastic gravitational wave background. This study investigates dipole-type statistical anisotropy in the primordial power spectrum within a phenomenological framework. We demonstrate that the primordial dipole induces both dipolar and quadrupolar anisotropies in the energy density spectrum of scalar-induced gravitational waves (SIGWs), without generating extra polarization modes. Based on this anisotropic spectrum, we derive the corresponding PTA overlap reduction functions (ORFs), which exhibit frequency dependence, with the anisotropies enhanced on small scales. Furthermore, owing to the non-uniform distribution of millisecond pulsars over the sky in current PTA dataset, the ORFs exhibit a morphology that explicitly depends on the preferred direction of the anisotropy. However, our bayesian analysis of the NANOGrav 15-year dataset still yields no significant evidence for a preferred direction and a weak upper limit on anisotropy amplitude $(g\lesssim0.5)$. This result arises because the observational frequency band lies below the spectral peak, where our models predict suppressed anisotropic contributions. This limitation highlights the potential of future PTA observations. Specifically, datasets with broader frequency coverage are expected to tighten constraints on dipole-type anisotropy.

Dark matter search strategies have started advancing towards the neutrino fog. In this regard, compact objects such as neutron stars have already demonstrated their ability in probing such low DM-nucleon cross-sections from dark matter induced effects. In the optically thin limit, effect of dark matter self-interaction becomes relevant and may assist the capture and thermalization of dark matter inside stars, imparting observable changes on neutron star temperatures. The resulting radiation although weak can be potentially detected by the James Webb Space Telescope and upcoming Thirty Meter Telescope and the European Extremely Large Telescope. Observation of cold neutron stars accompanied by advancements in direct detection probes would provide stringent constraints or a smoking-gun signature for dark matter self-interactions. The potential detection of a neutron star with surface temperatures $\sim (1000 - 1200)$ K in the optically thin limit can push the bounds on asymmetric dark matter self-interaction cross-section to approximately two orders of magnitude more stringent than the bullet cluster.

In this work, we systematically investigate the quantum-information diagnostics of cosmological perturbations with a nontrivial sound speed, utilizing a normalized open two-mode squeezed-state framework. Rather than introducing new observables, our analysis focuses on how a modified sound speed dynamically reshapes the Schrödinger evolution of the squeezing parameters ($r_k$ and $\phi_k$). We demonstrate how these dynamical changes are inherited by the reduced density matrix of the observable sector. By employing a sound-speed-resonance parametrization, we derive and evaluate the purity, von Neumann entropy, Rényi entropies, and logarithmic negativity. To overcome the intrinsic multiscale stiffness of the post-inflationary equations, we introduce a bounded variable $x = \tanh r_k$ as a partial regularization, which enables reliable numerical simulations exclusively within the inflationary regime. Our numerical results reveal that a nontrivial sound speed significantly suppresses the purity of the reduced state, indicating enhanced effective mixedness. Simultaneously, it strongly amplifies and modulates both the entropic and entanglement diagnostics. More precisely, a nontrivial sound speed postpones the onset of classicality by modulating the decoherence process. Ultimately, we show that a nontrivial sound speed leaves distinct and identifiable quantum-information signatures within the entanglement structure of the early universe.

This article opens new window to obtain accelerating scaling attractors without any need of dark energy. We study cosmological dynamics in a two-fluid system where pressureless dark matter (DM) undergoes adiabatic particle creation and exchanges energy with a barotropic fluid. Considering six widely used interaction prescriptions, we formulate the corresponding autonomous systems in a compact phase space and perform a unified dynamical analysis. We find that accelerating scaling attractors, namely late-time states where both fluids coexist with fixed energy fractions, arise only when the interaction is controlled by the DM density and energy flows from DM to the second fluid. Such attractors appear in the global and local DM-based interactions, and in the global mixed case, but are entirely absent when the interaction depends on the second fluid or on local mixed terms, which instead drive the universe to a DM-dominated accelerating phase. These results clarify the unique conditions under which matter creation can mimic dark-energy-like behaviour without introducing a dark-energy component.