Papers for Tuesday, Dec 02 2025

It is argued that the lens described in the paper commented upon has a focussing efficiency of less than 0.03% and an angular resolution for broadband radiation that is an order of magnitude worse that the diffraction limit. Furthermore incident radiation that is not focussed will lead to background fog in the image plane.

We explore a phenomenological model of dark energy in which space is treated as an elastic brane with uniform tension $T_s$ supporting a longitudinal phonon fluid. The brane tension gives a residual geometric contribution to the vacuum energy, while the phonon sector is described by three scalar fields $\phi^I$ and an invariant $b = \sqrt{\det B_{IJ}}$ entering an effective action $F(b)$. At the background level this reproduces a perfect fluid with energy density, pressure and bulk modulus set by two dimensionless parameters $\varepsilon$ and $\kappa$. These parameters fix the enthalpy and bulk modulus in units of the space tension and determine the phonon sound speed through $c_s^2 = \kappa/\varepsilon$. Dissipative effects are modeled by a bulk viscous pressure obeying a Maxwell type relaxation law with a characteristic time scale $\tau(H)$ that depends on the Hubble rate. Motivated by a Boltzmann suppressed scattering rate at a mass-gap scale $H_\star$, we adopt a simple ansatz for $H\tau(H)$ and obtain a compact expression for the effective dark energy equation of state $w_{\rm eff}(H)$. The viscous correction is transient, is most active around $H \sim H_\star$ and drives a temporary phantom deviation. At the same time, we also show that $w_{\rm eff}$ approaches $-1 + \varepsilon$ at very early and very late times. Using a flat $\Lambda$CDM background for $H(z)$, we scan the parameter space and identify a region with $\kappa \ \simeq \ \varepsilon \ \simeq \ 1/3$ and $H_\star/H_0 \ \simeq \ 2.1$. This choice gives an ultralight phonon with sound speed close to the causal limit and a viscous dark energy history $w_{\rm eff}(z)$ that closely tracks a DESI motivated Chevallier-Polarski-Linder parametrization over the redshift range most relevant for the DESI BAO measurements.

We develop an effective mapping for low-redshift photon propagation that captures the leading path-dependent deviations from the standard FLRW redshift. Instead of relying on exact integrations of the Sachs optical equations, we introduce a minimal deformation of the redshift relation z_eff(z) = z minus alpha times f(z), where alpha is a small amplitude and f(z) suppresses the correction at z greater than approximately 0.1. This mapping does not modify the background expansion but encapsulates the leading contribution of inhomogeneous tidal fields to the accumulated Sachs redshift. We derive the implications for the luminosity distance, the low-redshift Hubble relation, and the directional dependence associated with local structure. The framework provides a clean general-relativistic description of path-dependent redshift drift and yields concrete predictions for forthcoming low-redshift surveys.

In this Research Note we present new gravity-darkening exponents ($\beta$) for several stellar evolution models from the ZAMS up to the giant phase. The models were computed using the MESA code (version 7385) for the composition $X = 0.70$ and $Z = 0.02$, adopting A09 opacities and a mixing length parameter of $\alpha_{\rm MLT} = 1.84$. Results were calculated using the triangle strategy, for initial stellar masses of 1.0, 1.5, 2.0, 2.5, 3.0, 5.0, 7.0, 10.0, 15.0, and 20.0 $M_{\odot}$, and provide an update on previously available $\beta$ values that used older models.

Prompt cusps (PCs) form from the direct collapse of overdensities in the early Universe, reside at the center of every dark matter halo, and have density profiles steeper than $r^{-1}$ NFW cusps. Using a suite of high-resolution N-body simulations, we study the evolution of isolated halos in self-interacting dark matter (SIDM) with massive PCs embedded at their centers, a scenario that could be realized in certain SIDM models with light mediators that predict a small-scale suppression of the linear matter power spectrum. We track the evolution of three equally concentrated $10^7\,{\rm{M}}_\odot$ halos, hosting PCs of various total masses, and quantify how the presence of a PC affects the processes of core formation and collapse. Early in the core-formation phase, halos with more prominent PCs exhibit a delayed evolution by a factor of $\sim 2$ due to smaller velocity dispersion gradients in the inner region. During most of the core-collapse phase, the halo evolution becomes closely aligned in physical time, with appropriate rescaling of densities, radii, and velocity dispersions. The scale densities and radii preserve the virial mass of the original halos, but with increased concentration. Deviations occur at the late phase of core-collapse at the level of $\sim 5\%$ relative to the reference collapse track of an NFW halo. These deviations depend non-trivially on both the increased concentration incurred by the PCs, as well as the velocity dispersion (temperature) of the outer halo regions, which can inhibit or enhance the heat transfer process. Our simulations illustrate the complex interplay between the inner and outer halo profiles in determining the onset of core collapse and motivate future studies in the full cosmological context.

The core mass function (CMF) of prestellar cores is essential for understanding the initial conditions of star and cluster formation. However, the universality of the CMF and its relationship to the initial mass function (IMF) remain unclear. We study the CMF in the earliest stage of high-mass star formation using 461 prestellar core candidates and 254 protostellar cores as a part of the ALMA Survey of 70 \mu m Dark High-mass Clumps in Early Stages (ASHES). We find that prestellar core candidates tend to have lower masses than protostellar cores. We also find that the lifetime of prestellar cores is several times longer than the freefall time, although it approaches the freefall time as the core mass increases. The CMF, including both protostellar and prestellar cores, has a power-law slope of $-2.05\pm0.04$, shallower than Salpeter's IMF slope of -2.35. Conversely, the CMF of gravitationally bound, prestellar cores has a steeper slope ($-2.32\pm0.30$), indistinguishable from Salpeter's slope. This finding is consistent with observations in both low-mass star-forming regions and high-mass protoclusters, implying a universal core formation mechanism. The protostellar CMF with a larger maximum core mass can be reproduced by the prestellar CMF when an external gas infall is considered. The inferred mass infall rate is higher than the Bondi-Hoyle-Lyttleton accretion rate and follows a shallower mass dependence (smaller power-law index), more consistent with the tidal-lobe accretion. This may contribute to the evolution of CMFs seen in later stages.

Episodic star formation cycles in both high- and low-redshift galaxies have gained more and more evidence. This paper aims to understand the detailed physical processes behind such behaviors and investigate how such an episodic star-forming scenario can explain the scatter in star-formation rate (SFR) of star-forming main-sequence galaxies. This is achieved through tracing back in time the history of z=0 star-forming central galaxies in the TNG100 simulation over the past 7-8 Gyrs. As the first paper in this series, we provide an overview of the episodic star formation history. We find that two branches of star formation typically develop during each episode: while one branch happens in heavily metal-enriched gas in the centers of galaxies, a secondary branch starts in lower-metallicity regions at galaxy outskirts where fresh gas first arrives, and gradually progresses to inner regions of galaxies. Additionally, the temporal variation in the SFR at galaxy outskirts is more significant than that at centers. As a consequence, the metallicities in both gas and young stars exhibit remarkably different distributions between SFR peaks and valleys. The resulting temporal SFR fluctuation within individual galaxies has an average of ~ 0.2 dex, while the intrinsic differentiation between (the historical mean of) galaxies is ~ 0.15 dex. These two together can well account for the scatter in SFR of ~ 0.25 dex as observed for z=0 star-forming main-sequence galaxies.

Shocks in astrophysical transients are key sites of particle acceleration. If the shock upstream is optically thick, radiation smoothens the velocity discontinuity at the shock (radiation-mediated shocks). However, in mildly magnetized outflows, a collisionless subshock can form, enhancing the efficiency of particle acceleration. We solve the hydrodynamic equations of a steady-state, radiation-mediated shock together with the radiative transfer equations accounting for electron and proton acceleration. Our goal is to explore the impact of the magnetic field and non-thermal radiation on the shock structure and the resulting spectral distribution of photons. To this purpose, we assume a relativistic upstream fluid velocity ($\Gamma_u = 10$) and investigate shock configurations with variable upstream magnetization ($\sigma_u = 0$, $10^{-8}$, $10^{-4}$, $0.1$, and $0.3$). We find that synchrotron self-absorption alters the shock profile for $\sigma_u \gtrsim 10^{-8}$, with resulting changes up to $100\%$ in the bulk Lorentz factor at the shock; for $\sigma_u \gtrsim 0.1$, a prominent subshock forms. The spectral energy distributions of upstream- and downstream-going photons are also altered. Radiative processes linked to accelerated protons are responsible for a high-energy ($\gtrsim 10$ GeV) tail in the photon spectrum; however, the radiation flux and pressure are negligibly affected with consequent minor impact on the shock structure. Our work highlights the importance of coupling the shock hydrodynamics to the transport of photons, electrons, protons, and intermediate particles to forecast the multi-messenger emission from astrophysical transients.

Protoplanetary disks with inner dust cavities (often referred to as "transition disks") are potential signposts of planet formation. We use Gaia astrometry to search for planetary and stellar companions in a sample of 98 transition disks, assessing the occurrence rate of such companions and their potential influence on cavity formation. For the 98 Young Stellar Objects (YSOs), we compute Gaia proper motion anomalies which, together with the RUWE, identify companions with mass ratios $q \gtrsim 0.01$ at $\sim$0.1-30 au. We assess the impact of disk gravity, accretion, disk-scattered light, dippers, starspots, jets, and outflows on the measured proper motion anomalies, concluding that astrometric techniques such as the one of this work can be robustly applied to YSOs. Significant proper motion anomalies are found in 31 transition disks (32% of the sample), indicative of companions. We recover 85% of known companions within our sensitivity range. We model the semi-major axis and mass required for a companion to reproduce the observed astrometric signals. Most inferred companions have $M > 30$ M$\rm{_{J}}$, placing many within or near the stellar mass regime. Seven sources host companions compatible with a planetary mass ($M < 13$ M$\rm{_{J}}$, HD 100453, J04343128+1722201, J16102955-3922144, MHO6, MP Mus, PDS 70, and Sz 76). For the non-detections, we provide the companion masses and semi-major axes that can be excluded in future searches. About half (53%) of detected companions cannot be reconciled with having carved the observed dust cavities. We find that transition disks host as many companions within our sensitivity range as do randomly sampled groups of YSOs and main-sequence stars. If dust cavities are shaped by companions, such companions must reside at larger orbital separations than those of the companions detected here, and we predict them to be of planetary mass. [abridged]

The search for exomoons, or moons in other star systems, has attracted significant interest in recent years, driven both by advancements in detection sensitivity and by the expanding population of known exoplanets. The $\beta$ Pictoris system is a particularly favorable target, as its proximity and directly imaged planets allow for precise astrometric monitoring. We present astrometric constraints on the presence of binary planets and exomoons in the $\beta$ Pictoris system using archival observations from the GRAVITY interferometer and SPHERE instruments. We calculate these limits by modeling the motion of the two orbiting planets and introducing an additional perturbation to the model that simulates the astrometric motion caused by an exomoon orbiting the planet $\beta$ Pictoris b. We find that for short orbital periods ($\approx50$ days), a lunar companion is only allowed if its mass remains below $\approx 180~M_{\oplus}$ ($0.6~M_{\text{Jup}}$) at $3\sigma$ confidence. At intermediate periods near 300 days, we exclude moons more massive than $\approx 65~M_{\oplus}$ ($0.2~M_{\text{Jup}}$) at $3\sigma$ confidence. At longer orbital periods, we place the tightest constraints, ruling out any potential exomoon above $\approx 50~M_{\oplus}$ ($0.15~M_{\text{Jup}}$) at $700$ days and $\approx 30~M_{\oplus}$ ($0.1~M_{\text{Jup}}$) at $1,100$ days (both at $3\sigma$ confidence). These results place the first astrometric constraints on moons and binary planets in the $\beta$ Pictoris system and demonstrate the sensitivity of interferometric observations for exomoon studies.

A galaxy's metallicity and its relation to stellar mass encode the history of gas accretion, star formation, and outflows within cosmic ecosystems. We present new constraints on the low-mass end of the mass-metallicity relation (MZR) at $z\sim2-3$ from ultra-deep JWST/NIRSpec spectroscopy of seven continuum-faint galaxies in the Chemical Evolution Constrained using Ionized Lines in Interstellar Aurorae (CECILIA) Faint sample (Raptis et al. 2025). Our sample includes Ly$\alpha$-selected and other low-luminosity star-forming galaxies with stellar masses $\log(M_\star / M_\odot)\sim7.2-9.7$ and moderately faint rest-UV magnitudes ($-20.7 \lesssim M_{\rm UV} \lesssim -17.3$). Gas-phase oxygen abundances, calculated using empirical calibrations of [O III]/H$\beta$ together with [N II]/H$\alpha$ constraints, span $\sim0.04-0.5$ $Z_\odot$. We measure a steep MZR slope of $\gamma = 0.48 \pm 0.11$, suggesting a rapid increase in metal retention efficiency with mass, consistent with energy-driven outflows. Comparison with lower- and higher-redshift studies indicates an evolution in normalization from $z\sim0$ to $z\sim2$, reflecting less metal enrichment in early galaxies. We find no significant evolution in the MZR between $z\sim2$ and the Epoch of Reionization, suggesting that our galaxies may serve as useful analogs of reionization-era systems. Expanded samples and direct $T_e$-based abundance measurements will be crucial to fully trace the build-up of metals in low-mass galaxies during the peak epoch of cosmic star formation and to test the reliability of strong-line calibrations in these galaxies.

The temperature of baryons at the end of the cosmic dark ages can be inferred from observations of the 21-cm hyperfine transition in neutral hydrogen. Any energy injection from the dark sector can therefore be detected through these measurements. Dark compact objects and dark-matter substructures can modify the baryon temperature by transferring heat via dynamical friction. In this work, we evaluate the prospects for detecting dynamical friction-induced heating from dark compact objects with a mass in the range $10^2 M_{\odot}$ to $10^5 M_{\odot}$, as well as from axion minihalos, using upcoming 21-cm experiments. We find that both the 21-cm global signal and power-spectrum measurements will be sensitive to dark compact objects that constitute about 10% of the dark matter, and will substantially improve our sensitivity to axion-like particles with masses in the range $10^{-18}$ eV to $10^{-9}$ eV.

