Papers for Thursday, Apr 09 2026

The VASCO project has identified over 100,000 sub-second optical transients on photographic plates from the First Palomar Observatory Sky Survey (1949-1957), all predating artificial satellites. Cann (2026a) established that transient detection rates are dose-dependently suppressed during geomagnetic storms (Z = -3.391, p = 0.0007), ruling out emulsion defects and confirming the transients as real, magnetospherically coupled phenomena. Villarroel et al. (2022) constrained the source altitude to ~42,000 km (geosynchronous orbit) through an Earth-shadow deficit. This paper presents two results. First, a pre-registered empirical test reveals the full temporal recovery profile: transient rates remain suppressed at 55% of baseline during days 7-21 post-storm, then rise to 309% of baseline during days 25-45 (p = 0.00066, Wilcoxon rank-sum; all robustness checks significant). Combined with the dose-response staircase, the overall significance reaches 3.6-4.7 sigma (Fisher's method, range reflecting sensitivity to the independence assumption). Second, we propose a candidate physical mechanism: storm-enhanced electromagnetic trapping of charged micrometeoroid dust at L ~ 6.6, followed by aggregation of icy cometary grains under restored cold plasmaspheric conditions. A flux dilution analysis demonstrates that specular reflection from a partially reflective icy aggregate only 1-4 m in diameter suffices to produce the observed plate magnitude at 42,000 km. This mechanism connects the VASCO transients to independently observed magnetospheric dust swarms correlated with geomagnetic activity (Sommer 2024) and explains the extinction of the transient population following the onset of the space age. Multi-site replication is required to confirm these results.

Stellar astrophysics relies critically on accurate descriptions of the physical conditions inside stars. Traditional solvers such as \texttt{MESA} (Modules for Experiments in Stellar Astrophysics), which employ adaptive finite-difference methods, can become computationally expensive and challenging to scale for large stellar population synthesis ($>10^9$ stars). In this work, we present an self-supervised physics-informed neural network (PINN) framework that provides a mesh-free and fully differentiable approach to solving the stellar structure equations under hydrostatic and thermal equilibrium. The model takes as input the stellar boundary conditions (at the center and surface) together with the chemical composition, and learns continuous radial profiles for mass $M_r(r)$, pressure $P(r)$, density $\rho(r)$, temperature $T(r)$, and luminosity $L_r(r)$ by enforcing the governing structure equations through physics-based loss terms. To incorporate realistic microphysics, we introduce auxiliary neural networks that approximate the equation of state and opacity tables as smooth, differentiable functions of the local thermodynamic state. These surrogates replace traditional tabulated inputs and enable end-to-end training. Once trained for a given star, the model produces continuous solutions across the entire radial domain without requiring discretization or interpolation. Validation against benchmark \texttt{MESA} models across a range of stellar masses yields a Mean Relative Absolute Error of $3.06\%$ and an average $R^2$ score of $99.98\%$. To our knowledge, this is the first demonstration that the stellar structure equations can be solved in a fully self-supervised and data-free fashion employing PINNs. This work establishes a foundation for scalable, physics-informed emulation of stellar interiors and opens the door to future extensions toward time-dependent stellar evolution.

The frequency band between $0.1$ and $10\mathrm{Hz}$ remains largely unexplored in gravitational-wave astronomy due to strong seismic, Newtonian, and suspension thermal noise that limit ground-based detectors. The Cryogenic sub-Hz cROss torsion-bar detector with quantum NOn-demolition Speed meter (CHRONOS) is a novel detector concept designed to access this frequency range from the ground. CHRONOS combines cryogenic torsion-bar test masses with a triangular Sagnac interferometer implementing a speed-meter readout, which suppresses quantum radiation-pressure noise and enables quantum non-demolition measurements in the sub-Hz regime. The detector targets a strain sensitivity of $h \sim 10^{-18}\mathrm{Hz^{-1/2}}$ around $2\mathrm{Hz}$ and stochastic gravitational wave background of $\Omega_{GW} \sim 2\times 10^{-3}$ at $2\mathrm{Hz}$. This sensitivity opens a new observational window between space-based detectors such as LISA and ground-based interferometers, enabling observations of intermediate-mass black hole binaries, searches for stochastic gravitational-wave backgrounds, and tests of macroscopic quantum measurements.

Gas properties in the outer intracluster medium (ICM) are not well-constrained, as traditional probes lose sensitivity at Mpc distances. We show that the morphology of radio relics effectively encodes the power spectrum of the surrounding density fluctuations, and that they hence represent a new observational window. To demonstrate this, we use cosmologically motivated shock-tube simulations in which we systematically vary the coherence length, amplitude, and power-law slope of the upstream density power spectra. We then post-process our simulations with the cosmic ray electron spectral solver, Crest, thereby producing a suite of mock radio relics. We find that the downstream morphology of our simulated relics is independently sensitive to each of the aforementioned parameters. Specifically, we show that observed 'double strand' features can be formed by curved shock fronts in projection, and that the scale of these features maps directly to the fluctuation coherence length. Increasing the fluctuation amplitude, meanwhile, progressively lengthens the downstream extent of the relic, thus explaining why relics are observed to be broader than the idealised expectation. It also broadens the Mach number distribution across the shock, which simultaneously increases the integrated radio flux density and produces patchier emission. Finally, steepening the power-law slope makes 'double strand' features more likely, and additionally increases both the number of radio filaments oriented parallel to the shock front and their spacing. At higher Mach numbers, steepening the power-law slope can further lead to the production of curved radio filaments, which trace large eddies. We apply our analysis to the Toothbrush and Sausage relics, and find evidence for a typical fluctuation coherence length of ~500 kpc, a non-uniform amplitude, and power spectra that are steeper than a Kolmogorov-like scaling.

We use controlled N-body simulations to study the collisional exchange of energy between stars and dark matter in ultra-faint galaxies. We find that dynamical friction between stars and subsolar-mass dark matter particles results in the depletion of dark matter from the galaxies' centers, thereby transforming dark matter cusps into constant-density cores. The process is particularly effective in tidally limited galaxies with low stellar velocity dispersion. As high-mass stars sink toward the center of the dark matter halo, the dynamical-to-stellar mass ratio within the stellar half-light radius decreases monotonically. The stellar population of a dark matter-dominated galaxy is thereby compacted into a dense, baryon-dominated cluster, surrounded by a dark matter halo. Such a cluster would share the chemical composition of an ultra-faint galaxy, yet would be virtually dark matter-free within its half-light radius. We moreover find that the collisional cooling with dark matter particles provides an efficient pathway for the formation of stellar binaries in the contracting cluster. The contraction is eventually slowed down due to the decreasing central dark matter densities and the formation of stellar binaries. Our models highlight that the dynamical processes governing the faintest galaxies give rise to a rich phenomenology, blurring the line between the dynamics of globular clusters and galaxies.

We present JWST/MIRI mid-infrared integral field spectroscopy combined with ALMA CO(2-1) observations of the post-starburst galaxy ESO 420-G13, hosting a low-luminosity AGN. The unprecedented spatial and spectral resolution of MIRI enables a detailed study of the molecular and ionised gas kinematics, excitation, and energetics in the nuclear kiloparsec, revealing the impact of AGN feedback in a system with modest radiative output. Despite its faint radio and X-ray emission ($L_{2-10keV} \sim 10^{40}$ erg/s), ESO 420-G13 exhibits powerful kinetic feedback in the form of massive molecular and ionised gas outflows, with a total kinetic power of $\sim 1.5 \times 10^{41}$ erg/s. This corresponds to a jet-ISM coupling efficiency of ~3.8%, within the range observed in more powerful AGN. The feedback is driven by a previously undetected compact jet, traced by collimated coronal-line and extended X-ray emission to >870 pc from the nucleus. The interaction is strongest ~370 pc north of the nucleus, where a fast ionised gas stream emerges perpendicular to the jet axis, coinciding with a bend in the jet direction. Enhanced velocity dispersion in warm H2 surrounds this gas stream, consistent with an expanding molecular bubble. Massive molecular outflows are detected at its edges; the blueshifted outflow is devoid of CO emission, likely due to CO destruction in shocks or by cosmic rays from the jet-ISM interaction. About 5% of the central molecular reservoir has already been expelled, and the remaining gas is turbulent and warm, suggesting an ongoing phase of AGN-driven feedback in this post-starburst galaxy. Our results highlight the enormous potential of mid-IR imaging spectroscopy to uncover jet-driven feedback in low-luminosity AGN. Without the spatially resolved MIRI diagnostics, the kinetic power of the AGN in ESO 420-G13 and its role in shaping the host galaxy ISM would have remained hidden.

Neutron stars, at their cores, are highly dense and, thus, are expected to have a number of exotic processes. This includes a possible phase transition to deconfined quark matter at the core, leading to a hybrid star. The quark matter is expected to additionally be color superconducting. The physics of superconductivity plays an important role in understanding the high density matter in the interiors of neutron/hybrid stars. At their high densities, additionally, both proton superconductivity and neutron superfluidity are expected. We study the effect of superconducting (quark/proton) matter, along with the internal magnetic field, leading to pressure anisotropy within hybrid stars. We aim to probe the effect of superconductivity, especially from color superconducting quarks, in hybrid star structure. We propose new phenomenological model anisotropy profiles within a one-dimensional framework. We model quark matter using the vector interaction enhanced Bag model, and hadron matter with the DD2 equation of state. A Maxwell construction joins both phases. We further investigate the possible observational signatures of these hybrid stars. These include mass enhancement and continuous gravitational waves, possibly arising from the anisotropy induced deformation, helping us further constrain our model and its physical parameters.

We present a parameterized ("toy") model for shock interaction and $\gamma$-ray emission in classical novae, in which a white dwarf envelope of mass $M_{\rm env}$ is removed over a timescale $\tau$ (proportional to the nova speed class, $t_{2}$) in an outflow that accelerates on the same timescale to a terminal speed $v_{\rm f}$. Particle acceleration occurs at the reverse shock generated when the outflow collides with a thin, dense shell of slower material released earlier. Accelerated protons are then advected into the shell, where for typical ${ M_{\rm env}, \tau, \text{and } v_{\rm f}}$ they radiate in the calorimetric limit, consistent with correlated optical and $\gamma$-ray emission seen in well-sampled novae. The maximum proton energy, set by a Hillas-like argument, scales with the thickness of the hot post-shock region. Recent work shows turbulent mixing of hot post-shock gas with cooler dense gas may limit this thickness to $\lesssim 10^{-4}$ of the shock radius, explaining low X-ray luminosities. Using this empirically motivated thickness, and assuming efficient magnetic amplification, we predict maximum proton energies $E_{\rm max} \sim 10$ GeV, consistent with $\gamma$-ray spectra of Fermi-detected novae near optical peak ($\sim \tau$). However, as the shock and post-shock layer expand, $E_{\rm max}$ can grow to $\gtrsim 10$ TeV on timescales of a few $\tau$, enabling potential detection by atmospheric Cherenkov telescopes. We encourage TeV follow-up of Fermi-detected novae weeks to months after the optical/GeV peak and quantify the most promising events.

Discoveries of gamma-ray bursts (GRBs) have become commonplace in recent decades, totalling $\mathcal{O}(10^4)$ unique detections across various missions. However, there have been no confirmed discoveries of a gravitationally-lensed GRB, despite expected lensing rates of $\sim1$ in $10^{3}$. In light of this, we complete an archival search for lensed GRBs by cross-matching well-localised \emph{Swift}/XRT-detected bursts with a large all-sky sample of galaxy clusters as potential lenses. We find a total of 17 candidate lensed GRBs defined by a 2 arcminute search radius from a cluster in our sample. 14 of our candidates are either confirmed to be at higher redshifts than their cross-matched cluster, or are consistent with a higher redshift origin based on the Amati relation between $E_{p,i}$ and $E_{\rm iso}$ of GRBs, indicating they are, at some level, lensed by their nearby cluster. Using the Amati relation and the lens-GRB separation, we quantify the magnification experienced by each GRB. We find $\mu < 10$ for all except for one candidate, GRB~071031, which is consistent with $\mu > 10$, but is uncertain. Another candidate, GRB~050509B, does not have a directly measured redshift, but was previously assumed to be at the redshift of its nearby cluster, $z=0.225$. We produce a lens model of this cluster and show that GRB~050509B is consistent with $z>1$ and magnified by $\mu\simeq2-6$. We present these findings in anticipation of future lensed GRB discoveries enabled by facilities such as the Vera C. Rubin Observatory in the coming years.

We present the first implementation of an evolving dust grain size distribution (GSD) within a semi-analytic cosmological model (SAM) of galaxy evolution. This flexible model self-consistently accounts for stellar dust production, shattering, coagulation, accretion of gas-phase metals, and destruction in supernova-driven shocks and hot gas, successfully reproducing key observational constraints. The purpose of this paper is to present the key physical elements of this novel dust implementation in a SAM and to explore controlled numerical experiments to identify the mechanisms shaping the GSD and extinction law in galaxies. Our results show that the GSD evolves from a large-grain-dominated regime at high redshift to a flatter, MRN-like shape at low redshift. This transition occurs earlier for massive galaxies, at a characteristic metallicity determined by the galaxy depletion time. The resulting extinction curves show an increase of the UV/optical slope and a pronounced $2175$ A bump toward lower redshift, in good agreement with the extinction properties of the MW. Through numerical experiments, we find that once stars provide the initial reservoir of large grains, shattering and ISM accretion are the principal mechanisms driving the growth of small grains. When accretion is included, the model robustly reproduces the observed $z \approx 0$ dust masses, largely independent of the specific assumptions adopted for grain-size physics. The extinction properties of MW-like galaxies are also generally recovered, except in extreme cases, such as when grain velocities in turbulent media are assumed to be independent of grain size.

Understanding the multiplicity of stellar systems and the correlations between their hierarchical components provides crucial insights into star formation processes. If binary companions form independently in each component of a wide binary (WB), the fraction of quadruple systems, i.e., 2+2 configurations where both components are themselves close binaries (CBs), should equal the product of individual CB fractions. Using \textit{Gaia} DR3 radial velocity spectroscopy (RVS) data for WB systems, we measure the CB fraction $p$ and quadruple fraction $P_{2+2}$, suggesting an enhancement factor $\kappa = P_{2+2}/p^2 = 2.34_{-0.11}^{+0.12}$, significantly exceeding unity expected under a statistical model of independence. We confirm the significance of this excess by performing two sets of tests: (1) shuffling WB pairings while preserving the overall $\Delta G$ distribution shows no significant enhancement, ruling out selection effects; (2) simulations preserving the spectral type (temperature-dependent) CB fraction also yield the same null excess. When examined as a function of WB separation, the enhancement remains strong at separations $\leq 5\,000$ AU, but shows a decline towards unity at the widest separations ($\geq 10\,000$ AU). An independent proper motion anomaly (PMa) consistency check confirms the enhancement, suggesting a similar value. We further find that the enhancement declines with increasing peculiar velocity, suggesting that dynamical processing in older or dynamically hotter populations may transform 2+2 quadruples into triples over time. Our results provide strong evidence for correlated binary formation processes operating in WB systems.

We present the Dartmouth Stellar Evolution Emulator (DSEE), a flow-based stellar evolution model emulator trained on a comprehensive database comprising over eight million evolutionary tracks that vary across twenty input-physics dimensions and span broad ranges in mass and composition. DSEE learns phase-conditioned stellar state snapshots, unifying track and isochrone construction as marginals of one generative model. It delivers continuous interpolation across high-dimensional physics, probabilistic predictions with calibrated credible intervals, and orders-of-magnitude speedups over direct modeling. Validation against current stellar evolution models shows high fidelity across the HR diagrams, while distributional tests recover the full distributions obtained from brute-force Monte Carlo sampling. To broaden impact, DSEE is integrated into the open-source CONF1DENCE package, enabling fast, end-to-end creation of stellar tracks and isochrones. CONF1DENCE includes the ability to make uncertainty-aware age determinations for clusters taking into account observational effects. CONF1DENCE replaces bespoke, fixed-physics grids with a generative, physics-marginalized emulator, setting a practical new standard for stellar modeling and enabling survey-scale analyses with rigorous uncertainty.

Type Ia supernova (SN) magnitudes present correlations due to the fact that their peculiar velocities are sourced by the large-scale structure of the Universe. This effect can be used to constrain properties related to the distribution and growth of matter perturbations. We analyze both Pantheon+ and Dark Energy Survey (DES-Y5) SN catalogues in combination with CMB data from Planck PR4 to constrain $\sigma_8$ in $\Lambda$CDM, optionally including both curvature and a modified growth index $\gamma$. We show that SN and CMB datasets are highly complementary and capable of measuring $\sigma_8$, $\gamma$ and $\Omega_k$ simultaneously. Using only SN, we find $\sigma_8 = 0.73 \pm 0.22$ ($0.87 \pm 0.31$) for Pantheon+ (DES-Y5) in the base flat $\Lambda$CDM model. Interestingly, allowing for free $\gamma$ and $\Omega_k$, we find hints of positive curvature: $\Omega_k = -0.011 \pm 0.006$ $(-0.014 \pm 0.005)$, which exclude flatness at 2.2$\sigma$ (3.0$\sigma$), for the combination of CMB with Pantheon+ (DES-Y5). Such hints do not degrade if we also include a modified amplitude of CMB lensing, parametrized by $A_L$. We find that $\gamma = 0.519^{+0.061}_{-0.099}$ ($0.461^{+0.085}_{-0.069}$), which are consistent with the predictions of General Relativity. In terms of $f\sigma_8(z)$, we find $f\sigma_8(0.024)=0.461^{+0.066}_{-0.035}$ ($f\sigma_8(0.038) = 0.498^{+0.045}_{-0.050}$) for CMB + Pantheon+ (DES-Y5). Finally, the strong degeneracy between all three $\Omega_k$, $\gamma$ and $H_0$ results in a broader CMB $H_0$ posterior. However, if we include SH0ES $H_0$ data, which is in known strong tension with the CMB in flat $\Lambda$CDM, we find that the $H_0$ tension is recast in terms of a significantly negative curvature and suppressed growth of structures.

Here we investigate a possible presence of a fifth force at the Galactic Center (GC), and its potential influence on the stellar orbits around the central supermassive black hole of our Galaxy. For this purpose we simulated the stellar orbits in a Yukawa gravity model that predicts the emergence of a fifth force, and fitted them into the observed orbit of S2 star around Sgr A* at the GC. The fitting was performed using Markov chain Monte Carlo method which enabled us to constrain the parameters of Yukawa interaction describing the strength $\delta$ and the range $\lambda$ of a fifth force. We studied the following cases for a fifth force range $\lambda$, when it is: i) about a few hundred AU (i.e. deep inside the orbit of S2 star), ii) about a thousand AU (i.e. approximately the size of S2 star orbit), and iii) several thousand AU (i.e. much larger than the size of S2 star orbit). The obtained results showed that as the range $\lambda$ of a fifth force increases, its strength $\delta$ also increases and relative error $\Delta\delta/\delta$ decreases. The resulting fifth-force strengths in all three cases are respectively: $\delta\sim$ 0.005, 0.02 and 0.15. These results are consistent with the corresponding results of both our previous studies and those of other authors, regardless of the different Yukawa-like potentials used to model a fifth force. In addition, we also studied whether the possible small discrepancies from the prediction of General Relativity for the Schwarzschild precession of S2 star could be caused by a fifth force. For this purpose we used the $f_\mathrm{SP}$ parameter that was recently measured in the case of S2 star by GRAVITY Collaboration in 2020. We found that the obtained estimates in all three cases are compatible, within the error intervals, with the measured value of $f_\mathrm{SP} = 1.10\pm 0.19$.

