Papers for Wednesday, Aug 27 2025

With dayside temperatures elevated enough for all atmospheric constituents to be present in gas form, ultra-hot Jupiters offer a unique opportunity to probe the composition of giant planets. We aim to infer the composition and thermal structure of the dayside atmosphere of the ultra-hot Jupiter WASP-121b from two NIRISS$/$SOSS secondary eclipses observed as part of a full phase curve. We extract the eclipse spectrum of WASP-121b with two independent data reduction pipelines and analyse it using different atmospheric retrieval prescriptions to explore the effects of thermal dissociation, reflected light, and titanium condensation on the inferred atmospheric properties. We find that the observed dayside spectrum of WASP-121b is best fit by atmosphere models possessing a stratospheric inversion with temperatures reaching over 3000K, with spectral contributions from H2O, CO, VO, H-, and either TiO or reflected light. We measure the atmosphere of WASP-121b to be metal enriched (~10x stellar) but comparatively titanium poor (~1x stellar), potentially due to partial cold-trapping. The inferred C/O depends on model assumptions such as whether reflected light is included, ranging from being consistent with stellar if a geometric albedo of zero is assumed to being super-stellar for a freely fitted Ag = 0.16 +/- 0.02. The volatile-to-refractory ratio is measured to be consistent with the stellar value. We infer that WASP-121b has an atmosphere enriched in both volatile and refractory metals, but not in ultra-refractory titanium, suggesting the presence of a nightside cold-trap. Considering H2O dissociation is critical in free retrieval analyses, leading to order-of-magnitude differences in retrieved abundances for WASP-121b if neglected. Simple chemical equilibrium retrievals assuming that all species are governed by a single metallicity parameter drastically overpredict the TiO abundance.

KMT-2018-BLG-0029Lb and OGLE-2019-BLG-0960Lb were the lowest mass-ratio microlensing planets at the time of discovery. For both events, microlensing parallax measurements from the Spitzer Space Telescope implied lens systems that were more distant and massive than those inferred from the ground-based parallax. Here, we report on the detection of excess flux aligned to the event locations using Keck Adaptive Optics imaging, which is consistent with the expected brightness of main-sequence hosts under the ground-based parallax, but inconsistent with that predicted by Spitzer. Based on the excess flux, ground-based parallax, and angular Einstein radius, we determine KMT-2018-BLG-0029Lb to be a $4.2\pm0.5 M_\oplus$ planet orbiting a $0.70\pm0.07 M_\odot$ host at a projected separation of $3.1\pm0.3$ au, and OGLE-2019-BLG-0960Lb to be a $2.0\pm0.2 M_\oplus$ planet orbiting a $0.40\pm0.03 M_\odot$ host at a projected separation of $1.7\pm0.1$ au. We report on additional light-curve models for KMT-2018-BLG-0029 under the generalized inner-outer (offset) degeneracy, which were not reported in the original analysis. We point out inconsistencies in the inner/outer labeling of the degenerate models in the lens and source planes, and advocate for the lens-plane convention, which refers to the planet being closer or further to the host star compared to the image it perturbs. Lastly, we discuss the possibility of breaking this degeneracy via ground concurrent observations with the Roman Space Telescope.

Blazars are promising high-energy neutrino source candidates. Yet, leptohadronic models face challenges in describing neutrino emission within a viable energy budget, and their predictive power is limited by the commonly used single-zone approximation and the reliance on phenomenological parameters. In this work, we present a new leptohadronic model where a sub-Eddington jet evolves from magnetically- to kinetically dominated. A small fraction of the electrons and protons picked up by the jet are continuously accelerated to a power-law spectrum, estimated based on the local magnetic field strength, turbulence, and ambient density, for which we assume power-law profiles. The model parameters are thus directly tied to the jet physics and are comparable in number to typical single-zone models. We then numerically calculate the emission along the jet. Applying the model to the IceCube candidate TXS 0506+056, we find that protons are accelerated to EeV energies in the inner jet, producing a neutrino flux up to order 100 PeV that is consistent with the public 10 year IceCube point-source data. Proton emission at 0.1 pc describes the X-ray and gamma-ray data, while electron emission at the parsec scale describes the optical data. Protons carry a power of about 1% of the Eddington luminosity, showing that the model is energetically viable. The particle spectra follow $E^{-1.8}$, with diffusion scaling as $E^{0.3}$, ruling out Bohm-like diffusion. Additional particle injection near the broad line region can reproduce the 2017 flare associated to a high-energy neutrino. We also apply the model to blazar PKS 0605-085, which may be associated with a recent neutrino detected by KM3NeT above 100 PeV. The results suggest that blazars are efficient neutrino emitters at ultra-high energies, making them prime candidates for future experiments targeting this challenging energy range.

The faint radio-source population includes sources dominated both by star formation and active galactic nuclei (AGN), encoding the evolution of activity in the Universe. To investigate its nature, we probabilistically classified 4,471 radio sources at z < 0.947 using low-frequency radio data from the LoTSS Deep Fields alongside a multi-component model for nebular emission, sampled by spectra obtained with the Dark Energy Spectroscopic Instrument (DESI). This was done by combining three tools: (i) the identification of a radio excess, (ii) the BPT diagram, and (iii) a modified Mass Excitation diagram, alongside Monte Carlo methods to estimate the probability that each source is either a star-forming galaxy (SFG), a radio-quiet AGN (RQ AGN), or a high-\low-excitation radio galaxy (HERG or LERG). This approach extends the probabilistic classification framework of previous works by nearly doubling the redshift range, such that we can now probabilistically classify sources over the latter half of cosmic history. Often regarded as the 'gold standard' method, spectroscopic classifications allow us to evaluate the performance of other methods. Using a 90 per cent reliability threshold, we find reasonable overall agreement (~77 per cent) with state-of-the-art photometric classifications, but significant differences remain, including that we identify 2-5 times more RQ AGN. Furthermore, our high-confidence spectroscopic classifications show that radiatively-efficient and inefficient AGN exhibit clearly distinct Eddington-scaled accretion rate distributions, contrary to recent findings in the literature. Overall, our results highlight the need for new and forthcoming spectroscopic campaigns targeting radio sources, on the pathway to the SKA.

We present the discovery of deep, irregular, periodic transits towards the white dwarf ZTF$\,$J1944$+$4557 using follow-up time-series photometry and spectroscopy from Palomar, Keck, McDonald, Perkins, and Lowell observatories. We find a predominant period of 4.9704$\,$hr, consistent with an orbit near the Roche limit of the white dwarf, with individual dips over 30$\%$ deep and lasting between 15 and 40 minutes. Similar to the first known white dwarf with transiting debris, WD$\,$1145$+$017, the transit events are well-defined with prominent out-of-transit phases where the white dwarf appears unobscured. Spectroscopy concurrent with transit photometry reveals the average Ca$\,$K equivalent width remains constant in and out of transit. The broadening observed in several absorption features cannot be reproduced by synthetic photospheric models, suggesting the presence of circumstellar gas. Simultaneous $g+r$- and $g+i$-band light curves from the CHIMERA instrument reveal no color dependence to the transit depths, requiring transiting dust grains to have sizes $s \gtrsim0.2\,\mu$m. The transit morphologies appear to be constantly changing at a rate faster than the orbital period. Overall transit activity varies in the system, with transit features completely disappearing during the seven months between our 2023 and 2024 observing seasons and then reappearing in 2025$~$March, still repeating at 4.9704$\,$hr. Our observations of the complete cessation and resumption of transit activity provide a novel laboratory for constraining the evolution of disrupted debris and processes like disk exhaustion and replenishment timescales at white dwarfs.

We investigate the impact of radiation pressure on the circumbinary discs surrounding accreting massive black hole binaries (MBHBs) at milli-parsec separations, using 3D hyper-Lagrangian resolution hydrodynamic simulations. The circumbinary discs in our simulations evolve under an adiabatic equation of state. The gas temperature is therefore allowed to change through viscous heating, black-body cooling and self-gravity. We take a significant step further by including the contribution of radiation pressure in the simulations. We model binaries with a total mass of $10^6 \, M_{\odot}$, eccentricities $e=0,0.45,0.9$ and mass ratios $q= 0.7, 1$. We find that the radiation pressure significantly alters the vertical and thermal structure of the disc, resulting in a geometrically thinner, therefore colder configuration. This leads to a reduced accretion rate onto the binary and suppresses cavity eccentricity growth and precession in circular equal mass binaries. The binary eccentricity remains approximately constant, while the semi-major axis decreases over time due to net negative torque, regardless of the initial binary orbital parameters.

We analyze high-contrast, medium-spectral-resolution $H_{\rm \alpha}$ observations of the star AB Aurigae using the Very Large Telescope's Multi Unit Spectroscopic Explorer (MUSE). In multiple epochs, MUSE detects the AB Aur b protoplanet discovered from Subaru/SCExAO data in emission at wavelengths slightly blue-shifted from the $H_{\rm \alpha}$ line center (i.e. at 6558.88--6560.13 Å; $\sim$ -100 km s$^{-1}$) and in absorption at redshifted wavelengths (6562.8--6565.1 Å; $\sim$ 75 km s$^{-1}$). AB Aur b's $H_{\rm \alpha}$ spectrum is inconsistent with that of the host star or the average residual disk spectrum and is dissimilar to that of PDS 70 b and c. Instead, the spectrum's shape resembles that of an inverse P Cygni profile seen in some accreting T Tauri stars and interpreted as evidence of infalling cold gas from accretion, although we cannot formally rule out all other nonaccretion origins for AB Aur b's MUSE detection. AB Aurigae hosts only the second protoplanetary system detected in $H_{\rm \alpha}$ thus far and the first with a source showing a spectrum resembling an inverse P Cygni profile. Future modeling and new optical data will be needed to assess how much of AB Aur b's emission source(s) originates from protoplanet accretion reprocessed by the disk, a localized scattered-light feature with a unique $H_{\rm \alpha}$ profile, or another mechanism.

We investigate the possible presence of systematic rotation in the member galaxies of a sample of 17 nearby ($z<0.1$), rich (at least 80 identified members) Abell clusters. We also assess the extent to which low-number statistics may influence the recovery of the rotation parameters. Following the methods often used in the context of globular clusters and of clusters of galaxies, we estimate a representative value of the systematic rotation velocity and the position angle of the projected rotation axis for the set of spectroscopically confirmed member galaxies within 1.5 Mpc from the centre of each cluster. We study the robustness of our rotational velocity measurements as a function of the number of galaxies included in the analysis with a bootstrapping technique. Eight clusters with sufficiently abundant and regular data (A1367, A1650, A2029, A2065, A2142, A2199, A2255 and A2670) exhibit a significantly high rotational velocity, when compared to their velocity dispersion ($v_{rot}/\sigma\geq 0.15$). Interestingly, three of them (A1650, A2029 and A2199) are confirmed to be cool-core, relaxed clusters with no evidence of recent mergers, as suggested by X-ray observational data. We also find a general tendency to overestimate the value of $v_{rot}$ when the number of galaxies with measured velocities is reduced, for which we put forward an analytical justification.

