Locally authored papers of the past 5 days

Compilation of papers presented by the VERITAS Collaboration at the 39th International Cosmic Ray Conference (ICRC), held July 14 through July 24, 2025 in Geneva, Switzerland.

The JWST discovery of "Little Red Dots" (LRDs) has revealed a population of compact, red sources at $z \sim 5-10$ that likely host supermassive black holes (SMBHs). Recent observations of the gravitationally lensed LRD R2211-RX1 reveal century-scale photometric variability and a hysteresis loop in the luminosity-temperature plane, strongly suggesting that the optical emission originates from a pulsating, stellar-like photosphere rather than an accretion disk. This supports the "quasi-star" hypothesis, where a rapidly growing black hole seed is embedded within a massive, radiation-pressure supported envelope. In this work, we investigate the stability of these envelopes using the stellar evolution code MESA coupled with the non-adiabatic oscillation code GYRE. We identify a theoretical "Quasi-Star Instability Strip" with a blue edge at $T_{\mathrm{eff}} \approx 5000-5200$ K. Models hotter than this threshold are stable, consistent with the non-variable LRD R2211-RX2 ($T_{\mathrm{eff}} \approx 5000$ K), while cooler models are unstable to radial pulsations driven by the $\kappa$-mechanism in helium and hydrogen ionization zones. For quasi-star masses in the range $M_\star \sim 10^4-10^5 M_\odot$, we find that the unstable fundamental radial modes ($\ell =0$, n$_{\rm p}=1$) have periods in the range $\sim 20-180$ years. The first overtone ($\ell =0$, n$_{\rm p}=2$) is also unstable or marginally stable in some of our models, with typical pulsation timescales $\sim 10-30$ years. These oscillations match the co-moving frame variability timescale of RX1. We argue that these violent pulsations likely drive enhanced mass loss analogous to super-AGB winds, which could affect the duration of the quasi-star phase and regulate the final mass of the seeded black hole.

Mars lacks a global intrinsic dipole magnetic field, but its interaction with the solar wind generates a global induced magnetosphere. Until now, most studies have relied on single-spacecraft measurements, which could not simultaneously capture upstream solar wind conditions and the induced magnetic fields, thereby limiting our understanding of the system. Here, we statistically re-examine the properties of Mars' induced magnetic field by incorporating, for the first time, real-time upstream solar wind conditions from the coordinated MAVEN and Tianwen-1 observations. Our results are show that both solar wind dynamic pressure and the interplanetary magnetic field (IMF) magnitude enhance the strength of the induced magnetic field, but they exert opposite effects on the compression ratio: higher dynamic pressure strengthens compression, while stronger IMF weakens it. The induced field is stronger under quasi-perpendicular IMF conditions compared with quasi-parallel IMF, reflecting a stronger mass-loading effect. We further investigate the clock angle departures of the induced fields. They remain relatively small in the magnetosheath near the bow shock, increase gradually toward the induced magnetosphere, and become significantly larger within the induced magnetosphere. In addition, clock angle departures are strongly enhanced under quasi-parallel IMF conditions. Their dependence on upstream drivers further shows that, within the magnetosheath, clock angle departures are minimized under low dynamic pressure, high IMF magnitude, and low Alfven Mach number conditions. These results may enhance our understanding of solar wind interaction with Mars, and highlight the critical role of multi-point observations.

It is plausible that the dark matter particles have non-gravitational interactions among themselves. If such self interactions are large enough, they could leave an imprint on the morphology of galaxies. These effects can be studied with numerical simulations, which serve as the primary tool to predict the non-linear evolution of galactic structure. A standard assumption is that the course-grained phase-space distribution of the macroscopic simulation particles follows the same evolution equation as that of the fundamental dark matter particles. This Letter tests this assumption directly for the case of frequent dark matter scatterings, demonstrating that this is not generically true. Specifically, we develop a first-principles map from a microscopic particle physics description of self-interacting dark matter to a representation of macroscopic simulation particles for theories in the short-mean-free-path regime. Using this procedure, we show the emergence of an effective force between the simulation particles and derive their interaction cross section, which depends on the one from fundamental particle physics. This work provides the first explicit map from particle physics to simulation, which will facilitate exploring the phenomenological implications for galactic dynamics.

