Papers for Thursday, Oct 30 2025

We present the first stellar stream discovered with the Vera C. Rubin Observatory, around spiral galaxy M61 (NGC 4303) in Virgo First Look imaging. The stream is narrow, radially-oriented in projection, and ~50 kpc long. It has g-band surface brightness (SB) mu_g ~ 28 AB mag arcsec^-2, color g-z ~ 1.0, and stellar mass M_* ~ 2x10^8 M_Sun. This dwarf galaxy interaction may have provoked the M61 starburst, and foreshadows the bounty of accretion features expected in the ten-year Rubin Legacy Survey of Space and Time (LSST).

We combine deep photometric data in the COSMOS and XMM-LSS fields with high-resolution cosmological hydrodynamical simulations to explore two key questions: (1) how does the galaxy stellar mass function, particularly in the dwarf (Mstar < 10^9.5 MSun ) regime, vary with environment, defined as distance from the large-scale structure (LSS) traced by nodes and filaments in the cosmic web? (2) is there a generic 'missing dwarfs' problem in LambdaCDM predictions when all environments - and not just satellites around Milky Way like galaxies - are considered? The depth of the observational data used here enables us to construct complete, unbiased samples of galaxies, down to Mstar ~ 10^7 MSun and out to z ~ 0.4. Strong environmental differences are found for the galaxy stellar mass function when considering distance from LSS. As we move closer to LSS, the dwarf mass function becomes progressively flatter and the knee of the mass function shifts to larger stellar masses, both of which result in a higher ratio of massive to dwarf galaxies. While the stellar mass functions from the three simulations (NewHorizon, TNG50 and FIREbox) considered here do not completely agree across the dwarf regime, there is no evidence of a generic missing dwarfs problem in the context of LambdaCDM, akin to the results of recent work that demonstrates that there is no missing satellites problem around Galactic analogues.

We derive a simple, accurate, non-linear, global equation governing spiral density waves in thin, non-self-gravitating, inviscid accretion discs. These discs may have any slowly varying surface density or temperature profile. For specific 'self-similar' disc profiles, solutions to our equation match (novel) smooth non-linear exact spiral solutions derived via a separate method, which highlight that non-linear spiral waves need not shock. Indeed, at low amplitudes, we find that dispersion can overcome wave steepening, and may prevent the inner spiral wakes excited by low mass planets (below roughly 1% of a thermal mass) embedded in protoplanetary discs from shocking. At high amplitudes, we find a simple universal description of non-linear spiral waves with shocks, as well as caps on the possible amplitude and wave action flux of non-linear spirals both with and without shocks, depending on how many arms they have. We further find that highly non-linear spirals are far more loosely wound than their linear counterparts. These developments shed light on why two-armed spirals are prevalent across a range of astrophysical systems which don't necessarily possess an intrinsic twofold symmetry, and why they appear surprisingly loosely wound in observations. These results are supported by very high-resolution numerical simulations.

This study of the M81 complex and its Hubble flow delivers new and improved Tip of the Red Giant Branch (TRGB)-based distances for nine member galaxies, yielding a total of 58 galaxies with high-precision TRGB distances. With those, we perform a systematic analysis of the group's dynamics in the core and its embedding in the local cosmic environment. Our analysis confirms that the satellite galaxies of the M81 complex exhibit a flattened, planar distribution almost perpendicular to the supergalactic pole and thus aligned with a larger-scale filamentary structure in the Local Universe. We demonstrate that the properties of the group's barycentre are robustly constrained by the two brightest members, M81 and M82, and that correcting heliocentric velocities for the solar motion in the Local Group decreases the velocity dispersion of the group. Then applying minor and major infall models, we fit the local Hubble flow to constrain the Hubble Constant and the total mass of the M81 complex. The joint best-fit parameters from both models yield $H_0 = \left(63 \pm 6 \right)$ km/s/Mpc and total mass of $(2.28\pm 0.49) \times 10^{12} M_{\odot}$. We thus arrive at an increased mass estimate compared to prior work but reach a higher consistency with virial, $(2.74 \pm 0.36)\times 10^{12}\,M_\odot$, and projected-mass estimates, $(3.11 \pm 0.69)\times 10^{12} M_\odot$. Moreover, our $H_0$ estimate shows an agreement with Planck, consistent with other TRGB-based Local-Universe inferences of $H_0$ and still within a 2-$\sigma$ agreement with Cepheid-based Local-Universe probes.

We investigate accretion disk winds commonly observed in galactic black hole (BH) X-ray binaries (XRB), which manifest as blueshifted absorption features in X-ray spectra. We model these winds as ideal magnetohydrodynamic outflows of hot plasma driven by global magnetic fields threading the accretion disk around the BH. Using Monte Carlo simulations with MONACO, we solve three-dimensional radiative transfer equations to determine the large-scale ionization structure that produces the observed ionic column densities. Focusing on the high/soft state of the BH XRB, where disk emission provides the dominant source of ionizing X-rays, we calculated synthetic spectra showing resonance absorption and scattered emission from ions in various charge states. Our results demonstrate that systems viewed at high polar angles exhibit prominent multi-ion absorption lines with asymmetric profiles, accompanied by P-Cygni-like emission features that partially reproduce the characteristics seen in the observed spectra. This further implies that even a dense disk wind with a high polar angle is unlikely to be saturated due to effective scattering.

We present proper motion (PM) measurements for Draco II, an ultra-faint dwarf satellite of the Milky Way. These PMs are measured using two epochs of Hubble Space Telescope Advanced Camera for Surveys (HST/ACS) imaging separated by a 7 year time baseline. Measuring PMs of low-luminosity systems is difficult due to the low number of member stars, requiring a precise inertial reference frame. We construct reference frames using three different sets of external sources: 1) stars with Gaia DR3 data, 2) stationary background galaxies, and 3) a combination of the two. We show that all three reference frames give consistent PM results. We find that for this sparse, low-luminosity regime including background galaxies into the reference frame improves our measurement by up to $\sim2\times$ versus using only Gaia astrometric data. Using 301 background galaxies as a reference frame, we find that Draco II's systemic PM is $(\mu_{\alpha}^*, \mu_{\delta}) = (1.043\pm0.029, 0.879\pm0.028)$ mas/yr, which is the most precise measurement of the three we present in this paper.

The amount of surface water is thought to be critical for a planet's climate stability and thus habitability. However, the probability that a rocky planet may exhibit surface water at any point its evolution is dependent on multiple factors, such as the initial water mass, geochemical evolution, and interior composition. To date, studies have examined the influence of interior composition on the water inventory of the planet or how surface oceans may be impacted by planet topography individually. Here, we provide the first exploration on the impact of interior composition, topography, and planet radius on the water inventory of rocky planets using a sample of 689 rocky planets with spectroscopically derived stellar abundances from APOGEE and GALAH. We find that the oxidation state of the mantle (FeO content) significantly impacts the mantle water storage capacity and potential for surface flooding. For an FeO ~11 wt%, the water storage capacity of a 1 M$_\oplus$ is 2 times that of Earth, indicating that the oxidation state may reduce the amount of surface water. We quantify the impact of topography on seafloor pressures, showing that flat topographies are more likely to be flooded for all planet compositions and radii. We also find that Mars-like topographies are more likely to have seafloor pressures that may form high-pressure ice, reducing seafloor weathering. Thus, for the first time, we show that the composition and topography of the mantle influence the water inventory of rocky planets.

Pulsar Timing Array experiments are rapidly approaching the era of gravitational wave background anisotropy detection. The timing residuals of each pulsar are an integrated measure of the gravitational-wave power across all angular scales. However, due to the limited number of monitored pulsars, current analyses are only able to reconstruct the angular structure of the background at large scales. We show analytically that this mismatch between the integrated all-sky signal and the truncated reconstruction introduces a previously unaccounted source of systematic bias in anisotropic background map reconstruction. The source of this systematic error, that we call ''small-scale leakage'', is the intrinsic presence of unaccounted gravitational wave power at scales smaller than the reconstructed scales. This unmodeled power leaks into large-scale modes, artificially increasing the recovered value of the inferred angular power spectrum by at least one order of magnitude in a wide range of scales. Importantly, this effect is fundamentally independent of the geometry of the pulsar configuration, the anisotropy reconstruction method, the use of different regularization schemes, and the presence of pulsar noise. As the quality of pulsar timing array experiments improves, a robust understanding of small-scale leakage will become paramount for reliable detection and characterization of the gravitational wave background. Thus, the theoretical formalism developed here will be essential to estimate the magnitude of this systematic uncertainty in anisotropy searches.

The James Webb Space Telescope (JWST) North Ecliptic Pole (NEP) Time-Domain Field (TDF) has been monitored by NuSTAR and XMM-Newton with a regular cadence for five years starting in 2019. The survey has accumulated 3.5Ms of NuSTAR exposure and 228 ks quasi-simultaneous XMM-Newton observations covering 0.31 deg^2. This paper presents the results from the most recent two-years' 2Ms NuSTAR and 166 ks XMM observations in NuSTAR cycles 8 and 9. These observations reached a 20%-area flux of 2.20 x 10^-14 erg cm^-2 s^-1 in the 8-24 keV band. 75 NuSTAR sources and 274 XMM-Newton sources are detected at 99% reliability level. The logN-logS measured in cycles 8+9 are consistent with those measured in the previous cycle 5+6 NuSTAR NEP survey, but in a larger area (0.3 deg^2 compared with 0.19 deg^2). The slope of the cycles 8+9 8-24 keV logN-logS curve is flatter than other works ({\alpha}89 = 1.13 +/- 0.46), but is consistent with the Euclidean value of {\alpha} = 1.50. In addition, we found ~36% of the NuSTAR sources to be heavily obscured (NH >= 10^23 cm^-2). The Compton-thick (NH >= 10^24 cm^-2) (CT-) AGN fraction is 9+18-8% in the NEP-TDF, which is consistent with the measurements in previous surveys.

