Papers

We present a photometric and timing study of the near-period-gap SU~UMa-type dwarf nova YZ~Cancri, based on nearly continuous Transiting Exoplanet Survey Satellite (TESS) photometry and long-baseline ground-based observations. Our main observational result is that the positive-superhump (SH) waveform follows a repeatable, stage-dependent sequence during superoutbursts (SOs). In the two well-covered TESS SOs, and consistently in the long-baseline ground-based SO sample, the plateau waveform evolves from an early saw-tooth profile (ST), through a double-humped profile with unequal maxima (DHd), to a more symmetric double-humped profile with nearly similar maxima (DHs). The global and time-resolved Lomb–Scargle periodograms show that power near the positive-SH time scale and its harmonics is concentrated during SOs, whereas quiescent and normal-outburst intervals lack a comparably persistent SH-band signal. The dense TESS maxima–minima timing sequence shows different clock stability in different waveform stages: precursor modulations have slightly longer local periods, the middle-to-late DHs plateau is the most regular timing interval, and the post-plateau evolution is affected by phase switching and possible secondary/late-SH contamination. The TESS data also reveal profile-clock coupling, with SH amplitude and rise/decay durations evolving together with the timing residuals. As a secondary timing constraint, we obtain a common TESS-timing-based mean positive-SH period of $P'_{\rm sh}=0.09043(27)$~d, corresponding to a SH excess of $4.03(31)\%$ and an approximate mass ratio of $q=0.175(11)$. The repeatable ST–DHd–DHs sequence makes YZ~Cnc a useful system near the lower edge of the period gap, and may trace the growth, redistribution, stabilization, and decay of the light source associated with an eccentric, precessing accretion disk during SOs.

Radio surveys provide exponentially increasing amounts of observational data, which requires the development of novel data reduction and analysis methods. To this end, realistic simulations of observational data provide controlled environments for testing and developing new methods. Diffusion models excel at synthesizing realistic image data, but are limited by size constraints and high computational demands. We implemented a latent diffusion model (LDM) to synthesize realistic radio survey map cutouts. Further, we developed a technique for sampling radio survey maps of arbitrary size with the LDM. This is done by sampling in multiple sequential steps, whereby consistency is provided through context vectors and image inpainting. We trained our model on observations from the LOFAR telescope obtained from the third data release of the LOFAR Two-Metre Sky Survey. We implemented a LDM composed of a vector-quantized variational autoencoder and a diffusion model operating on the resulting latent representations. We added information on source locations, brightness and sizes by extracting parameters from the public source catalog and encoding them into a vector format passed as context to the LDM. We trained the model for image inpainting to facilitate continuous sequential sampling. We were able to generate realistic cutouts of 512 pixel side length, corresponding to ~5.7 arcmin, with precise control over source locations and brightness. Control over source extension was achieved at lower precision. Sampling of arbitrary-sized maps was made possible, showing realistic results with seamless combination of individual samples and limited only by occasional sampling artifacts.

As the Artemis program advances toward the lunar south pole and permanently shadowed regions (PSRs), understanding lunar charging is increasingly important for protecting astronauts, robotic systems, instruments, and infrastructure. Persistent darkness, cryogenic temperatures, low regolith conductivity, and long charge-relaxation times may allow energetic-particle-induced charge to accumulate beneath the surface. This Technical Memorandum addresses a key unresolved question: how is electrostatic charge generated, separated, retained, and accumulated within lunar regolith? Existing models predict that solar energetic particles and galactic cosmic rays may produce subsurface electric fields approaching dielectric breakdown thresholds, but the microscopic connection between incident particles and macroscopic volumetric charge sources remains unclear. Energy deposition alone does not determine retained charge. Incident particles and their secondary particles may stop, implant, backscatter, transmit, recombine, become trapped, or escape. This memorandum therefore defines the required microscopic input as the average signed retained-charge distribution per unit depth and incident particle, resolved by particle species and energy. Combined with incident flux and energy spectra, this response provides a depth-dependent volumetric charge-source rate that can be coupled with charge continuity, conduction, dielectric relaxation, and Poisson's equation. Key uncertainties include the effects of mineralogy, grain and pore geometry, temperature, and pre-existing potential on charge retention. Particle-resolved modeling and cryogenic high-vacuum irradiation experiments are needed to constrain these processes and assess subsurface electric fields, dielectric breakdown, dust transport, contamination, and charge-mitigation requirements for sustained lunar polar exploration.