The Sloan Digital Sky Survey IV DR17 Extended Baryon Oscillation Spectroscopic Survey (eBOSS) consists of 2,233,939 high-quality optical galaxy spectra obtained through 2" fibers, providing a rich spectroscopic resource for studying galaxy evolution across a broad redshift range. eBOSS was designed primarily for large-scale structure and BAO measurements and, as such, focused on galaxy redshifts, leaving much of the information contained in the spectra unexplored. In addition to the trove of spectra, the large number of repeat observations (197,521 duplicate spectra) enables evaluation of the survey's spectrophotometric quality. To unlock this potential, we introduce the eBOSS Data Analysis Pipeline (eBOSS-DAP), adapted from the MaNGA-DAP, which delivers uniform measurements of emission-line fluxes and equivalent widths, stellar and gas kinematics, continuum spectral indices, and stellar population fits. Using the eBOSS-DAP, we successfully analyze 1,899,553 high-quality galaxy spectra below a redshift of $z < 1.12$ to produce an extensive spectroscopic catalog for the eBOSS galaxy sample. We characterize the calibration performance, quantify the reliability of the derived measurements, and release a suite of data products that fully exploit the power of the eBOSS dataset. These catalogs open the door to a new generation of studies in galaxy evolution and cosmology.

Be stars are rapidly rotating main-sequence (MS) stars that play a crucial role in understanding stellar evolution and binary interactions. In this letter, we propose a new formation scenario for black hole (BH) + Be star binaries (hereafter BHBe binaries), where the Be star is produced through the Wind Roche Lobe Overflow (WRLOF) mechanism. Our analysis is based on numerical simulations of the WRLOF process in massive binaries, building upon recent theoretical work. We demonstrate that the WRLOF model can efficiently form BHBe binaries under reasonable assumptions on stellar wind velocities. Using rapid binary population synthesis, we estimate the population of such systems in the Milky Way, predicting approximately $\sim$ {1800-3200} currently existing BHBe binaries originating from the WRLOF channel. These systems are characterized by high eccentricities and exceptionally wide orbits, with typical orbital periods exceeding 1000 days and a peak distribution around $\sim$10000 days. Due to their long orbital separations, these BHBe binaries are promising targets for future detection via astrometric {and interferometric} observations.

Solar type II radio bursts are commonly attributed to coronal shocks driven by coronal mass ejections (CMEs). However, some metric type II bursts have occasionally been reported to occur in the absence of a CME and to be associated with weak solar activities. This study aims to identify the driver of the coronal shock in this kind of type II event. We investigate a high-frequency metric type II burst with clear band splitting, observed simultaneously by the Chashan Broadband Solar radio spectrograph and the Nançay Radioheliograph. It is associated with a C3.1-class flare and a small-scale jet, but without a detectable CME in the coronagraphs. The type II burst is preceded by multiple type III bursts, one of which exhibits characteristics of a type N burst. The type II burst source is associated with the jet-induced perturbation front propagating through nearby closed loops at a speed of $\sim$880 km s$^{-1}$, rather than the much slower jet front. This suggests that the disturbance initiated by the jet can convert to a shock wave within low Alfvénic coronal loops, providing the necessary conditions for electron acceleration and subsequent radio emission. Our findings offer new insights into the formation mechanism of high-frequency type II bursts associated with weak flares and jets.