During the 2026 ESO La Silla Observing school, about twenty students attended lectures and performed observations to learn various aspects of observational astronomy. The school, which took place during the first two weeks of February 2026, made use of EFOSC2/NTT and HARPS+NIRPS/3.6m. One of the groups was devoted to the study of binary stars. Several projects were considered and followed up by some of the six students in this group. The first subgroup used HARPS to study the Rossiter-McLaughlin effect in binary stars to infer the relative inclination of the rotation axis of the primary with respect to the orbital plane. A detailed study of the contact binary system HD 115264 led to the conclusion that the primary is well aligned, likely as a result of strong tidal forces within the binary. The second subgroup analysed blue straggler stars (BSS) in open clusters, using both HARPS and EFOSC2. With HARPS, they looked at some well-known long-period binary with the aim of determining their chemical abundances, thereby confirming their membership to the cluster, as well as looking for any chemical anomalies that might be explained by mass transfer. EFOSC2 was used to derive radial velocities of rapidly varying BSS. For one of them - the star Rediet - the students clearly detected and analysed the radial velocity variations due to the second overtone pulsation, thereby confirming its delta Scuti character. Finally, one student used EFOSC2 to study planetary nebulae (PN) - taking nice images of some of these intricate objects, as well as doing time-resolved photometry and spectra of some others. In one case, the binary nature of the central star of the PN was proven, confirming some previous estimates done with ZTF. Each subgroup was thus able to obtain useful research results, which we present hereafter.

Observed galactic cosmic ray intensity can be subjected to a transient decrease. These so-called Forbush decreases are driven by coronal mass ejection induced shockwaves in the heliosphere. By combining in situ measurements by space borne instruments with ground-based cosmic ray observations, we investigate the relationship between solar energetic particle flux, various solar activity indices, and intensity measurements of cosmic rays during such an event. We present cross-correlation study done using proton flux data from the SOHO/ERNE instrument, as well as data collected during some of the strongest Forbush decreases over the last two completed solar cycles by the network of neutron monitor detectors and different solar observatories. We have demonstrated connection between the shape of solar energetic particles fluence spectra and selected coronal mass ejection and Forbush decrease parameters, indicating that power exponents used to model these fluence spectra could be valuable new parameters in similar analysis of mentioned phenomena. They appear to be better predictor variables of Forbush decrease magnitude in interplanetary magnetic field than coronal mass ejection velocities.

Close-binary evolution couples Roche-lobe overflow (RLOF), common-envelope (CE) drag, stellar winds, magnetic braking, and gravitational-wave losses, exchanging mass and angular momentum while reshaping orbits and spins. We present interoperable effects in the REBOUNDx extension to REBOUND that embed these processes within high-accuracy N-body dynamics. The suite includes: a momentum-conserving RLOF operator with conservative and systemic channels and configurable specific-j loss; a CE drag model based on Mach-dependent dynamical friction with kick limiting; isotropic Reimers winds, Parker-type thermal winds, and Eddington-limited outflows powered by a parametric stellar-evolution module supplying mass-dependent R and L; magnetic braking via the Verbunt-Zwaan/Kawaler torque with a saturation-aware closed-form spin update; and post-Newtonian corrections 2PN point-mass and spin-spin; 2.5PN radiation reaction. Linear momentum is conserved for conservative transfer, a minimal corrective torque enforces angular-momentum consistency, and adaptive sub-stepping stabilizes evolution near contact. Inter-module flags coordinate wind/RLOF/CE activity. The unit-agnostic framework enables self-consistent, time-resolved studies of close binaries in isolated or dynamically rich settings. Multiple examples and comparisons against other codes are provided in the Appendix. The code is available at this https URL .

Over the past 30 years, thousands of exoplanets have been discovered, revealing detailed demographics of planets outside the Solar System. One of the most dramatic features of the planet radius distribution is the radius gap, a lack of planets between $\sim$1.8-2 $R_\oplus$. The radius gap is thought to mark the distinction between rocky and gas planets. Recent research has found that the radius gap may not be present in binary star systems. In past studies of planets in binary star systems, the common assumption has been that all of the planets are hosted by the primary star. In many cases, the radius of the planet would be significantly larger if it were orbiting the companion star, which could potentially affect the true radius distribution. It is possible to identify the host stars of planets through stellar density estimates obtained from transit fitting. Using this method, we made probabilistic estimates for the host stars of a sample of 15 transiting exoplanets across 10 binary star systems hosting either 1 or 2 planets, at least one of which would reside in the canonical radius gap if it was circumprimary. We found that 5 of the planets are highly likely to be circumprimary, while the remainder have ambiguous host stars. The lack of unambiguously circumsecondary planets is caused by physical and observational biases that favor circumprimary planets. Nonetheless, the summed posterior probabilities suggest that the canonical radius gap appears less vacant for planets in binaries.

Polarimetric observations at sub-GHz frequencies offer unique access to the magnetized universe through Faraday rotation and depolarization studies, but achieving reliable polarization calibration at these frequencies remains challenging. We report the identification and resolution of a systematic polarization calibration instability in the upgraded Giant Metrewave Radio Telescope (uGMRT) Band 4 (550--750\,MHz). Through diagnostic observations of multiple calibrators, we discovered that the cross-hand phase response varies with the fractional polarization of the observed source, violating the fundamental assumption of calibration transferability in radio interferometry. Systematic engineering tests traced this behaviour to the Quadrature Hybrid (QH) polarizer in the frontend signal chain. We conducted a controlled experiment in which the QH was bypassed in seven antennas, converting them to linear polarization feeds. The bypassed system shows dramatically improved performance: instrumental leakage reduced from 10--15\% to 2--5\%, residual leakage after calibration reduced from $\sim$0.5\% to less than $0.2\%$, and stable cross-hand phases independent of source polarization. For the polarized source DA\,240 (RM\,$=$\,3.3\,rad\,m$^{-2}$), the QH-bypassed system accurately recovers the expected $25^\circ$ polarization angle rotation across the band, which the with QH system fails to reproduce. These results establish that the QH polarizer is the dominant source of polarimetric instability in uGMRT Band\,4 and demonstrate that its removal enables reliable sub-GHz polarimetry. We recommend the linear feed configuration for science cases requiring accurate polarization angle and rotation measure measurements.

