Papers for Tuesday, May 19 2026

We describe Mosaic, a computational geometry method for constructing the surface mosaic induced when a Cartesian volume grid intersects a spherical shell. The motivating application is conservative coupling between data produced on rectangular grids and diagnostics or boundary conditions defined on spherical surfaces, as occurs in space-weather, magnetohydrodynamic, atmospheric, and geophysical models. The method identifies Cartesian cells that intersect the shell, constructs cell-sphere prepatches, splices those regions by the spherical colatitude grid, and then splices by azimuth to produce final patches indexed by the five-tuple (nx, ny, nz, ntheta, nphi). The implementation explicitly treats the polar singularity by separating polar-derived theta patches from ordinary phi splicing. A near-pole numerical failure mode, caused by linear interpolation in azimuth, is removed by computing exact intersections between great-circle boundary segments and meridian planes. The prepatch construction also handles several degeneracy cases that occur beyond the ordinary corner-straddling geometry, including doubly-crossing edges, lens-shaped prepatches, secondary closed loops, and re-entrant face arcs. For a representative nonuniform Cartesian test grid, the current implementation builds 3618 intersecting Cartesian cells, 3602 ordinary prepatches, 6476 theta patches, and 9714 final phi patches, with all 3602 ordinary cells built successfully and zero theta or phi splicing failures. The final phi-spliced patches close in normalized area to roundoff relative to their theta-spliced parents. The method is implemented in the Java/Maven application mdi-mosaic, which provides visualization, diagnostics, mouse-over patch inspection, and JSON export of final patch boundaries and normalized areas.

GW231123 is a short, massive binary-black-hole event whose source properties show strong waveform dependence. We use this event to test gravitational-wave polarization birefringence, modeled as a frequency-dependent rotation of the tensor-polarization basis. Instead of sampling a distance-normalized coefficient directly, we sample the band-differential rotation $\delta_{\rm br}=\Delta(448\,\mathrm{Hz})-\Delta(20\,\mathrm{Hz})$ with prior $[-\pi,\pi]$, and report the derived coefficient $\beta_{\rm br}^{\rm derived}$ for comparison with standard propagation parametrizations. We analyze three waveform families: IMRPhenomXPHM (XPHM), IMRPhenomXO4a (XO4a), and NRSur7dq4. The derived posteriors are consistent with the general relativity value, giving $90\%$ upper limits $|\beta_{\rm br}^{\rm derived}|_{90}=0.378,\,0.097,\,0.273$ for XPHM, XO4a, and NRSur7dq4, respectively. The directly sampled $\delta_{\rm br}$ posterior remains broad, with $|\delta_{\rm br}|_{90}\simeq2.8\,\mathrm{rad}$, so the accumulated rotation across the analysis band is weakly constrained. The Bayes factors are waveform dependent: $\ln\mathcal{B}_{\rm br/GR}=-1.26\pm0.30$, $+3.64\pm0.28$, and $-0.86\pm0.29$, respectively. We therefore find no waveform-robust evidence for parity-violating propagation. The positive XO4a result is better interpreted as a waveform-dependent birefringence-like response associated with the mass-ratio--distance--spin degeneracy of this short high-mass event.

Redshift drift is usually discussed for observers comoving with the cosmological background, but realistic observations are made by observers with nonzero peculiar motion. In this work, we calculate the expected redshift drift for tilted observers within the covariant \(1+3\) formalism. Starting from the exact redshift measured in the tilted frame, we derive the corresponding drift as an FLRW background term plus a directional correction driven by the observer's peculiar kinematics, encoded through peculiar expansion, projected shear, and projected acceleration along the line of sight. We analyse first an Einstein--de Sitter (EdS) background, which isolates the purely kinematic effect of tilt in the absence of background acceleration, and then extend the calculation to \(\Lambda\)CDM in order to quantify how the same anisotropic corrections deform the standard drift signal in the concordance model.

We use FLAMINGO to investigate why feedback from active galactic nuclei (AGN) significantly depletes gas in galaxy groups but is ineffective in clusters. We delineate three radial zones: an inner zone where AGN feedback heats halo gas via shocks; an intermediate buoyancy zone where the heated halo gas rises; and an outer zone where the outflow may stall in a termination shock. Heating in the inner zone self-limits because, once the gas is sufficiently hot, shocks become too weak to deposit further entropy. Consequently, outflows have a ceiling entropy value ($360\, {\rm keV\, cm^2}$) that is nearly independent of halo mass. These values (and trends with redshift and feedback variants) are explained using an argument based on the Rankine-Hugoniot relations. Outflows rise at fixed entropy through the buoyancy zone, escaping the halo if the ceiling value is sufficiently elevated over that of the inflowing gas. This condition is satisfied only for halo masses $M_{\rm 200m}<10^{13.7}\,{\rm M_\odot}$, because inflow entropy tracks the virial relation. Variants with stronger (or weaker) feedback have a higher (or lower) entropy ceiling and a correspondingly modified critical mass of $M_{\rm 200m}=10^{14.0}\,{\rm M_\odot}$ (or $10^{13.5}\,{\rm M_\odot}$). In clusters above the critical mass, the increased inflow entropy causes the outflow to stall and potentially shock at the 'splashback' radius. We derive an expression for the time evolution of the virial gas fraction, which shows how lingering gas is reincorporated as the halo virial radius expands. This effect dominates over outflows unless they rejoin the Hubble flow; as a result, virial gas fractions rise as a function of mass starting at $M_{\rm 200m} = 10^{13.0}\,{\rm M_\odot}$. These effects explain why groups have depleted gas, while clusters have close to the cosmic baryon fraction.

We present a theoretical model in which the recently detected over-massive black holes (OBHs), and possibly Little Red Dots (LRDs), arise during a halo-driven transient phase preceding the established coevolution of supermassive black holes (SMBHs) and their host galaxies. In this model, halo gravity drives an early phase of rapid black hole growth, leading to systems in high-redshift haloes that lie above the local scaling relations. As the halo evolves, a transition in halo thermodynamics leads to the onset of a hot, pressure-supported medium that suppresses accretion, reducing the black hole growth rate and driving the system toward the local black hole mass$-$stellar mass relation. LRDs may represent an observational manifestation of the rapid, halo-driven growth phase, while OBHs trace its direct mass signature. Our model thus provides a unified framework in which these systems form and evolve toward the regulated coevolution observed in the local Universe.

Stars in the initial and carbon-oxygen core mass ranges of $\sim140-260$ and $50-130$ M$_\odot$, respectively, with low metallicity are predicted to experience copious electron-positron pair production in their cores, leading to a runaway thermonuclear explosion that obliterates the entire star in a luminous and long-duration pair-instability supernova explosion. Some previous supernovae have been interpreted in this context but lack the full range of predicted properties. Here, we report detailed observations and modeling of the hydrogen-rich supernova 2023vbw, which exploded in a low-metallicity ($\sim0.1$ Z$_\odot$) environment in a dwarf star-forming galaxy at a redshift of $0.088$. Its light curve exhibits a luminous ($1.6\times10^{43}$ erg s$^{-1}$) and long-duration ($190$ days) main peak, resulting in a total radiated energy of $3\times10^{50}$ erg, more than an order of magnitude greater than canonical core-collapse supernovae. Semi-analytical light-curve modeling yields a blue supergiant-like progenitor with an ejecta mass of $170-350$ M$_\odot$, radioactive nickel mass of $1.2-1.6$ M$_\odot$, and explosion energy of $(6-13)\times10^{52}$ erg, well matched by pair-instability models. The early and late-phase light curve and spectra also show evidence for interaction of the supernova ejecta with an aspherical circumstellar medium. Discoveries of numerous such events with the upcoming Rubin Observatory and Roman Space Telescope will shed light on the deaths of the most massive stars in the Universe.

We propose that black hole feedback regulates gas in massive halos by establishing an entropy ceiling; the resulting buoyant gas migrates to the virial radius with no additional energy input required. The FLAMINGO simulations support this picture: at the virial radius, outflow entropy is mass-independent for isotropic thermal feedback but depends on the solid angle of directly heated gas for jet feedback. Above a critical halo mass $M_\rm{crit} \approx 10^{13.5\text{--}14}\, M_\odot$, virial shocks overwhelm the ceiling, predicting rejuvenation of star formation in the most massive galaxies, supported by new low-redshift evidence from star formation rates and morphologies.

Wide-angle outflows, or winds, are associated with a broad range of astrophysical systems, including protostars, massive stars, X-ray binaries, tidal disruption events (TDEs), luminous fast blue optical transients (LFBOTs), and starburst galaxies. When these winds first ``turn on," they inflate a ``bubble" into their surroundings, bounded by two shocks and a contact discontinuity, and evolve through distinct adiabatic phases prior to the onset of significant radiative cooling. For sufficiently overdense ejecta, the flow quickly relaxes into an interaction-dominated similarity state at early times and later enters an energy-conserving regime. We present a systematic study of these phases for adiabatic winds expanding into power-law density profiles $\rho \propto r^{-n}$ with $0 \leq n \leq 2$. Using analytic scalings together with one-dimensional shock-capturing hydrodynamic simulations, we quantify both the relaxation timescales and the accuracy with which the corresponding similarity solutions reproduce the fluid velocity, density, and pressure throughout the shocked bubble. We show that the interaction solutions are attained within only a few dynamical times and remain valid until the reverse-shocked shell is no longer thin relative to the forward-shocked shell, corresponding in practice to an instantaneous overdensity of order unity. For $n < 2$, the flow subsequently converges to the generalized energy-conserving scaling $R_s \propto t^{3/(5-n)}$, while the special case $n=2$ exhibits a single persistent similarity state. We discuss the durations and implications of these phases for stellar and galactic outflows, TDEs, and LFBOTs.

Detecting light reflected off the dayside of an exoplanet in high-resolution spectroscopic data has proved to be a notoriously difficult endeavour. Despite several attempts, the faint signal has yet to be detected. We present a new effort at finding reflection signatures and show how a strong rotational broadening of the reflected spectrum can complicate this objective. We introduce a new figure of merit that quantifies the favourability of different systems for a reflection study, the reflection spectroscopy metric. Applying this metric, we identify the KELT-9 system, which features a highly misaligned, rapidly rotating host star, as the target for a case study based on a spectroscopic time series obtained by CARMENES. We also perform an injection-recovery test to determine the detectability of the signal in our data and demonstrate its sensitivity to rotational line broadening. The search for a genuine reflection signal in our data resulted in a non-detection. The injection-recovery test puts this finding into context by revealing the critical importance of taking rotational broadening into account when dealing with systems featuring rapidly rotating stars and large spin-orbit misalignments. The case study presented here underscores the need to incorporate stellar rotation and spin-orbit misalignment into assessments of a given planet's favourability to reflection studies.

The latest Gaia data release provides low-resolution spectra for approximately 100 000 white dwarfs. Though useful for pre-classification, they lack the resolution required for accurate spectral type and parameter determination, motivating spectroscopic follow-up campaigns. In this work, we assess the reliability of machine-learning spectral classifications derived from Gaia spectra through comparison with medium-resolution spectroscopy, determine the nature of objects classified as "massive helium-rich (DB)" by automated methods, and characterise the properties of warm and hot DQ (carbon-dominated) white dwarfs, magnetic and metal-polluted objects. To do this, we observed 255 white dwarfs with the Gran Telescopio Canarias equipped with the OSIRIS instrument (R ~ 1000). Spectral types were assigned through visual inspection and compared with machine-learning classifications applied to Gaia spectra. Magnetic objects were identified via Zeeman splitting, and magnetic field strengths were estimated. We find machine-learning classifications are highly accurate (> 90% for spectral types in their training sets), despite the low resolution of Gaia spectra. We show "massive DBs" to be mostly magnetic white dwarfs and warm DQs, with only 5 of 112 observed (4.46%) confirmed as genuine DBs. Warm DQs are found along the Gaia Q branch and exhibit unusually high tangential speeds. We provide spectral classifications for 255 white dwarfs, demonstrate that Random Forest algorithms reliably classify low-resolution Gaia spectra into main spectral types, determine the nature of "massive DBs", and identify a large population of magnetic white dwarfs and carbon-rich objects. Several rare subtypes are identified, including 1 DAQ, 1 DQZA, 4 hot, 29 warm DQ stars, and 63 magnetic white dwarfs. The properties of warm DQs are consistent with previous studies, supporting their proposed origin as merger remnants.

We present late-time ($\delta t > 3000$ d) optical (Keck), X-ray (Chandra and NuSTAR), and radio (VLA, ALMA, and the uGMRT) observations of the Type Ic-BL SN 2012ap. Previous studies of this SN suggested that it stands out as a key example of a weak engine-driven explosion due to the lack of gamma-ray burst detection and a mildly relativistic ejecta. Recently, radio sky surveys revealed the rebrightening of the radio emission from this SN, highlighting the possibilities of a density enhancement at large radii or the existence of an off-axis relativistic jet. While the late-time optical spectra does not exhibit the broad emission lines seen in other interacting SNe, our analysis of the broadband radio and X-ray emission implies that both scenarios are plausible. If a density enhancement is responsible for the radio rebrightening, it has to result from a change in the mass-loss rate and/or wind velocity, possibly due to the transition of the progenitor from a red supergiant to a Wolf-Rayet star. If the late-time radio component is a result of an off-axis relativistic jet, we find that an energetic narrow jet viewed at $\theta_{\rm obs} \geq 80^{\circ}$ is needed. In this scenario, SN 2012ap is not a result of a weak engine-driven explosion, and, instead, it is similar to other GRBs. However, radio rebrightenings of Type Ic-BL SNe are not enough on their own to determine the existence of off-axis jets and our planned VLBA observation will help reveal the true nature of this SN.

We present a comprehensive characterisation of close stellar encounters in the solar vicinity, with a particular focus on placing the predicted fly-by of GJ~710 in context. This star will come extremely close ($0.0621$\,pc or $\sim10^4$\,AU) to the Solar System in approximately 1.3\,Myr. Using a linear motion approximation, we identified past and future close stellar encounters within 1\,pc of the Solar System, using a complete sample of nearby stars. We assessed the completeness of our dataset and applied corrections to the radial velocities, accounting for gravitational redshift and convective blueshift. Such effects can bias the measured velocities and affect the derived encounter parameters. Furthermore, we computed close encounters for all the stars in the Solar System vicinity to build a statistically significant sample of such events. We accounted for binary systems and common proper motion pairs, applying corrections to account for incompleteness at the edges of our time window. We derived reliable statistics for close stellar encounters of stars within 25\,pc of the Sun. We report a rate of encounters within 1\,pc and within 0.47\,Myr of $10.6 \pm 4.5$ per Myr and star, implying an average of one encounter every $95^{+71}_{-28}$\,kyr. Additionally, from the resulting distribution, we can evaluate the relative frequency of the upcoming GJ~710 fly-by, which is found to be quite rare, at a rate of one per $\sim50$\,Myr. This work provides new insights into the dynamic environment of the solar neighbourhood by quantifying the rates and distributions of close stellar encounter events.

Hierarchical triple evolution provides a promising alternative to isolated binary formation models for black holes (BHs) and neutron stars (NSs) with low-mass stellar companions. To search for tertiaries, we perform deep, adaptive optics-assisted, near-infrared imaging of five quiescent BH low-mass X-ray binaries (LMXBs), Gaia BH1, and twelve Gaia NSs. We detect several faint stars previously unresolved in survey imaging, but none are close enough to robustly rule out a chance alignment. To achieve high contrast sensitivity at close separations, we use the reference star differential imaging strategy with the Karhunen-Loéve Image Processing algorithm to model and subtract the point-spread function of each target. We identify tertiary candidates in the speckle-dominated regime, but injection-recovery tests suggest most 5$\sigma$ detections are likely artifacts. We derive $5\sigma$ contrast curves and convert these to limits on the mass of main sequence (MS) tertiaries and the effective temperature of white dwarf (WD) tertiaries consistent with a non-detection. We rule out plausible MS tertiaries and young, hot WD tertiaries at projected separations $\gtrsim 500$ au for the Gaia compact object binaries and $\gtrsim 2000$ au for the more distant BH LMXBs. While the recent discovery of a $1.2\,M_{\odot}$ tertiary to V404 Cygni supports triple formation scenarios for BH LMXBs, our results suggest such companions are relatively rare. Our observations remain consistent with intermediate-mass tertiaries that have since evolved into cool WDs, detectable with deeper JWST imaging. Follow-up observations are required to measure proper motions and confirm or rule out physical association of tertiary candidates.

We investigate the origin of the scatter in the supermassive black hole (BH) masses for galaxies in the L* regime, using the ARTEMIS and EAGLE simulations. By classifying galaxies based on their central BH / galaxy stellar masses ratios, we follow the evolution of galaxies with the lowest and highest such ratios (denoted LBH and HBH, respectively). We find that the properties of these two galaxy samples are comparable at z ~ 2 but diverge significantly towards lower redshifts. Galaxies with less massive BHs were able to maintain higher gas fractions and sustained star formation during their evolution, whereas those with more massive BHs formed stars earlier, grew BHs faster and experienced more efficient feedback and subsequent quenching. The simulations broadly match the observed scatter in the BH masses and galaxy morphologies in the L* regime and explain the origin of this scatter in terms of differences between merger histories. Galaxies with more active merger histories contain more massive BHs at present time and tend to be elliptical, while galaxies with more quiescent histories have less massive BHs and tend to be disc-like. Mergers enhance BH growth through higher gas accretion rates onto central regions and direct BH-BH mergers. However, these channels operate differently: in HBH galaxies, BHs grow primarily (~ 90%) by gas accretion, whereas in LBH they grow both through gas accretion (~ 60%) and BH-BH mergers (~40%). Our results suggest that the different BH masses in MW and M31 could be explained by differences in merger histories.

The relative rate of neutron stars and black holes produced by the collapse of massive stars is highly uncertain. We simulate the stellar collapse of $195$ progenitors with masses between $9\, M_\odot$ and $120\, M_\odot$, incorporating a schematic treatment of neutrino flavor conversion. We find that flavor transformation reshapes the explodability of massive stars-especially in the $16$-$30\, M_\odot$ mass range-and modifies the compact remnant mass distribution. Our findings identify neutrino flavor conversion as a fundamental ingredient in predicting neutron star and black hole populations, while naturally easing the red-supergiant and the supernova-rate problems, as well as reconciling theoretical expectations with the low-mass tail of the observed neutron star mass distribution.

The tight empirical correlation linking the stellar disk scale length $R_d$ to the dark matter scale radius $r_0$ has been proposed as possible evidence for a fundamental coupling between baryons and dark matter beyond gravity. We re-examine the physical origin of this relation using a sample of 31 galaxies from the NIHAO cosmological hydrodynamical simulations, which include no dark matter-baryon interactions beyond gravity and baryonic feedback processes. NIHAO naturally reproduces both the normalization and the small scatter of the observed $R_d-r_0$ relation at $z=0$, while yielding a slightly shallower slope. By tracking galaxies from $z=2$ to $z=0$, we identify three distinct evolutionary classes: systems undergoing disk expansion, contraction, and quasi-static evolution. Using a Bayesian hierarchical framework, we provide the first characterization of the cosmic evolution of the $R_d-r_0$ relation, tracing the evolution of its normalization, slope, and intrinsic scatter from $z=2$ to the present day. We find a mild decrease in normalization ($\sim0.07$ dex), a flattening of the slope from $\alpha \simeq 1.05$ to $\alpha \simeq 0.95$, and a weak decline in the intrinsic scatter toward lower redshift, suggesting that galaxies evolve preferentially along the relation while jointly re-balancing their stellar and dark matter scales. By comparing hydrodynamical simulations with their dark-matter-only counterparts, we isolate the impact of baryons and baryonic feedback on halo structure. Our results show that stellar feedback alone can reshape the central potential and naturally establish the observed coupling between luminous and dark matter, without requiring modifications to the dark sector.