From the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we identify 14 off-nuclear broad (FWHM>1000 km/s) Halpha and/or Hbeta emission line sources that indicate spatially offset active galactic nuclei (AGN) candidates. In addition to massive black holes (MBHs) in on-going galaxy mergers, this selection can also find MBHs that have been ejected from the host galaxy nucleus due to MBH binary coalescence and asymmetric gravitational wave emission or the dynamical `slingshot' mechanism. Recoiling/slingshot MBHs are predicted to affect co-evolution between MBHs and their host galaxies, and they are observational tracers of past binary MBH mergers and gravitational wave emission. This is the first systematic search through an integral field spectroscopy survey for ejected MBHs to enable uniform constraints on their surface densities. We find that 42% (6/14) have optical image counterparts consistent with galaxy stellar cores from infalling MBHs before the close binary MBH stage. The remaining 58% (8/14) have large broad line luminosities relative to their stellar core mass upper limits (~2 times larger than for central AGN), suggesting merger-driven MBH accretion enhancements or potentially ejected MBHs. The signatures of AGN-ionized narrow emission lines for recoil/slingshot candidates are weaker by 68%, which is consistent with the ejected MBH scenario. The broad line projected velocity offsets range from ~10-600 km/s and suggest motion within the host galaxy potentials. Finally, the implied recoiling MBH surface density upper limit is consistent with predictions that assume random spin orientations in MBH binaries.

We report results of Chandra X-ray observations of CI Tau, a young magnetically active classical T Tauri star for which previous studies have reported periodic variability attributed to a massive planet in a short-period orbit. CI Tau was clearly detected by Chandra in four separate observations acquired in late 2023. The X-ray emission was steady in the first two observations with a characteristic plasma temperature kT ~ 2 keV (~23 MK) and X-ray luminosity log Lx = 29.74 erg/s. During each of the last two observations obtained two weeks later the count rate increased slowly and the X-ray plasma temperature was much higher but remained nearly steady at kT ~ 4 - 5 keV (~46 - 58 MK) and peak luminosity log Lx = 30.5 erg/s. Such variable X-ray emission in T Tauri stars accompanied by high plasma temperatures is a signature of magnetic activity, consistent with the known presence of a strong magnetic field in CI Tau. We summarize the variable X-ray emission properties of CI Tau within the framework of T Tauri stars of similar mid-K spectral type, identify possible variability mechanisms, and assess the effects of stellar X-ray irradiation on the claimed planet.

It has recently been suggested that the typical separation between cores in molecular clouds dominated by turbulence is determined by the sonic scale, the size scale at which the turbulent velocity dispersion equals the sound speed. In this work, we test this hypothesis using a suite of turbulent simulations with Mach numbers $\mathcal{M}=4$ and 8, and three turbulent forcing wavenumbers ($k_{\rm for}=2,4$ and 8). Dense cores are identified through dendrogram analysis of column density maps, and their separations are compared to the sonic scale measured from velocity structure functions. We find no statistical correlation between the core separation and the sonic scale nor with the driving scale. Instead, for each run, the core separation spans the entire range of values between these two scales. Our results indicate that fragmentation in turbulence-dominated clouds is not governed by an intrinsic scale in the turbulent cascade. This finding calls into question the use of the sonic scale as a predictive quantity in star formation theories and cautions against interpreting observational core spacings as evidence for universal turbulent fragmentation physics.

Standard 1D stellar evolution model has poor descriptions of the near-surface layers of stars, and this can be improved by using the atmosphere model computed from 3D hydrodynamical simulations. In this work, we validated the model inferences of the 1D-3D coupling models using 18 well-studied stars from the Kepler LEGACY Sample. We compared our estimates of the fundamental parameters determined with other six pipelines and obtained good consistency. The results indicate that the 1D-3D coupling models can be applied to characterizing solar-like stars with confidence. Our analysis showed similar pattern for the surface term in stars with effective temperature range from ~5000 K to ~6400 K, suggesting that the surface term of the 1D-3D coupling models is simpler and easier to deal with than that of models using classical atmosphere.

Despite the fact that the terrestrial planets formed from the protoplanetary disk, their compositions show marked departures from that of solar nebula condensates. Metallic cores fix oxygen fugacities ($f$O$_2$s) of the planets to 5 (Mercury) to 1 log units below the iron-wüstite (IW) buffer, orders of magnitude higher than the nebular gas. Their oxidised character is coupled with a lack of volatile elements with respect to the solar nebula. Condensates from a solar gas at different temperatures ($T_0$) have Fe/O (by mass) of 0.93 ($T_0$ = 1250 K) to 0.81 ($T_0$ = 400 K), far lower than that of Earth (1.06). Because the reaction Fe(s) + H2O(g) = FeO(s) + H2(g) proceeds <600 K, temperatures at which most moderately volatile elements (MVEs) have condensed, oxidised planets should be volatile-rich, and vice-versa. That this is not observed suggests that planets did not accrete from equilibrium nebular condensates and/or underwent additional volatile depletion/$f$O$_2$ changes. Indeed, MVEs in small telluric bodies (Moon, Vesta) indicate near equilibrium evaporation/condensation at IW-1 and 1400-1800 K. Volatile-depleted elemental yet near-chondritic isotopes of larger telluric bodies (Earth, Mars) reflect mixing of bodies of variable volatile depletion, overprinted by volatile-undepleted material. From the Cr- and Ti isotopes in the BSE, such undepleted matter has been proposed to be CI chondrites. 6% CI added late to an enstatite chondrite-like proto-Earth would match the Earth. However, because Earth is an end-member in isotopic anomalies of heavier elements, no combination of existing meteorites alone can account for its chemical- and isotopic composition. Instead, the Earth is made partially or essentially entirely from an NC-like missing component. If so, the oxidised-, yet volatile-poor nature of inner solar system bodies, including Earth and Mars, is intrinsic to the NC reservoir.

We explore the imprints of fuzzy dark matter (FDM) on the redshifted 21~cm signal from the Cosmic Dawn and the Epoch of Reionization by employing the wavelet scattering transform (WST). FDM, composed of ultralight scalar particles with masses $m_{\mathrm{FDM}} \sim 10^{-22}\,\mathrm{eV}$, exhibits quantum pressure that suppresses the formation of small-scale structures below the de~Broglie wavelength, thereby delaying star formation and modifying the thermal history of the intergalactic medium. Using modified \texttt{21cmFAST} simulations that incorporate both linear and nonlinear effects of FDM on structure formation, we analyze the two-dimensional 21~cm brightness temperature fields through the first- and second-order WST coefficients. The first-order coefficients, $S_1(j)$, quantify scale-dependent variance analogous to the power spectrum, while the normalized second-order ratio $R(j_1,j_2)=S_2/S_1$ captures non-Gaussian cross-scale couplings. We find that low-order couplings, particularly between large and intermediate scales, are highly sensitive to the FDM particle mass and remain robust under SKA1-Low-like thermal noise. Quantitatively, the WST coefficients yield pairwise distances of $\Delta \simeq 225$ between CDM and FDM with $m_{\mathrm{FDM}}=10^{-22}\,\mathrm{eV}$, demonstrating that this framework can effectively discriminate between wave-like and cold dark matter scenarios even under realistic observational conditions. Our results establish the WST as a powerful, noise-tolerant statistical tool for probing the wave nature of dark matter through forthcoming 21~cm observations.

Exoplanet research is essential for understanding planetary formation and the potential for life beyond our solar system. The direct imaging method captures exoplanet light while minimizing light from the host star. This is conventionally achieved with a coronagraph, which allows detailed characterization of planetary atmospheres and features. The Photonic Lantern Nuller (PLN) is an innovative instrument designed for the direct detection of closely orbiting exoplanets within the inner working angle of standard coronagraphs. Unlike traditional coronagraphs, where the planet's signal is usually rotationally invariant, with the same point-spread-function at different position angles, and which also overlaps minimally with the residual stellar signal, data from a PLN consist of a one-dimensional collection of points that do not have rotational invariance and overlap significantly with the residual starlight arising from wavefront errors. Exploiting angular diversity to subtract these stellar residuals with the PLN thus requires adapting the Angular Differential Imaging (ADI) technique for use with non-rotationally invariant planet signals at close separations, where strong self-subtraction effects occur. We reformulate ADI using principal component analysis to develop a method to extract spatial parameters of exoplanets from simulated one-dimensional PLN data. We test two variations of ADI on simulated data and show that injecting an antiplanet signal before stellar estimation helps localize the planet due to self-subtraction at lower separations.

After decades of efforts, optical long-baseline interferometry has become a mainstream observational technique in terms of operation robustness and user friendliness. Interferometry has opened a new observational window, enabling (sub)au-scale resolution of massive stars and direct measurements of orbital parameters, wind structures, and magnetic phenomena. This paper reviews recent advances in interferometric studies of massive stars, focusing on multiplicity, magnetism, and stellar winds.

Weak stability boundary structures have been widely applied to the analysis on ballistic capture and the construction of low-energy transfers. The first step of this application is to compute/identify weak stability boundary structures. Conventional numerical and analytical methods cannot simultaneously achieve computational efficiency and identification precision. In this paper, we propose an efficient and precise method to identify weak stability boundary structures based on deep neural network. The geometric and dynamical properties of weak stability boundary structures are firstly analyzed, which provides further insights into the training of the deep neural network models. Then, the optimal hyperparameter combinations are determined by examining the identification precision of the trained deep neural network models. The performance of the models with the optimal hyperparameter combinations is further validated using the representative test datasets, achieving the precision of 97.26-99.91%. The trained models are also applied to constructing weak stability boundary structures.

Low-energy cosmic rays ($E\lesssim 1$ GeV) are responsible for the ionisation and heating of molecular clouds. While the role of supra-thermal electrons produced in the ionisation process in inducing excitation of the ambient gas (mostly molecular hydrogen) has been studied in detail, the role of primary cosmic-ray nuclei (protons and heavier nuclei) has been generally neglected. Here, we introduce, for the first time, cross sections for proton impact on H$_2$, calculated using the semi-classical implementation of the molecular convergent close-coupling method. Our findings show that proton-induced H$_2$ excitation is comparable in magnitude to that caused by electrons. We discuss the possible implications on the estimate of the cosmic-ray ionisation rate from observations in the near-infrared domain and on the cosmic-ray-induced H$_2$ ultraviolet luminescence. We also derive a new approximated analytical parameterisation of the spectrum of secondary electrons that can be easily incorporated in numerical codes.

We describe how the ESA Comet Interceptor mission, which is due to launch in 2028/29 to a yet-to-be-discovered target, can provide a conceptual basis for a future mission to visit an Interstellar Object. Comet Interceptor will wait in space until a suitable long period comet is discovered, allowing rapid response to perform a fast flyby of an object that will be in the inner Solar System for only a few years; an enhanced version of this concept could realistically provide the first in situ investigation of a visitor from another star system.

Complex organic molecules (COMs) in starless cores provide critical insights into the early stages of star formation and prebiotic chemistry. We present a chemical survey of 16 starless cores (including five prestellar cores) in the Orion A and B molecular clouds, targeting CH3OH, N2H+, CCS, and c-C3HD, using the Atacama Compact Array (ACA) and the Yebes 40-m telescope. CH3OH was detected toward all targets, confirming its ubiquity in starless cores, consistent with previous surveys in Taurus and Perseus. ACA imaging shows that CH3OH, CCS, and c-C3HD generally trace the outer layers of the dense cores outlined by N2H+, each exhibiting distinct spatial distributions. Meanwhile, Comparison with Yebes data reveals an extended, flattened CH3OH component. CCS and c-C3HD tend to be detected or non-detected together across cores, while cores near dust-rich regions on a large scale often lack both, suggesting environmental influences linked to the interstellar radiation field. Within individual cores, CCS typically resides in an outer layer relative to c-C3HD. Our findings underscore the importance of high-resolution studies for understanding the origins and spatial differentiation of COMs and carbon-chain molecules in cold, quiescent environments.

The Hydra I galaxy cluster (Abell 1060) is a nearby example of a low-temperature cluster that exhibits intermediate cool core and non-cool core properties. We aim to extend the characterization of the intracluster medium (ICM) properties at least until $R_{200}$ and study the correlation between the X-ray emission and non-thermal emission within $R=0.15$$R_{500}$, and optical/IR galaxy distribution beyond $R_{200}$. We used data from the first four SRG/eROSITA All-Sky Surveys and an archival Chandra observation to image the X-ray emission from Abell 1060. We also used multiwavelength data from TGSS (radio), 2MASS (IR), and NED (optical) to investigate the non-thermal emission, 2D galaxy distribution, and its redshift evolution, respectively. The surface brightness and spectral analyses are also extended until 3$R_{200}$ and $R_{200}$, respectively, following a detailed cosmic X-ray background (CXB) analysis. Our fully corrected eROSITA image showcases a relaxed ICM morphology within $R_{500}$. We detect two weak ICM shocks with Mach number $M\approx 1.5$ near the central galaxy NGC 3311 that coincide with diffuse radio emission along the line of sight. Furthermore, we detect multiple soft X-ray excesses with high spatial correlation with the 2D optical galaxy distribution beyond $R_{200}$. In particular, the excess in the north has a significance of $5.8\sigma$ above the local CXB level. This suggests that Abell 1060's outskirts are actively accreting baryons. We also estimate the average ICM temperature and metallicity of $\langle k_\mathrm{B}T \rangle=2.27\substack{+0.15\\-0.11}\thinspace$ keV and $\langle Z\rangle=0.19\substack{+0.05\\-0.03}\thinspace Z_\odot$, respectively, from the 0.2-0.5$R_{500}$ annulus. Overall, the temperature profile is broadly consistent with the average temperature profiles from hydrodynamical simulation and Suzaku between 0.28$R_{200}$ and $R_{200}$.

Planetary materials show systematic variations in their nucleosynthetic isotope compositions that resonate with orbital distance. The origin of this pattern remains debated, limiting how these isotopic signatures can be used to trace the precursors of terrestrial planets. Here we test the hypothesis that interstellar ices carried supernova-produced nuclides by searching for a supernova nucleosynthetic fingerprint in aqueous alteration minerals from carbonaceous and non-carbonaceous chondrite meteorites. We focus on zirconium, a refractory element that includes the neutron-rich isotope $^{96}$Zr formed in core-collapse supernovae. Leaching experiments reveal extreme $^{96}$Zr enrichments in alteration minerals, showing that they incorporated supernova material hosted in interstellar ices. We show that the Solar System's zirconium isotope variability reflects mixing between these ices and an ice-free rocky component. Finally, the presence of supernova nuclides in a volatile carrier supports models where the Solar System's nucleosynthetic variability was imparted by thermal processing of material in the protoplanetary disk and during planetary accretion.

The recent data release DR2 from the Dark Energy Spectroscopic Instrument (DESI) has reinforced compelling evidence supporting the dynamical nature of dark energy. In this respect, we here explore a two-scalar field cosmological model, dubbed Chiral-Quintom paradigm, that allows for a unified description of early- and late-time cosmic accelerations, namely inflation and dark energy, respectively. Moreover, we show that it provides a mechanism to cross the phantom divide without instabilities. To do so, we first focus on scenarios where the scalar fields evolve under an exponential potential, leading to distinct cosmological behaviors, including tracking solutions and slow-roll or hyperbolic inflation. Afterwards, by considering a nonlinear potential and a mixing in the kinetic sector, we show that the model can also describe the matter-dominated era, offering a potential unification of the dark sector. Accordingly, we place bounds on this double-field paradigm by employing recent low-redshift observational data, including supernovae from the Pantheon+ and Union3.0 compilations, observational Hubble parameter measurements, and baryon acoustic oscillation data from the latest DESI release. By means of the Akaike information criterion, the standard $\Lambda$CDM scenario and the $w_0w_a$CDM model are thus compared with the Chiral-Quintom approach, showing that, in principle, multiple-field dynamics does not seem to be ruled out, ultimately providing a flexible framework for describing late-time dynamics.

We present the high-resolution X-ray spectrum of the Z-source Cygnus X-2, obtained with X-Ray Imaging and Spectroscopy Mission (XRISM). The observations have enabled a precise characterization of the Fe-K emission lines from the accretion disk corona (ADC) and a possible detection of an ultra-fast outflow (UFO). The ADC component has at least two distinct regions. The Fe-K emission lines are remarkably broad, with a velocity dispersion of $\sigma_v=2300^{+900}_{-700}$~km~s$^{-1}$. This can be larger than what was observed in previous Chandra observations ($1120\pm870$~km~s$^{-1}$) and recent XRISM observations of the similar Z-source GX 340+0 (360--800~km~s$^{-1}$). Furthermore, we marginally detect a blueshifted absorption feature, which we identify as either \ion{Fe}{26} or \ion{Fe}{25}, with outflow velocities of $-0.29c$ or $-0.26c$, respectively. The mass loss rate is $\dot{M}\sim10^{-8}M_\odot\,\mathrm{yr}^{-1}$. Our findings suggest that the ADC in Cyg X-2 has a kinematically active environment with a high degree of turbulence and kinetic energy.

The Hilda population occupies the stable 3:2 mean-motion resonance of Jupiter and provides a window into Solar System evolution, including collisional processes. The NSF-DOE Vera C. Rubin Observatory will conduct the ten-year Legacy Survey of Space and Time (LSST). We present a simulation of Rubin's discovery of Hildas with the Sorcha (Merritt et al. 2025; Holman et al. 2025) survey simulator and the recovery of their light curves. We constructed a synthetic Hilda population model which includes distributions of orbital properties, sizes, collisional families and colors. We included two color classes corresponding to the Jupiter Trojan populations (Wong & Brown 2017). We applied three distinct populations of sinusoidal light-curves to this same orbit-size-color model: (1) a Gaussian kernel density estimate (KDE) fit to rotational periods and amplitudes from the Lightcurve Database (LCDB; Warner et al. 2009) (2) a super-fast rotator (SFR) population (0-3 hours) and (3) a super-slow rotator (SSR) population (100-1400 hours). Over the ten-year simulated survey, we predict LSST will discover ~33,400 Hildas, a fivefold increase over the known population. Using a multiband Lomb-Scargle Periodogram via Astropy (Price-Whelan et al. 2022) we confidently recover ~46.5% of Hildas in our LCDB-based population, higher than typical in observational searches. This suggests our light-curve population model may differ from the intrinsic population. We find strong biases in light-curve amplitude, with recovery efficiency dropping sharply below 0.1 magnitudes, while biases from rotational period are comparatively weak aside from cadence-related features such as LSST's ~36 minute revisit cadence. ...

The southern Galactic plane has been mapped at optical wavelengths and at under one-arcsecond angular resolution by the VST Photometric Ha Survey of the Galactic plane and bulge (VPHAS+). Anticipating the release of a uniform photometric calibration of the entire survey, we examine the properties of VPHAS+ ugriHa photometry of r < 19 mag. point sources in the third Galactic quadrant (longitudes 210^o to 260^o). We compare our interim calibration in gri with that of Pan-STARRS, the DECam Plane Survey (DECaPS-2) and Skymapper. We use the comparisons to identify small gri photometric offsets. Corrections to the Ha and u magnitude scales are determined via comparison with synthetic photometry. VPHAS+ and its northern counterpart, the INT Galactic Plane Survey (IGAPS), are shown to closely align, where they overlap across the celestial equator. Aided by Gaia Data Release 3, SIMBAD, and specialist catalogues, we present selections of: A stars; sub-luminous stars; intrinsically-red luminous stars; young stellar objects; emission-line and OB stars. Attention is drawn to stellar variability as a contaminant in selecting emission line objects via (r - Ha) excess. It is argued the (r - i, r - Ha) plane is the better choice for this selection than (g - i, r - Ha). Using A stars to map extinction, we trace the main run of dust obscuration, situated at mainly negative Galactic latitudes. Like the dust, OB and emission line stars are more frequent below the Galactic equator: at heliocentric distances of up to ~7 kpc, these stars' distribution fit in with the known warping of the Galactic plane. An overdensity of B stars, several degrees across and potentially in the Outer Arm, is found around (l,b) = (212.0, -0.6).

We summarize recent results of the observations of high (1 TeV-100 PeV) and ultrahigh ($\geq 100$ PeV) energy neutrinos, including the detection of a diffuse cosmic high-energy neutrino background, the identification of the first neutrino source candidates, and the observation of high-energy neutrinos from the Galactic plane. These findings open a new window to the universe by enabling the use of neutrinos to probe the cosmos that are otherwise inaccessible via photons. Although the origins of most detected neutrinos remain uncertain, we highlight several distinctive features of their sources that have emerged from current observations.

Gravitational lensing of background galaxies by intervening matter is a powerful probe of the cosmological model. In the era of Stage IV surveys, contamination from galaxies below the detection threshold has emerged as a significant source of bias. Adopting a noise-bias-free machine-learning method to estimate shear, we quantify the impact of faint galaxies for a Euclid-like survey. In our baseline simulations, faint blends induce a multiplicative shear bias of -0.008, well above Euclid's requirement. Similar to previous studies, we find that calibration simulations must include neighbouring galaxies to AB apparent magnitudes as faint as 27.0 (+2.1, -0.9) and within approximately 1.0 (+0.2, -0.2) arcsec of each bright sample galaxy (BSG; the galaxy for which shear is measured). By varying faint galaxy properties, we identify which ones significantly affect shear biases and quantify how well they must be constrained. Crucially, we find that biases not only depend on the mean projected faint-galaxy density and apparent-magnitude distribution across the sample, but also on how these quantities vary with the observed brightness of the BSG. Furthermore, biases are sensitive to radial and tangential alignments and positional anisotropy of faint galaxies relative to BSGs. By contrast, shear coherence between BSGs and faint galaxies, parallel orientation alignments, and variations in the faint galaxy size-magnitude relation have negligible impact within the parameter ranges explored. Our results guide calibration simulations and highlight the critical role of deep observations in measuring the properties of faint galaxies.

We present fast cadence and high resolution observations of flare ribbons from the Solar Orbiter Extreme Ultraviolet Imager (EUI). Utilizing the short-exposure observations from the EUI High Resolution Imager in EUV (HRIEUV), we find small-scale blob/bead-like kernel structures propagating within a hook at the end of a flare ribbon, during the impulsive phase of a C9.9-class solar flare. These bead structures are dynamic, with well-resolved spatial separations as low as ~420-840 kilometers (3-6 pixels) - below the observable limit of full-disk solar imagers. We analyze the evolution of the plane-of-sky apparent velocity and separation of the flare ribbon structures, finding evidence for multiple processes occurring simultaneously within the flare ribbon. These processes include - quasi-periodic pulsation (QPP)-like brightenings, slow back-and-forth zig-zag motions along the ribbon, rapid apparent motions along the ribbon (600+ km/s), and stationary blob-like structures. Finally, we conduct Fast Fourier Transform analysis and analyze the start times of exponential growth in the power spectrum at different spatial scales across the flare ribbon. Our analysis reveals that the ribbon beads form with a key spatial separation of 1.7-1.9 Mm, before developing into more complex structures at progressively larger and smaller spatial scales. This observation is consistent with predictions of the tearing mode instability.

We present the measurement of 26 elemental abundances of SMSS J022423.27$-$573705.1 (SMSS 0224$-$5737), an extremely metal-poor (EMP) star with a weak $r$-process signature. We report the measurements of N, O, V, Zn, and Ba, and the upper limits for Mo, Ru, Pd, Ag, and Eu for the first time. SMSS 0224$-$5737 exhibits low C abundance and high N and O abundances suggesting that C is converted to N by the enhanced mixing during the evolution. The abundance pattern up to the Fe-peak elements is generally in good agreement with the average abundance of EMP stars, although a notable feature is the high [Zn/Fe] ratio ([Zn/Fe] $= +0.88$). We confirm the enhancement of the first-peak neutron-capture elements (Sr, Y, and Zr) and determine a low Ba abundance [Ba/H] $= -5.25$, that is, [Ba/Fe] $= -1.45$. The extremely high ratio of [Zr/Ba] $= +2.60$ makes SMSS 0224$-$5737 the EMP star with the most pronounced weak $r$-process signature observed to date. The abundance pattern of the neutron-capture elements is compared with the yields from $r$-process nucleosynthesis models. The sharp decline in abundances beyond Zr disfavors neutron star merger or electron-capture supernova models, but are reproduced either by proto-neutron star wind models or by magneto-rotational supernova models. Considering the high [Zn/Fe] ratio, a magneto-rotational supernova is the most plausible origin of SMSS 0224$-$5737. This study demonstrates that the abundance measurements of both light and neutron-capture elements, even at low abundances, are crucial for unveiling the astrophysical sites of the weak $r$-process.

We determine the $z=6-8$ ultraviolet (UV) LF in the JWST UNCOVER field behind Abell 2744, through a depth-tied completeness model and source-plane selection taking multiple images into account. We compute the intrinsic $M_{UV}$ LFs and lens-dependent effective volumes $V_{eff,\ell}(M_{UV})$ for two of them (CATS, GLAFIC), and construct binned LFs with statistical (Gehrels) errors presented separately from the lens-model spread. When fitting each lens model independently with both the original and turnover-extended form (while keeping $M^*=-20.5$), we obtain decisive model selection in favour of a faint-end turnover. The peak of the turnover posteriors is at $\alpha \simeq -1$, $\beta \simeq 1.6$, and drive $M_T$ to the high side of our prior ($\sim -18$) and the difference in normalization between CATS-GLAFIC is absorbed into $\phi^*$. A similar magnification-bias test in the image plane gives a $B(\mu)$ trend matching to what is expected, which provides an internal consistency of our selection and volumes. Our results indicate an early suppression of the UV LF at $z=6-8$, and carry implications for the ionizing contribution of ultra-faint systems during reionization.

We present the results of the multi-epoch broadband spectral study of 1ES 2344+514 and study the evolution of physical parameters. We used nearly simultaneous data obtained from 2017 June 6 to 2022 August 6 (MJD 57910 -- 59797) in optical, UV, X-ray and $\gamma$-ray wavebands from various instruments including Swift-UVOT, Swift-XRT, NuSTAR, AstroSat (SXT and LAXPC), Fermi-LAT, and TeV flux from MAGIC. During 2017 July, 1ES 2344+514 appeared to be in the highest flaring state compared to other periods. We used the 0.5 -- 7.0 keV and 3.0 -- 20.0 keV data, respectively, from SXT and LAXPC of AstroSat and also 0.3 -- 8.0 keV and 3.0 -- 79.0 keV data, respectively, from Swift-XRT and NuSTAR. A joint fit between SXT and LAXPC, and between Swift-XRT and NuSTAR has been done for constraining the synchrotron peak. A clear shift in the synchrotron peak has been observed, which included 1ES 2344+514 in the HSP BL Lac family. A `harder-when-brighter' trend is observed in X-rays, and the opposite trend, i.e., `softer-when-brighter', is seen in the $\gamma$-rays. The multi-epoch broad-band spectral energy distributions (SEDs) of this source were built and studied to get an idea of the radiation processes. The SEDs were fitted using a steady-state leptonic one-zone synchrotron+SSC model, and the fitted parameters of the emission region are consistent with those of other TeV BL Lacs. In this study, we found a weak correlation tendency between bolometric luminosity and magnetic field (B), as well as between bolometric luminosity and the break Lorentz factor ($\gamma_{break}$).