We present a time-resolved, joint Swift-Fermi spectral study of GRB 241030A (z=1.411) that cleanly isolates the synchrotron origin of its prompt emission and favors a matter-dominated, internal-shock scenario. The light curve shows two episodes separated by a quiescent gap. Episode I (0-45 s) is well described by a single power law with photon index $\simeq -3/2$, consistent with the fast-cooling synchrotron slope below the peak. Episode II (100-200 s), exhibits two robust spectral breaks: a low-energy break at $E_{b}$$\sim$$2-3$ keV that remains nearly constant in time, and a spectral peak $E_{p}$ that tracks the flux within pulses but steps down between them. The photon indices below and above $E_{b}$ cluster around -2/3 and -3/2, respectively, as expected for fast-cooling synchrotron emission. The burst displays an unusually small (consistent with zero) spectral lag across GBM bands. At later times ($\geq 230$ s), the spectrum softens toward $\sim-2.7$, as expected when the observing band lies above both $\nu_m$ and $\nu_c$. These behaviors are difficult to reconcile with a globally magnetized outflow with a decaying field, which naturally produces hard-to-soft Ep evolution, growing $\nu_c$, and appreciable lags. By contrast, internal shocks with a roughly steady effective magnetic field and a time-variable minimum electron Lorentz factor (equivalently, e.g., a varying fraction of accelerated electrons simultaneously account for (i) the stable $E_{b}$, (ii) the intensity-tracking yet step-down $E_{p}$, (iii) the canonical -2/3 and -3/2 slopes, and (iv) the near-zero lag.

The initial-final mass relation (IFMR) links a star's birth mass to the mass of its white dwarf (WD) remnant, providing key constraints on stellar evolution. Open clusters offer the most straightforward way to empirically determine the IFMR, as their well-defined ages allow for direct progenitor lifetime estimates. We construct the most comprehensive open cluster WD IFMR to date by combining new spectroscopy of 22 WDs with an extensive literature review of WDs with strong cluster associations. To minimize systematics, we restrict our analysis to spectroscopically confirmed hydrogen-atmosphere (DA) WDs consistent with single-stellar origins. We separately analyze a subset with reliable Gaia-based astrometric membership assessments, as well as a full sample that adds WDs with strong cluster associations whose membership cannot be reliably assessed with Gaia. The Gaia-based sample includes 69 spectroscopically confirmed DA WDs, more than doubling the sample size of previous Gaia-based open cluster IFMRs. The full sample, which includes 53 additional literature WDs, increases the total number of cluster WDs by over $50\%$ relative to earlier works. We provide functional forms for both the Gaia-based and full-sample IFMRs. The Gaia-based result useful for $M_i \geq 2.67\,\mathrm{M}_\odot$ is $$M_f = \left[0.179 - 0.100 H(M_i-3.84\,\mathrm{M}_\odot) \right ] \times (M_i-3.84\,\mathrm{M}_\odot)+0.628\,\mathrm{M}_\odot$$ where $H(x)$ is the Heaviside step function. Comparing our IFMR to recent literature, we identify significant deviations from best-fit IFMRs derived from both Gaia-based volume limited samples of field WDs and double WD binaries, with the largest discrepancy occurring for initial masses of about $5\,\mathrm{M}_\odot$.

We analyze the Dark Energy Survey (DES) Year 3 data using predictions from the Effective Field Theory of Large-Scale Structure (EFTofLSS). Specifically, we fit three two-point observables (3$\times$2pt), galaxy clustering, galaxy-galaxy lensing, and cosmic shear, using the one-loop expressions for the projected angular correlation functions. We validate our pipeline against numerical simulations and we check for several internal consistencies before applying it to the observational data. Fixing the spectral tilt and the baryons abundance, we measure $S_8=0.833\pm 0.032$, $\Omega_m = 0.272\pm 0.022$, and $h = 0.773\pm 0.049$, to about $3.8\%$, $8.1\%$, and $6.3\%$, at $68\%$CL, respectively. While our results are consistent within $\lesssim 2\sigma$ with those from Planck and the BOSS full-shape analyses, they exhibit stronger tension with those from DES collaboration 3$\times$2pt analysis combined with a Big-Bang Nucleosynthesis prior, highlighting the impact of modeling, scale cuts, and choice of prior. The theory code and likelihood used for our analyses, \texttt{PyFowl}, is made publicly available.

The diminishing of star formation is accompanied by size differentiating, as quiescent galaxies are more compact than star-forming galaxies at fixed stellar mass. In order to understand how galaxy quenching is related to galaxy sizes, we performed a demographic study of 46 massive quiescent central galaxies with stellar mass from $10^{10.5}\rm M_\odot$ to $10^{11}\rm M_\odot$ in the TNG50 simulation. We found that, in addition to the triggering active galactic nucleus (AGN) feedback, galaxy size is also a major determinant of the quenching process, as small and compact galaxies are immediately quenched by the kinetic AGN feedback, while galaxies with large sizes are still active until strangulated by the cutoff of new gas replenishment. Further spatially resolved inspection reveals that this short and intense kinetic AGN feedback can only suppress the star formation within 1-2 kpc, resulting in this size-dependent effect of quenching. We also identify a long-term effect of a few Gyr timescale that the gas inflow rate is progressively suppressed after triggering kinetic feedback, which appears to effectively quench large galaxies entirely. We conclude that kinetic AGN feedback has two key roles in quenching: a short-term, intense effect that quenches the central 2 kpc region, and a long-term effect that suppresses the gas inflow rate and further quenches the entire galaxy.

Thanks to new advances in astronomical imaging, we can now routinely explore disk galaxy profiles about two magnitudes deeper than the data available 20 years ago. In this regard, it is an opportune time to reevaluate the past classifications of different surface brightness Types. In this paper, we explore the stellar mass profiles of a sample of disk galaxies with similar stellar masses ($\sim$10$^{10}$ M$_{\odot}$) using IAC Stripe82 Legacy Project data. We find that Type I, II-CT, and III-d galaxies exhibit edges at surface brightnesses $\mu_g$(R$_{edge}$)$\sim$26.5 mag/arcsec$^2$ and surface mass densities $\Sigma$(R$_{edge}$)$\sim$0.5-1 M$_{\odot}$/pc$^2$. These surface brightnesses were outside the range of typical SDSS images and, therefore, unstudied. The present data suggest (although the statistics are inconclusive) that the main difference between the previous profile Types is the presence of either a more or a less intense bulge. Half of the Type II-CT galaxies are Sc, whereas half of the Type III-d galaxies are Sb. Finally, Type III-d galaxies have low mass surface density disks ($\Sigma$($\Delta$ $R_{\mathrm{exp}}$)$\sim$0.5 M$_{\odot}$/pc$^2$) and are bluer $(g-r)_0$=0.32$\pm$0.02 than the other Types ($(g-r)_0$=0.42$\pm$0.02 for Type I and II-CT), suggesting that their disks probably formed later in time.

The lunar surface is a compelling location for large, distributed optical facilities, with significant advantages over orbital facilities for high spatial resolution astrophysics. The serious development of mission concepts is timely because of the confluence of multiple compelling factors. Lunar access technology is maturing rapidly, in the form of both US-based crewed and uncrewed landers, as well as international efforts. Associated with this has been a definitive maturation of astronomical optical interferometry technologies at Earth-based facilities over the past three decades, enabling exquisitely sharp views on the universe previously unattainable, though limited at present by the Earth's atmosphere. Importantly, the increasing knowledge and experience base about lunar surface operations indicates it is not just suitable, but highly attractive for lunar telescope arrays.

The study of solar active regions (ARs) is of central importance to a range of fundamental science, as well as the practical applications of space weather. Active region emergence and life cycles are two areas of particular interest, yet the lack of consistent full-Sun observations has made long-term studies of active regions difficult. Here, we present results from a study to identify and characterize long-lived active regions (LLARs), defined as those which were observed during at least two consecutive Carrington rotations and which did not undergo significant successive flux emergence once the decay phase began. Such active regions accounted for 13% of all NOAA-identified ARs between 2011 and 2019, and their distribution closely follows the annual sunspot number. This implies that LLARs are produced by the same basic driving processes as regular ARs. LLAR areas tend to be significantly larger and contain more magnetic flux compared to other ARs, but the two categories have similar magnetic complexity distributions. The most striking result, however, is that LLARs are 3-6 times more likely than other ARs to be the source of a solar flare of GOES class C or greater. This highlights the importance of studying what makes a LLAR and how to identify them at emergence with a view towards improved space weather forecasting. The further implications of these findings for AR heating spatial and temporal patterns will be explored in an upcoming study.

Neutrinos oscillate over cosmic distances. Using 11.4 years of IceCube data, the flavor composition of the all-sky neutrino flux from 5\,TeV--10\,PeV is studied. We report the first measurement down to the $\mathcal{O}$(TeV) scale using events classified into three flavor-dependent morphologies. The best fit flavor ratio is $f_e:f_{\mu}:f_{\tau}\,=\,0.30:0.37:0.33$, consistent with the standard three-flavor neutrino oscillation model. Each fraction is constrained to be $>0$ at $>$ 90\% confidence level, assuming a broken power law for cosmic neutrinos. We infer the flavor composition of cosmic neutrinos at their sources, and find production via neutron decay lies outside the 99\% confidence interval.

We present evidence of an extreme scattering event (ESE) toward PSR J2313+4253 using high-cadence observations taken with the Green Bank Observatory 20m telescope. The high density of observations in time allow for detailed tracking of the event. We observe a pair of caustic spikes along with the characteristic drop in scintillation bandwidth that is expected during an ESE. This pattern implies that the structures predominantly responsible for scattering occur at different distances than those from previous and subsequent epochs. A secondary spectrum processed during the event shows a detached feature similar to those found in double lensing events from previously observed ESEs. We measure this event as originating from a scattering region with a distance of 1.04(1) kpc, a transverse size of 15 AU, and a duration of approximately 220 days. We model the event using a Gaussian plasma lens, which, when compared to a pure-noise model, is 197 times more probable. These rare events provide opportunities to study the properties of small-scale structures in the ISM.

We present the first high-resolution X-ray spectrum of NGC 7213 obtained with XRISM/Resolve, supported by simultaneous XMM-Newton, NuSTAR, and SOAR optical data. The XRISM spectrum resolves the neutral Fe\,K$\alpha$ into two components: a narrow core ($\rm FWHM = 650_{-220}^{+240}\,\rm km\,s^{-1}$) consistent with emission at the dust sublimation radius, and a broader, asymmetric line best described by disk-like emission from $\sim 100\,\rm R_{g}$. This disk component mirrors the profile of the double-peaked H$\alpha$ line observed in the optical. In addition, we detect broadened Fe XXV and Fe XXVI emission lines whose inferred locations bridge the gap between the inner disk and the optical broad-line region. The weak narrow Fe K$\alpha$ equivalent width ($\rm EW = 32 \pm 6\,eV$) and absence of a Compton hump imply a low-covering-fraction, Compton-thin torus. Together, these results reveal a radially stratified structure in NGC 7213, spanning nearly four orders of magnitude in radius, and place the source in an intermediate accretion state ($\rm \lambda_{Edd} = 0.001-0.01$) where the inner disk and BLR remain, while the torus shows signs of dissipation.

Vela X-1 is a well-studied accreting X-ray pulsar, with a distinctive pulse profile that has been found to be very similar in different observations spread out over decades. On the other hand, significant variations down to the timescale of individual pulses have been observed. The physical mechanisms leading to the energy-resolved pulse profile and its variations are not fully understood. Long, uninterrupted observations of Vela X-1 with XMM-Newton in 2000, 2006 and 2019 at different orbital phases allow us to study variations of the pulse properties in the soft X-ray range on all timescales in detail. We aim to characterize and quantify the variations of pulse profiles and individual pulse cycles on all timescales probed, and to identify possible factors driving the observed variations on these timescales. We generated consistent energy-resolved pulse profiles for each observation, as well as profiles built from subsets of individual pulse cycles selected by time, flux, or similarity to the mean profiles. We identified five pulsed phases based on the profile morphology and hardness, and examined the relative contributions over time. To quantify short-timescale variability, we compared individual pulse cycles with synthetic light curves derived from scaled versions of the average profiles. The pulse profile of Vela X-1, when averaged over many pulse cycles, remains remarkably stable, as expected. The most prominent variations between epochs are attributable to changes in absorption. Residual systematic differences are primarily flux-dependent, with profiles showing less variability at higher flux levels. On shorter timescales, most individual pulse cycles resemble the average profile, even though significant, sporadic deviations are also present.

This study presents a novel multi-static radar technique for space debris characterisation using micro-Doppler signatures, developed within the Southern Hemisphere Asteroid Radar Programme (SHARP). The method employs C-band continuous waveforms transmitted from NASA's Deep Space Network (DSN) Canberra station, with reflected signals captured by distributed ground-based telescopes converted from astronomical to radar operations. From March 2021 to September 2025, we conducted systematic observations of 20 distinct rocket bodies at various orbital altitudes and object sizes. Micro-Doppler analysis of axial rotation signatures enabled extraction of rotation periods, object dimensions, surface characteristics, and mass distribution parameters with enhanced accuracy in orbital refinement. We implemented advanced imaging reconstruction techniques to generate morphology data of the targets. Results demonstrate successful determination of debris rotation periods with second-level precision, dimensional estimates within 10\% accuracy compared to known specifications, and improved orbital parameter determination reducing position uncertainties by up to 30\%. The technique provides a cost-effective approach for space situational awareness, leveraging existing astronomical infrastructure to enhance Australia's sovereign capabilities in tracking both near-Earth objects and space debris. This multi-static configuration offers significant advantages over traditional monostatic radar systems for geostationary and highly elliptical orbit debris monitoring.

Observations of circumbinary discs reveal inner cavities, with their shape and size varying strongly between different systems. The structure of the cavity is determined by the complex interplay between spirals induced by tidal forcing from the binary and the viscous and radiative damping of the spirals at the cavity edge. To fully understand what determines the properties of observed cavities, it is therefore necessary to capture the effect of radiative processes in modelling. To this end, we run 27 simulations of circumbinary discs in 2D using the PLUTO code. These simulations include various size scales, binary eccentricities and thermodynamic models. We find that the diverse cavity shapes are a natural outcome of the radially-varying cooling timescale, as different radiative processes mediate cooling at different disc size regimes. For binaries with separation of a few au, where the cooling timescale is comparable to the orbital timescale at the cavity edge, we recover much more circular cavities than for quickly- or slowly-cooling discs. Our results show that the cavity structure around several binary systems such as Cs Cha and GG Tau can be explained with one physical model, and highlight the importance of radiative cooling in modelling the dynamical evolution of circumbinary discs.

Supernova (SN) 1987A, with its long-rising ($\gtrsim$40~days) light curve, defines a rare subclass of type II SNe known as 1987A-like events. Representing only $\sim$1-3\% of all core-collapse SNe and often found in low-metallicity environments, their large diversity suggests a wide range of progenitor and explosion properties. This study aims to improve the understanding of 1987A-like SNe by characterizing their explosion parameters, including kinetic energy, ejected mass, progenitor radius at explosion, and synthesized $^{56}$Ni mass. Additionally, it seeks to identify systematic trends in both the physical properties and the observed features of these peculiar events. A new Bayesian parameter estimation method, based on our $^{56}$Ni-dependent analytical model for hydrogen-rich SNe, is applied to derive explosion parameters from the light curves and expansion velocities of one of the largest and most comprehensive 1987A-like SN samples to date. These data are measured through a consistent analysis of observations available in the literature. The analysis reveals a heterogeneous population that nevertheless clusters into two main groups: (i) lower-energy explosions with modest $^{56}$Ni yields ($\sim$0.07~M$_\odot$), similar to SN~1987A, and (ii) more energetic events (up to $\sim$5~foe) with larger nickel production and, in some cases, unusually extended progenitors. We confirm a robust correlation between $^{56}$Ni mass, peak luminosity, and explosion energy, as well as between ejecta mass and recombination timescale. An anti-correlation between Ba~II line strength and photospheric velocity indicates that stronger Ba~II absorptions in 1987A-like SNe arise from more compact, slowly expanding ejecta. Our study underscores the need to extend analytical frameworks to include additional power sources, enabling scalable and accurate modelling of the growing number of peculiar transients.

The origin of radio afterglows or delayed radio flares in tidal disruption events (TDEs) is not fully understood. They could be generated either by a forward shock (FS) propagating into diffuse circumnuclear medium (CNM), or a bow shock (BS) around a dense cloud, each of which is fundamentally different. To elucidate the distinctions between these two scenarios, we conducted two-fluid simulations incorporating relativistic electrons to investigate the spatial evolution of these electrons after being accelerated by shock. Based on their spatial distribution, we performed radiative transfer calculations to obtain the synchrotron spectra. In Paper I (Mou 2025), we reported the results for the FS scenario; in this article, we focus on the BS scenario. Compared to that from the FS, the radio emission from the BS exhibits a higher peak frequency, and its flux shows a much steeper rise and a more rapid decline. The radio flux from the BS also responds to fluctuations in the outflow. The combined effects of the BS and FS substantially alter radio spectra, causing significant deviations from the single-zone emission model, and in some cases producing double-peaked or flat-top features in spectra. This study highlights the importance of the BS, and inspires a novel approach for probing dense gas on sub-parsec scales in galactic nuclei by decomposing the BS radio spectrum to reveal the conditions of circumnuclear dense gas.

Interstellar comet 3I/ATLAS has been approaching its 2025 October 29 perihelion while opposite the Sun from Earth, hindering ground-based optical observations over the preceding month. However, this geometry placed the comet within the fields of view of several space-based solar coronagraphs and heliospheric imagers, enabling its continued observation during its final approach toward perihelion. We report photometry from STEREO-A's SECCHI HI1 and COR2, SOHO's LASCO C3, and GOES-19's CCOR-1 instruments in 2025 September--October, which show a rapid rise in the comet's brightness scaling with heliocentric distance r as r^(-7.5+/-1.0). CCOR-1 also resolves the comet as an extended source with an apparent coma ~4' in diameter. Furthermore, LASCO color photometry shows the comet to be distinctly bluer than the Sun, consistent with gas emission contributing a substantial fraction of the visible brightness near perihelion.

The 2025 June 01 Forbush Decrease in the terrestrial ground-level flux of cosmic ray secondaries was recorded by many cosmic ray systems. This was the deepest such decrease, from the quiescent value of the flux, which has been observed in the past two decades. It resulted from a complex series of solar events, none of which on its own reached the most extreme level. The extreme depth of this Decrease has enabled measurements of the flux reduction to be made, which would normally be severely limited by particle counting statistics. In particular, here we examine the Decrease phenomenon over a primary cosmic ray energy range which is rarely accessible, due to the low flux of high energy cosmic rays. This work considers data mainly from a muon telescope system which can respond to both unaccompanied muons and small cosmic ray air showers, providing data from GeV to mid-TeV energies, where the Forbush Decrease ceases to be statistically observable. This paper examines the depth of the flux Decrease, as a fraction below its quiescent value, over that primary energy range. The progressive development of internal time structure in the flux through the seven days of the phenomenon is also demonstrated.

Magnetic reconnection is a proposed mechanism for nanojets associated with coronal heating. We investigate the characteristics of reconnection-driven nanojets just before and during a prominence eruption using the High Resolution Imager (HRI) of the Extreme Ultraviolet imager (EUI) aboard Solar Orbiter during its perihelion on September 30, 2024. Extreme UV (EUV) images at unprecedented high spatial and temporal resolution from \hrieuv were analysed. The dimensions and propagation speeds of nanojets were estimated and used to estimate the kinetic energies. Nanojet activity was compared with GOES X-ray flux to assess its relation to flare evolution. The high spatial and temporal resolution in the EUV was found to be essential to fully capture the properties and numbers of reconnection nanojets. Approximately 120 nanojets were detected during the eruption, with 40 analysed in detail. Nanojets exhibited lengths of $200 - 5000$~km, widths of $200 - 500$~km and durations of $2-12$~s. Instant velocities ranged from 150 km~s$^{-1}- 600~$km~s$^{-1}$ with kinetic energies reaching $1.56\times10^{27}$~erg. These nanojets are faster, longer, more energetic and more numerous compared to previous studies. We also find clear signatures of acceleration and deceleration, reflecting magnetic tension release and reach of new equilibria. Reconnection events during the eruption were found to be more frequent and energetically intense. Pre-flare nanojet clustering indicates small-scale reconnection may precede large eruptive activity. These results suggest that nanojets also occur in fully ionised coronal plasma, playing a role in both quiescent and eruptive solar activity.

We investigate the magnetic field orientation and strength in the massive star-forming region G35.20-0.74N (G35), using polarized dust emission data obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) as part of the Magnetic fields in Massive star-forming Regions (MagMaR) survey. The G35 region shows a filamentary structure (a length of $\sim$0.1 pc) with six bright cores located along the filament's long axis. Magnetic field strengths across the G35 region range from 0.2 to 4.4 mG with a mean value of 0.8 $\pm$ 0.4 mG. The mass-to-flux ratio ($\lambda$) varies from 0.1 to 6.0 the critical value. The highest values are found locally around cores, whereas the remains of the filament are subcritical. A H$^{13}$CO$^+$ (3--2) velocity gradient of 29 km s$^{-1}$ pc$^{-1}$ is evident along the filament's long axis, aligned with the magnetic field direction. At larger scales ($\sim$0.1 pc), the magnetic field lines appear roughly perpendicular to the filament's long axis, in contrast to the smaller-scale structure ($\sim$0.003 pc) traced by ALMA. The magnetic field lines could be dragged along the filament as a result of the gas motion induced by the gravitational potential of the filament. Six cores in the filament have similar spacings between 0.02--0.04 pc. The initial filament fragmentation could have produced a core spacing of 0.06 pc, following filament fragmentation theory, and the current core spacing is the result of cores comoving with the gas along the filament. This core migration could occur in a few 10$^4$ years, consistent with high-mass star formation time scales.