Tidal interactions are one of the primary drivers of orbital evolution for massive planets with short orbital periods. Tidal dissipation within host stars can cause the orbits of such planets to decay. However, the mechanisms of tidal dissipation are difficult to probe. Generally, tidal dissipation is parameterized by the modified stellar tidal quality factor, or $Q_{*}^{'}$, but the lack of observational evidence of orbital decay to confirm dissipation theories has resulted in orders of magnitude of uncertainty in $Q_{*}^{'}$. We present a new transit timing analysis of 54 systems with varying stellar evolutionary states in an attempt to search for orbital decay across multiple stages of stellar evolution. For each system, we obtained mid-transit times from new TESS data and evaluated potential departures from a linear ephemeris using the Bayesian Information Criterion. We then determined tidal quality factors using widely tested theoretical relations. Of the systems studied, 25 showed evidence of a decrease in orbital period over time and 8 showed evidence of a decrease that was inconsistent with 0 by three standard deviations: CoRoT-2, TrES-5, WASP-4, WASP-12, WASP-19, WASP-45, WASP-99, and XO-3. However, the significance of some of these detections may be influenced by unreliable transit time measurements. Similarly, we see that the lower limit on $Q_{*}^{'}$ for evolved systems is marginally lower than it is for the sample of all systems, though both limits are orders of magnitude below the expected theoretical values for both samples. These new constraints on transit times and $Q_{*}^{'}$ values will help to narrow the search for orbital decay in the future and place important constraints on current theories of star-planet interactions.

The AT 2018cow-like fast blue optical transient AT2022tsd showed a large number of few-minute-duration, high-luminosity (~10^43 erg/s) flares. We present an intensive search for such flares from another 18cow-like event, AT2024wpp. We have used the Large Array Survey Telescope (LAST) to observe this transient between 28 and 74 days after the approximate time of zero flux. The target was observed for about 23 hours to a sensitivity that allows one to detect 3x10^42 erg/s flares at S/N>5. No optical flares have been found, suggesting a one-sided 2-sigma confidence upper limit of <0.02 on the flare's duty cycle, and flare rate lower than about 0.11/hr. These limits suggest that not all 18cow-like objects display a high rate of minute-timescale luminous flares. This can be explained either by diversity in the 18cow-like population or by viewing angle effects (e.g., beaming), or rather that the optical depth towards the central emitting region did not fall below unity during the particular search window.

We present a general-relativistic study of the distribution of proton superconductivity in strongly magnetized neutron stars (NSs), employing the $XNS$ code to solve the coupled Einstein-Maxwell equations. This work extends the previous investigation by us in a few key directions: (i) we explore equilibrium configurations with both toroidal and poloidal magnetic field geometries, (ii) we account for complex many-body effects incorporating microscopically derived pairing gaps, and (iii) our models utilize equation(s) of state (EoS) derived from microscopic many-body theory, including realistic two- and three-body nuclear forces, as well as from relativistic mean-field models. We compare superconducting topologies across our collection of EoS and explore the influences of magnetic field geometry in stellar models parameterized by central density. Our models confirm the absence of $S$-wave superconductivity in the inner core and, importantly, reveal that non-superconducting regions exhibit complex three-dimensional geometries: doughnut-shaped for toroidal fields and prolate-shaped for poloidal fields -- spatial structures that are inherently absent from one-dimensional analyses. We also compute magnetic deformations and ellipticities for several millisecond pulsars (MSPs), estimating their continuous gravitational wave (CGW) strain. While these MSPs remain undetectable by current detectors, next-generation instruments such as the Einstein Telescope and Cosmic Explorer may detect their signals, opening an observational window into internal superconductivity and internal magnetic field of NSs, as well as, the fundamental microphysics of dense matter.

The characteristic rapid rise and decline at optical wavelengths of a kilonova is the product of the low ejecta mass ($\lesssim 0.05 M_\odot$) and high ejecta velocity ($\gtrsim 0.1$c). We show that, even at very early times ($\lesssim 2$ days), regions of ejecta fall below critical density and temperature thresholds at which non-local thermodynamic equilibrium (NLTE) effects become important. Here, we present an approximate method for calculating the ionization state of the ejecta that accounts for the NLTE impact of high-energy electrons produced in the beta decay of freshly synthesized $r$-process elements. We find that incorporating ionization from high-energy electrons produces an ``inverted" and ``blended" ionization structure, where the most highly ionized species are located in the fastest moving homologous ejecta and multiple ionization states coexist. In radiation transport calculations, the higher degree of ionization reduces line blanketing in optical bands, leading to improved agreement with the light curve properties of AT\,2017gfo such as the duration, decay rates, brightness, and colors. Our quasi-NLTE implementation helps to alleviate tensions in kilonova modeling: for high-velocity ($\sim 0.3c$) ejecta components our models require less mass for a given peak brightness in optical bands, by as much as a factor of 3; our models can explain the presence of observed features associated to Sr II, W III, Se III, and Te III under conditions where LTE models would predict only neutral species; and we naturally predict the coexistence of species like Sr II and Ce III without the need for fine-tuning of the ejecta properties.

Few spectra of directly-imaged exoplanets have been obtained in the mid-infrared (> 3 $\mu$m). This region is particularly rich in molecular spectral signatures, whose measurements can help recover atmospheric parameters and provide a better understanding of giant planet formation and atmospheric dynamics. In the past years, exoplanet interferometry with the VLTI/GRAVITY instrument has provided medium-resolution spectra of a dozen substellar companions in the near infrared. The 100-meter interferometric baselines allow for the stellar and planetary signals to be efficiently disentangled at close angular separations (< 0.3''). We aim to extend this technique to the mid-infrared using MATISSE, the VLTI's mid-infrared spectro-interferometer. We take advantage of the fringe tracking and off-axis pointing capabilities recently brought by the GRA4MAT upgrade. Using this new mode, we observed the giant planet $\beta$ Pictoris b in L and M bands (2.75-5 $\mu$m) at a spectral resolution of 500. We developed a method to correct chromatic dispersion and non-common paths effects in the fringe phase and modelled the planet astrometry and stellar contamination. We obtained a high-signal-to-noise spectrum of $\beta$ Pictoris b, showing the planet continuum in L (for the first time) and M bands, which contains broad absorption features of H$_2$O and CO. In conjunction with a new GRAVITY spectrum, we modelled it with the ForMoSA nested sampling tool and the Exo-REM grid of atmospheric models, and found a solar carbon-to-oxygen ratio in the planet atmosphere. This study opens the way to the characterization of fainter and closer-in planets with MATISSE, which could complement the JWST at angular separations too close for it to obtain exoplanet spectra. Starting in 2025, the new adaptive optics system brought by the GRAVITY+ upgrade will further extend the detection limits of MATISSE.

Knowledge of the spatially resolved star formation history (SFH) of disk galaxies provides crucial insight into disk assembly, quenching, and chemical evolution. However, most reconstructions, both for the Milky Way and for external galaxies, implicitly assume that stars formed at their present-day radii. Using a range of zoom-in cosmological simulations, we show that stellar radial migration introduces strong and systematic biases in such SFH estimates. In the inner disk (R < h_d), early star formation is typically underestimated by 25-50% and late star formation overestimated, giving the misleading impression of prolonged, moderate activity. An exception occurs in the very central bin considered (~ 0.4h_d), which is consistently overestimated due to a net inflow of inward migrators. At intermediate radii and in the outer disk, migration drives the opposite trend: intermediate-age populations are overestimated by 100-200% as stars born in the inner disk migrate outward, whereas genuinely in-situ populations are underestimated by ~ 50% as they themselves continue to migrate. The net effect is that SFH peaks are suppressed and broadened, and the true rate of inside-out disk growth is systematically underestimated. These distortions affect all galaxies in our sample and have direct implications for interpreting spatially resolved SFHs from IFU surveys such as CALIFA and MaNGA, where present-day radii are often used as proxies for stellar birth sites. Correcting these biases will require accounting for disk mass, bar presence, disk kinematics and morphology, while recent birth-radius estimation techniques for Milky Way stars offer a promising path forward.

We present a systematic investigation of the evolution of the mass-metallicity relation (MZR) and fundamental metallicity relation (FMR) using uniform metallicity diagnostics across redshifts $z\sim0$ to $z\sim3.3$. We present new Keck/DEIMOS measurements of the [OII]$\lambda\lambda3726,3729$ emission line doublet for star-forming galaxies at $z\sim1.5$ with existing measurements of redder rest-optical lines from the MOSDEF survey. These new observations enable uniform estimation of the gas-phase oxygen abundance using ratios of the [OII], H$\beta$, and [OIII] lines for mass-binned samples of star-forming galaxies in 6 redshift bins, employing strong-line calibrations that account for the distinct interstellar medium ionization conditions at $z<1$ and $z>1$. We find that the low-mass power law slope of the MZR remains constant over this redshift range with a value of $\gamma=0.28\pm0.01$, implying the outflow metal loading factor ($\zeta_\text{out}=\frac{Z_{\text{out}}}{Z_{\text{ISM}}}\frac{\dot{M}_{\text{out}}}{\text{SFR}}$) scales approximately as $\rm \zeta_{out}\propto M_*^{-0.3}$ out to at least $z\sim3.3$. The normalization of the MZR at $10^{10}\ \text{M}_\odot$ decreases with increasing redshift at a rate of $d\log(\text{O/H})/dz =-0.11\pm0.01$ across the full redshift range. We find that any evolution of the FMR is smaller than 0.1 dex out to $z\sim3.3$. We compare to cosmological galaxy formation simulations, and find that IllustrisTNG matches our measured combination of a nearly-invariant MZR slope, rate of MZR normalization decrease, and constant or very weakly evolving FMR. This work provides the most detailed view of MZR and FMR evolution from the present day through Cosmic Noon with a fine time sampling of $1-3$ Gyr, setting a robust baseline for metallicity evolution studies at $z>4$ with JWST.

Dust plays a pivotal role in shaping the observed morphology of galaxies. While traditional cosmological simulations often assume a fixed dust-to-gas (DTG) or dust-to-metal (DTM) mass ratio to model dust effects, recent advancements have enabled on-the-fly (OTF) dust modeling that captures the spatial and temporal evolution of dust. In this work, we investigate the impact of dust modeling on galaxy morphology using the NewCluster simulation, which implements a detailed OTF dust model. We generate mock images of NewCluster galaxies under both OTF and fixed DTM models using the radiative transfer code SKIRT, and compare their morphology to JWST observations. We measure morphology indices and use the $G-M_{20}$ test to classify galaxies. We find that the OTF galaxy models exhibit brighter centers and more pronounced bulges than those of the fixed DTM models, resulting in a lower late-type galaxy (LTG) fraction, particularly at high redshifts. This central brightening is linked to a phenomenon we refer to as the DTM cavity, a localized depression in the DTM ratio driven by intense bulge starbursts. Our results highlight the importance of modeling dust evolution in a physically motivated manner, as fixed DTM models fail to capture key morphological features.

