2026-04-16 | | Total: 84
We present Visplot, a free, open-source, web-based tool for hardware-aware visibility analysis and heuristic scheduling of both sidereal and non-sidereal astronomical observations. Visplot computes visibility windows as finite unions of disjoint intervals by intersecting user-defined constraints. This framework natively incorporates celestial parameters (airmass, moon distance, twilight), mechanical telescope boundaries (altitude and hour-angle limits), and custom temporal restrictions defined in UTC or Local Sidereal Time, allowing for a high degree of scheduling flexibility. The scheduling engine combines deterministic pre-allocation for mandatory targets with a multi-objective heuristic optimization of the remaining target pool, balancing scientific priority, target urgency, altitude, and telescope slew overhead. Originally developed to address an operational need for flexible and lightweight scheduling support at the Nordic Optical Telescope (NOT) in La Palma, Visplot has been in continuous use since 2016. Its nearly decade-long operational history, together with routine use by astronomers at multiple observatories worldwide, demonstrates its practical value in real-world observational workflows. Its client-side, zero-installation architecture facilitates real-time schedule refinement, making it particularly suited for time-domain triggers (e.g., GRB/GW alerts) and geographically distributed remote observing. A user survey indicates that the tool significantly reduces the cognitive overhead of nightly planning while ensuring that generated schedules remain strictly within the mechanical and operational limits of the telescope hardware. Visplot provides a robust, lightweight alternative to monolithic scheduling suites, supporting the practical needs of modern PI-led observatories.
We present a blind time-delay cosmography measurement of the Hubble constant $H_0$ based on the quadruply imaged quasar SDSSJ1433+6007. Our analysis combines deep Hubble Space Telescope imaging, extended time-delay monitoring from the Wendelstein and Maidanak Observatories, and spatially resolved stellar kinematics from the Keck Cosmic Web Imager and Reionization Mapper. We build a robust lens model to reconstruct the mass distribution and high-signal-to-noise kinematic maps to break the mass-sheet degeneracy (MSD), explicitly accounting for the lens galaxy's oblateness, rotation, and anisotropy. Furthermore, we constrain the external convergence ($κ_{\rm ext}$) by characterizing the line-of-sight environment using wide-field photometry from the Dark Energy Spectroscopic Instrument (DESI) Legacy Survey data release 10. We incorporate these constraints into our joint lensing and dynamical model, running multiple iterations to estimate random and systematic uncertainties. Accounting for maximal flexibility of the mass-sheet transformation, and assuming a flat $Λ$CDM cosmology and an $Ω_{\rm m, 0}$ prior from DESI data release 2, we infer $H_0 = 73.2^{+4.8}_{-4.7}$ km s$^{-1}$ Mpc$^{-1}$ (a $6.5\%$ precision), and an internal mass-sheet parameter $λ_{\rm int}=1.12^{+0.05}_{-0.06}$. Notably, $λ_{\rm int}$ is $2σ$ away from unity for this system, highlighting the importance of treating it as a free parameter. Our $H_0$ measurement is consistent with the result from our 2025 milestone paper, and it will be included in our next hierarchical analysis to improve the overall precision. Moving forward, the comprehensive pipeline demonstrated herein establishes a robust framework that can be readily applied to future strongly lensed systems to further refine cosmological constraints.
We present high-resolution (0.94" $\approx$ 55 pc) ALMA CO(2-1) and 13CO(2-1) observations of the highly inclined (i~87.5 deg) galaxy NGC 4565 covering out to galactocentric radius Rgal > $\pm$ 17 kpc. The combination of sensitivity and resolution enables the detection of CO emission well into the HI-dominated outer disk while isolating individual molecular clouds across the full extent of the galaxy. Although often described as an edge-on Milky Way analog, the molecular gas profile of NGC 4565 has a central gap which is more similar to M31. The 13CO/12CO ratio remains consistent at 0.086 $\pm$ 0.009 from Rgal = 5-13 kpc. Based on fits to position-velocity slices, the molecular disk remains thin, with a FWHM scale height of 79.1 $\pm$ 1.6 pc measured from the vertical intensity profile with little evidence for vertical flaring. Molecular clouds in NGC 4565 show sizes, linewidths, and surface densities consistent with those found in similar environments in PHANGS-ALMA galaxies and in M31. We identify a prominent star-forming complex on the ring-an overdensity of molecular gas we term the East Ring Pileup. This feature hosts a compact, multiwavelength-bright region, which we call the Jewel. Effects of galaxy inclination on molecular cloud radius, velocity dispersion, surface density, and virial parameter appear as second-order effects that are strongest in velocity dispersion. At this resolution, GMCs are preferentially aligned with the disk of the galaxy and horizontally elongated by a factor of~2.
Protoplanetary disk ice lines shape a multitude of planet formation processes, setting the environmental composition through evolution. Ice line locations depend on molecular sublimation and deposition properties, but in dynamic disks where temperature and density structures change, so do the expected compositions of planets and planetesimals. In turbulent viscous disks with particle drift, thermal evolution, and desorption/adsorption, Price et al. 2021 demonstrated that the CO/H$_2$O ice ratio beyond the CO ice line can become enhanced by $\sim10\times$. We expand on their work by incorporating additional carbon, nitrogen, and oxygen species, more particle sizes, and a broader disk parameter exploration. We find that before $\sim0.5$Myr, volatile ices are enhanced relative to H$_2$O as the outer disk is desiccated by drift, while at later disk times outward advection and volatile deposition further increase relative volatile icy enhancements beyond the evolving critical disk radius. The outcome of these combined relative icy enhancement to H$_2$O mechanisms is solid C/O $\sim$ N/O $\sim1$ beyond the hypervolatile ice lines, much higher than expected in static disks. Hypervolatiles (N$_2$, CO, and CH$_4$) robustly increase to $\sim100\times$ across the explored parameter space, while mid-volatiles (CO$_2$ and NH$_3$) are sensitive to model choices, with enhancements ranging from $\sim2-50\times$. Together these results demonstrate that coupling disk dynamics with simple sublimation and deposition chemistry is fundamental to predicting grain, planetesimal, and planetary compositions, particularly the role of advection in redistributing volatiles across disk radii.
Star evolution models predict the lightest compact objects in the universe to have masses greater than that of the Sun. Nonetheless, alternative scenarios could lead to the formation of sub-solar mass (SSM) compact objects, such as primordial black holes (PBHs). The LIGO-Virgo-KAGRA Collaboration (LVK) has performed a search for gravitational-wave (GW) signals from compact binary coalescences (CBCs) including at least one SSM component during the first part of their fourth observing run (O4a), reporting no statistically significant candidate. This non-detection sets upper limits on the merger rate of such systems, which can be used to constrain PBH formation models and the fraction of dark matter (DM) in PBHs. For PBH binaries forming at late times, the fraction of DM in PBHs is constrained to be <= 1 for masses above 0.9 M_sun in the case of monochromatic mass functions. In the early-formation scenario, this fraction is limited to $\le 7\%$ at 1 $M_\odot$ and $\le 40\%$ at 0.35 $M_\odot$.
Filamentary structures in molecular clouds are thought to play a central role in star formation, yet the statistical distribution of their mass per unit length -- and any connection to the stellar initial mass function (IMF) -- remains poorly constrained observationally. Here we present a systematic analysis of the filament linear density function (FLDF) across seven nearby molecular clouds spanning $\sim$100--1000\,pc in distance and a wide range of star-forming environments, from quiescent to massive, using the multiscale extraction method $getsf$. The median linear densities of filaments increase approximately linearly with spatial scale, $\tildeΛ \propto Y$, and the fraction of supercritical filaments varies widely among clouds, from a few per cent to over 50%. Only when integrating over the full hierarchy of spatial scales do the combined linear densities across all clouds yield a FLDF that follows a power law ${\rm d}N/{\rm d}\logΛ\propto Λ^{-α}$ with $α\approx 1.30$--$1.34$, mirroring the Salpeter IMF slope of $1.35$. Our results suggest that the stellar mass spectrum is pre-encoded in the filamentary structure of molecular clouds, providing a direct observational link between the large-scale distribution of interstellar gas and the universal slope of the stellar IMF.
Using integral field spectroscopy from SDSS-IV MaNGA, we investigate the radial distributions of star formation rate (SFR) and gas-phase metallicity in spiral galaxies that reside in spiral-elliptical (S+E) pairs. Spirals in S+E pairs show suppressed central star formation and elevated metallicities, whereas spirals in spiral-spiral pairs exhibit centrally enhanced star formation and reduced metallicities. The degree of SFR suppression and metallicity enhancement in S+E pairs depends on the masses of the pair members. Spirals with more massive elliptical companions experience stronger star-formation suppression and larger increases in metallicity, while lower-mass spirals show more pronounced metallicity enhancement. In addition, within S+E systems, galaxies with asymmetric gas velocity fields display enhanced SFR and higher metallicities, whereas those with symmetric velocity fields exhibit clear central suppression. Based on these results, we infer that in S+E pairs, the spiral galaxy experiences suppressed gas accretion once it enters the hot circumgalactic medium of its early-type companion, which leads to the observed decline in star-formation activity. When a close encounter takes place, tidal perturbations can compress the remaining cold gas and trigger enhanced star formation, producing rapid chemical enrichment and the associated increase in metallicity.
We present BROOM, a new python package for the application of blind, minimum-variance component-separation techniques to microwave observations. The package enables the reconstruction of signals with known spectral energy distributions, such as the Cosmic Microwave Background (CMB), Sunyaev--Zeldovich distortions, or foreground moments, in both temperature and polarization through a suite of Internal Linear Combination (ILC) implementations, in the presence of astrophysical and instrumental contaminants. In addition, BROOM supports the blind reconstruction of coherent emission components with unknown covariance properties via a Generalized ILC (GILC) framework. Beyond component separation, the package provides tools to diagnose foreground complexity and to estimate residual contamination leaking into reconstructed maps across angular scales and sky regions. It also includes utilities to generate realistic microwave simulations for arbitrary CMB experiments and to compute angular power spectra of the resulting products. We present a comprehensive description and validation of the implemented pipelines in two representative experimental configurations: a full-sky satellite mission and a ground-based experiment. BROOM is publicly available, fully documented, and easily installable at https://github.com/alecarones/broom
An observation of a large prominence on the solar limb took place on February 15, 2022, by the Extreme Ultraviolet Imager (EUI) on board Solar Orbiter. We aim to determine the range of physical parameters of this prominence, such as temperature, radial velocity, and altitude, and examine how these parameters affect the formation of the 304~Å line of \ion{He}{2}, especially how collisional excitation and radiative processes contribute to line formation. After constraining these parameters, we generate 200 random models and compute the \ion{He}{2} 304~Å line profile. We present these results using parallel coordinate plots to explore how these parameters affect the results. This allows us to infer the key physical parameters that impact the formation of the \ion{He}{2} 304~Å line. This study demonstrates that column mass and the steepness of the temperature profile are key factors in the formation of the \ion{He}{2} 304~Å line during the solar prominence eruption on February 15, 2022. Radiative processes remain dominant in the formation of the \ion{He}{2} 304~Å line. These insights provide a foundation for future research and comparative studies.
The first dedicated observation of an off-limb prominence by the Spectral Imaging of the Coronal Environment (SPICE) instrument took place on April 15, 2023. We aim to create parameter maps on the prominence region, including temperature, pressure, and column mass, by studying the integrated intensity of the Lyman $β$ and Lyman $γ$ lines from SPICE data. After constraining the altitude and radial velocity in the prominence, we use a 1D non-LTE radiative transfer code to generate 1000 random models and compute the Lyman $β$ and Lyman $γ$ line profiles. The computed intensities are compared with observed integrated intensities from SPICE. Then, we create models which simultaneously give a reasonable match with the observed intensities in both lines. Unlike previous approaches, our method uses contribution functions to guide the optimisation of temperature and pressure profiles. Our approach enables a physically constrained and consistent match to both spectral lines. The method in this paper enables us to generate models from pixels on the prominence region and use this information to generate parameter maps. The results obtained have potential for future research.
The first dedicated observation of an off-limb prominence by Solar Orbiter took place on April 15, 2023. Our aim is to determine the range of different physical parameters of this prominence and to examine how these parameters affect the formation of the Lyman $β$ and Lyman $γ$ lines of hydrogen. We have found a way to refine key physical parameters by observational data. We will test the method by this prominence observation. We generate 200 random non-LTE models using these observational constraints and compute the Lyman $β$ line and the Lyman $γ$ line profiles. We use the Spectral Imaging of the Coronal Environment (SPICE) full-disk mosaic from November 13, 2023 to constrain the incident radiation. We present the parameters and results of 200 random models using parallel coordinate plots to explore how different parameters affect the results. This allows us to infer the key physical parameters (e.g., central pressure, column mass and temperature gradient) that impact the formation of the Lyman $β$ line and the Lyman $γ$ line in this observation.
The first dedicated observation of an off-limb prominence by the Spectral Imaging of the Coronal Environment (SPICE) instrument on board Solar Orbiter took place on April 15, 2023. Our aim is to provide an overview of the potentiality of the diagnostics using these data. We show that we can derive the changes in the physical parameters of the pre-eruptive and eruptive prominence using the Lyman lines. We investigate the integrated intensity and line widths of the Lyman $β$ and Lyman $γ$ lines, finding variations between the prominence, disk, and coronal regions. The results reflect dynamic changes in density, temperature, and optical thickness. We analyze the spatial and temporal evolution of the Lyman $β$ and Lyman $γ$ line profiles. Using a simple geometric model, we obtain the radial velocity of this prominence at the early phase of its eruption with GONG H$α$ images. This offers a way of calculating the radial velocity of an eruptive filament from a pair of 2D images. The result helps us understand the potential Doppler effect in line profiles. Overall, the spectral profiles indicate that the eruption enhances spatial and temporal variations in line intensity, reflecting dynamic changes in plasma conditions within the prominence. These findings highlight the diagnostic potential of SPICE observations, and future Non-LTE radiative transfer modeling will help to further constrain prominence plasma parameters.
We present multi-epoch high-resolution optical spectroscopy ($R \simeq 80{,}000$) of the emission-line object NaSt1 to test its proposed binary nature, along with long-term multiband photometry, mid-infrared spectroscopy, and spatially resolved integral field unit (IFU) spectroscopy to probe the circumstellar kinematics of the system. We detect two groups of 35 emission lines showing radial velocities (RVs) variation of the same period of 310 $\pm$ 6 d, but with opposing phase, which we associate with the optically thick wind of the stripped primary star and the wind-wind collision region with the companion star, providing a strong evidence for binarity. The RV and light curve (LC) periods are consistent within the uncertainties, ruling out ellipsoidal modulation, which would require an orbital period of about 620 d. The phase relationship between the RV and LC is inconsistent with stellar pulsations and supports a binary origin. We model the 1--5~$μ$m spectrum of NaSt1 and find two optically thin dust components: hot $T_{\rm h} \simeq 1230$ K, $M_{\rm h} \simeq 2 \times 10^{-10} M_{\odot}$ and warm $T_{\rm c} \simeq 660$ K, $M_{\rm c} \simeq 3 \times 10^{-8} M_{\odot}$. IFU spectroscopy spatially resolves the circumstellar medium in the [\ion{N}{2}] $\lambda6548$ and $\lambda6584$ emission lines, showing a deprojected expansion velocity of $\sim31$ km~s$^{-1}$, implying a dynamical age of $\sim40$ yr. This short timescale suggests that the nebula was produced by recent mass loss. The system may represent a Galactic analog of a massive binary undergoing a mass-loss process to become a stripped-envelope supernova progenitor.
We present a detailed investigation of the shock properties associated with solar energetic particle (SEP) events that exhibit a concave (``nose-like'') shape in their energy spectrogram, characterized by inverse velocity arrival (IVA) of the particles, where high-energy particles arrive later than mid-energy ones. Using measurements from Solar Orbiter and Parker Solar Probe between 2018 and 2025, we identify 26 such SEP events and reconstruct the observed shock fronts in three dimensions. We derive shock parameters along the magnetic field lines connected to each spacecraft using kinematic modeling and coronal magnetohydrodynamic simulations. Our analysis indicates that IVA-SEP events arise due to the spatial and temporal evolution of the shock properties and magnetic connectivity. In most cases analyzed here, the magnetic connectivity starts on the flanks of CME-driven shocks, where shocks tend to be weak, and shifts toward the shock apex sampling stronger portions of the shock front. This evolution of the shock properties at the connected field lines likely leads to the delayed arrival of high-energy particles and the progressive hardening of the SEP energy spectrum, observed in some of the events. We find a correlation between the transition energy at which the IVA begins and the shock speed along the connected field lines, consistent with expectations from time-dependent diffusive shock acceleration. Our results underscore the importance of the evolving shock properties, magnetic connectivity, and instrumental sensitivity in shaping SEP intensity profiles and the formation of IVA signatures.
Gravitational waves (GWs) from the galactic double white dwarf (DWD) systems are one of the primary targets for upcoming space-based detectors. Due to their vast abundance and widespread distribution throughout the Galactic disk and bulge, these systems may provide a high-statistical population for probing GW microlensing effects induced by Galactic compact objects. To evaluate the detectability of such effects, in this work we simulate the four-year observation of DWD systems by Taiji, in the form of a second-generation Time Delay Interferometry (TDI) data stream. Within a Bayesian inference framework, we estimate parameters for lensed GWs from DWD systems for different values of the lens parameters, including the lens mass $M_\mathrm{L}\in [10, 10^6]$\,M$_\odot$, the effective velocity $v_\mathrm{eff}\in [50, 500]$\,km/s and the initial separation $L\in [R_\mathrm{E}, 3R_\mathrm{E}]$, and obtain the uncertainties of the corresponding parameters. These results characterize the capability of future Taiji observations to probe such systems. We further employ the Bayesian model selection framework to distinguish between lensed and unlensed scenarios, and investigate the impacts of three key physical parameters of the lens system: $M_\mathrm{L}$, $v_\mathrm{eff}$, and $L$ on distinguishing lensing events. Our results show that when $M_\mathrm{L}$ is below $10^5$\,M$_\odot$ or $L\geq3R_\mathrm{E}$, it is not possible to distinguish between lensed and unlensed models. For $v_\mathrm{eff}$, although the Bayes factor decreases as $v_\mathrm{eff}$ decreases, the lensed and unlensed models can still be distinguished within our parameter range.
We investigate the structure of the optical F-corona (inner zodiacal light) using a publicly released wide-field image of a total solar eclipse obtained during the Artemis II crewed lunar flyby. In this image, the solar disk is fully occulted by the Moon, providing a rare view of diffuse circumsolar emission over a wide angular extent. Although the dataset is a rendered RGB JPEG image without full photometric calibration, field stars show that the instrumental response is approximately linear in the unsaturated regime, allowing analysis of the relative morphology of the F-corona. The observed F-corona exhibits a flattened, nearly elliptical morphology aligned with the ecliptic plane, with flattening indices of 0.52-0.59 across the three color channels. The east-west intensity distribution is approximately symmetric, while a modest north-south asymmetry is present. Radial intensity profiles along ecliptic longitude and latitude are well described by power laws in solar elongation and are broadly consistent with previous observations by STEREO/HI-1A and LASCO. Comparison with the ZodiSURF zodiacal light model indicates that the observed F-corona is more extended than model predictions, particularly in the north-south direction. Within the model framework, we explore this discrepancy through an empirical, phenomenological modification of the number density distribution of the interplanetary dust. Within an empirical, phenomenological framework, a shallower radial number-density dependence with a power-law index of ~0.7 improves the match to the observed relative profiles.
Radio astronomy is uniquely coupled to exponential trends in computation because the optics (cross-correlation, beamforming, and imaging) and spectrometry (i.e. channelization) can now be done digitally. Inexpensive analog-to-digital converters (ADCs) can sample signals from large numbers of antennas and graphics processing units (GPUs) allow us to coherently process wide-field radio data in real time, motivating large-$N$ aperture arrays at moderate cost. We describe the 256-antenna Coherent All-Sky Monitor (CASM-256), a dense aperture array operating at 375-500\,MHz, currently being deployed at the Owens Valley Radio Observatory (OVRO) in Big Pine, California. The large field-of-view (FoV$\sim10^4$\,deg$^2$) and point-source sensitivity of CASM-256 will allow it to detect local Universe fast radio bursts (FRBs). The nearby sample is ideal for unveiling the physical origin of FRBs, measuring the baryonic content of nearby galaxy halos, and discovering prompt multi-wavelength and multi-messenger counterparts to FRBs. CASM will search for fast transients in the Milky Way such as FRB analogs, pulsar giant pulses, and the new source class known as long-period radio transients. We describe the instrument and present on-sky data from the first two dozen antennas, including an operational real-time GPU based FRB search pipeline. We emphasize the scalability of the concept and describe paths to a future CASM array with tens of thousands of antennas that could detect one million FRBs.
The KM3NeT Collaboration reported the detection of a neutrino, designated as KM3-230213A, with a reconstructed energy peaking at 220 PeV and equatorial coordinates (J2000) of RA=$94.3\degree$ and Dec=$-7.8\degree$. As the highest-energy neutrino event documented to date, its astrophysical origin remains unascertained. Prior preliminary investigations have probed potential associations between this neutrino event and gamma-ray bursts (GRBs), factoring in the possibility of Lorentz invariance violation (LV). In this study, we perform a comprehensive analysis to explore correlations between KM3-230213A and all viable GRBs. We explicitly account for the angular uncertainties intrinsic to both the neutrino event and the respective GRBs. Our analysis identifies a larger set of correlated GRBs. For each associated GRB, we compute the LV scale, integrating uncertainties from redshift measurements and neutrino energy determinations to enhance the robustness of our findings.
Early observations from the James Webb Space Telescope (JWST) have revealed an overabundance of massive high-redshift galaxies, raising the question of whether this points to new physics beyond $Λ$CDM, or an enhanced formation efficiency of massive stars. We revisit this issue going beyond earlier analyses based on direct comparisons to theoretical bounds at a fixed cosmology, by performing a full Bayesian analysis of the most extreme galaxies in the CEERS imaging and FRESCO spectroscopic samples, jointly constraining cosmological parameters and the baryon-to-star conversion efficiency $ε$. We do so not only within the spatially flat $Λ$CDM model, but also in models where the dark energy equation of state $w$ and/or the spatial curvature parameter $Ω_K$ are allowed to vary, carefully discussing the impact of both $w$ and $Ω_K$ on the cumulative comoving stellar mass density. Within the flat $Λ$CDM model, once cosmological parameters are marginalized over, the CEERS sample provides a weak $2σ$ lower limit of $ε\gtrsim 0.07$, compatible with astrophysical expectations. In contrast, the FRESCO sample requires $ε\gtrsim 0.5$ at $2σ$, with values $ε\lesssim 0.2$ disfavored at $>5σ$. These results do not qualitatively change when we allow $w$ and/or $Ω_K$ to vary, with no evidence for deviations from $w=-1$ or $Ω_K=0$. Our results therefore suggest that the origin of the ``JWST tension'' is unlikely to be cosmological, but lies in the astrophysics of galaxy formation.
Direct observational constraints on the earliest, stellar-wind-dominated phases of galactic outflows remain scarce. We present medium-resolution VLT/X-shooter spectroscopy of six Type I superluminous supernova (SLSN-I) host galaxies at z = 0.15-0.51, exploiting the bright SLSN continua as single, down-the-barrel probes of the host interstellar medium. From nebular emission lines we derive dust-corrected star-formation rates as low as 0.06-0.44 Msun yr$^{-1}$, and gas-phase metallicities in the extremely metal-poor regime (less than nine percent solar). Moreover, all hosts exhibit narrow, blueshifted Mg II 2796, 2803 absorption, indicative of the presence of low-ionization outflows along the line of sight. Voigt modeling of the Mg II absorption yields maximum outflow velocities of $v_{max}$ = 37-104 km s$^{-1}$, placing these galaxies systematically below the empirical $v_{max}$-SFR relations for more evolved galaxies of similar SFR. Given the short lifetimes of the SLSN massive progenitors, we argue that these outflows must originate from preceding stellar wind episodes. Assuming a constant-velocity outflow over 3 Myr and spherical symmetry, we infer wind masses M$_{wind}$ = (0.02-1.0) $\times 10^6$ Msun and mass-outflow rates $\dot{M}_{wind} = 0.01-0.33$ Msun yr$^{-1}$, corresponding to mass-loading factors $η<1$. These results indicate that, during the first few Myr of a burst, stellar winds and radiation pressure alone drive slow and weak outflows in low-mass systems, prior to the onset of dominant supernova feedback. Our work provides one of the first empirical constraints on early feedback phases relevant for high-redshift galaxies, and for time-dependent implementations of stellar feedback in galaxy formation simulations.
Cosmological parameter inference from galaxy clustering relies critically on accurate estimates of the covariance and precision matrices. These are often obtained from a limited number of mock catalogs, introducing noise and bias in the precision matrix when the data-vector dimension becomes comparable to the number of available realizations. We present the first application of the Rotational Invariant Estimator (RIE) to the large-scale clustering of galaxies, benchmarking it against the standard sample covariance and the non-linear shrinkage estimator NERCOME for both the two-point correlation function (2PCF) and power spectrum. Using controlled synthetic data sets with analytically known covariance matrices, we estimate the covariance with all three methods across a range of mock-to-dimension ratios $q = N/D$ and data-vector sizes $D$. We then perform Bayesian inference with an EFT-based model and quantify each estimator through the Figure of Bias (FoB) and Figure of Merit (FoM). After correction for finite-$N$ effects, the sample covariance recovers unbiased average uncertainty volumes but suffers from growing best-fit scatter and bias at small $q$ due to the Dodelson--Schneider effect. Both NERCOME and RIE substantially reduce these stochastic shifts; however, the uncertainties they assign are probe-dependent. In configuration space, both estimators can yield overly tight constraints, with a bias that grows with $D$. In Fourier space, RIE delivers markedly improved best-fit stability with only mild FoM bias, whereas NERCOME tends to overestimate the constraining power. Among the estimators tested, RIE emerges as the most effective at stabilizing best-fit recovery, particularly in Fourier space, where it closely reproduces the reference posteriors even when the number of mocks barely exceeds the data-vector dimension.
Protostellar outflows display wide-angle winds and collimated jets, the magnetocentrifugal launching of which enables accretion onto the protostar. The majority of the outflow mass is likely ejected or entrained molecular H$_2$, which can now be studied in unprecedented detail with JWST. Using JWST MIRI/MRS observations towards 13 single and 20 multiple Class 0 and I protostars, we investigate the nature and evolution of the H$_2$ wind and jet morphology, mass outflow rate, and velocity and temperature structure. We construct line flux and velocity maps of the H$_2$ S(1) and S(7) lines as well as the sub-mm CO traced by ALMA. Low-$J$ ($J\le4$) H$_2$ transitions trace extended wide-angle, low-velocity (0-20 km s$^{-1}$) winds within the contours of the low-velocity ($< 30$ km s$^{-1}$) sub-mm CO emission, while high-$J$ ($J >5$) transitions are associated with shocks and knots. In Class 0 sources with a known high-velocity ($> 30$ km s$^{-1}$) molecular CO or SiO jet, higher H$_2$ velocities are found along the jet axis. The opening angle of the wind traced by the H$_2$ S(1) line broadens from $\sim20^\circ$ to $\sim90^\circ$ through the Class 0 to Class I stage. Near the base of each blue-shifted outflow lobe, we extract representative spectra, where rotation diagram fitting of the H$_2$ lines is combined with the outflow width and H$_2$ line velocity to measure the mass-loss rates. The rotation diagrams show a warm $\sim 600$ K, component with two orders of magnitude more mass than the hot, 1500-3000 K component. The H$_2$ outflow mass-loss rates decline by two orders of magnitude from the Class 0 to Class II stage and are correlated with bolometric luminosity. The declining warm H$_2$ mass loss rates and increasing opening angles from the Class 0 to I stages, and the absence of H$_2$ jets in the Class I sources, are consistent with the predictions of MHD disk wind models.
It is difficult to implement numerical simulations on a region extending from the vicinity of a black hole to the Bondi radius. Most previous numerical simulations have primarily concentrated on the region close to the black hole. They found that strong winds can be generated in the hot accretion flows near the black hole, and that radiative cooling significantly affects the strength of these winds. However, the effects of radiative cooling on the production and properties of winds around the Bondi radius remain unclear. In this paper, we perform two-dimensional magnetohydrodynamic simulations to study the impact of radiative cooling on the dynamics and wind production in hot accretion flows around the sub-Bondi radius. As the increase of mass accretion rate, radiative cooling gradually becomes strong, resulting in a reduction in the thickness of the accretion disk (defined as the accretion flows within the density scale height). Based on the Høiland criterion, we find that within the accretion disk, the region of convective stability accounts for \sim 55\% to 62\%, and therefore the accretion flows are marginally stable in convective stability. In the runs with weak radiative cooling, the winds play a significant role in the inward decrease of the mass inflow rate. In the runs with strong radiative cooling, the mass outflow rate of winds is significantly reduced and then the inward decrease of the mass inflow rate is mainly attributed to turbulence driven by magnetorotational instability and convection. Radiative cooling has the potential to suppress accretion processes and reduce the power of winds.
We study the optical depth to reionization in a cosmological setting that includes backreaction from matter inhomogeneities, using the Buchert averaging formalism. We construct a spacetime model consisting of multiple inhomogeneous domains, hereafter referred to as the backreaction model, characterized by a set of parameters. We first examine how these parameters influence the computation of the optical depth to reionization, $τ_{reion}$. Next, we carry out a Markov Chain Monte Carlo (MCMC) analysis based on the PantheonPlus+SH0ES Type Ia supernova sample to infer the best-fit values of the model parameters, and then use these to evaluate $τ_{reion}$. We obtain $τ_{reion} = 0.0581^{+0.0105}_{-0.0096}$ (68$\%$ confidence limits). This result indicates that, when PantheonPlus+SH0ES data are used to constrain the model parameters, our backreaction model yields a value of $τ_{reion}$ that aligns more closely with observational estimates than the value predicted by the standard cosmological model. We further demonstrate that the backreaction model leads to a modest reduction of the Hubble tension, while avoiding the need for exotic or non-standard physics.
In the past few decades, large universal structures have been found that challenge the homogeneity and isotropy expected in standard cosmological models. This study examines burst clustering in both galactic hemispheres using a recently developed methodology, using spheres in 3D space for testing regularities. Using our new method in both hemisphere we find only one deviation from isotropy. A small one in the Southern and a huge one in the Northern hemisphere. This itself suggests that the two deviations do not likely to come from statistical fluctuation. The northern huge group contains app. 125 gamma-ray bursts (GRBs) corresponds with the so - called Hercules-Corona Borealis Great Wall. The southern group contains 4-5 GRBs locating very close to each other. Two of them (GRB050822 and GRB050318) are close not just in redshift and the angular position but they are very close in observing time (5 months). We concluded that the third important result of this work that this method could not find other overdensity deviation from homogeneity in the GRBs spatial distribution. We have shown that the large-scale density increase in the spatial distribution of gamma-ray bursts does not necessarily violate the cosmological principle.