GRS 1915$+$105 was the stellar-mass black hole that best reproduced key phenomena that are also observed in Type-1 active galactic nuclei. In recent years, however, it has evolved to resemble a Type-2 or Compton-thick AGN. Herein, we report on the first XRISM observation of GRS 1915$+$105. The high-resolution Resolve calorimeter spectrum reveals that a sub-Eddington central engine is covered by a layer of warm, Compton-thick gas. With the obscuration acting as a coronagraph, numerous strong, narrow emission lines from He-like and H-like charge states of Si, S, Ar, Ca, Cr, Mn, Fe, and Ni dominate the spectrum. Radiative recombination continuum (RRC) features are also observed, signaling that much of the emitting gas is photoionized. The line spectrum can be fit by three photoionized emission zones, with broadening and bulk velocities suggestive of an origin in the outer disk atmosphere and/or a slow wind at $r \simeq 10^{6}~GM/c^{2}$. The Fe XXV He-$\alpha$ and Fe XXVI Ly-$\alpha$ lines have a broad base that may indicate some emission from $r \sim 3\times 10^{3}~GM/c^{2}$. These results broadly support a picture wherein the current state in GRS 1915$+$105 is due to obscuration by the irradiated outer disk. This could arise through disk thickening if the Eddington fraction is higher than inferred, but it is more likely due to a warped, precessing disk that has brought the outer disk into the line of sight. We discuss the strengths and weaknesses of this interpretation and our modeling, and possible explanations of some potentially novel spectral features.

Detecting gaseous debris disks around white dwarfs offers a unique window into the ultimate fate of planetary systems and the composition of accreted planetary material. Here we present a systematic search for such disks through the Ca II infrared triplet using the Dark Energy Spectroscopic Instrument (DESI) Early Data Release. From a parent sample of 2706 spectroscopically confirmed white dwarfs, we identify 22 candidate systems showing tentative emission-line features, which corresponds to a raw occurrence rate of 0.81%, more than ten times higher than previous estimates. The detected emission lines are predominantly weak and require confirmation by follow-up observations. Three of these candidates also exhibit infrared excess in WISE photometry, suggesting a possible coexistence of gas and dust. However, the high candidate rate indicates that most are likely false positives due to telluric residuals or unresolved binaries. This work demonstrates the potential of DESI spectra for blind searches of rare circumstellar phenomena. The recently released DESI DR1, with its substantially larger spectroscopic sample, will enable searches for more gaseous disks and provide better insights into their occurrence and nature.

We investigate the long-term variability of the known Changing Look Active Galactic Nuclei (CL AGN) Mrk 1018, whose second change we discovered as part of the Close AGN Reference Survey (CARS). Collating over a hundred years worth of photometry from scanned photographic plates and five modern surveys we find a historic outburst between ~1935-1960, with variation in Johnson B magnitude of ~0.8 that is consistent with Mrk 1018's brightness before and after its latest changing look event in the early 2010s. Using the combined modern and historic data, a Generalised Lomb-Scargle suggests broad feature with P = 29-47 years. Its width and stability across tests, as well as the turn-on speed and bright phase duration of the historic event suggests a timescale associated with long-term modulation, such as via rapid flickering in the accretion rate caused by the Chaotic Cold Accretion model rather than a strictly periodic CL mechanism driving changes in Mrk 1018. We also use the modern photometry to constrain Mrk 1018's latest turn-off duration to less than ~1.9 years, providing further support for a CL mechanism with rapid transition timescales, such as a changing mode of accretion.

Accretion near black holes (BHs) is multidimensional, with turbulence, resistivity, thermal transport, and dust dynamics all playing essential roles. In cold accretion discs (ADs) or the region of an AD where magnetic fields (MFs) are negligible (or absent), hydrodynamic (HD) turbulence is probably dominating. However, Magneto-rotational instability (MRI) is the primary cause of turbulence in ADs. Significant velocity variations and rapid pressure changes are characteristics of turbulent flows, which allow better mixing and more angular momentum (AM) and energy transfer. Also, in accretion flow (AF), the interaction between turbulence and resistivity determines the efficiency of energy dissipation and heat transfer. Radiation, convection, and thermal conduction (TC) are the heat transport modes existing in AFs, where TC enables energy transfer in accreting materials via heat flux. Moreover, convection, also generated by turbulence, significantly impacts the stability of the AD and its vertical structure. The disc may be affected by radiation from the AD surrounding the BH when X-ray emission occurs. The emission from the disc is also affected by dust particles. Dust grains near BH are exposed to high temperatures and intense radiation, which might affect the flow characteristics, as seen in Active Galactic Nuclei (AGN). This chapter highlights the combined effect of turbulence, resistivity, transport mechanisms, and dust particles on BH AF. Future studies in this field must thoroughly investigate how dust, transport mechanisms, and turbulence interact in the BH accretion system.