Fixed angular apertures sample different physical regions of nearby galaxies and can therefore bias optical emission-line diagnostics. We use 20,545 high-confidence emission-line galaxies observed by both the Dark Energy Spectroscopic Instrument (DESI) and the Sloan Digital Sky Survey (SDSS) as a two-fibre experiment, in which the DESI 1.5 arcsec fibre is nested within the SDSS 3 arcsec fibre. To isolate aperture effects from line-measurement systematics, we refit every spectrum using the same eMILES stellar-continuum and emission-line model while preserving the wavelength-dependent line-spread function of each spectrum. This common refit yields 20,200 quality-controlled matched galaxies. Relative to SDSS, the smaller DESI aperture produces small but highly significant offsets: $\Delta\log({\rm N2})=+0.0067\pm0.0002$, $\Delta\log({\rm S2})=-0.0148\pm0.0003$, $\Delta\log({\rm O3})=-0.0246\pm0.0006$, $\Delta{\rm O3N2}=-0.0327\pm0.0007$, and $\Delta\log({\rm H}\alpha/{\rm H}\beta)=-0.0195\pm0.0003$, where $\Delta$ denotes DESI minus SDSS. These offsets persist in Baldwin–Phillips–Terlevich-selected star-forming galaxies and in a near-concentric subset. Their amplitudes vary across the redshift distribution and retain a measurable dependence on apparent fibre coverage after redshift is controlled, showing that redshift and relative aperture coverage jointly modulate the aperture terms. We conclude that, for high-S/N emission-line galaxies, DESI–SDSS comparisons carry measurable, diagnostic-dependent fibre-aperture terms even when continuum modelling, emission-line fitting, and spectral-resolution treatment are held fixed.

We introduce DustRover, an open-source Python package designed to model dust extinction curves using multi-wavelength regimes. DustRover supports optical/infrared spectroscopy, broadband photometry, and X-ray data within a unified framework, enabling robust analysis of dust extinction curves in distant galaxies. The software features a modular architecture that strictly decouples data management, dust models, and Markov Chain Monte Carlo (MCMC) based statistical fitting. The package is highly user-oriented, using simple YAML configuration files and offering built-in visualisation tools for SEDs and extinction curves.

The Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope produces high-dimensional data-cubes where even minor astrometric misalignments between multiple exposures can introduce severe artifacts and compromise data quality. While automated pipelines provide a crucial baseline, subtle offsets require human intervention. To address this, we introduce PERRY, an interactive, Python-based Graphical User Interface (GUI) explicitly designed for the visual inspection and manual fine-tuning of MUSE image alignments. Playing an essential role of finding manual offset adjustments for data-cubes within the Pythonic AAO Reduction Environment (PARE), PERRY ingests automatically pre-aligned data, provides synchronized visual and quantitative diagnostics, and enables real-time coordinate translation and rotation updates. The resulting precise adjustments are preserved in reproducible configuration files, maximizing the scientific return and efficiency of advanced data reduction pipelines.

In many scientific disciplines, weak signals of interest are obscured by dominant nuisance signals that are several orders of magnitude stronger. Recovering these weak signals requires subtracting the dominant ones; however, this calibration process inherently distorts or partially suppresses the underlying signal of interest. To address this problem, we propose the use of normalizing flow models to reconstruct calibration-affected weak signals. By leveraging the statistical invariance of the target signals and assuming minimal initial suppression, our framework effectively recovers the lost signal components. We provide a comprehensive theoretical overview of this normalizing flow-based recovery method and demonstrate its efficacy using simulated data.

We investigate bulk viscosity (BV) effects on post-merger gravitational wave (PMGW) signals from binary neutron star (BNS) mergers and their detectability using realistic equation of state (EoS) posteriors constrained by recent progress in gravitational-wave (GW) and X-ray observations, as well as neutron skin thickness measurements. A linear fitting formula based on recent BNS simulations with BV is used to estimate the peak frequency shift of the PMGWs caused by BV effects. Subsequently, we evaluate this peak frequency shift for our set of EoS posterior samples. We use the Fisher information matrix and the dataset of simulated observable events from the ET and CE detector network in an optimistic scenario to estimate the measurement accuracy of the peak frequency. We find that BV effects on PMGWs may be detectable only for EoS models with large values of the symmetry energy slope ($L_{\rm sym}$), as favored by neutron-skin thickness measurements. Thus, the detection of BV effects on PMGWs itself can provide abundant information about the symmetry energy and its slope. However, BV effects on PMGWs are generally weak and may be observed only in some extreme and optimistic cases. Therefore, observations of PMGWs and their analyses may not be able to provide very accurate measurements of the symmetry energy through BV effects.