We test whether the Efron-Petrosian (E-P) method can recover the inverse-square law dependence of the radio pulsar flux, using a synthetic catalog generated according to the specifications of the Parkes multi-beam survey using the {\tt PsrPopPy} software. We find that the E-P method cannot reproduce the inverse-square law, except over a narrow range of flux thresholds and even here we don't get pristine agreement. The main reason for the deviation is that the synthetic radio pulsar catalog is truncated based on a cut on the pulsar signal to noise ratio (SNR), which has a non-linear dependence on the flux along with plenty of scatter. We show that the disagreement is exacerbated as we raise the SNR threshold. We then demonstrate that if we create a synthetic catalog based on a flux cut (instead of an SNR-based threshold), we can recover the true distance exponent, with an accuracy ranging from pristine agreement to within $\pm 1 \sigma$ depending on the chosen flux threshold.

An asteroseismic analysis has revealed a magnetic field in the deep interior of a slowly-rotating main-sequence F star KIC9244992, which was observed by the Kepler spacecraft for four years. The star shows clear asymmetry of frequency splittings of high-order dipolar gravity modes, which cannot be explained by rotation alone, but are fully consistent with a model with rotation, a magnetic field and a discontinuous structure (glitch). Careful examination of the frequency dependence of the asymmetry allows us to put constraints on not only the radial component of the magnetic field, but also its azimuthal (toroidal) component. The lower bounds of the root-mean-squares of the radial and azimuthal components in the radiative region within 50 per cent in radius, which have the highest sensitivity in the layers just outside the convective core with a steep gradient of chemical compositions, are estimated to be $\mathsf{B}_{\mathrm{r}}^{\min}=3.5 \pm 0.1$ kG and $\mathsf{B}_{\phi}^{\min}=92 \pm 7$ kG, respectively. The much stronger azimuthal component than the radial one is consistent with the significant contribution of the differential rotation although the star has almost uniform rotation at present. The estimated field strengths are too strong to be explained by dynamo mechanisms in the radiative zone associated with the magnetic Tayler instability. The aspherical glitch is found to be located in the innermost radiative layers where there is a steep gradient of chemical composition. The first detection of magnetic fields in the deep interior of a main-sequence star sheds new light on the problem of stellar magnetism, for which there remain many uncertainties.

Green Bean is a rare type of galaxy which represents a short-lived phase in the life cycle of active galactic nuclei (AGN), characterised by large-scale, powerful ionised clouds in the circumgalactic medium. Recent studies demonstrate that these extended ionised structures may reflect fading signatures of past AGN activity, often manifested in the form of large-scale ionisation cones. The analysis of their observational properties provides unique constraints on AGN lifetimes, feedback mechanisms, and transitions between radiative and kinetic modes of activity. In this paper we announce the first results of the project dedicated to the long-slit spectroscopic and scanning Fabry-Perot interferometric observations of Green Bean galaxies at the Russian 6-m telescope with SCORPIO-2 multi-mode instrument. We describe the data reduction and spectral fitting procedures that allow one to characterise ionisation conditions in extended gaseous regions of the galaxy SDSSJ095100.54+051026.7.

We study the multiplicity of galaxies in the largest sample of red-Herschel sources ($S_{250 \mu \mathrm{m}} < S_{350 \mu \mathrm{m}} < S_{500 \mu \mathrm{m}}$) using archival ALMA observations. Out of 2416 fields with ALMA detections (from a total of 3,089 analyzed maps), we identify 474 multiple systems within a radius of 16 arcsec (equivalent to the 500 $\mu$m Herschel beam-size): 420 doubles, 51 triples, and 3 quadruples. In each case the brightest source contributes, on average, 64, 48, and 42 per cent of the total flux in double, triple, and quadruple systems. The average combined ALMA flux density of the sources in double systems is comparable to that of the two brightest components within triple and quadruple systems. Non-parametric tests suggest that only a small fraction of the double systems ($\lesssim13$ per cent) are comprised of sources with compatible redshifts, while 47-67 per cent of triple and quadruple fields contain at least one potentially associated pair. Simulations using a mock catalogue of dusty star-forming galaxies suggest that 32 per cent of the double systems are likely physically associated ($\Delta z < 0.01$, i.e. $\lesssim$10 comoving Mpc at $z = 3$) and, while only 8 per cent of the triple and none of the quadruple systems meet this criterion, $\sim$ 70 per cent of them include at least one likely associated pair. Our results suggest that enhanced star formation rates in submillimetre galaxies are primarily driven by internal processes rather than large-scale interactions. This study also provides a catalog of potential overdensities for follow-up observations, offering insights into proto-cluster formation and evolution.

We present an empirical investigation into the relationship between galactic age and dark matter content across a sample of 16 nearby, well-resolved spiral galaxies. Using raw rotation curve data from IOA Tokyo's publicly available repository, we model each galaxy's mass distribution via a three-component decomposition (Hernquist bulge, exponential disk, and a Navarro-Frenk-White (NFW) dark matter halo) fit using Monte Carlo simulations. The onset of dark matter dominance was identified using the NFW scale radius, beyond which we computed the total enclosed mass via Keplerian dynamics. I-band luminosities for these regions were estimated using a calibrated Tully-Fisher relation, yielding precise mass-to-light (M/L) ratios. We further calculated dark matter mass and density using NFW profile equations, and galaxy ages were retrieved through an extensive literature survey of stellar population studies. Our analysis reveals strong positive correlations between galactic age and both dark matter mass (Pearson $r \approx 0.91$; Spearman $\rho \approx 0.93$) and density (Pearson $r \approx 0.91$; Spearman $\rho \approx 0.91$), as well as M/L ratios, suggesting a robust link between evolutionary history and dark matter build-up. These findings are in quantitative agreement with predictions from large-scale cosmological simulations that incorporate assembly bias and smooth accretion-dominated growth, reinforcing the view that older galaxies, having formed earlier in high-density peaks, have accumulated significantly more dark matter over cosmic time. Our results offer observational evidence for time-dependent dark matter assembly and establish galactic age as a meaningful tracer of halo evolution.

We investigated the shadow and optical appearance of Bardeen black hole (BH) immersed in perfect-fluid dark matter (PFDM). Using the EHT data, we find that the DM parameter is restricted to a narrow allowed range, confined to values of $\mathcal{O}(10^{-1}-10^{-2})$. Additionally, we observed that increasing the DM parameter substantially enlarges the photon sphere, the critical impact parameter, and thus the shadow radius, whereas the magnetic charge $g$ produces only negligible corrections. The DM component also modifies the optical appearance of the accretion disk: higher $b$ systematically suppresses the observed radiation intensity and reduces image brightness, while changes in $g$ yield only marginal effects. Subsequently, using the Novikov-Thorne thin disk model, we analyze primary and secondary images, observed flux, and redshift patterns, showing that PFDM noticeably alters image size and brightness whereas magnetic charge has a negligible influence. No blueshifted regions appear in any configuration. These results highlight several promising observational signatures for testing DM environments and nonlinear electrodynamic effects around BHs.

Based on the Zwicky Transient Facility (ZTF), we selected 10 blazars as our sample sources. Among these, we found four blazars (J 0923.5+4125, J 1221.3+3010, J 1503.5+4759, and J 1652.7+4024) showing possible indications of quasi periodic oscillations (QPOs) modulation. We conducted a detailed analysis of their optical light curves (g- and r-bands) over the past five years using the root mean square (RMS)-Flux relation, flux distribution, and QPO detection methods to investigate their variability characteristics. A linear RMS-Flux relation is present in both bands, and their flux distributions follow a log-normal form. This suggests that optical variability may arise from multiplicative, nonlinear processes across different timescales and flux states. Further QPO analysis using the weighted wavelet Z-transform (WWZ), Lomb-Scargle periodogram (LSP), and autoregressive integrated moving average (ARIMA) methods identified candidate periodic signals in four blazars. J 0923.5+4125 (period $\sim$ 205 days) and J 1221.3+3010 ($\sim$ 630 days) show local significances of $\sim 3 \sigma$, whereas J 1503.5+4759 ($\sim$ 38.5 days) and J 1652.7+4024 ($\sim$ 48 days) reach $\sim 4 \sigma$. After accounting for the look-elsewhere effect, the global significances for J 1503.5+4759 in the g- and r-bands are $\sim 2.7 \sigma$, while for J 1652.7+4024 they are approximately $\sim 2.5 \sigma$ in both bands. These two blazars warrant further monitoring and investigation.