We use weak gravitational lensing measurements from Year 3 Dark Energy Survey data to calibrate the masses of 443 galaxy clusters selected via the Sunyaev-Zel'dovich effect from Atacama Cosmology Telescope Data Release 5 maps of the cosmic microwave background. We incorporate redshift and SZ measurements for individual clusters into a hierarchical model for the stacked lensing signals and perform Bayesian analyses to constrain the hydrostatic mass bias of the clusters. Our treatment of systematic uncertainties includes a prescription for measuring and accounting for the weak lensing boost factor, consideration of a miscentering effect, as well as marginalization over uncertainties in the source galaxy photometric redshift distributions and shear calibration. The resultant constraints on the normalization of the mass-observable relation have a precision of approximately 7\%, with the mean WL halo mass of $M_{\rm 500c} = 5.4 \times 10^{14} M_{\odot}$. We measure the bias between the true cluster mass and the mass estimated from the SZ signal based on an X-ray--calibrated scaling relation assuming hydrostatic equilibrium, to be $1-b = 0.75^{+0.04}_{-0.06}$ over the full sample. When splitting the clusters into high ($z$=0.43-0.70) and low ($z$=0.15-0.43) redshift bins, we measure $1-b = 0.58^{+0.06}_{-0.05}$ and $0.82^{+0.07}_{-0.07}$, respectively. When introducing additional freedom in redshift and mass to the hydrostatic bias model, we find that $1-b$ decreases with redshift (with the power law of $-2.0^{+0.7}_{-0.4}$, 99.95\% confidence), consistent with findings from other recent studies, while we do not find any significant trend in mass. We also demonstrate that our result is robust against various systematics. The weak-lensing mass calibration presented in this study will be a useful tool for using the ACT clusters as probes of astrophysics and cosmology.

Improved measurements of $B$-modes in the cosmic microwave background can be obtained through accurate calibration of the orientation of detector antennas as projected onto the sky. Miscalibration of the detector polarization angle leads to a leakage of $E$-modes into $B$-modes, which can bias the detection of the latter. To achieve a $\sigma(r)$ of 0.003, the Simons Observatory Small Aperture Telescopes are required to calibrate the global polarization angle on the sky with an accuracy ${\lesssim}0.1^\circ$. We demonstrate a fully remote-controllable calibration system using a ``sparse wire grid," which injects a rotatable linear polarized signal across the telescope's focal plane. This calibration system is installed and operational on a Small Aperture Telescope at its observing site at the Parque Astronómico in the Atacama desert in Chile. We developed a pipeline for the detector polarization angle calibration, and demonstrate it using initial data for 93~GHz and 145~GHz frequency bands. The observed distribution of detector polarization angles is in agreement with the instrument design. Statistical uncertainties for the relatively calibrated polarization angles are $0.02^\circ$ and $0.03^\circ$ at 93~GHz and 145~GHz, respectively. Systematic uncertainty was evaluated to be $0.08^\circ$ at the hardware development and fabrication stage. Their sum in quadrature is less than $0.1^\circ$.

Quasi-periodic MicroPulses (QMP) are quasi-periodic microstructural features manifested in individual pulsar radio pulses, the study of which is crucial for understanding pulsar radiation mechanisms. Manual identification of QMP in large-scale pulsar single-pulse datasets remains highly inefficient. To address this, we propose a Dual-Stage Residual Network (DSR) that achieves automated QMP detection in FAST observational data through joint analysis of single-pulse profiles and their Amplitude Distribution Profiles (ADP), defined as the power spectra of the autocorrelation function derivatives of the microstructure residuals. The model was trained on PSR B1933+16 data from 2019 (10,486 single pulses) and evaluated on manually annotated PSR B1933+16 data from 2020 (9,657 single pulses). DSR achieved 96.10\% recall and 95.85\% precision on the test set. This approach provides an automated pipeline for large-scale, reproducible QMP identification and establishes the foundation for in-depth investigation of their physical mechanisms.