We report VLT spectroscopy of the interstellar comet 3I/ATLAS (C/2025 N1) from $r_{\rm h}\!\simeq\!4.4$ to $2.85$ au using X-shooter (300-550 nm, $R\!\simeq\!3000$) and UVES (optical, $R\!\simeq\!35k-80k$). The coma is dust-dominated with a fairly constant red optical continuum slope ($\sim$21-22\%/1000Å). At $r_{\rm h}\!\simeq\!3.17$ au we derive $3\sigma$ limits of $Q({\rm OH})<7.76\times10^{23}\ {\rm s^{-1}}$, but find no indications for [O I], C$_2$, C$_3$ or NH$_2$. We report detection of CN emission and also detect numerous Ni I lines while Fe I remains undetected, potentially implying efficiently released gas-phase Ni. From our latest X-shooter measurements conducted on 2025-08-21 ($r_{\rm h} = 2.85$\,au) we measure production rates of $\log~Q(\mathrm{CN}) = 23.61\pm 0.05$ molecules s$^{-1}$ and $\log~Q$(Ni) $= 22.67\pm0.07$ atoms s$^{-1}$, and characterize their evolution as the comet approaches perihelion. We observe a steep heliocentric-distance scaling for the production rates $Q(\mathrm{Ni}) \propto r_h^{-8.43 \pm 0.79}$ and for $Q(\mathrm{CN}) \propto r_h^{-9.38 \pm 1.2}$, and predict a Ni-CO$_{(2)}$ correlation if the Ni I emission is driven by the carbonyl formation channel. Energetic considerations of activation barriers show that this behavior is inconsistent with direct sublimation of canonical metal/sulfide phases and instead favors low-activation-energy release from dust, e.g. photon-stimulated desorption or mild thermolysis of metalated organics or Ni-rich nanophases, possibly including Ni-carbonyl-like complexes. These hypotheses are testable with future coordinated ground-based and space-based monitoring as 3I becomes more active during its continued passage through the solar system.

Wolf-Rayet (WR) stars are the evolved descendants of the most massive stars and show emission-line dominated spectra formed in their powerful stellar winds. Marking the final evolution stage before core collapse, the standard picture of WR stars has been that they evolve through three well-defined spectral subtypes known as WN, WC, and WO. Here, we present a detailed analysis of five objects that defy this scheme, demonstrating that WR stars can also evolve directly from the WN to the WO stage. Our study reveals that this direct transition is connected to low metallicity and weaker winds. The WN/WO stars and their immediate WN precursors are hot and emit a high flux of photons capable of fully ionizing helium. The existence of these stages unveil that high mass stars which manage to shed off their outer hydrogen layers in a low-metallicity environment can spend a considerable fraction of their lifetime in a stage that is difficult to detect in integrated stellar populations, but at the same time yields hard ionizing flux. The identification of the WN to WO evolution path for massive stars has significant implications for understanding the chemical enrichment and ionizing feedback in star-forming galaxies, in particular at earlier cosmic times.

The escape of Lyman continuum (LyC) radiation from early galaxies transformed the intergalactic medium (IGM) and is intimately connected to the fueling and feedback processes that regulate galaxy evolution. IGM attenuation interferes with high-redshift LyC observations, but growing samples of LyC observations at z<0.1 are revealing the properties of LyC-emitting galaxies. Along with multi-wavelength observations of nearby LyC-emitting candidates, cosmological simulations, and simulations of LyC escape from star-forming clouds, recent studies are providing insights into the physics of LyC escape and the possible characteristics of the galaxies that reionized the universe. Here, I review progress in LyC detections, the inferred indirect signatures of LyC escape and their application to high redshift, and our current understanding of the physical conditions that lead to high LyC escape. These findings include: LyC-emitting populations are diverse, and multiple factors correlate with LyC escape, particularly neutral gas absorption, dust attenuation, nebular ionization, and concentrated star formation. Radiative feedback plays a critical role in the youngest starbursts with the highest LyC escape fractions, but mechanical feedback may also contribute. Further research is needed to clarify the timing and role of different feedback mechanisms and to connect local LyC-production sites with the broader interstellar medium. Indirect LyC diagnostics show promise, but we need to understand whether and how the properties of LyC-emitting galaxies evolve from low to high redshift.

Elements heavier than hydrogen and helium, collectively termed metals, were created inside stars and dispersed through space at the final stages of stellar evolution. The relative amounts of different isotopes (variants of the same element with different masses) in stellar atmospheres provide clues about how our galaxy evolved chemically over billions of years. M dwarfs are small, cool, long-lived stars that comprise three-quarters of all stars in our galaxy. Their spectra exhibit rich fingerprints of their composition, making them potential tracers of chemical evolution. We measure rare carbon and oxygen isotopes in 32 nearby M dwarfs spanning a range of metallicities using high-resolution infrared spectroscopy. We find that stars with higher metal content have lower 12C/13C ratios, indicating they formed from material progressively enriched in 13C over time. This pattern is consistent with models where novae eruptions contributed significant amounts of 13C to the interstellar medium over the past few billion years. Our measurements of 16O/18O ratio match theoretical predictions and suggest that metal-rich stars reach 16O/18O ratios significantly lower than the Sun. These results establish M dwarfs as tracers of chemical enrichment throughout cosmic history.

If dark matter (DM) exists in halos around rotating neutron stars (NSs), it will be essential to understand the effects of rotation on the distribution of DM and baryonic matter (BM) in the stars to interpret observations. In this work, we construct rapidly rotating dark matter admixed neutron stars (DANS) with DM halos using the two-fluid approximation, where the BM and DM interact only through gravity. Our goal is to describe rapidly rotating millisecond-period DANS spun up by the accretion of BM from a zero angular momentum state. We extend the Rapidly Rotating Neutron Star (RNS) code to compute axisymmetric configurations in which the BM rotates rigidly while the DM remains torque-free and differentially rotates through the frame-dragging of spacetime. For the first time, we examine in detail local and global definitions of mass in general relativity for two-fluid systems, showing how their differences affect the interpretation of baryonic and dark component masses. We compute energy density and frame-dragging frequency profiles for DANS with three different characteristic DM halos. We demonstrate that rapid BM rotation reduces DM halo sizes if central energy densities are kept constant between non-rotating and rotating models. We also construct sequences of DANS to create mass and radius curves and compare rotating and non-rotating cases. Finally, we quantify deviations in the spacetime metric outside the baryonic surfaces of these sequences of stars caused by the DM halos. We hypothesize that the size of this quantity could indicate whether a DM halo will significantly impact X-ray pulse profile modeling. These results provide a framework for assessing the observational consequences of DM halos around rapidly rotating NSs.

Active galactic nuclei (AGNs) often exhibit broad-line regions (BLRs), populated by high-velocity clouds in Keplerian orbits around the central supermassive black hole (SMBH) at subparsec scales. During episodes of intense super-Eddington accretion, the disk can launch a powerful radiation-driven wind that overtakes the BLR clouds, forming bowshocks. Two shocks arise: one into the wind and another into the cloud. If adiabatic, electrons and protons are efficiently accelerated via Fermi processes to relativistic energies. In dense winds, the resulting high-energy photons are absorbed and reprocessed within the photosphere, while neutrinos from inelastic $pp$ collisions escape. We explore the potential of super-accreting AGNs as neutrino sources and propose a new class of emitter: an AGN without jets or gamma-ray counterparts, but with a strong opaque disk wind. As a case study, we consider a SMBH with $M_{\rm BH}=10^6,M_{\odot}$ and accretion rates consistent with tidal disruption events (TDEs). We compute the main cooling processes for relativistic particles and show that super-Eddington SMBHs can produce detectable neutrino fluxes with only weak electromagnetic signatures. Such fluxes may be observable by IceCube-Gen2 in nearby galaxies with a high BLR cloud filling factor. For more massive black holes, detection remains possible with moderate filling factors if the source is close, or at larger distances if the filling factor is high. Our model thus provides a plausible scenario for extragalactic neutrino sources, where both flux and timescale are determined by the number of orbiting clouds and the duration of the super-accreting phase.

Wide separation gas giant planets present a challenge to current planet formation theories, and the detection and characterisation of these systems enables us to constrain their formation pathways. The WIde Separation Planets In Time (WISPIT) survey aims to detect and characterise wide separation planetary-mass companions over a range of ages from <5 to 20 Myr around solar-type host stars at distances of 75-500 (median: 140) parsecs. The WISPIT survey carries out two 5 minute H-band exposures with the VLT/SPHERE instrument and IRDIS camera, separated by at least six months to identify co-moving companions via proper motion analysis. These two H-band observations in combination with a follow-up Ks-band observation were used to determine the colour-magnitude of the co-moving companions and to derive their masses by comparing to AMES-COND and AMES-DUSTY evolutionary tracks. We report the discovery of WISPIT 1b and WISPIT 1c, two gas giant exoplanets that are co-moving with the stellar binary WISPIT 1, which itself consists of a K4 star and M5.5 star in a multi-decadal orbit. The planets are at projected separations of 338 au and 840 au and have masses of 10 Mj and 4 Mj respectively. We identified two common proper motion planetary companions to a (previously unknown) stellar binary with a Sun-like primary. These targets are ideal for follow up characterisation with both ground and space-based telescopes. Monitoring of the orbit with the GRAVITY interferometer will place constraints on their eccentricity, and spectroscopic characterisation will identify the composition and metallicity, providing information on their formation pathways.

Cosmography is a model-independent phenomenological approach to observational cosmology, relying on Taylor series expansions of physical quantities as a function of the cosmological redshift or other analogous variables. A recent work developed the formalism for a cosmographic analysis of astrophysical and local measurements of the fine-structure constant, $\alpha$, and provided first constraints on the corresponding parameters. Here we update the earlier work, both by including more recent measurements of $\alpha$, and by extending it to other fundamental dimensionless couplings. We find no statistically significant evidence for such variations, and place stringent constraints on the first two terms of all these cosmographic series: the linear coefficient in the $\alpha$ series is constrained to parts per billion level, thanks to recently improved atomic clock constraints, while the other coefficients are constrained to parts per million level. Additionally, we use the same data to place cosmographic constraints on models from a broad class of Grand Unified Theories in which varying fundamental constants occur, and highlight how future data can provide a discriminating cosmographic test between freezing and thawing dark energy models.

On 2017-09-20 we observed GJ 4334, an M5V dwarf rotating with a period of 23.5 days, simultaneously with both the Space Telescope Imaging Spectrograph aboard Hubble (1160 -- 1710 Angstroms) and the Dual Imaging Spectrograph mounted on the 3.5m telescope at Apache Point Observatory (3750 -- 5050; 5800 -- 6950 Angstroms) as part of a larger survey of intermediately active M dwarfs. GJ 4334 flared during the observation, starting with a rise in the flux of optical chromospheric emission lines, followed by the rapid rise and decay of multiple far-ultraviolet emission lines formed in the transition region, followed by the slow decay of the optical lines. We find significant broadening and asymmetries in the optical emission lines that are potentially from bulk plasma motion, a post-flare elevated flux in both the optical and far-ultraviolet, and trends in the rise and decay timescales of the Balmer series such that higher-order lines rise earlier and decay faster than lower-order lines. The equivalent durations of the flare in individual lines range from 800 -- 3e4 seconds, mapping to flare energies of 1e28 -- 3e29 erg for each line. To contextualize GJ 4334's flare behavior we measure and compare its optical flare frequency distribution with TESS to EV Lacertae, a similar mass but faster rotating M dwarf, and find that GJ 4334 has an excess of large flares relative to the power-law established by the majority of its smaller flares. This dataset is a rare opportunity to characterize flares near a critical transition in stellar magnetic activity.

