2026-11-06 | | Total: 106
We present BeAR, the Bern Atmospheric Retrieval code and successor of Helios-r, and its new high-resolution cross-correlation framework for emission and transmission spectra. BeAR performs joint retrievals by adding the log-likelihoods of high- and low-resolution data computed from the same model atmosphere. The framework implements established cross-correlation log-likelihoods, with the line-depth scaling either fitted analytically or sampled as a free parameter, their per-pixel generalization, fast model filtering, and model re-injection adapted to this filtering. The forward model and likelihood run entirely on a graphics processing unit (GPU). One forward model takes about 1 ms on an NVIDIA GH200, allowing millions of models per minute on a GPU cluster. In separate retrievals of Gemini South/IGRINS and JWST/NIRSpec dayside spectra of the hot Jupiter WASP-77Ab, BeAR reproduces published abundances, C/O ratios, and metallicities to within 1.5 sigma, the preference for 13CO, and the thermal structure and line contrast to within 1 sigma. The joint retrieval confirms a sub-solar metallicity and a near-solar C/O ratio. A full IGRINS retrieval with 17 free parameters takes 27 to 42 minutes on a single NVIDIA GeForce RTX 5080 or about 14 minutes on eight RTX 4090 GPUs. Injection--recovery tests show that the posteriors underestimate the uncertainties of abundance ratios by about a factor of two.
Observations and theory increasingly indicate that protoplanetary disks are weakly turbulent, implying that magnetized disk winds may dominate their angular momentum transport and long-term evolution. Yet dust evolution in wind-dominated disks remains poorly understood, despite its importance for planet formation. In particular, whether wind-dominated disk evolution can create regions with enhanced dust-to-gas mass ratios that enable planetesimal formation via the streaming instability is unclear. We build a parametric hybrid disk evolution model that includes turbulent accretion, wind-driven accretion, and wind-driven mass loss to investigate dust evolution in magnetized protoplanetary disks. We employ this model to carry out a parameter study in which the initial total accretion rate, the wind-driven mass loss rate, and the turbulence strength are varied. We find that although wind-dominated disks are more favorable to dust growth and concentration than turbulence-dominated disks, conditions suitable for the streaming instability arise only within the inner few au and only under the most extreme conditions: very weak turbulence ($α_\text{turb} = 10^{-5}$) where magnetized winds drive most of the accretion ($\dot{M}_\text{acc,tot} \simeq \dot{M}_\text{wind,acc}$), and wind-driven mass loss comparable to the accretion rate ($\dot{M}_\text{wind,outflow} \simeq \dot{M}_\text{acc,wind}$). Our results therefore imply that, beyond the inner few au, the parametric wind-driven evolution models considered in the literature do not by themselves produce conditions suitable for planetesimal formation via the streaming instability, motivating more physically complete models in which spontaneous, long-lived dust traps can arise.
The dispersion measures (DMs) of fast radio bursts (FRBs) trace ionised baryons along the line of sight, probing baryonic feedback. Two complementary observables can be constructed from an FRB catalogue: the DM-redshift (Macquart) relation of host-identified FRBs, and the angular power spectrum of the DM field of the full catalogue. Their joint analysis requires their cross-covariance. We derive it within the halo model and show it is a super-sample effect: only the monopole of an individual DM correlates with the angular power spectrum, yielding an exact, model-independent bound on the correlation coefficient. We validate the calculation against log-normal forward simulations of $10^4$ FRBs. For individual DMs, the correlation with the bandpowers is at the per cent level. Neglecting the cross-covariance changes the width of the $χ^2$ distribution by only 0.8 per cent, so the two observables can be treated as effectively independent for wide-area surveys. In a joint forecast with an LSST-like weak lensing survey, adding the Macquart relation tightens constraints on the physical baryon density by 32 per cent and on the mean host contribution by a factor of 4.1, while improving baryonification parameters by only about 5 per cent, indirectly through the host contribution. Finally, we show that the redshift distribution of FRBs without identified hosts, which enters every DM power spectrum, can be reconstructed by inverting their DMs with the jointly fitted mean relation without biasing the inferred parameters, whereas a 20 per cent error in the mean relation biases the DM power spectrum by up to 2$σ$ per band for $5\times 10^5$ FRBs.
Thermal radiation transport (TRT) is computationally challenging because the photon specific intensity depends simultaneously on space, angle, and frequency, producing a high-dimensional state space whose cost grows rapidly with resolution. Tensor-network methods, such as the tensor train (TT) decomposition, can mitigate this curse of dimensionality when the solution admits low-rank structure, but localized nonseparable structure can drive substantial rank growth. Adaptive mesh refinement (AMR) provides a complementary means of making high spatial resolution tractable by concentrating resolution on features of interest that would otherwise be prohibitively expensive to resolve. When tensor representations are constructed locally on AMR blocks, this decomposition offers an additional advantage: nonseparable structure can remain confined to the blocks in which it occurs rather than driving rank growth throughout the domain. This paper introduces a framework for combining block-adaptive mesh refinement with block-local tensor-train representations of the photon specific intensity. The approach couples spatial adaptivity with localized low-rank compression, allowing fine-scale features to be resolved while limiting the impact of locally complex structure on tensor ranks elsewhere in the domain. We develop conservative ghost-cell transfer, coarse-fine flux correction, and refinement criteria that balance spatial discretization and tensor-rounding errors. Numerical experiments demonstrate that the method preserves low-rank structure across much of the domain while resolving localized features with AMR, achieving compression factors exceeding $\sim$10$^3$. The resulting framework provides a practical foundation for applying tensor-network methods to high-dimensional TRT on the multiscale, adaptively refined meshes encountered in radiation hydrodynamics simulations.
We present an analysis of JADES-GS-BD-5, an extremely cold Y dwarf identified in the JADES/GOODS-South field. Using multi-epoch JWST/NIRCam imaging spanning a $\simeq2$-yr baseline, we calibrate the JWST astrometry and obtain the first robust parallax measurement for a brown dwarf discovered in a JWST extragalactic deep field. We measure a proper motion of $μ_{\rm tot}=49.8^{+0.3}_{-0.3}$ mas yr$^{-1}$ and a parallax of $\varpi=15.7^{+1.4}_{-1.3}$ mas, corresponding to a geometric distance of $d=64^{+6}_{-5}$ pc, $\sim3\times$ further than that of the previously known most distant Y dwarfs with parallax-measured distances. We hence determine a tangential velocity of $v_{\rm tan}=15.0^{+1.4}_{-1.3}$ km s$^{-1}$, consistent with a Galactic thin-disk object. We find no evidence for binary acceleration, significant photometric variability, or a plausible host. The $J$-band spectral morphology of JADES-GS-BD-5 is consistent with existing Y-dwarf templates, while its extremely faint luminosity and overall near-infrared spectral profile confirm that it is one of the coldest dwarfs discovered to date. We therefore adopt a classification of $>$Y1, consistent with a NIRCam colour-based proxy of cY1.5. The robust parallax-determined distance allows us to determine its bolometric luminosity, $\log(L_{\rm bol}/L_{\odot})=-7.08^{+0.09}_{-0.09}$. Evolutionary modelling places JADES-GS-BD-5 in the planetary-mass regime and makes it a rare benchmark for the study of cold giant-planet atmospheres. Current atmospheric models fail to reproduce its SED in detail, pointing to the importance of disequilibrium chemistry and vertical mixing. Finally, we also report additional proper motion measurements for eight of the other nine known JADES Y-dwarf candidates, highlighting JWST extragalactic surveys as a powerful route to building a statistically useful, parallax-calibrated sample of Y dwarfs.
We investigate how spatial resolution affects emission-line diagnostics in the interstellar medium (ISM) of the Small and Large Magellanic Clouds (SMC and LMC) using SDSS Local Volume Mapper (LVM) data obtained at and degraded to 19, 42, 170 pc (for the SMC) and 17, 68, 272 pc (for the LMC), respectively. We compare the distributions of strong emission-line ratios through 1D and 2D Kolmogorov Smirnov tests. We deduced that diagnostics involving ionization and density sensitive ratios such as [OIII]λ5007/H\b{eta}, [SII]λλ6717+6731/Hα and [NII]λ6584/Hα are the most strongly impacted by resolution loss in this order. The SMC shows significant deviations even from intermediate scales (42 pc), whereas for the LMC, the coarsest resolution (272 pc) remains relatively stable, indicating that the SMC is intrinsically more affected by spatial mixing. Resolution degradation produces a marked homogenization effect: extreme values present at the finest scales are progressively diluted, causing diagnostic diagrams to converge toward their mean loci. Diffuse ionized gas is particularly sensitive to this mixing, while Hα selected ionized regions exhibit more robust behavior across resolutions. The sensitivity of diagnostic diagrams involving emission-line ratios, equivalent widths, and kinematic properties varies with both spatial resolution and galaxy, with the SMC generally showing stronger changes than the LMC. Our results demonstrate that sub-50 pc resolution is required to resolve the detailed structure of low-metallicity galaxies, and they highlight at least the need for careful interpretation of line-ratio diagnostics in extragalactic studies where coarse spatial resolution is unavoidable.
In hierarchical galaxy formation, supermassive black hole binaries are a ubiquitous consequence of frequent galaxy mergers, yet subparsec binaries remain difficult to be confirmed. Periodic flux modulations may arise from the orbital motion of black holes in close binaries, while relativistic jet precession can provide an additional source of periodic variability. Such periodic signatures therefore motivate us to identify close binary candidates. Here, we report an optical periodicity of $269^{+13}_{-14}$ days in both Zwicky Transient Facility g and r bands of the blazar PKS 0301-243, with corrected global significances of $2.75σ$ and $2.18σ$, respectively. This signal persists over 7 cycles across a baseline of ~2000 days. We interpret this period as a possible modulation associated with the orbital motion of a supermassive black hole binary. A longer optical periodicity of $889^{+47}_{-46}$ days is identified in Katzman Automatic Imaging Telescope, with corrected global significance of $1.97σ$. The longer periodicity spans approximately 5 cycles over a baseline of ~4500 days, potentially associated with jet precession. Both periodicities are independently identified using multiple methods. A supermassive black hole binary incorporating a relativistic jet and a circumbinary disc could reproduce the observed light curve ensemble. For a total mass of $10^{8.95} M_\odot$ and a mass ratio q = 0.1, the derived milli-parsec separation and a short gravitational wave driven coalescence time, $t_c = 1122^{+144}_{-144}$ years, imply that PKS 0301-243 may be a rare source of nano-hertz gravitational waves and a potential multimessenger laboratory for the terminal stages of supermassive black hole coalescence.
We generalize a previous formulation for the estimation of polarization cross-talk induced by atmospheric seeing to the case of image jitter that must be described as a non-stationary, and generally non-isotropic, random process. We demonstrate how the newly derived formalism converges towards the previous one in the limit of stationary random image jitter. In the case of non-stationary processes, we compare the predictions of the two formalisms, and illustrate the salient differences. For the sake of demonstration, we adopt for our analysis a specific example of power-spectral density derived from a circulant auto-correlation function, which describes a class of stationary random processes periodic over the length of the time series, and which we successively modified by the insertion of off-diagonal power spectral features at different frequencies in order to test the extended modeling framework for the estimation of polarization cross-talk.
Space-based astrometric missions such as Hipparcos and Gaia rely on one-dimensional scanning strategies that couple the measured signals with both observation times and scanning geometry. This coupling can introduce spurious periodicities and complicate their interpretation. We present a non-parametric extension of the phase-distance correlation (PDC) periodogram that quantifies and reduces dependence on scan geometry at the period-search stage. Using Hipparcos data, we identified recurrent periodicities associated with characteristic scanning-law timescales and showed that they are strongly suppressed by the scan-angle-controlled statistic. We then applied the method to Gaia Data Release 4 (DR4) prerelease epoch astrometry, where a recurrent ~54-day feature is suppressed across multiple sources while independently known orbital signals remain detectable. The method provides a diagnostic of scan-angle dependence and a way to mitigate its impact on period searches. These results are particularly relevant for Gaia DR4, where public epoch astrometry will enable large-scale searches for companions ranging from exoplanets to dormant black holes.
Previous studies combined rotation periods derived from Kepler photometry, stellar parameters from previous Gaia Data Releases, and planet hosts from the NASA Exoplanet Archive to conclude that statistically planet-hosting stars rotate slower than stars without detected planets. Here, we revisit these studies by using similar Kepler rotational period sample, but with an updated Gaia DR3-based stellar parameter catalogue and an updated exoplanet catalogue. We apply minimal sample cuts in evolutionary state, metallicity, and distance, and then model the stellar rotation period as a function of effective temperature alone. For each star, we derive the ratio of its observed period to the temperature-expected period, and compare the ratio distributions for the planet-host and non-host stars. For periods shorter than 20 d, the planet-host distribution of this ratio is shifted towards longer periods; the period-ratio median is ~15% larger for hosts than for non-hosts (p_KS = 2.4 x 10^-4). We also examine the recently published TESS All-Sky Rotation Survey (TARS) and find a similar trend. However, the strongly structured period distribution precludes a robust comparison between the host and the non-host stars. Therefore, the updated analysis, using more robust stellar and planetary parameters, recovers the previously reported association, strengthening the case for a genuine connection between stellar rotation and planet occurrence.
The vast majority of cosmological simulations employ periodic boundary conditions to achieve statistical homogeneity and to mimic the infinite spatial extent of our Universe. In a cubic periodic domain the Green function breaks full rotational invariance and introduces a geometry-induced anisotropy in the gravitational interaction. We quantify how this intrinsic anisotropy affects cosmic structure growth and investigate the extent to which it biases the commonly used large-scale structure and halo-based observables in dark-matter-only simulations. We analyse an ensemble of 498 cosmological dark-matter-only $N$-body simulations in cubic $T^3$ topological manifolds, evolved to $z=0$, spanning box sizes from $L=6.25\,\mathrm{Mpc}$ to $L=100\,\mathrm{Mpc}$ at fixed mass and force resolution. We compare anisotropy-sensitive statistics across the suite, including the matter power spectrum, halo abundance, merger histories, shapes and alignments, and filament orientations. Finite periodic boundaries imprint orientation-dependent signatures aligned with the box geometry. These signatures propagate beyond the largest modes and can induce systematic shifts in both field-level and halo-level statistics. The detectable effects propagate to scales much smaller than the box size. Our results demonstrate that periodic-box anisotropy is a non-negligible source of systematics for small simulation volumes, including those commonly employed in modern hydrodynamical simulations.
Reliable collisional rate coefficients are crucial for interpreting the rotational emission spectra of molecules in interstellar clouds. In this work, we construct a two-dimensional potential energy surface (PES) for the $\mathrm{HC}_5\mathrm{NH}^+$--He system at the CCSD(T)-F12B/cc-pVTZ-F12 level of theory. Rotationally inelastic cross sections are calculated using the quantum close-coupling (CC) method, and the corresponding collisional rate coefficients are obtained by thermal averaging over a Maxwell--Boltzmann distribution. We use these collisional data in a non-local thermodynamic equilibrium (non-LTE) radiative transfer model to investigate the excitation of $\mathrm{HC}_5\mathrm{NH}^+$ in the cold dark cloud TMC-1. A global $χ^2$ minimization analysis of seven observed lines yields best-fitting values of $n(\mathrm{H}_2)=1.63\times10^{5}$\,cm$^{-3}$ and $N(\mathrm{H}_2)=5.09\times10^{11}$\,cm$^{-2}$, with reduced $χ^2_{\min}=0.87$. The modelled line intensities are in good agreement with the observations. These results provide insight into the collisional excitation and radiative properties of $\mathrm{HC}_5\mathrm{NH}^+$ under the physical conditions prevailing in TMC-1.
This study reports the meter-wavelength radio imaging of multiple type-II bursts occurred on 01 October 2011 observed by LEARMONTH radio spectrograph in the frequency range of 180 - 25 MHz, with two coronal shocks from two coronal mass ejections (CMEs) from the nearby active regions lying adjacent to each other. The two CMEs are well observed by STEREO and LASCO FOV. The eruption is associated with the multi-peaked soft X-ray flare erupted from N16W45 (first and second C3 and C6 class flares) and from N09W04 (third one an M Class flare). The first type-II burst starts during the impulsive phase of the first flare and the second type-II burst starts during the peak time of the third flare. The drift rate of the second type-II burst is larger than that of the first type-II burst. The EUV images at AIA 211 reveal the shock-like EUV wave structure along the north-western limb around 09:06 UT, the first type-II appearing around the same time (09:06 UT) confirms the existence of the shock wave and as an outward propagating wave along the disk around 09:55 UT, the second type-II appearing around the same time (09:55 UT) confirms the existence of the shock wave. According to the simultaneous EUV and radio imaging data the radio sources appears along the nose-front of the shock like EUV wave structure for the first type-II burst and for the second type-II burst they are located along the shock flank. The study shows the possibility of two coronal shocks from two CMEs generating the multiple type-II radio bursts. The second event is strongly affected by the projection and scattering effects, the corrections are achieved by applying spherical cone model and 3D Monte Carlo ray-tracing code for turbulent solar environment at fundamental frequency.
Superbubbles are vast hot cavities carved out by the stellar winds and supernovae of clusters of high-mass stars. We investigate G92.4+5.0 (the CTB 102 superbubble), a Perseus Arm superbubble associated with the HII region CTB 102. We use new and archival radio observations (total intensity, rotation measure, and HI spectral line) from the Dominion Radio Astrophysical Observatory Synthesis Telescope, combined with single-antenna data from the Effelsberg Medium Latitude Survey, the Global Magneto-Ionic Medium Survey, and the HI4PI survey to provide data on large structures. The upper region of the superbubble is faint in HI. In order to detect its full extent, we supplement our observations with H-alpha data from the Wisconsin H-Alpha Mapper Survey and the Ziegenbalg Survey, and observe H-alpha emission in the upper region. We measure the CTB 102 superbubble to be at least ~170 pc wide and ~420 pc tall, with an expansion rate of ~12 km/s in the lower regions of the superbubble and ~22.5 km/s in the upper regions. The elongated structure and increasing expansion rate with height are consistent with the evolutionary properties of a superbubble described by the Kompaneets model. We fit the Kompaneets model to the CTB 102 superbubble and obtain an estimated age of 1-6 Myr. This supports a possible scenario in which the stellar cluster currently ionizing CTB 102 is the expansion energy source for the superbubble. Alternatively, if the superbubble expanded under the power of a previous generation of stars, the current CTB 102 stellar cluster may have formed in the walls of the superbubble.
Increased computing power and more efficient codes have allowed cosmological simulations of galaxy formation to grow massively in scale: the analysis of their outputs has consequently become a big data computational challenge. We present Octavius, a high-performance, hybrid-parallel simulation analysis toolkit. Octavius ingests a simulation snapshot with associated halo catalogue, locates galaxies with a 6D friends-of-friends algorithm, then produces an HDF5 catalogue of properties for haloes, subhaloes, and galaxies. The analysis is fast and memory efficient: we compare the performance against the widely used predecessor toolkit, Caesar, and show Octavius delivers order-of-magnitude improvements in speed with a factor of >5 reduction in peak memory. We demonstrate near-linear strong scaling of compute and linear peak memory scaling with the number of halo-assigned particles. Compute times scale with the superlinear computational load of larger haloes; we note the pipeline is presently I/O-limited. The code is written in Python and uses a declarative, plugin-based architecture so as to be easily extensible. Octavius can currently analyse SIMBA, SWIFT-KIARA, SWIFT-EAGLE, SWIFT-COLIBRE, and TNG snapshots; AHF, HBT-HERONS, and SUBFIND halo catalogues are supported. Properties computed include masses, galaxy photometry with dust attenuation, and kinematic, radial, morphological, and virial quantities. Output catalogues are easily accessed through a provided user API. The code is open sourced at https://github.com/jp-duminy/octavius
Double-periodic variables (DPVs) are interacting binaries with an orbital modulation and a longer cycle of uncertain origin. I searched for long-cycle-dependent eclipse-center variations in Magellanic DPVs and tested whether changing disk brightness can produce analogous shifts with fixed conjunctions. I analyzed Optical Gravitational Lensing Experiment (OGLE) photometry for 109 eclipsing-DPV candidates. In separate long-phase intervals, I measured centers of eclipses by profile fitting to local flux minima ($Q$) and by shifting a fixed eclipse template ($T$). A single synthetic semidetached DPV provided a separate numerical experiment. Of 49 curated systems with sufficient timing data, 27 show at least one $Q$ or $T$ modulation under the adopted procedure. Some selections depend on the fitting choice, but ten systems retain a selected $Q$ curve with all tested local functions and windows. Three have both estimators selected in both eclipses, without establishing whether conjunction times change. Using local color interpolation, 14 of 16 systems measured in both eclipses are redder at the deeper, phase-0 eclipse. In the artificial system, fixed conjunctions coexist with large photometric shifts and quadrature-brightness changes, although the large tilted-disk cases do not reproduce the selected observed crossing phases. Long-cycle-dependent photometric eclipse-center variations occur in a substantial fraction of the evaluable, curated sample, not necessarily of all DPVs. Circumstellar brightness changes can produce such shifts, but the single artificial system neither fits an observed DPV nor establishes disk tilt as the cause. Genuine conjunction-time changes remain possible in individual objects.
We present results from a stellar occultation by the asteroid (623) Chimaera on October 21st, 2025, when it passed in front of V magnitude 13.86 star UCAC4 626-004193. Chimaera is the largest asteroid in the Chimaerid asteroid family and will be subject to a spacecraft flyby in the 2030s. Eighteen telescopes were deployed, spanning the United States from Arizona to North Carolina, incrementally spaced by $\sim4.4$\,km with respect to the cross-track of the shadow. In total, sixteen telescopes collected useful data. This multi-chord observation, when compiled into a single sky-plane plot, suggests a profile which is well-described by an elongated elliptical fit, with a semi-major axis $a=26.2\pm0.8$ km and a semi-minor axis $b=16.0\pm1.0$ km. The projected area from ellipse fitting is 1316 km$^2$, and a limb fit is calculated to enclose 1312 km$^2$. Photometry collected in the Gaia G band over a seven-month period coinciding with the occultation corresponds to a geometric albedo $p_G=0.059\pm0.005$; comparing with archival occultation observations and direct observations gives $p_G=0.051\pm0.004$ and $p_V=0.051\pm0.004$. We compare occultation sky-plane projections with 3-D modeling from inversion techniques and demonstrate that the models can accurately predict the orientation of Chimaera at the time of each occultation.
Direct imaging of Earth-like exoplanets requires strong suppression of stellar light at small angular separations, achievable with extreme adaptive optics (XAO). The pyramid wavefront sensor (PWFS) is widely used for its sensitivity but suffers from nonlinearities and limited dynamic range. Modulation extends the range but reduces sensitivity and limits loop speed. We investigate an unmodulated infrared PWFS (IR-umPWFS), where the dynamic range scales with wavelength, removing the need for modulation and enabling faster operation. We adopt a 3-sided PWFS, which eliminates the roof defect of the classical 4-sided design. We assess its feasibility and performance through end-to-end simulations, lab tests, and on-sky observations on the PAPYRUS bench at the Observatoire de Haute-Provence (OHP), and compare it with a visible modulated PWFS (VIS-mPWFS). Simulations show reliable operation, reaching Strehl ratios up to ~60% at 1400 nm for seeing up to 4''. Infrared operation improves linearity, reduces optical-gain variations and modal cross-talk. Laboratory results validate the simulations and confirm the pyramid quality. On sky, stable closed-loop operation was achieved over three nights. Compared with the VIS-mPWFS at 500 Hz, the IR-umPWFS at 1000 Hz showed gains of ~1.63 in Strehl ratio and ~1.26 in encircled energy, and a ~2.03-fold reduction in residual wavefront error. The IR-umPWFS remained robust to tip-tilt offsets up to 3 lambda/D at 1400 nm. It is a strong candidate for future XAO systems: operating in the infrared extends the linear range sufficiently to remove the need for modulation, enabling higher loop speeds no longer limited by the modulation mirror, which contributed to the improved on-sky performance. The OHP turbulence (D/r0 ~ 17 at 1400 nm) is comparable to typical VLT conditions (D/r0 ~ 20), supporting the relevance of these results for VLT/RISTRETTO.
We present a multi-probe observational analysis of a clustering tachyonic Dirac-Born-Infeld~(DBI) dark energy model. Its noncanonical kinetic term gives a subluminal sound speed, so dark energy perturbations can grow on subhorizon scales. Working in a multifluid framework that tracks both the background expansion and linear perturbations, we constrain the model with nested sampling, using Planck 2018 compressed CMB distance priors, DESI DR2 BAO, cosmic chronometers, $fσ_8$ growth data, and three Type Ia supernova compilations: Pantheon$+$, DES-SN5YR, and Union3. The standard cosmological parameters agree with $Λ$CDM within $2σ$, with DBI preferring a slightly lower $H_0$ and higher $Ω_{m0}$, with the largest shift being $\simeq1.8σ$ in $H_0$ for Pantheon$+$. The Bayesian evidence favours DBI in every compilation, moderately for Pantheon$+$ ($Δ\ln\mathcal{Z}\simeq3.4$) and inconclusively for DES-SN5YR and Union3 ($Δ\ln\mathcal{Z}\simeq0.6$-$0.7$). The field thaws naturally. Dark energy perturbations stay several orders of magnitude below matter perturbations, yet the noncanonical sound speed still alters the Hubble friction term, raising $S_8$ by only $\approx0.5σ$ relative to $Λ$CDM. Four independent diagnostics, namely the CPL parametrization, $Om(z)$, the statefinder hierarchy, and the growth index, all point to genuinely evolving dark energy rather than a cosmological constant. The thermal history is left intact, with negligible early dark energy. However, one limitation remains. Because $w_φ=\dotφ^2-1\geq-1$, the phantom regime is out of reach, which is in line with the no-go theorem for a single minimally coupled scalar field, so the model cannot follow the DESI DR2 preference for an effective equation of state crossing $w=-1$.
We study global atmosphere models of AF supergiants. In agreement with previous studies, we found super-Eddington layers in the photospheres of these stars, which lead to density inversions. Our tests showed that the super-Eddington layers are caused by hydrogen ionization. The analysis of radial perturbations has not revealed any instability of such layers, but we show that super-Eddington layers in an optically thick environment are prone to instabilities due to nonradial perturbations. We predict the wind parameters of AF supergiants from our global wind code with a comoving-frame radiative force. Despite the existence of density inversions, the predicted mass-loss rates are significantly lower than those estimated from observations. Theoretical values approach the observational estimates only for the most luminous supergiants, which are close to the Eddington limit due to light scattering on free electrons. Part of the discrepancy between the mass-loss rates could be attributed to the influence of the clumping on the observational properties. Low mass-loss rates imply that the winds of AF supergiants may be strongly influenced by the magnetic fields, suggesting an explanation for the observed spectral variability of these stars.
Large spectroscopic surveys have measured atmospheric parameters for millions of Milky Way stars, but differences in instruments, wavelength ranges, analysis methods, and assumptions introduce survey-dependent systematics that complicate their joint use. We aim to place effective temperature, surface gravity, and metallicity estimates from major spectroscopic surveys onto a common, externally validated scale, and to provide a homogenised catalogue for large-scale Galactic studies. We construct the calibration hierarchically. We first define an internal reference scale from the best-populated and best-overlapped surveys (APOGEE, GALAH, and LAMOST LRS), and model survey-to-reference corrections and their uncertainties with machine-learning regression. The scale is then extended to LAMOST MRS, Gaia RVS, RAVE, DESI, Gaia-ESO, and SEGUE, and anchored to high-quality measurements from the PASTEL database. The calibration reduces survey-to-survey offsets and improves consistency in independent overlap samples: independent calibration of GALAH and LAMOST LRS temperatures onto the APOGEE scale reduces their median offset from 44.5 K to 6.2 K. The final catalogue contains about 10.13 million stars, including 6.10 million in the recommended high-quality subset and 3.69 million in the metallicity-optimised subset. For the high-quality subset, Gaia benchmark stars give median residuals of 30 K in $T_{\rm eff}$, 0.03 dex in $\log g$, and $-0.02$ dex in [Fe/H]. The [Fe/H]-optimised subset gives median residuals of 0.003 dex against Gaia benchmark stars, 0.05 dex for globular clusters, and $-0.02$ dex for open clusters. As an example application, we confirm the presence of metal-poor stars on thin-disc-like orbits, including stars with [Fe/H] $< -2$, showing that the catalogue can be used to study rare Galactic populations on a consistent multi-survey scale.
The atmospheric boundary layer (ABL) presents turbulence regimes that diverge from the classical Kolmogorov 1941 (K41) power-law spectrum widely assumed in optical atmospheric profiling, yet turbulence characterization for astronomical applications remains predominantly based on K41 assumptions. Using high-frequency sonic anemometer data from site-testing campaigns at the Giant Magellan Telescope (GMT) and Thirty Meter Telescope (TMT) sites in Chile, we apply detrended fluctuation analysis (DFA)---a technique from fractal time-series analysis---to identify the inertial sub-range and evaluate statistical isotropy in the ABL, comparing our approach with the Fourier-based algorithm for robust identification of the inertial subrange (ARIIS). Strict Kolmogorov conditions---simultaneous isotropy and $k^{-5/3}$ scaling across all fields---are satisfied in fewer than 7 per cent of observations. Temperature fields systematically decouple from velocity components under stable stratification, exhibiting Bolgiano-Obukhov scaling ($k^{-7/5}$) during nighttime hours, indicating active rather than passive scalar behaviour. DFA-detected regime transitions correlate with characteristic turbulence scales, providing objective, model-free boundaries for the inertial subrange. Comparative $C_n^2$ estimation reveals that structure function methods assuming K41 scaling diverge significantly from flux-based parameterizations under stable conditions. The rarity of true Kolmogorov conditions at these sites necessitates stratification-aware methodologies for accurate optical turbulence characterization at next-generation extremely large telescope facilities.
Magnetic-field studies of star-forming regions are complicated by overlapping structures along the line of sight. We investigate the large-scale magnetic field toward NGC~1333 in Perseus, accounting for foreground polarization through stellar distances and optical polarimetry. Our wide-field $R$-band survey covers nearly one square degree and provides polarization measurements for approximately 840 stars, with distances from Gaia DR3 and StarHorse. We used polarization-distance trends and Stokes parameters to identify line-of-sight components and estimate foreground polarization. After subtracting it in Stokes space, we characterized the field using circular statistics and the angular dispersion function (ADF), comparing it with molecular-line and complementary polarimetric observations. The data reveal a polarized foreground at 160-300 pc and a cloud component beyond 300 pc. After correction, the plane-of-sky field rotates smoothly from $126.3°\pm0.8°$ in the northwest to $162.2°\pm0.6°$ in the southeast, remaining well ordered in both regions. The ADF gives turbulent-to-ordered field ratios of $0.054^{+0.025}_{-0.022}$ and $0.095^{+0.021}_{-0.047}$, corresponding to plane-of-sky strengths of $59^{+44}_{-23}$ and $46^{+50}_{-13}~μ\mathrm{G}$. Residual foreground fluctuations may affect these estimates. CO data show multiple velocity components with spatially varying contributions. Corrected optical polarization agrees with the large-scale orientation traced by Planck, while near-infrared and JCMT/POL-2 data reveal localized deviations in denser structures. Our results indicate a coherent cloud-scale magnetic field with a smoothly changing orientation, alongside localized variations in dense structures. They demonstrate the value of combining stellar distances and wide-field optical polarimetry to disentangle foreground contributions.
The Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) is a radio interferometer array designed for surveying fast radio bursts (FRBs) in the nearby Universe with a wide instantaneous field of view (FoV) of $\approx$6000 deg$^2$ and precisely localizing them with sub-arcsecond angular resolution. The first phase of BURSTT, BURSTT-256, consists of a main antenna array comprising 256 log-periodic dipole antennas (LPDAs) observing at 400-800 MHz for real-time beamforming and FRB search, currently with 256 beams covering $\approx$ 3200 deg$^2$ FoV, together with several smaller outrigger stations with 16-64 antennas and separated by $\mathcal{O}(10^2 - 10^3)$ km baselines for joint observations to localize FRBs using very-long baseline interferometry. This work presents the design, deployment, and commissioning of the BURSTT-256 main array, including site selection, radio-frequency (RF) background characterization, receiver antenna array and RF front-end, as well as real-time, factorized digital beamforming backend that currently generates 256 beams. The system's on-sky performance, sensitivity, and beamforming fidelity have been calibrated through observations of bright celestial calibrators. The antenna array in the initial phase using commercial LPDAs has a system equivalent flux density of $\approx$11 kJy, and is expected to achieve $\approx$4 kJy using custom-design LPDAs after optimization in late 2026. BURSTT-256 serves as a critical pathfinder for future expansion of the array to 2048 elements and beyond, operating at 300-700 MHz with $\approx$6000 deg$^2$ FoV. It establishes the framework required to discover and localize nearby FRBs with sub-arcsecond precision. By substantially increasing the number of well-localized nearby FRBs, BURSTT will facilitate multi-wavelength and multi-messenger counterpart identification, helping to unveil the elusive origin of FRBs.
Flexion, the second-order gravitational lensing distortion, probes small-scale mass distribution and complement the cosmic shear measurement. However, its detection is highly challenging due to the faint signal and the contamination arising from the intrinsic galaxy morphology. In this Letter, we present GradFlex, a new method that extracts flexion from the gradient field of galaxy surface brightness distributions without assuming a specific point-spread function model or galaxy morphology. Using surface-brightness conservation, we relate the pre- and post-lensing gradient fields and derive flexion estimators from ensemble-averaged statistics of the corresponding spins. Tests on numerical simulations across diverse galaxy morphologies and multiple noise realizations demonstrate its reliability under realistic observing conditions. Combined with shear measurements, it could tighten constraints on cluster density profiles and enable more detailed studies of substructure than shear alone.