Gravitational-wave observations are playing an instrumental role in understanding the population of binary compact objects in the Universe. These observations have begun to hint at the mass distribution of binary black holes (BBHs), with tentative evidence for features in the mass distribution beyond a simple power-law. Such features, hence, can be connected with different formation scenarios of BBHs and lead to important astrophysical conclusions. However, it is crucial to understand whether these features are truly astrophysical or connected with any unknown systematics. We show in this work that waveform modelling uncertainties can significantly distort inferred features in the BBH mass distribution, which can be more pronounced than the statistical uncertainty, even with the current generation detectors, which can peak close to the lower edge of the pair instability supernovae (PISN) mass gap, and also can impact the slope of the power-law distribution. So, in order to have a confirmed detection of astrophysical features in the BBH mass distribution and connecting them with BBH formation channels, it is important to consider waveform systematics in the astrophysical population analysis. We show the typical scaling of the systematic error and discuss a few avenues to mitigate this effect for robust measurements in the future.

In this work we investigated methods for the accurate and efficient incorporation of photometrically classified supernovae into cosmological analyses, and to assess the impact of the additional uncertainty associated with this procedure on the ability of Type Ia supernovae (SNeIa) tests to place constraints on cosmological models. We proposed a simplified likelihood, in which the contamination is described as a redshift dependent change in the mean of the usually assumed Gaussian distribution, and we tested this hypothesis against the usual two-component approach, based on the BEAMS framework. Using the latest version of the DES supernova sample, dubbed DES-Dovekie, we compared the results when using type probabilities from different classifiers, such as SNIRF and SCONE, and applying different cuts on these probabilities. We show that the new model is strongly favored by the Bayes factor, when compared with the current one, for all configurations, allowing an improvement on the constraining power of photometric supernova data.

Context. Type Ibn supernovae (SNe) are characterized by narrow helium emission lines arising from ejecta-circumstellar medium interaction, yet their progenitors remain debated, with both massive Wolf-Rayet stars and low-mass helium stars in binaries proposed. Aims. We aim to directly identify the progenitor of the Type Ibn SN 2023fyq and to characterize its environment in order to constrain the progenitor's nature and evolutionary channel. Methods. We search for the SN progenitor based on pre-explosion and late-time HST and JWST images and derive its properties by fitting the spectral energy distribution. We investigate the SN environment by probing the stars, dust, ionized gas and molecular gas with a multi-wavelength dataset including HST and JWST imaging, VLT/MUSE integral-field-unit spectroscopy and ALMA CO (2--1) radio interferometry. Results. We discover a pre-explosion source at the SN position, which is consistent with a hot ($T>$15000 K) and luminous (log($L$/$L_\odot$) $\gtrsim$ 5.5) SN progenitor and a possible host star cluster. The progenitor is confirmed to have disappeared after explosion. Analysis of the SN environment implies that the progenitor likely has an age of log($t$/yr) = 7.1--7.2. These phenomena disfavor a very massive single-star progenitor and instead support a binary scenario involving a low-mass helium star and a compact object; the observed progenitor emission likely arises from binary interaction that began at least $\sim$12 yr before the explosion. Conclusions. SN 2023fyq is the first Type Ibn SN with a directly detected progenitor and a possible host star cluster. It adds to the diversity of Type Ibn SNe in terms of their progenitor channels and mass-loss mechanisms.

We present VLA observations of 16 Type IIn and Type II-L supernovae (SNe IIn and SNe II-L) at ~1000--7000 days after explosion, probing circumstellar matter (CSM) at distances >10^16 cm from the progenitor corresponding to mass-loss over hundreds to thousands of years before core collapse. We detect radio emission from four SNe (1998S, 2005ip, 2008fq, and PTF11iqb) with the remaining 12 yielding upper limits of nu L_nu < 10^35--10^36 erg s^-1 at 3--11 GHz. The detected sources span approximately two orders of magnitude in radio luminosity, reflecting a wide range of CSM densities and pre-explosion mass-loss histories. All detected sources exhibit steep spectral indices (~<-0.4) consistent with optically-thin synchrotron emission, and the spectral evolution supports internal free-free absorption as the dominant absorption mechanism at these late epochs. We infer progenitor mass-loss rates of \dot{M}/v_w ~< 10^-5--10^-3 Msun yr^-1/(1000 km s^-1), with the most radio-luminous objects requiring sustained mass-loss over hundreds to thousands of years. The detection of the intermediate SN IIn/SN II-L object PTF11iqb at luminosities between classical SNe IIn and SNe II-L supports a continuum between these subtypes in terms of CSM interaction strength. Our limits further suggest that SNe IIn and SNe II-L are not separated by long-term mass-loss rate at the radii probed here, but chiefly by the presence and strength of dense circumstellar material immediately before explosion. At epochs >5000 days, some SNe (e.g., 1979C and 1986J) maintain nearly constant radio luminosity while others decline rapidly, suggesting that the most radio-luminous SNe IIn arise from progenitors with sustained mass-loss extending >10^4 yr before explosion.

The sky-averaged redshifted 21-cm signal from Cosmic Dawn is expected to provide a unique view of the first compact objects. However, its measurement remains daunting. Difficulties are driven by the large dynamic contrast between the intervening foregrounds and the signal-of-interest, which places extremely high demands on instrumental calibration and data quality measures. The ongoing debate within the field concerning the evidence of a potential first detection by the EDGES experiment highlights the need for a more robust set of analysis methods and tools that are reliable and accessible. In this paper, we detail for the first time the precise calibration and analysis methodology adopted in previous EDGES data releases. These methods are presented in the context of a new open-source end-to-end analysis and simulation package for 21-cm global signal experiments that both formalizes these methods and provides general tools for the broader community. Finally, we describe the raw data used in previous EDGES papers and release these data publicly for extended scrutiny.

Most main sequence stars, unlike our Sun, belong to multiple systems with two or more stars. How and when these multiples come together and become bound is uncertain, since the earliest stages of star formation are difficult to resolve. We analyze simulations of star cluster formation in Milky Way-like conditions, including all key physics and stellar feedback mechanisms, to understand how multiple systems form. We show that $\approx 70-80\%$ of binaries are gravitationally bound from the moment the second star forms. Binaries evolve and accrete together, which will affect their planetary systems and chemical evolution. Half of the binaries are disrupted by the end of the star-formation epoch, such that $\approx40\%$ of the final single stars belonged to a multiple at some point, with implications for the stellar initial mass function. Formation in multiples is the dominant mode of star formation, accounting for at least 57% of stars.

Accreting millisecond X-ray pulsars (AMXPs) are transient X-ray sources likely to be in the final evolutionary phases of low-mass X-ray binaries (LMXBs). In this work, we have investigated the torque and X-ray luminosity variations of five AMXPs during outburst and quiescence, using a model previously employed to explain some typical behaviours of LMXBs. Most of these systems spin up in the outburst state and spin down in quiescence, while they slow down on the long term. We have obtained reasonable results with the model in agreement with these observations. We have found that the torques produced by the disc-magnetosphere interaction, the magnetic dipole radiation, and the mass accretion are compatible with the X-ray luminosity and rotational properties of the sources in their outburst and quiescent states, and during the resultant long-term evolution. Nevertheless, our results are not sufficient to rule out the spin-down contribution of the gravitational radiation torques due to significant timing noise and uncertainties about the bolometric corrections for X-ray luminosity during the outbursts of some sources.

The first detected member of a new astronomical class is often not representative of the underlying population, but instead reflects the selection effects of the observing technique that found it. We apply this idea to the first remote detection of biosignatures with two leading near future strategies: JWST transmission spectroscopy and HWO reflected light direct imaging. Using the known signal scalings of the two methods together with a simple detectability model, we show how a rare but observationally favored planet class can dominate early detections even when it is intrinsically uncommon. For JWST, an early biosignature detection is most likely to arise from a detectability favored outlier, such as a sub-Neptune or other atmosphere rich planet around a nearby M dwarf, rather than from a true Earth analog. For HWO, the situation is subtler. Among accessible habitable-zone targets around FGK-type stars, differences in maximum observable distance and hence in effective survey volume may be smaller than in the JWST case, weakening the volume bias. At the same time, stellar-type-dependent photochemistry can alter biosignature abundances, so the first HWO biosignature may emerge from a balance between photochemical enhancement and geometric accessibility. Nevertheless, within the accessible sample, planets with stronger biosignature features and higher reflected light contrast may still be favored in early detections. A first HWO biosignature could be a selection favored outlier and should not be assumed to represent inhabited rocky planets in general. Crucially, the longest lived biosphere on a planet is not necessarily its most spectrally detectable one. If the first detection turns out to be an outlier, that may still suggest that a more broader range of habitable environments awaits discovery.

We investigate how cosmic ray (CR) transport in molecular clouds and their substructures can be probed using multi-wavelength observations. The detailed microphysics regulating the penetration and coupling of CRs in dense molecular structures is unsettled. Self-generated turbulence can produce scattering and diffusive transport, while ion-neutral damping in cold, dense gas promotes ballistic CR propagation. We construct a self-consistent framework for CR transport and interactions in magnetized molecular clouds, considering three limiting propagation scenarios: ballistic transport, diffusion, and a hybrid configuration with a diffusive envelope and quasi-ballistic core. By forward-modeling pion-decay $\gamma$-ray emissivities, CR-driven ionization-rate profiles, and electron synchrotron emission in the hard X-ray band, we connect GeV attenuation and propagation signatures to independent diagnostics of secondary production and CR penetration. As an illustrative example, we apply our framework to the Taurus molecular cloud complex and selected embedded clumps. We show that CR scattering may be substantially enhanced on clump scales, with inferred CR diffusion coefficients suppressed relative to canonical interstellar medium (ISM) values at GeV energies. In this interpretation, CRs are more closely coupled with dense gas in the ISM, and a diffusive envelope boosts the effective gas column density encountered by the CRs. This increases the hadronic interaction rate in the cloud. In turn, the secondary CR electron injection is also increased, and CR ionization rates are elevated at higher densities. We show that a hard X-ray synchrotron emission component is also generated, which may be detectable with near-future facilities. Finally, we discuss how future $\gamma$-ray, X-ray, and ionization constraints will provide firm tests of CR propagation theories in molecular cloud environments.

Helium stars stripped of their hydrogen envelopes represent pivotal phases in binary evolution, yet their origins, particularly within the intermediate-mass range of $2-8\, M_{\odot}$, still remain poorly understood. This population bridges the gap between low-mass hot subdwarfs and massive Wolf-Rayet stars, but has remained largely unobserved. In this study, we employ binary population synthesis to systematically investigate the formation and properties of intermediate-mass helium stars (IMHeS) across various galactic metallicities. Our results indicate that metallicity and common-envelope ejection efficiency are the dominant factors shaping the IMHeS population. We estimate that several thousand IMHeS exist in the Milky Way, with several hundred more in the Magellanic Clouds. The vast majority of IMHeS reside in binaries, with fewer than $10\%$ appearing as single stars. Among IMHeS binaries, $\gtrsim 50\%$ are expected to have main-sequence companions, and the remainder host compact companions (including helium stars, white dwarfs, neutron stars, or black holes). The former systems form mainly through stable mass transfer, whereas the latter arise predominantly from common envelope evolution. Our work provides quantitative predictions for the populations of these elusive stars formed through binary interactions and offers guidance for future observational searches.

Wide-field slitless spectroscopic galaxy surveys, such as the one performed by the upcoming Chinese Space Station Survey Telescope (CSST), are crucial for precision cosmology but present formidable data analysis challenges. Because spectra are dispersed directly onto the detector, they are convolved with the 2-dimensional (2D) spatial morphology, which complicates wavelength calibration and consequently degrades the fidelity of subsequent 1-dimensional (1D) spectral extraction. To overcome these limitations, we present a deep learning framework that extracts redshifts directly from 2D slitless spectral images, bypassing 1D extraction entirely. We construct a realistic mock dataset for the CSST $GV$ and $GI$ band using high-resolution images from HSC-SSP PDR3 and spectral energy distributions (SEDs) from DESI DR1. A Bayesian convolutional neural network implemented by Monte Carlo dropout is employed to map the 2D spectral images to redshift estimations while simultaneously quantifying uncertainties. We find that our model can achieve a precision $\sigma_{\rm NMAD}=0.0104$ and mean uncertainty $\langle E / (1 + z_{\rm true}) \rangle=0.0155$ for sources with ${\rm SNR}_{GI}\geq1$. For sources with ${\rm SNR}_{GI}$ higher than 3.0, 5.0 and 10.0, $\sigma_{\rm NMAD}$ can achieve 0.0047, 0.0037 and 0.0024 respectively, matching the redshift precision requirements for studies such as BAO using the CSST slitless spectroscopic surveys. Furthermore, by utilizing spatial augmentations, the network demonstrates resilience to wavelength calibration errors. This work provides a novel and robust pathway for data analysis of next-generation slitless spectroscopic galaxy surveys.

We detail the spectral analysis of a 170 ks XRISM Resolve observation of the core of Cygnus A. The high spectral resolution of Resolve have enabled us to probe the inner accretion region of Cygnus A by analyzing the 6.4 keV Fe K$\alpha$ line complex. We find that it consists of two Keplerian broadened components. (1) A broad component with a velocity dispersion of $3400^{+800}_{-600}$ km s$^{-1}$ and (2) a narrow component of $440^{+60}_{-50}$ km s$^{-1}$. For an inclination of $50^{\circ}-85^{\circ}$, constrained by VLBI, we find that the broad component arises from a distance of $\sim 0.1-0.17$ pc ($800-1400$ gravitational radii) and the narrow component from $\sim 6-10$ pc ($50,000-80,000$ gravitational radii) from the central black hole depending on the inclination angle. Our result suggests that the origin of the broad component is consistent with the broad line region and the narrow component from the torus of Cygnus A. We also find a potential emission line possibly from intermediate ionized Fe XVII with a very low dispersion ($<80$ km s$^{-1}$) that originates from either the outer edge of the torus or the narrow line region. Finally, we find that the Fe K edge is redshifted compared to the Fe K$\alpha$ line components, suggesting a line of sight bulk velocity of $470 \pm 100$ km s$^{-1}$. Such a shift may be due to an inflowing wind or relative motion between the two components originating from the near and far side of an inflowing torus, respectively.

We present near-infrared (NIR) and optical observations of the highly reddened and moderately energetic Type Ib supernova (SN) 2023dbc, {\bf covering a period from} 2 to 136 days after the explosion. By comparing its color {\bf evolution}, specifically in $r-JHK_{\mathrm{s}}$ and $i-JHK_{\mathrm{s}}$, with those of broad-lined Type Ic (Ic-BL) and Type IIb SNe, we estimate a significant extinction of $A_{V}=4.1\pm0.1$\,mag toward the SN. The extinction-corrected peak absolute magnitudes are $M_{J} = -16.8\pm0.2$\,mag, $M_{H} = -16.8\pm0.2$\,mag, and $M_{K_{\mathrm{s}}} = -17.0\pm0.2$\,mag. The SN {\bf exhibited} an $r$-band rise time of 14.9 days. The spectra {\bf display} broad features {\bf indicative of} high expansion velocities; the He~{\sc i} line velocity was measured at $16,000\,\mathrm{km\,s^{-1}}$ at $t=-4$\,d. Its spectral profile is broader than {\bf those} of typical moderately energetic Type Ib SNe, {\bf yet narrower than those of Type Ic-BL SNe, placing it in an intermediate category}. Based on the light-curve timescale and velocity, we estimate {\bf a} kinetic energy of $E_k = (4.1\pm0.7) \times 10^{51}$\,erg, {\bf an} ejecta mass of $M_{\mathrm{ej}} = 2.3\pm0.7\,M_{\odot}$, and a radioactive $^{56}\mathrm{Ni}$ mass of $(3.8\pm0.1) \times 10^{-2}\,M_{\odot}$. {\bf An} analysis using a two-component model suggests a steep density profile in the outer layer {\bf contrasted with} a dense inner core, {\bf which implies} ejecta asphericity. The low $^{56}\mathrm{Ni}$ mass is consistent with a partial fallback scenario. We conclude that SN 2023dbc originated from an aspherical explosion with partial core fallback, {\bf arising} from a progenitor ($M_{\mathrm{ini}} \simeq 15\,M_{\odot}$) that had retained {\bf its} helium envelope {\bf within} a binary system.

We present the second data release (DR2) of Spectra and Polarisation in Cutouts of Extragalactic sources from RACS (SPICE-RACS). SPICE-RACS DR2 is derived from the third low-band epoch of the Rapid ASKAP Continuum Survey (RACS-low3) and covers the entire sky from the South celestial pole up to a declination of $+49^\circ$; approximately 87.5% of the celestial sphere. We produce 'cutout' spectral cubes in Stokes $I$, $Q$, $U$ around 4 million radio sources and extract spectra towards 5 million radio components. Across our observed band of 799.5--1087.5 MHz we find an $rms$ noise of $\sim200\mu\mathrm{Jy/PSF}$, an angular resolution of $\sim15''$, and residual wide-field instrumental polarisation on the order of 0.1%. After de-duplication, our polarisation catalogue contains the detection of $2.5\times10^5$ ($3.4\times10^5$) Faraday rotation measures (RM) for components with a linearly polarised signal above $8\sigma$ ($6\sigma$). This places SPICE-RACS DR2 as the largest single RM catalogue ever produced by nearly an order of magnitude; the number of RMs in our catalogue alone is $\sim5$ times larger than every previous RM catalogue combined. Our resulting RM grid has an areal density of $6.7^{+1.8}_{-1.7}\mathrm{deg}^{-2}$, providing an effective 'resolution' of $\sim23'$, and reveals striking features across the sky. The broad-band RMs have a median uncertainty of $\sim2\ \mathrm{rad\ m}^{-2}$, and include complexity metrics and information from the time domain. The breadth and quality of the SPICE-RACS DR2 dataset will enable a new generation of RM science. Further, SPICE-RACS will provide an ideal reference for forthcoming deep polarisation surveys such as the ASKAP POSSUM survey. All of our data products are publicly available on the CSIRO Data Access Portal (DAP) and the CSIRO ASKAP Science Data Archive (CASDA).

Magnetic fields are expected to permeate the circumgalactic medium (CGM) of galaxies, yet direct constraints at high redshift remain limited by the lack of high-quality Faraday rotation measure (RM) data. Using the RMs from SPICE-RACS DR2 combined with the DESI DR1 quasar catalogue, we compile the largest sample to date of 2483 quasar sightlines with associated RMs, including 612 with intervening Mg II absorbers tracing foreground galaxies and 1871 control sightlines without Mg II absorbers. After subtracting the Galactic RM contribution and restricting the analysis to sightlines with low Milky Way HI column density and H$\alpha$ intensity, we obtain a foreground-cleaned sample of 757 quasars (191 Mg II / 566 control) spanning redshifts $0.13<z<3.45$. In this foreground-cleaned sample, Mg II sightlines exhibit a $4.5\sigma$ excess in the residual RM dispersion of $4.13 \pm 0.91~\mathrm{rad\,m^{-2}}$ relative to the control sample, at a median absorber redshift of $z\sim1.14$. This implies model-dependent CGM magnetic field strengths of $\sim0.4 - 0.8\, \mu$G over projected radii of $20 - 150$ kpc. This indicates that substantial CGM magnetisation was already established by $z\sim1$, enabling new constraints on the growth and amplification of magnetic fields in galaxy halos over cosmic time.