Core-collapse supernovae provide a unique opportunity to probe axions because they can be a copious source of the particles. It has recently been proposed that axion helioscopes can be used for the direct search for supernova axions if a supernova event appears within a few hundred parsecs. However, the event number of supernova axions has been estimated only within the post-process framework. In this study, we perform long-term supernova simulations for a 9.6M_sun star coupled with the axion emission to reevaluate the event number of axions detected by the helioscopes. We find that the additional cooling induced by the axion emission can significantly decrease the temperature in the proto-neutron star. As a result, the axion luminosity and hence the axion event number are reduced, compared with the result obtained through post-processing. Our result indicates that the nonlinear feedback of the axion emission is an essential factor to predict the axion detectability, and underscores the need for systematic simulation studies across various progenitor models.

Noncircular motions have been observed across various spatial scales in disk galaxies, yet the physical properties of the gas involved in these motions remain poorly constrained. Using data from 19 galaxies from the PHANGS-MUSE sample, we investigated the prevalence of noncircular flows at spatial resolutions of tens of parsecs. We developed a new tool for 3D kinematic modelling of data cubes and applied it to the PHANGS-MUSE H$\alpha$ spectral lines to recover the underlying circular, noncircular motions, as well as the intrinsic velocity dispersion in these objects. The PHANGS-MUSE galaxies exhibit rotation supported disks with $V_\mathrm{rot}/\sigma_\mathrm{intrin}$ ratios $\gtrsim$ 5. Our analysis revealed ionized gas exhibiting noncircular motions at different amplitudes, with low velocity amplitudes of about $5\mathrm{km\,s^{-1}}$ associated with the axisymmetric rotation component, deviations of $\sim10\mathrm{km\,s^{-1}}$ primarily linked to interarm and spiral arms, and larger deviations ($>20 \mathrm{km\,s^{-1}}$), found in the central and bar regions. We found that the velocity dispersion and the strength of ionization correlate with the amplitude of noncircular motions, suggesting that the underlying dynamics of the warm gas are closely tied to its physical properties.

The magnetic field is integral to our understanding of the formation and dynamical evolution of molecular clouds and star formation within. We present a polarimetric survey of 17 massive protostellar cluster forming clumps, covered in 34 pointings in the 230-GHz window using the Atacama Large Millimeter/submillimeter Array (ALMA). The two array configurations, C43-1 and C43-4, probe linearly polarized dust emission, hence the plane-of-the-sky orientation of magnetic fields, at resolutions of 1\arcsec\ and 0\arcsec.4 that correspond to approximately 0.01pc core and $10^3$ au envelope scales, respectively. The relative orientations (ROs) of the magnetic field probed at two spatial scales are analyzed for the entire protostellar cluster sample and for a subset of objects in NGC 6334. We found a bimodal distribution of ROs with peaks at 0° (parallel) and 90°(orthogonal) for the entire sample combined as well as for NGC 6334. We investigate the physical origin of this bimodal distribution through a projected Rayleigh statistic (PRS) analysis in relation to column densities and local gravity in NGC 6334. We found an excess of parallel magnetic fields at column densities $> 10^{23}$ \cmm. The underlying cause of the RO distribution of the magnetic field is gravitational collapse at higher gas densities, which drags and reorients the magnetic field as shown in the alignment between the magnetic field and the direction of gravitational forces. The distribution of ROs observed here is consistent with the evolution of relative orientations of an initially sub-Alvénic cloud that becomes magnetically super-critical and super-Alvénic as the cloud collapses to form stars.

Differential heating and radiation on asymmetric asteroids can cause measurable changes in their rotation rates and spin axes, known as the YORP effect. In binary systems, such radiation-driven torques can change the mutual asteroid orbits, termed the binary YORP or BYORP effect. To study how binary asteroid shapes and thermophysical properties affect surface temperatures and BYORP, we developed a new 3D thermophysical model which balances insolation, 1D conduction, visible light reflection, and mutual heating through scattered infrared radiation. Using 3D ray tracing, we include eclipses, shadowing from horizons and topography, and mutual radiation exchange between the primary and secondary asteroids. We perform global modeling of the binary asteroid (175706) 1996 FG3, a Janus mission target. At perihelion, we find that the 1996 FG3 system experiences temperatures between 100 and 475 K. We find that eclipses and thermal inertia can alter secondary surface temperatures by up to 14%, with a mean difference due to radiation from the primary of just over 1%. We also present a model for calculating the BYORP effect using binary thermophysical model results. This model compares well to analytical approximations of the BYORP coefficient B, and suggests that thermal effects like eclipses and thermal inertia can reduce torque in the 1996 FG3 system and alter the BYORP coefficient by up to several percent. For 1996 FG3, eclipses alter B by approximately 7%, resulting in a lower torque on the secondary. Though small, in the absence of tidal effects this would reduce the contraction of the semimajor axis by about 20 meters over 10,000 years. Our findings suggest that thermal effects can alter temperatures and BYORP calculations sufficiently that they should be included when modeling binaries. The relative importance of each effect is predicted to vary with the properties of the studied system.

We present a study of bright, young $\delta$ Scuti stars near the zero-age main sequence using TESS light curves and Gaia DR3 data. From a sample of 2041 stars with G<7 and $G_{BP}-G_{RP}$ colour in the range 0-0.6, we identified 444 $\delta$ Scuti pulsators. We measured a pulsator fraction of ~70% in the middle of the instability strip, tapering off towards the edges. A period - luminosity diagram reveals a concentration along the fundamental mode and overtone ridges. We addressed sample completeness and identified low-frequency pulsators possibly exhibiting mixed modes. Cross-matching with nearby young associations showed that 63 $\delta$ Scuti stars from our sample are association members.

We propose to search for a faint yet distinguishable contribution to the cosmic infrared background (CIB) spectrum arising from the radiative decay of the cosmic neutrino background (C$\nu$B). In the Standard Model of particle physics, neutrino decay is highly suppressed, with a predicted lifetime on the order of $10^{43}$ years. However, non-standard models suggest the possibility of significantly shorter lifetimes, ranging from $10^{12}$ to $10^{17}$ years. Observations to date, however, only provide a lower limit of approximately $10^{12}$ years for the neutrino lifetime. In PRIMA's low-resolution mode ($R \sim 100$), a diffuse background analysis, combined with the removal of point sources associated with known galaxies in a wide-field ($\sim 1$ square degree) spectroscopic survey covering in the 24-240 $\mu$m range could facilitate a search for neutrino decay lifetimes up to $O(10^{15})$ years. The expected signal from C$\nu$B decay at 50 $\mu$m for a neutrino lifetime of $O(10^{15})$ years is more than an order of magnitude fainter than the CIB and over three orders of magnitude fainter than the zodiacal emission foreground. Taking advantage of the characteristic spectral features of C$\nu$B decay, this level of sensitivity can be achieved with approximately 100 hours of total observation time, based on estimated surface brightness sensitivity. Searching for neutrino decay at this sensitivity could place strong constraints on several non-standard theories. A positive detection would provide compelling evidence for non-standard contributions to neutrino decay and could directly reveal the cosmic neutrino background. Furthermore, the decay photon spectrum could offer insights into the absolute mass of neutrinos, and key cosmological parameters.

Bibliometric methods provide valuable tools for assessing scientific productivity and impact across disciplines, yet their application in astronomy journals remains relatively limited. This study conducts a bibliometric analysis of Japanese astronomy publications before and after the commissioning of the Subaru Telescope, a major national investment in observational infrastructure. Using data from Scopus and SciVal, we examine peer-reviewed journal articles published between 1996 and 2007 by authors affiliated with Japanese institutions, focusing on field-normalized citation indicators such as the Field-Weighted Citation Impact (FWCI) and the share of publications in the top 10% most cited globally. Subaru Telescope-based publications are identified through cross-referencing with official telescope publication lists and are compared against national and global benchmarks. The results show that Subaru Telescope-based publications, while accounting for less than 10% of Japan's total scholarly output in astronomy, consistently achieved FWCI values above 2.0 and a significantly higher proportion of highly cited papers. This indicates that the Subaru Telescope substantially enhanced Japan's research visibility and impact, especially during its early operational years. This study demonstrates the utility of bibliometric evaluation in capturing the academic return of large-scale research facilities and contributes to broader discussions on research infrastructure in astronomy.

X-ray binaries play a significant role in the thermal and ionization history of galaxies. Their X-ray luminosity can shed light on galactic star formation rates and histories. Compact objects are also crucial in the evolution of gravitational wave progenitors. Here we present the results from our work to extend the binary population and spectral synthesis (BPASS) code suite to incorporate X-ray emission onto compact remnants in binary systems. We self-consistently model the accretion disc for each interacting binary system in a grid of stellar evolution models and then combine these to obtain the total X-ray spectra for stellar populations over a range of ages and metallicities. Crucially, these are estimated using the same stellar models as those used for modelling the stellar spectral energy distribution. We utilise first principle equations to calculate the X-ray binary (XRB) evolution, luminosity and spectral energy densities of individual accreting compact objects. Population synthesis using observationally motivated values for R_inner (the accretion disc inner truncation radius) reproduces the observed X-ray number evolution in the Small Magellanic Cloud and the inferred X-ray flux evolution for M51, validating our models. Using these models, we explore the implications of a self-consistent stellar and XRB emission population synthesis for ionizing photon production, the XRB dependence on metallicity and, for XRBs as a potential source of nebular He II emission seen in the spectra of high redshift galaxies. We conclude that XRBs contribute towards powering nebular He II emission without causing significant overestimates of hydrogen ionization.

The solar corona is much hotter than the photosphere and chromosphere, but the physical mechanism responsible for heating the coronal plasma remains unidentified yet. The thermal microwave emission, which is produced in strong magnetic field above sunspots, is a promising but barely exploited tool for studying the coronal magnetic field and plasma. We analyzed the microwave observations of eight solar active regions obtained with the Siberian Radioheliograph in years 2022-2024 in the frequency range of 6-12 GHz. We produced synthetic microwave images based on various coronal heating models, and determined the model parameters that provided the best agreement with the observations. The observations and simulations strongly favour either a steady-state (continuous) plasma heating process, or high-frequency heating by small energy release events with a short cadence. The average magnetic field strength in a coronal loop was found to decrease with the loop length, following a scaling law with the most probable index of about -0.55. In the majority of cases, the estimated volumetric heating rate was weakly dependent on the magnetic field strength, and decreased with the coronal loop length following a scaling law with the index of about -2.5. Among the known theoretical heating mechanisms, the model based on wave transmission or reflection in coronal loops acting as resonance cavities was found to provide the best agreement with the observations. The obtained results did not demonstrate a significant dependence on the emission frequency in the considered range.