Searching for fleeting radio transients like fast radio bursts (FRBs) with wide-field radio telescopes has become a common challenge in data-intensive science. Conventional algorithms normally cost enormous time to seek candidates by finding the correct dispersion measures, of which the process is so-called dedispersion. Here we present a novel scheme to identify FRB signals from raw data without dedispersion using Machine Learning (ML). Under the data environment for multibeam receivers, we train the EfficientNet model and achieve both exceeding 92% accuracy and precision in FRB recognition. We find that the searching efficiency can be significantly enhanced without the procedure of dedispersion compared with conventional softwares like TransientX and presto. Specifically, the impact of radio frequency interference (RFI) for single-beam and multibeam data has been investigated, and we find ML can naturally mitigate RFI under the multibeam environment. Finally, we validate the trained model on actual data from the current FRB surveys carried out by the Five-hundred-meter Aperture Spherical radio Telescope, which provides considerable potential for real implementation in the future.

We present the design, implementation, and operation of a pyramidal horn radio telescope built for detecting the Galactic 21 cm neutral hydrogen line emission. The system employs an SDR-based pipeline to obtain drift-scan observations, which were calibrated and processed to generate HI sky maps, a Galactic rotation curve and spiral arm features. This demonstrates that this low-cost system is effective both for educational purposes and scientific exploration of Galactic structure at radio frequencies.

A deep understanding of the life-cycle of galaxies, particularly those of high mass, requires clarifying the mechanisms that regulate star formation (SF) and its abrupt shutdown (quenching), often capable of stopping SF rates of hundreds of solar masses per year. What initially triggers quenching, and what sustains the quiescent state thereafter, especially given the frequent presence of large gas reservoirs or even massive gas inflows, are unsolved key issues. Ultimately, the crucial connection between the galaxy life-cycle and the surrounding Intergalactic (IGM) and Circumgalactic (CGM) Medium remains largely unclear. Addressing these issues requires studying star formation, chemical enrichment, and quenching homogeneously up to high redshift. The upcoming AO-assisted Extremely Large Telescope (ELT), will deliver sharper and deeper data than the JWST. SHARP is a concept study for a near-IR (0.95-2.45 mu) spectrograph designed to fully exploit the capabilities of ELT. Designed for multi-object slit spectroscopy and multi-Integral Field spectroscopy, SHARP points to achieve angular resolutions (~30 mas) far superior to NIRSpec at JWST(100 mas) to decipher and reconstruct the life-cycle oa galaxies.

HESS J1857+026 remains a mysterious gamma-ray emitter since its discovery in 2008. Despite the disclosure of a nearby pulsar and multiple studies in the high-energy (HE, E > 100 MeV) and very-high-energy (VHE, E > 100 GeV) regimes, there have been no confirmed counterparts (e.g., an SNR shell or other extended structure) in X-ray or other wavelengths. We present the result of our study of the VHE emission of HESS~J1857+026 with VERITAS as part of a multi-wavelength investigation to uncover its emission mechanisms. Our result confirms the extended nature of the source and we characterize its spectral and morphological features in the VHE band. Using the morphology of the source revealed in our analysis, we also explore the underlying transport process of a possible electron population in a leptonic PWN scenario for the gamma-ray emission.

We consider various possible consequences of time-varying dark energy due to a quintessence scalar field whose energy density is partially converted to particles as the field evolves down its potential. This particle production acts as a source of thermal friction on the field that can make it difficult to distinguish whether dark energy is due to a radiating field rolling down a steep potential, a purely self-interacting field moving down a flatter potential, or a cosmological constant. By reducing the acceleration of the scalar field, thermal friction increases the amount of accelerated expansion and can cause a sizable bump in the quintessence equation of state. We take special interest in the case where a steep potential rapidly changes from positive to negative as the field evolves, resulting in the end of cosmic expansion and the beginning of contraction. Even in this case, we find that thermal friction lengthens the period of accelerated expansion and consequently delays the end of cosmic expansion, making it challenging to detect the impending transition to contraction using conventional cosmological tests. However, particle production can also provide alternative avenues for detection by generating a background of thermal dark radiation, partly comprised of neutrinos or other particles, whose energy density exceeds the remnant photon energy density.