We compute the thermal conductivity and thermoelectric power (thermopower) of the inner crust of compact stars across a broad temperature-density domain relevant for proto-neutron stars, binary neutron-star mergers, and accreting neutron stars. The analysis covers the transition from a semi-degenerate to a highly degenerate electron gas and assumes temperatures above the melting threshold of the nuclear lattice, such that nuclei form a liquid. The transport coefficients are obtained by solving the Boltzmann kinetic equation in the relaxation-time approximation, fully incorporating the anisotropies generated by non-quantizing magnetic fields. Electron scattering rates include (i) dynamical screening of the electron-ion interaction in the hard-thermal-loop approximation of QED, (ii) ion-ion correlations within a one-component plasma, and (iii) finite nuclear-size effects. As an additional refinement, we evaluate electron-neutron scattering induced by the coupling of electrons to the anomalous magnetic moment of free neutrons; this contribution is found to be subdominant throughout the parameter range explored. To assess the sensitivity of transport coefficients to the underlying microphysics, we perform calculations for several inner-crust compositions obtained from different nuclear interactions and many-body methods. Across most of the crust, variations in relaxation times and in the components of the anisotropic thermal-conductivity and thermopower tensors reach up to factors 3-4 and 1.5-2, respectively, with the exception of the region where pasta phases are expected. These results provide updated, composition-dependent microphysical inputs for dissipative magneto-hydrodynamic simulations of warm neutron stars and post-merger remnants, where anisotropic heat and charge transport are of critical importance.

The search for potentially habitable exoplanets is a primary objective in modern astrophysics, yet the vast number of candidates discovered by missions like Kepler and TESS presents a significant challenge for detailed follow-up characterization. An efficient and reliable method for prioritizing the most promising targets is therefore essential. In this paper, we propose a novel first-order filter for identifying potentially habitable worlds based on a simple geometric ratio: the orbital semi-major axis to the stellar diameter ($d/D_s$). Using data from the NASA Exoplanet Archive, we demonstrate that the ideal value for this ratio is not constant, but is dependent on the host star's spectral class. We establish a tiered framework of ideal ratios, beginning with $\approx 108$ for G-type stars (anchored by the Earth-Sun system), and decreasing by a factor of two for K-type ($\approx 54$) and M-type ($\approx 27$) stars, respectively. Our analysis reveals a strong correlation, showing that exoplanets whose $d/D_s$ ratios are close to these empirically derived values consistently exhibit high Earth Similarity Index (ESI) scores. We propose that these tiered ratios represent "Habitability Main Sequences," analogous to the Hertzsprung-Russell diagram for stars, providing a valuable and straightforward tool for the astronomical community to rapidly screen large datasets and efficiently shortlist high-priority candidates for further investigation with next-generation observatories.

The rapid expansion of satellite constellations is transforming the radio-frequency environment around the Earth. At the same time, radio astronomy is entering a new era of sensitivity and survey capability, requiring unprecedented control of interference. This primer introduces satellite operators, engineers, spectrum managers and policy makers to the basic concepts of radio astronomy, explains why the discipline is uniquely vulnerable to interference, and outlines the regulatory and practical tools available to manage coexistence.

A future measurement of Lense-Thirring (LT) precession using a binary radio pulsar is expected to yield the pulsar's moment of inertia ($I_{\rm p}$). However, most of the known pulsar-binary systems expected to provide this opportunity will exhibit linear variations in the orbital elements due to LT precession that are difficult to separate from variations induced by other mechanisms. In this work, we demonstrate that the pulsar-timing signature of LT precession for an arbitrary orbital orientation produces nonlinear orbital variations; if detected, these nonlinear variations provide a means to constrain $I_{\rm p}$ without the need for timing-independent measurements of orbital geometry or distance. We show through simulations that these signatures are indeed detectable in pulsar-binary systems with an appropriate combination of timing precision and orbital geometry. Our simulations also show that nonlinear orbital variations from LT precession are expected to be detectable in PSR J1757$-$1854 after only 15 yr of dedicated timing.

The high-precision photometry from NASA's Kepler and TESS missions has revolutionized exoplanet detection, enabling the discovery of over 5500 confirmed exoplanets via the transit method and around 10000 additional candidates awaiting validation. However, confirming these candidates as true planets demands meticulous vetting and follow-up observations, which hampers the discovery of exoplanets in large-scale datasets. To address this challenge, we developed a machine learning framework trained on Kepler's catalog of confirmed exoplanets and false positives to accurately identify true planetary candidates. Our model uses transit properties, planetary characteristics, and host stellar parameters as training features. The optimized model achieved 83.9% accuracy in cross-validation. When applied to 3987 TESS candidates with complete observational data, the model identified 1595 new high-confidence planets and correctly recovered 86% (358/418) of all previously confirmed TESS exoplanets in a blinded validation test. Our analysis revealed 100 previously unrecognized multi-planet systems, including five systems--that host habitable-zone exoplanets. Additionally, we identified 15 more planets within the habitable zone of a single system, suggesting strong potential for liquid water stability under conservative planetary albedo assumptions. This work demonstrates that machine learning can accelerate exoplanet validation while maintaining scientific rigor. Our modular design enables direct adaptation to future photometric missions like PLATO or Earth 2.0.

We study the influence of positrons on the outer crusts of neutron stars and the interiors of white dwarfs, introducing them as a novel component in both the composition of matter and in transport processes. We solve a system of coupled Boltzmann kinetic equations for the electron and positron distribution functions in the relaxation-time approximation, taking into account electron-ion, positron-ion, and electron-positron collisions. The relevant scattering matrix elements are calculated from one-plasmon exchange diagrams, with in-medium polarization tensors derived within hard-thermal-loop effective theory. Numerical results are obtained for matter composed of carbon nuclei. We find that the conductivity rises with temperature, following a power law sigma proportional to the 4th power of T in the semi-degenerate regime and sigma proportional to T in the nondegenerate regime, due to the intense creation of thermal electron-positron pairs and the resulting collisions among them. These results highlight the importance of including positrons in the transport properties of heated, dense astrophysical plasmas.

Sunspots and solar flares are two different manifestations of magnetic activity on the surface of the Sun. On the Sun, flares typically occur close to spots. In this paper we test this the connection between spots and flares on other stars. We detect 218,386 stellar flares on 14,163 spotted stars using a new algorithm called \textsc{toffee}. Inhomogeneous spot distributions mean that as stars rotate they become brighter when less spots are facing the observer, and dimmer when more spots are facing the observer. We determine that flares occur when the star is brighter $49.97\pm 0.21\%$ of the time, i.e. there is an equal preference for the flares to occur when the star is relatively bright or dim. We therefore find no evidence for a correlation between flare rate and spot occurrence, contrary to what is seen on the Sun.

1H 0323+342 is the nearest gamma-ray narrow-line Seyfert 1 galaxy (z=0.063). Its X-ray spectrum (0.3-10 keV) is characterised by significant spectral variability observed by many authors, with a backbone with photon index ~2 occasionally superimposed by a hard tail. This spectral variability has been interpreted as the interplay between the X-ray corona and the relativistic jet. The X-ray fluxes in the 0.3-10 keV energy band are generally around ~10^-11 erg cm^-2 s^-1, making it easier to get sufficient statistics even with short exposures. Here I present a reanalysis of all the available X-ray observations with Swift (181 obs), XMM-Newton (7 obs), Chandra (1 obs), and Suzaku (2 obs) performed between 2006 and 2025. Possible interpretations are proposed and discussed.

We present narrowband observations of the Fe XIV (530.3 nm), Fe X (637.4 nm), and Fe XI (789.2 nm) coronal emission lines from the 2023 April 20 Total Solar Eclipse in Australia. We deployed pairs of telescopes for each emission line that were equipped with narrowband filters centered on, and several nanometers away from, the center wavelengths of the lines. The secondary continuum telescopes were used to measure and remove the combined continuum K- (electron) and F- (dust) corona, which dominate coronal emission at optical and infrared wavelengths. Significant emission was detected from all three lines from 1.03 solar radii (Rs) continuously outward to at least 6 Rs. The brightness of the lines and continuum are absolutely calibrated to the solar disk, and are validated by a comparison with LASCO-C2 observations made at the same time. Using these observations, we inferred the line emission ratios resolved throughout the middle-corona (defined as 1.5 to 6 Rs) for the first time. These line ratios are a probe of the electron temperature, which have important implications for constraining models of coronal heating and the characterization of solar wind formation, yet these emission lines have scarcely been quantified beyond 3 Rs in the corona. This study demonstrates the enduring potential of eclipse observations for coronal physics and suggests that future spacecraft missions could observe these lines farther out than has been attempted previously.

White dwarfs (WDs) are the final fate of about 97\% of the stars in our galaxy, making them vital tracers of stellar history. A fraction of WDs exist in cataclysmic variable (CV) systems, accreting matter from a nearby companion star. A subset of CVs undergo episodic rapid mass transfer, termed dwarf novae (DNe) outbursts. Some accreting WDs exhibit near sinusoidal photometric variations, interpreted as $g$ mode pulsations. However, identifying pulsation modes in accreting WDs remains challenging due to the paucity of available observed modes. In this work, we present a comprehensive computation of the observable $g$ mode frequencies across a range of WD parameters, varying the WD mass, size of the newly accreted layer and core temperature. We also introduce a novel method for mode identification based on the time evolution of pulsation periods following an accretion episode. Our mode identification method does not rely on the direct detection of the consecutive radial mode orders, frequently required in isolated WDs. Moreover, this work improves upon our previous WD modeling efforts. We use a more realistic core temperature in addition to thermohaline mixing and element diffusion enabled during the accretion phase.

We present a generalized elliptical cylindrical flux rope model for interplanetary coronal mass ejections (ICMEs) that allows for a non zero poloidal component in the internal magnetic field. We introduce a two step reconstruction algorithm that decouples the geometric configuration from the magnetic field fitting in order to improve numerical stability and physical consistency. Applied to Parker Solar Probe data, the new radial poloidal model preserves the global flux rope geometry while achieving a substantially better fit to the internal magnetic field than the traditional radial model, offering a more accurate and realistic description of ICME structures.

On April 17, 2024, the third successful Hi-C sounding rocket flight, Hi-C Flare, recorded coronal images in Fe XXI 129 A emission from 11 MK plasma during the post-maximum phase of an M1.6-class solar flare, achieving unprecedented spatial (~300 km) and temporal (1.3 s) resolutions. The flare started at 21:55 UT, peaked at 22:08 UT, and lasted ~40 minutes. Hi-C observed for over five minutes (22:15:45 to 22:21:25), starting roughly eight minutes after flare maximum. A sudden compact bright burst - 875 +/- 25 km wide, lasting 90 +/- 1.3 s, displaying a plane-of-sky motion of ~50 km/s toward the loop apex, and splitting into two toward the end - occurs near the foot of some post-flare loops. Its size and brightness are reminiscent of flare-ribbon kernels during a flare's rapid rise phase, kernels marking sites of sudden heating and hot plasma upflow, making its occurrence during the late phase surprising. Such isolated brightenings in a flare's post-maximum phase are rare, and have not been previously reported. The kernel was detected in all SDO/AIA channels. Its 1600 A light curve peaked ~50 s earlier than its 131 A light curve, similar to that of flare-ribbon kernels, albeit with a smaller delay of ~25 s, during the impulsive phase of the flare. In SDO/HMI magnetograms, the kernel sits in unipolar positive magnetic flux near an embedded clump of negative flux. Although localized magnetic reconnection within the kernel (a microflare) cannot be ruled out for its cause, the observations favor the localized brightening being an isolated, exceptionally late flare-ribbon kernel, resulting from an exceptionally late burst of the flare's coronal reconnection.

We present an investigation into the fate of disrupting near-Sun comet C/2019 Y4 (ATLAS). Imaging observations with the Lowell Discovery Telescope (LDT), obtained 3--5 months after the reported disruption and the last sighting, constrained the primary component C/2019 Y4-B to an absolute magnitude of $H>20.5$ or a diameter of $D<0.5$~km (assuming a geometric albedo of 0.04) at an outbound heliocentric distance of 2.71~au. A search of shallower data obtained from the Zwicky Transient Facility (ZTF) survey, conducted when components A and B receded from 2.5 to 2.9~au from the Sun, yielded no detections, suggesting that the fragments had either ceased to exist or were not substantially active during this period. The uncertain fate of C/2019 Y4 highlights the challenge in understanding the state of presumably disrupted comets. Our review of six other presumably disrupted long-period comets reveals that negative confirmations were not reported for half of them, leaving their disruption presumptions inconclusive.

The nature of the supernova remnants (SNRs) 3C 397 and W49B has long been a subject of debate, with prior studies offering conflicting interpretations between thermonuclear and core-collapse scenarios. To help settle this debate, we present a systematic, spatially resolved, spectroscopic analysis of both remnants using XMM-Newton. By applying multi-component thermal models, we derive key physical properties including elemental abundances, ejecta temperatures, ambient densities, and explosion energetics. We compare the inferred metal abundance ratios to a wide range of core-collapse and thermonuclear nucleosynthesis models, including new models whose explosion energies differ from the canonical value of $10^{51}$ ergs. We find that the observed Fe/Si and Ca/Si ratios in both SNRs are best matched by certain thermonuclear models. However, no model fully reproduces the complete set of observed abundance patterns. In 3C 397, high Fe enrichment and spatial abundance variations suggest interaction with a dense progenitor environment, and W49B's composition is overall consistent with a thermonuclear origin; however both require a low energy ($\sim 10^{50}$ erg) supernova explosion. We additionally map the Fe K$\alpha$ line centroid energies and find a spread, with W49B falling within the core-collapse region -- highlighting both environmental complexity and the limitations of this diagnostic for supernova classification. Our results highlight the need for caution in relying on any single diagnostic or nucleosynthesis model for supernova typing, underscore the need for improved nucleosynthesis models, and motivate future high-resolution, high-throughput observations.

We present the application of the data-driven branch of the MURaM code to the extensively studied flare-productive active region 11158. We refine the hybrid model strategy, which was described in the earlier paper of this series, to model the emergence of the active region during 4 solar days starting shortly before 2011 February 11 and the eruption of an X2.2 flare on February 15. After 4 days of evolution, a major eruption of a magnetic flux rope occurs in the simulation at approximately 3 hours (3\% difference) before the real flare. The eruption leads to magnetic reconnection that contributes to bulk heating in the chromosphere and corona. The deposition of flare energy in the chromosphere causes strong condensations and evaporations, which fill hot post-flare loops and bright flare ribbons that exhibit separation and extension similar to the observed ribbon evolution. The synthesized soft X-ray flux corresponds to X class, which is close to the real event. The upward eruption of the flux rope leads to a piston-driven shock and horizontal expansion that exerts a strong downward impact on the lower atmosphere and generate an apparently fast-propagating chromospheric Moreton wave. We conclude that the data-driven radiative simulation of this active region can reproduce the key observational results of the real flare and demonstrate the great potential of this method for studying solar eruptions in a realistic corona environment.

Interchange reconnection is believed to play a significant role in the production of solar jets and solar wind. However, the dynamics of interchange reconnection in the low corona might be more complex than recognized before in higher temporal and spatial resolutions. Using unprecedentedly high-resolution observations from the Extreme Ultraviolet Imager (EUI) onboard the Solar Orbiter, we analyze the dynamics of interchange reconnection in a small-scale fan-spine-like topology. Interchange reconnection that continuously occurs around the multi-null points of the fan-spine-like system exhibits a quasi-periodicity of ~200 s, nearly covering the entire evolution of this system. Continuous evolution and reversal of multiple current sheets are observed over time near the null point. These results reveal that the dynamics of interchange reconnection are likely modulated by the emerging magnetic structures, such as mini-filaments and emerging arcades. Moreover, a curtain-like feature with a width of 1.7 Mm is also observed near the interchange reconnection region and persistently generates outflows, which is similar to the separatrix curtain reported in the pseudo-streamer structure. This study not only demonstrates the complex and variable reconnection dynamics of interchange reconnection within small-scale fan-spine topology but also provides insights into the self-similarity of magnetic field configurations across multiple temporal and spatial scales.

We present new radio-continuum observations of the Large Magellanic Cloud (LMC) supernova remnant (SNR) N206, which we give the name "Goat's Eye". Goat's Eye contains an interior radio structure that is likely a pulsar wind nebula (PWN), which we analyse in further detail. We use new radio observations from the Australia Telescope Compact Array (ATCA) telescope, as well as several archival radio observations, to calculate spectral indices, and find a steep spectral index for the whole SNR ($\alpha = -0.60\pm0.02$), and a flatter spectral index for the PWN ($\alpha = -0.16\pm0.03$). We also measure the polarisation and magnetic field properties of the PWN. Previously reported as a linear structure, the new observations show an unusual "zig-zag"-like structure, visible in radio-continuum total intensity, linear polarisation, and magnetic field orientations. The origin of this zig-zag structure is unclear, but we propose some origin scenarios that will require further observations to differentiate between.