The detection of gravitational waves from merging black holes with masses $\sim\,80-150\,\mathrm{M_\odot}$ suggests that some proportion of black hole binary systems form hierarchically in dense astrophysical environments, as most stellar evolution models cannot explain the origin of these massive black holes through isolated binary evolution. A significant fraction of such mergers could occur in Active Galactic Nuclei disks (AGN), however connecting individual black hole mergers to host galaxies is a challenging endeavor due to large localization uncertainties. We assess the feasibility of determining the fraction of hierarchically merging black hole binaries by computing the angular cross-correlation between gravitational wave localization posteriors and galaxy catalog skymaps. We forecast when the clustering of gravitational wave sky localizations can be measured accurately enough to distinguish the AGN origin scenario from hierarchical mergers in galaxies that do not host AGN. We find that if the observed merging population is dominated by binaries formed dynamically in AGN, then this could be determined with $\mathcal{O}(5000)$ mergers detected at the sensitivity that is projected for the upcoming A\# gravitational wave detectors.

The new 1.5~m telescope (AZ1500) is operating at the National Astronomical Observatory Rozhen, Bulgaria. This paper gives an overview of the telescope and presents a snapshot of the current performance. Science observations are under way, and we give brief highlights from a number of programs that have been enabled.

Are the most massive objects in the Universe today the direct descendants of the most massive objects at higher redshift? We address this question by tracing the evolutionary histories of haloes in the MultiDark Planck2 simulation. By following the 100 most massive halos at $z = 0$ across cosmic time, we find that only 40\% of them were among the largest 100 halos at $z = 1$. This suggests that many of today's most massive clusters were not the most dominant structures at earlier times, while some of the most massive objects at high redshift do not remain in the top mass ranks at later epochs. The hierarchical nature of structure formation predicts that, on average, massive haloes grow over time, with their abundance in comoving space decreasing rapidly at higher redshifts. However, individual clusters exhibit diverse evolutionary paths: some undergo early rapid growth, while others experience steady accretion or significant merger-driven mass changes. A key assumption in self-similar models of cluster evolution is that the most massive objects maintain their rank in the mass hierarchy across cosmic time. In this work, we test this assumption by constructing a mass-complete sample of haloes within the $(1 h^{-1}{\rm Gpc})^3$ volume of MultiDark and analysing when clusters enter and exit a high-mass-selected sample. Our results demonstrate that cluster selections must be carefully constructed, as significant numbers of objects can enter and leave the sample over time. These findings have important implications for observational cluster selection and comparisons between simulations and surveys, especially at high redshift.

(Abridged) Quasar studies with Herschel/SPIRE often report host luminosities ranging from 10^{12} to 10^{14} Lsun, suggestive of star formation rates (SFRs) of up to several thousand Msun/yr. Due to the limited spatial resolution of SPIRE, it is uncertain whether the far-infrared (FIR) emission originates from the quasar itself, nearby sources, or unrelated sources within the SPIRE beam. High-resolution observations at wavelengths close to the SPIRE coverage are needed to pinpoint the true source of the FIR emission. In this work we unambiguously identify the ALMA Band 7 counterparts of a statistical sample of 152 FIR-bright SDSS quasars and estimate the multiplicity rates among these systems. Based on the multiplicities, we assess the importance of mergers as triggers for concomitant accretion onto supermassive black holes (SMBHs) and extreme star formation. In ~60% of cases, the submm emission originates from a single counterpart within the SPIRE beam, centred on the optical coordinates of the quasar. The multiplicity rate increases by a factor of ~2.5 between redshifts 1 and 2.5. The incidence of multiplicities is consistent among broad absorption line (BAL) quasars and non-BAL quasars. The multiplicities observed in a fraction of the sample indicate that, while mergers enhance gas inflow efficiency, there must be viable alternatives for driving synchronous SMBH growth and intense star formation in isolated systems. We report the serendipitous detection of two CO(6-5) and three CO(7-6) transitions out of the eight such transitions expected based on the spectral setup and the redshifts of the objects in the sample. Higher transitions are not detected, indicating that the quasars are not exciting sufficiently the gas in their hosts. Finally, we also detect a potential emission of H2O, HCN (10-9) or a combination of both in the spectrum of a quasar at redshift 1.67.

We report the discovery of an exceptionally short burst recurrence time in the well-known clocked burster GS~1826$-$238, observed with the CubeSat X-ray observatory NinjaSat. The source had remained in a persistent soft spectral state since its hard-to-soft transition in 2015 July until a soft-to-hard transition occurred in 2025 May. On 2025 June 23, NinjaSat began monitoring GS~1826$-$238 in the hard state and continued until the source returned to a steady soft state. During this period, we detected 19 X-ray bursts: 14 during the hard state, 4 in the transitional state, and 1 in the soft state. In the hard state, we identified a new clocked bursting epoch, during which the burst recurrence time remained highly stable and unprecedentedly short among the clocked bursting phases of GS~1826$-$238, with $t_{\rm rec} = 1.603 \pm 0.040$~hr~($1\sigma$ error). Furthermore, this recurrence time deviates significantly by approximately~37\% from the previously established empirical relation $t_{\rm rec} \propto F_{\rm bol}^{-1.05}$, where $F_{\rm bol}$ is the bolometric persistent flux. These results suggest that the ignition conditions for X-ray bursts differed from those in earlier epochs. We explore possible explanations for this discrepancy and propose a scenario in which a smaller fraction of the neutron star surface was involved in fuel accumulation and nuclear burning. This situation leads to a high local accretion rate compared to earlier epochs, which naturally results in the exceptionally short recurrence time, as well as the observed reductions of approximately 38\% in blackbody normalization (proportional to the emitting area) and 29\% in burst fluence.

Context: Cosmic rays drive several key processes for the chemistry and dynamical evolution of star-forming regions. Their effect is quantified mainly by means of the cosmic-ray ionisation rate $\zeta_2$. Aims: We aim to obtain a sample of $\zeta_2$ measurements in 20 low-mass starless cores embedded in different parental clouds, to assess the average level of ionisation in this kind of sources and to investigate the role of the environment in this context. The warmest clouds in our sample are Ophiuchus and Corona Australis, where star formation activity is higher than in the Taurus cloud and the other isolated cores we targeted. Methods: We compute $\zeta_2$ using an analytical method based on the {column density} of ortho-$\rm H_2D^+$, the CO abundance, and the deuteration level of HCO$^+$. To estimate these quantities, we analysed new, high-sensitivity molecular line observations obtained with the Atacama Pathfinder EXperiment (APEX) single-dish telescope and archival continuum data from Herschel. Results: We report $\zeta_2$ estimates in 17 cores in our sample and provide upper limits on the three remaining sources. The values span almost two orders of magnitude, from $1.3 \times 10^{-18}\, \rm s^{-1}$ to $8.5 \times 10^{-17}\, \rm s^{-1}$. Conclusions: We find no significant correlation between $\zeta_2$ and the core's column densities $N\rm (H_2)$. On the contrary, we find a positive correlation between $\zeta_2$ and the cores' temperature, estimated via Herschel data: cores embedded in warmer environments present higher ionisation levels. The warmest clouds in our sample are Ophiuchus and Corona Australis, where star formation activity is higher than in the other clouds we targeted. The higher ionisation rates in these regions support the scenario that low-mass protostars in the vicinity of our targeted cores contribute to the re-acceleration of local cosmic rays.

The shape and dynamics of coronal mass ejections (CMEs) vary significantly based on the instrument and wavelength used. This has led to significant debate about the proper definitions of CME/shock fronts, pile-up/compression regions, and cores observations in projection in optically thin vs. optically thick emission. Here we present an observational analysis of the evolving shape and kinematics of a large-scale CME that occurred on May 7, 2021 on the eastern limb of the Sun as seen from 1 au. The eruption was observed continuously, consecutively by the Atmospheric Imaging Assembly (AIA) telescope suite on the Solar Dynamics Observatory (SDO), the ground-based COronal Solar Magnetism Observatory (COSMO) K-coronagraph (K-Cor) on Mauna Loa, and the C2 and C3 tele- scopes of the Large Angle Solar Coronagraph (LASCO) on the Solar and Heliospheric Observatory (SoHO). We apply the updated multi-instument version of the recently de- veloped Wavetrack Python suite for automated detection and tracking of coronal eruptive features to evaluate and compare the evolving shape of the CME front as it propagated from the solar surface out to 20 solar radii. Our tool allows tracking features beyond just the leading edge and is an important step towards semi-automatic manufacturing of train- ing sets for training data-driven image segmentation models for solar imaging. Our findings confirm the expected strong connection between EUV waves and CMEs. Our novel, de- tailed analysis sheds observational light on the details of EUV wave-shock-CME relations that is lacking for the gap region between the low and middle corona.

Superluminous supernovae (SLSNe) are a distinct class of stellar explosions, exhibiting peak luminosities 10-100 times brighter than those of normal SNe. Their extreme luminosities cannot be explained by the radioactive decay of $^{56}\mathrm{Ni}$ and its daughter $^{56}\mathrm{Co}$ alone. Consequently, models invoking newly formed millisecond magnetars have been widely proposed, capable of supplying additional energy through magnetic dipole radiation. For these rapidly rotating magnetars, however, gravitational-wave (GW) emission may also contribute significantly to the spin-down, particularly during their early evolutionary stages. While high-energy photons initially remain trapped within the optically thick ejecta, they will eventually escape as the ejecta becomes transparent during the expansion, thereby influencing the late-time lightcurve. In this work, we adopt an analytical framework to systematically explore the combined effects of GW emission and high-energy leakage on the lightcurve of SLSNe. Compared to scenarios that neglect these processes, we find that for magnetars with initial spin periods of millisecond, the combined influence suppresses early-time luminosities but enhances late-time emission. We further investigate the effects of the neutron-star equation of state to the lightcurve, GW emission efficiency, ejecta mass, and other relevant quantities. Our results highlight the complex interplay between GW-driven spin-down and radiative transport in shaping the observable features of SLSNe, offering new insights into diagnosing the nature of their central engines.

Atmospheres of transiting exoplanets can be studied spectroscopically using space-based or ground-based observations. Each has its own strengths and weaknesses, so there are benefits to both approaches. This is especially true for challenging targets such as cooler, smaller exoplanets whose atmospheres likely contain many molecular species and cloud decks. We aim to study the atmosphere of the warm Neptune-like exoplanet WASP-107 b (Teq~740 K). Several molecular species have been detected in this exoplanet in recent space-based JWST studies, and we aim to confirm and expand upon these detections using ground-based VLT, evaluating how well our findings agree with previously retrieved atmospheric parameters. We observe two transits of WASP-107 b with VLT/CRIRES+ and create cross-correlation templates of the target atmosphere based on retrieval results from JWST studies. We create different templates to investigate the impact of varying volume mixing ratios of species and inclusion or exclusion of clouds. Considering this target's observational challenges, we create simulated observations prior to evaluating real data to assess expected detection significances. We report detections of two molecular species, CO (~6 S/N) and H2O (~4.5 S/N). This confirms previous space-based detections and demonstrates, for the first time, the capability of VLT/CRIRES+ to detect species in targets cooler than hot Jupiters using transmission spectroscopy. We show our analysis is sensitive to cloud inclusion, but less so to different volume mixing ratios. Interestingly, our detection deviates from its expected location in our Kp-vsys diagrams, and we speculate on possible reasons for this. We demonstrate that the error budget for relatively cooler exoplanets is severely reduced in comparison to hotter exoplanets, and underline need for further work in context of high-resolution spectroscopy.