This paper presents almanac and ephemeris message representation for lunar navigation satellites supporting the Lunar Augmented Navigation System (LANS). The proposed method combines osculating orbital elements, Chebyshev polynomials, and Fourier series to efficiently represent lunar satellite trajectories subject to complex perturbations from lunar gravity and third-body effects. For the ephemeris, a hybrid Chebyshev--Fourier formulation improves fitting accuracy over long arcs while maintaining message compactness under the data-size constraint of the LunaNet Interoperability Specification. For the almanac, a compact low-order polynomial and Fourier model is introduced to capture mid-term orbital variations over a 15-day fitting arc. The approach is validated for multiple orbit regimes, including 30-hour, 24-hour, and 12-hour elliptical lunar frozen orbits (ELFOs) and a 6-hour polar orbit. Results show that the proposed framework achieves sub-meter position and sub-millimeter-per-second velocity fitting errors within the 900-bit limit for 6-hour ephemeris arcs, and almanac fitting accuracy sufficient for reliable satellite-visibility identification in warm-start operations.

Observational studies of HII region-molecular cloud interactions constrain models of feedback and quantify its impact on the surrounding environment. A recent hypothesis proposes that a characteristic spectral signature in ground state hyperfine lines of hydroxyl (OH) -- the OH flip -- may trace gas that is dynamically interacting with an expanding HII region, offering a new means of probing such interactions. We explore this hypothesis using dedicated Jansky Very Large Array (VLA) observations of three Galactic HII regions, G049.205$-$0.343, G034.256+0.145 and G024.471+0.492, in 1--2 GHz continuum emission, all four 18-cm ground-state OH lines, and multiple hydrogen radio recombination lines. A Gaussian decomposition of the molecular gas data reveals complex OH emission and absorption across our targets. We detect the OH flip towards two of our sources, G049.205$-$0.343 and G034.256+0.145, finding agreement between key predictions of flip hypothesis and the observed multi-wavelength spectra, kinematics and morphology. Specifically, we demonstrate a strong spatial and kinematic association between the OH flip and the ionized gas of the HII regions -- the first time this has been demonstrated for resolved sources -- and evidence from $^{13}$CO(1--0) data that the expected OH component originates from the non-disturbed gas of the parent cloud. While we detect no flip in G024.471+0.492, we do find evidence of interacting molecular gas traced by OH, providing further support for OH's ability to trace HII region-molecular cloud interactions.

This study investigates the origin of gas fueling secondary star formation, i.e., rejuvenation in nearby galaxies. From the MaNGA IFU survey, we use stellar absorption features D$_n$4000 and H$\delta_A$ to identify regions that started the rejuvenation within the last $\sim$200~Myr. We compare the gas-phase metallicity, metallicity gradients, environments, and H\Romannum{1} gas fractions of the rejuvenating galaxies (RJGs) to controlled star-forming and quiescent galaxy samples. We demonstrate that, for the majority of RJGs, the rejuvenating gas is originally in the galaxy rather than accreted gas. The evidence includes: (1) gas metallicities consistent with the mass-metallicity relation of SF galaxies; (2) metallicity gradients that are not flattened, arguing against radial inflows; (3) gas velocities in rejuvenating regions consistent with their surroundings, and (4) high H\Romannum{1} gas fractions comparable to SF galaxies, indicating a pre-existing reservoir. Furthermore, we find no evidence that the rejuvenating events are triggered by tidal interactions with neighbors. While internal processes appear to dominate, we also present a clear example of rejuvenation triggered by gas accretion. The galaxy MaNGA 12080-12705 hosts a low-metallicity, kinematically distinct star-forming region in an overall old, massive galaxy, providing unambiguous evidence of an external origin, such as accretion or a minor merger. Our analysis demonstrates that using D$_n$4000 and EW(H$\delta_A$) provides a reliable way to identify current rejuvenation events in large spectroscopic surveys. The method will enable statistical studies to understand rejuvenation across cosmic time.

We have observed the low-mass molecular cloud core G204.4-11.3A2-NE (G204NE) in the direction of Orion B giant molecular cloud with the Atacama Large Millimeter/submillimeter Array in Band 6. The $\rm 1.3\ mm$ continuum images and visibilities unveil a compact central structure with a radius of $\sim$12 au, while showing no signature of binarity down to 18 au. The bolometric temperature and luminosity of this source are derived to be ${\sim}$33 K and $\sim$1.15 $~L_{\odot}$, respectively. Chemical stratification is observed in dense gas tracers, with C$^{18}$O emission peaking at the continuum position surrounded by the spatially extended emission of N$_2$D$^+$ and DCO$^+$. This implies that the core is in a very early evolutionary stage in which $\rm CO$ depletion occurs in most regions except for a small area heated by the central source. The envelope kinematics indicates a rotating and infalling structure with a central protostar mass of 0.08-0.1 $M_{\odot}$. The protostar drives a collimated outflow traced by CO, SiO, SO, and H$_2$CO, with misaligned blueshifted and redshifted lobes exhibiting a pair of bow-like patterns. High-velocity jets, extending up to 720 au, are detected in CO, SiO, and SO lines. The jet launching region is likely within twice of the dust sublimation zone. The absence of a binary signature suggests the outflows and jets are driven by a single protostar, although a close binary cannot be ruled out. The observed deflection of the outflows and jet is likely due to turbulent accretion in a moderately magnetized core.

X-ray observations of the tidal disruption event (TDE) candidate AT 2019avd show drastic variabilities in flux and spectral shape over hundreds of days, providing clues on the accretion disc-corona evolution. We utilize a disc-corona model, in which a fraction of the gravitational energy released in the disc is transported into the hot corona above/below. Some soft photons emitted from the disc are upscattered to X-ray photons by the hot electrons in the optically thin corona. By fitting the NICER observations of AT 2019avd during epochs when the spectra exhibit significant hardening, we derive the evolution of the mass accretion rate, $\dot{m}$, and the coronal energy fraction, $f$. Our results show that $f$ decreases with increasing $\dot{m}$, which is qualitatively consistent with that observed in active galactic nuclei (AGNs), while the slope of this source, $f\propto \dot{m}^{-0.30}$, is much shallower than that of AGNs. We also find that the non-thermal X-ray spectrum in this source is significantly softer than those typically seen in AGNs and black-hole X-ray binaries. We argue that these quantitative differences can be a powerful diagnostic of the underlying magnetic turbulence, which may imply a stronger magnetic field within the TDE accretion disc than that in typical AGNs. It is also found that the evolution of the fitted neutral hydrogen column density follows a similar pattern to that of the accretion rate evolution, which may reflect the accumulation of absorbing material originating from the inflowing streams of stellar debris and/or other related sources.

We propose a dual-resolution prescription meant for the Boltzmann neutrino transport, in which the advection and collision terms are calculated with different angular resolutions in momentum space. The purpose is to address the issue of the low resolution that afflicts the $S_N$ method in the multi-dimensional neutrino transport simulations for core-collapse supernovae. We handle with a high resolution the advection term alone, assuming that the collision term does not require such high resolutions. To confirm this surmise as well as our new conversion scheme, from low- to high-angular resolutions and vice versa, we run a couple of experimental one-zone (in space) simulations. Neutrino scatterings on nucleons are considered with small recoils fully taken into account whereas the advection term is replaced by the angle- and energy-dependent source terms that are designed to mimic the results of a Boltzmann simulation, inducing the anisotropy in momentum space. For the conversion from a low-resolution distribution function to a high-resolution one, we employ a polynomial interpolations in the zenith and azimuth directions separately with the number conservation and continuity (and periodicity only in the azimuth direction). We find that this dual-resolution scheme works well and that the current angular resolution employed in the canonical supernova simulations with our Boltzmann solver or a bit better in the $\phi_\nu$ direction will be sufficient for the collision terms if they are coupled with the advection terms calculated with a high-angular resolution via this prescription.

We show that the merger of oscillons results in a broad spectrum of the oscillon mass. A huge number of oscillon samples obtained from numerical lattice simulations reveal that the oscillon mass distribution has an exponential tail in a heavy-mass region. This enables us to infer the fractional abundance of heavy oscillons. Using the criterion for the primordial black hole (PBH) formation from the oscillon collapse obtained in previous studies, we estimate the abundance of PBHs and conclude that a sizable number of PBHs can be produced from oscillons. It can be an alternative PBH formation mechanism without employing the tuning of the inflaton potential to enhance the small-scale density fluctuations in the conventional PBH formation scenario.

Planetary rings are not only ubiquitous around the giant planets in the outer Solar System, but have also been discovered around several small distant bodies. In contrast, no rings have been observed around any inner Solar System objects. To constrain the dynamical origin of this ringed-versus-ringless dichotomy, we employ a numerically cross-checked analytical model of gigayear-scale Poynting-Robertson (PR) drag due to the solar flux acting on an isolated particle, expressed as a function of the host body's heliocentric distance \(\,a_{\mathrm{pla}}\) and the particle radius \(\,r_{\mathrm{par}}\). Here we show that, in the absence of additional perturbations, PR drag alone can explain the observed ring architecture of the Solar System: outer planets and Centaurs/TNOs are able to retain rings for the age of the Solar System, whereas any rings around the inner planets are removed on much shorter timescales. Because the PR-drag lifetime scales steeply with heliocentric distance \(\bigl(\tau_{\mathrm{decay}}\propto a_{\mathrm{pla}}^{2} \,r_{\mathrm{par}}\bigr)\), we predict that forthcoming surveys will reveal an ever-growing population of ring-bearing bodies in the distant Solar System.

The Compton Spectrometer and Imager (COSI) is a Compton telescope designed to survey the 0.2 - 5 MeV sky, consisting of a compact array of cross-strip germanium detectors. It is planned to be launched in 2027 into an equatorial low-Earth (530 km) orbit with a prime mission duration of 2 years. The observation of MeV gamma rays is dominated by background, mostly from extragalactic and atmospheric photon but also from the activation of the detector materials induced by cosmic-ray interactions. Thus, background simulation and identification are crucial for the data analysis. In this work we perform Monte Carlo simulations of the background for the first 3 months in orbit, and we extrapolate the results to 2 years in orbit, in order to determine the build-up of the activation due to long-lived isotopes. We determine the rates of events induced by the background that are reconstructed as Compton events in the simulated COSI data. We find that the extragalactic background photons dominate at low energies (<660 keV), while delayed activation from cosmic-ray primaries (proton/alpha) and albedo photons dominate at higher energies. As part of this work, a comparison at low latitude (<1 deg) between recent measurement of the SAA by the High-Energy Particle Detector (HEPD-01) on board the China Seismo-Electromagnetic Satellite (CSES-01) and the AP9/AE9 model has been made, showing an overestimation of the flux by a factor 9 by the model. The systematic uncertainties associated with these components are quantified. This work marks a major step forward in estimating and understanding the expected background rates for the COSI satellite mission.

Circumbinary configurations offer a test of planet formation in an altered environment, where the inner binary has perturbed a protoplanetary disc. Comparisons of the physical and orbital parameters between the circumbinary planet population and the population of exoplanets orbiting single stars will reveal how these disc perturbations affect the assembly of planets. Circumbinary exoplanets detected thus far typically have masses $< 3 \,\rm M_{jup}$ raising the question of whether high-mass circumbinary planets are possible, and also whether population features such as the brown dwarf desert would appear in circumbinary configurations like for single star systems. Here, we report observations taken with the SOPHIE high-resolution spectrograph. These observations reveal an $m_{\rm b}\,\sin i_{\rm b} = 20.9 \,\rm M_{jup}$ outer companion, on an eccentric ($e = 0.43$), $1800\,\rm d$ orbit, which we call BEBOP-4 (AB) b. Using dynamical arguments we constrain the true mass $m_{\rm b}< 26.3 \,\rm M_{jup}$. The inner binary's two eclipsing stellar components have masses $M_{\rm A} = 1.51\,\rm M_\odot$, and $M_{\rm B} = 0.46\,\rm M_\odot$. Their orbital period is $72\,\rm d$, and their eccentricity is $0.27$. This system contains the longest period binary surveyed by the BEBOP project. BEBOP-4b is expected to be detectable using Gaia DR4 single epoch astrometric measurements. Despite a large period ratio of $\sim 25:1$, the substantial eccentricities of both orbits mean that the outer orbit is on the edge of orbital stability, and located in between two destabilising secular resonances. Should the outer companion survive, the BEBOP-4 system appears like a precursor to several post-common envelope binaries exhibiting eclipse timing variations where very massive circumbinary companions have been proposed.