We present updated forecasts for the scientific performance of the degree-scale (0.5 deg FWHM at 93 GHz), deep-field survey to be conducted by the Simons Observatory (SO). By 2027, the SO Small Aperture Telescope (SAT) complement will be doubled from three to six telescopes, including a doubling of the detector count in the 93 GHz and 145 GHz channels to 48,160 detectors. Combined with a planned extension of the survey duration to 2035, this expansion will significantly enhance SO's search for a $B$-mode signal in the polarisation of the cosmic microwave background, a potential signature of gravitational waves produced in the very early Universe. Assuming a $1/f$ noise model with knee multipole $\ell_{\rm knee} = 50$ and a moderately complex model for Galactic foregrounds, we forecast a $1\sigma$ (or 68% confidence level) constraint on the tensor-to-scalar ratio $r$ of $\sigma_r = 1.2\times10^{-3}$, assuming no primordial $B$-modes are present. This forecast assumes that 70% of the $B$-mode lensing signal can ultimately be removed using high resolution observations from the SO Large Aperture Telescope (LAT) and overlapping large-scale structure surveys. For more optimistic assumptions regarding foregrounds and noise, and assuming the same level of delensing, this forecast constraint improves to $\sigma_r = 7\times10^{-4}$. These forecasts represent a major improvement in SO's constraining power, being a factor of around 2.5 times better than what could be achieved with the originally planned campaign, which assumed the existing three SATs would conduct a five-year survey.

The observation of filamentary X-ray structures near bow-shock pulsar wind nebulae (PWNe) -- such as the Guitar, Lighthouse, and PSR J2030$+$4415 nebulae -- and of slow-diffusion regions around pulsars like Geminga, Monogem, and PSR J0622$+$3749, challenges the standard picture of cosmic-ray transport in the interstellar medium, implying a diffusion coefficient two orders of magnitude smaller than the Galactic average. The suppressed diffusion can be attributed to self-generated magnetic turbulence, driven -- via the non-resonant streaming instability -- by electron--positron pairs escaping the PWNe. This instability requires a net current, yet the beam of escaping pairs is expected to be charge-neutral. We show that a charge-neutral pair beam propagating through an electron--proton plasma can spontaneously generate a net current. Using fully kinetic two- and three-dimensional particle-in-cell simulations with realistic mass ratio, we find that beam electrons get focused into self-generated magnetic filaments produced by the nonlinear evolution of the Weibel instability, while beam positrons remain unconfined. The resulting net (positron) current drives the non-resonant streaming instability, further amplifying the magnetic field. This mechanism provides a pathway for the onset of charge asymmetries in initially charge-neutral pair beams and for the growth of magnetic fluctuations that efficiently scatter the beam particles, with implications for the formation of X-ray filaments and, more broadly, for particle self-confinement in TeV halos around PWNe.

Thermal images at different wavelengths probe varying subsurface depths of planetary bodies, and therefore can inform us about their compositions, thermophysical properties, and impact histories. We identified six archival observations of Callisto obtained by the Atacama Large Millimeter/submillimeter Array (ALMA) between 2012 July 17 and 2012 November 4 at wavelengths of 0.43-0.47 mm (701.9-641.5 GHz). These wavelengths are shorter than those of nearly all other Callisto ALMA data and are sensitive to subsurface emission at depths (the upper ~cm) between those sounded by millimeter and infrared observations. We estimate the disk-averaged brightness temperature as 133$\pm$15 K, and use a thermophysical mixture model to find that Callisto's thermal emission is best fit by a ~50-50% two-component thermal inertia mixture of $\Gamma_{\text{low}}$~15-40 and $\Gamma_{\text{high}}$~1200-2000 J m$^{-2}$ K$^{-1}$ s$^{-1/2}$, consistent with recent ALMA observations of Callisto at longer wavelengths. Finally, we present several previously unpublished Galileo Photopolarimeter-Radiometer (PPR) observations of Callisto and derive thermal inertia and spectral emissivity maps using the same model. Altogether, these ALMA and PPR maps improve our understanding of the thermal properties and spatial distribution of Callisto's shallow subsurface regolith, and demonstrate the value of ALMA flux density calibrator data for extending frequency coverage of existing science data.