Context: Observations suggest that magnetic fields of disk-bearing stars may have non-dipolar configurations. However, the influence of these configurations on magnetospheric accretion remains poorly understood. Aims: We aim to simulate magnetospheric accretion incorporating non-dipolar and strong magnetic field. Our model is informed by observations of IRAS 17449+2320, a post-merger belonging to the group of FS CMa stars, which indicate a dominant dipolar magnetic field with an additional quadrupole component. Methods: Using the PLUTO code, we conduct 2.5-D non-ideal viscous-resistive magnetohydrodynamical (MHD) simulations of star-disk magnetospheric interactions. We consider a thin accretion disk and strong stellar magnetic field ($B_\star= 6.2\mathrm{kG}$) under four configurations: pure dipole, pure quadrupole, dipole plus quadrupole, and dipole plus octupole. In the latter two cases, different magnetic polar strength ratios are explored. Results: For asymmetric magnetic field configurations, we find that accretion exhibits funnel streams below the midplane, indicating the dominance of the quadrupolar and octupolar components. In contrast, in dipolar configurations, we observe the formation of two symmetrical funnels with respect to the midplane. However, in the quadrupolar configuration, accretion is entirely confined to the disk midplane forming a cone-like pattern that leads to disk widening. Remarkably, the presence of a quadrupolar component gives rise to highly asymmetric substructures in the corona region. Conclusions: Multipolar stellar magnetic fields drive non-uniform accretion and lead to asymmetric density distributions in both the disk and corona. These results resemble observed features of some FS CMa post-mergers and Herbig Ae/Be stars, highlighting the critical role of magnetic field complexity in shaping circumstellar environments.

We report the discovery of diffuse radio emission within SNR G338.3-0.0 using new MeerKAT observations at 816 MHz and 1.4 GHz. The radio emission spatially overlaps with the X-ray pulsar wind nebula (PWN) powered by PSR J1640-4631 and the GeV/TeV gamma-ray source HESS J1640-465. The morphology of this radio emission is centrally peaked and its extent is well-contained within the SNR shell. A lack of mid- and far-infrared counterparts and the absence of catalogued H II regions argues against a thermal origin, while the morphology and radial profile are suggestive of a PWN origin powered by PSR J1640-4631. Under this assumption, we use a one-zone, time dependant model to reproduce the size and broadband (radio, X-ray, and gamma-rays) spectral energy distribution of the PWN. The modelling and broadband properties of this PWN suggests it is currently interacting with the reverse shock within its host SNR. This evolutionary stage is associated with particles escaping the PWN and entering the ISM, suggesting this object may be an important source of Galactic PeV e+/e-

Excellent (<25 mas) H$_{\alpha}$ images of the star TYC 5709-354-1 led to the discovery of a rare H$_{\alpha}$ protoplanet. This star was discovered by the WISPIT survey to have a large multi-ring transitional disk, and is hereafter WISPIT 2. Our H$_{\alpha}$ images of 2025, April 13 and April 16 discovered an accreting (H$_{\alpha}$ in emission) protoplanet: WISPIT 2b (r=309.43$\pm$1.56 mas; (~54 au deprojected), PA=242.21$\pm$0.41 degrees) likely clearing a dust-free gap between the two brightest dust rings in the transitional disk. Our SNR=12.5 detection gave an H$_{\alpha}$ ASDI contrast of (6.5$\pm$0.5)x10$^{-4}$ and a H$_{\alpha}$ line flux of (1.29$\pm$0.28)x10$^{-15}$ erg/s/cm$^2$. We also present L' photometry from LBT/LMIRcam of the planet (L'=15.30$\pm$0.05 mag) which, when coupled with an age of 5.1$^{+2.4}_{-1.3}$ Myr, yields a planet mass estimate of 5.3$\pm$1.0 Mjup from the DUSTY evolutionary models. WISPIT 2b is accreting at 2.25$^{-0.17}_{+3.75}$x10$^{-12}$ Msun/yr. WISPIT 2b is very similar to the other H$_{\alpha}$ protoplanets in terms of mass, age, flux, and accretion rate. The inclination of the system (${\it i}$=44 degrees) is also, surprisingly, very similar to the other known H$\alpha$ protoplanet systems which all cluster from 37$\leq{\it i}\leq$52 degrees. We argue this clustering has only a ~1.0% (2.6 sigma) probability of occurring randomly, and so we speculate that magnetospherical accretion might have a preferred inclination range (~37-52 degrees) for the direct (cloud free, low extinction) line of sight to the H-alpha line formation/shock region. We also find at 110mas (~15au deprojected) a close companion candidate (CC1) which may be consistent with an inner dusty 9$\pm$4 Mjup planet.

In the past decades several thousand exoplanet systems have been discovered around evolved, main-sequence stars, revealing a wide diversity in their architectures. To understand how the planet formation process can lead to vastly different outcomes in system architecture we have to study the starting conditions of planet formation within the disks around young stars. In this study we are presenting high resolution direct imaging observations with VLT/SPHERE of the young ($\sim$5 Myr), nearby ($\sim$133 pc), solar-analog designated as WISPIT 2($=$ TYC~5709-354-1). These observations were taken as part of our survey program that explores the formation and orbital evolution of wide-separation gas giants. WISPIT 2 was observed in four independent epochs using polarized light and total intensity observations. They reveal for the first time an extended (380 au) disk in scattered light with a multi-ringed sub-structure. We directly detect a young proto-planet WISPIT 2b, embedded in a disk gap and show that it is co-moving with its host star. Multiple SPHERE epochs demonstrate that it shows orbital motion consistent with Keplerian motion in the observed disk gap. Our $H$ and $K_s$-band photometric data are consistent with thermal emission from a young planet. By comparison with planet evolutionary models, we find a mass of the planet of $4.9^{+0.9}_{-0.6}$ Jupiter masses. This mass is also consistent with the width of the observed disk gap, retrieved from hydrodynamic models. WISPIT 2b is the first unambiguous planet detection in a multi-ringed disk, making the WISPIT 2 system the ideal laboratory to study planet-disk interaction and subsequent evolution.

In the search for life in the Universe, molecular oxygen (O$_2$) combined with a reducing species, such as methane (CH$_4$), is considered a promising disequilibrium biosignature. In cases where it would be difficult or impossible to detect O$_2$ (e.g., mid-IR or low O$_2$ levels), it has been suggested that ozone (O$_3$), the photochemical product of O$_2$, could be used as a proxy for determining the abundance of O$_2$. As the O$_2$-O$_3$ relationship is nonlinear, the goal of this series of papers is to explore how it would change for different host stars and atmospheric compositions and learning how to use O$_3$ to infer O$_2$. We used photochemistry and climate modeling to further explore the O$_2$-O$_3$ relationship by modeling Earth-like planets with the present atmospheric level (PAL) of O$_2$ between 0.01% to 150% along with high and low CH$_4$ abundances of 1000% and 10% PAL, respectively. Methane is of interest not only because it is a biosignature, but also the source of hydrogen atoms for hydrogen oxide (HO$_x$) which destroys O$_3$ through catalytic cycles and acts as a catalyst for the smog mechanism of O$_3$ formation in the lower atmosphere. We find varying CH$_4$ causes changes to the O$_2$-O$_3$ relationship in ways that are highly dependent on both host star and O$_2$ abundance. A striking result for high CH$_4$ models in high O$_2$ atmospheres around hotter hosts is that enough CH$_4$ is efficiently converted into H$_2$O to significantly impact stratospheric temperatures, and therefore the formation/destruction rates of O$_3$. Changes in HO$_x$ also influenced the HO$_x$ catalytic cycle and smog O$_3$, causing variations in harmful UV reaching the surface as well as changes in the 9.6~$\mu$m O$_3$ feature in emission spectra. This demonstrates the need to explore the O$_2$-O$_3$ relationship in order to use O$_3$ as a reliable proxy for O$_2$ in future observations.

Cosmological measurements have revealed tensions within the standard $\Lambda$CDM model, notably discrepancies in the Hubble constant and $S_8$ parameter. Additionally, recent observations suggested evidence of a possible non-flat spatial curvature and an anomalous CMB lensing amplitude that beyond the $\Lambda$CDM framework. In this study, we explore whether introducing a variation in the electron mass $m_e$, allowing non-zero spatial curvature $\Omega_K$, and a free lensing amplitude $A_{\rm lens}$ can resolve these persistent tensions. Without relying on any local Universe measurements, we obtain $H_0 = 69.61^{+0.60}_{-0.55} \rm \, km \, s^{-1} \, Mpc^{-1}$ and $S_8= 0.808\pm0.012$, with $\Delta m_e / m_e = 0.0109^{+0.0066}_{-0.0062}$ and $A_{\rm lens} = 1.030^{+0.039}_{-0.037}$, both of which exceed the $\Lambda$CDM expectations. We find no indication of spatial curvature deviating from flatness, even when including the Cosmic Chronometers and SNe Ia samples. Our results point to a promising direction for cosmological models to reconcile these discrepancies, although more precise data from future experiments will be necessary to clarify these modifications.

Motivated by recent results from the DESI collaboration, we explore two classes of quintessence models that can give rise to crossing of the dark energy equation of state through the ``phantom divide'' $w=-1$. These are models with Lagrangians that involve higher powers of the kinetic energy $\dot\phi^2$, or where the dark matter (DM) mass is a function of $\phi$. Both have similar features with respect to the reconstructed redshift-dependent $w(z)$: moderate tuning of parameters is required to achieve the desired shape, and it is difficult or impossible for $w(z)$ to continue evolving smoothly as $z$ becomes large. Nevertheless, they give a strong improvement over $\Lambda$CDM in fitting the data. We point out that models of coupled dark matter and dark energy that cross the phantom divide are under pressure from constraints on long-range DM forces. They rule out the simplest renormalizable coupling of scalar DM to quintessence, but leave the fermionic case marginally allowed, while exponentially coupled models are safe from current constraints.

We study the galaxy bispectrum multipoles in the Hu-Sawicki $f(R)$ gravity model, where a scalar degree of freedom mediates a fifth force that is screened in high-density environments. The model is specified by $f_{R0}$, the present-day background value of the scalar field, which controls the strength of deviations from General Relativity (GR). Using perturbation theory, we compute the redshift-space galaxy bispectrum with the full scale- and time-dependent second-order kernels, incorporating corrections from the scale-dependent growth rate and nonlinear screening. Expanding the bispectrum in spherical harmonics, we analyze the sensitivity of the multipoles to modified gravity and forecast their detectability in a \textit{Euclid}-like survey. The monopole ($B_0^0$) and quadrupole ($B_2^0$) show the strongest signatures, with relative deviations of $2\%$-$8\%$ at $z=0.7$ and $k_1\simeq0.3\,h\,{\rm Mpc}^{-1}$ for $f_{R0}=10^{-5}$. Higher multipoles provide weaker but complementary signals. Finger-of-God damping and galaxy bias significantly modulate the results, requiring joint modeling of these nuisance parameters. For \textit{Euclid}, we forecast signal-to-noise ratios up to $\sim30$ for the monopole and $\sim15$ for the quadrupole. These results demonstrate that bispectrum multipoles are a powerful probe of gravity, capable of breaking degeneracies with bias and velocity effects and strengthening constraints on deviations from $\Lambda$CDM.