Neutron star mergers are multimessenger sources whose dynamics and signals depend critically on neutrinos and their flavor transformations. We investigate whether fast and collisional neutrino flavor instabilities (FFIs and CFIs) arise in a GW170817-like post-merger accretion disk, and how they develop and relax, by performing global and local classical and quantum-kinetic simulations that resolve anisotropies and inhomogeneities in the full six-dimensional phase space. In the accretion disk, the neutrino radiation field naturally develops electron-lepton-number crossings through the interplay between the more isotropic electron neutrino field and the more anisotropic electron antineutrino field. The neutrino field in the disk is also unstable to CFI, although on longer timescales than the FFI. Using local, multi-energy quantum-kinetic calculations at selected points, we find that the growth of unstable modes is well-predicted by a fully anisotropic linear stability analysis and the flavor transformation increases the heavy lepton neutrino fluxes. CFI likewise enhances heavy-flavor fluxes, shows significant impacts from the growth of multi-energy anisotropic modes, and breaks the symmetry of the heavy-flavor sector by raising the average energy of heavy-flavor antineutrinos above that of heavy-flavor neutrinos. However, the CFI remains subdominant to the FFI in most of the disk. In our global quantum-kinetic simulations with an attenuated Hamiltonian, flavor coherence develops primarily in the polar regions. Because the attenuation causes advection to outpace the growth of the instabilities, coherence and flavor conversion remain artificially suppressed within the disk. These results emphasize the resolution and scaling requirements for future global simulations that capture instability growth, saturation, and advection simultaneously.

Ultracompact binary systems, consisting of two compact objects in an orbit $\lesssim 0.5 R_\odot$, should exhibit measurable rates of orbital period change ($\dot{P} \neq 0$) due to the emission of gravitational waves (GWs). Measurements of \pdot\ have so far been limited to the shortest-period ultracompact binaries ($\lesssim 20$\,min). Among the AM\,CVn-type subclass, several works have proposed the presence of extra angular momentum loss beyond GW emission, with magnetic braking being a widely discussed mechanism. If present, this magnetic braking would dominate the angular momentum loss of AM\,CVn-type binaries with orbital periods $\gtrsim 30$\,min. In this work, we present a long-term eclipse timing study of two AM\,CVn-type binaries, YZ\,LMi and Gaia14aae, with respective orbital periods of 28.3\,min and 49.7\,min and continuous observations since 2006 and 2015. Both systems show $\dot{P}$ consistent with zero within $2\sigma$. Their $3\sigma$ upper limits are $1.1 \times 10^{-13}\,{\rm s \, s}^{-1}$ and $9.7 \times 10^{-14}\,{\rm s \, s}^{-1}$ respectively. These non-detections are most simply explained by a scenario in which secular angular momentum loss is not substantially stronger than GW emission at all orbital periods, but is combined with deviations from the secular $\dot{P}$ whose timescales span decades but whose amplitude is $\lesssim 10^{-13}\,{\rm s \, s}^{-1}$. %, orders of magnitude smaller than the eclipse timing variations seen in hydrogen-dominated cataclysmic variables. Our non-detections of $\dot{P}$ represent a limit on the strength of any enhanced angular momentum loss beyond pure GW emission.

Star formation occurs within dense regions of giant molecular clouds (GMCs), however, exactly how gas collects and evolves to form individual stars and what role dense cores play remains unclear. We use the Lagrangian cell information in the STARFORGE simulation suite to track star-forming gas in three GMCs with varying magnetic field strengths. We find that, once a protostar forms, the lifetime of the unaccreted gas correlates with the final stellar mass, where low-mass stars ($M_*$ < 0.5 M$_\odot$) accrete for 0.5-0.6 Myr from a relatively local reservoir of gas, and high-mass stars ($M_*$ > 2 M$_\odot$) accrete over 3.3-4.7 Myr from a much larger volume. Although the protostellar accretion time increases weakly with magnetic field strength, the accreting gas radii, velocity dispersions, virial parameters, and magnetic energy ratios are largely insensitive to the global cloud properties. At the time of protostar formation, the unaccreted gas exhibits linewidth-size and mass-size relations characteristic of turbulently regulated, isothermal dense cores, following $\sigma_v \propto R^{1.0-1.1}$ and $M \propto R^{0.47-0.55}$, respectively. Low- and intermediate-mass stars undergo relatively continuous accretion and their accretion histories are well-fit by either isothermal sphere, turbulent core, or competitive accretion models, where no one model fits all masses. However, many high-mass stars experience intermittent accretion and their accretion histories are not well-fit by any of these models. While the distribution of accreting gas is more extended than typically-defined dense cores, the physical properties and structure of the star-forming gas resemble those of observed cores and are largely regulated by turbulence and feedback.

Magnetars are usually interpreted as highly magnetized neutron stars, yet a small subset of low spin-down sources has motivated alternative scenarios involving highly magnetized white dwarfs. We test whether the observed magnetar sample is consistent with a single neutron-star population or whether the data favor an additional compact-object channel. We combine exploratory machine-learning diagnostics with hierarchical Bayesian population modeling. First, we apply principal component analysis and K-means clustering in $(P,\dot{P},L_X)$ space, and then train a Random Forest classifier with leave-one-out cross-validation to identify the observables driving the empirical split. We subsequently construct a hierarchical Bayesian mixture model that links spin parameters to magnetic-field distributions through covariate-dependent mixing fractions. Posterior inference is performed with Hamiltonian Monte Carlo, and predictive performance is assessed with Pareto-smoothed importance sampling leave-one-out cross-validation. The exploratory analysis reveals a reproducible sub-structure: the Random Forest reaches $>95\%$ LOOCV accuracy, with $L_X$, $\dot{P}$, and $kT$ emerging as the dominant predictors. However, the Bayesian comparison shows no statistically significant preference for a two-population model. Instead, a few low spin-down sources receive intermediate posterior membership probabilities, indicating that they are better interpreted as transitional or outlying objects than as members of a clearly distinct class. Overall, current data do not require a separate white-dwarf magnetar population. The main result is therefore conservative but strong: the observed sample is adequately described by a predominantly neutron-star population, while still allowing physically interesting deviations in specific sources.

Close quasar pairs are rare products of galaxy mergers in which both supermassive black holes (SMBHs) are actively accreting, offering strong constraints on merger-driven active galactic nuclei (AGN) evolution. Identifying close quasar pairs at $z\gtrsim4$ is challenging due to the declining quasar number density in the early Universe. Here we report the confirmation of a close quasar pair at $z=5.7$, J2037--4537, utilizing high-resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations. The quasar host galaxies exhibit tidal disturbed features in both the far-infrared continuum emission and the {\cii} line emission, ruling out the doubly-imaged lensed quasar scenario. The two quasar hosts are massive $(M_\text{dyn}\gtrsim10^{10}M_\odot)$ and star-forming (SFR $\gtrsim500 M_\odot~ \mathrm{yr^{-1}}$). The confirmation of J2037--4537 puts a lower limit on the quasar pair fraction at $5.5<z<6$, $F_\text{pair}>1.2\%$, which is much higher than the quasar pair fraction at $z\lesssim4$. J2037--4537 is expected to form a gravitationally-bound SMBH binary within $\lesssim2$ Gyr. The elevated quasar pair fraction at $z>5.5$, as indicated by J2037--4537, likely contributes to the high gravitational-wave background reported by recent Pulsar Timing Array experiments.

We utilize JWST/MIRI IFU observations from the Galaxy Activity, Torus and Outflow Survey (GATOS) to investigate the diverse range of ionized outflow rates of obscured AGN with similar bolometric luminosity and explore potential associations with AGN feedback. We explore spatial correlations between ionized emission potentially associated with fast shocks ([Fe II]5.34{\mu}m) and the excitation of H2. We further constrain our investigation to the inner 400 pc (the nuclear and circumnuclear regions r < 200 pc), and estimate the excitation temperature and column density of H2 assuming local thermodynamic equilibrium (LTE) and using the S(1) to S(8) rotational H2 emission lines visible to JWST/MIRI spectroscopy. We report the molecular gas temperature of the deprojected 400 pc nuclear region to correlate with the ionized mass outflow rate. We also observe the stronger degree of spatial correlation between [Fe II]5.34um emission and H2 gas temperature. We observe regions of enhanced [Fe II]5.34{\mu}m / [Ar II]6.99{\mu}m spatially coincident with the ionization cones of objects with higher ionized outflow rate and [Fe II]5.34{\mu}m / [Ar II]6.99{\mu}m in the deprojected 400 pc nuclear region to scale positively with both ionized outflow rate and estimated molecular gas temperature. We do not observe the estimated jet cavity power within the central 400 pc to strongly correlate with the ionized mass outflow rate or molecular gas temperature of the nuclear region. We take the preceding observations to suggest a higher degree of interaction between AGN outflows and the circumnuclear disk.

We have evaluated fluctuations in neutrino fluxes caused by solar gravity (g) modes based on the analysis of linear adiabatic oscillation of a spherically symmetric star. We find that the first-order fluctuation is zero due to geometrical cancellation. We still find that the second-order fluctuation is non-zero, which consists of time-varying and non-time-varying components. The amplitude of the time-varying component is small (${\sim} 10^{-9}$ in relative difference, in the case of $\mathrm{^{8}B}$ neutrino) and well below the detection limits of the current neutrino detectors, when we assume the g-mode amplitude parameter $A_{n \ell}$ to be $10^{-5}$, which corresponds to the assumed maximum relative temperature perturbation inside the Sun. Thus, it is at the moment fair to say that detecting individual solar g-modes via the solar neutrino flux measurement is almost impossible. However, the net increase in the mean neutrino flux that originates from the non-time-varying component could be non-negligible. In particular, since $A_{n \ell}$ may be related to convection amplitude, which could change in accordance with the solar magnetic activity, the total net increase in the neutrino flux, which is proportional to $A_{n \ell}^2$, should also change with the solar activity cycle. Such a long-period variation~(${\sim} 11$~years) in the neutrino flux could thus be interpreted as evidence for a bunch of solar g-modes. Comparison of the theoretical prediction with the solar neutrino measurements by, e.g., Super-Kamiokande, may have a potential to put constraints on the theory of the excitation mechanism of solar g-modes.

Long secondary periods (LSPs) in luminous red giants remain the only major class of long-period stellar variability without a secure physical origin. Competing hypotheses include binaries with dusty companions and oscillatory convective dipole modes. We identify the physical and geometric conditions under which oscillatory convective dipole modes produce distinctive harmonic signatures that contrast with those expected from binary systems, and apply this diagnostic to a filtered subset of the OGLE-III LSP sample to identify examples consistent with oscillatory convective dipole modes. We model the geometric flux modulation from oscillatory convective dipole modes and map the range of inclinations, temperature amplitudes, and observing wavelengths for which harmonic features are observable. Using OGLE-III I-band light curves, we require statistically significant power at both sequence D and its harmonic, keeping a filtered sample of 249 stars (2.1\% of the ridge-selected sample). We apply iterative Lomb-Scargle and weighted Fourier decomposition to isolate the fundamental and harmonic components. The relative phase ($\Delta\phi$) between these distinguishes secondary maxima predicted by an inclined dipole from secondary minima caused by eclipsing or ellipsoidal binary this http URL majority of high amplitude stars in the filtered subset show $\Delta\phi$ consistent with secondary minima produced by binary systems. However, a small but statistically non-negligible subset exhibits $\Delta\phi$ consistent with secondary maxima that are difficult to reconcile by eclipsing or ellipsoidal binaries, and instead match the geometric predictions for highly inclined, non-rotating oscillatory convective dipole modes with temperature amplitudes consistent with published models.