Strong magnetic fields and plasmas are intrinsically linked in both terrestrial laboratory experiments and in space phenomena. One of the most profound consequences of that is the change in relationship between the frequency and the wave number of electromagnetic waves propagating in plasma in the presence of such magnetic fields when compared to the case without these fields. Furthermore, magnetic fields alter electromagnetic wave interaction with relativistic plasma waves, resulting in different outcomes for particle and radiation generation. For a relativistic plasma wave-based photon acceleration this leads to an increased frequency gain, and, thus, potentially to higher efficiency. The influence of a magnetic field leads to quantitative and qualitative change in the properties of photon acceleration, amplifying the increase in the electromagnetic wave frequency.

The spin-orbit tilt angles $\theta_{1(2)}$ of merging stellar-mass black holes provide key insights into their astrophysical origin. The LIGO, Virgo, and KAGRA Collaborations (2025a, arXiv:2508.18083) report that the spin-orbit tilt distribution of mergers in the latest Gravitational-Wave Transient Catalog 4.0 exhibits a global peak at near-perpendicular directions $\cos\theta_{1(2)}\approx0$. Here, we recover this feature using hierarchical Bayesian inference with parametric models that are tailored to enhance the diagnostic power about astrophysical formation channels. We find that the spin distribution of the low-mass bulk of the binary black hole merger population $(m_1\lesssim 44.3^{+8.7}_{-4.6}\,\rm M_\odot)$ can be well-modelled by a dominant Gaussian component that peaks at $\cos\theta_{1(2)}\approx0$, possibly mixed with a subdominant isotropic component. Models that include a component with spins preferentially aligned with the orbit are disfavoured by current data (with Bayes factors $|\Delta\ln\mathcal{B}|\approx1$ to $3$) and constrain its contribution to be small ($\xi\sim\mathcal{O}(1)\,\%$). If these findings are reinforced by more detections, they would challenge any major contribution from the traditional isolated-binary formation scenario yielding closely aligned spins. Instead, the dominant component with near-perpendicular spins qualitatively matches expectations from the evolution of isolated massive stellar triples in the galactic field, where the Lidov-Kozai effect naturally produces a unique overabundance of mergers with $\cos\theta_{1(2)}\approx0$.

Accurate redshift estimates are a critical requirement for weak lensing surveys and one of the main uncertainties in constraints on dark energy and large-scale cosmic structure. In this paper, we study the potential to calibrate photometric redshift (photo-z) distributions for gravitational lensing using the Dark Energy Spectroscopic Instrument (DESI). Since beginning its science operations in 2021, DESI has collected more than 50 million redshifts, adding about one million monthly. In addition to its large-scale structure samples, DESI has also acquired over 256k high-quality spectroscopic redshifts (spec-zs) in the COSMOS and XMM and VVDS fields. This is already a factor of 3 larger than previous spec-z calibration compilations in these two regions. Here, we explore calibrating photo-zs for the subset of KiDS-1000 galaxies that fall into joint self-organizing map (SOM) cells overlapping the DESI COSMOS footprint using the DESI COSMOS observations. Estimating the redshift distribution in KiDS-1000 with the new DESI data, we find broad consistency with previously published results while also detecting differences in the mean redshift in some tomographic bins with an average shifts of Delta Mean(z) = -0.028 in the mean and Delta Median(z) = +0.011 in the median across tomographic bins. However, we also find that incompleteness per SOM cell, i.e., groups of galaxies with similar colors and magnitudes, can modify n(z) distributions. Finally, we comment on the fact that larger photometric catalogs, aligned with the DESI COSMOS and DESI XMM and VVDS footprints, would be needed to fully exploit the DESI dataset and would extend the coverage to nearly eight times the area of existing 9-band photometry.