We investigated the impact of the number of binary neutron star merger events and neutron star (NS) mass distributions on constraining the equation of state (EoS), tidal deformability and radius of NSs, as well as the nuclear parameters, using binary neutron star inspiral gravitational wave signals with third-generation detectors. We generate simulated gravitational wave signals and compute the Fisher information matrix for each event after the number of events and mass distribution models are given. The covariance of EoS parameters is obtained by holding the non-EoS waveform parameters fixed at their injected values and accumulating the Fisher matrix over all events. Finally, the posterior samples of EoS, tidal deformability, radius and nuclear parameters are generated based on this covariance matrix. We find that larger number of events lead to tighter constraints due to more observed events and data accumulation, as expected given the increased number of detections. Meanwhile, we note that the mass distribution plays a more important role in constraining the EoS. We compare a realistic bimodal Gaussian distribution, a uniform distribution and another uniform including sub-solar mass NSs. The results show that the uniform distribution yields tighter constraints than bimodal models because it includes more low-mass and massive NSs. This also implies that events near $1.4\,\Msun$ provide partly redundant information about the EoS. Additionally, including sub-solar mass NSs can further improve the constraints by significantly reducing the EoS uncertainties at sub-saturation densities and inner-crust, highlighting the importance of sub-saturation density EoS.

We present a multi-wavelength characterization of IPHAS J190812.63+045728.1. Independently discovered by us after being serendipitously observed by XMM-Newton, followed up with optical photometry with SPARC4/OPD and spectroscopy with GMOS/Gemini South, and exploring archival TESS and ZTF observations, we reveal the system to be an eclipsing accreting white dwarf (WD) system with an orbital period of 5.073 h. The deep optical and flat-bottomed X-ray eclipses completely occult the WD and the accretion disk, providing a model-independent, uncontaminated view of the donor star which we classify as an M2.3 dwarf with a mass of $\sim 0.50\,M_\odot$. The system is a modestly luminous (L$_{x;0.3-10 keV} \sim$ 1.7 $\times$ 10$^{32}$ erg s$^{-1}$) hard thermal (consistent with a shock-heated plasma cooling down from $k$T $\sim$ 26 keV) X-ray source. While these X-ray properties are typical of low luminosity intermediate polars (IPs), long-term archival photometry reveals frequent, large-amplitude outbursts ($\Delta m \sim 2$ mag in the ZTF $zr$ band) that are the typical characteristics of dwarf novae, and would be unusual for IPs. The weakness of the HeII $\lambda$4686 emission line does not allow us to decide between the IP and the dwarf nova scenarios. On the other hand, the optical and X-ray light curves do not exhibit periodic variations caused by the WD rotation. If the system is indeed an IP, the accretion geometry is likely unusual while the dwarf nova-like outbursts indicating disk instabilities would make IPHAS J190812.63+045728.1 a benchmark for testing disk instability models in the presence of magnetically truncated accretion disks. Otherwise, the system is a dwarf nova with X-ray properties at the extreme end of its class.

Binary black holes (BBHs) emit gravitational waves (GWs) as they inspiral towards merger. These GWs can be gravitationally lensed by large-scale structure along the line of sight, potentially creating multiple images of the same source with fixed time delays determined by the lensing geometry. As the BBHs inspiral, the GW frequency increases, leading to successive constructive and destructive interference between the multiple images. BBHs also have spins $\mathbf{S}_i$ that may be misaligned with their orbital angular momentum $\mathbf{L}$. As the BBHs inspiral, these misaligned spins cause $\mathbf{L}$ to precess about the total angular momentum $\mathbf{J}$, modulating the GW emission similar to pulsar emission resulting from a misaligned jet rotating in and out of the line of sight. We investigate the ability of a single L-shaped GW detector to distinguish between these two sources of modulation. We find that precessional modulation can mimic the lensing modulation between two images with comparable magnifications when the time delay between the images is short enough that fewer than three interference fringes occur during the time the GW signal spends in the sensitivity band of the detector. As strong lensing is rare for GW sources at moderate redshift while misaligned spins are common for BBHs produced in certain formation channels, ruling out precessional modulation is essential to identifying genuinely lensed systems.