We present Egent, an autonomous agent that combines classical multi-Voigt profile fitting with large language model (LLM) visual inspection and iterative refinement. The fitting engine is built from scratch with minimal dependencies, creating an ecosystem where the LLM can reason about fits through function calls-adjusting wavelength windows, adding blend components, modifying continuum treatment, and flagging problematic cases. Egent operates directly on raw flux spectra without requiring pre-normalized continua. We validate against manual measurements from human experts using 18,615 lines from the C3PO program across 84 Magellan/MIKE spectra at SNR~50-250. We find per-spectrum systematic offsets between Egent and expert measurements, likely arising from differences in global continuum placement prior to manual fitting; after accounting for these offsets, the agreement is 5-7 milliangstrom. The LLM's primary role is quality control: it confirms good fits (~60-65% of lines are LLM-refined and accepted), flags problematic cases (~10-20%), and occasionally rescues edge cases where tool use improves fits. Agreement between GPT-5 and GPT-5-mini confirms reproducibility, with GPT-5-mini enabling low-cost analysis at ~200 lines per US dollar. Every fit stores complete Voigt parameters, continuum coefficients, and LLM reasoning chains, enabling exact reconstruction without re-running. Egent compresses what traditionally requires months of expert effort into days of automated analysis, enabling survey-scale EW measurement. We provide open-source code at this https URL, including a web interface for drag-and-drop analysis and a local LLM backend for fully offline operation on consumer hardware.

To investigate whether apparently non-repeating Fast Radio Bursts (FRBs) are truly one-off transients, we conducted systematic follow-up observations of 27 out of 81 non-repeating FRBs identified in the Parkes Transient Database. Using 59.0 hours of data from the Parkes Ultra-Wideband Low (UWL) receiver and 6.3 hours from the Five-hundred-meter Aperture Spherical Telescope (FAST) 19-beam receiver, we searched for repeated bursts from these sources. No additional bursts were detected from any of the 27 FRBs. Combining these non-detections with prior archival observations, we derived stringent upper limits on their repetition rates above 1 Jy under two statistical models: Poisson process constraints range from $\sim10^{-3.5}$ to $10^{-1.9}\,\mathrm{h^{-1}}$, while Weibull process constraints range from $\sim10^{-3.4}$ to $10^{-1.5}\,\mathrm{h^{-1}}$. These limits are approximately an order of magnitude stricter than those reported in previous studies. By applying consistent observational setups and analytical methodologies across all sources, the derived rate limits converge to a narrow, well-defined range. This suggests that these FRBs form a relatively homogeneous population with extremely low intrinsic activity rates.

In the wideband paradigm of pulsar timing, the time of arrival of a pulsar pulse is measured simultaneously with the corresponding dispersion measure from a frequency-resolved integrated pulse profile. We present a new method for performing wideband measurements that rigorously accounts for measurement noise. We demonstrate this method using observations of PSR J2124$-$3358 made as part of the Indian Pulsar Timing Array experiment using the upgraded Giant Metre-wave Radio Telescope, and show that our method produces more realistic measurement uncertainty estimates compared to the existing wideband measurement method.

We use several smoothed particle hydrodynamics+N-body models as part of the GASTRO library to study the role of high-density star-forming clumpy regions and a single merger on the formation of the $\alpha$-rich and $\alpha$-poor populations in the disk galaxies. These experiments are tailored to mimic what is expected to be the Gaia-Sausage/Enceladus (GSE) accretion event, which occurred circa 10 Gyr ago in the Milky Way (MW). We find that either an early clumpy phase or a retrograde merger significantly reduces the star formation rate (SFR) of the disk, giving rise to a chemical bimodality qualitatively similar to the MW's. The decrease of the SFR as the cause of the chemical bimodality is consistent with previous idealized and cosmological simulations. On the other hand, a prograde radial merger does not significantly modify the SFR of the disk, resulting in no clear chemical bimodality. We further show that stars originating from the inner regions ($R_{form}<4$ kpc) do not create the disk's chemical bimodality, although they can enhance it. Finally, only the models with an early clumpy phase can produce a significant fraction of old, age $>11$ Gyr, $\alpha-$poor stars with disk-like orbits, similar to what has been recently observed in the MW. Our results strengthen the case of clumpy disky galaxies observed at redshift $z\approx 1-2$ as likely progenitors of our Galaxy.

Using K giants from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) Milky Way (MW) Survey, we measure the shape, orientation, radial profile, and density anisotropies of the MW stellar halo over 8 kpc$<r_\mathrm{GC}<200$ kpc. We identify a triaxial stellar halo (axes ratio $10:8:7$), 43 degrees tilted from the disk, showing two break radii at $\sim16$ kpc and $\sim76$ kpc, likely associated with Gaia-Sausage/Enceladus (GSE) and Large Magellanic Cloud (LMC), respectively. The inner stellar halo ($<30$ kpc) is oblate and aligned with the disk, whereas the outer stellar halo becomes prolate and perpendicular to the disk, consistent with the Vast Polar Structure of MW satellites. The twisted halo may arise from the disk-halo angular momentum shift triggered by the infall of a massive satellite. The anisotropic density distribution of the stellar halo is also measured, with successful re-identification of the Hercules-Aquila Cloud South/North (HAC-N/-S) and Virgo overdensities (VOD). Break radii are found at 15/30 kpc for VOD/HAC-N(-S). We identify the LMC transient density wake with a break radius at 60 kpc in the Pisces overdensity region. We also find new observational evidence of the LMC collective density wake, by showing a break radius at $\sim$100 kpc in the northern Galactic cap with a clear density peak at 90 kpc. In the end, we found that more metal-poor halo stars are more radially extended. Our results provide important clues to the assembly and evolution of the MW stellar halo under the standard cosmic structure formation framework.

We perform a joint analysis of the Milky Way (MW) and Andromeda (M31) satellite populations to constrain the properties of fuzzy dark matter (FDM) and thermal-relic warm dark matter (WDM). We combine MW satellite observations from the Dark Energy Survey (DES) and Pan-STARRS1 (PS1) with M31 satellite data from the Pan-Andromeda Archaeological Survey (PAndAS), and model the corresponding observable satellite populations using the galaxy--halo connection model together with the appropriate selection functions. Uncertainties in the virial masses of the MW and M31 are incorporated through host-mass priors that linearly scale the relevant model parameters, allowing us to infer the full posterior distributions of all parameters. For the FDM case, we obtain $m_{\mathrm{FDM}} > 1.75 \times 10^{-20}~\mathrm{eV}$ (95% CL) and $m_{\mathrm{FDM}} > 1.41 \times 10^{-20}~\mathrm{eV}$ (20:1 posterior ratio). For thermal-relic WDM, we find $m_{\mathrm{WDM}} > 6.22~\mathrm{keV}$ (95% CL) and $m_{\mathrm{WDM}} > 5.75~\mathrm{keV}$ (20:1 posterior ratio). These results represent a moderate improvement over MW-only constraints, and provide the strongest constraints to date on the FDM and WDM derived from satellite galaxy populations in the Local Group.

Double white dwarf binaries (DWDBs) with Helium components are progenitors to AM CVn binary systems. Their Galactic production rate may be given by the number of stably mass transferring DWDBs in the Milky Way. The theoretical criteria for stable mass transfer in DWDBs is calculated assuming that component white dwarfs are completely cold and degenerate. Respective fractions of surviving AM CVn and DWDB which merge are then calculated by applying this criteria to population synthesis estimates for Galactic DWDB. However, emerging observations of the local DWDB population suggest that Helium white dwarf (He WD) components are typically hot, and only partially degenerate when they begin mass transferring. Using recent numerical simulations of He WD donors in DWDBs, we qualitatively describe a temperature dependent stable mass transfer criteria for Galactic DWDBs. Mass transfer is even more stable than previously thought, or equivalently, DWDB mergers are even rarer. Realistic finite temperature treatments will deepen the dearth in observed AM CVn binaries compared to DWDB merger products.

The rate of observed gravitational waves (GWs) from neutron star-neutron star (NS-NS) mergers detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO) indicates the existence of more than one short gamma-ray bursts (sGRBs) similar to GW/sGRB 170817A within the total gamma-ray bursts (GRBs) recorded by satellite detectors such as BATSE, Fermi-Gamma-ray Burst Monitor (Fermi-GBM), and Swift-Burst Alert Telescope (Swift-BAT). We investigated sGRBs in the Fermi-GBM dataset based on their MeV-GeV $\gamma$-ray emission features, to identify sGRBs similar to sGRB 170817A. Any addition of such events can impact the rate of NS-NS CBC events observed by LIGO. SGRB 170817A exhibits two distinct emission components: a non-thermal peak and a thermal component. We adopted a multifaceted approach to identify analogous sGRBs, which involved computing the hardness ratios $HR_{1}$ and $HR_{2}$ and then clustering them via the K-means algorithm. Our further studies reveal the presence of eight such events in Fermi-GBM data, which will enable us to calculate the rate of electromagnetic (EM) counterparts associated with LIGO GW events (GW+sGRB events) across all observing runs. Giving an estimation, by the end of the $O_4$ LIGO run, there could be nearly 5 GW+sGRB events. Deviation from this number may raise concerns about our understanding of the evolution of such events over distance.

The Sun produces a steady signal of high-energy gamma rays through interactions of Galactic cosmic rays (GCRs) with its atmosphere. Observations with Fermi-LAT and HAWC have revealed a gamma-ray flux significantly higher than early theoretical predictions, with unexpected temporal and spectral features that suggest a crucial role of the solar magnetic field. In this work, we model GCR-induced gamma-ray emission at the solar disk using the CRPropa framework with realistic hadronic interactions, chromospheric density profiles, and several magnetic field configurations over the solar cycle. This allows us to quantify the gamma-ray emission of the entire solar disk for different phases of the solar activity cycle and we present, for the first time, maps of the production locations of gamma rays on the solar surface. We consider both mono-energetic and realistic power-law injection spectra in a simplified dipole-quadrupole current sheet model and potential-field source surface (PFSS) extrapolations for Carrington rotations during solar maximum and minimum. Our results show that magnetic mirroring and large-scale field topology strongly affect the spectral shape and spatial distribution of the emission, with slightly enhanced fluxes predicted at solar minimum. While our simulated baseline fluxes remain below observations, additional effects, such as heavier nuclei, Parker-field mirroring, and deeper atmospheric interactions, could result in further enhancements of fluxes closer to observational values. Hadronic interactions do not only produce gamma rays but also neutrinos. We estimate the expected neutrino flux from the Sun based on our gamma-ray predictions. We find that the expected flux is slightly below current upper limits from IceCube.

We report the discovery of 45 compact hub-filament systems (HFSs; median size $\sim$2.4 pc) in infrared-dark clouds (IRDCs) in the W33 complex, located at the junction of the Scutum and Norma spiral arms. Using {\it Spitzer} 8 and 24 micron, and unWISE 12 $\mu$m images, HFSs are identified as regions where three or more filaments converge onto a central hub, appearing as absorption features toward IRDCs. In each IRDC, HFSs mainly lie at the intersections of elongated substructures, associated with groups of protostars and lacking radio continuum emission. Minimum Spanning Tree (MST) analysis shows that protostars are closely associated with the HFSs, with protostellar core separations of $\leq$ 0.7 pc, indicating strong clustering within fragmented structures. The HFSs form two main groupings spanning 10--15 pc, with member separations of 1--3.3 pc. Around 65\% are tightly clustered ($<$ 2 pc), exhibiting rich small-scale structures and emphasizing the uniqueness of the complex. MST analysis of ALMAGAL 1.38 mm continuum cores -- predominantly low-mass and embedded in ten HFSs -- reveals a median core separation of $\sim$0.03 pc. The protostellar spacing ($\sim$0.7 pc) significantly exceeds the thermal Jeans length ($\sim$0.08 pc for temperature $\sim$18 K and density $\sim$10$^{5}$ cm$^{-3}$), whereas the core spacing is smaller than the Jeans length, suggesting that thermal fragmentation may influence core formation but alone cannot explain the larger-scale protostellar distribution. All these findings together support a picture in which fragments of clouds/filaments form clumps hosting compact HFSs that facilitate efficient and clustered star formation, often yielding massive stars.

The interactions between stars and their orbiting planets, driven by forces such as stellar radiation and gravity, play an essential role in shaping exoplanetary atmospheres and gas-rich debris discs. One way to look into the composition of these environments is to observe how they can contaminate the stellar photospheres. For that, we examine how stellar radiation pressure and gravity influence atomic species and analyse their effects across various stellar effective temperatures. Using the radiative-to-gravitational force ratio, we determined the atomic movement direction and assessed the velocity boost imparted to neutral atoms escaping from exoplanet atmospheres or debris discs. Incorporating the solar far ultraviolet/extreme ultraviolet spectrum to address flux discrepancies in the {\sc{atlas9}} model, we find that radiation affects atoms differently according to their ionisation state, with highly ionised species less affected by stellar radiation. Our results conclude that the stars most suitable for observing stellar contamination are those between 6,500 and 8,000 K, with neutral noble gases and ionised iron-peak elements as the most likely contaminants.

In space plasmas, the rarity of collisions leads to complex structures in the velocity space where a turbulent cascade of the velocity distribution function fluctuations is thought to occur. Previous studies have explored this phenomenon using the Hermite decomposition of the ion velocity distribution function (VDF) in both magnetosheath data and numerical simulations. In this work, we investigate the Hermite spectrum of the ion VDFs measured by Parker Solar Probe in the inner heliosphere. We analyze a superalfvénic stream at a radial distances of $R \approx 28 R_{sun}$ and a subalfvénic at $R \approx 11 R_{sun}$, the former characterized by a prevalence of VDFs with suprathermal beams (also known as hammerhead). The Hermite analysis is also compared with various proxies of energization and dissipation, in order to establish a connection between turbulent cascades in real space and those in the velocity space. A qualitative agreement between the energization proxies and the Hermite analysis is observed. The results are suggestive of the presence of a dual cascade in real and velocity space.

Covariant density functionals have been successfully applied to the description of finite nuclei and dense nuclear matter. These functionals are often constructed by introducing density dependence into the nucleon-meson couplings, typically through functions that depend only on the vector, i.e., proper baryon density. In this work, we employ a Bayesian framework to investigate how different parametrizations, characterized by distinct functional forms and by their dependencies on vector and scalar densities, affect the properties of dense matter and compact stars. Our analysis demonstrates that although all considered parametrizations yield broadly comparable inferences, the differences in the equation of state and the symmetry energy remain significant at suprasaturation densities, reflecting the sensitivity to the chosen functional form of the density dependence. We find that allowing the nuclear saturation properties in the isoscalar channel, including the skewness coefficient $Q_{sat}$, to be freely adjusted provides adequate flexibility for the current modeling of nuclear and neutron star matter. In contrast, the isovector channel requires further refinement, with freedom extended at least up to the curvature coefficient $K_{sym}$ to capture variations in the symmetry energy and particle composition at high densities. This work advances prior studies by implementing a rational-function parametrization of the density dependence, informed and constrained by multimessenger astrophysical observations.

In high-contrast imaging, a novel detection algorithm for angular differential imaging (ADI) sequences has recently been introduced: the Regime Switching Model (RSM). In this study, we apply the RSM algorithm to analyze the F150 sample from the SHINE high-contrast imaging survey carried out with VLT/SPHERE, aiming to enhance detection limits and identify new exoplanet candidates. Additionally, we investigate how environmental conditions influence post-processed noise distributions and detection thresholds. We generate detection maps and contrast curves for 213 observations in the F150 SHINE sample using the RSM algorithm. A clustering approach based on environmental parameters is used to group observations with similar noise characteristics. We propose two methods for defining radial detection thresholds in the RSM maps: fitting a log-normal distribution to the post-processed noise and maximizing the F1 score. We also assess the performance of various combinations of post-processing techniques within the RSM framework to identify optimal configurations. This study demonstrates the utility of clustering based on observational parameters, effectively distinguishing features like wind-driven halos and low-wind effects. Detection thresholds vary significantly across clusters, differing by up to a factor of 10, highlighting the importance of considering observational environments. Log-normal thresholds provide conservative, noise-aware limits, while F1 score-based thresholds offer observation-specific results, both showing compatibility overall. RSM improves detection limits by an average factor of two at 1arcsec and five at inner working angles compared to standard PCA processing. This study reports more than 30 newly detected signals, including one promising candidate awaiting second-epoch confirmation.

The Hubble Missing Globular Cluster Survey has secured high quality astro photometric data in two bands for 34 clusters never observed with HST. When combined with Gaia positional measurements, this data set enables the investigation of the bulk motion and the internal kinematics of these poorly studied clusters to an unprecedented level of detail. Focusing on the case of NGC 6749, we here showcase how the combined Gaia HST proper motions have a quality sufficient to accurately assess the cluster stellar membership, determine its absolute proper motion with a precision superior to Gaia, and to investigate its kinematic profile for the first time. Proper motions are determined using the public code GAIAHUB, which for NGC 6749 combines data sets separated in time by 8 years. The resulting measurements improve the precision of Gaia proper motions by a factor of 10 at the faint end, and enable recovering the proper motion for 662 stars for which Gaia could only measure the positions. These proper motions are efficient in decontaminating the colour magnitude diagram of NGC 6749, and make it possible to compare the efficacy of a method of statistical decontamination that relies only on the photometric information extracted from the HST parallel fields. Finally, using the sample of best measured proper motions we determine the velocity dispersion and anisotropy profiles of NGC 6749, that reveal an isotropic behaviour in the cluster inner regions and a slight radial anisotropy outside 1.5 half light radii. The proper motions and the code to statistically decontaminate the clusters color magnitude diagram are made available as public products of the survey.