We identify shell-like tidal structures in flattened haloes that appear stream-like under different projections. This projection dependence demonstrates how changes in the host halo directly impact the formation and classification of tidal debris, highlighting the challenges of relying solely on visual inspection. To address this, we employ our clustering-based classification framework to systematically categorise tidally disrupted satellites into stream-like and shell-like structures. Our host consists of a static three-component MW model with flattening introduced along the z-axis NFW dark halo. We consider three halo shapes: spherical q = 1, extremely oblate q = 0.5, and prolate q = 1.5. We evolve three subhalo types: a highly radial massive subhalo favouring shell formation, an eccentric orbit leading to stream formation, and an intermediate orbit. We first classify the tidal structures visually using face-on and edge-on density projections of the 3D position distribution. This reveals shell-like and stream-like formations across face-on projections, while edge-on views lead to contrary classifications in some cases. To resolve these ambiguities, we apply the classification method developed in our earlier work, analysing structures in ordered density, radial, and energy-angle space. We further investigate the spatial dispersion of stream-like structures and the rate at which core density reduces as the flattening parameter varies. Our results demonstrate that halo shape variations affect tidal debris formation and classification, as well as the spatial dispersion and core density evolution of streams. These findings offer new insights into the role of dark matter halo geometry in shaping tidal structure formation and its contribution to hierarchical galaxy formation and evolution.

Gamma-ray bursts (GRBs) have traditionally been classified by their prompt emission duration and spectral hardness, with short GRBs (sGRB; $\lesssim2 \ \rm{s}$) originating from compact object mergers and long GRBs (LGRB; $\gtrsim2 \ \rm{s}$) from massive star core-collapse. Recent kilonova (KN) associations with long-duration GRBs have challenged this standard picture. We analyze the host galaxies of nine GRBs with associated kilonova candidates at $z<0.6$, including five sGRB-KNe and four LGRB-KNe. Using both parametric and non-parametric modeling of the host light distributions, we investigate the progenitor environments of these events and test whether their hosts show evidence for recent galaxy interactions that could favor dynamical formation channels or isolated pathways following merger-driven star formation episodes for neutron star binaries. We find that five of the nine hosts display tidal features that show they have likely undergone recent mergers, suggesting that merger-driven, dynamical formation pathways may contribute in some systems. We find no clear morphological distinction between sGRB-KN and LGRB-KN hosts as both populations span a wide range of morphologies, including ellipticals, spirals, and interacting systems with tidal features. Multi-Sérsic modeling of the host light profiles further shows that host-normalized offsets inferred from single-Sérsic fits can be overestimated when the transient is associated with a specific subcomponent of a complex host light profile. These results highlight the importance of decomposing host morphology into physically relevant components when interpreting GRB environments and galactocentric offsets.

We present XRISM/Resolve high-resolution X-ray spectroscopy of the prototypical magnetic cataclysmic variable AM Herculis. All satellite lines of highly ionized Fe are fully resolved. Lighter element lines (Si, S, Ca) show 2 - 3 eV widths consistent with purely thermal broadening, while the broader 6 - 7 eV Fe lines require additional bulk Doppler broadening. Spin-phase-resolved modulations are clearly detected in the Fe XXV and Fe XXVI lines, with semi-amplitudes of $81.8\pm6$ km s$^{-1}$ and $132.5\pm9$ km s$^{-1}$, and mean velocities of $143.6\pm6$ km s$^{-1}$ and $225.6\pm8$ km s$^{-1}$, respectively. After removing these bulk Doppler shifts, we obtain intrinsic Doppler widths of $5.23_{-0.15}^{+0.16}$ eV for Fe XXV and $6.23_{-0.18}^{+0.19}$ eV for Fe XXVI, directly revealing gradients of bulk velocity and temperature in the cooling-flow plasma. We additionally examined the resonance anisotropy predicted by Terada et al. (1999, 2001): the equivalent widths of the Fe XXV and Fe XXVI resonance lines increase at the pole-on phase by factors of 1.30 - 1.35, in positive correlation with their oscillator strengths. Combining XRISM with simultaneous NuSTAR data and PSAC/MCVSPEC plasma models, we derive a self-consistent shock temperature of $24.0\pm0.1$ keV and shock velocity of $1,116\pm2$ km s$^{-1}$. Radiative transfer simulations of the resonance lines further constrain the shock density to about $(5 - 6)\times10^{15}$ cm$^{-3}$, providing a new density diagnostic for accretion columns. The resulting accretion column geometry has a height of 200 - 300 km and a radius of 200 - 400 km.

Stellar obliquity serves as a key diagnostic for tracing the dynamical evolution of bound systems-from giant planets and brown dwarfs to stellar binaries-revealing whether these diverse populations share analogous histories. Here, we report the first obliquity measurement for a double M dwarf system, determined via the Rossiter-McLaughlin effect. The spin axis of the primary star, TOI-5375 ($M_\ast=0.62\pm0.02\,M_\odot$), is well aligned with the orbit of its low-mass stellar companion ($M_c=84.8\pm1.5\, M_J$, $\rm P=1.72\,days$) with a projected obliquity of $\lambda=-13.5_{-13.8}^{+12.4}\,^{\circ}$ and a true 3D obliquity of $\psi=37.5_{-13.4}^{+10.6}\,^{\circ}$. The result indicates that the system either formed with a primordially aligned configuration or has undergone tidal realignment. We further investigate obliquity patterns across giant planets, brown dwarfs and binary stars. It turns out that a few obliquity trends observed in giant planets also tentatively exhibit in the latter two higher-mass populations: 1) well-aligned orbits are preferentially found around cooler host stars ($T_{\rm eff}\leq 6250\,K$); 2) wide-orbit ($a/R_\ast\geq 10$) companions are predominantly aligned; 3) no significant correlation shows up between obliquity and orbital eccentricity in any of the companion classes. By modeling $|\lambda|$ with a two-component Gaussian distribution, we find that the low-$|\lambda|$ components of binary stars and brown dwarfs are more concentrated near zero than giant planets while the high-$|\lambda|$ components of brown dwarfs and binaries remain unclear due to the small sample size.

We present UMI (Unified Median Iterative), a novel robust location estimator for detrending photometric time series in exoplanet transit surveys. UMI modifies the standard Tukey bisquare M-estimator with two innovations: (1) an asymmetric weight function that penalizes downward deviations (transit dips) more aggressively than upward ones, exploiting the physical constraint that transits are always below the stellar continuum, and (2) an upper-RMS scale estimator computed from above-median residuals only, ensuring that transit dips never contaminate the noise estimate. Implemented as a fused HIP/CUDA GPU kernel, UMI achieves 69x faster detrending (3.4 ms vs 234 ms per star) and 37x faster full pipeline throughput compared to the standard wotan's biweight. On 1000 real TESS stars, UMI reduces the median per-star depth recovery error at 0.1% transit depth from 20.5 (biweight) to 15.8%, a 23% improvement. On Kepler, where lower photometric noise allows the asymmetry to work more effectively, the improvement grows to 71% (4.2% vs 14.6%). Validated on 802 confirmed exoplanets across both missions, UMI recovers more planets than biweight, Welsch, and Savitzky-Golay combined. The tool is publicly available as pip install torchflat.

We present an end-to-end, iterative pipeline for efficient identification of strong galaxy--galaxy lensing systems, applied to the Euclid Q1 imaging data. Starting from VIS catalogues, we reject point sources, apply a magnitude cut (I$_E$ $\leq$ 24) on deflectors, and run a pixel-level artefact/noise filter to build 96 $\times$ 96 pix cutouts; VIS+NISP colour composites are constructed with a VIS-anchored luminance scheme that preserves VIS morphology and NISP colour contrast. A VIS-only seed classifier supplies clear positives and typical impostors, from which we curate a morphology-balanced negative set and augment scarce positives. Among the six CNNs studied initially, a modified VGG16 (GlobalAveragePooling + 256/128 dense layers with the last nine layers trainable) performs best; the training set grows from 27 seed lenses (augmented to 1809) plus 2000 negatives to a colour dataset of 30,686 images. After three rounds of iterative fine-tuning, human grading of the top 4000 candidates ranked by the final model yields 441 Grade A/B candidate lensing systems, including 311 overlapping with the existing Q1 strong-lens catalogue, and 130 additional A/B candidates (9 As and 121 Bs) not previously reported. Independently, the model recovers 740 out of 905 (81.8%) candidate Q1 lenses within its top 20,000 predictions, considering off-centred samples. Candidates span I$_E$ $\simeq$ 17--24 AB mag (median 21.3 AB mag) and are redder in Y$_E$--H$_E$ than the parent population, consistent with massive early-type deflectors. Each training iteration required a week for a small team, and the approach easily scales to future Euclid releases; future work will calibrate the selection function via lens injection, extend recall through uncertainty-aware active learning, explore multi-scale or attention-based neural networks with fast post-hoc vetters that incorporate lens models into the classification.

We discover an intrinsic degeneracy in the semi-analytic two-spot model for parameter inference in thermal X-ray pulse-profile modeling. Although this degeneracy exists in our simplified model with two small circular hot spots and without Doppler effects, it causes the likelihood surface of parameter inference based on ST-U and more complex models with Doppler effects to have multi-modal structures. Consequently, the posterior surface may also exhibit multi-modal structures if there is insufficient prior knowledge of parameters. Because of this, the inferred value of the neutron star radius can be biased even by $30\%$. This finding also provides a promising way to explain the multi-modal structures discovered in the evaluation of recovery performance using synthetic pulse-profiles that mimic the PSR J0030+0451 pulse-profiles~\citep{Vinciguerra2023,Vinciguerra2024}. Our work may have profound implications for the reanalysis of NICER data and the analysis of upcoming eXTP data.

Total solar eclipses are not only astronomical spectacles but also great astrophysical laboratories. Their historical records are particularly helpful for assessing the past variability of the Earth's rotation speed. Chinese records played a key role for such analyses. However, Chinese eclipse records from the Míng period have not been used for {\Delta}T reconstructions, partially because most of the contemporaneous eclipse reports are found not in official histories but in local treatises. This study examines eclipse records in the (quasi-)contemporaneous local treatises, concentrating on what explicitly mentioned eclipse totality on the day of a total solar eclipse and what were compiled during the Míng Dynasty. On their basis, our study revised the {\Delta}T constraint in 1361 to -408 s =< {\Delta}T =< 601 s and set new {\Delta}T constraints of 277 s =< {\Delta}T =< 890 s in 1514, -328 s =< {\Delta}T =< 332 s in 1542, and -1762 s =< {\Delta}T =< 1091 s in 1575, respectively. We also revised most of the existing {\Delta}T constraints in the 14th to 16th centuries, using the ephemeris data of the NASA JPL DE 441. Overall, our {\Delta}T constraints generally tighten the {\Delta}T variations more than what M+21 fit for their {\Delta}T spline curve, requiring downward modification and upward modifications for the {\Delta}T reconstructions around 1361 and 1542, respectively. Our results suggest that the {\Delta}T decrease between 1514 and 1567 was slightly steeper than previously considered.

We report results from the simultaneous XRISM (183 ks) and NuSTAR (62 ks) observations of the Seyfert-2 galaxy NGC 4388. This AGN has the brightest Fe K$\alpha$ line among Compton-thin, obscured sources. To model the reflection continuum and fluorescent lines, we employ an updated version of XCLUMPY and a broad line region model with a disk-like geometry. The profile of the neutral Fe-K fluorescent line is well described as the sum of three components convolved with Gaussians with FWHM values of $\sim 290\ \mathrm{km\ s^{-1}}$, $\sim 1470\ \mathrm{km\ s^{-1}}$, and $\sim 11100\ \mathrm{km\ s^{-1}}$. These line widths correspond to radii of 1.5 pc, 0.060 pc, and $1.0\times10^{-3}$ pc by assuming Keplerian motion, which we interpret as the dusty torus, its inner edge region, and the BLR, respectively. The data suggest that the Fe K$\alpha$ BLR component is larger than that of H$\alpha$ (FWHM of 4500 $\mathrm{km\ s^{-1}}$) in the polarized optical spectrum, implying that the velocity field of the BLR is dominated by that parallel to the equatorial plane. In addition, Fe XXVI Ly$\alpha$ and Fe XXV absorption lines are detected, characterized by $\log{\xi} \sim 3.50~\mathrm{erg\ cm\ s^{-1}}$, $\log{N_{\mathrm{H}}} \sim 22.1~\mathrm{cm^{-2}}$, $v_{\mathrm{out}} \sim 40\ \mathrm{km\ s^{-1}}$, and $\sigma_v \sim 160\ \mathrm{km\ s^{-1}}$. We infer that the absorber is gravitationally bound and is possibly associated with a failed wind, consistent with a radiation-driven fountain flow.