Constraining the plasma properties and energetics of coronal mass ejections (CMEs) in the inner corona is essential for understanding their early evolution, yet remains challenging because of limited observations. We investigate the mass and density evolution, and energy partitioning of two CMEs observed on 2025 September 14 and 16, and assess the potential of coordinated Proba-3/ASPIICS and Aditya-L1/VELC observations to constrain CME energetics in the low corona. We present the first coordinated observations of CMEs obtained simultaneously by Proba-3/ASPIICS and Aditya-L1/VELC. Using ASPIICS white-light observations, we estimate the CME mass, volume, number density, and the evolution of kinetic, thermal, and magnetic energies. Magnetic energies are estimated from observed CME properties using observationally constrained, physically motivated assumptions. Simultaneous VELC Fe XIV 5303 Angstorm observations provide independent estimates of the emission measure, electron number density, thermal energy, and CME lateral extent. The two events exhibit markedly different energy partitioning in the low corona. For the 14 September CME, the kinetic and magnetic energies are comparable, while the thermal energy remains nearly two orders of magnitude smaller, indicating limited plasma heating. In contrast, the 16 September CME exhibits a substantial thermal-energy enhancement, with thermal energy eventually becoming comparable to kinetic energy. For both CMEs, the estimated magnetic energy remains comparable to or exceeds the kinetic energy over the observed height. Our results demonstrate the scientific potential of synergetic ASPIICS and VELC observations for constraining CME mass, density, and energetics in the inner corona, providing new observational constraints on the early evolution of CMEs.

Aims. The objective of this work is to adapt a deep neural network model to perform galaxy morphological classification trained from crowd annotations, considering the training scheme, the agreement between the annotators, and the hierarchy. Methods. We use Galaxy Zoo 1 as our experimental testbed and trained a convolutional neural network (CNN) for the automatic classification of galaxies' morphologies. We analyze the impact of the following aspects on the classification accuracy and training efficiency: (i) Training only the last layer vs. training all the network; (ii) Classification with only the CNN vs. considering the hierarchy; (iii) Comparing the models trained with different amounts of data and levels of agreement between the annotators; (iv) Training by stages, transferring knowledge from one model to another; and (v) Combining several models as an ensemble. Results From the experiments, we derive the following results: (i) Training all the layers in the network significantly improves the accuracy (10% increase in exact match), compared to training only the last layer; (ii) There is a tradeoff between the amount of data and the level of agreement between the annotators used for training; (iii) Using the hierarchy can improve accuracy when the amount of training data is reduced; (iv) Training by stages through transfer learning (curriculum learning) produces higher accuracy for limited data; (v) Ensembles can improve accuracy; (vi) Models achieve a low accuracy for the most difficult cases, but, if we consider hierarchical measures, we can derive useful results for upper levels in the hierarchy. An accuracy above 99% is achieved when training all layers of the network and considering a high agreement between the annotators. Conclusions. Training deep learning models from crowd annotations involves additional challenges than learning from hard annotations.

The aim of this work is to predict the opacity of plasmas under stellar conditions. We focus on iron and nickel, as these elements have been extensively investigated both theoretically and experimentally. In certain regimes, notably under non-local thermodynamic equilibrium, calculating the spectral opacity can be computationally demanding. To mitigate this difficulty, we employ deep learning models to predict accurate opacities while achieving a substantial reduction in computational cost. Specifically, we develop and train a hybrid model that combines a convolutional neural network with a multilayer perceptron. The proposed approach provides both the spectral opacity over the 0–10,000 eV energy range, and the corresponding mean opacities for large sets of temperatures and mass densities, in close agreement with standard computational methods. The agreement is very satisfactory, except in the L- and M-shell energy ranges for nickel at high temperatures, where a very large number of emerging transitions leads to substantial discrepancies.

We present interferometric observations with the Northern Extended Millimeter Array (NOEMA) of NGC 1333 IRAS 4B and 4B', two young stellar objects (YSOs) within the Perseus molecular cloud. YSOs, especially Class 0 protostars, are rich in gas-phase interstellar complex organic molecules (iCOMs), but relatively little is known about the spatial distribution of these important prebiotic molecules. Various molecules were imaged within IRAS 4B and the surrounding area in the 145 GHz spectral range, revealing hundreds of molecular transitions within the warm inner envelope of IRAS 4B. This work provides derived physical parameters (rotational temperature, column density, velocity shift with respect to local standard of rest, and full-width at half-maximum spectral line width) for 11 molecules toward IRAS 4B including methanol (CH$_3$OH and isotopologues), methyl formate (HCOOCH$_3$), dimethyl ether (CH$_3$OCH$_3$), acetaldehyde (CH$_3$CHO), ethanol (C$_2$H$_5$OH), glycolaldehyde (CH$_2$OHCHO), acetone (CH$_3$COCH$_3$), and isocyanic acid (HNCO). Out of these molecules, 10 are iCOMs (excluding HNCO) that trace the inner envelope of IRAS 4B and mark it as a rich source of organic material. There is a strong spatial correlation between the continuum emission of IRAS 4B and iCOM molecular emission and among various iCOMs with each other. Physical parameter and integrated intensity maps show that IRAS 4B is a binary protostar system that drives a jet in a north-south orientation and an outflow in a northwest-southeast orientation. IRAS 4B' also drives a west-east outflow but shows no evidence of molecular emission cospatial with its continuum emission.