Context. Current exoplanet surveys are focused on detecting small exoplanets orbiting in the liquid-water habitable zones of their host stars. Despite the recent significant advancements in instrumental developments, the current limitation in detecting these exoplanets is the intrinsic variability of the host star itself. Aims. Our aim is to use the full CARMENES guaranteed time observations (GTO) data set spanning more than 8 years of observations of over 350 stars to investigate the wavelength dependence of high-precision radial velocities (RV), as stellar activity features should exhibit a wavelength dependence while the RV variation due to an orbiting planet will be wavelength independent. Methods. We use the chromatic index (CRX) to quantify the slope of the measured RVs as a function of logarithmic wavelength. We investigate the dependence of the CRX in the full CARMENES GTO sample on 24 stellar activity indices in the visible and near-infrared channels of the CARMENES spectrograph and each star's stellar parameters. We also present an updated convective turnover time scaling for the calculation of the stellar Rossby number for M dwarfs. Results. Our results show that approximately 17\% of GTO stars show a strong or a moderate correlation between CRX and RV. We can improve the measured RVs by a factor of up to nearly 4 in rms by subtracting the RV predicted by the CRX-RV correlation from the measured RVs. Mid M dwarfs with moderate rotational velocities, moderate CRX-gradients and quasi-stable activity features have the best rms improvement factors. Conclusions. We conclude that the CRX is a powerful diagnostic in mitigation of stellar activity and the search for low mass rocky planets.

We review major achievements in our understanding of multiphase coronal plasma, where cool-dense and hot-tenuous matter coexists, brought about by advances in modeling and theory, inspired by observations. We give an overview of models that self-consistently form solar (or stellar) prominences and filaments, or (postflare) coronal rain, and clarify how these different phenomena share a common physical origin, relating radiative losses and coronal heating. While we do not fully understand the coronal heating, multi-dimensional models of solar prominence and rain formation demonstrate how thermal instability triggers condensations, and how their morphology may reveal aspects of the applied heating at play. We emphasize how the many pathways to linear instability due to combined ingredients of heat-loss, gravity, flows, and magnetic topologies are all involved in the resulting nonlinear magnetohydrodynamics. We provide some challenges to future model efforts, especially concerning prominence fine structure, internal dynamics, and their overall lifecycle.

Recent observations with the James Webb Space Telescope (JWST) have begun to reveal a surprising morphological diversity in galaxies within the first billion years after the Big Bang, including indications of structural maturity previously thought to arise much later. These findings raise fundamental questions about when and how well-known structural components of galaxy morphology, such as bulges and disks, first emerged. However, directly identifying and resolving such structures at z $>$ 6 remains challenging due to limited spatial resolution and sensitivity. In this work, we present a clear and robust morphological analysis of a sample of 190 galaxies at z $>=$ 6, demonstrating that distinct bulge and disk components were already beginning to emerge during this early epoch. Using multi-component light profile fitting, we model the radial brightness distributions of a subset (20) of galaxies with an inner spheroidal (Sersic) component and an underlying exponential disk. These systems exhibit high bulge-to-total (B/T) light ratios (~ 0.47) and central stellar mass surface densities (~ 2.82*10$^{8}$ M$_{sun}$ kpc$^{-2}$ ) - values close to those of nearby quiescent galaxies. Combined with their intense central star formation rate surface densities (~ 1.26*10$^{1}$ M$_{sun}$ yr$^{-1}$ kpc$^{-2}$ ), our results indicate a rapid building of inner stellar mass and bulge assembly within these young systems. We propose that these early bulge-disk galaxies represent progenitors of massive star-forming and quiescent systems observed at lower redshifts. Their subsequent evolution may proceed through physical processes such as disk growth, compaction, quenching, or bulge-disk co-evolution, driven by both internal dynamics and external interactions.

The effective redshift distribution $n(z)$ of galaxies is a critical component in the study of weak gravitational lensing. Here, we introduce a new method for determining $n(z)$ for weak lensing surveys based on high-quality redshifts and neural network-based importance weights. Additionally, we present the first unified photometric redshift calibration of the three leading stage-III weak lensing surveys, the Dark Energy Survey (DES), the Hyper Suprime-Cam (HSC) survey and the Kilo-Degree Survey (KiDS), with state-of-the-art spectroscopic data from the Dark Energy Spectroscopic Instrument (DESI). We verify our method using a new, data-driven approach and obtain $n(z)$ constraints with statistical uncertainties of order $\sigma_{\bar z} \sim 0.01$ and smaller. Our analysis is largely independent of previous photometric redshift calibrations and, thus, provides an important cross-check in light of recent cosmological tensions. Overall, we find excellent agreement with previously published results on the DES Y3 and HSC Y1 data sets while there are some differences on the mean redshift with respect to the previously published KiDS-1000 results. We attribute the latter to mismatches in photometric noise properties in the COSMOS field compared to the wider KiDS SOM-gold catalog. At the same time, the new $n(z)$ estimates for KiDS do not significantly change estimates of cosmic structure growth from cosmic shear. Finally, we discuss how our method can be applied to future weak lensing calibrations with DESI data.

We present a statistical framework that establishes an accelerating star formation scenario for dense clumps using ATLASGAL and ALMAGAL samples. By employing the cumulative distribution function of dust temperature as a monotonic evolutionary indicator, we linearize clump evolution into a normalized timescale, enabling direct comparison across different samples. The virial mass of clumps increases exponentially with this normalized time, revealing an accelerating buildup of star-forming gas within protoclusters. The evolution of the maximum core mass further shows that the growth timescales of protoclusters and their embedded most massive protostars are comparable, implying a self-similar acceleration of star formation from the stellar to the protocluster scale. This unified framework naturally reproduces the observed evolution of luminosity, the core mass function, the mass growth of the most massive protostars, and the dense gas star formation law on clump scales, establishing a coherent picture of accelerating star formation across scales.

Brown dwarfs display Jupiter-like auroral phenomena, such as rotationally modulated electron cyclotron maser radio emission. Radio observations of cyclotron maser emission can be used to measure their magnetic field strength, topology, and to deduce the presence of magnetically interacting exoplanets. Observations of the coldest brown dwarfs (spectral types T and Y) are especially intriguing, as their magnetospheric phenomena could closely resemble those of gas-giant exoplanets. Here we report observations made over ten epochs, amounting to 44 hours, of WISEP J101905.63+652954.2 (J1019+65 hereinafter) using the LOFAR telescope between 120 and 168 MHz. J1019+65 is a methane dwarf binary (T5.5+T7) whose radio emission was originally detected in a single-epoch LOFAR observation to be highly circular polarised and rotationally modulated at $\approx 3$h. Unexpectedly, our long-term monitoring reveals an additional periodic signature at $\approx 0.787$h. We consider several explanations for the second period and suggest that it could be the rotationally modulated emission of the second brown dwarf in the binary, although follow-up infrared observations are necessary to confirm this hypothesis. In addition, the data also allow us to statistically estimate the duty cycle and observed radio-loud fraction of the 120-168\,MHz cyclotron emission from methane dwarfs to be $\langle D \rangle = 0.030^{+0.034}_{-0.030}$ and $F^{'}_{\rm radio} = 0.088^{+0.168}_{-0.088}$ respectively.

The distribution of the coronal electron density at different distances from the Sun strongly influences the physical processes in the solar corona and is therefore a very important topic in solar physics. Most methods, including radio observations, used for estimating coronal electron density were not fully validated due to the absence of in situ observations closer to the Sun. Consequently, space weather forecasting models that simulate coronal density lacked proper validation. Newly available PSP in situ observations at distances close to the Sun provide an opportunity to study plasma properties near the Sun and to compare observational and modeling results. This work studies type III bursts, estimates their propagation path, and validates coronal electron density obtained from radio, in situ observations, and modeling with EUHFORIA. Type III bursts observed during the second PSP perihelion are analyzed using radio triangulation and modeling. We determine 3D positions of radio sources and use EUHFORIA to estimate electron densities at various locations. The electron densities derived from radio observations and EUHFORIA modeling are inter-validated with in situ PSP measurements. We studied 11 type III bursts during the second PSP perihelion, with radio triangulation showing propagation paths southward from the solar ecliptic plane. Radio source sizes ranged from 0.5 to 40 deg (0.5 to 25 Rs) with no clear frequency dependence, indicating that scattering of radio waves was not very significant. Comparison of electron densities from radio triangulation, PSP data, and EUHFORIA modeling showed a large range of values, influenced by different propagation paths and model limitations. Despite these variations, EUHFORIA identified high-density regions along type III burst paths.

We present the analysis of 900 ks high-resolution RGS XMM-Newton observations of the nearby Seyfert galaxy Mrk 509 with the use of a self-consistent warm absorber (WA) model. We perform a physically motivated approach to the modeling of WA as a stratified medium in constant total pressure (CTP) regime. Powerful outflows are fundamental ingredients of any active galactic nuclei (AGN) structure. They may significantly affect the cosmological environment of their host galaxy. High-resolution X-ray data are best suited for outflow's studies, and the observed absorption lines on heavy elements are evidence of the physical properties of an absorbing gas. Our models allow us to fit continuum shapes bounded together with the line profiles, which gives additional constraints on the gas structure of WA in this source. In this work, we benchmark and test the CTP model on the soft X-ray spectrum of Mrk 509. A grid of synthetic absorbed spectra was computed with the photoionization code Titan assuming that the system was under CTP. As an illuminating spectral energy distribution (SED), we used the most actual multiwavelength observations available for Mrk 509. We apply these models to the high-resolution spectrum of the WA in the Mrk 509, taking into account cold/warm/hot Galactic absorption on the way to the observer. CTP gas with $\log \xi_{0} \sim 1.9$, defined on the cloud surface, fits the data well. A higher ionization component is needed for Ne X absorption. The best-fit model is optically thin with $\log N_{\rm H }= 20.456 \pm 0.016$. The lines are non-saturated, and the CTP spectral fit aligns with previous analyses of Mrk 509 with a constant density WA. The model constrains the gas density, placing the WA cloud at 0.02 pc, consistent with the inner broad line region and the thickening region of the accretion disk.

The Early Dark Energy (EDE) model has been proposed as a candidate mechanism to generate cosmic birefringence through a Chern-Simons coupling between a dynamical scalar field and the cosmic microwave background (CMB) photon. Such birefringence induces a nonzero cross-correlation between the CMB $E$- and $B$-modes, providing a direct observational signature of parity violation. Recent measurements of the $EB$ and $TB$ power spectra, however, cannot yet unambiguously separate instrumental miscalibration ($\alpha$) from a true cosmic-rotation angle ($\beta$). For this reason, we perform a model-independent analysis in terms of the total effective rotation angle $\alpha+\beta$. We analyze the latest $EB$ and $TB$ measurements from the SPIDER, Planck, and ACT experiments and derive constraints on the Chern-Simons coupling constant $gM_{Pl}$ and on the polarization rotation angle $\alpha+\beta$. We find that the coupling $gM_{Pl}$ is not compatible with the SPIDER data, while it provides reasonable fits to the Planck and ACT measurements. The fits for $\alpha+\beta$ prefer a value larger than zero: when combined, Planck+ACT yield a detection significance of approximately 7$\sigma$. We also find that ACT data alone do not provide sufficiently tight constraints on either $gM_{Pl}$ or $\alpha+\beta$, whereas the combination Planck+ACT improves the statistical consistency of ACT's high-$\ell$ results and leads to a better PTE for those measurements.

The astrophysical origin of the rapid neutron-capture process (r-process), which produces about half of the elements heavier than iron, remains uncertain. The oldest, most metal-poor stars preserve the chemical signatures of early nucleosynthesis events and can reveal the nature of the r-process sites. We present a homogeneous chemical abundance analysis of ten r-process-enhanced, metal-poor stars that show strong enrichment in r-process elements with minimal contamination from other nucleosynthetic sources. Using high-resolution, high signal-to-noise spectra, we examined over 1400 absorption lines per star through equivalent width measurements and spectral synthesis under one-dimensional LTE assumptions with the MOOG radiative transfer code. Abundances for 54 chemical species were derived, including 29 neutron-capture elements spanning the full r-process pattern. We quantified the cosmic scatter of elemental ratios relative to Zr (light) and Eu (heavy) and found remarkably small dispersions for the rare-earth and third-peak elements, {\sigma}[La/Eu] = 0.08 dex and {\sigma}[Os/Eu] = 0.11 dex, while the light-to-heavy ratio shows slightly larger variation, {\sigma}[Zr/Eu] = 0.18 dex. A kinematic study indicates that the stars likely originated from ten distinct progenitor systems, allowing us to probe the intrinsic variation between independent r-process events. These results imply that the main r-process operates under highly uniform conditions across diverse astrophysical sites.

The Polarimetric and Helioseismic Imager (SO/PHI) onboard Solar Orbiter is a spectropolarimeter scanning the Fe I line at 617.3 nm, providing data of the solar photosphere. The same line is sampled by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) and many other on-ground instruments. In this paper, we aim at assessing the consistency between line-of-sight (LoS) velocity measurements from the two instruments. Reliable measurements of up and down flows from SO/PHI are crucial and unique when Solar Orbiter is facing the far side of the Sun. Also, a combination of measurements from two vantage points to study horizontal flows must rely on consistent observations. For this purpose, we compare the LoS velocity measured by SO/PHI's High Resolution Telescope (SO/PHI-HRT) and SDO/HMI on 29 March 2023, when Solar Orbiter was crossing the Sun-Earth line at a distance of 0.39 au from the Sun. Because such co-alignments are rare, this configuration offered an almost unique opportunity to directly compare data products from both telescopes. The data are aligned and remapped to allow a pixel-by-pixel comparison of the whole time series of 4 hours length. Temporal and spatial variations are considered for a direct combination of the measurements. The LoS velocity distributions are evaluated and a clear linear relation is found between the two instruments with a slope of 0.94 and a correlation of 90%. We find that the signals form at similar heights, with a separation of 10$\pm$12 km, which is larger than previous estimates. A close-up look at the penumbra of a sunspot and its Evershed flow is presented. We conclude that the signals inferred by SO/PHI-HRT and SDO/HMI show very good agreement and high correlation when instrumental effects and large-scale velocities on the Sun are properly accounted for.