While collisionless plasmas are ubiquitously present near astrophysical compact objects, the impact that their composition has on the high-energy emission is presently unknown. We present the first investigation of particle-acceleration mechanisms in kinetic, special-relativistic turbulence, modeling electrons, positrons, and protons with realistic mass ratios. Under global charge neutrality, we introduce a positron fraction and cover regimes ranging from an electron-proton plasma over to pair-dominated plasmas. Using a novel generalized Ohm's law for multispecies relativistic plasmas, we analyze particle acceleration due to electric fields in reconnection events that self-consistently emerge from turbulence. We demonstrate, for the first time, that energization occurs at reconnection current sheets driven by the divergence of the relativistic pressure tensor, which locally aligns with the particle velocity and leads to an efficient energy transfer. The imbalance between electrons and positrons systematically favors electron acceleration, highlighting the necessity of realistic multispecies modeling to capture the nonthermal contributions in accretion flows and relativistic jets from black holes.

Short and long-short gamma-ray bursts (GRBs) are widely believed to be powered by neutron star mergers. In this work, we calculate local rate of such GRBs and find a relatively high value of $\sim 786-2468~{\rm Gpc^{-3}~yr^{-1}}$ when including the very narrow collimation event GRB 061201. Considering that its redshift is not very reliable, after excluding this event, the rate is $\sim 195-666~{\rm Gpc^{-3}~yr^{-1}}$. We also calculate the electromagnetically (EM) bright neutron star merger rate inferred from the LIGO/Virgo/KAGRA observations up to the end of the first epoch of the O4 run, and derive a rate of $\sim 66-347~{\rm Gpc^{-3}~yr^{-1}}$. This rate is somewhat lower than the value obtained from the GRBs, even after excluding GRB 061201. The non-detection of any viable EM bright merger in the O4b and O4c observing runs favors an even lower rate, which starts to challenge the neutron star merger origin of the short and long-short GRBs and may suggest additional contribution from the mergers of other compact object (like the neutron star-white dwarf) binaries, as speculated initially by King et al. (2007) in interpreting the long-short event GRB 060614.

The study of planetary habitability beyond Earth remains a central and challenging project in planetary science. Analysis of large volumes of planetary data from space missions such as CoRoT, Kepler, and JWST is directed ultimately at finding a planet similar to Earth, the Earth's twin, and answering the question of potential exo-habitability. The Earth Similarity Index (ESI) is a first step in this quest, ranging from 1 (Earth) to 0 (totally dissimilar to Earth). To identify planets that may be habitable to the extreme forms of life, we introduce the Mars Similarity Index (MSI). However, extreme forms of life have also been hypothesized under specific conditions in the upper atmosphere of Venus, motivating comparative habitability studies beyond Earth and Mars. The Venus Similarity Index (VSI), introduced here, is defined as the geometric mean of radius, density, escape velocity, and surface temperature, normalized in Venus units (VU). VSI values range from 0 (complete dissimilarity) to 1 (maximum similarity). The VSI provides a comparative framework for identifying Venus-like planetary environments within exoplanet populations. To explore habitability evolution, we further introduce the Ancient Venus Similarity Index (AVSI) and the Future Earth Similarity Index (FESI) to examine early Venusian conditions relative to ancient Earth and to assess potential future evolutionary pathways for Earth-like planets.

Due to its brightness and proximity, the radio galaxy 3C 84 (optical counterpart NGC 1275 in the Perseus cluster) has been the target of extensive studies investigating the central parsec region of its active galactic nucleus. In 2003, its most recent active phase resulted in a plasma ejection visible in the southern jet, which presented a unique opportunity to study jet formation and evolution at high angular resolution with very long baseline interferometry (VLBI). We aim to study the morphology, evolution, and spectral properties of the restarted jet three years after the first ultra-high angular resolution observations with the RadioAstron space-VLBI satellite in September 2013. To study 3C 84, we used space-VLBI observations carried out in September 2016 at 22 GHz with a global VLBI network and the 10 m Spektr-R radio telescope in orbit as well as quasi-simultaneous multifrequency observations at 4.8, 8, 15, and 43 GHz from the Very Long Baseline Array, including the Effelsberg 100 m telescope. We present the 22 GHz RadioAstron image of 3C 84 from 2016, which reveals the source's central region at a 58 microarcsecond effective resolution. During the three years that elapsed between the first and second space-VLBI observations, the source underwent significant morphological changes. We confirm the existence of the limb-brightened jet and counter-jet reported earlier as well as a flip in the position of the hotspot discovered recently via VLBI monitoring at 43 GHz. Based on measuring the collimation profile, we find that it has evolved from being quasi-cylindrical to parabolic. This is most likely the result of the decreased pressure of the mini-cocoon, which was inflated by the jet and contains hot gas that cannot confine the jet efficiently as it propagates further away from the core. Finally, we also constrained the magnetic field strength in the core region and the hotspot.

We present the first observational test of the hybrid ring strategy, a general coordinated signaling scheme proposed by Seto (2025), which provides a practical Schelling-point realization for interstellar signaling. We use the exceptionally bright GRB 221009A as the anchoring flash for the scheme, together with the accurately measured distance to the Galactic center. This combination provides a high-precision relation linking sky position to a tightly constrained arrival-time window. TESS observed the region around the GRB nearly continuously for ~50 days in 2024, providing survey light curves that enable a direct test of this scheme with sharply predicted arrival-time windows of $\sim$3.4 days. Among 58 carefully selected stars, we identify two that show noticeable single-time-bin brightenings inside their predicted windows (where each time bin corresponds to a 200 s integrated TESS exposure). In both cases the brightenings coincide with excursions in at least one nearby star and are therefore most consistent with instrumental origins. This test demonstrates that the hybrid ring strategy is practical with existing survey data and could serve as a promising basis for future technosignature searches.

M dwarfs are the most abundant stars in the Galaxy and exhibit diverse magnetic behaviours. Understanding their large-scale magnetic fields is essential to study stellar dynamos and assess the impact of magnetic activity on planetary environments, yet their magnetic properties and long-term variability remain poorly characterised. We aim to characterise the large-scale magnetic fields of 6 M dwarfs in the southern PLATO field, with rotation periods from ~1 to 17 days and masses between 0.26 and 0.64 Msun. Five stars are partially convective, one fully convective, extending the mass-rotation diagram to previously unsampled regions. We analysed TESS light curves to determine accurate rotation periods and optimise phase coverage for spectropolarimetric observations. SPIRou data were reduced to obtain LSD profiles and longitudinal field measurements, while synthetic spectra fitting yielded small-scale field strengths. ZDI was applied to reconstruct large-scale magnetic topologies. We report a wide diversity of magnetic topologies among the 6 M dwarfs, with 3 main results: (1) Rapidly-rotating (Prot < 2 d) early M dwarfs can generate dipole-dominated fields of moderate intensity, similar to less massive mid-M dwarfs; (2) rapidly-rotating mid-M dwarfs can generate non-axisymmetric large-scale fields with a significant toroidal component; (3) a moderately-rotating (Prot ~ 17 d) early M dwarf shows a surprisingly weak large-scale field. Our findings highlight the diversity of magnetic configurations, including in previously unexplored regions. Long-term monitoring is crucial to distinguish persistent features from variability-driven excursions and to characterise the evolution of surface magnetic fields. Complementary PLATO photometry, including flare and spot-induced variability analyses, will be essential to link surface activity with magnetic properties.

Millisecond magnetars, rapidly rotating neutron stars with ultra-strong magnetic fields, have long been proposed as central engines of gamma-ray bursts (GRBs). For GRBs produced by neutron star mergers, the survival of a long-lived magnetar remnant remains uncertain, as the merger remnant may rapidly collapse into a black hole. In GRB 230307A, multiwavelength observations together with a previously reported 909-Hz periodic signal consistent with millisecond spin in its prompt emission provide strong evidence that such a post-merger magnetar may power the burst. Here we report the discovery of a quasiperiodic modulation with a characteristic period of 4.5 s in the spectral evolution of GRB 230307A, detected consistently across multiple gamma-ray instruments. The modulation is manifested as a coherent, energy-dependent variation of the spectral shape, with the strongest signature in the evolution of the peak energy. Within the magnetar-engine framework, such a low-frequency modulation can be interpreted as a manifestation of large-scale periodic variations associated with the central engine. If interpreted in terms of free precession, the observed timescale implies a stellar ellipticity of $\epsilon \gtrsim 2.4 \times 10^{-4}$, corresponding to an internal magnetic field strength of $B_t \gtrsim 1.6 \times 10^{16}$ G, alongside a dipole field of $B_p \approx 5.6 \times 10^{15}$ G inferred from the early X-ray emission. These results suggest that such systems may provide potential sources of post-merger gravitational waves (GWs), motivating targeted searches following GRB triggers.

Forbush decreases (FDs) are short-term reductions in galactic cosmic ray flux caused by interplanetary disturbances. During some interplanetary coronal mass ejection (ICME) events, neutron monitor (NM) data also contain variations produced by geomagnetic storms. Earlier studies emphasized apparent effects near 10~GV, but storm-time changes in geomagnetic cutoff rigidity can either increase or decrease the ground-level count rate. Using a recently published hourly proton flux reconstructed from NM data for May 2011 through October 2019, the interval covered by the published AMS daily proton fluxes, we show that these localized anomalies can extend to lower rigidities and reach 1~GV in some events. Such effects can bias the rigidity dependence inferred from NM-based hourly proton spectra during disturbed intervals. Because AMS measures proton rigidity directly in space, its daily proton spectrum is not affected by cutoff variations at ground stations and provides a stable reference. We therefore use AMS to constrain corrections for selected events. The correction removes localized anomalies while preserving the broader FD evolution, and for a representative ICME event it brings the corrected daily averages closer to the AMS measurements. Our results show that short-timescale cosmic ray variability during FDs reflects both heliospheric modulation and storm-time changes in geomagnetic shielding.

The discovery of the phase space spirals in the Solar neighborhood in Gaia Data Release 2 has prompted various attempts to understand their origin. A source of bending waves, which has been neglected as a cause of the phase spiral, is irregular gas inflow along the warp. We aim to study whether perturbations by the gas warp could induce phase spirals. Accounting for this additional formation scenario for phase spirals could improve our current understanding of the perturbation history of the Milky Way disc. We use two N-body + SPH (Smooth Particle Hydrodynamics) simulations of an isolated galaxy to search for, and study, warp-induced phase spirals. We study the emergence and propagation of the detected phase spirals using Fourier decomposition. We detect strong one-armed phase spirals in the warped simulation. These phase spirals are prevalent and persist over ~10 Gyr. The morphology of these phase spirals varies with location and evolves with time. In particular, the emergence rate of the phase spiral evolves with the gas inflow at the outer disc and the bending wave amplitude, indicating that these phase spirals are a record of warp-induced bending waves. We find that these phase spirals can reach amplitudes comparable to those in the Gaia DR3. We only detect weak and stochastically distributed phase spirals in an unwarped control simulation. We conclude that phase spirals can be induced by the irregular gas accretion along the warp. These phase spirals occur globally and are long-lived.

We present a low-redshift-agnostic compression of anisotropic baryon acoustic oscillation (BAO) distances to constrain the normalized dark-energy density evolution, $X(z)\equiv \rho_{\rm DE}(z)/\rho_{\rm DE}(0)$, above the lowest BAO redshift node $z_1$. Standard BAO summaries include the transverse comoving distance $D_{\rm M}/r_{\rm d}$, which depends on the integral of $H^{-1}(z)$ from $z=0$ to $z$ and therefore mixes the expansion history at $z<z_1$ with the higher-redshift signal. We instead replace the set $D_{\rm M}(z_i)/r_{\rm d}$ by adjacent increments $\Delta D_{\rm M}(z_i,z_{i+1})/r_{\rm d}$ while retaining the radial distances $D_{\rm H}(z_i)/r_{\rm d}$. The mapping is linear, so the covariance propagates exactly. This compression intentionally removes one absolute transverse-distance mode, namely the additive contribution to $D_{\rm M}/r_{\rm d}$ below the first BAO node, and preserves the remaining information relevant to reconstructing the expansion history above $z_1$. Applied to DESI DR1 and DR2 anisotropic BAO measurements, the method yields almost uncorrelated constraints on piecewise-constant interval parameters $X_j$. In this sense, the compressed likelihood provides a conservative band-power-like estimate of dark-energy evolution: each interval is constrained mainly by BAO information from its own redshift range, while one nonlocal transverse mode and stronger global assumptions are deliberately projected out or marginalized over. Because our baseline analysis also marginalizes over bin-local matter-density and distance-scale parameters with broad external priors, the resulting $X_j$ constraints should be interpreted as a low-redshift-agnostic BAO baseline rather than as a fully prior-free reconstruction. All bins are consistent with $X=1$ within current uncertainties.