Recent observations of the local Little Red Dot (LRD) analog SDSS~J1025+1402 have revealed that the optical continuum is polarized at the ~1.5% level, while broad Halpha has lower polarization (~0.7%) and a different polarization angle. We explain these properties in a scenario in which LRDs are dust-obscured Little Blue Dots (LBDs) viewed at high inclination and powered by super-Eddington accretion. In this framework, electron scattering in the geometrically thick inner accretion flow produces intrinsic continuum polarization, while unpolarized emission from the outer thin disk dilutes the signal, particularly at optical wavelengths. A circumnuclear dust screen contributes dichroic polarization to both the continuum and broad-line emission. Broad Halpha lacks the intrinsic disk component and therefore has lower polarization and a different position angle. We show that this model quantitatively reproduces the observed polarization properties, although degeneracies between accretion rate and inclination remain. These results illustrate the potential of spectropolarimetry as a probe of accretion-flow geometry and orientation in the little-dot population, complementary to spectral-energy-distribution fitting and line diagnostics.

Continuous gravitational waves are long-lived signals emitted by spinning neutron stars (NSs) or boson clouds around black holes. While the Milky Way is expected to host $\mathrm{O}(10^{8}-10^{9})$ NSs, only $\mathrm{O}(10^{3})$ are currently identified through electromagnetic observations. This large observational gap motivates searches based on alternative messengers. In particular, all-sky searches for continuous gravitational waves (GWs) offer a unique opportunity to detect NSs that are electromagnetically silent or otherwise undiscovered. By exploring a broad region of parameter space without any a priori source information, these searches can probe the vast hidden NS population of our Galaxy. In this work, we present a novel way to include a Neural Network (NN) classifier into an all-sky search strategy for isolated NSs. Starting from candidates produced by the FrequencyHough pipeline running on data from ground-based detectors such as LIGO and Virgo, the proposed method improves sensitivity without increasing computational cost and can be parallelized across multiple detectors to enhance detection probability and reduce false alarms. We analyze data from the third observing run (O3) — April 1, 2019 to March 27, 2020 — and focus on the frequency range [129, 229] Hz and spin-down interval [$-2.5\cdot10^{-9}, 1.0\cdot10^{-9}$] Hz/s. The model is trained on noise constructed to be statistically consistent with real O3 data and tested on real O3 data. Results show robust performance in distinguishing signal from noise, even at low strain, and successful identification of hardware injections outside the training spin-down range.

T Coronae Borealis is the most compelling anticipated Galactic recurrent nova and one of the best opportunities to test whether nova eruptions can become genuine multi-messenger transients. Recent observational work on classical novae has shown that the GeV $\gamma$-ray luminosity correlates strongly with the differential velocity between interacting outflows, suggesting that shock power depends strongly on velocity contrast. We apply this empirical framework to T CrB and compare it with the 2021 eruption of RS Ophiuchi, the only nova securely detected at TeV energies. We argue that T CrB may lie at the extreme high-$\Delta v$ end of the nova population and that its strongest shock should arise at the external interaction between the nova ejecta and the slow circumbinary medium. This makes the external shock the most likely site of any TeV emission, while the GeV signal may contain contributions from both internal, if present, and external shocks. Because T CrB is substantially closer than RS Oph, it is an especially favorable target for TeV detection. The neutrino case is promising but more uncertain because it depends sensitively on the density and structure of the circumbinary material. We also revisit the recurrence time of T CrB. Using the observed scatter of better-sampled recurrent novae, we infer a representative population-scale recurrence of $T_{\rm rec}\simeq80\pm12$ yr. Taken together, its proximity and potentially extreme shock velocities make T CrB perhaps the best opportunity for a nova to become the first securely established multi-messenger nova.

The recent discovery of S301 provides an unique opportunity to investigate dark matter-baryon interaction in the innermost region of the Milky Way. S301 follows a highly eccentric orbit around Sgr~A*, and its pericenter distance is only around 12 AU with velocity about $25000\,{\rm km/s}$. These extreme properties make S301 highly sensitive to dark matter-baryon scattering in the center of galaxy, where the dark matter density may be strongly enhanced. We have calculated the orbit-averaged energy deposited rate in S301 via dark matter-proton scattering, and show that over a broad dark matter parameter space, cross sections satisfying the existing bounds from terrestrial experiments and astronomical observations can inject sufficient energy to modify the luminosity and evolution of S301. Our results establish S301 as an exceptionally sensitive probe of dark matter-baryon interactions and open a new avenue for probing dark matter interactions in high-density and high-velocity astrophysical environments.