We present a fast method for obtaining fully analytical approximations for gravitational waveforms produced by merging of neutron stars and/or black holes for the earliest stages of the merger process. The obtained analytical formula is compared with numerical calculations, its accuracy and limits of applicability are evaluated. Our results may be useful not only for the earliest evalution of properties of the nature of binary system in gravitational-wave detectors but also will give early alerts for gamma-ray, optical and neutrino observatories.

The scope of this paper is focused on the question if gas is transported to dense clumps inside high-mass star-forming regions through filamentary structures on a scale of several parsecs to the sub-parsec scale. We quantify the gas flows from a scale of up to several parsecs down to the sub-parsec scale along filamentary structures. In this work the accretion flow mechanisms based on gas kinematic data in the three high-mass star-forming regions G75.78, IRAS21078+5211 and NGC7538 are studied with data obtained from the IRAM 30 m telescope. The analysis is carried out using the surface density derived from 1.2 mm continuum emission and velocity differences estimated from HCO+ and H13CO+ molecular line data.

Recent observations suggest that planetary engulfment by a giant star may produce radiation that resembles subluminous red novae. We present three-dimensional hydrodynamical simulations of the interaction between an eccentric $5 \,M_J$ giant planet and its $1\,M_\odot$ red-giant host star. The planet's pericenter is initially $60\%$ of the stellar radius and is fully engulfed after tens of orbits. Once inside the stellar envelope, the planet generates pressure disturbances that steepen into shocks, ejecting material from the envelope. We use post-processing to calculate the light curves produced by planetary engulfment. We find that the hot stellar ejecta enhances the stellar luminosity by several orders of magnitude. A prolonged hydrogen recombination plateau appears when the ejecta cools to about $10^4\,\rm{K}$. The late-time rapid dimming of the light curve follows dust formation, which obscures the radiation. For planets with lower eccentricity, the orbital decay proceeds more slowly, although the observable properties remain similar.

A review of publications devoted to the study of the characteristics of the Radcliffe wave has been given. The advent of mass measurements of radial velocities of stars has recently led to a number of interesting results obtained from the analysis of spatial velocities of stars and open star clusters. An important place in the study has been given to issues related to the clarification of the direct or indirect influence of magnetic fields on the process of formation of the Radcliffe wave. The hypothesis of Parker instability of the galactic magnetic field as one of the reasons for the formation of wave-type inhomogeneities in the galactic disk has been discussed.

The LIGO-Virgo-KAGRA collaborations have reported gravitational wave signals from more than 150 binary black holes in the fourth catalog (GWTC-4). Here, we investigate the population properties of these binary black holes using the mixture-model framework Vamana. We present one-dimensional distributions of masses and spins, explore their correlations, and examine their evolution with redshift. These features may reflect astrophysical processes associated with binary black hole formation channels, although most remain poorly constrained. A notable feature is a peak near $10M_\odot$ in the primary mass and $8M_\odot$ in the chirp mass. Additionally, the primary and secondary masses correlate uniquely, producing pronounced peaks in the chirp mass around $14M_\odot$ and $27M_\odot$. The three peaks are roughly separated by a factor of two. A simple explanation for such well-placed peaks is a hierarchical merger scenario, in which the first peak arises from mergers of black holes of stellar origin, and higher-mass peaks arise from repeated mergers of black holes from lower-mass peaks. Although most binaries do not exhibit the high spins and characteristic mass ratios expected from hierarchical mergers, those that do are associated with the peaks observed in the chirp mass distribution.

Aims. eROSITA on board of the Spectrum Roentgen Gamma (SRG) spacecraft performed its first X-ray all-sky survey (eRASS1) between December 2019 and June 2020. It detected about 930000 sources, providing us with an unprecedented opportunity for a detailed blazar census. We present the properties of blazars and blazar candidates in eRASS1 and the compilation of the eROSITA blazar catalog. Methods. We compile a list of blazar and blazar candidates from the literature and match it with the eRASS1 catalog, creating the Blazars in eRASS1 (BlazEr1) catalog. For sources with more than 50 counts we obtain their X-ray spectral properties. We compile multiwavelength data from the radio to the gamma-ray regimes for all sources, including multiwavelength spectral indices and redshifts. The full catalog is available online. Results. We present the BlazEr1 catalog, containing 5865 sources, of which 2106 are associated with confirmed blazars. For 3668 sources, eROSITA provides the first X-ray data. The contamination from non-blazar sources of the entire sample is less than 11%. Most candidates exhibit properties typical for blazars. We present properties of the entire X-ray detected blazar population, including the distributions of X-ray luminosities and photon indices, multiwavelength properties, and the blazar log N-log S distribution. Our catalog provides follow up targets, such as potential MeV and TeV blazars. Conclusions. The BlazEr1 catalog provides a compilation of X-ray detected blazars and blazar candidates. The catalog serves as a starting point for exploiting further eROSITA surveys using the same methodology, enabling us to study the X-ray variability and a large number of spectral energy distributions of blazars in the future.

We combine weak lensing, galaxy clustering, cosmic microwave background (CMB) lensing, and their cross-correlations (so-called 6$\times$2pt) to constrain cosmology and baryonic feedback scenarios using data from the Dark Energy Survey (DES) Y3 Maglim catalog and the Planck satellite PR4 data release. We include all data points in the DES Y3 cosmic shear two-point correlation function (2PCF) down to 2.$^\prime$5 and model baryonic feedback processes via principal components (PCs) that are constructed from the ANTILLES simulations. We find a tight correlation between the amplitude of the first PC $Q_1$ and mean normalized baryon mass fraction $\bar{Y_\mathrm{b}}=\bar{f}_\mathrm{b}/(\Omega_\mathrm{b}/\Omega_\mathrm{m})$ from the ANTILLES simulations and employ an independent $\bar{Y_\mathrm{b}}$ measurement from Akino et al. (2022) as a prior of $Q_1$. We train a neural network $6\times2$pt emulator to boost the analysis speed by $\mathcal{O}(10^3)$, which enables us to run an impressive number of simulated analyses to validate our analysis against various systematics. For our 6$\times$2pt analysis, we find $S_8=0.8073\pm0.0094$ when including a $Q_1$ prior from $\bar{Y_\mathrm{b}}$ observations. This level of cosmological constraining power allows us to put tight constraints on the strength of baryonic feedback. We find $Q_1=0.025^{+0.024}_{-0.029}$ for our 6$\times$2pt analysis and $Q_1=0.043\pm{0.016}$ when combining with external information from Planck, ACT, DESI. All these results indicate weak feedback, e.g., the tensions to Illustris ($Q_1=0.095$) and OWLS AGN T8.7 ($Q_1=0.137$) are 2.9$\sigma$-3.3$\sigma$ and 4.7$\sigma$-5.9$\sigma$, respectively.

Tidal dissipation in star-planet systems occurs through various mechanisms, including the precessional instability. This is an instability of laminar flows (``Poincaré flows") forced by axial precession of a rotating, oblate, spin-orbit misaligned fluid planet or star, which excites inertial waves in convective regions if the dimensionless precession rate (``Poincaré number" $\mathrm{Po}$) is sufficiently large. We constrain the contribution of the precessional instability to tidal dissipation and heat transport, using Cartesian hydrodynamical simulations in a small patch of a planet, and study its interaction with turbulent convection, modelled as rotating Rayleigh-Bénard convection. The precessional instability without convection results in laminar flow at low values and turbulent flow at sufficiently high values of $\mathrm{Po}$. The associated tidal dissipation rate scales as $\mathrm{Po}^2$ and $\mathrm{Po}^3$ in each regime, respectively. With convection, the Poincaré number at which turbulent flow is achieved shifts to lower values for stronger convective driving. Convective motions also act on large-scale tidal flows like an effective viscosity, resulting in continuous tidal dissipation (scaling as $\mathrm{Po}^2$), which obfuscates or suppresses tidal dissipation due to precessional instability. The effective viscosities obtained agree with scaling laws previously derived using (rotating) mixing-length theory. By evaluating our scaling laws using interior models of Hot Jupiters, we find that the precessional instability is significantly more efficient than the effective viscosity of convection. The former drives alignment in 1 Gyr for a Jupiter-like planet orbiting within 23 days. Linearly excited inertial waves can be even more effective for wider orbits, aligning spins for orbits within 53-142 days.