The orbital decay of binary systems is a critical process for understanding the evolution of massive binary black holes (MBBHs) and binary star formation. Performing high-resolution three-dimensional magnetohydrodynamic (MHD) simulations, we investigate a binary system that accretes gas from an infalling envelope analogous to the collapse of molecular clouds in the context of binary star formation. Our simulations reveal the presence of outflows/jets launched from both the circumstellar (mini) disks and the circumbinary disk (CBD). The magneto-rotational instability is also excited within the CBD. These magnetic processes efficiently extract orbital angular momentum from the binary and thus drive orbital decay, while a purely hydrodynamical model exhibits orbital expansion. The decay rate reaches $\sim 0.3-0.7\%$ per orbital period, depending on the initial magnetic field strength. By appropriately scaling these numerical results, we propose a new mechanism for MBBHs mergers within a Hubble time, overcoming the bottlenecks encountered at separations near the final parsec scales. Additionally, we present a formation scenario for close twin binary star systems, emphasizing the significant role of magnetic processes in driving orbital evolution across various astrophysical systems.

The approximately flat outer parts of spiral galaxy rotation curves are commonly interpreted as evidence for a discrepancy between the observed baryonic mass and the dynamical mass inferred from the measured orbital velocities. In most standard analyses, this discrepancy is quantified using $v^2(R)=GM(<R)/R$, which is exact only under spherical symmetry. However, spiral galaxies are flattened disk systems, for which mass exterior to the galactocentric radius under consideration can contribute non-negligibly to the gravitational field. We introduce the Lost and Found (LF) model, a geometrically consistent Newtonian framework based on direct full-disk gravitational integration and a parametrized representation of the disk surface density. In this approach, the gravitational field is computed without imposing spherical symmetry, and the disk mass distribution is represented by two exponential components with a smooth outer truncation. We apply the LF model to a heterogeneous sample of disk galaxies spanning a broad range of masses and radial extents. The model reproduces the main observed features of the rotation curves, including the inner rise and the approximately flat outer behavior, without explicitly invoking a dark matter halo or modifying Newtonian gravity. Across the sample, the LF-inferred mass scales nearly linearly with the conventional dynamical mass, with a characteristic reduction factor $\eta_{LF}$ ~ 0.67. These results indicate that part of the inferred mass discrepancy may arise from the geometric treatment of gravitation in disk galaxies, and motivate a reassessment of mass inference in non-spherical systems.

Blazars, the jet dominated class of AGN comprising flat spectrum radio quasars (FSRQs) and BL Lac objects (BL Lacs) are now increasingly identified as potential sources of high energy neutrinos. Such neutrino blazars are ideal targets to investigate the high energy emission processes and to understand their role as neutrino sources. We report results on four candidate neutrino blazars, PKS 0446+112, TXS 0506+056, PKS 1424$-$418 and PKS 1502+106. We carried out $\gamma$-ray spectral and timing analysis on three time periods that comprise a quiescent epoch, an epoch that corresponds to neutrino detection and a flaring epoch. We also carried out modeling of the broadband pectral energy distribution (SED) on those three epochs. We found that the $\gamma$-ray spectra of the BL Lac TXS 0506+056 can be adequately described by a power-law, while the spectra of the other three FSRQs require a log-parabola model. On shorter timescales, we observed flux variability with doubling/halving timescales of 4.70 hrs, 9.24 hrs, 30.76 hrs and 15.42 hrs for PKS 0446+112, TXS 0506+056, PKS 1424$-$418 and PKS 1502+106, respectively. The SEDs of most of the epochs for the sources are well explained by a leptonic scenario. However, the quiescent epoch of PKS 1502+106 and the neutrino-emission epoch of PKS 0446+112 required an additional hadronic component to reproduce the observed SEDs. Our analysis reveals a complex interplay of leptonic and hadronic processes. While certain neutrino-associated epochs align with a leptonic model, others necessitate a hadronic component to explain the emission features.

We analyze three NuSTAR observations and two NICER observations of the transient accreting millisecond pulsar SAX J1808.4-3658 in the hard spectral state during its most recent outbursts in 2022 and 2025. The spectral analysis of the persistent emission shows that the continuum is well described by an absorbed thermal Comptonization model with a high plasma temperature of ~25-90 keV. A prominent iron emission line around 5-8 keV and a Compton hump around 15-30 keV have been detected from all NuSTAR observations, indicating the reflection of the hard X-ray photon from the accretion disk. We employ the relativistic reflection model relxillCP to describe the reflection phenomena. The spectral fit of three NuSTAR observations shows that the inner disk radius moves outward, the Comptonized thermal emission decreases in flux, the mass accretion rate decreases, and the disk becomes less ionized as we proceed from the 2022 to the 2025 observations. Reflection studies also reveal a moderate inclination of the source within ~30-50 degrees. During the 2025 September observation, the inner radius of the disk is significantly truncated (~23R_g), and the corresponding magnetospheric radius is comprehensively larger than the disk's co-rotation radius, suggesting a hint of the transition to the propeller regime. Although the disk is truncated at the larger radius, accreted material is still reaching the surface of the neutron star, which is confirmed through the detection of a Type-I X-ray burst during this NuSTAR observation. The spectral analysis of the burst suggests helium burning at a low ignition depth.

We perform a comprehensive observational test of a canonical quintessence model driven by an exponential potential, motivated by its emergence in higher-dimensional theories, string-inspired scenarios, and modified gravity. Using a Markov Chain Monte Carlo framework, we constrain the model with the latest high-precision observational datasets including Cosmic Chronometers, Baryon Acoustic Oscillation, Pantheon+, and DES-SN5YR Type Ia Supernovae. The combined data significantly tighten the parameter bounds on (H0, Omega_m0, eta0, gamma) and yield predictions for the Hubble parameter H(z), the distance modulus mu(z), and the scaled comoving angular diameter distance that remain in excellent agreement with observations and closely follow the LCDM baseline. An information-theoretic model comparison using the Akaike Information Criterion shows that the exponential quintessence model remains statistically comparable with LCDM despite having additional parameters. The model successfully reproduces the transition from matter domination to late-time acceleration, maintains w_tot > -1, and provides an age of the universe consistent with Planck 2018. Statefinder diagnostics indicate trajectories approaching the LCDM fixed point with small deviations, and energy condition analysis confirms physical viability, with only the Strong Energy Condition violated at late times as required for acceleration.

Interactions between stellar and planetary magnetic fields are expected to produce observable radio and optical signals modulated by their orbital periods, but direct detections remain elusive. We analyze 17 years of spectroscopic data of the GJ 436 system. This M2.5 V star hosts a transiting Neptune-sized planet in a close-in, inclined orbit. The data shows repeated enhancements of the stellar chromospheric activity at approximately the same phase of its 8-year activity cycle modulated by a combination of the planet's orbital period and the stellar rotation. We interpret this modulation as star-planet interaction. We propose a new geometrical model to interpret these signals, then, estimate the power of the interaction and, from models, estimate the magnetic field of GJ 436 b to be between 6 and 110G. This finding opens new pathways to detect star-planet interactions and to investigate planetary magnetic fields and their implications on atmospheric retention and detectability.

We examine the relationship between the mass of present-day central supermassive black holes (SMBHs, $M_{\rm BH}$), and the stellar mass ($M_{\star}$) and halo mass ($M_{200}$) of their host galaxies in the EAGLE simulation, and find that scatter about these relations correlates with both halo structure and galaxy morphology. EAGLE reproduces the observed $M_{\rm BH}$-$M_{\star}$ relation, including (qualitatively) its dependence on morphology: at fixed $M_{\star}$, disc-dominated galaxies host less massive SMBHs than ellipticals. We show that $M_{\rm BH}$ correlates with $M_{200}$, as expected if SMBHs are regulated by processes acting on the scale of the host dark matter halo, but exhibits a tighter correlation with the halo binding energy ($E_{\rm bind}$), signalling that this quantity, which encodes information about both halo mass and halo structure, is more fundamental to $M_{\rm BH}$. As with $M_{\rm BH}$-$M_{\star}$, scatter about the $M_{\rm BH}$-$E_{\rm bind}$ relation is strongly correlated with morphology. Gas in the central few parsecs of galaxies with present-day discs retains strong rotational support as the galaxy grows, inhibiting inward transport and precluding periods of rapid SMBH growth by gas accretion. In galaxies destined to be present-day ellipticals, however, this rotational support is disrupted, enabling gas to be funnelled onto the central SMBH, triggering rapid growth. Evolution of the mass fraction of stars formed ex-situ indicates that this disruption is caused by galaxy-galaxy interactions and mergers. Our findings corroborate the conclusion of recent studies, based on controlled simulations of an ~$L^{\star}$ galaxy, that prolonged secular galaxy evolution inhibits central SMBH growth.

Extinction maps are essential for tracing interstellar dust and enabling accurate stellar population studies in galaxies. Here, a high-resolution extinction distribution of nearby galaxy M33 is constructed by fitting multiband color indexes of the individually resolved red giant branch (RGB) stars from the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) survey. Achieving an angular resolution of approximately 6$^{\prime\prime}$ ($\sim$ 24.4 pc), the extinction map reveals the intricate and heterogeneous distribution of dust throughout the entire disk of M33, with distinct delineation of spiral arms, inter-arm regions, and compact dust clouds. In addition, it exhibits strong spatial correspondence with the distributions of total hydrogen, H I, and CO, underscoring the reliability of the extinction map for tracing both diffuse and dense components of the interstellar medium. The derived $V$-band extinction reaches up to 2.5 mag per pixel, with a mean value of about 1.05 mag. Beyond providing new insights into the dust structure of M33, the extinction map offers a robust foundation for accurate extinction corrections and will support future studies, including upcoming observations with the Chinese Space Station Telescope.

Dust in the diffuse interstellar medium remains incompletely understood with regard to the structure, composition, size distribution, and alignment properties of the grains. Joint observations of reddening, starlight polarisation spectra, and polarised dust emission for individual sightlines provide essential constraints on such properties. We study a far-UV selected sample of 96 reddening curves, for which optical linear polarisation spectra were obtained with FORS at the VLT as part of the Large Interstellar Polarisation Survey (LIPS). Starlight polarisation spectra for 60 stars are presented in this work. These data are combined with Gaia distance estimates and Planck thermal dust emission. A three-component dust model is made publicly available. It consists of nanoparticles, amorphous grains, and micrometre-sized dust agglomerates, varying axial ratios, porosities, sizes, element abundances, and alignment efficiencies to match the observations. The diversity of reddening and polarisation spectra is well reproduced by prolate grains with typical axial ratios of two, a porosity of 10%, and high alignment efficiencies. Such efficiencies can be achieved with radiative torque alignment theory (RAT), but not with imperfect Davis-Greenstein (IDG) alignment, except when adjusting the magnetic-field orientation to maximise the polarisation. Micrometre-sized dust contributes wavelength-independent grey extinction in the optical, accounts for about one-third of the visual extinction, and carries one-third of the dust mass. A follow-up submillimetre survey with high-resolution polarimetry will further constrain grain shapes and alignment physics.

Dust torus plays the key role in determining active galactic nuclei (AGN) observational appearance. Here, the scenario of accretion regulated central dust torus is tested for the first time in the individual changing-look AGN (CLAGN) SDSS J1011+5442. Through the dependence of broad H$\alpha$ luminosity on continuum luminosity, the scenario of moving dust clouds can be ruled out in SDSS J1011+5442. Meanwhile, virial BH mass in the bright state is consistent with the M-sigma relation determined mass, indicating the virialization assumptions efficient in central BLRs. However, the virial BH mass determined in the dim state is 60 times smaller than the M-sigma relation determined value. The contrary properties of broad H$\alpha$ in different states can be naturally explained by the scenario of accretion regulated dust torus. Below a critical Eddington ratio, opening angle of dust torus declines with increasing accretion rate, leading to only outer part of central BLRs for broad H$\alpha$ with smaller line widths detected in the dim state but all the BLRs detected in the bright state. The results in this manuscript not only indicate properties of central dust torus having apparent effects on variability properties of CLAGN, but also indicate that studying CLAGN could provide further clues to check dynamical evolving models for dust torus in AGN.

Debris disks are created and sculpted by planetary bodies in the orbital space they share. The properties of these disks, including mass, orbital extent, and morphology, can be indicators of their planetary shepherds. Recently, T. Pearce and collaborators placed limits on the masses and orbits of hypothetical planets around 178 stars with resolved debris disks. We consider 176 of these stars, all the objects that have astrometric data in the Gaia Data Release 3 archive, to assess planet detection from astrometry. Our analysis begins with a set of stellar hosts of known exoplanets, selected to roughly match the parallax, apparent magnitude, and color of the 176 debris disk systems. We confirm that Gaia's ruwe parameter, a measure of the quality of astrometric fitting to a linear drift model, is sensitive to the presence of massive companions, even planetary ones. Guided by ruwe and a metric derived from a machine-learning algorithm trained on Gaia parameters from the exoplanetary host data set, we identify promising stars with debris disks that may host as-yet-undiscovered planets. These stars will be compelling subjects for time-series analyses with Gaia Data Release 4.