We present a novel model for inelastic dark matter in a simple extension of the Standard Model (SM) with a local $U(1)'$ baryon flavor symmetry, such as $U(1)_{B_i-B_j}$ with $i\neq j$. In this framework, we realize a realistic flavor structure for the mixing Yukawa couplings due to higher dimensional operators or renormalizable couplings between the SM quarks and extra vector-like quarks that are anomaly-free by themselves. We regard a Dirac singlet fermion as being dark matter and obtaining a small mass splitting for inelastic dark matter, only after the $U(1)'$ symmetry is broken spontaneously. We show that there is a consistent parameter space for explaining both the LZ nuclear recoil event and the correct relic density for dark matter and the existing collider and indirect detection bounds are satisfied at the same time. In this case, we predict a relatively light $Z'$ boson with $10-20\,{\rm GeV}$ mass and a weak gauge coupling, suggesting a dedicated search at colliders. The future data with higher recoil energies at direct detection experiments will play a decisive role of discriminating between the scenarios of light and heavy dark matter.

Deep-space exploration beyond the solar system requires eliminating chemical propellant mass penalties to achieve relativistic flight velocities (0.166c at 180 s, reaching the mission target of 0.20c at 227 s). This study presents a 3D multi-physics numerical framework for a 4.0-meter circular lightsail propelled by a 100 GW ground laser array, coupling 3D Maxwell FDTD wave optics, non-linear membrane mechanics, and Stefan-Boltzmann thermal radiation in ANSYS Mechanical APDL and Ansys Lumerical. A four-level grid convergence study establishes numerical independence with an ASME GCI_21 of 0.13%, resolving peak membrane stresses of 530.88 MPa with a 3.77x safety factor against stoichiometric Si3N4 tensile failure. With optical absorption constrained to 10 ppm (A = 1.0 x 10^-5), the steady-state core temperature stabilizes at 923.02 K (0.44% deviation from radiation theory), maintaining a +1,247 K margin below sublimation, while fundamental drumhead modal resonance (7.92 Hz) provides a 7.92x safety buffer against laser jitter. The electrodynamic radiation pressure formulation is cross-verified against published flight telemetry from JAXA IKAROS and NASA LightSail 2 within 0.12% and 2.13%, confirming classical momentum transfer modeling across solar and beamed propulsion regimes.

We present an implementation of an HLLD approximate Riemann solver for the AthenaK astrophysics code with support for full general relativistic magnetohydrodynamics via a tetrad frame transformation. Our implementation uses an initial guess for the HLLD iterative solve which eliminates the need for an additional conserved-to-primitive inversion, which greatly accelerates performance without affecting accuracy. Additionally, by coupling the method with a first-order flux correction, we are able to use the method reliably even when the magnetization exceeds $10^4$, which we achieve in a SANE accretion disk. Our SANE disk shows that HLLD leads to a more strongly magnetized funnel and more accurate horizon fluxes when compared with HLLE. We further apply the new HLLD implementation to an equal-mass binary neutron star merger. For our long-lived remnant, HLLD enhances the magnetic shear stresses in the outer layers and leads to weaker differential rotation. However, due to weaker gravitational wave emissions in the post-merger phase, the remnant is consistently less compact while producing more dynamical ejecta and a more massive disk. The cost of this new solver is relatively modest thanks to the improved initial guess: our accretion disk tests are only ${\sim}10-25\%$ slower than HLLE, and for our binary neutron star runs with a microphysical equation of state, we find that HLLD is only ${\sim}3\%$ slower than HLLE across all runs.

The recently launched Nancy Grace Roman Space Telescope (Roman) will soon begin observations that promise to revolutionize our understanding of exoplanet demographics. Roman's Galactic Bulge Time-Domain Survey (GBTDS) will observe the Galactic Bulge via transit and microlensing techniques, providing a unique opportunity to study hot Jupiters as well as long-period/free-floating exoplanets in a dense and kinematically hot stellar environment. In this Letter, we quantify how repeated, impulsive stellar flybys truncate the outskirts of planetary systems in the Bulge. We find that flybys result in ejections on $\sim$Gyr timescales for wide orbits (${\gtrsim} 250\ \text{au}$) around low-mass hosts. Provided these wide orbits are populated, these ejections likely contribute to the population of free-floating planets probed by Roman. Surviving planets outside of ${\sim} 10$ au are placed on eccentric and misaligned orbits. We show that, though systems suffer many close encounters, the evolution is dominated by the single strongest perturbation; planetary orbits (even those that evade ejection) undergo superdiffusion (i.e. a Lévy flight). We comment on the implications of our results for the origin of hot Jupiters amenable to $\textit{Roman}$'s transit search. Our work serves as a quantitative baseline from which future dynamical studies can build to better understand the evolution of planetary systems in the Galactic Bulge.