Polarization of the cosmic microwave background (CMB) brings out information not only on the early universe but also on the late-time large-scale structure via weak gravitational lensing. Here, we show that circular polarization is induced in principle from CMB temperature fluctuations when the optical Magnus effect is incorporated into gravitational lensing. This is a consequence of the transverse shift of a trajectory of light depending on its helicity that requires right-handed and left-handed components at the same observation point to be sourced from different points of the surface of last scattering. Whereas the resulting circular polarization is found far beyond the scope of current detection, our work establishes the optical Magnus effect on gravitational lensing as a new fundamental mechanism to produce circular polarization of CMB.

We carry out a Bayesian model comparison analysis to determine the dark matter model that best describes the Milky Way rotation curve, using four recent compilations of rotation-curve velocities with the Gaia data. We compare these data using three different baryonic models along with seven different dark matter models. Within the adopted modeling framework and Gaia-based rotation curve datasets, we find that the Einasto profile provides the preferred phenomenological fit compared to the NFW profile} across most combinations of datasets and baryonic models. We also find that the cored dark profiles are preferred over the (cuspy) NFW profile. We also test MOND using three different interpolating functions and find that, within the implementations considered, the models provide poorer fits than the preferred dark matter profiles. Finally, among the different baryonic models considered, none is decisively favored over the others.

Red novae or luminous red novae are a class of optical transients that have emerged over the past two decades. They occupy an intermediate luminosity regime between classical novae and supernovae and are characterized by cool, slowly expanding ejecta and a pronounced evolution toward red, dust-enshrouded remnants. These events are now widely interpreted as the outcome of binary coalescence involving non-compact stars, providing a rare opportunity to directly observe the dynamical phases of stellar mergers and their immediate aftermath. Observational studies of red novae provide a glimpse into the still poorly understood physics of unstable mass transfer and common-envelope evolution in binary stars, responsible for the formation of high-energy astrophysical phenomena, compact binary systems, and gravitational wave sources. In this review, we synthesize current observational knowledge of red novae, including their outburst properties, population characteristics, and long-term remnants. Observations of light curves, spectra, and circumstellar environments reveal a complex interplay between mass ejection, collisions, radiative processes, and dust formation. Archival detections of red novae progenitors show a diversity of systems, ranging from low-mass contact binaries to massive evolved stars, with a notable representation of post-main-sequence stars. We examine current efforts to predict red nova outbursts and establish robust event rates, both of which remain challenging. The growing sample of extragalactic transients suggests that the brightest red novae may be even more frequent than core-collapse SNe in the local Universe, underscoring their importance for binary evolution and stellar population studies. Finally, we outline future prospects, including the impact of large-scale time-domain surveys and the potential connection between stellar mergers and gravitational-wave sources.

We analyze the line ratio of the $^{13}$CO (2-1) to $^{12}$CO (2-1) rotational transitions observed from new ALMA observations of 100 Giant Molecular Clouds (GMCs) that span the Large Magellanic Cloud. We measure a median line ratio of $^{13}\mathrm{CO}(2{-}1)/^{12}\mathrm{CO}(2{-}1) = 0.078$ with $68\%$ of the sample falling between 0.058 and 0.107. A regression analysis confirms a nearly linear relationship across two orders of magnitude in line luminosity. Moreover, we find that the inclusion of $(L_{\text{FIR}})$ from Young Stellar Objects as a predictor variable of the line ratio significantly improves the quality of the fit, with clouds hosting IR-bright YSOs having relatively brighter $^{13}$CO emission. This analysis indicates that active star forming molecular clouds have different internal conditions than more quiescent clouds.

JWST/NIRCam observations in the North Ecliptic Pole Time Domain Field (TDF) identify 4.4 micron counterparts for 206 of 211 radio sources with S(3 GHz) \gapprox 5 micro-Jy in a 65arcmin^2 field. One of the remaining radio sources is likely to be a radio lobe of a nearby Seyfert galaxy, and the four radio sources without counterparts could be spurious. All but five counterparts are brighter than magnitude 23.5 AB at 4.4 micron. A simple position match with radius 0.3 arcsec would have identified 198 of the counterparts but only in a 4.4 micron catalog created with aggressive deblending of multiple peaks within an object's brightness distribution into distinct catalog sources. The properties of the radio-host galaxies are mostly consistent with those found in Paper 1: the median redshift is 1.14, and the radio emission, calculated taking into account the non-linear dependence of radio luminosity on star-formation rate, is consistent with a star formation origin in ~79% of the sample. For the other ~21%, the radio flux could come from star formation hidden behind dust or from an active galactic nucleus. One difference from other studies of radio-source counterparts is that 66% of the radio hosts show at least one indication of an AGN's presence. The presence of AGN and of hidden star formation could be elucidated by monitoring for source variability, and the TDF is the field most suited to such observations.

Magnetars are among the most extreme laboratories in the universe, harboring surface magnetic fields reaching $10^{15}$~G. At these supercritical scales, Maxwell's linear electrodynamics is superseded by Nonlinear Electrodynamics (NLED). While vacuum birefringence has provided initial observational evidence for these effects, its broader impact on photon propagation remains largely unexplored. In this work, we demonstrate that NLED significantly alters photon propagation in the vicinity of magnetars, deviating light from standard null-geodesics. We estimate that neglecting these corrections leads to relative errors in inferred stellar radii by means of ray-tracing techniques of approximately $10\%$. Furthermore, we find that NLED induces a systematic minimal travel-time delay of approximately $350~n$s, a value that already far exceeds the $100$~ns temporal resolution of missions like NICER. These results are critical for the interpretation of X-ray pulse profiles from current and future observatories, such as eXTP, which rely on high-precision light-bending and timing models to determine neutron-star masses and radii. Finally, our results underscore the role of magnetars as a vital window into the physics of superdense matter and supercritical fields, and we briefly highlight other astrophysical observables--such as glitches and antiglitches--that may be affected by NLED.

The Systematic KMTNet Planetary Anomaly Search series was conducted using the KMTNet data archived from $2016$ to $2019$. From this first phase of the series, we reported a total of $50$ planetary systems hidden in the data archive, which represent about $35\%$ of the total microlensing planets discovered from $2016$ to $2019$, demonstrating that this semi-machine-based search is a crucial channel for building a complete microlensing planet sample. We continue this series for $2021$ and beyond to expand the microlensing planet sample. In this work for the $2021$ KMTNet high-cadence fields (Prime fields), we find seven hidden planetary systems and three planet candidates. These new planets represent about $33\%$ of the total microlensing planets discovered within the Prime fields observed during the $2021$ bulge season. While the by-eye search is the primary channel for detecting microlensing planets (i.e., two-thirds of microlensing planet discoveries), this work clearly shows that a systematic search series is still necessary for constructing a complete microlensing planet sample. Such a sample is essential for conducting unbiased statistical studies of planet demographics in our Galaxy. Datasets for all the events used for analyses in this work are publicly available.

Even though the prompt emission of gamma-ray bursts (GRBs) is highly beamed, high-latitude emission still produces a distinct light curve break after the intrinsic emission ceases and the edge of the jet comes into view. This curvature effect offers a direct probe of the jet structure during the prompt phase. To uncover the geometric structure of the GRB jet encoded in the prompt light-curve evolution, we develop a numerical model that calculates synchrotron light curves from structured jets to interpret the observed break. We apply this model to the prompt emission of GRB 230307A, which displays a rare late-time break. Our analysis demonstrates that simple spherical outflow and top-hat jet models are inadequate to reproduce the light curve. Instead, the observations are best described by a power-law wing jet with a uniform core ($\theta_{\rm core}=0.0147$ rad) and a surrounding power-law wing. Our results demonstrate that the break in late-time prompt emission can be a powerful diagnostic of GRB jet structure.

Reflected sunlight from a solar-system body produces a flux at Earth that scales as the heliocentric distance to the negative fourth power, whereas self-luminous emission scales as the negative second power. This difference defines the $\alpha$-slope test proposed by Loeb & Turner (2012), a photometric technosignature diagnostic applicable to any solar-system body observed at multiple heliocentric distances. Of 22 (observatory $\times$ band) analysis bins for Pluto in the Minor Planet Center (MPC) archive, none recovers the reflected-sunlight flux--distance slope predicted when photometry is restricted to a single instrument and band. The archive cannot cleanly execute the $\alpha$-slope test on the brightest, most-observed trans-Neptunian object. We formalize a six-criterion eligibility pipeline (Q1--Q6) for the Loeb & Turner technosignature test and apply it to every numbered TNO. Of 8,557 candidate bins (KBO $\times$ observatory $\times$ band), 1,089 pass Q1--Q3 and 186 additionally pass Q4--Q6, splitting into 53 consistent with reflected sunlight ($\alpha = -4$), 24 with self-luminous emission ($\alpha = -2$), and 109 anomalous. The anomalous bins exhibit slopes steeper than $\alpha = -4$ or shallower than $-2$, consistent with uncorrected per-instrument calibration offsets rather than any single physical mechanism. All 24 self-luminous-like bins originate from Pan-STARRS PS1/PS2; no other observatory contributes any. This indicates a per-instrument calibration systematic. The Rubin Observatory's ten-year survey will deliver uniform single-instrument calibration on a tenfold larger sample and either resolve the test at $>10\sigma$ on hundreds of TNOs or, by reproducing the same clustering, falsify the calibration-systematic interpretation.

Certain solar flares exhibit a distinctive candle-flame or cusp-shaped feature above the bright flare arcade visible in extreme ultraviolet (EUV) and X-ray channels sensitive to high-temperature plasma. The presence of a cusp-like structure is generally regarded as a key piece of morphological evidence for magnetic reconnection to power explosive energy release in solar flares. In addition, downward-propagating plasma flows above the flare arcade have often been interpreted as outflows driven by magnetic reconnection. However, the relationship between the observed candle-flame-shaped morphology and the underlying magnetic field geometry for reconnection remains unclear. Likewise, the observed speed of the plasma downflows has been found to be too low compared to the upstream Alfvén speed predicted by reconnection theories. With the help of a recently developed three-dimensional magnetohydrodynamics (MHD) model, we examine the locations where magnetic topology changes from antiparallel to closed (Y-points) in a candle-flame-shaped flare, compare the observational emission features with synthetic EUV images generated from the model, and analyze their time evolutions. We also investigate the role of projection effects and line-of-sight integration in the measurements of plasma downflow speeds. Our analysis reveals that the Y-points do not necessarily coincide with the apparent cusp tip. Also, the apparent speeds of the supra-arcade downflows, as derived from tracks in the time-distance plots, underestimate the true Alfvén speeds in the reconnection inflow region by at least a factor of two up to an order of magnitude.

The one-parameter fuzzy dark matter (FDM) model has faced increasingly stringent constraints from both Lyman-$\alpha$ forest observations and local measurements of dwarf galaxies. A natural extension to mitigate these limits is the inclusion of FDM self-interactions. In this study, we derive constraints in the two-dimensional parameter space $(m_a, f_a)$ using the dark matter density profile inferred from a Jeans analysis of the stellar kinematics in the dwarf galaxy Leo II, which has previously been employed to constrain non-interacting FDM. We find that, for a fixed particle mass $m_a$, attractive (repulsive) self-interaction leads to a more concentrated (more diffuse) FDM density profile relative to the non-interacting case, thereby improving (worsening) agreement with the Jeans analysis results. Our results indicate that, for either attractive or repulsive SI with strength $f_a^{-1}\lesssim 10^{-14}\,\mathrm{GeV}^{-1}$, the $95\%$ confidence-level lower limits on $m_a$ lies within the range $(1-10)\times10^{-22}\,\mathrm{eV}$, although the precise bounds depend to some extent on the statistical method employed. This analysis simultaneously constrains the two parameters $(m_a, f_a)$ without relying on assumptions about cosmological or galaxy evolution histories, and thus offers a complementary probe to existing constraints.

We consider predictions for the MeV gamma-ray spectrum emitted by the $\beta$ decays of freshly synthesized isotopes from a neutron star merger at timescales of relevance for post-merger (days) and remnant (years) emission. We develop a search algorithm to identify observable spectral peaks and then determine if a specific isotope has a dominant emission line producing the spectral feature. We predict emission spectra using nucleosynthesis calculations which consider nuclear models with distinct masses, $\beta$-decays, and fission properties as well as variations on main ($A>130$) and weak ($A<130$) $r$-process astrophysical conditions. We tabulate all lines from decaying isotopes that our procedure identifies and provide the predicted range in time over which each line could be visible. We find that Rh-106 presents a unique opportunity to distinguish between main and weak $r$-process emission, as our calculated spectrum above $\sim 1$ MeV for an event dominated by the weak $r$ process is identical to the Rh-106 emission spectrum from $\sim$ 0.2 to $\sim$17 years. We further find emission from species such as Hf-181, Ta-182, Ta-184, and Re-188 offers the potential to be able to distinguish between nuclear models. We investigate whether the 2.6 MeV strong gamma-ray line from Tl-208 is predicted to be robustly observable across calculation variations on both timescale of days and years. We find Tl-208 to consistently shine through on the order of years, though it can face competition from Ga-72 and La-140 at early times ($\sim$ days). We additionally highlight numerous isotopes of interest for observation and nuclear experiment.

We present resolved HI observations of six dwarf galaxies drawn from a sample of baryon-dominated dwarf galaxy (BDDG) candidates previously identified using global HI spectra from ALFALFA and optical inclinations from SDSS, both of which suffer from systematic uncertainties in irregular dwarf galaxies. Using uGMRT interferometric observations, we obtain high-resolution HI cubes that enable more reliable determination of their geometry, circular velocity, and dynamical mass. We find that optical axial ratios systematically underestimate true disc thickness, inflating inclinations and underestimating rotation velocities in earlier work. Our HI-derived axial ratios and kinematic position angles yield larger inclination corrections and hence larger dynamical masses. Four of these galaxies, UGC 6438, UGC 7983, AGC 191707, and AGC 733302, appear dark-matter deficient. The latter three of these four exhibit high baryon enhancement efficiency factor (ratio of baryon mass accumulated by a halo to the maximum expected value for its halo mass) exceeding 50%, with AGC 191707 appearing formally super-efficient. Only UGC 9500 and AGC 220901 are consistent with being dark-matter dominated. Two of these high-efficiency dwarf galaxies lie in relatively isolated environments, showing no clear signatures of tidal disturbance or stripping, making their dark-matter deficiency difficult to reconcile with standard $\Lambda CDM$ expectations for low-mass halos. Our results underscore the importance of resolved HI kinematics in confirming genuine BDDGs and suggest that more such systems may exist. Identifying a larger sample is essential for assessing their implications for baryon-halo coupling and structure formation within the $\Lambda CDM$ paradigm.

We report optical and near-infrared photometry of a WZ Sge-type dwarf nova, TCP J23580961+5502508, during its 2022 superoutburst, obtained using the 1.5-m Kanata telescope. Our observation detected early superhumps on three consecutive nights which included the rising phase toward the peak of the outburst. The early superhumps exhibited a profile dominated by the primary maximum during the first two days, while a prominent secondary minimum appeared on the third day. We reconstructed the structure of the accretion disk from multi-color light curves using the early superhump mapping. The accretion disk has a prominent flaring structure on the leading side of the disk during the first two days. An additional flaring structure emerged on the opposite side on the third day, forming a two-armed pattern that can be interpreted within the framework of the 2:1 resonance model. The reconstructed disk structure in the first two days suggests the presence of an additional mechanism operating during the initial stage of early superhump development.

Fast blue optical transients (FBOTs) are luminous, rapidly evolving explosions whose radio emission provides a sensitive probe of shock interaction and the circumstellar material (CSM) surrounding the progenitor. However, the origin of their diverse radio light-curve morphologies, especially the very steep post-peak declines seen in several well-sampled events, remains unclear. We present a forward-shock synchrotron model in which mildly relativistic ejecta interact with a dense but radially confined CSM. The CSM is described by a broken power-law density profile, and the radio emission is modeled by including both synchrotron self-absorption and external free-free absorption. Applying this framework to multi-frequency radio observations of a representative sample of FBOTs, we show that their radio diversity can be explained by shock propagation through a finite CSM shell. The early radio evolution is regulated by absorption, while the rapid post-peak fading marks the forward shock's transition from the dense inner CSM into a more tenuous outer environment. The inferred shock velocities are trans-relativistic, $v_{\rm sh}\sim0.1$--$0.5c$. The radio-emitting CSM requires high mass-loading rates, $\dot{M}\sim10^{-4}$--$10^{-3}\,M_{\odot}\,{\rm yr}^{-1}$, but modest total CSM masses, $M_{\rm CSM}\sim10^{-4}$--$10^{-2}\,M_{\odot}$. These properties point to brief episodes of enhanced mass loss in the final years to decades before explosion, rather than long-lived steady winds. Our results provide a dynamically consistent interpretation of FBOT radio emission and establish radio light curves as a diagnostic of the immediate pre-explosion mass-loss history of FBOT progenitors.

The shape of the power spectrum of the black hole low-mass X-ray binary evolves systematically over different spectral states during an outburst. Therefore, the power colours (ratio of the variability amplitude at different frequency ranges) and the hue parameter, quantifying the power spectral shape, can be utilised to identify the spectral states of the system. We present the comprehensive power colour analysis and subsequent identification of spectral states using the entire archival data (2016-2024) from AstroSat. We detected 29 QPOs (quasi-periodic oscillations), along with several associated harmonics and shoulders, and investigated their properties as a function of hue. We examined the evolution of the QPO RMS variability and time lag, along with hue and QPO frequency. We report a sign change in the average QPO time lag around the QPO frequency of ~2 Hz for high inclination sources, during the hard-to-hard intermediate state transition. At lower frequency, the hard lags showed an increasing trend reaching up to ~36 ms, but the soft lags above 2 Hz remained confined within ~10 ms, suggesting an evolution to a compact corona. Conversely, for low inclination sources, no such transition was found. Furthermore, for high inclination sources, the harmonic lag remains unaffected during state transition, in contrast to the QPO lag behaviour. Our results are consistent with a transition from an elongated jet-like corona to a compact corona and reveal vital clues about the dynamical evolution of the corona and disk.

Solar flares are intense eruptive events caused by the rapid release of magnetic energy, often impacting Earth's space environment through electromagnetic radiation and high-energy particles. Accurate flare prediction is critical for space weather forecasting. However, many existing deep learning approaches often rely on single-modal inputs or shallow feature fusion, limiting their ability to capture complementary information. In this study, we propose a dual-branch multimodal fusion deep learning model for predicting 24-hour solar flares. The model integrates magnetograms and magnetic parameters through cross-attention mechanisms, followed by cross-scale interactions at the feature level to enhance multi-scale representation. It is designed to perform both binary prediction of $\geqslant$ C-class flares and multi-class classification of C, M, and X-class flares. To ensure rigorous evaluation, we employ a stratified group five-fold cross-validation scheme to preserve class representativeness and adopt a splitting-before-sampling strategy based on NOAA active region numbers to prevent data leakage. Experimental results show that the model achieves a TSS of 0.661 and an HSS of 0.658 for binary $\geqslant$ C-class prediction, while notably attaining a TSS of 0.780 and an HSS of 0.775 for X-class flares in the multi-class task. Compared with existing approaches, the model demonstrates superior performance in predicting intense X-class flares, effectively suppresses the false alarm rate, and exhibits strong generalization capability.