Understanding the magnetic field environment around Mars and its response to upstream solar wind conditions provide key insights into the processes driving atmospheric ion escape. To date, global models of Martian induced magnetosphere have been exclusively physics-based, relying on computationally intensive simulations. For the first time, we develop a data-driven model of the Martian induced magnetospheric magnetic field using Physics-Informed Neural Network (PINN) combined with MAVEN observations and physical laws. Trained under varying solar wind conditions, including B_IMF, P_SW, and {\theta}_cone, the data-driven model accurately reconstructs the three-dimensional magnetic field configuration and its variability in response to upstream solar wind drivers. Based on the PINN results, we identify key dependencies of magnetic field configuration on solar wind parameters, including the hemispheric asymmetries of the draped field line strength in the Mars-Solar-Electric coordinates. These findings demonstrate the capability of PINNs to reconstruct complex magnetic field structures in the Martian induced magnetosphere, thereby offering a promising tool for advancing studies of solar wind-Mars interactions.

We investigated whether "velocity excess" in circumstellar maser lines can diagnose the earliest evolutionary phases of Water Fountains (WFs). Here we define "velocity excess" as maser emission (e.g., H$_2$O 22.235 GHz or OH 1665/1667 MHz) detected at velocities outside the velocity range of the OH 1612 MHz line, which traces the terminal expansion velocity of a spherical circumstellar envelope (CSE). Such excess velocities serve as an indicator of gas motions deviating from spherical expansion and may signify the onset of asymmetric outflows. Based on recent studies (Fan et al.2024; Xie et al.2025), we analyzed 17 sources showing velocity excess and fitted their infrared spectral energy distributions (SEDs) with the one-dimensional radiative transfer code DUSTY. Seven sources are well reproduced, implying outer CSEs that remain nearly spherical despite inner asymmetries. Among these, five exhibit single-peaked, AGB-like SEDs and two show double-peaked, post-AGB-like profiles. IRAS variability indices and NEOWISE-R light curves reveal pulsations (~600-1000 days) in three sources, supporting their AGB classification. Considering the magnitude of the velocity excess, two objects-IRAS 19229+1708 and IRAS 19052+0922-may represent the earliest or incipient WF phase, in which asymmetric outflows are beginning to emerge within otherwise spherical envelopes. These results support a morphological sequence in which bipolar jets and tori arise first in the central regions while the outer CSE remains spherical, and they show that selecting WF candidates via velocity excess effectively identifies objects at the onset of jet formation and early morphological transformation.

We explore a novel analysis framework for parameter inference with large-scale CMB polarization data. Our method uses simulation-based inference combined with the needlet internal linear combination (NILC) algorithm and cross-correlation-based statistics to compress the data into a vector that is robust to model misspecification and small enough to be amenable to neural posterior estimation with normalizing flows. By leveraging this compressed data representation, our method enables the robust use of the anisotropic and non-Gaussian information in the foreground fields to more accurately separate the CMB polarization signal from these contaminants. Using an idealized ground-based experimental setup inspired by the Simons Observatory Small Aperture Telescopes, we demonstrate improved statistical constraining power for the tensor-to-scalar ratio $r$ compared to the (constrained) NILC algorithm and improved robustness to complex foregrounds compared to other techniques in the literature. Trained on a relatively simple semi-analytical foreground model, the method yields unbiased $r$ results across a range of PySM Galactic foreground simulations, including the high-complexity d12 model, for which we obtain $r=(1.09 \pm 0.27)\cdot 10^{-2}$ for input $r=0.01$ and sky fraction $f_{\mathrm{sky}} = 0.21$. We thus demonstrate the feasibility and advantages of a complete, maps-to-parameters, simulation-based analysis of large-scale CMB polarization for current ground-based observatories.