The Magellanic Stream (MS) is believed to have formed from gas removed from the Large and Small Magellanic Clouds through tidal forces and hydrodynamic interactions with the Milky Way's gaseous halo. It provides an important laboratory for studying how stripped gas fragments, mixes, and evolves in a circumgalactic environment. In this work, we present HI observations of a $3.8^\circ\times2.2^\circ$ field in the MS IV region, using the data from the Commensal Radio Astronomy FasT Survey (CRAFTS). The total HI mass in the analyzed field is $\simeq 5.3 \times 10^{6} (d/120\,{\rm kpc})^{2}\,M_{\odot}$, where $d$ is the distance to the MS IV gas. The data resolve the emission into three coherent filamentary HI structures with related but distinct velocity trends. To characterize the HI structures and study potential multiphase gas, we adopt a Gaussian decomposition procedure to identify and reconstruct sources. Our results indicate that the field is dominated by one major filamentary HI complex, together with several smaller kinematic clump-like components. Among these identified sources, only one source likely shows multiple velocity components. Because several sources appear spatially overlapped in projection, we further examine their apparent overlap regions using position-velocity (P-V) diagrams. P-V diagrams across the apparent overlap region show no clear intermediate-velocity bridge or V-shaped structure, favoring line-of-sight projection over direct cloud-cloud collision. These results demonstrate the value of deep, high-angular-resolution HI observations for resolving faint emission, compact morphologies, and kinematic structure in the MS.

Orbit-averaged equations are widely used to model the gravitational-wave-driven evolution of compact binaries. Their validity relies on the radiation-reaction timescale being much longer than the orbital period, an assumption that can break down for highly eccentric systems. In this work, we investigate the impact of orbit averaging on merger-time estimates by comparing orbit-averaged and non-orbit-averaged prescriptions. We apply our analysis to two different astrophysically motivated populations of black hole binaries in hierarchical triple systems, where von Zeipel-Lidov-Kozai oscillations can drive the inner binary eccentricities close to unity. We find that merger-time estimates obtained with orbit-averaged equations remain remarkably accurate across a wide range of eccentricities, even when the orbital timescale exceeds the radiation-reaction timescale. Significant discrepancies arise only for the most extreme and compact binaries, simultaneously characterized by very high eccentricities and small semi-latus recta. In this regime, the merger time becomes strongly dependent on the initial orbital phase, with different phases yielding predictions that can differ by several orders of magnitude. We also compare our results with an orbit-averaged merger-time prescription specifically developed for hierarchical triple systems and find substantial differences between this prescription and both the generic orbit-averaged and non-orbit-averaged calculations in the high-eccentricity regime. While these results demonstrate the robustness of orbit-averaged merger-time estimates for most systems, they also highlight the potential importance of non-orbit-averaged effects when merger times are used to infer compact-binary delay-time distributions and formation channels.

We iteratively model the GWTC-5.0 binary-black-hole census with three components. First, a physically normalized population from Rapster globular-cluster simulations spans cluster mass, metallicity, formation redshift, compactness, and natal black-hole spin. Hierarchical mergers reproduce higher-mass events and their larger effective-spin dispersion, provided black holes are born with zero natal spin. Second, a phenomenological field (isolated-binary) channel supplies the low-mass, preferentially aligned population. Its inclusion flattens the cluster compactness likelihood, permitting ordinary dense globular-cluster birth radii without requiring nuclear-cluster-like conditions. Finally, residual tension at $15$–$30 M_\odot$ motivates an intermediate-mass isotropic component representing field remnants reprocessed in clusters. This is a hypothesis, not a detection: current numerical support precludes a reliable evidence comparison. Physical normalization converts the high-mass rate into $\hat{f}_{\rm GC}\simeq0.39\%$, with local rates $R_{\rm cl}\simeq9.1$ and $R_{\rm field}\simeq16.0$ Gpc$^{-3}$ yr$^{-1}$. The model predicts linked mass-spectrum breaks near $35$ and $70 M_\odot$, a $q\simeq0.5$ feature from first-plus-second-generation pairings, and a symmetric effective-spin distribution that broadens sharply above $45 M_\odot$. These correlated, mass-resolved predictions can be tested as the gravitational-wave census grows.