Hub-filament systems (HFSs) play an important role in the formation of massive stars and star clusters. Although the velocity structures along dense filaments have been studied, the gas kinematics in the low density inter-filament regions has not been investigated. We use $^{13}$CO ($J$ = 1--0) and C$^{18}$O ($J$ = 1--0) observations obtained with the Nobeyama 45 m telescope to study the gas dynamics towards the Monoceros R2 (Mon R2) HFS. From the $^{13}$CO and C$^{18}$O data, tracing low- and high-density gas, respectively, we identify velocity coherent structures and divide them into filaments (Fs) and inter-filamentary regions (IFs). We estimate velocity gradients ($\Delta v$) and mass accretion rates ($\dot{M}$) along ($\parallel$) and across ($\perp$) the Fs and IFs. The mean ratio of $\dot{M}_\parallel$ to $\dot{M}_\perp$ in Fs is 6.8, while that in IFs is 1.5. These results show that the overall gas within both Fs and IFs flows directly into the hub and the gas flows faster along the Fs than the IFs. In addition, we found that at least 30\% of the gas mass in the IFs may flow towards the Fs replenishing the latter with new matter. Our study reveals the importance of considering the total gas mass reservoir, both low- and high-density, infalling into the hub and promoting the formation of massive stars, which are preferentially located in the hub of Mon R2.

The coalescence of binary neutron stars in the GW170817 event led to the generation of gravitational waves, accompanied by the electromagnetic counterpart known as a kilonova (KN). Since then, it has been a prime topic of interest, as it has provided much insight into multi-messenger astronomy. Apart from existing methods for parameter estimation, we propose an alternative technique for it, utilizing the strength and flexibility of a conditional variational autoencoder. Publicly available light curves are used as training data, conditioning on the corresponding physical parameters for a chosen model; after training, we carry out rapid parameter inferences. As this approach approximates the likelihood through variational inference, it yields results more efficiently. Through this innovative approach, we demonstrated that the total time, from training to parameter inference, is under $\approx3$h. We showed that for a given KN light curve, we can rapidly perform parameter inference based on the required model.

We investigate the $\gamma$-ray and X-ray population properties of changing-look blazars (CLBs) using sources from the Fourth \textit{Fermi} LAT Source Catalog Data Release 4 (4FGL-DR4) together with X-ray information from the Living \textit{Swift} XRT Point Source (LSXPS) catalog. The CLB sample is compared with large populations of confirmed BL Lac objects (BLLs) and flat-spectrum radio quasars (FSRQs) using spectral, variability, and broadband properties. In the $\gamma$-ray parameter space, CLBs mainly occupy intermediate and overlap regions between the BLL and FSRQ populations. However, the centroid locations in different parameter planes, along with the PCA and UMAP projections, show that the CLB population lies closer to the FSRQ region. The X-ray analysis also shows a similar behavior, where the overall distribution of CLBs in the X-ray parameter space is found to be nearer to FSRQs than to BLLs. In addition, the X-ray/$\gamma$-ray coupling relations and random-forest classification probabilities are consistent with this trend. Overall, the results suggest that CLBs form a transitional population between the two main blazar subclasses while retaining characteristics closer to the FSRQ population.

Using high-resolution photospheric and chromospheric observations taken by the Goode Solar Telescope (GST), we studied two groups of bidirectional plasma jets occurring in active region NOAA 13110. Supplementary observations are also provided by Solar Dynamics Observatory's (SDO) Atmospheric Imaging Assembly (AIA) and Helioseismic and Magnetic Imager (HMI). From the photospheric observations and magnetograms, the two successive bidirectional plasma jets were initially located in the vicinity of the polarity inversion lines and at the intersection of the umbra and penumbra of the sunspot, followed by magnetic flux emergence and cancellation. As the cool filamentary threads are continuously emerging from the lower chromosphere and interact with overlying horizontal magnetic loops, it leads to the bidirectional jets, erupting outflow plasmoids, and heating coronal magnetic loops. We find that the bidirectional jets extended from the central excitation location in opposite directions, at the speed of about dozens of km s$^{-1}$. For the first jet, the initial brightening first appears in 304 angstroms, about 30 s earlier than the H$\alpha$ observations, indicating that magnetic reconnection takes place in the transition region. While the initial reconnection for the second jet occurs simultaneously in H$\alpha$ and 304 angstroms, showing the recurrent eruptions. These observations confirm that the bidirectional plasma jets can be generated by magnetic reconnection between the rising filamentary threads or material and the overlying horizontal magnetic loops. Our results provide new insights into the generation of the bidirectional plasma jets and reconnection-based coronal heating.

The Chinese Space Station Survey Telescope (CSST) aims to map the universe across an unprecedented dynamic range of stellar densities, spanning from extragalactic voids to the crowded Galactic center (e.g. a few stars and galaxies in the voids and $>10^5$ stars per detector in Galactic center). However, processing such heterogeneous data with a general source extraction pipeline introduces significant systematic uncertainties, standard algorithms exhibit poor accuracy in crowded fields and suffer from increased astrometric uncertainty in void regions. To mitigate these systematics, we propose a hierarchical, two-stage deep learning model for adaptive data reduction. The first stage ('classification') employs a ResNet-34 model to classify images into six discrete density categories, achieving $98.83\%$ in global accuracy. This classification acts as a critical decision gate, ensuring high calibration accuracy in the crowded fields. In the second stage ('regression'), a ResNet-50 regression model predicts the bright stars ($<23.5$ mag) in the field, which is essential for astrometric calibration, achieving a mean absolute error (MAE) of 0.0824 dex. By decoupling density characterization from source extraction, our model ensures that photometric and astrometric algorithms are optimally matched to the stellar density environment, thereby enhancing the fidelity and homogeneity of CSST as well as future large sky survey data products.

We present a detailed spectro-temporal study of the Be X-ray binary pulsar $1A$ $1118-61$ during its brightest recorded outburst in 2026, using \textit{Swift} and \textit{NuSTAR} observations, and compare its properties with the 2009 outburst. Coherent pulsations at $\sim400$ s are detected throughout the outburst, with pulse profiles evolving strongly with energy and luminosity, indicating changes in emission geometry. A transient quasi-periodic oscillation (QPO) at $\sim$0.11 Hz is observed during the rising phase. The luminosity dependence of the QPO frequency during the current and previous outbursts suggests an origin associated with instabilities near the magnetospheric radius. The broadband spectra are well described by thermal Comptonization and show clear spectral hardening at higher luminosities. A cyclotron line is detected in the two \textit{NuSTAR} observations, with its energy remaining nearly constant despite a factor of $\sim25$ change in luminosity. Long-term monitoring reveals that the 2026 outburst is systematically harder and brighter, suggesting significant difference in the accretion structure and emission regions between the two outbursts.

We perform three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations of a near-maximally spinning black hole (spin parameter, a = 0.998) with varying initial magnetic field geometries, systematically exploring the parameter space connecting magnetically arrested disk (MAD), intermediate (INT), and standard and normal evolution (SANE) accretion states. The magnetic flux threading the black hole horizon emerges as the fundamental state variable controlling jet efficiency, flow magnetization, and radiative output across all three states. We introduce complementary diagnostics-broadband spectral energy distributions spanning radio through hard X-ray frequencies and time-resolved X-ray light curves-that together connect simulation dynamics directly to multiwavelength observables. The radiative output follows a clear MAD > INT > SANE hierarchy in time-averaged luminosity, mean X-ray emission, as well as variability. Furthermore, MAD exhibits the highest fractional variability through quasi-periodic magnetic flux eruption events, and INT and SANE show moderate variability driven by episodic reconnection and stochastic MRI turbulence, respectively. Scaling to GRS 1915+105, Cyg X-1, and HLX-1, we demonstrate that all twelve temporal classes of GRS 1915+105 map naturally onto our three magnetic states, Cyg X-1's persistent hard state is reproduced by a sustained INT configuration, and HLX-1's extreme luminosities arise through efficient Blandford-Znajek extraction in MAD states scaled to higher black hole mass.

Understanding galaxy morphology evolution across cosmic time requires models that can generate realistic galaxy populations conditioned on redshift. In this work, we study efficient redshift-conditioned generative modeling for astrophysical image synthesis using diffusion models and pixel-MeanFlow. We first review the connections between score-based diffusion models, Flow Matching, one-step generative models, and modern diffusion samplers. We then evaluate DDPM, DDIM, DEIS-AB2, DPM++2M, and one-step pixel-MeanFlow on the GalaxiesML-64 dataset using morphology-based metrics, including ellipticity, semi-major axis, Sérsic index, and isophotal area. Our results show a clear accuracy-efficiency trade-off: standard DDPM sampling achieves the best distributional fidelity but requires high computational cost, while second-order samplers substantially improve efficiency over DDIM. Pixel-MeanFlow enables single-step generation and achieves competitive performance on several morphology statistics, though it remains weaker than many-step DDPM for fine-grained structure. Our results demonstrate that one-step generative models can recover key galaxy morphology statistics at orders-of-magnitude lower computational cost, opening a path toward efficient conditional simulators for large cosmological surveys and simulation-based scientific inference.

X-ray surveys of active galactic nuclei (AGNs) provide direct constraints on the properties of individual AGNs, such as their emission, obscuration, and accretion rate. Previous AGN population synthesis models have not addressed such properties self-consistently. Here, we use a simulation-based inference (SBI) approach to constrain the geometrical and physical properties of the AGN population. We perform numerical simulations with our ray-tracing code, RefleX, which allows the self-consistent modelling of the X-ray emission of AGNs with flexible circumnuclear and source geometries. We create our synthetic population by sampling the intrinsic active black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) of local AGNs, and we construct a geometry based on the radiation-regulated model, along with Eddington-ratio-dependent emission spectra. Using the RefleX-simulated emission of the AGN population, we aim to simultaneously reproduce the cosmic X-ray background (CXB), differential AGN number counts, and several observed absorption properties of local AGNs, such as the fraction of $N_\mathrm{H}$ in bins of log($N_\mathrm{H}$), the Compton-thick fraction as a function of limiting flux, and the number of obscured and unobscured AGNs as a function of Eddington ratio. With this approach, we test the consistency of the radiation-regulated model with a very comprehensive set of X-ray observables, while constraining the size and density of the dusty torus and the evolution of the local AGN population. We derive an intrinsic Compton-thick fraction of $40\pm3$%, and find that a simple evolutionary prescription controlling the active fraction of supermassive black holes is sufficient for our synthetic population to reproduce the CXB.

Neutrinos and their weak interactions play a vital role in the physics of core-collapse supernovae and binary neutron star mergers. Their description within astrophysical simulations, including the weak rates, is of pivotal importance not only for the prediction of accurate neutrino fluxes and spectra, including the associated conditions relevant to nucleosynthesis, neutrinos are also responsible for heating and cooling of the stellar plasma as well as the transport of lepton number and entropy. In the present article, we develop an essential improvement of the description of the underlying nuclear medium, necessary for the calculations of charged-current weak rates, with the inclusion of explicitly momentum-dependent nuclear interactions. To this end, we introduce the relativistic Hartree-Fock (RHF) approach and the associated momentum-dependent nucleon self-energies. We discuss the resulting neutrino and antineutrino opacities and find large discrepancies comparing the weak rates at the RHF level with those of commonly used relativistic mean-field (RMF) models; in particular, we observe a substantial shift of previously reported large medium-dependent modifications associated with the RMF approach.

We report the detection of diffuse gamma-ray emission toward the young massive star cluster Berkeley 87 using Fermi data. The emission has an angular extension of 0.36 degree and a photon index of 2.68. The hadronic scenario is favored given the dense gas and the cluster's strong stellar winds.

We present detailed analyses of three anomalous microlensing events--KMT-2021-BLG-0209, KMT-2021-BLG-0901, and OGLE-2025-BLG-0356--identified from a systematic re-examination of KMTNet light curves for which previous modeling attempts failed or left persistent residuals. Although all three events show caustic-related features consistent with binary-lens microlensing, we find that their full light-curve structures can be described by four-body configurations that required four-body configurations involving a binary lens and a binary source. In KMT-2021-BLG-0209, weak caustic-exit residuals arise from a faint companion source undergoing an additional caustic interaction. In KMT-2021-BLG-0901, a late-time re-brightening is produced when the secondary source encounters the resonant caustic long after the primary. For OGLE-2025-BLG-0356, we test the degeneracy between 3L1S and 2L2S interpretations of a short isolated anomaly and find that the 2L2S model provides a significantly better fit. Source colors and magnitudes indicate binary sources composed of (G8V, M3V), (G8V, K2V), and (G6V, G8V) stars for the three events, respectively. Bayesian inference suggests that the lenses are predominantly low-mass binaries, including one system (KMT-2021-BLG-0901) with a companion consistent with a brown dwarf. These events add to the growing sample of well-characterized 2L2S systems and underscore the importance of systematically testing complex models, particularly in anticipation of the high-precision microlensing data expected from the Roman Space Telescope survey.

The human presence in outer space is undergoing a transition from one in which nation states are the dominant actors to an emerging status quo in which states merely supervise the activities of private entities. Such largely commercial ventures include extracting natural materials from the Moon, a celestial body of great cultural and spiritual reverence for some Indigenous societies. However, the existing international legal framework governing activities in space focuses on its "exploration and use", centered in a Western worldview that attaches to a past history of colonialism. While that framework, articulated in the Outer Space Treaty (OST), claims to guarantee that outer space will remain "the province of all [hu]mankind," only entities with significant political power have to date decided the limits of the acceptable uses of space. This paper examines the historical record for clues about how states, private actors and Indigenous societies might interact in the future on matters of outer space governance to achieve more just ends. It analyzes a key case study: the dispute between the NASA and the Diné people of the American Southwest over the launches of human cremated remains to the Moon in 1998 and 2024, acts the Diné president called "deeply disturbing and unacceptable to our people and many other tribal nations." In a future in which space becomes increasingly commercialized and entities like NASA transform into mere contract-administering agencies, it is unclear how an impending, exploitative human presence on the Moon can simultaneously honor Indigenous rights and perspectives on lunar issues. The presentation concludes that best practices for future engagement with the Moon must transcend the mere "due regard" and "international consultations" required by the OST in favor of arrangements where participants co-create a human future in outer space.

We present high-resolution magnetic field maps of the M17 SW molecular cloud using JCMT 850 $\mu$m dust polarization at a scale of 14$''$. The magnetic field exhibits a distinct arc-like structure that encircles three dense clumps (C1, C2, and C3). By combining polarization data with ammonia line observations, the plane-of-sky magnetic field strength, measured using the Skalidis-Tassis method to minimize angle dispersion errors, ranges from 0.1 to 2.4 mG (mean: 0.54 mG). Energy budget analysis reveals a hierarchy dominated by gravity ($e_G \approx 10^{-7.8}$ erg cm$^{-3}$), which exceeds both magnetic ($e_B \approx 10^{-8.3}$ erg cm$^{-3}$) and turbulent ($e_k \approx 10^{-8.7}$ erg cm$^{-3}$) energies. Since all three energy densities lie within one order of magnitude, gravitational dominance acts primarily as the global driver, while the system remains in a state of near-equipartition. Structurally, the northeastern boundary shows magnetic field lines perpendicular to the shock front, consistent with compression from the adjacent HII region. Within the cloud, magnetic field lines generally align with gravity to assist collapse, but turn perpendicular to gravity within curved accretion bridges. This configuration provides support against radial collapse while guiding gas flow. Kinematic evidence suggests that these channels transport material from Clump C3 onto the massive Clump C2. Star formation in M17 SW is globally driven by gravity but locally regulated by the magnetic field structure.

IGR J20084+3221 is a previously unclassified Galactic source first detected by INTEGRAL. Chandra observations led to possible classifications of either a magnetic Cataclysmic Variable (mCV) or high mass X-ray binary (HMXB) based on the hardness of its spectrum. Here, we report follow-up observations taken by XMM-Newton, NuSTAR, and the Hale Telescope at Palomar Observatory. Based on these observations, we conclude that IGR J20084+3221 is most likely an Intermediate Polar (IP) type mCV. Timing analysis of the X-ray data found a significant peak period of $P=635.0\pm0.4$ s, which we interpret to be the spin period of the white dwarf (WD). The X-ray spectrum is well fit to an absorbed Bremsstrahlung model with components accounting for partial covering, reflection, and a fluorescent Fe-line, all typical for an IP. The optical spectrum shows clear emission lines, consistent with emission dominated by an accretion disk. We find counterparts to the source across the optical and infrared (IR) bands, and, despite uncertainties in the distance and extinction, we estimate that the source is too faint in the IR to be an HMXB. Given the evidence pointing towards an IP classification, we fit the X-ray spectrum to a post-shock region model where we find a WD mass of $M=1.09^{+0.12}_{-0.11}\mathrm{M}_{\odot}$, larger than the average mass for a WD in an mCV.

Quintessence with a potential motivated by supergravity (SUGRA) exhibits several intriguing features. Depending on its initial conditions, it can behave either as dynamical dark energy or effectively as a cosmological constant. Moreover, when quintessence is coupled to dark matter, the effective dark-energy equation of state can cross the phantom divide. In this paper, we test both coupled and uncoupled SUGRA quintessence models using DESI BAO, DES-Dovekie SNIa, and Planck CMB data. We find that current observations strongly favor a coupling between dark energy and dark matter, with the coupling parameter deviating from zero at more than $4\sigma$. The data also favor the branch of coupled SUGRA quintessence in which the energy transfer between the two dark sectors changes sign, leading to a crossing of the phantom divide by the effective dark-energy equation of state. Interestingly, this coupled SUGRA branch is statistically indistinguishable from dark energy described by the CPL parametrization, with only a very small difference in $\chi^2_\mathrm{min}$. Our results suggest that coupled quintessence with a SUGRA potential provides a field-theoretic realization of the evolving dark energy behavior favored by the latest observations.

During the performance verification phase, XRISM observed the M82 galaxy for a net exposure of 207.7 ks, with the ultraluminous X-ray pulsar (ULXP) X-2 included in the field of view. A pulsation search identified a candidate signal with a period close to the previously known value, 1.38727 s, at a significance of $3.15\sigma$ based on Monte Carlo simulations. Using this candidate period, phase-resolved spectral analysis with the high spectral resolution of Resolve was performed. The spectra suggest that, if the candidate pulsation is real, the Fe K$\alpha$ emission line in the pulse peak phase has a larger width ($36^{+60}_{-13}$ eV) than that in the remaining phase at a significance exceeding $3\sigma$. This suggests that at least a fraction of the Fe K$\alpha$ emission is associated with the ULXP system. The observed width corresponds to a velocity dispersion of $(1.7^{+2.8}_{-0.6})\times10^3$ km s$^{-1}$, which is too large to be explained by motions in the companion star atmosphere. The rise time of the pulsation constrains the line-emitting region to be smaller than $6.3\times10^4$ km, suggesting an origin in the accretion flow. This work demonstrates the capability of XRISM Resolve for pulsation-resolved high-resolution spectroscopy of ULX pulsars.