Aims: We aim to characterize the physical and morphological properties of the binary system AFGL 4106, composed of two evolved massive stars. Understanding its mass-loss processes and circumstellar environment offers insight into the late stages of stellar evolution in massive binary systems. Methods: We obtained high-angular--resolution, high-contrast imaging using VLT/SPHERE with ZIMPOL (optical) and IRDIS (near-infrared) across multiple filters. We used aperture photometry to extract the spectral energy distributions (SEDs) of each star, and applied radiative transfer modelling to study the system and its surrounding dusty environment. Results: The observations resolve both components of the binary and unveil a complex, dusty nebula featuring asymmetric structures and cavities. SED fitting yields stellar temperatures of T$_1 = 6723\pm196$ K and T$_2 = 3394\pm264$ K, along with bolometric luminosities of L$_1 = (7.9 \pm 0.18) \times 10^4$ L$_\odot$ and L$_2 = (3.8 \pm 0.11) \times 10^4$ L$_\odot$. These values support the classification of the primary as being in a post-red supergiant (post-RSG) phase and the secondary as an active red supergiant (RSG). The luminosity ratio, combined with the inferred radii, indicates that both stars are at close yet distinct stages of their evolution. The binary is surrounded by an extended shell whose asymmetric morphology and large-scale features suggest interaction with the stellar winds and interstellar medium (ISM), and possibly the presence of a third, undetected companion. Conclusions: These observations provide the first resolved view of AFGL 4106's system and its dusty envelope. Our analysis sets constraints on the physical properties and evolutionary status of the system. This work contributes to understanding mass-loss processes in massive binaries and the shaping of nebulae around evolved stars.

Ultra-faint dwarf galaxies (UFDs) are among the oldest and most metal-poor stellar systems in the Universe. Their metallicity distribution encodes the fossil record of the earliest star formation, feedback, and chemical enrichment, providing crucial tests of models of the first stars, galaxy assembly, and dark matter halos. However, due to their faint luminosities and the limited number of bright giants, spectroscopic studies of UFDs typically probe only small stellar samples. Here, we present an analysis of multi-epoch Hubble Space Telescope and James Webb Space Telescope observations of the UFD Bootes I. Using deep color-magnitude diagram in the F606W and F322W2 bands, extending from the subgiant branch to the M-dwarfs, and stellar proper motions to identify likely members, we obtained an unprecedentedly clean census of the system. The exquisite quality of the diagram, combined with the sensitivity of M-dwarf colors to metallicity, allowed us to constrain the metallicity distribution in a large stellar sample. As a first step, we derived the binary fraction in Bootes I. This is crucial, since binaries can bias kinematic mass estimates, affect stellar population analyses, and shape the photometric signatures used to infer metallicity. We find that 20$\pm$2% of stellar systems in Bootes I are binaries with mass ratios larger than 0.4, corresponding to a total binary fraction of $\sim$30%. This value is comparable to the binary fractions observed in globular clusters of similar stellar mass, suggesting that the presence of dark matter does not significantly affect the binary properties of Bootes I. We then exploited the metallicity sensitivity of M-dwarf colors to derive the metallicity distribution function. We find that most of the stars $\sim$85% have [Fe/H]<-2, and that roughly $\sim$17% have [Fe/H]<-3.

The well-observed type IIP SN 2024bch with the short plateau is shown to be an outcome of the red supergiant explosion with the presupernova mass of 14-15 Msun, the explosion energy of 2x10^{51} erg, and presupernova radius of 1250 Rsun. The early gamma-ray escape demonstrated by the radioactive tail is due to the large Ni-56 extension up to 7400 km/s. The early-time spectral evolution indicates the presence of the circumstellar dense confined envelope with the mass of 0.003-0.006 Msun within 6x10^{14} cm. The deceleration of the outermost ejecta implies the wind with the mass-loss rate of about 6x10^{-4} Msun/yr. The inferred mass-loss rate is by one-two order larger compared to most of type IIP supernovae, but comparable to the wind of type IIL SN 1998S. The asymmetry of the broad H-alpha component on day 144 powered by the circumstellar interaction is the outcome of the Thomson scattering and absorption in the Paschen continuum in the unshocked ejecta.

Accurate forecasts of Coronal Mass Ejection (CME) arrival times and impact geometry remain a major challenge for space-weather operations. Coronagraph-based techniques typically achieve mean absolute errors of order ten hours, while in situ measurements at L1 provide excellent magnetic-field information but only tens of minutes of warning. In this work we introduce X-CME, a framework that links in situ flux-rope reconstructions at intermediate heliocentric distances with a physics-based CME propagation model. The internal magnetic structure is obtained with an elliptical cylindrical, radial poloidal flux-rope model and embedded into a tapered torus CME geometry. The subsequent propagation is computed by solving a drag-based equation of motion in a Parker solar-wind background, including gravitational deceleration, self-similar expansion of the cross section, and an explicit calculation of the CME wetted area and swept area in the ecliptic plane. We apply X-CME to two events observed by the Parker Solar Probe and Solar Orbiter spacecraft, respectively, and propagate the reconstructed structures to Earth and Mars. For both cases, the model reproduces the observed in situ signatures at L1 and predicts the CME arrival time at Earth with errors of a few hours (typically about 2-4 hours), while correctly distinguishing between central encounters and glancing blows. These results demonstrate that combining intermediate-distance magnetic reconstructions with a geometrically consistent propagation scheme can substantially improve CME arrival-time forecasts and impact assessment in the inner heliosphere.

Modeling how supermassive black holes co-evolve with their host galaxies is notoriously hard because the relevant physics spans nine orders of magnitude in scale-from milliparsecs to megaparsecs--making end-to-end first-principles simulation infeasible. To characterize the feedback from the small scales, existing methods employ a static subgrid scheme or one based on theoretical guesses, which usually struggle to capture the time variability and derive physically faithful results. Neural operators are a class of machine learning models that achieve significant speed-up in simulating complex dynamics. We introduce a neural-operator-based ''subgrid black hole'' that learns the small-scale local dynamics and embeds it within the direct multi-level simulations. Trained on small-domain (general relativistic) magnetohydrodynamic data, the model predicts the unresolved dynamics needed to supply boundary conditions and fluxes at coarser levels across timesteps, enabling stable long-horizon rollouts without hand-crafted closures. Thanks to the great speedup in fine-scale evolution, our approach for the first time captures intrinsic variability in accretion-driven feedback, allowing dynamic coupling between the central black hole and galaxy-scale gas. This work reframes subgrid modeling in computational astrophysics with scale separation and provides a scalable path toward data-driven closures for a broad class of systems with central accretors.

We present 22 sets of light curves and one radial velocity curve for the W UMa-type total eclipse contact binary system V2790 Ori, derived by combining all available public photometric data, the photometric data in previous studies, and our own spectroscopic and decade-long photometric observations. Our simultaneous analysis of the light curves and radial velocity curve shows that V2790 Ori is a W-subtype contact binary with a mass ratio of $q = 0.322(\pm0.001)$ and a shallow contact degree of $14.8(\pm0.6)\%$. The orbital period analysis based on 445 eclipsing minima reveals a secular decrease at a rate of $\dot P = -3.18 (\pm 0.75) \times 10^{-8}\mathrm{d~yr^{-1}}$, superimposed with a cyclic variation with an amplitude of $A = 8.98 (\pm 2.19) \times 10^{-4}~\mathrm{d}$ and a period of $ P_3 = 7.44 (\pm 0.52)~\mathrm{yr}$. The secular decrease is caused by AML via magnetic braking, while the cyclic period variation is explained by the light-travel time effect due to a third body, which is likely to be a brown dwarf. Furthermore, our analysis indicates a mass transfer from the more massive component to the less massive one at a rate of $1.22(\pm0.29) \times 10^{-8}~\mathrm{M_{\odot}~yr^{-1}}$. A model with a cool spot on each component was adopted to fit the O'Connell effect observed in the light curves. Since the O'Connell effect varies over time, we identified a solar-like magnetic activity cycle with a period of approximately 5.4 yr by analyzing the magnitude difference ($\Delta m$) at the two light maxima and the O'Connell effect ratio. In addition, evolutionary analysis suggests that V2790 Ori is a newly formed contact binary that evolved from a detached phase into the present contact configuration.

Young, massive stellar clusters offer a prime setting to explore brown dwarf (BD) formation under high densities and intense UV radiation. Trumpler 14 (Tr 14), a 1 Myr-old cluster located at a distance of 2.4 kpc, and particularly rich in O-type stars, is an ideal target for such a study. Our goal is to measure the initial mass function (IMF) in the young massive, high UV flux cluster. We present the deepest study to date of the IMF in Tr 14 based on GeMS/GSAOI imaging. We construct the IMF using both the Besancon Galactic model and an observational control field from VISTA for background correction. Completeness was assessed using artificial-star tests and applied to the IMF derivation. We estimate the IMF down to the 20% completeness limit found at 0.01 MSun. Using the control field-based IMF as our primary result, we find a slope of alpha=0.14+-0.19 for masses between 0.01-0.2 MSun, and alpha=1.72+-0.04 for 0.2-4.5 MSun. The low-mass slope is largely influenced by the incompleteness-affected lowest bin; excluding it brings our results into agreement with those from other young clusters. The resulting median for the star-to-BD ratio in the 0.03-1 MSun mass range is 4.0, with a 95% confidence interval of 2.8-5.8. Our analysis reveals that Tr 14 hosts a relatively flat substellar IMF, but this is strongly influenced by the lowest-mass bin, which may be affected by incompleteness. When that bin is excluded, the IMF becomes consistent with those of other regions. The star-to-BD ratio falls within the usually observed 3-6 range, indicating that brown dwarfs with masses above 0.03 MSun form with similar efficiency across environments. However, the relative lack of objects below this threshold suggests that the presence of an environment with both high stellar density and FUV flux may play a role in shaping the IMF by suppressing the formation of BDs at masses < 0.03 MSun.

We develop an accurate and computationally efficient emulator to model the gravitational lensing magnification probability distribution function (PDF), enabling robust cosmological inference of point sources such as supernovae and gravitational-wave observations. We construct a pipeline utilizing cosmological $N$-body simulations, creating past light cones to compute convergence and shear maps. Principal Component Analysis (PCA) is employed for dimensionality reduction, followed by an eXtreme Gradient Boosting (XGBoost) machine learning model to interpolate magnification PDFs across a broad cosmological parameter space ($\Omega_m$, $\sigma_8$, $w$, $h$) and redshift range ($0.2 \le z \le 6$). We identify the optimal number of PCA components to balance accuracy and stability. Our emulator, publicly released as ace_lensing, accurately reproduces lensing PDFs with a median Kullback-Leibler divergence of $0.007$. Validation on the test set confirmed that the model reliably reproduces the detailed shapes and statistical properties of the PDFs across the explored parameter range, showing no significant degradation for specific parameter combinations or redshifts. Future work will focus on incorporating baryonic physics through hydrodynamical simulations and expanding the training set to further enhance model accuracy and generalizability.

An empirical baryon density profile can be included in the Schrödinger-Poisson (SP) equations to influence the fuzzy dark matter (FDM) soliton formation. However, to probe the effects of baryon on the other dynamical evolutions of the FDM soliton, its equation of motion (EoM) inside the corresponding FDM soliton is needed. In this paper, given an empirical baryon density profile, we first provide the cylinderical symmetric FDM soliton solution about the FDM density and the total potential of FDM and baryon. Then, we build an analytical baryon EoM from the obtained FDM density and total potential by machine learning. Finally, we check that this baryon EoM works as well as an empirical baryon density profile for the FDM soliton formation, with the fractional errors $\lesssim0.04$. It should also work well for some other simple FDM soliton evolutions.

Ultra-high-energy photons play an important role in probing astrophysical models and beyond-Standard-Model scenarios. We report updated limits on the diffuse photon flux using Telescope Array's Surface Detector data collected over 14 years of operation. Our method employs a neural network classifier to effectively distinguish between proton-induced and photon-induced events. The input data include both reconstructed composition-sensitive parameters and raw time-resolved signals registered by the Surface Detector stations. To mitigate biases from Monte Carlo simulations, we fine-tune the network with a subset of experimental data. The number of observed photon candidates is found to be consistent with the expected hadronic background, yielding upper limits on photon flux $\Phi_\gamma(E_\gamma > 10^{19} \text{eV}) < 2.3 \cdot 10^{-3} $, and $\Phi_\gamma(E_\gamma > 10^{20} \text{eV}) < 3.0 \cdot 10^{-4} $ $ (\text{km}^2 \cdot \text{sr} \cdot \text{yr})^{-1} $.

We present a rigorous proof establishing the mathematical equivalence between two independent criteria for the marginal stability of multi-fluid relativistic stars: the dynamical criterion based on the vanishing of the fundamental radial pulsation mode's eigenfrequency, and the static criterion derived from the geometric alignment of mass and particle number gradients in the parameter space. Leveraging this equivalence, we introduce a powerful and computationally efficient framework as an upgraded version of the critical curve method, to systematically map the stability boundaries for multi-fluid mixed stars across the entire parameter space of central pressures. Our analysis, applied to a variety of nuclear and dark matter equations of state, reveals the existence of stable region in the observable mass-radius diagram. By resolving degeneracies with 3-dimensional Mass-Radius-Pressure diagrams, we provide a complete topological view of the ensemble. This work supplies a robust theoretical foundation for interpreting multi-messenger astronomical observations and constraining the properties of dark matter.

Concerning the stability of two-fluid star models, we prove the rigorous equivalence of two independent determining methods for mixed stars, after a brief review of the hybrid star case. Our derivations apply to general multi-fluid cases, and here we take dark matter admixed neutron star models for example, demonstrating a stability boundary distinct from the single-fluid case. Stable configurations form a surface in the three-dimensional parameter space of (either) central pressure, mass, and radius, which yields a group containing stable mixed stars. This group includes twin stars with identical masses and radii but different interior structures. These results can help interpret compact star observations and constrain dark matter properties through astrophysics.

Active region NOAA 14274 produced some of the strongest flares of Solar Cycle 25, including the X1.2 and X5.1 flares on 10 and 11 November 2025, respectively. We present the first large mosaic of speckle-restored images obtained with the improved High-resolution Fast Imager (HiFI+) at the 1.5-meter GREGOR solar telescope at the Observatorio del Teide in Izaña, Tenerife, Spain. The observations were obtained approximately 30 minutes before the onset of the X1.2 flare. The active region exhibited strongly curved penumbral filaments, sunspot rotation, and shear motions along the polarity inversion line (PIL), which led to a highly stressed magnetic field configuration that stored sufficient energy to release multiple M- and X-class flares. The first flare signatures appeared as small-scale brightenings, each with a width of a few tenths of an arcsecond, that trace penumbral filaments in the trailing sunspot.

By combining high-resolution observations from JWST and HST, we have measured the stellar masses, star formation rates (SFRs), and multi-wavelength morphologies of galaxies in the CANDELS fields. Furthermore, based on rest-frame 1 $\mu$m morphologies, we have derived spatially resolved stellar mass and SFR surface density ($\Sigma_*$ and $\Sigma_{\rm SFR}$) profiles for 46,313 galaxies with reliable structural measurements at $0<z<4$ and $\log(M_\ast /M_{\odot})>8$, and provide the corresponding catalogue. For star-forming galaxies (SFGs), our results show excellent consistency with previous studies in terms of the star formation main sequence and the size-mass relation, demonstrating the robustness of our stellar mass and SFR measurements. For spatially resolved profiles, we find that at higher redshifts ($z>2.5$), the median radial profile of $\Sigma_{\rm SFR}$ is nearly parallel to but slightly steeper than that of $\Sigma_*$. This results in mildly negative gradients in the specific SFR (sSFR) profiles across all stellar mass bins considered. These findings indicate that galaxies at $z>2.5$ cannot grow in size via only in-situ star formation, challenging the understanding of galaxy size evolution beyond the cosmic noon. In contrast, at $z<2.0$, the sSFR profiles transition to exhibit more and more positive gradients at lower redshifts, consistent with an inside-out growth scenario where star formation preferentially expands the galactic outskirts.