The detection of the redshifted 21 cm signal from neutral hydrogen during the epoch of reionization remains one of the central challenges in contemporary cosmology. Although upcoming radio interferometers are expected to measure 21 cm fluctuations at a few percent precision, an independent confirmation is essential to ensure robustness against instrumental systematics and foreground contamination. In this work, we propose a novel real-space stacking technique to detect the 21 cm signal, guided by the presence of dark gaps in quasar absorption spectra during the late stages of reionization ($z \sim 5.0$-$6.1$). Using a semi-numerical reionization model, we show that stacking 21 cm images on long dark gaps (length, $L > 32~h^{-1}\mathrm{cMpc}$ and threshold flux, $F_{\mathrm{th}} < 0.005$) should reveal a strong emission signal with a peak mean-subtracted 21 cm brightness temperature of $\delta \tilde{T}_b \sim 2-3~ mK$ for a volume-weighted average neutral fraction of $x_{\mathrm{HI}} \sim 0.1$ at $z = 5.9$. We further extend our analysis to realistic reionization lightcones, incorporating thermal noise and foreground avoidance effects, generated for two different reionization histories. In an optimistic foreground avoidance scenario (wedge field of view $\theta < 20^\circ$), a sample of roughly 100 dark gaps is sufficient to achieve a combined root-sum-square signal-to-noise ratio exceeding 10 for both reionization models, while a more conservative wedge avoidance may require more dark gaps for a solid detection. Given the large number of quasars known at $z > 6$, a detection of this stacked signal may thus already be within reach of existing optical (i.e., VLT, Keck) and 21 cm observatories (i.e., HERA, SKA). This synergistic approach offers a promising pathway towards a direct and statistically robust detection of the cosmological 21 cm signal.

Space-based gravitational-wave (GW) data may contain transient signals whose waveform morphologies are not known in advance. Extracting these signals is important for characterizing new sources and mitigating instrumental anomalies. However, existing deep neural network (DNN)-based extraction approaches rely on clean training targets and waveform-class-specific examples, which may limit their applicability when the transient morphology is not specified in advance. This work develops a Noise2Noise (N2N)-inspired self-supervised framework that learns from noisy observations without clean training targets and requires no transient-specific waveform templates at inference. A single model trained on noisy massive black-hole binary (MBHB) observations provides high-overlap MBHB recovery and recovers the dominant morphologies of instrumental glitches and other GW transient signals such as cosmic string bursts in source-confused test data. Beyond waveform recovery, when independent information identifies a candidate transient as instrumental, its extracted waveform can be subtracted from the data without an anomaly-specific template. In a simulated continuous data stream, this procedure substantially suppresses the injected-anomaly power within the glitch-dominated frequency band. These results support the use of self-supervised extraction for initial waveform estimation of candidate transients with unknown morphologies, enabling subsequent characterization and, where appropriate, conditional subtraction of instrumental anomalies.

The hydromagnetic structure of magnetically confined mountains on accreting white dwarfs is computed, alongside the effects of accretion on the reduction of the star's magnetic field. The equilibrium structure of the mountain is obtained by numerically solving the Grad-Shafranov equation with a self-consistent scheme that enforces magnetic flux freezing. For characteristic system parameters, it is shown that magnetic field lines are significantly deformed and dragged equatorwards by mountains with masses $\sim10^{-3}M_\odot$ that have a characteristic height of $\sim10^5\,\mathrm{cm}$. A mathematical relation is obtained for the maximum mass that a white dwarf can confine magnetically given the star's magnetic field and temperature. It is found that the accretion buries the star's magnetic field by reducing it by up to $10\%$ of its pre-accretion value. The star's magnetic field remains dominantly dipolar, but with added multipolar components.

We present a homogeneous HST/ACS F555W/F814W time-series census of variable stars in two M101 disk fields, using a uniform multi-band reduction and reproducible variable-selection workflow for nearby SN-host galaxies. Supplementary archival F555W epochs are used only to refine selected long-period solutions. Candidates are selected with Lomb-Scargle period searches, Fourier light-curve modeling, two-band consistency statistics, artificial-star-test photometric corrections, and mode-aware Cepheid refinement. The final catalog contains 1417 variable-star candidates, including 1102 secure fundamental-mode and 78 secure first-overtone classical Cepheids. Of these secure Cepheids, 420 are newly identified relative to the Shappee & Stanek catalog. The catalog also contains luminous candidates outside the classical Cepheid locus. We release the catalog and source-level two-band light curves for studies of Cepheid physics, stellar variability, and distance-scale applications.