Turbulence is a ubiquitous process that transfers energy across many spatial and temporal scales, thereby influencing particle transport and heating. Recent progress has improved our understanding of the anisotropy of turbulence with respect to the mean magnetic field; however, its exact form and implications for magnetic topology and energy transfer remain unclear. In this Letter, we investigate the nature of magnetic anisotropy in compressible magnetohydrodynamic (MHD) turbulence within low-$\beta$ solar wind using Cluster spacecraft measurements. By decomposing small-amplitude fluctuations into Alfvén and compressible modes, we reveal that the anisotropy is strongly mode dependent: quasi-parallel (`slab') energy contains both Alfvén and compressible modes, whereas quasi-perpendicular (`two-dimensional'; 2D) energy is almost purely Alfvénic, a feature closely linked to collisionless damping of compressible modes. These findings elucidate the physical origin of the long-standing `slab+2D' empirical model and offer a new perspective on the turbulence cascade across the full three-dimensional wavevector space.

Context. Bulges, located at the central regions of galaxies, are complex structures, expected to be shaped by the physical processes involved in the assembly history of their host galaxy, such as gravitational collapse, mergers, interactions, and bars. As a consequence a variety of bulges with distinct morphology and chemistry could be produced. Aim. We aim at exploring the existence of a stellar mass-metallicity relation of bulges, MZ*R, and analyze the possible imprint of characteristics features by accretion and migration of stars, which could store information on their assembly histories. Methods. We use 44 central galaxies from the CIELO cosmological simulations. Their stellar masses are within the range of [10^7.6, 10^10.6] Msun. We decomposed the galaxy into bulge and disk using the circularity and binding energies. We track the stellar populations in bulges back in time to their birth location, classifying them as bulge-born in-situ, and disk-born stars and accreted. Results. We find that most of the stars in our bulges are formed in-situ, but 33% of our bulges show a non-negligible contribution of stellar accretion from satellites, which could add to about 35% of the population. The accreted material is generally contribute by two or three satellites at most. In some bulges, we also find up to a 32% of stars that migrated from the disk due to secular evolution, with a median of 10%. Regardless of the formation histories, we found a clear MZ*R for bulges, which is more enriched by about 0.4 dex than the corresponding relation of the disk components, and about 0.15 dex more enriched than the galaxy MZ*R. We find evidence that the dispersion in the bulge MZ*R is influenced by both stellar accretion from satellites and migration from the disk, such that, at a fixed bulge mass, bulges with higher fraction of accreted and migrated stars tend to be less metal-rich... (abridged)

We investigate the properties of extreme emission line galaxies (EELGs) at $z=4-9$ and their role in reionization. Compact, low-mass galaxies with intense optical emission lines are linked to elevated specific star formation rates (sSFRs) and recent bursts of star formation. Feedback in these systems may enable the leakage of ionizing radiation into the intergalactic medium. Using JWST/NIRSpec spectroscopy from the CAPERS, CEERS, and RUBIES surveys, we compile 160 NIRCam-selected EELGs in the EGS field. These galaxies show extreme rest-frame equivalent widths (EWs), with a median EW([O III]+H${\beta}$)=1616Å and EW(H${\alpha}$)=763Å. They are low-mass (median log(M$_{\star}$/M$_{\odot}$)=8.26) with high sSFRs (median 43 Gyr$^{-1}$), above the $z\sim6$ main sequence. UV slopes are diverse, with a mean $\beta=-2.0$, and only 7% have extremely blue continua ($\beta<-2.6$). Emission-line diagnostics suggest stellar populations as the primary ionizing source, although an AGN fraction of 14% is found. These galaxies are efficient ionizing photon producers, with mean log($\xi_{\rm ion}$ [Hz erg$^{-1}$])=25.34, exceeding typical values at similar redshifts. Escape fractions, however, are heterogeneous: 9% of EELGs show escape fractions $>$10% for both Ly${\alpha}$ and LyC photons, while 82% lack detectable Ly${\alpha}$ emission. The median inferred LyC escape fraction is modest (4.4%) but enhanced in super-Eddington systems with sSFR >25 Gyr$^{-1}$. The galaxies are extremely compact, with a median effective radius of 0.49 kpc, and exhibit a recent star-formation burst. Our analysis indicates that sSFR and star-formation rate surface density are the primary drivers of their extreme emission line strengths.

Ultracompact Galactic binary systems (UCBs) emit low-frequency gravitational waves (GWs). The emission of GWs is causing these systems to lose angular momentum, which is detectable by observing an orbital period decay. ZTFJ213056.71+442046.5 (ZTFJ2130) is an UCB with a period of 39.3401(1) minutes consisting of a Roche lobe-filling hot subdwarf and a white dwarf companion. We attempt to measure the orbital decay rate $\dot{P}$ caused by GW emission of ZTFJ2130 and predict the expected GW signal for LISA. High-speed photometry was conducted using the FLI Kepler KL4040FI CMOS camera, mounted to the 1.2-meter Oskar Lühning telescope at the Hamburg Observatory as well as the Hamamatsu ORCA-Quest 2 qCMOS camera at the 1.23-meter telescope at CAHA in Spain. ZTFJ2130 was observed on six nights between August 2024 and September 2025. The obtained lightcurves combined with previous high-cadence observations were used to conduct an O-C timing analysis. Additionally, we employed the LISA data analysis tool ldasoft to model the expected GW data. We measure a period change of $(-1.97\pm0.05)\times10^{-12}\,\mathrm{ss^{-1}}$. Assuming only GW emission, this result was used to calculate a chirp mass of $(0.408\pm0.006)\,\mathrm{M_{\odot}}$. From ldasoft we predict that LISA will be able to measure the chirp mass with an uncertainty of 5%. We measure $\dot{P}$ with an uncertainty of only 2% and show that modern (q)CMOS detectors are well suited to provide precise timing measurements, enabling the measurement of the orbital decay of UCBs with high precision with modest size telescopes. The derived orbital decay is fully consistent with predictions from spectral and lightcurve modeling. We show that future observations with LISA can potentially provide a deviation from only gravitational wave effects, e.g. due to accretion, if the effect is sufficiently large.

The expansion history of the Universe prior to Big Bang Nucleosynthesis (BBN) remains largely unconstrained. The high-energy post-inflationary era may involve multiple distinct epochs, each characterized by a different equation of state (EoS). A key prediction of inflation is the generation of tensor perturbations that later manifest as a stochastic background of primordial gravitational waves (GWs). The large-scale amplitude and small-scale spectral tilt ($n_{\rm GW}$) of these GWs encode the inflationary energy scale and the subsequent expansion history, respectively. A soft post-inflationary EoS ($w<1/3$) yields red-tilted GW spectra ($n_{\rm GW}<0$), while a stiff EoS ($w>1/3$) results in a blue-tilt ($n_{\rm GW}>0$). In our previous work [arXiv:2407.07956], we developed an analytical framework for computing the GW spectral energy density, $\Omega_{\rm GW}(f)$, for multiple post-inflationary transitions ($w_1 \to w_2 \to \cdots \to w_n \to 1/3$), focusing on the parameter space relevant for future GW observations. In this paper, we extend that framework to systematically investigate the $morphological~diversity$ of inflationary GW spectra generated by multi-epoch post-inflationary histories. Remaining model agnostic, we demonstrate that a wide variety of spectral shapes, ranging from convex and concave monotonic profiles to multi-peaked non-monotonic spectra, can naturally emerge depending on the sequence and duration of these epochs. We also introduce GWInSpect, a publicly available Python package that computes $\Omega_{\rm GW}(f)$ for arbitrary sequences of EoS transitions, providing a practical tool to study the pre-BBN expansion history of the Universe.

Powerful, large-scale relativistic jets are usually associated with massive, old elliptical galaxies. This paradigm has recently been challenged by the identification of narrow-line Seyfert 1 (NLSy1) galaxies, thought to be young active galactic nuclei with low-mass black holes, harboring relativistic jets. Among them, sources hosting $>$100 kpc radio jets are extremely rare. Here, we report the discovery of large-scale, double-lobed radio structures in 33 NLSy1s with the projected linear size of at least 100 kpc from a recently published catalog of 22656 NLSy1 galaxies. These 33 include 29 confirmed double-lobed sources and 4 candidates whose radio structure requires further study. We suggest that their low black hole masses are unlikely to be due to their small angles of inclination to the line of sight. These enigmatic sources were identified by examining the radio observations taken with the Faint Images of the Radio Sky at Twenty centimeters, Very Large Array Sky Survey, Low Frequency Array, and Rapid ASKAP Continuum Survey. Among them, the NLSy1 source J1318+2626 ($z=0.62$) is found to host a radio jet with the projected linear size of 4.3 Mpc, making it the only NLSy1 galaxy hosting a Mpc-scale radio jet known as of today. We conclude that future observations of NLSy1 sources with the next generation of sensitive telescopes may reveal a much larger population harboring large-scale jets, thus providing crucial clues on their origin, propagation, and interaction with the ambient environment.

A recent asteroseismic analysis suggests that Kepler-56 -- a planet-hosting red giant -- exhibits a unique spin structure: (1) the spin axes of the core and envelope are misaligned; and (2) the envelope rotates approximately an order of magnitude faster than typical red giants. In this paper, we investigate a feasible scenario to reproduce this spin structure by estimating the amount of the angular momentum (AM) supply from the planets through the simplified calculation of the time evolution of AM. As a result, unless the tidal efficiency is extremely high, we show that the tidal interactions between the known close-in planets (Kepler-56 b and c) are insufficient to supply the AM required to accelerate Kepler-56 from the spin rate observed in typical red giants. We also show that the engulfment of a hot Jupiter can be expected to provide sufficient AM supply for the acceleration ant that the mass and orbit of the engulfed hot Jupiter are constrained by a mass of 0.5-2 Jupiter masses and an orbital period of 1-6 days. On the other hand, if Kepler 56 was already rapidly spinning before entering the RG stage and requires no acceleration, the obliquity damping by the known close-in planets can reproduce the spin structure of Kepler-56. Even in such cases, planetary engulfment during the MS stage might be involved in achieving rapid spin before the tidal alignment. These discussions demonstrate the importance of Kepler-56 as a candidate for planetary engulfment that may leave traces of its spin structure.

We develop a second-order cosmological perturbation theory on a background geometry expressed in terms of light-cone coordinates, extending the first-order analyses available in the literature. In particular, we investigate the gauge transformations of second-order perturbative quantities on the light-cone and establish their connection with standard perturbation theory. Through a consistent matching procedure, we identify the second-order gauge fixing that corresponds to the non-linear Geodesic Light-Cone gauge within standard perturbation theory, known as the Observational Synchronous Gauge. We then emphasize its conceptual similarities and differences wrt the standard Synchronous Gauge. Finally, within this new perturbative framework, and adopting a fully gauge-invariant approach, we compute the luminosity distance-redshift relation up to second order with anisotropic stress as seen by a free-falling observer. Remarkably, we show how divergences at the observer position can be eliminated in a completely model independent way. These results validate our perturbative framework and establish it as a novel formalism for evaluating cosmological observables at second order.

We present a simulation-based forward-modeling framework for cosmological inference from optical galaxy-cluster samples, and apply it to the abundance and weak-lensing signals of DES-Y1 redMaPPer clusters. The model embeds cosmology-dependent optical selection using a counts-in-cylinders approach, while also accounting for cluster miscentering and baryonic feedback in lensing. Applied to DES-Y1, and assuming a flat $\Lambda$CDM cosmology, we obtain $\Omega_m=0.254^{+0.026}_{-0.020}$ and $\sigma_8=0.826^{+0.030}_{-0.034}$, consistent with a broad suite of low-redshift structure measurements, including recent full-shape analyses, the DES/KiDS/HSC 3$\times$2 results, and most cluster-abundance studies. Our results are also consistent with \textit{Planck}, with the difference being significant at $2.58\sigma$. These results establish simulation-based forward-modeling of cluster abundances as a promising new tool for precision cosmology with Stage~IV survey data.

We present analysis of JWST NIRCam and NIRSpec observations of the galaxy SXDF-NB1006-2 at $z = 7.212$, as part of the Reionization and the ISM/Stellar Origins with JWST and ALMA (RIOJA) project. We derive the physical properties by conducting spectral energy distribution (SED) fitting, revealing that our target is a young (age $\sim2$ Myr) starburst galaxy with intense radiation field. We detect multiple nebular emission lines from NIRSpec IFS data. We identify a robust broad component of [O III]$\lambda5008$ emission, indicating the presence of ionized gas outflows. The derived gas depletion time of a few hundred Myr implies that our target could be one of the progenitors of massive quiescent galaxies at $z\sim4-5$ identified by recent JWST observations. The spatial distribution of optical and far-infrared (FIR) [O III] emission lines differs in morphology, likely resulting from different critical densities and inhomogeneous density distributions within the galaxy. Potential old stellar populations may be necessary to account for the derived metallicity of $\sim0.2\,\rm{Z}_\odot$, and their presence can be confirmed by future MIRI observations. Including our target, star-forming galaxies at $z>6$ detected by ALMA are generally very young but more massive and brighter in UV than galaxies identified by only JWST. The ALMA-detected galaxies may also have a steeper mass-metallicity relation. These findings suggest that the ALMA-detected galaxies may have experienced more efficient mass assembly processes in their evolutionary pathways.

The radiation-dominated polar emitting region of an accreting X-ray pulsar is simulated numerically in the framework of a three-dimensional (geometrically two-dimensional) model. The radiative transfer within the emitting region and the structure of the latter are calculated with the use of the self-consistent algorithm developed earlier. The magnetic scattering cross sections dependent on the photon energy and polarization have been incorporated. Second-order bulk Comptonization over entire emitting region, induced Compton scattering, the switching of the polarization modes, free-free processes, the cyclotron emission because of electron-proton collisions, and a realistic shape of the accretion channel have been taken into account. The case of a dipole magnetic field is considered. It is shown that the induced Compton effect can play a notable role in establishing the electron temperature in the post-shock zone. Within the model shock wave, a higher electron temperature is achieved than in the post-shock zone by means of the bulk-heating mechanism. The photons, gained the energy in the shock wave and above it due to bulk motion effects and the thermal Doppler effect, are responsible for the formation of high-energy regions in the emergent continuum of the polarization modes.

We present an investigation of quantum-corrected black hole spacetimes coupled with clouds of strings, examining two distinct theoretical models that incorporate quantum gravitational effects through different implementations of correction terms. Our study explores the geodesic structure, focusing on photon sphere properties, black hole shadows, and innermost stable circular orbits of test particles around these exotic geometries. The analysis reveals fundamental modifications to particle trajectories that distinguish quantum-corrected solutions from their classical counterparts, with observable implications for high-energy astrophysics. We investigate quasi-periodic oscillations arising from test particle motion, deriving frequency relationships that could serve as observational probes of quantum gravity effects in accreting black hole systems. Through rigorous gravitational lensing analysis using the Gauss-Bonnet theorem, we calculate weak-field deflection angles and identify distinctive signatures that enable discrimination between the two quantum correction models. The gravitational lensing study reveals opposite dependencies on quantum parameters between the models, providing unambiguous observational discriminants. Additionally, we analyze the thermodynamic properties including temperature, entropy, and heat capacity modifications, exploring topological characteristics and phase transition behavior in these quantum-corrected systems. The investigation reveals that the two models exhibit contrasting behaviors across multiple observational channels, from gravitational lensing deflection angles to quasi-periodic oscillation frequencies, providing remarkable discriminants for testing quantum gravity theories.

The Australasian strewn field covers more than 15% of Earth's surface, consisting of tektites and microtektites. Australasian tektites from Southeast Asia and Australia, as well as microtektites recovered from deep sea sediments and Antarctica, are established to be derived from upper continental crust sediments. An expedition to retrieve remnants of bolide CNEOS 2014 January 8 (IM1), held in the Pacific Ocean, was in proximity to the known extent of the Australasian strewn field, and yielded "BeLaU"-spherules, whose compositions did not match most well-studied solar system material. We therefore report precise and comprehensive elemental data for Australasian tektites to compare their elemental abundances to those of microtektites from deep sea sediments, and BeLaU. Our findings corroborate previous studies that Australasian tektites and microtektites closely resemble the elemental abundance patterns of the upper continental crust. Meanwhile, the elemental patterns of the BeLaU-spherules are distinct from the Australasian tektite/microtektite compositions.

Detection of sub-Hz gravitational waves is of significant importance for astrophysics. It enables the observation of intermediate-mass black hole mergers, the issuance of early alerts for gravitational-wave events, and the exploration of the stochastic gravitational-wave background. The Cryogenic sub-Hz cROss torsion-bar detector with quantum NOn-demolition Speed meter (CHRONOS) is a proposed gravitational-wave detector based on a Sagnac speed-meter topology that uses torsion bars as test masses. Its prototype design aims to achieve a strain sensitivity of $3 \times 10^{-18}~\mathrm{Hz}^{-1/2}$ at 1~Hz and thus enable the detection of $\mathcal{O}(10^4),M_\odot$ intermediate-mass black hole mergers at 100~Mpc with a signal-to-noise ratio of 3. We show that the torsion-bar-based speed meter can suppress noise originating from laser intensity fluctuations by canceling the net torque on the bar and by using a balanced homodyne readout. We further present, for the first time, an analytic intensity-noise model for a gravitational-wave detector employing a torsion-bar Sagnac speed-meter configuration. Using this model, we evaluate the expected performance of a 2.5~m arm-length CHRONOS prototype. The projected laser-intensity noise is $2.9 \times 10^{-20}~\mathrm{Hz}^{-1/2}$ at 1~Hz, which is sufficiently low to allow the detection of binary intermediate-mass black hole mergers.

We present the optical design and sensitivity modeling of the 2.5 m Cryogenic sub-Hz cROss torsion-bar detector with quantum NOn-demolition Speed meter (CHRONOS), a triangular Sagnac speed-meter interferometer incorporating power- and signal-recycling techniques. Using ABCD-matrix analysis and Finesse3 simulations, we show that stable eigenmodes are obtained with optimized mirror curvatures and focal placements, achieving mode-matching efficiencies above 99.5 %. The resulting configuration reaches a quantum-noise-limited strain sensitivity of $h \simeq 3\times10^{-18},\mathrm{Hz^{-1/2}}$ at 1 Hz, with a ring-cavity finesse $\mathcal{F}\simeq3.1\times10^{4}$ and round-trip Gouy phase $\psi_{\mathrm{rt}}\approx153^{\circ}$. The power-recycling cavity detuning ($\phi_p=-85^{\circ}$) dominates the low-frequency quantum noise, while the signal-recycling cavity detuning ($\phi_s=0^{\circ}$) mainly introduces a uniform quadrature rotation. The optimal homodyne angle ($\zeta_{\mathrm{opt}}\simeq46^{\circ}$) balances shot-noise and radiation-pressure effects to give the best sensitivity near 1 Hz. Assuming end-mirror reflectivity $R_{\mathrm{ETM}}=99.9999\%$ under cryogenic operation at 10 K, CHRONOS can achieve quantum-noise-limited performance on a laboratory scale and serve as a testbed for future long-baseline, cryogenic interferometers probing sub-hertz gravitational waves.

While the 14N(p,g)15O reaction plays a key role in the hydrogen-burning processes in various stellar conditions, its reaction rate is not known with sufficient precision. Therefore, the first scientific project at the recently launched Bellotti Ion Beam Facility of the Laboratori Nazionali del Gran Sasso was the measurement of the 14N(p,g)15O reaction cross section in the proton energy range between 250 and 1500 keV. In this paper, the experimental techniques are summarized with special emphasis on the description of solid state nitrogen target production and characterization. The first results of the reaction yield measured at 55 deg detection angle are also presented.

The chart of nuclei could be enlarged with a branch describing neutron stars that are huge nuclei of a few solar masses held together by gravity force and sustained by the pressure due to the degenerate Fermi sea. We contend in this manuscript that yet another branch could be added: objects with a large weak charge, with masses around $10^{-3}$ solar masses and having radii of a few meters, very compact, only slightly larger than their Schwarzchild radius, and sustained by the pressure generated by the weak force due to $Z$ exchange. This interaction, insignificant in normal neutron stars, could become dominant when ultrahigh densities are reached due to the action of gravity and lead to stable configurations if the appropriate conditions are met. They would constitute a physical realization of the equation of state proposed by Zeldovich some decades ago.

We study hierarchical triple systems formed by a compact binary orbiting a supermassive black hole (SMBH), focusing on the role of relativistic magnetic tidal interactions. Extending previous analyses of precession resonances to 0.5 post-Newtonian order, we incorporate quadrupolar magnetic tidal moments, which have no Newtonian counterpart. Averaging the Hamiltonian over the inner orbit, we find that magnetic tides introduce new resonances absent at lower order, leading to additional eccentricity excitations and significantly modifying the binary's long-term evolution. Numerical solutions of the Lagrange Planetary Equations confirm these analytical predictions and reveal how resonance strength depends on orbital eccentricity and inclination. The resulting dynamics accelerates the binary merger and imprint distinctive signatures on gravitational waves, potentially observable by LISA. Our findings identify magnetic tidal coupling as a novel strong-gravity effect and establish its importance for the resonant dynamics of compact-object binaries near SMBHs.

While photonic lanterns efficiently and uniquely map a set of input modes to single-mode outputs (or vice versa), the optical mode transfer matrix of any particular fabricated device cannot be constrained at the design stage due to manufacturing imperfections. Accurate knowledge of the mapping enables complex sensing or beam control applications that leverage multimode conversion. In this work, we present a characterisation system to directly measure the electric field from a photonic lantern using digital off-axis holography, following its evolution over a 73 nm range near 1550 nm and in two orthogonal, linear polarisations. We provide the first multi-wavelength, polarisation decomposed characterisation of the principal modes of a photonic lantern. Performance of our testbed is validated on a single-mode fibre then harnessed to characterise a 19-port, multicore fibre fed photonic lantern. We uncover the typical wavelength scale at which the modal mapping evolves and measure the relative dispersion in the device, finding significant differences with idealised simulations. In addition to detailing the system, we also share the empirical mode transfer matrices, enabling future work in astrophotonic design, computational imaging, device fabrication feedback loops and beam shaping.

The stochastic gravitational-wave background (SGWB) from eccentric stellar-mass binary black holes (SBBHs) holds crucial clues to their origins. For the first time, we employ a Bayesian framework to assess the detectability and distinguishing features of such an SGWB with spaceborne detectors, while accounting for contamination from the Galactic foreground. Our analysis covers eccentric SBBHs from three formation channels: isolated binary evolution, dynamical assembly in globular clusters (GCs), and in active galactic nuclei (AGNs). We find that TianQin, LISA, and Taiji can detect the SGWBs from both isolated and GC-formed SBBHs after 4 years of operation, with the corresponding signal-to-noise ratios of around 10, 60, and 170. However, these backgrounds are spectrally degenerate with a strictly power-law SGWB. Furthermore, highly eccentric SBBHs formed in AGNs yield an SGWB marked by a spectral turnover and sharp decline. While this feature lowers the signal-to-noise ratio by approximately an order of magnitude, it can enable a clear distinction from the strictly power-law background using LISA and Taiji.

We investigate how matter density distributions affect thin-wall bubble formation in the asymmetron mechanism, a scalar-tensor theory with a universal coupling to matter and explicit symmetry-breaking, and analyse the stability of its metastable state. We show that the screening mechanism of the asymmetron inside dense objects induces a surface tension associated with the boundary of the screening object, leading to a richer class of bubble solutions than the standard Coleman-Callan bulk nucleation. These boundary surface tensions are used to modify the Nambu-Goto action for instantons, allowing for the computation of the corresponding Euclidean action for bubbles nucleating on flat planes, as well as on concave and convex cylindrical surfaces. We find that the smallest Euclidean action occurs for bubbles nucleating along the edge of a concave spherical surface. Comparing this edge nucleation channel with the bulk one, we determine the maximum curvature radius for which concave edge nucleation is preferred. Since the maximum radius of curvature is exponentially suppressed by the action of a bulk bubble, we find that within the regime of the instanton approximation, edge nucleation is always preferred. This is largely due to the weak couplings of the asymmetron. We apply these findings to determine the maximum curvature radius of a cosmic void and discuss how our results affect the seeding of $N$-body simulations of asymmetron domains, showing that domain wall nucleation preferentially occurs at the edges of cosmological voids. We also demonstrate that the presence of a homogeneous gas around the dense substrates reduces the maximum curvature radius, enabling bulk bubbles to form preferentially as the asymmetron undergoes a density-driven phase transition.

We revisited the scenario of a resonant enhancement in the oscillation probability due to an interaction potential between neutrinos and dark matter with the novelty of the inclusion of the generalized uncertainty principle. It is shown that a new resonant conversion appears at higher energies. This effect could be tested with future neutrino data as new dips in the ultrahigh energy neutrino flux.

We investigate the cross-correlation between astrometric and timing-residual observables for distant sources, such as pulsars and galaxies, and equivalent observables for nearby solar system bodies. Using the unified spin-weighted formalism introduced in [1], we derive the angular correlation functions-generalised Hellings-Downs curves-that describe the response of these mixed observables to a stochastic, unpolarised gravitational-wave background (SGWB). We compute the expected signal-to-noise ratio (SNR) and sensitivity for such measurements, focusing on cross-correlations between pulsar timing array (PTA) redshift signals and astrometric or distortion (shimmering) effects induced in solar system objects such as asteroids. Although the current astrometric precision of asteroid tracking does not yet provide competitive constraints relative to PTA-only surveys, the method offers a complementary probe with enhanced sensitivity at higher frequencies. Future wide-field surveys capable of sub-milliarcsecond precision could make this approach a viable tool for detecting or constraining the SGWB. A key advantage of the technique is its reduced susceptibility to correlated systematics across different measurement domains, providing an independent cross-check of PTA detections and a potential observational bridge between PTA and LISA frequency bands.