The nature of the sub-Neptune K2-18b is debated between Hycean and mini-Neptune interpretations. We test whether self-consistent Hycean atmospheres are compatible with current JWST transmission spectra by combining one-dimensional photochemical modelling, radiative--convective equilibrium calculations, and forward modelling of transmission spectra. We assume H2-CH4-H2O atmospheres over a liquid ocean, compute altitude-dependent abundances with a 1D photochemical model, and couple them to P-T profiles that avoid runaway greenhouse states. Using the CH4-dominated 2.8-4.0 $\mu$m band, we constrain wavelength-independent offsets between NIRISS SOSS and NIRSpec G395H for multiple reductions, and then scan grids of CO and CO2 scaling factors, weighted by the CH4-band offset posteriors, to evaluate oxidised-carbon abundances consistent with the 4-5 $\mu$m region. Radiative--convective calculations further map pressures and albedos that yield non-runaway climates. Over a wide range of temperatures and pressures, liquid oceans can exist, and Hycean models with a 1 bar H2 envelope, percent-level CH4 and CO, and CO2 buffered at $\sim 10^{-3}$-$10^{-2}$ reproduce the NIRISS and NIRSpec spectra from 0.8 to 5.2 $\mu$m without invoking DMS or other additional species. Our photochemical simulations show that H2-CH4-H2O networks generically drive CO to mixing ratios of order 1-2 %. Mass-balance arguments imply that a $\sim$1 bar H2 envelope with percent-level CH4 requires interior replenishment on gigayear timescales, and the resulting vertical gradients naturally generate flat, CH4-dominated plateaux in transmission. While mini-Neptune scenarios remain viable, our results show that Hycean configurations are likewise consistent with the data, and current CO and CO2 constraints alone are not yet sufficient to rule out Hycean interpretations of K2-18b.

We investigate the effects of source confusion expected in observations with GREX-PLUS, a JAXA L-class space infrared telescope mission candidate with a wide-field infrared camera covering 2-8 um with a field of view of 0.50 deg$^2$. For the deep imaging band near 4 um, we calculate the GREX-PLUS PSF and ghost based on the latest optical design, and consider two representative imaging performance cases with PSF FWHM values of 0.9 and 1.2 arcsec. We construct simulated GREX-PLUS images at different depths by convolving JWST NIRCam imaging data from JADES, GLASS, CEERS, and COSMOS-Web with the PSF+ghost kernel. Comparing the limiting magnitudes estimated from random aperture photometry using the same aperture sizes, we find that the simulated GREX-PLUS images are shallower than the original JWST images, with larger deviations for deeper original JWST images. This likely reflects unresolved faint sources and extended PSF+ghost wings from bright sources, which elevate background fluctuations in blank regions. Nevertheless, the limiting magnitudes continue to improve with increasing integration time down to ~27 mag, without a clear plateau at depths comparable to the planned GREX-PLUS deep survey, although the improvement becomes progressively less efficient toward longer integrations. Based on Monte Carlo simulations, we estimate detection completeness and correct the number counts for magnitude bias and incompleteness, finding that confusion-induced blending can reduce the completeness even at magnitudes well above the nominal 5-sigma depth. The completeness-corrected number counts agree well with the JWST-based number counts down to around the detection limit. Overall, our results suggest that statistical studies of faint galaxies remain feasible for GREX-PLUS; however, survey planning should account for less efficient depth improvement toward longer integrations due to source confusion.

The rapid development of mega-constellations in low Earth orbit (LEO) severely impacts ground-based optical astronomical observations. By combining WorldWide Telescope (WWT) simulations with 2019 and 2023 observational data from the Xinglong Observatory 80-cm telescope and 2023 data from the Mini-SiTian (MST), we find that satellite visibility increases with deployment, particularly during the summer. For the 80-cm telescope, the fraction of images containing satellite trails increased from an average of 0.34% in 2019 to 0.7% in 2023; meanwhile, for the MST in 2023, the fraction rose from 5% in January to 12% by December, peaking at 19% in the summer. Through stratified analysis of solar elevation and local time, we find that observations during twilight and summer are particularly susceptible to satellite trail interference. Photometric analysis reveals that the interference intensity increases for fainter sources and those closer to the trails. Furthermore, a comparative analysis across different seeing conditions shows that the deviation of median standardized residuals ({\sigma}) is significantly greater under poor seeing than under good seeing conditions.

Recent progress in gravitational wave observations has positioned Pulsar Timing Arrays as a key tool for detecting the stochastic gravitational wave background in the nanohertz band. It is widely believed that this background is primarily attributed to the cosmic ensemble of inspiraling supermassive black hole binaries. While traditional analyses have predominantly focused on the spectral amplitude and frequency dependence of the gravitational wave background, higher-order statistics such as variance, skewness, and kurtosis could potentially be useful for extracting further physical information. However, these statistical moments are known to diverge when the redshift integration is extended down to z=0. In this study, we propose a strategy to resolve this issue by introducing a physically motivated lower integration limit, z_min, defined by the sensitivity for detecting individual sources. Since higher-order statistics are primarily determined by local sources, we may adopt the lowest-order approximation with respect to redshift in their computations. Under this approximation, we demonstrate that all higher-order statistics beyond the expectation value depend on the mass function only through a weighted average of the chirp mass, <\mathcal{M}^{10/3}>, irrespective of the redshift evolution model. We show that the ratio of the variance to the expectation value provides information on <\mathcal{M}^{10/3}>/<\mathcal{M}^{5/3}> independently of the total number of mergers. We also find a consistency relation between the kurtosis and the squared skewness, paving the way for testing the binary-origin hypothesis of the gravitational wave background. Our findings demonstrate that higher-order statistics provide a new window for interpreting the gravitational wave background, offering a methodology to break existing degeneracies and refine our understanding of the mass function.

Low-luminosity active galactic nuclei (LLAGN; $L_{\rm bol} < 10^{42}$~erg~s$^{-1}$) may comprise a significant fraction of the local AGN population, yet their weak emission makes them difficult to detect. In this paper, we analyse 38 LLAGN identified from a 15~GHz sub-arcsecond Very Large Array survey and assess the effectiveness of X-ray, optical, and infrared wavelengths in identifying LLAGN. We found that optical emission-line diagnostics recovered $84.2^{+15.8}_{-22.9}$\% (32/38) of the sample, X-rays detected $63.2^{+25.7}_{-19.6}$\% (24/38), and infrared methods only identified $13.2^{+14.5}_{-8.0}$\% (5/38), reflecting limited X-ray sensitivity, weak or absent optical lines, and strong host galaxy contamination in the infrared. Compared to \textit{Swift}--BAT AGN, our LLAGN are $\sim$ 4.1~dex fainter in bolometric luminosity (log $L_{\mathrm{bol}} \approx$ 39.3 - 41.9 erg s$^{-1}$), host smaller black holes ($\sim$0.7~dex lower), and accrete at much lower rates (log $\lambda_{\text{Edd}} \approx$ -6.5 to -1.3, i.e., $\sim$ 4.2~dex lower). Host galaxies span a broad range of morphologies, from disk- to bulge-dominated, with a subset exhibiting prominent bulges, potentially representing systems where nuclear activity has faded while the bulge remains dominant. LLAGN also reside in galaxies with lower stellar masses ($\sim$0.3~dex) and suppressed star formation rates ($\sim$0.5~dex) relative to \textit{Swift}--BAT AGN. Overall, LLAGN in our sample systematically host smaller and weakly accreting black holes, residing in galaxies with diverse morphologies, but lower stellar masses and reduced star formation activity, demonstrating the connection between low-level black hole accretion and host galaxy properties in the local Universe.

We report the first detection of a low-frequency quasiperiodic oscillation (QPO) extending above 250 keV in the black hole X-ray binary candidate Swift J1727.8-1613 using Insight-HXMT observations during its 2023 outburst. Swift J1727.8-1613 is one of the brightest X-ray transients discovered and presents a valuable opportunity for studying high-energy properties of QPOs. Owing to the large effective area of Insight-HXMT in hard X-ray, our observations indicate a remarkably strong QPO signal in the power spectrum above 100 keV. We utilize advanced Hilbert-Huang transform techniques to analyze phase-folded light curves across a wide range of energy bands, observing significant QPOs from 100 to 300 keV in the NaI and CsI detectors, respectively. The detection of QPO profiles above 250 keV can achieve significance levels of ${\sim} 8.9{\sigma}$ for the NaI detector and ${\sim} 5.7{\sigma}$ for the CsI detector. Our results indicate a decrease in QPO fractional rms above 100 keV and an increased soft phase lag with energy, suggesting a geometric origin for the QPOs, likely linked to the precession of a small-scale jet.

(Abridged) We aim to explore the parameters that influence the mid-infrared emission of C$_2$H$_2$ and H$_2$O, and if the spread observed in $F\rm{_{C_2H_2}}$/$F\rm{_{H_2O}}$ is tracing a variation of the C/O ratio. Our work is based on the DALI 2D thermochemical model to predict spectra readily comparable to JWST/MIRI observations. To robustly model organics in inner disks, several improvements have been made: (1) carbon chemistry adapted for warm environments, (2) updated UV shielding treatment, and (3) mutual line overlap in the raytracing. We are able to reproduce the observed C$_2$H$_2$ fluxes of T Tauri disks with a solar C/O ratio. Acetylene abundance is primarily set by a balance between formation initiated by CO dissociation by X-rays and destruction of carbon chains by atomic oxygen, the latter being generated by X-ray-induced destruction of H$_2$O and CO. The water UV shielding and hot temperatures of the inner disk also favor acetylene formation, as they prevent the destruction of carbon chains and allow overcoming activation barriers of reactions with H$_2$. C$_2$H$_2$ and H$_2$O emissions are not only sensitive to the C/O ratio but also to the total O/H elemental abundance, supporting recent claims. In particular, we find that enhanced O/H reduces acetylene emission due to an excess of atomic oxygen. $F_{\rm{C_2H_2}}$/$F_{\rm{H_2O}}$ is thus a promising tracer of the elemental composition of inner disks. Still, the dust size distribution also plays a key role in this line flux ratio. We find that increasing the abundance of small grains relative to large grains favors C$_2$H$_2$ flux over H$_2$O flux. Grain depletion does not affect the line flux ratio as previously suggested by observational works. A preliminary comparison with published JWST observations indicates a gas-phase C/O ratio below unity and suggests that enhanced O/H ratios may be common in T Tauri disks.

LHAASO J0248$+$6021, a possible $\gamma$-ray pulsar halo associated with PSR J0248$+$6021 (J0248), exhibits a highly curved spectrum as revealed by LHAASO and Fermi-LAT measurements. We propose a direct interpretation of this large curvature: the energy-loss bump in the parent electron spectrum has not yet significantly departed from the high-energy cutoff. This requires either that the ambient magnetic field strength $B$ around J0248 be lower than the typical value in the interstellar medium, or that the electron injection age be significantly shorter than the pulsar characteristic age. For the much older Geminga pulsar, the expected energy-loss bump in its $\gamma$-ray halo spectrum has shifted below $100\ \text{GeV}$, in excellent agreement with Fermi-LAT measurements. Thus, the broadband spectra of young and old pulsar halos find a unified interpretation in the picture of a time-dependent energy-loss bump. Meanwhile, the spectral measurements of LHAASO J0248$+$6021 only constrain the combination of $B$ and electron injection age. The uncertainty in $B$ leads to an order-of-magnitude variation in the fitted diffusion coefficient. Future X-ray observations are expected to break the degeneracies.

The interaction between the fast-moving ejecta and the interstellar medium can produce long-lasting radio signals after binary neutron star mergers. Searching for such radio signals is a way to test the central engine of kilonovae and short gamma-ray bursts. With a magnetar as the central engine, the spin-down energy powers the main thermal and late-time radio emissions of the kilonova. However, both the thermal and radio emissions are strongly affected by the ejecta distribution, e.g., the two-component ``blue" and ``red" emissions of AT 2017gfo corresponding to the GW 170817 event. In this study, we investigate the distribution of the merger ejecta, analyzing several possible anisotropic distributions and demonstrating their impacts on the emission properties, particularly the late-time radio light curves. Under a bipolar and equatorial ejecta configuration, corresponding to the wind and dynamical components of the merger ejecta, the late-time radio light curves reveal distinct two-peak features, which are consistent with the main thermal light curves. The anisotropic distribution of the ejecta intrinsically connects the main thermal and late-time radio emissions, forming a self-consistent evolutionary picture. A combined analysis of the main thermal and late-time radio emissions provides a way to constrain the geometry of the merger ejecta and to probe the properties of the central engine. Furthermore, using the fitting parameters from the main thermal emission of AT 2017gfo, we calculate the corresponding potential late-time radio light curves. The results show that, under typical parameters, the non-detection of radio signals in observations is consistent with the theoretical expectation.

(abridged) Elemental abundances of FGK stars can be derived routinely from high-resolution optical spectra, but this remains considerably more difficult for cooler stars. Machine-learning methods offer a practical route to infer otherwise inaccessible abundances from more widely available stellar data. We use a large APOGEE DR17 sample of red giant stars as the main training set and an independent HARPS sample of nearby FGK dwarfs for external validation. We benchmark several machine-learning regressors, optimise the strongest models, and analyse feature importance using gain-based metrics, permutation importance, single-feature models, and SHAP values. We also explored the prediction of C and O from Mg, Si, and [Fe/H], and derived simple empirical relations between selected abundance ratios (Fe/Si, Mg/Si, C/O, and Fe/O) and metallicity. Kinematic information alone recovers only a limited fraction of the variance in stellar metallicity, with a clear performance ceiling at RMSE $\sim$0.20 dex. The most informative predictor is the maximum vertical orbital excursion, $Z_{\max}$, followed by radial orbital parameters. When [Fe/H] is combined with kinematic information, the abundances of C, O, Mg, and Si are predicted significantly more accurately than with the baseline approximation $\mathrm{[X/H]}=\mathrm{[Fe/H]}$. In contrast, when predicting C and O from Mg, Si, and [Fe/H], most of the predictive power is already contained in the elemental abundances themselves, with Mg being the dominant contributor, and the addition of kinematic information provides little improvement. The trained models reproduce the main abundance trends associated with Galactic chemical evolution. We find that the slopes of the relations between Fe/Si, Mg/Si, C/O, and Fe/O and metallicity differ slightly between the HARPS and APOGEE samples, with fractional differences generally below 17\%.

M dwarfs host strong magnetic fields that can be measured using several complementary techniques. However, the impact of key methodological choices on Zeeman broadening diagnostics has not been systematically quantified. We assess the reliability of different approaches for inferring magnetic fields in M dwarfs using synthetic Stokes $I$ spectra of Ti I lines generated from three-dimensional magnetohydrodynamic simulations combined with MARCS model atmospheres. Synthetic observations were produced for different surface magnetic field strengths, projected rotational velocities, and inclination angles. Zeeman broadening and intensification were analysed using polarised radiative transfer calculations coupled with Markov chain Monte Carlo inference. We evaluated several statistical criteria (BIC, AIC, and WAIC) for selecting the number of magnetic filling factors and compared alternative strategies for treating line strengths. In most cases, BIC, AIC, and WAIC favoured the same number of components. For some active and rapidly rotating models, BIC and AIC selected fewer components and yielded lower field estimates, while WAIC generally produced closer agreement with the input field strengths. Treating individual line strengths as free parameters underestimated the field strength by 30--50%, whereas fitting a joint element abundance together with continuum scaling recovered the input field more reliably. Our results show that Zeeman broadening diagnostics can robustly recover magnetic fields in M dwarfs, but their accuracy depends strongly on methodological choices.

Red quasars (rQSOs) have been recognized as a short-lived, early stage in the evolutionary cycle of Active Galactic Nuclei (AGN), with fundamental differences in their intrinsic properties compared to blue quasars (bQSOs). In this work, we present the first large X-ray sample of 380 rQSOs, selected from the eROSITA/SDSS-V collaboration, providing uniform X-ray detection with optical spectroscopy accros half the sky, in the German hemisphere of eROSITA. We combine X-ray imaging, optical spectroscopy, and multi-wavelength photometry to fully probe the accretion, absorption and host properties of rQSOs. Independent Component Analysis is used to reconstruct optical spectra in a data-driven and non-parametric approach, while accounting for dust reddening and host contamination. rQSOs are intrinsically X-ray weak compared to bQSOs, with a higher fraction found at low X-ray luminosities (over 50$\%$ of the rQSO sample have $L_X < 10^{43.5} \, \rm erg \, s^{-1}$). We investigate the relative X-ray strength of rQSOs by measuring the spectral slope indicator $\alpha_{OX}$. Despite their suppressed X-ray emission, rQSOs are not optically faint, but show low $\alpha_{OX}$ values, indicating weak X-ray emission relative to their bright optical continua. X-ray spectral measurements reveal large gas column densities relative to optical reddening due to dust, implying that X-ray absorption could arise from dust-free gas close to the supermassive Black Hole (BH) rather than a classical dusty torus, while the dust responsible for optical reddening likely resides on larger host-galaxy scales or is associated with dusty gas carried in disc winds. rQSOs trace a phase of suppressed BH assembly relative to stellar mass growth, suggesting that they represent a distinct evolutionary stage where BH accretion is suppressed while the host galaxy continues to grow.

Time delays are a key observable in strong gravitational lensing systems. Their theoretical expression is usually written as a sum of a geometrical delay and a Shapiro delay, with cosmology entering through angular diameter distances and a redshift prefactor. In this work we derive this structure from the exact null geodesics of the Schwarzschild-de Sitter metric. The standard formula is recovered as the leading term in a small-angle expansion, and we identify the first correction to the usual geometrical-plus-Shapiro split. Such correction does not introduce any new cosmological dependence: it corresponds instead to a higher-order correction intrinsic to the Schwarzschild part of the metric. As a consequence, up to this order, the cosmological constant enters only through the unlensed angular diameter distances and the unlensed lens-redshift prefactor.

This paper presents a comprehensive study of quasar photometric classification and redshift estimation using machine learning techniques. We cross-matched photometric data from the Dark Energy Survey Data Release 2 (DES DR2) with spectroscopic classifications from the Sloan Digital Sky Survey Data Release 16 (SDSS DR16), yielding an initial sample of 168,738 point-like objects. Using a K-Nearest Neighbors (KNN) algorithm with PSF magnitudes in the $g$, $r$, $i$, and $z$ bands, we achieved high-precision quasar/galaxy classification against stellar contaminants, reaching a recall of 0.77 at 0.99 precision. Photometric redshifts were subsequently estimated using a hybrid machine learning approach combining a Boosted Decision Tree from ANNz and a Decision Tree Regressor from scikit-learn. The resulting catalog spans redshifts from $z \approx 0.5$ to $z > 3$, with a distinct population recovered at $z \approx 4$. A stacked outlier classifier was developed to mitigate catastrophic redshift errors. The full photometric redshift sample contains 872,372 objects and remains reliable for cosmological applications at $z \approx 4$. The cleaned catalog contains 675,683 objects and is suitable for large-scale structure studies in the range $0 < z < 3$. This robustly characterized quasar catalog provides a valuable resource for future cosmological investigations.

Evolutionary models of massive stars are quasi-exclusively computed using an Eddington gray atmosphere. This approximation does not accurately describe the complex physical phenomena occurring in the atmosphere of massive stars. We aim to include state-of-the-art atmosphere models in the evolution computations of massive stars and test how the Wolf-Rayet phase is impacted. We computed the evolution of Galactic massive stars with the code STAREVOL. During the advanced phases of evolution, we applied outer boundary conditions interpolated within a grid of CMFGEN model atmospheres at each time step. The effective temperature and effective gravity were extracted from the atmosphere models. We then compared the resulting evolutionary tracks with classical calculations assuming Eddington gray atmospheres. We find that including detailed model atmospheres has a significant impact on the effective radius and temperature of the models during the later stages of the evolution. The effective temperatures of the evolution models computed with detailed model atmospheres are greatly reduced and in better agreement with observations of Wolf-Rayet stars. On the other hand, the internal structure of the models is barely affected by the choice of the atmosphere. We show that applying post-processing corrections on effective temperature and gravity is a method equivalent to our direct inclusion of atmosphere models in evolutionary calculations. The inclusion of detailed atmosphere models in the computation of evolutionary models is necessary to correctly reproduce the position of evolved massive stars in the Hertzsprung-Russell diagram. However, this has no impact on the internal and chemical evolution.