Much of a solar flare's energy is thought to be released in the continuum. The optical continuum (white light) is of special interest due to the ability to observe it from the ground. We aim to investigate the prevalence of white-light (WL) emissions in simulations of purely electron beam-driven solar flares, what determines the occurrence of these enhancements, and the underlying causes. We utilized the F-CHROMA grid of flare simulations created using the radiative hydrodynamics code RADYN. We probed the spectral index, total energy, and low-energy cutoff to draw conclusions about their relationships to the white-light intensity. Furthermore, we calculated the 6684 Å continuum intensities, the Balmer, and the Paschen ratios. Finally, we analyzed two particular cases, one with high 6684 Å intensity and one with a large Balmer ratio, to determine the dominant mechanisms in these simulations. 33 of the 84 flares included in the F-CHROMA grid show white-light intensity enhancements that exceed 0.1% relative to the pre-flare level. We conclude that, with the parameters presented in the F-CHROMA grid, purely electron beam-driven simulations of solar flares are not able to reproduce observed WL enhancements, as the maximum enhancements in the grid are below 4%. The total energy (which is correlated with the maximum beam flux) is the main factor for deciding whether excess white-light emissions will be detectable. There is a linear relationship between the Balmer (and Paschen) ratio and the relative continuum increase. Both case studies show that during the time of maximum WL excess, hydrogen ionization and subsequent recombination in an optically thin medium is the dominant mechanism for WL continuum emission enhancements. Increased H$^-$ emission in the photosphere as a result of radiative backwarming becomes dominant during the declining phase of WL emissions in both case studies.

With the giant exoplanet occurrence rate peaking around stars of 1.5-2 solar masses, there is strong motivation to characterize the disks that set their formation conditions. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) allow us to investigate both the availability of different molecules in disks and infer the radial distribution of elemental abundances, enabling us to make connections to exoplanet systems. Here we present a survey of six transition disks around young F-, A-, and B-type stars using ALMA. We find 13C18O, CS, SO, and H2CO in all six systems, as well as ten additional molecules in a subset of disks, including detections of H2S, 33SO, and CH3OCH3. Using these data, and literature data where available, we construct the first comprehensive picture of Herbig disk chemistry. We find clear correlations between molecular tracers of C/O>1 environments (e.g., CS, C2H) and disk mass, as traced by C18O line flux. In contrast, tracers of C/O<1 environments (e.g., SO, CH3OH) do not show significant correlations with disk mass. Interestingly, these molecules are relatively brighter in lower-mass disks, with their presence primarily linked to disks with central cavities and spirals. Finally, we show that the observed chemical diversity seen across Herbig disks leads to varying C/O regimes at the orbital radii of candidate proto-planets identified within these disks. When comparing these inferred disk C/O ratios with those measured for directly imaged exoplanets, we find a notable overlap and show that low C/O gas is common on 10's of au scales in Herbig disks.

We present the continuation of a systematic search for new southern Galactic symbiotic stars, selecting candidates from the SuperCOSMOS H$\alpha$ Survey and 2MASS. Follow-up spectroscopy with the Southern African Large Telescope (SALT) was used to confirm their symbiotic nature and to characterize the cool and hot components of the full sample, including systems from earlier work. We report 14 newly confirmed bona fide symbiotic stars and identify 6 additional strong candidates. Photometric variability was examined using our data and archival light curves from multiple all-sky surveys. Most systems are variable, with the majority showing periodic modulation consistent with orbital motion or pulsations. Possible photometric orbital periods are reported for 19 confirmed and 3 candidate systems, pending spectroscopic confirmation. Eight objects exhibit signs of outburst activity. In one of the systems, multiple brightenings occur at similar orbital phases, closely resembling the evolution of FN Sgr, a symbiotic binary with a magnetic white dwarf. The peculiar variability of another symbiotic star is best explained by dust-obscuration events. These results expand the census of Galactic symbiotic stars.

On 9 October 2022, the brightest gamma-ray burst (GRB) ever recorded (GRB 221009A) was detected. Its remarkably bright emission, partially due to its close distance to Earth ($z\sim0.15$), makes this GRB a unique event. The outstanding characteristics of GRB 221009A, including the TeV detection by the LHAASO experiment, triggered deep follow-up observations of the source across all wavebands, including very-high-energy gamma rays with the first Large-Sized Telescope (LST-1) of the future Cherenkov Telescope Array Observatory. LST-1 observations started about one day after the onset of the prompt emission, under strong moonlight conditions. This resulted in a hint of a signal with a statistical significance of about 4$\sigma$. The monitoring of this source continued until the end of November 2022. This constitutes the deepest observation campaign performed on a GRB with LST-1. Here we show the results of this follow-up campaign.

We introduce a new method to infer the posterior distribution for planet occurrence rates from radial-velocity (RV) observations. The approach combines posterior samples from the analysis of individual RV datasets of several stars, using importance sampling to reweight them appropriately. This eliminates the need for injection-recovery tests to compute detection limits and avoids the explicit definition of a detection threshold. We validate the method on simulated RV datasets and show that it yields unbiased estimates of the occurrence rate in different regions, with increasing precision as more stars are included in the analysis.

Sub-Neptunes represent the largest exoplanet demographic, yet their bulk compositions remain poorly understood. Recent studies suggested that only very cold planets, such as Uranus and Neptune, could experience stratification of volatiles in their envelopes, implying that the envelopes of warmer sub-Neptunes instead have fully-miscible compositions. Here, we present ATHENAIA, an interior-atmosphere composition inference framework we leverage to assess the potential for water-hydrogen demixing on the $T_{\mathrm{eq}}=350$ K planet TOI-270 d, and more broadly for warm sub-Neptunes, using radiative-convective atmosphere models coupled to interior models. We find that the higher temperatures at which hydrogen and water demix in water-rich environments, combined with the shallower adiabatic gradients of water-rich envelopes, open a window for demixing on sub-Neptunes with bulk envelope metallicities of $\sim 100$ to $700\times$ solar, compatible with TOI-270 d. Demixing is easier to achieve on more massive and colder planets, but still broadly affects warm (330 to 500 K) metal-rich sub-Neptunes. Therefore, combining atmosphere metallicities with models of fully-miscible envelopes may lead to underestimated bulk envelope metallicities and mass fractions. Further, our modeling of TOI-270 d's envelope and interior reveals that, for a typical internal energy budget $T_\mathrm{int}$ of 25 K, the envelope-mantle boundary conditions likely preclude the presence of a molten magma ocean. This work encourages a reconsideration of the current paradigm for linking sub-Neptune atmospheres to their interiors and motivates further evolutionary modeling describing the onset of metallicity gradients in sub-Neptune envelopes.

We report the first observations in a rare family of class II methanol maser transitions in both CH$_3$OH and $^{13}$CH$_3$OH toward three southern high-mass star formation regions, along with the first maser detected in the $^{13}$CH$_3$OH line. The $8_2 \rightarrow 9_1 A^{-}$ methanol transition was observed in both CH$_3$OH and $^{13}$CH$_3$OH (at 28.9 GHz and 41.9 GHz respectively) toward three sources; G358.93-0.03, NGC6334I and G345.01+1.79, all of which are star formation regions with recent maser flaring events. We report the first maser detection of the 41.9 GHz line in $^{13}$CH$_3$OH toward G358.93-0.03 and the first confirmed maser detection of the 28.9 GHz line in CH$_3$OH toward NGC6334I. Additionally we report a maser detection of the 28.9 GHz line in CH$_3$OH toward G358.93-0.03, meaning that with our detection of the 41.9 GHz line, this is the first isotopic detection of these lines toward G358.93-0.03. The newly detected maser transitions are associated with the primary millimetre continuum sources (MM1) in both G358.93-0.03 and NGC6334I, within the varying positional uncertainties.

The fine structures of the solar chromosphere, driven by photospheric motions, play a crucial role in the dynamics of solar magnetic fields. Many have been already identified such as fibrils, filament feet, and arch filament systems. Still, high resolution observations show a wealth of structures that remain elusive. We have observed a puzzling, unprecedented chromospheric fibril singularity in close vicinity of a blow-out solar jet and a flaring loop. We aim to understand the magnetic nature of this singularity and the cause of its activity using coordinated high- resolution multi-wavelengths observations. We aligned datasets from Solar Orbiter, SST, IRIS, and SDO. We re-projected the Solar Orbiter datasets to match the perspective of the Earth-based instruments. We performed potential field extrapolations from Solar Orbiter/PHI data. We analysed the spatial and temporal evolution of the plasma structures and their link with the surface magnetic field. This leads us to derive a model and scenario for the observed structures which we explain in a general schematic representation. We have discovered a new feature, a singularity in the chromospheric fibril pattern. It is formed in a weak magnetic field corridor between two flux concentrations of equal sign, at the base of a vertically inverted-Y shape field line pattern. In this specific case some activity develops along the structure. Firstly a flaring loop at one end, secondly a blow-out jet at the other end, where a coronal null-point was present and associated with a chromospheric saddle point being located onto the fibril singularity. The observations sugge

We analyzed archival Herschel observations of water vapor emission toward the Horsehead photon dominated region (PDR), along with supporting ground-based and airborne observations of CO isotopologues and fine structure lines of ionized and atomic carbon to determine the distribution and abundance of water vapor in this low-UV illumination PDR. Water emission in the Horsehead nebula is very weak and, surprisingly, extends outward beyond other PDR tracers such as $^{12}$CO or [CI] 609 $\mu$m, reaching as far out as [CII] 158 $\mu$m. We model the observations using a newly developed PDR wrapper that takes into account the geometry of this region. PDR modeling of the molecular and atomic lines studied here provides strong constraints on the thermal pressure, but not on the UV illumination. Maximum model line intensities %typically agree to within ~40\% with the observations. and spatial profiles are well reproduced, except for CO isotopologues, where the increase on the illuminated side of the PDR is steeper than observed. Water vapor abundance in the model reaches $3.6 \times 10^{-7}$ at $A_V \sim 3$ mag. However, the ground state $o$-H$_2$O 557 GHz line is systematically overestimated by the models by at least a factor of 7 for any values of the model parameters. This line has a very high optical depth and the emergent line intensity is sensitive to radiative transfer effects such as line scattering by water molecules in a low-density halo surrounding the dense PDR and the assumed microturbulent line width. A more accurate model of the water surface chemistry is required.

We present the characterization of the TOI-6041 system, a bright ($V = 9.84 \pm 0.03$) G7-type star hosting at least two planets. The inner planet, TOI-6041b, is a warm Neptune with a radius of $4.55^{+0.18}_{-0.17}\,R_\oplus$, initially identified as a single-transit event in \textit{TESS} photometry. Subsequent observations with \textit{TESS} and \textit{CHEOPS} revealed additional transits, enabling the determination of its $26.04945^{+0.00033}_{-0.00034}$~d orbital period and the detection of significant transit timing variations (TTVs), exhibiting a peak-to-peak amplitude of about 1~hour. Radial velocity (RV) measurements obtained with the APF spectrographs allow us to place a $3\sigma$ upper mass limit of $28.9\,M_\oplus$ on TOI-6041b. In addition, the RV data reveal a second companion, TOI-6041c, on an 88~d orbit, with a minimum mass of $0.25\,M_{\mathrm{Jup}}$. A preliminary TTV analysis suggests that the observed variations could be caused by gravitational perturbations from planet c; however, reproducing the observed amplitudes requires a relatively high eccentricity of about 0.3 for planet c. Our dynamical stability analysis indicates that such a configuration is dynamically viable and places a $1\sigma$ upper limit on the mass of TOI-6041c at $0.8\,M_{\mathrm{Jup}}$. An alternative is the presence of a third, low-mass planet located between planets b and c, or on an inner orbit relative to planet b -- particularly near a mean-motion resonance with planet b -- which could account for the observed variations. These findings remain tentative, and further RV and photometric observations are essential to better constrain the mass of planet b and to refine the TTV modeling, thereby improving our understanding of the system's dynamical architecture.

The upcoming era of large scale, high cadence astronomical surveys demands efficient and robust methods for time series analysis. ARIMA models provide a versatile parametric description of stochastic variability in this context. However, their practical use is limited by the challenge of selecting optimal model orders while avoiding overfitting. We present a novel solution to this problem using a Bayesian framework for time series modelling in astronomy by combining Autoregressive Integrated Moving Average (ARIMA) models with the Nested Sampling algorithm. Our method yields Bayesian evidences for model comparison and also incorporates an intrinsic Occam's penalty for unnecessary model complexity. A vectorized ARIMA Nested Sampling framework is implemented allowing us to perform model selection across grids of Autoregressive (AR) and Moving Average (MA) orders, with efficient inference of selected model parameters. The method is validated on simulated and real astronomical time series, including the yearly sunspots number record, Kepler Lightcurves data of the red giant KIC 12008916, and TESS photometry of the exoplanet host star Ross 176. In all cases, the algorithm correctly identified the true or best-fitting model while simultaneously yielding well constrained posterior distributions for the model parameters. Our results demonstrate that Nested Sampling offers a potentially rigorous alternative to autoregressive model selection in astronomical time series analysis.

The current non-linearity correction for the Wide Field Camera 3 Infrared (WFC3/IR) channel is based on ground-based data acquired during WFC3's Thermal Vacuum 3 (TV3) testing campaign. In the current reference file, the correction coefficients derived for each pixel are averaged over each of the four detector quadrants. In this work, we compute a new pixel-based non-linearity correction using in-flight calibration observations with the internal tungsten lamp flats acquired between 2011 and 2013. We derive the new correction coefficients by fitting a third-order polynomial to the accumulated signal ``up-the-ramp" for each pixel. Approximately 2\% of IR detector pixels are flagged as bad, and a solution cannot be computed. For these, we use the quadrant averages of the new correction coefficients. An accompanying report (\cite{huynh2025}) provides detailed testing results using both internal flats and external science targets acquired in a variety of observing modes. The report highlights improvements in photometry derived from ``\this http URL" data products calibrated using the new reference file, with the largest improvement for pixels with fluence levels approaching the full well limit of $\sim$80,000 $e^- $. A new NLINFILE reference file was delivered in October 2025 and will be used to reprocess all WFC3/IR imaging and grism observations in the MAST archive.

Radio astronomy is facing critical challenges due to an ever-increasing human-made signal density filling up the radio spectrum. With the rise of satellites, mobile networks, and other wireless technologies, radio telescopes are struggling with radio frequency interference (RFI), which can masquerade, block or distort astronomical signals. In this chapter, we explain where RFI comes from, how it affects observations, and discuss different ways to reduce or remove interference. The techniques presented here reflect the state of the art in real-time RFI mitigation at the time of publication and include methods such as filtering, digital processing, and optimal scheduling. The proposed catalogue also explores new ideas like satellite avoidance through scheduling, the use of intelligent surfaces to block interference, and advanced computer algorithms to clean up data. The chapter also highlights the need for strong cooperation between astronomers and spectrum regulators to protect radio frequencies for future discoveries. By combining technical solutions and better policies, we can help ensure that radio astronomy continues to provide important insights into the universe.

Dark acoustic oscillations (DAOs) in the matter power spectrum can arise in many different dark sector models, and can imprint on a variety of cosmological observables. In this work we use measurements of the galactic UV luminosity function (UVLF) at high redshifts to constrain the dark acoustic oscillation feature at small scales in a model-agnostic way. We introduce a phenomenological transfer function model for a dark sector with a species undergoing DAOs which can accommodate sub-dominant dark matter abundances, and obtain constraints on its parameters. In order to predict the UVLF, we employ an Extended Press-Schechter formalism which we calibrate using N-body simulations with initial conditions featuring DAOs. Using measurements from the Hubble Space Telescope, James Webb Space Telescope, Subaru Telescope, and Canada-France-Hawaii Telescope, we constrain the wave number of the first DAO peak to be at $k \gtrsim 50\ h/\mathrm{Mpc}$, unless the fraction of dark matter undergoing DAOs is less than $0.07$.

We present NIRSpec/IFS observations of a rest-frame UV-bright, massive ($M_* \sim 10^{10}$ M$_\odot$, $z_{AB}=20.5$) galaxy highly magnified by gravitational-lensing observed just after the end of the epoch of reionization ($z=5.04$, $\bar{\mu}\sim90$). With JWST accessing the restframe UV and optical spectrum of this galaxy with high fidelity, we classify this UV-bright galaxy as post-starburst in nature -- due to weak/absent emission lines and strong absorption features -- making this an example of a new class of UV-bright but significantly quenched galaxies being discovered in this epoch. With a median $E(B-V)=0.44\pm0.14$, we identify the presence of stellar absorption across the arc both in Balmer lines and the MgII doublet, indicative of older stellar populations dominated by A stars (and potentially B stars). Using spatially-resolved maps of rest-optical strong emission lines, we find a heterogeneous distribution of nebular metallicities across the arc, potentially hinting at different enrichment processes. With a low median lensing-corrected H$\alpha$ star formation rate of SFR$_{H\alpha} = 0.024 \pm 0.001$ M$_\odot$ yr$^{-1}$, we find in the most "star-forming" clumps indications of lower ionization (log$_{10}$U $\sim -3.2$), lower nebular metallicities (12+log$_{10}$O/H $\lesssim$ 8.3), and hints of higher densities that suggest a possible recent infall of more pristine (low metallicity) gas onto the galaxy. Investigating the regions with no detectable H$\beta$ emission, we find (for the first time at $z>5$) signatures of diffuse ionized gas (DIG). Separating DIG from HII regions within a galaxy has predominantly been demonstrated at lower redshifts, where such spatial resolution allows clear separation of such regions -- highlighting the immense power of gravitational lensing to enable studies at the smallest spatial scales at cosmic dawn.