We present a JWST/NIRSpec-IFU study of metallicity gradients in seven low-metallicity systems at $z=7.2-9.5$. The main sample spans stellar masses of $\rm \log(M_*/M_{\odot}) \sim 7.8-9.5$, star formation rates (SFRs) of $\rm \log(\text{SFR} / M_{\odot} \text{yr}^{-1}) \sim 0.5-2.5$, and gas-phase metallicities of $4\%-15 \%~Z_\odot$. Within our sample, we also identify three low-metallicity satellite galaxies associated with two of our sources, providing a rare view of early-epoch interactions. The three satellites exhibit even more primordial properties, with metallicity $3\% -4\% ~Z_\odot$ and low star-formation activity ($\rm \log(\text{SFR} / M_{\odot} \text{yr}^{-1}) \sim -0.5$ to $-0.9$). We find that our galaxies, and especially the satellites, are significantly offset from the local Fundamental Metallicity Relation (FMR), with deviations reaching $\Delta \text{FMR} \approx -0.9$ dex. This indicates that these galaxies are likely experiencing strong accretion of pristine gas. Overall, we observe a large scatter in radial metallicity gradients, ranging from positive to negative with an average metallicity gradient of $\rm -0.02 \pm 0.04 \ dex \ kpc^{-1}$. Flat gradients are found in systems with confirmed satellites, suggesting that tidal interactions and mergers drive the radial mixing necessary to homogenise the interstellar medium. The (tentative) presence of an AGN in two of our sources suggests that strong feedback may also be responsible for the observed flat gradients. Conversely, the detection of a positive gradient in one source points toward a direct funnelling of metal-poor gas inflow into the central region of the galaxy. These results show that galaxies in the first billion years grow through diverse, episodic processes, suggesting that early evolution is characterised by structural variety rather than a single, predictable path.

Understanding the formation and evolution of late type galaxies (LTG) requires deep imaging for tracing the faintest stellar components in their outskirts. Despite their crucial role in the buildup of stellar mass, these low surface brightness (LSB) features remain largely unexplored due to observational limitations. The VST-SMASH is designed to fill this gap, providing deep, wide field optical imaging for a volume limited sample of nearby LTG, overlapping with the Euclid Wide Survey in the South. This paper aims to introduce the VST-SMASH survey and showcase its scientific potential through the analysis of IC 5332, a LTG observed in the g, r, and i bands. The main goal is to demonstrate the depth, quality, and diagnostic power of the dataset in tracing LSB features and structural components in galactic outskirts. We carried out detailed surface photometry of IC 5332 to extract radial surface brightness and color profiles down to LSB regime. We performed multicomponent Sersic decompositions and constructed stellar mass surface density profiles. We identified and characterized faint stellar streams, estimating their colors and comparing them with adjacent galactic regions. While the internal (1Reff) negative colour gradients can be explained by dissipative collapses and SN outflows, the color profiles at larger radii reveal a significant gradient toward redder colors, consistent with the presence of accreted populations in the outskirts. We also find bluer r - i, which could be explained by strong Ha emission. These findings support a scenario of ongoing stellar mass assembly through accretion and highlight the capability of VST-SMASH to uncover faint structures in nearby galaxies.

We present LightCurveLynx, a flexible and extensible software framework for end-to-end forward modeling time-domain light curves. Given the growing need for realistic simulations in the time-domain astronomy community, LightCurveLynx is designed to support a wide range of applications, including the development and validation of analysis pipelines, the optimization of survey strategies, and simulation-based inference studies. Realistic simulations can be generated from real survey metadata, forecasted survey plans, or user-defined mock survey strategies. We demonstrate the functionality of LightCurveLynx by generating a realistic simulation of Type Ia supernovae that is representative of the ZTF SN Ia Data Release 2 dataset and perform extensive comparisons between the simulated and observed samples to validate the software. The simulation shows excellent agreement with the data in parameter distributions (with the Kullback-Leibler divergence values around 0.01-0.02) and in noise properties. The Hubble diagram generated from the simulation also indicates that the sample is complete up to redshift 0.06, which is consistent with previous studies. Our results confirm that LightCurveLynx is robust, accurate, and ready for community use and contribution.

Little Red Dots (LRDs) have emerged as a key population linked to early black hole growth, yet photometric selections have predominantly targeted only the most extreme red systems, thereby shaping our current understanding of this new population of objects. In this work, we deliberately explore a broad range of optical redness while enforcing stringent compactness and visual inspection to ensure robustness and minimize contamination. Leveraging the depth and multiwavelength coverage of the JWST Advanced Deep Extragalactic Survey (JADES) data in the GOODS-North and GOODS-South fields, we construct the largest photometric census of LRDs to date in these fields, comprising 412 sources over $z\approx2\text{--}11$ across $\approx349.6$ arcmin$^2$. We show that classic extreme color cuts isolate only a minor fraction of this population ($\lesssim25\%$), while the majority of LRDs span a broader, largely unexplored parameter space. We quantify how selection strategies impact UV and optical luminosity functions and number density evolution, finding that current demographic trends of LRDs are strongly driven by selection biases and further limited by incomplete identification at both high and low redshift. Spectroscopically confirmed LRDs reveal a continuous range of spectral shapes, consistent with varying Active Galactic Nucleus (AGN) and host contributions in agreement with recent findings. Our results demonstrate that commonly adopted, purity-driven selections bias current demographic constraints toward the most extreme systems, potentially misrepresenting the diversity and evolution of the LRD population. Accounting for these selection effects is essential for interpreting LRDs and their role in early black hole growth.

We use a high signal-to-noise MIRI-LRS spectrum (5.45 - 11 $\mu$m, R$_\lambda$ $\sim100$) of COCONUTS-2~b revealing prominent molecular features of H$_2$O, CH$_4$ and NH$_3$. This dataset is combined with spectra from Gemini/FLAMINGOS-2 and JWST/NIRSpec (G395H), as well as photometry from WISE and Spitzer, resulting in almost continuous wavelength coverage from 1 to 15 $\mu$m. We analyze the data using five grids of self-consistent atmospheric models, spanning a wide range of T$_\text{eff}$, log(g), and [M/H]. We also investigate the use of Gaussian Processes to account for correlated noise either caused by the spectrograph or by systematic departures of models in the inversion framework. All models manage to fit the overall combined observations but predict fainter flux in Y- and N-bands. Classical model comparison suggests that the ATMO2020++ synthetic specra (with and without PH$_3$) are statistically preferred. Fitting for the correlated noise of the three spectroscopic instruments, ATMO2020++ models yields constraints consistent with previous studies and evolutionary models predictions: T$_\text{eff}$ $=496^{+5}_{-3}$ K, log(g) $=4.30^{+0.04}_{-0.02}$ dex, [M/H] $=-0.02^{+0.03}_{-0.02}$ dex, and R $=1.03^{+0.01}_{-0.02}$ R$_\text{jup}$. The extended wavelength coverage provided by MIRI (accounting for 41% of the bolometric flux) completes the SED, yielding a precise luminosity estimation of log(L/L$_{\odot}$) $=-6.166\pm0.002$ dex. Combined with a previous estimate of the system age ($414\pm23$ Myr), cooling models predict a mass of M $=7.3\pm0.3$ M$_\text{jup}$.

We investigate observational tests that can be used to distinguish between broad classes of cosmological models. This is achieved using curvature-consistency tests of the Friedmann-Lemaître-Robertson-Walker (FLRW) models, which we investigate in two scenarios where they can be violated; (i) when the optical properties of the cosmology deviate from the expectations of FLRW, and also (ii) when the large-scale expansion of the cosmology is different from FLRW. We identify useful ways to determine the properties of these alternative scenarios in terms of the violation of the curvature-consistency tests, and propose a new null test that can be used to isolate cosmologies with non-FLRW observational relations. The characteristic signatures we find can be used, together with the results of recent and upcoming cosmological observations, to probe and/or rule out large classes of cosmological models. This becomes an increasingly important task as the number of proposals in the literature increases, as cosmologists attempt to explain tensions, anomalies, and the dark sector of the Universe. Our approach provides a clear route for telling apart these different proposals, and offers a new opportunity for using precision cosmological data to efficiently discriminate between cosmological models.

HATS-75 b is one of the recently discovered Giant Exoplanets orbiting M-dwarf Stars (GEMS) with a transmission spectrum shaped by both its atmosphere and the active stellar surface it transits. As part of a JWST program studying 7 GEMS, we observed three transits of HATS-75 b with the NIRSpec PRISM instrument (0.6-5.3 um). The planet's spectra exhibit a slightly larger transit depth at shorter wavelengths, indicative of hazes or stellar contamination due to stellar heterogeneities outside the transit chord, i.e., the transit light source (TLS) effect. While both a hazy atmospheric model or TLS model can replicate the transmission spectrum, independent evidence (.e.g, stellar rotation, spot-crossing events) favors a model that includes contamination from unocculted starspots and faculae. Within this stellar heterogeneity / TLS-based framework, atmospheric retrievals yield remarkably low atmospheric metallicity (log[M/H]=-1.74^{+0.92}_{-0.76}) and super-solar carbon-to-oxygen (C/O=1.04^{+0.40}_{-0.09}), which paired with a best-fit interior model with bulk metallicity of Z_p=0.20+/-0.04, implies poor vertical mixing within the planet. Retrievals also detect robust absorption signatures of CH4, CO, and CO2. We obtain only an upper limit for H2O, consistent with its atmospheric spectral features being masked by stellar contamination. These results underscore the importance of accounting for stellar heterogeneity when interpreting exoplanet transmission spectra and highlight HATS-75 b as a significant asset to our understanding of giant exoplanets around M-dwarfs with JWST.

The formation of late-type galaxies has traditionally been described via two pathways: one producing pressure-supported classical bulges, the other rotationally supported pseudo-bulges. Early studies relied on photometric decompositions assuming an exponential disk extrapolated inwards. Recent high-resolution observations, however, reveal a far more complex landscape in disk galaxy centres. We investigated the morphology of central stellar components in intermediate-to-massive spiral galaxies, focusing on disentangling cold, warm, and hot orbital contributions, critically reassessing the standard approach of extrapolating the exponential disk profile inwards. We developed GLANCE (Galactic archaeoLogy via chronochemicAl and dyNamiCal modElling), a tool for photometric, chronochemical, and dynamical galaxy analysis, applied to 8 high-resolution MUSE galaxies to derive stellar population properties and decompose orbits into cold, warm, hot, and counter-rotating (CR) components. We uncovered remarkable structural diversity in the dynamically cold central component: one galaxy displays an exponential profile throughout, while the majority exhibit either a pronounced central drop resembling a doughnut-shaped structure or a compact inner disk significantly steeper than the outer disk. Most galaxies hosting nuclear disks are classifiable as classical bulges - hot, old, red, high bulge-to-total ratio - contrasting with galaxies showing a central cold-component deficit. Beyond the bulge, cold plus warm orbit contributions remain below the total, indicating non-negligible hot or CR orbits with Sérsic indexes consistently above unity. These results highlight the composite nature of disk galaxy centres and the need for decomposition methods that avoid extrapolating the outer disk inwards, requiring large IFS samples across a broad mass range, complemented by simulations such as IllustrisTNG50.

Multi-dimensional unified model atmospheres with winds of massive stars have so far been studied under the assumption of equal flux, Planck, and energy weighted mean opacities, which effectively means these models have been in local thermodynamic equilibrium (LTE). Although LTE may be a valid approximation in deeper atmospheric layers, it breaks down in the extended outflowing parts. As such, the opacities governing the heating and cooling of the gas are neither the same nor equal to flux-mean opacity in those regions. We present an approximate NLTE procedure that accounts for scattering in the computation of energy and Planck-mean opacity from a multitude of spectral lines in an accelerating medium. The formalism evaluates the opacities using Sobolev escape probabilities and effective thermalization parameters from a line database consisting of ~4 million spectral lines. RHD simulations are calculated as before with a hybrid opacity scheme combining Rosseland means with line opacities in an accelerating medium. Due to their high velocity dispersion, upon interaction, they produce localized shock fronts with the gas temperature exceeding the photon temperature. Due to improved treatment of heating and cooling in outflowing parts, the radiation and gas temperatures in the wind are no longer the same, as was the case in previous multi-dimensional simulations. Instead, gas gets heated at shock fronts, but due to strong radiative cooling remains localized. The net result is a multi-component wind structure not only in density and velocity, but also in temperature. This likely has important consequences for the formation and interpretation of observed O-type star wind spectra.