We present a comprehensive phase-space analysis of a quadratic dark energy model where the pressure includes a nonlinear term proportional to the square of the energy density. This minimal extension beyond the $\Lambda$CDM framework introduces a dynamical parameter $\eta(z)$ that governs transitions between different cosmological regimes. Through dynamical systems theory, we identify critical points and their stability properties, revealing that negative $\eta$ values drive the system toward stable phantom attractors (sinks), while positive values correspond to unstable repellers (sources). The model exhibits a distinctive asymptotic approach to the phantom divide ($w_{\rm eff}=-1$) from both quintessence and phantom sides without actual crossing, providing a non-crossing alternative to the phantom-crossing behavior preferred by recent DESI DR2 constraints. Our analysis shows that stable phantom attractors produce enhanced Hubble expansion rates and more pronounced late-time acceleration, features that can be compared with recent DESI observations suggesting evolving dark energy.

We investigate the existence and stability of highly-compact sub-stellar objects composed of strange quark matter (SQM), focusing on finite-size strangelets with baryon number $A \leq 100$. Motivated by the emergence of mass--radius outliers in the \textit{Gaia} DR3 era, we employ a Bayesian exploration of the MIT bag-model parameter space, explicitly accounting for finite-size surface and curvature contributions that become relevant at low baryon number. Enforcing the bulk absolute-stability requirement for SQM ($E/A < 930~\mathrm{MeV}$), we find that self-gravitating equilibrium sequences are confined to the sub-stellar regime, with typical masses $M \simeq 10^{-2}$--$10^{-1}\,M_{\odot}$ and characteristic radii of order $10^{3}$--$10^{4}$ km. We further show that rapid rotation, treated through a self-consistent framework that incorporates relativistic thermodynamics, can substantially inflate the equatorial radius and extend the accessible mass--radius domain. While rotation does not eliminate the intrinsic high-density compactness of these configurations, it shifts the most extended models closer to the observational parameter space of massive exoplanets. A comparison with objects from the NASA Exoplanet Archive reveals a pronounced density gap separating standard atomic-matter planets and brown dwarfs from the strangelet-rich branch predicted here. We conclude that light strangelets cannot account for solar-mass white dwarfs, but they robustly predict a previously unexplored population of ultra-compact sub-stellar objects, offering testable targets for future microlensing searches and high-cadence photometric surveys.

Chemical equilibrium calculations are a key ingredient for modelling and interpreting spectroscopic observations of (exo)planets, brown dwarfs, cool stars, and protoplanetary disks. As these applications increasingly probe non-solar elemental abundances and previously underrepresented elements, equilibrium chemistry solvers must be both numerically robust and capable of handling complex chemical systems. Here we present FastChem 4, a major update to the open-source FastChem equilibrium chemistry code. We extend the gas-phase solver with a multidimensional Newton-method that mitigates the slow convergence previously encountered for strongly non-solar elemental abundances. We further reformulate the gas-phase equations in logarithmic element densities, removing the dependence on quad-precision arithmetic and allowing FastChem to be applied at low temperatures on any platform supporting double precision. The condensate solver is upgraded with adaptive Levenberg-Marquardt regularisation, a perturbed-Hessian fallback, and a combined gas-condensate Newton solver. These changes lead to a strong increase in computational performance and stability. The thermochemical data is expanded using thermochemical data from the NIST-JANAF tables and the Barin compilation, and now comprises 800 gas-phase molecules and ions and 511 condensates spanning 44 elements. We apply the updated code to a wide pressure-temperature grid for both solar and carbon-rich (C/O = 2) elemental compositions. The resulting grids reproduce the classical solar-composition condensation sequence and reveal the marked shifts that occur under carbon-rich conditions. We also find that silicon monoxide is stable as a condensate over a limited pressure-temperature range, consistent with recent JWST observations of brown dwarfs. FastChem 4 is released under the GPLv3 licence, together with a pre-compiled Python package.

Observations of diffuse gamma rays above hundreds of TeV from the Galactic disk provide strong evidence for the existence of PeV cosmic-ray accelerators--so-called PeVatrons--in the Galaxy. However, mounting observations have ruled out most supernova remnants as likely PeVatron candidates, suggesting instead that multiple populations of cosmic-ray sources exist in the Galaxy. Recently, the LHAASO collaboration reported the detection of ultra-high-energy gamma rays from microquasars, establishing that the black holes in these systems, which accrete matter from companion stars, are powerful PeV particle accelerators. In this work, we propose a two-component source model to explain the observed cosmic-ray spectra and composition up to the PeV range. Below approximately 100 TeV, supernova remnants serve as the dominant sources; above this energy, microquasars are considered the primary candidate population. Within this scenario, the assumption of a charge-dependent cutoff well accounts for the latest measurements, including the proton and helium spectra up to the PeV range, the energy-dependent composition, and the all-particle spectrum. In contrast, the nuclei-dependent cutoff hypothesis is ruled out by the data.

Spectra of 6.7 GHz methanol masers from 21 pointings of known star-forming regions are reported. The C-band observations, using the Green Bank Telescope in full Stokes mode, have measured how polarization properties vary across the maser profiles in each spectrum and vary between different epochs of observation. Two-thirds of the sources are observed to have 6.7 GHz methanol masers, including one new detection (G240.316+0.071). Linear polarization is in the range 0 to 15% and circular polarization 0 to $\pm$9%, in line with previously reported values. The only instances where polarization is not observed is when these polarization ranges are below 3$\sigma$ detection limits. Zeeman splitting is observed in several sources, with splitting values derived from velocity separation between RCP and LCP components via Gaussian fitting. These values are seen to change with time and appear to correspond to changes in the linear and circular polarization. The polarization properties varying across the spectra and changing with time are most likely due to variations in the magnetic fields.

Star-planet magnetic interactions (SPMIs) can transfer energy from an exoplanet to its host star via Alfvén waves when the planet orbits within a sub-Alfvénic stellar wind. Similar conditions were encountered by the sungrazing comet Lovejoy as it passed through the solar corona in December 2011. The possibility that comets could trigger solar activity via magnetic interactions has never been investigated. The aim of this paper is to quantify the energetics of such interaction and to assess if enough energy could be deposited to either form a hotspot or trigger eruptions on the Sun. We used the magnetohydrodynamic WindPredict-AW model to reconstruct the coronal magnetic field and solar wind conditions along the comet's orbit, and determined the magnetic footpoints that connected the comet to the solar surface. By estimating the travel time of hypothetical Alfvén waves, we identified a brightening event observed by EUVI aboard STEREO-A that is spatially and temporally consistent with comet Lovejoy's passage. We then computed the SPMI power for all magnetic field lines anchored within 5 degrees of the brightening using SPMI power scaling laws. The resulting power distribution spans approximately 10$^{14}$-10$^{16}$ W. In comparison, the estimated radiative power of the brightening event in the STEREO-A/EUVI 195$Å$ channel is approximately 10$^{17}$ W. We find that comet Lovejoy does not generate sufficient SPMI power to be the energy source of the observed brightening intensity. However, it may still act as a perturbation to existing magnetic structures, and lead to the triggering of solar flares. Confirming this hypothesis would require additional observations that do not exist for this event, which make future passages of sungrazing comets valuable opportunities to study SPMI processes within our solar system.

Compact steep spectrum (CSS) sources generally show weak Doppler boosting, yet some exceptions show multi-year-scale radio flux variability and high-energy activity. Since 2022, the CSS quasar 3C 138 has been in a radio high state accompanied by multiple gamma-ray outbursts, offering unique opportunities to study changes in jet physical conditions. We estimated the synchrotron self-absorption (SSA) magnetic field ($B_{\rm SSA}$) in the SSA core of 3C 138 during its high state and compared it with the equipartition magnetic field ($B_{\rm eq}$) to assess the core field environment. Using extended Korean Very long-baseline interferometry Network (KVN) data at 22, 43, 86, and 129 GHz (2024-2025), we calibrated the visibilities and modeled resolved components with circular Gaussians. A single-zone SSA model fitted to the core spectrum provided the turnover frequency and peak flux density, from which we estimated the $B_{\rm SSA}$ and $B_{\rm eq}$. We used Very Large Array and Atacama Large Millimeter/submillimeter Array data to constrain the broadband spectra with the same model. The KVN SSA core shows a turnover at about 33 GHz and a peak flux of about 1.45 Jy. The inferred $B_{\rm SSA}$ is far below equipartition, with $B_{\rm SSA}/B_{\rm eq}\approx0.05$. The flux variability of 3C 138 is driven by a compact, particle-dominated core. Shock-driven particle injection in the inner jet could account for the core brightening and the production of X-ray/gamma-ray emissions through an inverse-Compton process without requiring extreme relativistic beaming effects.