In interacting dark energy (DE) and dark matter (DM) scenarios, the interaction function typically includes a coupling parameter $\xi$ that quantifies the strength of energy exchange between the dark sectors. While $\xi$ is often assumed to be constant, there is no fundamental reason to exclude a time-dependent coupling, which could provide a more general and realistic description of dark sector dynamics. In this work, we study two widely used interacting models involving pressureless DM and DE, where the coupling parameter is allowed to vary with the scale factor $a$. Specifically, we consider two parametrizations: $\xi(a) = \xi_0 + \xi_a (1-a)$ and $\xi(a) = \xi_0 \left(1 + \frac{1-a}{a^2 + (1-a)^2} \right)$, and constrain them using the latest cosmological observations, including Planck 2018 CMB data, DESI DR2 BAO measurements, and multiple Type Ia supernovae samples. Our results show that one scenario yields evidence for a non-zero interaction at more than 95\% confidence level, while the remaining cases indicate at most mild or inconclusive signs of interaction. These findings highlight the potential of variable coupling models and the importance of continued investigation into the nature of the dark sectors.

Since its observation in 2019, the first image of a super-massive black hole using Very Long Baseline Interferometry (VLBI) with an Earth-scale baseline has generated much scientific and public interest, including the possible extension of the baseline into space to obtain higher image resolution. Operating one or more VLBI nodes in space will require the use of frequency standards that are space qualified, greatly reducing the number of options available. The coherence function C(T) is the metric usually used to determine the viability of a frequency standard. Here we show that C(T) is a useful but not sufficient metric for gauging frequency standard performance in VLBI and instead derive an expression for the clock-limited VLBI visibility S/N. We evaluate this expression for real frequency standards and find only the Ultra-Stable Oscillator (USO) and hydrogen maser to be viable for upcoming high-frequency VLBI with the USO only useful for very limited integration times (30s at 90 GHz, 10s at 230 GHz, 5s at 345 GHz, and not viable at 630 GHz). The maser extends these, but may have prohivitive size for a space mission. We also evaluate emerging frequency standard technologies and find the Optical Local Oscillator portion of optical clocks to be very promising (conservatively >100s at 90 GHz, 60s at 230 GHz, 40s at 345 GHz, and 22s at 630 GHz) when accounting for both performance and potential operation in space.

Neutral hydrogen (HI) 21-cm Intensity Mapping (IM) holds the potential to map the large-scale structures in the Universe over a wide redshift range $(z \lesssim 5.5)$, measure cosmological parameters, and shed light on the nature of dark energy. In addition, the signal is also sensitive to how the HI is distributed among the dark matter haloes, this being quantified through the HIHM relation, which relates the HI mass to the halo mass. In this work, we investigate whether measurements of the 21-cm power spectrum (PS) and bispectrum (BS) at large scales can be used to estimate the HIHM relation, which quantifies the HI distribution at small scales. As a proof of concept, we consider the simulated 21-cm IM signal at $z=1$. We find that the measured 21-cm PS and BS at large scales $(k \le k_{ul} = 0.32 \, {\rm Mpc}^{-1})$ are well modeled using perturbation theory, with only two free parameters namely $[\Omega_{\rm HI} b_1]$ and $\gamma = b_2/b_1$. Combining the measured 21-cm PS and BS with an independent measurement of $\Omega_{\rm HI} $, we show that it is possible to estimate the three parameters that quantify the HIHM relation. We expect observational estimates of the HIHM relation to shed light on galaxy formation and the evolution of the ISM. Our preliminary analysis ignores redshift space distortion and the system noise in IM observations, which we plan to address in future work.

Stars and planets form in collapsing clouds of gas and dust. The presence of dust grains and their local distribution play a significant role throughout the protostellar sequence, from the thermodynamics and the chemistry of molecular clouds to the opacity of collapsing protostellar cores and the coupling between the gas and the magnetic field and down to planet formation in young and evolved disks. We aim to simulate the dynamics of the dust, considering the whole range of grain sizes, from few nanometers to millimeters. We implemented a neutral pressureless multifluid that samples the dust size distribution in the RAMSES code. This multifluid is dynamically coupled to the gas via a drag source term and self-gravity, relying on the Eulerian approach. We designed a Riemann solver for the gas and dust mixture that prevents unphysical dust-to-gas ratio variations for well-coupled grains. We illustrated the capacities of the code by performing simulations of a protostellar collapse down to the formation of a first hydrostatic core, both for small and large dust grains. Grains over 100 microns significantly decouple from the gas. The spatial maps and the probability density functions indicate that dust enrichment within the first hydrostatic core and in some locations of the envelope increases as a function of the grain size and the level of initial turbulence. Thanks to the novel Riemann solver, we recovered the terminal velocity regime, even at low resolution. Moreover, we successfully extended it to regimes where the grain inertia matters. The multifluid module performs the coupling between the dust and the gas self-consistently all through the dynamical scales. The dust enrichment in the first hydrostatic core and the envelope have been revised here, assuming the initial turbulence and grain sizes. This enables us to probe new potential conditions for planet formation.

Massive quiescent galaxies at high redshift show significantly more compact morphology than their local counterparts. To examine their internal structure across a wide redshift range and investigate potential redshift dependence, we performed spatially resolved SED fitting using pixedfit software on massive $(\log(M_*/M_\odot)\sim11)$ quiescent galaxies at $0<z<4$ with public James Webb Space Telescope and Hubble Space Telescope imaging data from the Public Release Imaging for Extragalactic Research and the Cosmic Evolution Early Release Science Survey. We find that at $z \sim 3.5$, the half-mass radius is about 5.4 times smaller than at $z \sim 0.5$. This growth is driven by stellar mass buildup in the outskirts ($r > 4$ kpc), while the central regions ($r \sim 1$ kpc) remain largely unchanged, with stellar mass surface density similar to local quiescent galaxies. The estimated star formation rates are too low to explain the stellar mass growth, indicating an additional stellar mass accumulation process, such as mergers, is necessary. We parameterize the size-mass relation of the most massive galaxies in our sample as $\log(R_{e,mass}) \propto \alpha \log(M_*)$, and find $\alpha = 2.67^{+1.14}_{-1.17}$ for $z\lessapprox2$, consistent with growth dominated by minor mergers, and $\alpha = 0.91^{+0.20}_{-0.16}$ for $z\gtrapprox2$, consistent with growth dominated by major mergers. These results indicate that massive quiescent galaxies originate from compact quenched systems and grow through combinations of minor and major mergers.

We introduce our new program to develop two-dimensional MEMS arrays of individually addressable micro-mirrors (''Micro-Mirror Devices'', MMDs) specifically optimized for astronomy, multi-slit spectroscopy in particular. After reviewing the main characteristics and performance of the currently available options, Micro Shutter Arrays by NASA/Goddard and Digital Micromirror Devices by Texas Instruments, we present our planned first generation/baseline devices with 30 micron x 30 miron pixel size arranged in a 1K x 1K format with tilt angle 15 degrees. Our goal is to bring to maturity a technology capable of delivering arrays of 2K x 2K element of 100 micron x 100 micron, buttable on two sides to achieve even larger formats. In additions to MEMS design, we will develop the associated device packaging and electronic control circuitry leveraging on the extensive expertise gained in the last 30+ years by leading experts from digital imaging industry.

This paper reports the observations of two coronal shocks from two Coronal Mass Ejections (CMEs) for the Successive type II Solar radio bursts observed on 02 May 2021 in the frequency range of 80 - 1 MHz with the time interval of ~ 20 minutes between them. Both the bursts show clear band splitting features in the harmonic band. The estimated heights for the source of the first type II burst lies in the range of 2.06 - 2.93 R0 with the average speeds of 601 + or - 76, 700 + or - 91 and 783 + or - 105 km/s for 2 X, 3.5 X and 5 X Saito electron density models, and the heights for the source of the second type II burst lies in the range of 2.24 - 3.83 R0 with the average speeds of 1063 + or - 113, 1287 + or - 145 and 1478 + or - 172 km/s. The successive CMEs are observed by the twin Solar Terrestrial Relations Observatory STEREO-A between 11:20 - 12:21 UT in the Extreme Ultra Violet Imager (EUVI) and in the Internally Occulting Refractive Coronagraphs (COR1) FOV, the two coronal shocks are generated by the two successive CMEs observed at (11:30) 11:26 UT and (11:55:00) 11:56 UT according to ST-A (EUVI) COR1 observations and most likely released from the same active region. The average speeds of CMEs at COR1 FOV are about 574 = or - 64 km/s and 595 + or - 82 km/s. The simultaneous observations of the EUV structures and the radio bursts, their coinciding height-time further confirms that the successive CMEs are responsible for the successive shocks and their related radio bursts in the corona. The observed band-splitting in the successive type-II radio bursts provides the compression ratios of 1.26 and 1.45 respectively. Therefore, these observations confirms the presence of shock waves in the corona.

We analyze binary black hole (BBH) mergers from the latest Gravitational Wave Transient Catalog (GWTC-3) using a flexible, non-parametric framework to infer the underlying black hole mass distribution. Our model employs Gaussian Processes (GPs) with astrophysically motivated priors to represent the mass distribution, assuming both component masses are independently drawn from a common mass function. This approach enables us to capture complex features in the population without relying on rigid parametric forms. Motivated by predictions from binary stellar evolution, we focus on the presence of a mass gap in the range $40-120~ M_\odot$, attributed to the pair-instability supernova (PISN) and pulsational PISN (PPISN) processes. Using our GP-based model, we find for the first time, strong evidence for the onset of this mass gap, locating its lower edge at approximately $45-60 ~M_\odot$. We observe a suppression in the BBH merger rate when comparing the latest constraints against our GP model, by a factor of $10-60$ at $\sim60~M_\odot$, corresponding to the minimum of the mass gap. Additionally, we find evidence of a subpopulation of mergers populating the mass gap around $70 ~M_\odot$, which we argue is due to hierarchical mergers, as well as of a feature in the $40-50~ M_\odot$ range, albeit with low significance, which we attribute to the predicted PPISN build up. We discuss the astrophysical implications of this result in light of GW231123, a recently reported BBH merger with a total mass potentially above the PISN gap, suggesting the need for revised models of massive stellar evolution or alternative formation channels.

Primordial Black Holes (PBHs) are compelling candidates for explaining the present-day relic abundance of cold dark matter (CDM), yet their formation typically requires finely tuned early-universe dynamics. In this work, we propose a novel PBH formation mechanism within a well-established magnetogenesis framework. This scenario simultaneously accounts for the large-scale magnetic fields observed today and generates an enhanced curvature power spectrum at intermediate scales, leading to PBH formation with masses that can survive until the present epoch. We identify a narrow reheating temperature range, $10^5\,\mathrm{GeV} \leq T_{re} \leq 3\times 10^5\,\mathrm{GeV}$, within which the resulting PBHs can constitute the entirety of the observed CDM abundance. Furthermore, our model predicts a stochastic gravitational wave (GW) background as a byproduct of the PBH formation process. Remarkably, the predicted GW signal lies within the sensitivity reach of upcoming space-based interferometers, such as the LISA, DECIGO, or SKA mission, offering a direct observational probe of this PBH generation mechanism.