Detections of long-period giant exoplanets will expand dramatically with Gaia Data Release 4 (DR4), but interpreting these signals will require care. We derive the astrometric resoeccentric degeneracy: an astrometric analogue of the well-known radial velocity degeneracy in which a single eccentric planet can mimic two circular planets near a 2:1 period ratio. To first order in eccentricity, the sky-projected motion of a single eccentric orbit decomposes into a fundamental mode and first harmonic with an amplitude proportional to that eccentricity. A pair of coplanar, circular planets in a 2:1 orbital resonance produces the same harmonic structure: the outer planet sets the fundamental mode, while the inner planet supplies an apparent first harmonic. We present a mapping between the harmonic amplitudes and effective eccentricity ($e_\mathrm{eff}$) of a single planet that mimics a 2:1 configuration, demonstrating that $e_\mathrm{eff} = \, 2^{1/3}(M_{p,2}/M_{p,1})$, the masses of the inner and outer planets, respectively. Using simulated Gaia data we show that (1) coplanar 2:1 systems are statistically indistinguishable from a single eccentric planet and (2) mutual inclination can break this degeneracy. This bias favors detecting mutually inclined systems, often fingerprints of a dynamically hot history -- traces for processes such as planet-planet scattering or secular chaos. Determining the planetary architectures in which this degeneracy holds will be essential for measuring cool-giant occurrence rates with Gaia and for inferring their dynamical evolution histories.

We investigate a static, spherically symmetric black hole solution surrounded by a cloud of strings and immersed in an electromagnetic universe. By deriving the event horizon from the lapse function, we demonstrate that both the string cloud parameter and the electromagnetic background parameter significantly modify the horizon radius compared to the Schwarzschild case. Consequently, thermodynamic quantities-including the Hawking temperature, Bekenstein-Hawking entropy, and heat capacity-become explicit functions of these additional parameters, with the heat capacity exhibiting divergences that signal phase transitions. We analyze the motion of massive test particles in this spacetime, deriving the effective potential and calculating the innermost stable circular orbit radius, which governs the inner edge of accretion disks and influences orbital stability. Scalar perturbations are examined through the associated effective potential, and quasinormal mode frequencies are computed using the sixth-order WKB approximation; the negative imaginary parts confirm the stability of the black hole under such perturbations. We also study the photon sphere structure, black hole shadow radius, and photon trajectories, showing how the interplay between string clouds and the electromagnetic background shapes the optical properties of this spacetime. Finally, we investigate weak gravitational lensing phenomena by deriving the deflection angle for both massive particles and photons using the Gauss-Bonnet theorem applied to the optical geometry. The results exhibit notable deviations from the Schwarzschild geometry, with the string cloud enhancing the deflection through a $(1-\alpha)^{-1}$ factor while the electromagnetic parameter introduces competing corrections at second order.

We develop tools of Bayesian inference on the moduli space of Calabi--Yau (CY) manifolds. We sample from the invariant Weil--Petersson (WP) measure using Markov Chain Monte Carlo and normalising flows on \Kahler moduli space with dimension up to $h^{1,1}=30$, and present results on the spectrum of the CY volume and properties of divisors when the measure is restricted in physically meaningful ways. We furthermore present a theory-informed prior on axion masses and decay constants $(m_a,f_a)$ marginalised over the WP measure for all inequivalent CYs constructable from the Kreuzer--Skarke database with $h^{1,1}\leq 5$. We then impose likelihoods based on axion physics. We demonstrate how detection of a relatively heavy QCD axion at small $h^{1,1}$, e.g. by ADMX, provides detailed information about CY geometry and topology. Finally, we compute a full forward model incorporating likelihoods from the cosmic microwave background and Lyman-alpha forest and find the maximum posterior probability region on the moduli space of a given CY favoured by a resolution of the tension in these data by an ultralight axion composing $\mathcal{O}(1\%)$ of the dark matter. This demonstration serves as a blueprint for future statistical analyses within string phenomenology.

We show how a mass term for gravitational axions (''gravi-axions'') with a Chern-Simons coupling to gravity can naturally arise due to non-perturbative contributions from Euclidean wormholes, breaking the continuous shift symmetry of the standard theory. The induced mass can be generated in a cosmologically relevant range to be the dark matter or dark energy of the universe for a reasonable and well-motivated range of the symmetry breaking scale. Upon generating the gravi-axion mass term, we discuss the cosmology of the theory. We find that the gravi-axion can be produced to be the dominant dark matter component via misalignment or gravitational particle production, and that in a different regime, the gravi-axion can act as dynamical dark energy. We discuss gravi-axion decay into gravitons as a potential observational window for this theory. Finally, we find that, for late-time, astrophysical compact objects, cosmologically relevant gravi-axions behave as minimally-coupled, massive scalar fields.

Neutrino self-interaction beyond the Standard Model is well motivated by the nonzero masses of neutrinos, which are the only known particles guaranteed to have new physics. Cosmic messengers, especially neutrinos, play a central role in probing new physics, as they provide experimental conditions far beyond the reach of laboratories and serve as the link between laboratory fundamental-physics discoveries and their roles in the Universe, where many new physics motivations originate. In this work, we propose a novel probe of neutrino self-interactions through ultra-high-energy neutrinos scattering off the cosmic neutrino background when the lightest neutrino species remains relativistic today. This allows us to ``Widen the Resonance'' of such scattering. Meanwhile, we also provide a semi-analytic framework for cosmogenic UHE neutrino production, avoiding computationally intensive simulations and yielding results precise enough for BSM studies. The widened resonance enables future ultrahigh-energy neutrino telescopes, in particular GRAND, to probe mediator masses from MeV to GeV, reaching couplings down to $g \sim 10^{-3}$ -- up to two orders of magnitude beyond current bounds. Our results enhance the discovery potential of $\nu$SI in the high-mass regime, potentially offering crucial insights into the connections between the neutrino sector and dark sector.

Vickers et al. (2025) established that 58.20% of astrobiology experts believe intelligent extraterrestrial life likely exists, providing the first empirical baseline for public comparison. We surveyed 6,114 highly educated and scientifically engaged individuals (77.60% bachelor's degree+; 67.99% high-to-very-high scientific engagement) to assess their beliefs about extraterrestrial intelligent life existence: (1) personal beliefs, (2) perceived social circle beliefs, (3) perceived expert beliefs, and (4) responses to expert consensus revelation. Results showed 95.01% believed extraterrestrial intelligent life exists, with 62.59% holding definitive rather than probable convictions. Participants exhibited massive pluralistic ignorance, a 'cosmic closet', underestimating social circle beliefs by 46.07 percentage points despite near-universal personal conviction. Participants also exhibited a novel 'conviction intensity gap': while overestimating expert belief prevalence (67.63% vs. 58.20% actual), they underestimated expert conviction strength, perceiving only 21.10% as holding definitive beliefs. Experimental revelation of actual consensus (N = 5,106; 83.51% passed manipulation check) produced negligible personal belief change (d = -0.11) and small social belief change (d = 0.14). These findings demonstrate that consensus misperception operates along two dimensions, prevalence and intensity, and that even scientifically engaged audiences resist belief revision via expert consensus information.

For several years, students at an art college, working with NASA astronomers, have produced animations inspired by research on black holes, dark matter and more. They can be whimsical or poetic but still constrained by scientific rigour. The animations are used for scientific outreach and are freely available. Our program received a positive assessment through an evaluation we undertook. We are now planning a mobile STEAM exhibition to engage teenagers from underrepresented communities who may not typically consider STE(A)M for their studies. "Science anxiety" has been reported to be a significant barrier to learning. Mixing animation with astronomy can stimulate interest in STEAM, making science engaging in an unconventional way. One component would be activities where participants create artistic responses to astronomy. We undertook a workshop at a local city-run school, specialising in the arts for ages 14-17, to brainstorm the art/science activities. There we gave short scientific presentations leading to art activities: a giant colouring wall with projected celestial phenomena, a stop-motion station, and colouring images of comet 67P to produce an animation. Surveys before and after the activities showed positive responses. The hand of the artist has long been an important concept in animation (Crafton 1991). In a film entitled "The Movements of the Universe", this concept is adapted to the hands of scientists. Combining animation, filmed interviews at NASA (including a Nobel prize winner), and the scientists' hands, bring unexpected feelings of dream and humour to the audience. In this paper we explore three different viewpoints of these activities from a scientist, an animator, and an animation student.

While direct numerical simulations (DNS) are the most accurate method for studying turbulence, their large computational cost restricts their use to idealized configurations and to Reynolds numbers well below those found in practical systems. A recently proposed method, Large/Small Eddy Simulation (L/SES), aims to overcome this limitation while still providing the solution fidelity comparable to that of DNS. L/SES represents a pair of coupled calculations: a lower-fidelity Large Eddy Simulation (LES), which captures the large-scale flow structure, and a high-fidelity Small-Eddy Simulation (SES) targeting a sub-region of interest of the LES, in which the small-scale dynamics is fully resolved. In this study, we demonstrate the accuracy and performance of L/SES in large Reynolds-number homogeneous isotropic turbulence (HIT) up to Taylor-scale Reynolds number approximately 600. Turbulence properties obtained with L/SES are shown to be in close agreement with the literature, both in terms of global characteristics, such as kinetic energy spectra and dissipative anomaly, as well as small-scale properties, such as higher-order moments of the velocity gradients up to the 10th order and probability density functions of the intermittent quantities. Also using simulations of HIT, we systematically investigate various method parameters and determine their optimal converged values. Finally, we discuss the computational cost of L/SES and demonstrate that it is approximately 3 orders of magnitude lower than for a traditional DNS at the highest Reynolds number considered here. This highlights the potential of L/SES as a discovery tool, which brings high-fidelity simulations of realistic flows into the realm of feasibility.

We introduce a minimal Dark Standard Model (DSM) consisting of a single spin-0 particle with dark $U(1)$ gauge symmetry, and completely decoupled from the visible sector. Characterized only by the scalar mass $\mu$ and the dark charge $q$, this framework naturally gives rise to a rich phenomenology, including stable solitonic configurations that behave as dark "mini-MACHOs". We numerically build and evolve these gauged scalar-field solitons, derive their mass-radius relations, and identify a critical charge beyond which no gravitationally bound configurations exist. By combining these results with microlensing surveys that exclude compact objects heavier than the asteroid-mass scale ($M\lesssim 10^{-11}M_\odot$), we obtain the constraint $\mu\gtrsim 10\,\rm eV$ for viable configurations, depending on $q$. Our results represent a step forward in showing that purely gravitational observations can constrain the internal parameters of a dark gauge sector, and provide a framework for exploring broader DSM scenarios through future probes such as gravitational wave detections.

Source extraction is crucial in analyzing data from next-generation, large-scale sky surveys in radio bands, such as the Square Kilometre Array (SKA). Several source extraction programs, including SoFiA and Aegean, have been developed to address this challenge. However, finding optimal parameter configurations when applying these programs to real observations is non-trivial. For example, the outcomes of SoFiA intensely depend on several key parameters across its preconditioning, source-finding, and reliability-filtering modules. To address this issue, we propose a framework to automatically optimize these parameters using an AI agent based on a state-of-the-art reinforcement learning (RL) algorithm, i.e., Soft Actor-Critic (SAC). The SKA Science Data Challenge 2 (SDC2) dataset is utilized to assess the feasibility and reliability of this framework. The AI agent interacts with the environment by adjusting parameters based on the feedback from the SDC2 score defined by the SDC2 Team, progressively learning to select parameter sets that yield improved performance. After sufficient training, the AI agent can automatically identify an optimal parameter configuration that outperform the benchmark set by Team SoFiA within only 100 evaluation steps and with reduced time consumption. Our approach could address similar problems requiring complex parameter tuning, beyond radio band surveys and source extraction. Yet, high-quality training sets containing representative observations and catalogs of ground truth are essential.

Horizon-scale images of black holes provide a potential probe of fundamental physics, including tests of gravity and black hole hair. To assess the impact of scalar hair on accretion-flow imaging self-consistently, we construct an analytical model of a geometrically thick, magnetized disk around a rotating hairy black hole in Horndeski theory and analyze its 230 GHz image morphology. We find that scalar hair modestly alters the inflow and magnetic-field structure but strengthens gravitational redshift, markedly reducing the total flux and lensed ring brightness through relativistic transfer and spectral-shift effects. Moreover, we highlight a previously unexplored channel: the maximum interferometric diameter of the first photon ring responds strongly to the hair parameter but shows little dependence on accretion-flow details, making it a promising observable for constraining black-hole hair with future space-based interferometry.

The connection between turbulence and solar-wind acceleration, long known in space physics, is further developed in this Letter by establishing a robust empirical law that relates the bulk-flow speed to the magnetohydrodynamic-scale fluctuation energy in the plasma. The model is based on analysis of twenty-five years of near-Earth observations by NASA's Advanced Composition Explorer. It provides a simple way to estimate turbulence energy from low-resolution speed data -- a practical approach that may be of utility when high-resolution measurements or advanced turbulence models are unavailable. Potential heliospheric applications include space-weather forecasting operations, remote imaging datasets, and energetic-particle transport models that require turbulence amplitudes to specify diffusion parameters.

We propose a resonant annihilation as a way to reconcile the WIMP annihilation cross sections in a recently reported gamma ray signal from the Milky Way halo with that for the freeze-out and the upper limit from dwarf galaxies. We perform a simple model-independent analysis based on this hypothesis and determine the required parameters. We also present a simple particle-physics model that can accommodate them.

We present a dedicated search for gravitational-wave backgrounds with nontrivial peak structures using data from the first three and the initial part of the fourth observing runs of the LIGO-Virgo-KAGRA network. The analysis is motivated by a variety of early-Universe models characterized by signals with multiple peaks. We introduce a model independent parameterization of double-peaked spectra based on the superposition of two normalized broken power laws and perform a Bayesian inference study using the LIGO-Virgo-KAGRA isotropic cross-correlation data. While no statistically significant evidence for a multi-peak background is found, the analysis provides constraints on the inter-peak slopes in correlation with the signal amplitude. These results exhibit LIGO-Virgo-KAGRA's ability to probe signals beyond a single peak structure and establish a foundation for future targeted searches for nontrivial gravitational waves background spectral shapes in future observing runs and the advanced detector era.

We take a fresh look at the viability of physically realistic extremal black holes within our (non-supersymmetric) low energy physics. By incorporating prefactors and volume effects, we show that Schwinger discharge in charge neutral environments is far more efficient than commonly assumed. Using ionization estimates for neutral hydrogen, we obtain a new and robust lower bound on the mass of an extremal electrically charged black hole, exceeding $10^{14} M_\odot$. For magnetic black holes, we compute the Lee-Nair-Weinberg instability and revisit early universe pair creation rates, including singular instantons that substantially enhance production, to demonstrate that the extreme charges required for stability are cosmologically implausible. Finally, we suggest that an extremal Kerr black hole could shed angular momentum via superradiant scattering from the stochastic gravitational wave background. Taken together, our results provide a unified picture that extremal black holes of any type are unlikely to persist in our universe.

Tidal Love numbers of black holes, zero in classical general relativity for Kerr black holes in vacuum, become non-vanishing in the presence of exotic matter or in alternative theories of gravity, making them a powerful probe of fundamental physics. The gravitational-wave event GW250114, observed with an unprecedented signal-to-noise ratio, provides a unique opportunity to test this prediction. By analyzing this event, we conclude that the data is consistent with the binary black hole hypothesis, and we place a 90\% upper limit on the effective tidal deformability of $\tilde{\Lambda} < 34.8$. These bounds imply that any environment surrounding the black holes must contribute less than $\sim 7\times 10^{-3}$ of their mass, and they rule out some models of boson stars. Our findings provide the strongest observational constraints yet on black hole tidal deformability and show that the data remain fully consistent with the Kerr black hole prediction of vanishing tidal Love numbers.

We investigate the quantum-to-classical transition of primordial perturbations within a two-field inflationary framework where an adiabatic mode interacts with an entropic environment. In the case of a massive entropic environment, the attractor Bunch--Davies vacuum plays a special role: it is the only state that can undergo full recoherence, whereas all excited initial states exhibit persistent loss of purity. To characterize this behavior, we parameterize excited Gaussian initial states by their Bogoliubov coefficients and compute the purity and Rényi-2 entropy of the reduced adiabatic state as information-theoretic indicators of decoherence dynamics. We find that excited states display \emph{purity-freezing} at a non-zero plateau, where residual quantum correlations persist indefinitely, a qualitative departure from the complete recoherence observed for the Bunch--Davies vacuum. This sensitivity to initial conditions highlights the non-generic nature of full recoherence in the quantum-to-classical transition of inflationary perturbations.