[Abridged] In slow roll inflation, the power spectrum (PS) of primary gravitational waves (PGWs) generated from the quantum vacuum rises as $k^2$ over wave numbers $k$ which never leave the Hubble radius. In fact, over such small scales, the PS exhibits a similar behavior at any time after inflation. In a recent work, we had argued that the PS of PGWs has to be regularized to truncate the unphysical quadratic rise at large wave numbers. Assuming instantaneous transitions from inflation to the epochs of radiation and matter domination, we had shown that the regularized PS oscillates with a constant amplitude about a vanishing mean over small scales during these epochs. We had also smoothed the transition (actually, the `effective potential' governing the equation of motion of GWs) from inflation to radiation domination using a linear function and evaluated the regularized PS of PGWs post inflation. In such a case, we had shown that, over small scales, while the regularized PS continues to oscillate about zero, its amplitude decreases as $k^{-1}$. In this work, using the Born approximation, we examine the behavior of the regularized PS of PGWs over small scales when they are evolved through smoother and smoother transitions from inflation to the epochs of radiation and matter domination. We illustrate that, at small scales or high frequencies, the suppression in the regularized PS of PGWs occurs more and more sharply as the transition is smoothed further and further. With the help of examples, we also show that, in the case of transitions described by infinitely differentiable `effective potentials', the regularized PS of PGWs exhibits an exponential suppression on small scales. We argue that the observation of the exponential drop in the PS of PGWs can help us determine the energy scale and the time of the end of inflation. We clarify related issues and discuss the wider implications.

The geometry and morphology of symbiotic stars are conducive to exhibit a variety of scattering phenomena. The prominent among them is the Raman scattering of O VI doublet $\lambda \lambda$ 1032,1038 angstrom, which often show strongly polarized features in the visible spectrum. Similar Raman scattering of Ly$\beta$ photons has also been predicted to occur in symbiotic stars, though with fewer detections and with weak polarization amplitudes. Here, we present the discovery of strong variable polarization in the H$\alpha$ profile of a recently established symbiotic system Y Gem, over a period of nearly 22 months. This is, most likely, a very rare detection of the strongly polarized Raman scattered Ly$\beta$ photons, falling at the H$\alpha$ emission. Monte-Carlo simulations have been conducted to confirm the underlying Raman scattering process causing the polarized line profile, and a simple orbital model is constructed with typical parameters available in the recent literature along with a complementary low-resolution spectroscopic data. These simulations and models are then used to validate the observed polarization variation of H$\alpha$ at different orbital phases corresponding to the epochs of observations. The possibility of such strong variable H$\alpha$ polarization, being caused by Raman scattering of Ly$\beta$, would thus open up avenues of exploring such effects in various other astrophysical situations having similar morphology.

The redshifted 21-cm signal from neutral hydrogen (HI) in the intergalactic medium (IGM) is a powerful probe of the Epoch of Reionization (EoR). Owing to the complex growth and morphology of ionized regions, the 21-cm brightness-temperature field becomes strongly non-Gaussian during the EoR, limiting the information captured by the standard power spectrum alone. While higher-order statistics such as the bispectrum can recover part of this information, they are computationally expensive and often less straightforward to interpret. In this work, we investigate marked statistics as an alternative framework for characterizing the EoR 21-cm signal. Using semi-numerical 21-cm simulations, we introduce a set of EoR-tailored mark functions, inspired by and extending existing marked-statistics ideas, and study both the power spectrum of the mark and that of the marked field. We show that suitably chosen marks can selectively enhance contributions from different IGM environments, capture additional non-Gaussian information beyond the standard power spectrum, and improve the statistical constraining power on EoR model parameters in a Fisher-matrix analysis. These results demonstrate that Fourier-space marked statistics, including both the power spectrum of the mark and that of the marked field, provide a computationally simple and flexible extension of standard two-point statistics for extracting astrophysical information from the EoR 21-cm signal.

We present three-dimensional hydrodynamic and magnetohydrodynamic core-collapse supernova simulations of a rapidly rotating, high-compactness $39 M_\odot$ progenitor to investigate the roles of rotation and magnetic fields in shock revival and outflow morphology. This study is designed to separate neutrino-driven expansion, rotation-induced deformation, and magnetically aided polar outflow within the same progenitor. We evolve three models: a non-rotating hydrodynamic baseline, a rotating hydrodynamic model, and a rotating magnetized model. All three models reach runaway shock expansion within the simulated interval, but with markedly different morphologies and timescales. The magnetized model revives first and develops the clearest bipolar outflow. The rotating non-magnetized model undergoes the latest shock revival and remains comparatively compact at the end of the simulation. The non-rotating model also undergoes shock revival, but subsequently collapses to a black hole about one second after core bounce. In the magnetized model, Maxwell stresses redistribute angular momentum and extract energy from the differential rotation of the protoneutron star, reducing the inner-core spin and helping channel rotational free energy into the emerging polar outflow. Neutrino emission provides an additional, though smaller, angular-momentum sink in both rotating models. We find that rapid rotation and strong magnetic fields can launch an early magnetically aided polar outflow in 3D, while the resulting dynamics remain intrinsically non-axisymmetric. In this extreme progenitor, rotation also provides significant support against prompt black-hole formation, although the longer-term remnant stability remains uncertain beyond the simulated interval.

Dust properties, such as mass and porosity, impact planet formation directly. Understanding the time evolution of dust distribution across multiple properties requires numerical computation. However, available ways to calculate the multi-component coagulation-fragmentation are highly time-consuming. This study aims to develop a fast and accurate algorithm for multi-component coagulation. We assumed that two pairs of colliding aggregates reproduce a similar outcome if the dust properties are similar, and that the ratio of dust properties in logarithmic space gives the similarity as a "distance". These assumptions enable us to apply the tree algorithm, which groups distant bins and calculates interactions together, to coagulation. The algorithm reduces the computational complexity from $O (N^{2d})$ to $O (d N^d \log N)$, considering $N$ bins per $d$ components. We tested the algorithm by comparing it with the conventional direct method for cases where analytic solutions are known. We measured the dependencies of the wall-clock time, $L_2$ error in the distribution, and relative error of the total mass, on the $d, N$, opening angle $\theta_c$, and maximum dust distribution width after coagulation $k_c$. The algorithms are found to calculate coagulation consistently. For $d=1$, the tree method is faster than the direct method for a specific range of parameters. For $d=2$, however, the tree method is faster for all parameter regions surveyed, speeding it up by tens of times. Increasing $N$ and decreasing $\theta_c$ or $k_c$ made it slower and more accurate. Additionally, using a small $k_c$ performs worse than when using a large $k_c$, suggesting that limiting $k_c$ is unnecessary. We present a fast tree algorithm for the multi-component coagulation equation. It will enable us to evolve the multi-component dust distribution, such as in mass-porosity space, in protoplanetary disks.

Photometric variability in young stellar objects (YSOs) provides critical insight into the mechanisms of mass accretion, disk evolution, and circumstellar extinction in early stellar evolution. We present an analysis of day-timescale optical variability in the Orion A central region using two-night 7-Dimensional Telescope (7DT) medium-band photometry obtained on March 23 and 24, 2024. The 7DT observations provide optical spectral sampling with 16 medium-band filters spanning 400--825 nm, enabling direct two-epoch comparisons. To remove satellite-trail contamination, we used an SSIM-based ResNet classifier (accuracy 0.97; F1 = 0.93) to exclude affected exposures. Subsequent photometry and two-epoch variability measurements yielded a working sample of 769 YSO candidates, among which we identified 110 variables ($\sim$14\%), including seven extreme cases with $|\Delta m_\lambda|>0.5$ mag. To describe the wavelength dependence of the variability, we compared five simple templates: extinction-like changes ($R_V =$ 3.1 and 5.5), a gray (wavelength-independent) change, and two spot-like toy models (hot and cold) implemented as two-temperature surface mixtures. The best-fit results are dominated by spot-like templates (37 cold-spot and 22 hot-spot objects), with 37 sources best matched by extinction-like templates and 14 by the gray template. The m650 excess fraction is higher in the hot-spot and gray templates than in the others. This could be compatible with more frequent line/veiling-related contributions in those groups, although the m650 excess is not a direct accretion diagnostic.

Aims. X1908+075 is a supergiant X-ray binary (SgXB) composed of an evolved OB star and a neutron star (NS) in a 4.4-day orbit. We aim to characterize the stellar wind, constrain the system geometry, and investigate the origin and variability of the Fe K$\alpha$ fluorescence line. Methods. We analyzed three Chandra HETGS observations at different orbital phases. The continuum was modeled with a bulk motion Comptonization (Bmc) model including partial-covering absorption. We performed a blind line search with Monte Carlo simulations accounting for the look-elsewhere effect. The orbital modulation of $N_{\rm H}$ was modeled through particle swarm optimization to constrain the inclination and donor mass-loss rate. Results. Fe K$\alpha$ emission is detected in all observations and remains significant after correcting for multiple trials (global $p < 0.005$). Tentative features include a Compton shoulder and Fe K$\beta$ emission, indicating dense reprocessing material. No highly ionized Fe XXV or Fe XXVI lines are confirmed. The Fe K$\alpha$ flux correlates with the continuum flux, while its equivalent width anticorrelates with both $N_{\rm H}$ and orbital phase, opposite to the canonical curve of growth seen in many HMXBs. Line broadening reaches velocities of $\sim3000$ km s$^{-1}$. Modeling of the $N_{\rm H}$ variability yields an inclination $i = 46 \pm 3^\circ$ and a donor mass-loss rate $\dot{M}_{\rm W} = (9.1 \pm 1.6) \times 10^{-7}\,M_\odot\,{\rm yr}^{-1}$. Conclusions. X1908+075 is a classical wind-fed SgXB in the direct accretion regime. The observed continuum and line variability reflect the complex density and ionization structure of the stellar wind.

Stellar mergers produce explosive outflows that serve as transient sources of IR line luminosity and inject mechanical energy early into the natal molecular cloud. We present the first velocity-resolved maps of the [O I] 63 and 145 um fine-structure line emission from the wide-angle outflow in Orion BN/KL, the nearest explosive outflow. The data were obtained with SOFIA and include sensitive [C II] 158 um and OH maps. They allowed us to disentangle the quiescent cloud gas from the outflow, traced by a broader [O I] component with a line FWHM of about 20-30 km/s and exhibiting a spatial distribution similar to that of the shock-excited H2 emission seen with JWST. The OH 119 um line shows a prominent P-Cygni profile covering 160 km/s, similar to the very broad CO lines. The total [O I] 63 and 145 line luminosity is remarkably high, 86.5 L_sun, comparable to the H2 and CO line luminosities, implying an outflow mass-loss rate of (9.1+/-2.6)x10^-3 M_sun/yr and a mass of 3.3-5.9 M_sun. The [O I] 63 / 145 and [O I] 63 / [C II] 158 intensity ratios reach very high values in the line wings (20-30 and 40-60, respectively), exceeding those found in PDRs. These ratios are consistent with the presence of dense (10^5 to 10^6 cm^-3 ) and warm (~500 K) postshock gas. We analyzed the fine-structure line-wing intensities using magnetized shock models that include UV irradiation, to which the [C II] 158 line intensity is particularly sensitive. We find that the [O I] and [C II] intensities are consistent with emission from dissociative J-type shocks with velocities of 30-40 km/s and preshock gas densities of a few 10^4 cm^-3, illuminated by external UV radiation generated by surrounding fast shocks and possibly by massive (proto)stars in the region. We also report a broad [O I] emission feature around the BN star, which we attribute to an unresolved outflow or wind bow shock.

The use of methods that investigate the value of the Hubble constant H$_0$ in different patches (60 $ or 90$^ size) across the sky to probe the statistical isotropy of the Universe using large SNe Ia databases has led to contradictory claims of either anisotropy or isotropy. The anisotropy directions vary amongst research works. The objective of this paper is to clarify the abovementioned claims and study the lack of basis for depicting directions of anisotropy with the present SNe Ia samples. We explain the type of limitation embedded in the SN Ia lightcurve method to determine the isotropy of H_0 and the corresponding consequences. The widely used analysis through the Region Fitting and the Hemisphere Comparison methods is done here using the Pantheon+ database, simulating 2000 distinct directions in the sky within a Bayesian Markov Chain Monte Carlo approach. We also study a smaller SNe Ia database, the Carnegie Supernova Project sample, leading to a similar kind of result as that from the Pantheon+ sample. We investigate the validity of the directions found for anisotropy within these analyses. We have found that within the tests used here, the Region Fitting method and the Hemisphere Comparison method, one can not determine with robustness the direction of an anisotropy of H$_0$ using the present SNe Ia large data samples. This is intrinsic to the way H$_0$ is obtained with the SN Ia lightcurve method. Achieving robust constraints will require a quite uniform sky coverage from larger SNe Ia samples with improved systematics.

Systematic observations of the Sun are performed at Meudon observatory since 1908 under the form of monochromatic images. However, a major technical improvement occurred in 2017; since this date, spectroscopic datacubes are obtained daily with the spectroheliograph and a fast CCD camera. Line profiles of H$\alpha$ (6562.8 A), CaII K (3933.7 A), CaII H (3968.5 A) and H$\epsilon$ (3970.1 A in the wing of CaII H) are recorded over the full solar disk under the form of 3D FITS files (x, y, $\lambda$). The optical spectral resolution is 0.15 A for Calcium and H$\epsilon$ (0.093 A/pixel), and 0.25 A for H$\alpha$ (0.155 A/pixel); the spatial sampling is about 1 arc sec (the usual seeing is 2 arcsec). Datacubes are freely available since July 2017 in raw TIF (level 0) or processed FITS (level 1) format. Access to observations and typical line profiles associated to solar features are presented.

Hyperactive repeaters provide a unique window into the evolving environments and energy budgets of fast radio burst (FRB) sources, though they may not be representative of the FRB population in general. High-cadence observations are key to capturing the rarest and most energetic bursts, which occur only once per hundreds to thousands of hours. Here we present an unprecedented $4{,}200$-hour observing campaign targeting FRB 20240114A as part of the HyperFlash and ÉCLAT FRB monitoring programs. Over $806$ days, we detected $178$ high-energy ($\sim$$10^{40-42}$ erg) bursts with HyperFlash, which together amount to $4.4 \times 10^{42}$ erg of released radio energy (assuming isotropic emission and 1-GHz emission bandwidth). The cumulative energy of the HyperFlash bursts is about twice that of $\sim$$11{,}000$ lower-energy bursts detected with FAST, emphasising the significant role that the highest-energy bursts play in depleting the central engine's stored energy. In fact, the single most brilliant burst from our sample, which we term the STROOP, contributes roughly $1/3$ of all the energy we measure, and is at the maximum energy seen in studies of both repeating and apparently one-off FRBs alike. We also find a break in the burst energy distribution at $\sim$$2\times10^{40}$ erg and a linear dispersion measure (DM) increase of $+0.96 \pm 0.06$ pc cm$^{-3}$ over a period of $318$ days. We discuss these findings in the context of a magnetar source model and highlight comparisons with the energetics of intermediate and giant X-ray/$\gamma$-ray flares from Galactic sources.

NectarCAM is a Cherenkov camera designed to detect gamma rays with energies between 80 GeV and 50 TeV. It will equip nine medium-sized telescopes (MSTs) of the Cherenkov Telescope Array Observatory. NectarCAM consists of 1855 pixels distributed over 265 modules. Each pixel consists of a photomultiplier tube that is connected to a NECTAr3 chip. This NECTAr3 chip contains a 12-bit digitizer with a GHz sampling rate, and has a typical readout deadtime of ${\sim}0.7$ $\mu$s. In these proceedings, we highlight the performance of the NectarCAM in terms of time resolution and charge resolution. We also present the latest calibration techniques that were recently implemented for the camera. Finally, we briefly present the current status and plans of the NectarCAM production; the first production-line NectarCAM will be ready for shipment by Summer 2026, and it is planned to equip one of the MST pathfinders of CTAO.

We reconstruct the primordial power spectrum $P_{\mathcal{R}}(k)$ from the BOSS DR 12 LRG and eBOSS DR 16 QSO catalogs with a non-parametric Bayesian method. The $P_{\mathcal{R}}(k)$ is reconstructed by linearly interpolating $N$ knots in the $\{ \log k, \log P_{\mathcal{R}}(k) \}$ plane. We use a parametric model to describe the galaxy power spectra of the BOSS+eBOSS catalogs, assuming any power-law deviations and BAO contributions separately from the matter power spectrum template, composed of seven parameters $\Theta_{\text{model}}$. This parametric model enables us to reconstruct $P_{\mathcal{R}}(k)$ at non-linear scales, reaching $k = 0.3 \text{ h} \text{ Mpc}^{-1}$. The method is validated by applying it to different Primordial Features (PF) templates and by recovering the input power law of \textsc{MD-Patchy} and \textsc{EZmock} mock catalogs, representative of the BOSS and eBOSS data. These mocks provide additional information on $\Theta_{\text{model}}$, enabling us to impose Gaussian correlated priors on $\Theta_{\text{model}}$. This prior set allows us to reconstruct $P_{\mathcal{R}}(k)$ more precisely and to alleviate the degeneracies between the model and knot parameters. The results for both individual and combined $z$-bins and galactic caps of the BOSS and eBOSS catalogs are consistent, showing no evidence of the presence of PF in $P_{\mathcal{R}}(k)$ and pointing to a quasi-scale-invariant power law as the preferred model for $P_{\mathcal{R}}(k)$, as predicted by most slow-roll inflationary models. With a different prior set that relaxes the Gaussian constraints on $\Theta_{\text{model}}$ and imposes Planck-based priors on the extreme knots, the results also favor the power law. From the knot reconstructions, we robustly constrain the spectral index $n_s = 0.976 \pm 0.021$, compatible with the Planck value.

We present a new method for fast computation of the realization-dependent bias, a major computational bottleneck in measurements of the cosmic microwave background (CMB) lensing power spectrum. The method accelerates the bias calculation by differencing two correlated estimates: one based on fully realistic masked simulations and the other on isotropic simulations, for which the bias is analytically tractable. We show that our algorithm reduces the total computational cost of a lensing power spectrum measurement by approximately a factor of five for Atacama Cosmology Telescope- or Simons Observatory-like noise levels, or by a factor of three if current anisotropic filtering methods are left unchanged. Owing to its simplicity, the method can be readily implemented in existing CMB lensing analysis pipelines.

We present new astrophysical parameters for the open cluster NGC~6793 based on new CCD $UBV(RI)_{KC}$ photometry. We derived a reddening of $E(B-V) = 0.24 \pm 0.02$~mag and a heavy element abundance of $Z = 0.024$ ($[Fe/H] = +0.20$~dex). Padova isochrone fitting to the $V \times (B-V)$ colour-magnitude diagram yields an intermediate age of $525 \pm 51$~Myr and a distance modulus of $\mu = 8.80 \pm 0.05$~mag, corresponding to a distance of $d = 575 \pm 58$~pc from the Sun. The core radius of NGC~6793 appears to be shrinking due to advanced dynamical evolution ($\log\tau_{2} = 1.13$), driven by mass segregation and the evaporation of low-mass stars from the central region. The ratios of core to half-mass radius ($R_{c}/R_{h}$) and half-mass to Jacobi radius ($R_{h}/R_{J}$) indicate that the cluster's evolution is governed by the combined effects of internal two-body relaxation, mass segregation, and external tidal perturbations. The ratio $R_{t}/R_{J} = 0.99$ suggests that the cluster is currently in a tidally filling state. The parameter pairs ($t_{diss}/t_{rlx_{1}} = 40$, $\log R_{J}/R_{c} = 0.72$) and ($R_{h}/R_{J} = 0.38$, $\log\rho_{amb} = -0.88$) place NGC~6793 among the relatively compact clusters within $R_{GC} < 7.9$~kpc. This implies a compact internal structure that is stable against the combined effects of two-body encounters and tidal heating. Given its current state, NGC~6793 will likely dissolve and disperse before entering the final contraction phase ($R_{4}$ regime).

We describe the design, hardware integration, and calibration performance of the Wide-Field Imager (WFI) instruments for the Polarimeter to Unify the Corona and Heliosphere (PUNCH) mission. The WFI instruments are a trio of visible-light heliospheric imagers that, together, view the outer corona and solar wind from under 3.5° to over 47° from the Sun, via sunlight that is Thomson-scattered from free electrons. In flight, the WFIs are arranged so that their collective fields of view form an approximately symmetric trefoil on the sky, comprising three circular-truncated square fields spaced 120° apart in position angle. The WFIs work with the NFI instrument, described elsewhere, to implement the full PUNCH field spanning all solar position angles, at elongations from 1.5° to 47° from disk center. WFI is implemented using dioptric (lens) optics and deep multi-stage baffles that attenuate solar, planetary, and lunar stray light sufficiently for ground processing to reveal the faint signal for the primary science. WFI measures both total brightness (tB) and polarized brightness (pB), via an on-board polarizing filter wheel (PFW) and charge-coupled device (CCD) camera that share a common design with those of the NFI instrument.

Recent DESI BAO measurements, combined with CMB data, reveal a tension within the $\Lambda$CDM model that leads to a discrepancy between cosmological and laboratory bounds on the summed neutrino mass. We show that a recoupled interacting radiation component can alleviate this cosmological tension, as well as the one with neutrino oscillation experiments. Sterile neutrinos interacting through a light pseudoscalar mediator provide a concrete realization of this scenario. The resulting interacting fluid modifies the CMB phenomenology, lowers the preferred matter density, and improves the consistency between CMB and DESI BAO measurements. Combining CMB with DESI DR2 BAO measurements, we find a $2.7 \sigma$ preference for a nonzero interacting sterile neutrino component, $\nsp=0.253 \pm 0.094$, corresponding to an improvement $\Delta \chi^2=-8.98$ relative to $\Lambda$CDM. The model also reduces the tension with the SH0ES determination of the Hubble constant to the $2.4\sigma$ level.

Habitable Worlds Observatory (HWO) presents a key challenge for technology development in the coming years, requiring a $>$ $100\times$ more stable system than \textit{JWST}. WaveDriver is a concept for a laser guide star spacecraft coupled to an adaptive optics (AO) system onboard HWO that would enable HWO to reach its picometer-level wavefront stability requirements while relaxing other HWO subsystem requirements. At LLNL and UCSC we are revisiting the concept initially proposed by Douglas et al.\ (2019). We present key results key initial results from the first phase of our project, including (1) AO control developments, including with Linear Quadratic Gaussian control and machine learning, (2) AO wavefront sensor (WFS) trade study simulations, and (3) simulations of a photonic lantern natural guide star WFS. A key finding from our work is that WaveDriver could be needed to enable HWO's primary mirror segment stability and/or low order wavefront stability requirements.

Pulse profile modelling via relativistic ray-tracing can constrain the system parameters of neutron stars, notably their mass and radius. Among these objects, accreting millisecond pulsars (AMPs) are promising targets, because they are bright in X-rays and their potentially polarized radiation can lead to complementary constraints on the emission geometry. We perform combined analysis of NICER and IXPE observations of the recently discovered the 448-Hz pulsar SRGA J144459.2-604207, with IXPE providing X-ray polarization information. NICER and IXPE jointly favour a large mass and radius for our best-fitting model, for which the neutron star has two independent hotspots. The primary hotspot is centered near the northern rotational pole, the secondary in the southern hemisphere, and the observer inclination is in the range 50-75 degrees. The primary hotspot is large (up to half the surface area) and contributes the majority of the non-pulsed X-rays, while the secondary is hotter and the major contributor to the overall pulse profile shape. However, many parameters are inferred to be near the prior bounds, which could indicate that the model does not adequately account for important physics. Furthermore, we tested several different methodologies for joint analysis of the two data sets: the results are sensitive to the method used, something that merits further study with synthetic data. In the future, we expect simultaneously recorded data will lead to improved parameter constraints, especially when multi-band and polarized data are combined.

Orbital eccentricity in gravitational-wave signals from merging compact object binaries is a powerful indicator of their formation channel. Several binary black hole mergers and a neutron star--black hole merger have been reported to exhibit signs of eccentricity, but which events are identified and the significance of the eccentricity differs between studies. Measurements of eccentricity can change depending on the choice of prior. The choice of prior is subtle: eccentricity is commonly measured at an arbitrary reference frequency, which varies from study to study. We use the candidate eccentric neutron star--black hole merger GW200105_162426 as a case study, employing a range of priors and reference frequencies, and find the results to be strongly prior-driven. We show that the varied results reported across different studies can be partially reconciled by accounting for the evolution of eccentricity with reference frequency. In order to make conclusive statements about eccentricity, we propose a detection statistic that does not depend on reference frequency, and which marginalises over astrophysically-motivated distributions in eccentricity. Using this detection statistic, we find reduced support for the eccentric hypothesis for GW200105_162426: we obtain a natural log Bayes factor ln B $\leq$ 0.9 comparing the eccentric, aligned-spin hypothesis to the quasi-circular, precessing hypothesis. Our results cast doubt on the eccentric interpretation of GW200105_162426 and underscore the importance of modelling the astrophysical distributions of eccentricity in nature.

Recent cosmological tests have discovered a fresh new set of anomalies in the large-scale isotropy of the universe. Motivated thus by the numerous pieces of evidence for large-scale cosmic isotropy violation with the advent of the 'precision cosmology' era, we are led to explore the viability of anisotropic Thurston geometries, described in William Thurston's geometrization conjecture. In this work, we examine the coherent temperature and polarization signals generated in the CMB sky by such geometries. We begin with introducing Thurston spacetimes as our background model and the formalism we use to obtain the patterns. We then construct a set of transfer equations relative to a given background and solve them for each spacetime geometry. We finally discuss the role of spatial curvature in these FLRW limiting models along with their underlying geometry, and attempt to establish some general results on the symmetries of the patterns produced by their time evolution in terms of the Stokes parameters P, Q, U and V. We show the evolution of temperature and polarization amplitudes in terms of such Stokes parameters at different timestamps and attempt to isolate individual Thurston geometries.

Continuous gravitational waves, i.e., persistent and nearly-monochromatic signals emitted by asymmetric spinning neutron stars, remain elusive. Searches for these signals from unknown binary systems are the most computationally challenging, but they are essential, given that binary accretion provides a natural mechanism for creating the required asymmetry, and around half of the known pulsars rotating above 25 Hz are part of a binary system. Here we report on a search of a large uncharted parameter-space region: for the first time we cover gravitational-wave frequencies above 520 Hz (from 50 to 1000 Hz), and, for the first time with advanced detectors, orbital periods lower than 3 days are explored. No signal is detected, and we set the most stringent constraints to date on the amplitude of signals of this kind. Our results exclude with $95\%$ confidence neutron stars within 100 pc and rotating faster than $\sim$ 495 Hz from having ellipticities above $5.2 \times 10^{-8}$. Within the same distance our results also exclude r-mode amplitudes above $1.5 \times 10^{-6}$ for stars rotating faster than $\sim$ 740 Hz.

We study the timelike and null geodesic structure of a static, spherically symmetric black hole sourced by a Kalb--Ramond (KR) field coupled to nonlinear electrodynamics (NED). The geometry is characterized by the mass $M$, the magnetic monopole charge $q$, and the Lorentz-violating parameters $(\gamma,\lambda)$. Closed-form expressions are derived for the effective potential, as well as the specific energy and angular momentum of massive particles on circular orbits. We further analyze the photon sphere, black hole shadow, and the Lyapunov exponent associated with unstable null circular geodesics. The latter determines the eikonal quasinormal-mode frequencies through $\omega_{\rm eik}=(\ell+1/2)\,\Omega_c-i(n+1/2)\,|\lambda_L|$. The shadow radius is compared with the Event Horizon Telescope (EHT) observations of M87$^\ast$ and Sgr~A$^\ast$, allowing us to identify the viable region in the $(q,\gamma)$ parameter space. Finally, we compute the Hawking temperature, horizon area, and the Gray--Visser sparsity parameter. We demonstrate that the combined effects of the KR field and magnetic monopole charge increase the sparsity parameter from the Schwarzschild value $16\pi^3 \simeq 496$ to nearly $1.7\times10^3$. This indicates a significantly sparser Hawking cascade compared to the Schwarzschild case, while the photon ring remains consistent with the EHT $1\sigma$ observational bounds across most of the physically allowed parameter range.

We propose a natural resolution to the fine-tuning problem inherent in the freeze-in dark matter paradigm by embedding a sterile singlet within a 3-3-1 electroweak extension. By imposing an exact $Z_{13}$ discrete gauge symmetry, we formally suppress all low-dimensional portals to ensure that the dark sector communicates with the Standard Model (SM) exclusively through a dimension-six operator. This theoretical structure allows the extraordinarily small coupling required for dark matter production to emerge naturally from the profound hierarchy between the electroweak scale and the ultra-high Peccei-Quinn symmetry breaking scale. Detailed numerical integration of the Boltzmann equations demonstrates that the sterile singlet can be produced via the infrared freeze-in mechanism to match the observed relic abundance of $\Omega_S h^2 = 0.12$. The resulting keV-scale warm dark matter candidate remains consistent with stringent Lyman-alpha forest constraints while offering a viable solution to galactic-scale discrepancies such as the cusp-core and missing satellites problems. Ultimately, this framework provides a self-consistent unification of dark matter genesis and the strong CP solution that is completely independent of ad hoc parameter adjustments.

Diffusion and flow-based models are ubiquitously used for generative modelling and density estimation. They admit a deterministic probability flow ordinary differential equation (PF-ODE), analogous to continuous normalizing flows (CNFs), which describes the transport of the probability mass. Obtaining the likelihood from these models is of interest to many workflows, especially Bayesian analysis, and requires solving the trace of the Jacobian to compute the divergence of the learned PF-ODE, which is either $\mathcal{O}(D^2)$ to compute exactly or $\mathcal{O}(D)$ with a noisy estimate. We introduce StAD, a new distillation method to predict and learn the divergence of the PF-ODE using the Langevin-Stein operator without ever computing the Jacobian. We show that our method is competitive with the Hutchinson and Hutch++ on CIFAR-10, ImageNet and other density estimation tasks, consistently improving the variance and speed of the likelihood predictions compared to the Hutchinson. We additionally show our method will generalize to a varied class of generative models, and show that under some regularity conditions these learned vector fields can be made to satisfy the Stein class.

We present a full 3D numerical evolution code to study neutron stars in massive-scalar-tensor theories. The code is embedded in the Einstein Toolkit framework and its implementation constitutes a modified version of the Baumgarte-Shapiro-Shibata-Nakamura formalism with an additional nonminimally coupled scalar field. The approach we follow preserves the standard hydrodynamic evolution for matter fields, allowing eventually for a straightforward inclusion of more microphysical effects and better flexibility. Using this code, we examine the gravitational collapse of rapidly rotating, scalarized neutron stars to a black hole by exploring the influence of the scalar field on the dynamical features of the process and on the gravitational-wave emission. We find that for the configurations studied in this work, there is an observational degeneracy in the tensorial gravitational-wave emission between collapsing scalarized stars and their counterparts in general relativity. However, this degeneracy can be broken through the emission of scalar radiation, which carries an energy of ~10^-3 M_sun c^2. This is orders of magnitude higher than the quadrupolar emission (~10^-7 M_sun c^2) and might be used as an observational probe of modified gravity. We also find that rapid rotation can enhance this signal, since fast rotating stars can sustain larger scalar field amplitudes.

We present a bidirectional internal squeezing scheme for gravitational-wave detectors and show that it saturates the lowest known lower bounds on quantum noise from internal optical dissipation. The scheme uses two optical parametric amplification stages inside the signal-extraction cavity that act on intra-cavity fields propagating in opposite directions. Thereby, most vacuum fields entering the interferometer are squeezed, while the signal and internal vacuum fields are amplified so that loss in the readout path adds no further noise. We show that the resulting signal-referred quantum noise spectral density is independent of the arm-cavity input and signal-extraction transmissivities at high frequencies, opening design freedom to mitigate technical constraints and radiation-pressure noise. We derive these results analytically, compare them with other internal squeezing and amplification schemes, and validate the full quantum-noise spectrum through numerical simulations. We also assess realistic implementations, including dissipation mechanisms and transverse mode mismatch introduced by the scheme, and find that 'mode healing' in the signal-extraction cavity can suppress mismatch losses. These results identify bidirectional internal squeezing as a possible upgrade path for gravitational-wave observatories such as LIGO, and the scheme may also benefit future observatories and other interferometry experiments.

In this study, we precisely evaluated the feasibility of the chaotic inflation model within a non-canonical kinetic framework. By applying the slow-roll approximation and imposing constraints on the equilateral non-Gaussianity $f_{\rm NL}^{\rm equil}$, we imposed constraints on the feasible range of the potential index $n$. We established physical bounds for the non-canonical parameter $\alpha$. To obtain precise parameter constraints, we solved the primordial perturbation equations numerically and conducted a rigorous MCMC analysis by using a comprehensive joint P-ACT-LB-BK18 dataset. For these potentials $n=1/3$, $2/3$, and $1$, our results respectively tightly limit $\alpha$ to the levels of $8.8^{+1.6}_{-2.8}$, $11.7^{+1.7}_{-2.6}$, and $16.4^{+3.7}_{-7.0}$, within the corresponding $1\sigma$ confidence intervals. Meanwhile, the required number of $e$-foldings naturally converges to $N \simeq 54$, without the need for fine-tuning. These findings confirm that non-standard mechanisms can resurrect excluded chaotic inflation models within the $1\sigma$ allowed regions of high-precision cosmological data.

In this paper, we investigate background-hierarchy bounds in Type~3 of New General Relativity (NGR). These bounds arise when the contribution associated with the evolution of the background spacetime exceeds that of the quadratic kinetic term in the perturbed Lagrangian. Type~3 of NGR has two free parameters and preserves diffeomorphism invariance and spatial rotations, while breaking Lorentz-boost invariance. We first review Type~3 and identify preferable gauge choices for metric-affine gauge theories of gravity with Weitzenböck connection, including NGR, from the viewpoint of symmetry in both Dirac--Bergmann analysis and linear perturbation theory. We then revisit the perturbative analysis of Type~3 and show that the propagating modes are correctly identified even when the perturbed Lagrangian is not written solely in terms of gauge-invariant variables. Finally, we derive the background-hierarchy bounds for the scalar, transverse-vector, and tensor modes around a flat FLRW background, and identify the region of parameter space in which the linear perturbation theory of Type~3 remains viable for cosmological applications.

We investigate the warm inflationary scenario in the Weyl geometric gravity theory, in which the action is constructed by adding matter to the simplest conformally invariant gravitational action in Weyl geometry. The $\tilde{R}^2$ theory can be formulated equivalently as a linear theory supplemented by an additional scalar degree of freedom originating from higher-order curvature terms, with the equations of motion obtained via variational methods. We investigate the cosmological implications of the theory by considering the warm inflationary scenario of the early evolution of the Universe, in which radiation, the inflaton field, and the Weyl vector coexist. We consider the widely studied linear dissipation coefficient model along with a quartic potential, and investigate the influence of the Weyl vector term on the dynamics. We have performed numerical computations for different coupling models, and we have successfully developed a warm inflationary model in which the Universe transitions naturally from an inflationary epoch to a radiation-dominated era. The relevant cosmological observables have been calculated and compared with the latest observational constraints from the ACT data.

This paper describes the setup and the results of the direct on-wafer measurements of a FET noise parameters obtained with a source-pull method at temperatures down to T=4K and in the 5-12 GHz frequency range. The setup consists of a cryostat with wafer probes, two reflectometers, a programmable impedance generator, wideband isolators and bias tees and low noise preamplifier, all cooled to cryogenic temperatures, allowing to perform a full vector error-corrected wafer-level measurements of the discrete transistors and amplifier dies. The setup and its calibration procedure are designed in a such way that allows simultaneous calibration, S-parameters, noise parameters and I-V curve measurements of several FETs all in one cooldown. Using the described setup we perform first measurements of 14nm FinFETs and also measure noise parameters of an LNA based on these FETs. Resulting noise temperature values are compared against those obtained using independent and alternative measurement techniques.

Second-generation Starlink Direct-to-Cell (DTC) satellites carry an additional payload for direct cellular phone connectivity whose unintended electromagnetic radiation (UEMR) at sub-300 MHz frequencies has not been individually characterised. We reanalyse 112,534 detections from 1,806 Starlink satellites observed with the Engineering Development Array version 2 (EDA2) at 21 frequencies between 72.685 and 234.375 MHz (Grigg et al. 2025), separating 175 DTC and 1,623 Ku-only v2-Mini comparison satellites via the McDowell General Catalogue (McDowell 2020). DTC satellites emit a range-corrected flux density 1.45x that of the Ku-only comparison (Cliff's delta = +0.30, p = 2.6e-11). At 230.469 MHz the XX detection fraction reaches 0.811 against a 0.481 baseline (p ~ 1e-274), and 11 of 21 frequency channels show Benjamini-Hochberg-significant polarisation anomalies. The DTC population is brighter in eclipse than in sunlight (illuminated/eclipsed flux density ratio 0.47) while the Ku-only comparison shows the opposite sense (1.18); the reversal persists across altitude, sub-satellite latitude, frequency, and launch-epoch matching. The reversal strongly disfavours UEMR mechanisms that scale monotonically with instantaneous solar photocurrent and favours an active on-board source whose effective duty cycle is larger at lower equilibrium temperature. Within the 230.469 MHz coarse channel, fine-channel inspection isolates the excess to a single ~24 kHz bin near 230.627 MHz, tail-driven and absent at five control channels. Three falsifiable mechanism-discrimination tests show this feature is not coincident with the LOFAR-resolved Bassa et al. (2024) clock fundamentals, is unresolved at the EDA2 24 kHz resolution, and is heterogeneously expressed across the v2-Mini fleet rather than driven by a few permanently bright units or by uniform thermal scaling.

Light, weakly coupled sectors can retain information about the cosmological background in which they are produced. We study light axions produced during reheating and their contribution to dark radiation, $\Delta N_{\rm eff}$. We develop a shift-symmetric EFT in which an inflaton-dependent axion kinetic term systematically organizes the leading production channels. The same kinetic function generates both direct inflaton decay and inflaton annihilation from the oscillating inflaton background. Direct decay is described by an invisible inflaton branching fraction, while annihilation is a genuinely reheating-sensitive source controlled by a coherent combination of Wilson coefficients. We derive the contribution to $\Delta N_{\rm eff}$ from both channels and show that they scale oppositely with the reheating temperature: the decay contribution falls as $T_{\rm rh}^{-2}$, whereas the annihilation contribution grows approximately as $T_{\rm rh}^{4/3}$. Their crossing is missed by treatments that keep only one production channel. We translate current and projected $\Delta N_{\rm eff}$ sensitivities into constraints on the Wilson coefficients of the kinetic function, obtaining a two-dimensional EFT map of axion dark radiation from reheating. This map can imply both lower and upper bounds on the reheating temperature, showing that light axion relics can turn dark radiation measurements into constraints on reheating.

Ultralight bosons with self-interactions, such as axion-like particles, can form astrophysical Bose-Einstein condensates around stars or compact objects, often referred to as gravitational atoms. In this work, we adopt a recently proposed dynamical formation mechanism for these halos and estimate their impact on extreme- and intermediate-mass-ratio inspirals when present around the primary black hole. We show that, for signal-to-noise ratios $\lesssim 100$, LISA can distinguish gravitational waveforms from binaries embedded in such halo overdensities. Our analysis indicates that LISA can probe boson masses $m_\mathrm{dm}\sim10^{-17}$-$10^{-15}\,\mathrm{eV}$ and decay constants $f_a\sim10^{10}$-$3.2 \times 10^{12}\,\mathrm{GeV}$ using binaries with total masses $M\sim10^4$-$10^5\,M_\odot$, assuming conservative DM densities consistent with the central values of Navarro-Frenk-White profiles. Allowing for higher background densities and different extreme-mass-ratio configurations further extends the accessible parameter space. Moreover, we find that for a binary configuration with $M\sim10^4\,M_\odot$, $\rho_\mathrm{dm} = 10^3\,\mathrm{GeV/cm^3}$, and signal-to-noise ratio $\mathrm{SNR}\sim20$, a particle mass of $m_\mathrm{dm}\sim2.5\times10^{-16}\,\mathrm{eV}$ and decay constant $f_a\sim6.3\times10^{10}\,\mathrm{GeV}$ maximize the dephasing due to dynamical friction, enabling the recovery of the particle parameters at the percent level. These results demonstrate that LISA can place constraints on axion-like particle masses and self-interactions without requiring additional couplings to Standard Model fields.

In this work we address the possibility of having a positive running of the spectral index in inflationary theories. The recent ACT data indicate mildly that the running of the spectral index might be positive, and several other physical indications point out this possibility. If the running of the spectral index is confirmed to be positive by future cosmic microwave background experiments, this can rule out quite popular inflationary scenarios. We investigate how it is possible to obtain a positive running of the spectral index in the context of minimally coupled scalar field gravity and modified gravity. For the modified gravity we choose two mainstream and of string origin candidate theories, $F(R)$ gravity and Einstein-Gauss-Bonnet gravity. In the case of scalar field inflation and $F(R)$ gravity inflation, we demonstrate the difficulties for obtaining a positive running of the spectral index for a viable inflationary regime, so scalar theories and $F(R)$ gravity are mostly compatible with the Planck data. But nuanced scalar field scenarios can be compatible with the ACT data and produce a positive running of the spectral index. In the context of Einstein-Gauss-Bonnet theories which are compatible with the GW170817 event, the running of the spectral index can easily be positive while in parallel having a viable inflationary era.

Kerr quasinormal modes near algebraically special frequencies exhibit anomalous behavior, including apparent bifurcation, disappearance, and a nonsmooth connection to the Schwarzschild limit, which has remained puzzling for decades. Tracking poles and zeros of Green-function building blocks across different Riemann sheets, we show that the bifurcation is due to an avoided crossing accompanied by resonant excitation, while the disappearance is due to pole skipping caused by cancellation of a quasinormal-mode pole by a Matsubara-mode zero. This resolves the physical origin of these long-standing anomalies.

The unique resonance features in the $^{12}$C+$^{12}$C fusion reaction lead to significant fluctuations in the branching ratio $R_{p/\alpha}=\sigma_p/\sigma_\alpha$, making it difficult to determine the $R_{p/\alpha}$ at astrophysical energies. By combining Hauser--Feshbach statistical-model calculations with constraints from direct charged-particle and gamma-ray measurements, we investigate the energy dependence of the averaged $R_{p/\alpha}$ and predict its behavior within the Gamow window. Owing to the strong energy dependence of $R_{p/\alpha}$, the corresponding reaction-rate ratios, $\langle \sigma v \rangle_p / \langle \sigma v \rangle_\alpha$, during core and shell carbon burning are determined to be 0.29, 0.45, and 0.52 at $T_9 = 0.5$, 1.0, and 1.2, respectively, significantly lower than the widely adopted CF88 constant value of 0.78. The implications of the revised $\langle \sigma v \rangle_p / \langle \sigma v \rangle_\alpha$ ratio for stellar nucleosynthesis and white-dwarf evolution are also discussed.

This work investigates the impact of seasonal variations in seismic noise on the low-frequency performance of the Einstein Telescope (ET) at the Sardinia candidate site, focusing on implications for compact binary coalescence observations. Using seismic data collected between 2022 and 2025 in deep boreholes, we characterize monthly noise variations and identify representative best and worst case scenarios, corresponding to July and December. The measured seismic spectra are used to estimate the Newtonian noise contribution in the 2-10 Hz band and to derive modified ET sensitivity curves. These are implemented in a simulation framework to evaluate their effect on the signal-to-noise ratio (SNR) of binary neutron star and intermediate mass black hole signals, assuming the triangular ET configuration. We find that the low seismic noise of the Sardinia site results in only minor seasonal variations in detector sensitivity. The corresponding impact on SNR is limited to a few percent, even without including Newtonian noise mitigation. These results indicate that seasonal environmental fluctuation have a minor effect on the early inspired detectability of compact binaries, confirming the suitability of the Sardinia site for achieving ET low-frequency sensitivity goals.

We study scalar field theories that break diffeomorphism invariance down to the subgroup of transverse diffeomorphisms through the matter sector in cosmological backgrounds. We focus on single- and multi-field models and develop the corresponding cosmological perturbation theory. We analyze the different contributions to the pressure perturbation, discussing the adiabaticity and the effects in the perturbation coefficients of the interactions that arise in the multi-field case as a consequence of the symmetry breaking. We also consider the stability of the perturbations in terms of the effective speed of sound and present particular models that could be of phenomenological interest.

Gravitational-wave observations provide a powerful probe of compact objects and strong-field gravity. In this work, we investigate the detectability of binaries containing (sub-)solar-mass black holes and superspinars with current and next-generation gravitational-wave observatories. Such objects may arise from primordial formation channels or from more exotic high-energy scenarios, and their detection would provide important insights into the population of low-mass compact objects and the physics of extreme gravitational fields. We model the gravitational-wave signals using the frequency-domain post-Newtonian inspiral waveform model TaylorF2, and truncate the signal at the innermost stable circular orbit (ISCO) to avoid contamination from the post-inspiral regime. We assess the observability of these systems using the sensitivities of current detectors such as Advanced LIGO and upcoming third-generation observatories including the Einstein Telescope and Cosmic Explorer. Our results show that while current detectors have limited reach for very low-mass binaries, third-generation observatories can enhance both detection capability and parameter-estimation precision. Their improved strain sensitivity and extended low-frequency coverage allow these observatories to track the inspiral phase over a substantially larger number of gravitational-wave cycles. As a result, they achieve considerably higher signal-to-noise ratios and provide dramatically improved constraints on binary parameters. In particular, it is possible to measure the primary spin parameter with precision $\Delta \chi_{1z}~\sim~10^{-4}-10^{-3}$, potentially allowing clear observational discrimination between near-extremal black holes and superspinars in the mass range $0.1~M_\odot-2~M_\odot$ and with signal-to-noise ratio of $\sim 100-350$.

We investigate Bayesian chiral effective field theory ($\chi$EFT) uncertainties, which assign a statistical interpretation to equation of state (EOS) distributions near nuclear saturation density, n$_0$, as well as constraints from perturbative quantum chromodynamics (pQCD) to Bayesian EOS inference from LIGO/Virgo, NICER and pulsar mass observations. The tails of the $\chi$EFT uncertainties allow for broader pressure ranges in our priors, but large parts of these are excluded by the astrophysical observations, so that the EOS and the resulting mass-radius posteriors are still very consistent with our earlier work. Within our broad prior ranges, we observe a clear stiffening of the EOS at $n \gtrsim 3 n_0$. Moreover, the impact of the pQCD constraints on the posterior EOS and mass-radius range is negligible due to the astrophysics constraints. Exploiting the strong correlation between pure neutron matter and matter in beta equilibrium, we infer the symmetry energy slope parameter $L$ from astrophysics. For the $68\%$ credible interval, we obtain $L=42.6-52$ MeV and $L=44.2-56.7$ MeV using piecewise-polytrope and speed-of-sound high-density extensions, respectively. The $L$ posterior is mainly driven by the combination of GW170817 LIGO/Virgo and PSR J0740+6620, PSR J0437-4715, and PSR J0614-3329 NICER observations.

Inverse problems are ubiquitous in modern scientific studies and involve recovering an underlying signal from noisy observations often transformed by a measurement operator. These problems are frequently ill-posed, particularly in imaging, leading to multiple plausible solutions and considerable uncertainty in reconstructed images. In fields like the physical and biological sciences, accurate uncertainty quantification (UQ) is critical for trustworthy scientific analyses and confident diagnoses. Current UQ methods for imaging often fall short; they can be inaccurate, or require unavailable or difficult-to-acquire ground truth data for calibration, which can introduce hidden biases due to distribution shifts between calibration and observed data. We introduce a UQ approach that leverages equivariant bootstrapping to generate heuristic coverages by exploiting data symmetries. We then refine these coverages through a conformal prediction calibration step, while crucially employing a self-supervised approach to avoid the need for ground truth calibration data. We demonstrate this method with weak lensing mass-mapping, where we aim to reconstruct the convergence field from shear measurements of distant galaxies weakly-lensed by gravitational fields. Mass-mapping in particular benefits from the self-supervised approach, as simulating calibration data is expensive and relies on specific cosmological models that could introduce biases in downstream cosmological inference tasks.

The exponential growth of scientific submissions has strained the peer review system. Despite the rapidly expanding global pool of researchers, this unprecedented scale has rendered the previous approach of manual expert identification unfeasible. Therefore, institutions have naturally turned to Large Language Models (LLMs) to automate intricate processes like expert reviewer identification. However, the reliability of these new models in accurately identifying domain experts lacks rigorous evaluation. We conduct a comprehensive empirical evaluation of statistical and AI-driven expertise identification methodologies to benchmark their reliability and limitations. Framing expert identification as an information retrieval problem, we utilize the distributed peer review system of a major international astronomical observatory, where proposal authorship serves as our proxy ground truth for domain expertise. Evaluating six retrieval methodologies utilized across observatories and computer science conferences, we demonstrate that traditional statistical representations outperform generative AI. Specifically, Term Frequency-Inverse Document Frequency successfully identified a labeled expert within the top 25 recommendations 79.5% of the time, compared to 51.5% for GPT-4o mini. Our results highlight that distinguishing subfield expertise requires fine-grained vocabulary, which is obscured by the semantic smoothing in generative methods. By establishing a rigorous evaluation framework for automated peer review, we demonstrate that transparent and reproducible statistical representations still outperform computationally expensive LLMs in specialized scientific tasks.