The James Webb Space Telescope is characterising the atmospheres of sub-Neptunes. The presence of magma oceans on sub-Neptunes is expected to strongly alter the chemistry of their envelopes (100 bar-100 kbar) and atmospheres (1 mbar-100 bar). At the magma ocean-envelope boundary (MEB, >10 kbar), gas properties deviate from ideality, yet the effects of real gas behaviour on chemical equilibria remain underexplored. Here, we compute equilibrium between magma-gas and gas-gas reactions using real gas equations of state in the H-He-C-N-O-Si system for TOI-421b, a canonical hot sub-Neptune potentially hosting a magma ocean. We find that H and N are the most soluble in magma, followed by He and C. We fit real gas equations of state to experimental data on SiH$_4$, and show that, for a fully molten mantle, SiH$_4$ dominates at the MEB under accreted gas metallicity of 1$\times$ solar, but is supplanted by CH$_4$ at 100$\times$ solar. Lower mantle melt fractions lower both magma-derived Si abundances in the envelope and the solubility of H and He in magma, yielding H$_2$- and He-rich envelopes. Projecting equilibrium chemistry through the atmosphere (1 mbar-100 bar), we find that condensation of $\beta$-quartz 'clouds' depletes Si-bearing gases, although SiH$_4$ remains abundant at solar metallicity. SiH$_4$/CH$_4$ ratios increase with mantle melt fraction and decrease with metallicity. These effects also deplete H$_2$ via CH$_4$ and SiH$_4$ formation and the dissolution of H$_2$ into magma. The competition between SiH$_4$ and CH$_4$ presents a diagnostic of metallicity and magma oceans. The corollary is that H$_2$- and He-rich, SiH$_4$- and CH$_4$-poor (<5%) atmospheres may indicate a limited role or absence of magma oceans on sub-Neptunes.

We investigate a qualitatively new regime of inflationary models with small and rapid oscillations in the potential--resonant non-Gaussianity. In contrast to the standard scenario, where most of the observable information is encoded in the power spectrum, in this regime the oscillatory signal predominantly appears in higher-order correlation functions with large $n$. This behavior emerges when the oscillation frequency $\omega$ exceeds the naive cutoff of the theory, $4\pi f$. However, as noted by Hook and Rattazzi [2306.12489], the actual cutoff is somewhat higher--though only logarithmically--when the amplitude of the oscillations is small. We identify a phenomenologically relevant window in which $n$-point functions with $3 \lesssim n \lesssim 9$ are potentially observable. In this regime, the signal exhibits 350 - 1000 oscillations per decade in $k$.

The Vera C. Rubin Observatory's Legacy Survey of Space and Time will increase interstellar object detection rates to one every few months, substantially elevating the probability of identifying objects with characteristics suggesting artificial origin. Despite this imminent capability, no evidence-based crisis communication framework exists for managing potential technosignature discoveries. We present the Adaptive Communication Framework (ACF), a theoretically grounded protocol that integrates crisis communication theories with the SPECtrum of Rhetoric Intelligences model to address diverse cognitive, social and emotional processing styles. Through analysis of communication failures during COVID-19, Fukushima, and asteroid 99942 Apophis (2004 MN4), we identify critical gaps in managing scientific uncertainty under public scrutiny. The ACF provides graduated protocols calibrated to the Loeb Scale for Interstellar Object Significance, offering specific messaging strategies across four rhetoric intelligence channels (Systematic, Practical, Emotional, Creative) for each evidential level group. Recognizing that Artificial Intelligence (AI) systems will mediate public understanding, the framework incorporates safeguards against synthetic media manipulation and algorithmic misinformation through pre-positioned content seeding and deepfake detection protocols. Our theoretical framework reveals that successful communication of paradigm-shifting discoveries requires simultaneous activation of multiple cognitive channels, with messages adapted to cultural contexts and uncertainty levels. This framework provides essential theoretical groundwork for ensuring humanity's potentially most transformative discovery unfolds through understanding rather than chaos.

We present the first unified constraints on a broad class of extended dark matter compact objects (EDCOs) from interstellar gas heating. These include axion stars, Q-balls, axion miniclusters, dark fermion stars and primordial black holes surrounded by dark matter halos, which arise in a wide range of theories beyond the Standard Model. As such massive objects traverse the interstellar medium, their gravitational influence generates wakes and, if sufficiently compact, drives accretion flows that heat gas in their vicinity. Our general framework extends standard dynamical friction treatments by incorporating finite-size effects, internal density profiles, gas penetration through objects, and criteria for accretion disk formation. We perform detailed numerical calculations of wake formation and gas heating and apply our results to the Leo T dwarf galaxy, establishing new constraints on the dark matter fraction in EDCOs heavier than a solar mass spanning several orders of magnitude in both mass and abundance.

Scalar fields with non-trivial kinetic term derived from a nonlinear sigma model are motivated by UV completions of gravity such as string theory. We discuss the $\mathrm{SL}(2,\mathbb{R})$ and $\mathrm{O}(3)$ sigma models with interacting potentials and simulate their full nonlinear dynamics on black hole spacetimes. We study the properties of the field as a function of the curvature of the sigma model with respect to the free massive scalar case. In the accretion process, the $\mathrm{SL}(2,\mathbb{R})$ model behaves as a self-interacting field with attractive interaction, while the $\mathrm{O}(3)$ one exhibits a repulsive phenomenology. In the case of a binary black hole system, these models cause a positive or negative dephasing of the gravitational waveform, respectively, when compared to the non-interacting case and as long as the field's Compton wavelength is larger than the binary separation. We observe no bosenova emission, which suggests that nonlinearities tend to suppress this phenomenon. Our results highlight how the kinetic and potential terms are both relevant in determining the field dynamics.

The Earth's Moon presents a uniquely advantageous environment for detecting astrophysical gravitational waves (GWs) in the frequency range of millihertz to decihertz. Unlike Terrestrial GW detectors, the quiet seismic environment of the Moon does not impede detection in this band; in fact the ground motions of the Moon will be excited by GWs, making the Moon a resonant amplifier at low frequencies. The Laser Interferometer Lunar Antenna (LILA) mission aims to be limited by thermal Brownian noise in its optics across most target frequencies. By taking advantage of the lunar normal mode resonances, we show that the first phase of the mission, LILA Pioneer, achieves the GW sensitivity required to study astrophysical sources through the millihertz to decihertz range. The advanced phase of the mission, LILA Horizon, would increase GW sensitivity to the cosmological horizon in this band.

We investigate how screening mechanisms, reconciling light scalar fields driving cosmic acceleration with local fifth force constraints, can be probed via their impact on non-local quantum correlations between entangled spin pairs, whose evolution on a curved background is affected by General Relativity (GR) and screened modified gravity effects. We consider a gedankenexperiment featuring a pair of massive, spin-1/2 particles orbiting the Earth, evaluating their non-local correlations through spin observables associated to the Clauser-Horne-Shimony-Holt (CHSH) inequality. Using a general formalism developed earlier for curved space-time spin evolution, we compute the effects of screening on the CHSH inequality, finding its degree of violation to be suppressed relative to the flat space-time case. Applying this formalism to the chameleon, symmetron, and dilaton mechanisms, we identify currently unconstrained regions of parameter space where the screening contribution is comparable to that of GR. While detecting these effects will be challenging, our work provides a proof-of-principle for testing screened dark energy through quantum non-locality.

Neural density estimation has seen widespread applications in the gravitational-wave (GW) data analysis, which enables real-time parameter estimation for compact binary coalescences and enhances rapid inference for subsequent analysis such as population inference. In this work, we explore the application of using the Kolmogorov-Arnold network (KAN) to construct efficient and interpretable neural density estimators for lightweight posterior construction of GW catalogs. By replacing conventional activation functions with learnable splines, KAN achieves superior interpretability, higher accuracy, and greater parameter efficiency on related scientific tasks. Leveraging this feature, we propose a KAN-based neural density estimator, which ingests megabyte-scale GW posterior samples and compresses them into model weights of tens of kilobytes. Subsequently, analytic expressions requiring only several kilobytes can be further distilled from these neural network weights with minimal accuracy trade-off. In practice, GW posterior samples with fidelity can be regenerated rapidly using the model weights or analytic expressions for subsequent analysis. Our lightweight posterior construction strategy is expected to facilitate user-level data storage and transmission, paving a path for efficient analysis of numerous GW events in the next-generation GW detectors.

Ultralight bosons can be excited around spinning black holes via black hole superradiance. These boson clouds may play an important role in the orbital evolution of binary black holes. In this work, we investigate the formation and evolution of common envelopes of ultralight boson clouds in comparable mass-ratio black hole binaries. We describe the cloud evolution using gravitational molecular eigenstates and analyze the possible level transitions during orbital decay, as well as the impact on orbital dynamics. We find that the cloud can generally lead to eccentricity growth. In particular, the eccentricity may vary significantly during level transition, leaving an eccentricity of ${\cal O}(0.1)$ within the detection band of ground-based gravitational wave detectors.

The eclipsing binary 2M1510 AB made of two brown dwarfs, the supermassive black hole M87* and our solar system may all have in common the presence of a massive pointlike distant companion perturbing in some ways their inner dynamics. That is, respectively, a circumbinary exoplanet whose orbital plane should be perpendicular to that of its host inner binary, an intermediate mass black hole potentially causing the observed jet precession, and a still unseen planet, known as Telisto or Planet Nine, which may, perhaps, have been recently spotted in archival images although in a quite different orbit with respect to that originally expected. To this aim, the long-term rates of change of all the Keplerian orbital elements of the relative motion of a gravitationally bound two-body system acted upon by a massive remote perturber are analytically worked out in full generality and presented in a way facilitating a straightforward application to those and other potentially interesting astronomical scenarios. For each of the considered cases, exclusion plots in their parameter spaces are provided.

We consider the impact of the first string corrections of minimally coupled single scalar field theory on inflationary dynamics. Specifically we consider separately the string corrections $\sim \alpha'\xi(\phi)c_2\,\left( \partial_{\mu}\phi \partial^{\mu}\phi\right)^2$ and $\sim \alpha'c \square \phi \partial_{\mu}\phi \partial^{\mu}\phi$, where $\alpha'$ is the square of the string scale. Our aim is to develop a theory which is self consistent in the sense that the field equations reproduce themselves in the slow-roll approximation. Such a requirement for the theory with $\sim \alpha'\xi(\phi) c_2\left( \partial_{\mu}\phi \partial^{\mu}\phi\right)^2$ resulted to a trivial non-minimal coupling function $\xi(\phi)$, however a self-consistent framework emerged from the theory with correction term $\sim \alpha' c \square \phi \partial_{\mu}\phi \partial^{\mu}\phi$. The resulting theory can easily be worked out analytically and we obtained an inflationary theory that can easily be fitted with the Atacama Cosmology Telescope constraints on the scalar spectral index and the updated Planck constraints on the tensor-to-scalar ratio.