The Lorentz force exerted by a primordial magnetic field (PMF) on the coupled baryon-dark matter system may enhance total matter power at small scales after recombination. In the post-reionization (post-EoR) era, a weakly scale-dependent PMF of sub-nG strength is thus expected to influence the Lyman-$\alpha$ (Ly$\alpha$) power spectrum, the 21 cm power spectrum, and the Ly$\alpha$-21 cm cross-spectrum at scales $k\gtrsim 1\:h/\textrm{Mpc}$. We investigate the prospects of constraining the PMF sector via these three cosmological observables, by employing SNR estimation and Fisher forecast on the PMF amplitude $B_0$ and spectral index $n_{\rm B}$, for a next-generation DESI-like spectroscopic survey and two upcoming 21 cm facilities, namely SKA1-Mid and PUMA. Our results indicate the possibility of constraining both PMF parameters with $\lesssim10\%$ relative errors through the uncontaminated 21 cm auto-spectrum as well as the Ly$\alpha$-21 cm cross-spectrum probed with the DESI-like+SKA1-Mid combination. Indicatively, the Ly$\alpha$-21 cm cross-correlation via DESI-like+SKA1-Mid is predicted to constrain a fiducial scenario $B_0=0.8$ nG and $n_{\rm B}=-2.9$ with $1\sigma$ errors $\Delta B_0\approx 0.07$ nG and $\Delta n_{\rm B}\approx0.02$. The DESI-like+PUMA setup is predicted to fare relatively worse due to its restriction to larger scales, resulting in comparatively one order of magnitude relaxed error bounds for similar fiducials. Since the Ly$\alpha$-21 cm cross-signal is expected to be largely insensitive to foreground contamination (unlike the 21 cm auto-spectrum), it may serve as an optimal foreground-immune post-EoR probe to constrain a weakly scale-dependent sub-nG PMF via future DESI-like+SKA1-Mid observations.

We present a framework to compute non-Gaussian likelihoods for two-point correlation functions. The non-Gaussianity is most pronounced on large scales that will be well-measured by stage-IV weak-lensing surveys. We show how such a multivariate likelihood can be constructed and efficiently evaluated using a copula approach by incorporating exact one-dimensional marginals and a dependence structure derived from the exact multivariate likelihood. The copula likelihood is found to be in better agreement with simulated sampling distributions of correlation functions than Gaussian likelihoods, particularly on large scales. We furthermore investigate the effect of the non-Gaussian copula likelihood on posterior inference, including sampling the full parameter space of contemporary weak-lensing analyses. We find parameter shifts in $S_8$ on the order of one standard deviation for $1 \ 000 \ \mathrm{deg}^2$ surveys but negligible shifts for areas of $10 \ 000 \ \mathrm{deg}^2$, suggesting Gaussian likelihoods are sufficient for stage-IV surveys, though results depend on the detailed mask geometry and data-vector structure.

The Friedmann cosmological solution of the standard Einstein gravitational field equation has a curvature singularity at a moment in time known as the Big Bang. It has been suggested that this Big Bang curvature singularity can be eliminated by use of a degenerate spacetime metric. This proposal was the main topic of our talk at the Workshop, but, here, we also discuss the possible appearance of CPT-conjugated worlds and the conjectured relevance of an extended version of Einstein's field equation.

Since its humble origins, humans have left imprints on the face of the planet. From the profound transformation unleashed by the Neolithic Revolution, about 12000 years ago, till the present, humans have reshaped the planet significantly. From the second half of the XX century, the impact on the atmosphere, biosphere, cryosphere, hydrosphere and upper lithosphere is so overwhelming that a new geological age, the Anthropocene, was proposed to consider the extent of these transformations. However, despite the ubiquitous nature of the changes in course, the International Union of Geological Sciences rejected in March 2024 formalizing the Anthropocene as a new geological epoch. This controversial decision implies that geologists are not quite convinced that human activities have reached the level of an encompassing new geological age. Nevertheless, it is beyond any doubt that there is no single spot on the planet where the signs of the transformations ensued by the human activities are not felt. Furthermore, the interconnection of the human activities has reached a level of entanglement that it makes the Anthropocene an inescapable feature of our present and immediate future. Thus, more important than framing our present condition in a way that it can be recognised by geologists in the future, is the understanding that by its very nature, the Anthropocene is a condition that is continuously being reshaped to the point that we should instead regard our time as a Transientocene, a time of significant and multidimensional transformations.

We examine several observable optical properties of a Skyrmion black hole (BH), focusing on the photon sphere, BH shadow, and photon trajectories. The Skyrme term, along with other geometric parameters of the spacetime, determines the photon sphere location and shapes the resulting BH shadow. Parameter variations produce observable departures from standard BH geometries, offering potential signatures of nonlinear field effects. We also analyze the sparsity of Hawking radiation and the associated energy emission spectra, showing how these quantities respond to the Skyrme coupling and background parameters. Our findings illuminate the connection between nonlinear field contributions and BH optics, with implications for observational and theoretical studies of modified gravity scenarios.

We present a geometric framework for adaptive whitening in gravitational-wave detectors, reformulating the problem from a sequence of spectral factorizations to parallel transport on a principal bundle. We identify the whitening filter as a section over the manifold of power spectra and derive the minimum-phase connection as the unique geometric structure that enforces signal causality while preserving signal-to-noise ratio. This construction yields a rigorous definition of geometric drift, a coordinate-independent scalar measuring the intrinsic instability of the detector noise floor. The central result is the flatness theorem, which proves that the curvature of the connection vanishes for scalar fields. This establishes a holonomic update law, guaranteeing that the optimal filter correction is path-independent and determined solely by the instantaneous noise state, free from geometric phase or hysteresis. This framework unifies the static theory of Wiener-Hopf factorization with the dynamic requirements of real-time control, providing a rigorous certification for the stability of zero-latency calibration routines and establishing a foundation for gauge-theoretic signal processing (GTSP) in next-generation detector networks.

Neutrino telescopes employ multi-level reconstruction chains, where computationally expensive high-quality reconstructions are applied only to events that survive initial quality cuts based on fast, coarse directional estimates. Currently, event selection between reconstruction levels is source-agnostic, giving no priority to events from directions of known neutrino source candidates. We propose a simple modification to inter-level event selection: preferentially retain events whose early-level reconstruction places them within an angular tolerance of pre-specified candidate source directions from established multi-messenger catalogs, while continuing to subsample remaining events at the baseline rate. Using a realistic two-level detector model with energy-dependent angular resolution, we show that this source-informed selection can improve median point source sensitivity by factors of $\sim 2$--$3$ compared to uniform subsampling, with the improvement depending on the baseline selection efficiency, angular tolerance, and correlation between reconstruction qualities at different levels. For catalogs of $\mathcal{O}(100)$ sources, the additional computational overhead is modest ($\sim 7$--$14\%$). This approach offers a path to substantially enhance the discovery potential of current and future neutrino telescopes without requiring new detector capabilities.

The soliton-like coherent structure (SCS), which has been verified to exist in both transitional and turbulent boundary layers1-4, still poses a challenge in the understanding of its formation and behavior. In our previous study (Niu et al.5), the SCS was also found to exist in the transitional wake flow behind a sphere. In present study, the formation and evolution of the SCS is further investigated at four Reynolds numbers by numerical simulation. The results show that at the early stage of the turbulence transition, the SCS appears as a form of wave packet during the Tollmien-Schlichting (T-S) wave stage. With the increase of the Reynolds number, the SCS reaches its maximum amplitude downstream where the velocity discontinuity occurs. This position is located after the breakdown of the T-S wave and the three-dimensional structure is formed. Then, the SCS conserves its shape and amplitude over a long distance downstream. The relationships among the SCS, the spikes, the vortex structures, and the high-shear layers are further analyzed. It is found that the SCS in the wake flow has similarities to the phenomena observed in boundary layer flows during the turbulent transition. The vortex structures and high-shear layers mostly wrap around the border of the SCS. The vortex structure is considered to be as a consequence of the development of the SCS rather than its cause.

This proceedings contribution provides an overview of Loop Blow-up Inflation and updates its observational predictions and their comparison with the latest CMB and BAO data from combined analyses of SPT, Planck, ACT, and BICEP/Keck, as well as ACT DR6 constraints on extra dark radiation. It is based on work originally published in arXiv:2403.04831, carried out in collaboration with L. Brunelli, M. Cicoli, A. Hebecker, and R. Kuespert, and presented at the 2025 Workshop on Quantum Gravity and Strings. We focus on string loop corrections to the Kähler potential, long regarded as a potential threat to blow-up inflation in the Large Volume Scenario. We argue that these corrections, previously assumed avoidable, are in fact generically present and qualitatively alter the original non-perturbative picture: they invalidate slow-roll near the minimum and instead generate a new slow-roll regime at larger field values, where the scalar potential transitions from an exponential to a power-law plateau. This leads to modified inflationary dynamics and distinct cosmological predictions, including an increased tensor-to-scalar ratio. We confront all three SM-location scenarios with the latest constraints on $n_s$, $r$, and $\Delta N_{\rm eff}$. The tighter bound on extra dark radiation requires an updated Giudice-Masiero coefficient in Scenario III, yielding revised predictions presented here for the first time. All scenarios remain consistent with recent observations, with the ACT+DESI combination yielding near-perfect agreement in $n_s$ for vanishing extra dark radiation at $0.03\sigma$ deviation. We also comment on subleading loop corrections, which improve robustness by reducing the field value required for slow roll. These results highlight that string loop effects, rather than being merely detrimental, can play a constructive role in realising viable inflation in string compactifications.

Optically levitated sensors inside a Fabry-Pérot cavity have been proposed for high-frequency gravitational-wave (GW) detection, though their configuration for gravitational wave sensitivity exhibits counterintuitive features. We provide a new detailed general relativistic derivation of the interaction between a gravitational wave and a levitated object in an optical cavity, demonstrating gauge independence of the observable response. We find a strong asymmetric dependence of the strain signal on trap position, maximized when the sensor is located near the input mirror, in agreement with previous results. A key new result of this work is the consequence of this asymmetry on the noise coupling: the coupling of input-mirror displacements to the strain signal can be highly suppressed relative to that of end-mirror displacements and common-mode mirror motion. These results clarify the physical origin of the gravitational wave interaction with such a sensor and establish crucial design principles for optical levitation based high-frequency GW detectors.

Non-thermal fermions can be produced non-perturbatively in the early universe during coherent oscillations of a scalar field. We explore fermion production in $\lambda\phi^{4}$ inflation through this mechanism and analyze the momentum spectrum of the fermions produced, which depends on a coupling parameter $q$. For $q \gtrsim 0.01$, the main contribution to the total number density comes from an approximately half-filled Fermi sphere as a result of non-adiabaticity. For $q\lesssim 0.01$, we find that the major contributions instead come from resonance peaks at higher momentum values. We find a simple relation to predict the momentum values corresponding to resonance peaks for any $q$. We also obtain analytic power-law approximations for the total number density of fermions and find that it is proportional to $q^{1/2}$ for $q\lesssim 0.01$ and proportional to $q^{3/4}$ for $q\gtrsim 10$. If fermions produced by this mechanism make up the entirety of dark matter, we estimate lower bounds on their mass.

Gravitational effective theories are essential for characterizing the space of deviations from General Relativity (GR). Testing these theories against fundamental principles, such as causality and unitarity, can yield constraints on the corresponding parameters. In this paper, we perform such an analysis on the very interesting dynamical Chern-Simons (dCS) gravity. This is a parity violating correction to GR wherein a new scalar field couples to the Pontryagin density $^*R\,R$. It has generated significant interest, including possible new gravitational wave shapes for LIGO/Virgo and new phenomena from cosmic inflation. In this work, we begin by deriving the dispersion relation and wave packet speed on top of a gravitational wave background in dCS gravity. This alters the corresponding Shapiro time delay (which we compute to second order), potentially giving superluminality. Causality then demands a bound on the dCS coupling constant, which we find to be moderately sharper, but compatible with, standard estimates. We then examine a UV completion in the form of a set of $N$ fermions with a (pseudo) Yukawa coupling. By imposing perturbativity and a gravitational species bound, we find that the dCS coupling constant is constrained significantly more, depending on the choice of scale of the species bound. We also identify higher order operators generated from the UV completion. Overall, we find that any dCS corrections to gravitational dynamics should likely be very small on macroscopic systems of observational interest, such as in late-time merging black holes.