Mapping the local and distant Universe is key to our understanding of it. For decades, the Sloan Digital Sky Survey (SDSS) has made a concerted effort to map millions of celestial objects to constrain the physical processes that govern our Universe. The most recent and fifth generation of SDSS (SDSS-V) is organized into three scientific ``mappers". Milky Way Mapper (MWM) that aims to chart the various components of the Milky Way and constrain its formation and assembly, Black Hole Mapper (BHM), which focuses on understanding supermassive black holes in distant galaxies across the Universe, and Local Volume Mapper (LVM), which uses integral field spectroscopy to map the ionized interstellar medium in the local group. This paper describes and outlines the scope and content for the nineteenth data release (DR19) of SDSS and the most substantial to date in SDSS-V. DR19 is the first to contain data from all three mappers. Additionally, we also describe nine value added catalogs (VACs) that enhance the science that can be conducted with the SDSS-V data. Finally, we discuss how to access SDSS DR19 and provide illustrative examples and tutorials.

Cometary activity from interstellar objects provides a unique window into the environs of other stellar systems. We report blue-sensitive integral field unit spectroscopy of the interstellar object 3I/ATLAS from the Keck-II-mounted Keck Cosmic Web Imager on August 24, 2025 UT. We confirm previously reported CN and Ni outgassing, and present, for the first time, the radial profiles of Ni and CN emission in 3I/ATLAS. We find a characteristic $e$-folding radius of $593.7\pm14.8$ km for Ni and $841.0\pm15.4$ km for CN; this suggests that the Ni emission is more centrally concentrated in the nucleus of the comet and favors hypotheses involving easily dissociated species such as metal carbonyls or metal-polycyclic-aromatic-hydrocarbon molecules. Additional integral field spectroscopy after perihelion will offer a continued opportunity to determine the evolution of the radial distributions of species in interstellar comet 3I/ATLAS.

Stellar surface inhomogeneities such as spots and faculae introduce Doppler variations that challenge exoplanet detection via the radial velocity method. While their impact on disc-integrated spectra is well established, detailed studies of the underlying local line profiles have so far been limited to the Sun. We present an observational campaign targeting the active star WASP-85 A during transits of its hot Jupiter companion. The transits span two stellar rotation periods, allowing us to probe the evolution of active regions. From ground-based photometry we identify seven active regions, six containing dark spots. Using simultaneous ESPRESSO transit spectroscopy, we spatially resolve these regions on the stellar surface by using the planet as a probe. We detect significant bisector shape changes, line broadening, and net redshifts during spot occultations, with velocity shifts of 108-333 m/s (mean uncertainty 50 m/s). The observed broadening is consistent with the Zeeman effect, implying magnetic field strengths (Stokes $I$) $B$ = 2.7-4.4 kG (mean uncertainty 0.6 kG), comparable to solar umbrae. Combined with our photometric spot model, this yields lower limits to the disc-integrated field $Bf = 16 \pm 3$ G and $61 \pm 9$ G for the two hemispheres probed -- at least three times higher than Sun-as-a-star values. We also measure centre-to-limb variations in FWHM, line depth, equivalent width, and convective blueshift, which broadly agree with solar observations and 3D MHD models. This work demonstrates a new way to characterise the surfaces of exoplanet host stars, paving the way for future analyses incorporating synthetic line profiles from 3D MHD simulations.

We present observations of comet 67P/Churyumov-Gerasimenko during its 2021/22 apparition, aiming to investigate its dust and gas environment and compare the results with those obtained in 2015/16 using the same telescope. Quasi-simultaneous photometric, spectroscopic, and polarimetric observations were carried out at the 6-m BTA SAO telescope. The comet was observed on 6 October 2021, 31 days before perihelion, with \textit{g}-SDSS and \textit{r}-SDSS filters, and on 6 February 2022, 96 days after perihelion, using narrowband cometary filters: BC ($\lambda4450/62$~Å), RC ($\lambda6839/96$~Å), and CN ($\lambda3870/58$~Å). These were complemented by images from the 2-m Liverpool Telescope (La Palma). On 6 October 2021, a sunward jet and long dust tail were detected. By 6 February 2022, the dust coma morphology had changed noticeably, revealing a bright sunward neckline structure superimposed on the projected dust tail, along with two jets at position angles of 133$^{\circ}$ and 193$^{\circ}$. Spectra showed strong CN emission, with relatively weak C$_2$, C$_3$ and NH$_2$ emissions. The dust production rate $Af\rho$ did not exceed 200~cm (uncorrected for phase angle) in both epochs. An unusual CN coma morphology was observed, with evidence of an additional CN source associated with dust jets. Geometric modeling of the jets' dynamics indicated an active area at latitude $-70^{\circ} \pm 4^{\circ}$ with a jet opening angle of $20^{\circ} \pm 6^{\circ}$ on 6 October 2021, and two active areas at latitudes $-58^{\circ} \pm 5^{\circ}$ and $-53^{\circ} \pm 10^{\circ}$, separated by longitude $150^{\circ} \pm 20^{\circ}$, producing the observed jets on 6 February 2022. The average particle velocity in the jets was about $0.32 \pm 0.04$~km~s$^{-1}$.

We present new VLTI/GRAVITY astrometry and updated orbit fits for the directly imaged companions YSES 1 b and HR 2562 B, substellar objects straddling the planet-brown dwarf boundary. Using high-precision astrometry, radial velocity (RV) data, and proper motions, we derive revised orbital parameters with orbitize! arXiv:1910.01756. For YSES 1 b, the inclusion of GRAVITY astrometry and a relative radial velocity measurement from arXiv:2409.16660 overcomes the traditional challenge of constraining eccentricities for distant companions, enabling the first orbit fit and yielding a constrained eccentricity of 0.44 (0.20). This represents the first full orbit fit for the system. Additionally, we calculate a median line-of-sight stellar obliquity of 12 (+11, -8) degrees, providing further insight into the system's dynamical architecture. For HR 2562 B, our analysis agrees with arXiv:2302.04893, confirming a low-eccentricity orbit (0.34 (0.20)) and an inclination of 87 (1) degrees. We find HR 2562 B's orbit to be nearly coplanar with the debris disk, with a mutual inclination of 3.7 (0.3) degrees. For both YSES 1 b and HR 2562 B the lower eccentricities favor an in situ formation scenario over extreme scattering or cloud fragmentation.

We report initial observations aimed at the characterization of a third interstellar object. This object, 3I/ATLAS or C/2025 N1 (ATLAS), was discovered on 2025 July 1 UT and has an orbital eccentricity of $e\sim6.1$, perihelion of $q\sim 1.36$ au, inclination of $\sim175^\circ$, and hyperbolic velocity of $V_\infty\sim 58$ km s$^{-1}$. We report deep stacked images obtained using the Canada-France-Hawaii Telescope and the Very Large Telescope that resolve a compact coma. Using images obtained from several smaller ground-based telescopes, we find minimal light curve variation for the object over a $\sim4$ day time span. The visible/near-infrared spectral slope of the object is 17.1$\pm$0.2 %/100 nm, comparable to other interstellar objects and primitive solar system small bodies (comets and D-type asteroids). 3I/ATLAS will be observable through early September 2025, then unobservable by Earth-based observatories near perihelion due to low solar elongation. It will be observable again from the ground in late November 2025. Although this limitation unfortunately prohibits detailed Earth-based observations at perihelion when the activity of 3I/ATLAS is likely to peak, spacecraft at Mars could be used to make valuable observations at this time.

LTT-9779 b is an ultra-hot Neptune (Rp ~ 4.7 Re, Mp ~ 29 Me) orbiting its Sun-like host star in just 19 hours, placing it deep within the "hot Neptune desert," where Neptunian planets are seldom found. We present new JWST NIRSpec G395H phase-curve observations that probe its atmospheric composition in unprecedented detail. At near-infrared wavelengths, which penetrate the high-altitude clouds inferred from previous NIRISS/SOSS spectra, thermal emission reveals a carbon-rich atmosphere with opacity dominated by carbon monoxide (CO) and carbon dioxide (CO2). Both species are detected at all orbital phases, with retrieved mixing ratios of 10^-1 for CO and 10^-4 for CO2, indicating a globally well-mixed reservoir of carbon-bearing gases. We also moderately detect water vapor (H2O) and tentatively detect sulfur dioxide (SO2), providing insight into its chemistry and possible photochemical production under intense stellar irradiation. From these detections we infer a carbon-to-oxygen ratio near unity (C/O ~ 1) and a metallicity exceeding 500X Solar, consistent with equilibrium chemistry predictions for high-temperature atmospheres. This enrichment raises the mean molecular weight, reducing atmospheric escape, and likely helps LTT-9779 b retain a substantial atmosphere despite extreme irradiation. Our findings show that LTT-9779 b survives where few planets can, maintaining a carbon-rich atmosphere in a region where hot Neptune-class worlds are expected to evaporate. This makes LTT-9779 b a valuable laboratory for studying atmospheric escape and chemical processes under extreme conditions, offering new insight into the survival of planets in the hot Neptune desert.

We apply the exponential operator method to derive the propagator for a fermion immersed within a rigidly rotating environment with cylindrical geometry. Given that the rotation axis provides a preferred direction, Lorentz symmetry is lost and the general solution is not translationally invariant in the radial coordinate. However, under the approximation that the fermion is completely dragged by the vortical motion, valid for large angular velocities, translation invariance is recovered. The propagator can then be written in momentum space. The result is suited to be used applying ordinary Feynman rules for perturbative calculations in momentum space.

This is the second paper in the HOWLS (higher-order weak lensing statistics) series exploring the usage of non-Gaussian statistics for cosmology inference within \textit{Euclid}. With respect to our first paper, we develop a full tomographic analysis based on realistic photometric redshifts which allows us to derive Fisher forecasts in the ($\sigma_8$, $w_0$) plane for a \textit{Euclid}-like data release 1 (DR1) setup. We find that the 5 higher-order statistics (HOSs) that satisfy the Gaussian likelihood assumption of the Fisher formalism (1-point probability distribution function, $\ell$1-norm, peak counts, Minkowski functionals, and Betti numbers) each outperform the shear 2-point correlation functions by a factor $2.5$ on the $w_0$ forecasts, with only marginal improvement when used in combination with 2-point estimators, suggesting that every HOS is able to retrieve both the non-Gaussian and Gaussian information of the matter density field. The similar performance of the different estimators\inlinecomment{, with a slight preference for Minkowski functionals and 1-point probability distribution function,} is explained by a homogeneous use of multi-scale and tomographic information, optimized to lower computational costs. These results hold for the $3$ mass mapping techniques of the \textit{Euclid} pipeline: aperture mass, Kaiser--Squires, and Kaiser--Squires plus, and are unaffected by the application of realistic star masks. Finally, we explore the use of HOSs with the Bernardeau--Nishimichi--Taruya (BNT) nulling scheme approach, finding promising results towards applying physical scale cuts to HOSs.

We report VLT spectroscopy of the interstellar comet 3I/ATLAS (C/2025 N1) from $r_{\rm h}\!\simeq\!4.4$ to $2.85$ au using X-shooter (300-550 nm, $R\!\simeq\!3000$) and UVES (optical, $R\!\simeq\!35k-80k$). The coma is dust-dominated with a fairly constant red optical continuum slope ($\sim$21-22\%/1000Å). At $r_{\rm h}\!\simeq\!3.17$ au we derive $3\sigma$ limits of Q({\rm OH})<7.76\times10^{23}\ {\rm s^{-1}}, but find no indications for [O I], C$_2$, C$_3$ or NH$_2$. We report detection of CN emission and also detect numerous Ni I lines while Fe I remains undetected, potentially implying efficiently released gas-phase Ni. From our latest X-shooter measurements conducted on 2025-08-21 ($r_{\rm h} = 2.85$\,au) we measure production rates of $\log~Q(\mathrm{CN}) = 23.61\pm 0.05$ molecules s$^{-1}$ and $\log~Q$(Ni) $= 22.67\pm0.07$ atoms s$^{-1}$, and characterize their evolution as the comet approaches perihelion. We observe a steep heliocentric-distance scaling for the production rates $Q(\mathrm{Ni}) \propto r_h^{-8.43 \pm 0.79}$ and for $Q(\mathrm{CN}) \propto r_h^{-9.38 \pm 1.2}$, and predict a Ni-CO$_{(2)}$ correlation if the Ni I emission is driven by the carbonyl formation channel. Energetic considerations of activation barriers show that this behavior is inconsistent with direct sublimation of canonical metal/sulfide phases and instead favors low-activation-energy release from dust, e.g. photon-stimulated desorption or mild thermolysis of metalated organics or Ni-rich nanophases, possibly including Ni-carbonyl-like complexes. These hypotheses are testable with future coordinated ground-based and space-based monitoring as 3I becomes more active during its continued passage through the solar system.

We present a statistical study of spatially resolved chemical enrichment in 18 main-sequence galaxies at $z=4$--6, observed with \jwst/NIRSpec IFU as part of the ALPINE-CRISTAL-\jwst\ survey. Performing an optimized reduction and calibration procedure, including local background subtraction, light-leakage masking, stripe removal, and astrometry refinement, we achieve robust emission-line mapping on kiloparsec scales. Although line-ratio distributions vary across galaxies in our sample, we generally find mild central enhancements in [O\,\textsc{iii}]/H$\beta$, [O\,\textsc{ii}]/[O\,\textsc{iii}], [S\,\textsc{ii}]$_{6732}$/[S\,\textsc{ii}]$_{6718}$, H$\alpha$/H$\beta$, and $L_{\rm H\alpha}/L_{\rm UV}$, consistent with elevated electron density, dust obscuration, and bursty star formation accompanied by reduced metallicity and ionization parameter. These features point to inside-out growth fueled by recent inflows of pristine gas. Nevertheless, the median metallicity gradient is nearly flat over a few kpc scale, $\Delta \log({\rm O/H}) = 0.02 \pm 0.01$ dex kpc$^{-1}$, implying efficient chemical mixing through inflows, outflows, and mergers. From pixel-by-pixel stellar and emission-line characterizations, we further investigate the resolved Fundamental Metallicity Relation (rFMR). Metallicity is described by a fundamental plane with stellar mass and SFR surface densities, but with a stronger dependence on $\Sigma_{\rm SFR}$ than seen in local galaxies. Our results indicate that the regulatory processes linking star formation, gas flows, and metal enrichment were already vigorous $\sim$1 Gyr after the Big Bang, producing the nearly flat metallicity gradient and a stronger coupling between star formation and metallicity than observed in evolved systems in the local universe.

We present new ALMA [OIII]$_{88}$ observations of eight previously [CII]$_{158}$-detected galaxies at $6.8 \lesssim z \lesssim 7.7$. Six of our targets -- the primary sample -- are massive, UV-luminous galaxies drawn from the REBELS survey, while the remaining two are UV-fainter galaxies that were previously serendipitously detected through their luminous [CII] lines in the REBELS fields. We detect [OIII]$_{88}$ emission in all eight galaxies at $6.2 - 17.7\sigma$ significance, and find them to be consistent with the local dwarf galaxy relation between $L_\mathrm{[OIII]}$ and star formation rate. Our sample spans [OIII]/[CII] $\approx 1.9 - 9.6$, which is typical for the high-redshift galaxy population. Five of the primary targets benefit from JWST/NIRSpec observations, enabling a direct comparison of the [OIII]/[CII] ratio against rest-optical ISM diagnostics. We supplement our high-redshift sample with eleven $z\approx6-14$ galaxies in the literature for which similar ALMA and JWST observations are available, and furthermore compare to the [OIII]/[CII] ratios measured for local dwarf galaxies. We find that, at fixed metallicity and ionization parameter, $z>6$ galaxies show elevated [OIII]/[CII] ratios compared to local dwarfs. Instead, we find that a large [OIII]$_{4959,5007}$+H$\beta$ equivalent width -- a proxy for burstiness -- is the main driver of the high [OIII]/[CII] ratios seen in the early Universe, which is primarily due to [CII] being suppressed in bursty galaxies. Given the apparent validity of the [OIII]$_{88}$-SFR relation across most of cosmic time, as well as the abundance of young, bursty galaxies at high redshift, [OIII]$_{88}$ is set to remain a powerful ISM tracer at the cosmic dawn.

The infall of the Large Magellanic Cloud (LMC) into the Milky Way's halo impacts the distribution of stars and dark matter in our Galaxy. Mapping the observational consequences of this encounter can inform us about the properties of both galaxies, details of their interaction, and possibly distinguish between different dark matter models. N-body simulations predict a localized overdensity trailing the LMC's orbit both in baryonic and dark matter, known as the wake. We collected wide-field, deep near-infrared, and optical photometry using VIRCAM and DECam across four fields along the expected wake, covering the sky region expected to span most of its predicted density contrast. We identify over 400 stars comprising two different tracers - near main sequence turn-off stars and red giants - that map the halo between 60-100 kpc, deriving stellar halo densities as a function of sky position and Galactocentric radius. We detect (1) a break in the halo radial density profile at 70 kpc not seen in Northern halo studies, and (2) a clear halo overdensity starting also at 70 kpc, with density contrast increasing steadily toward the expected current location of the wake. If this overdensity is the LMC wake, its peak density contrast is as pronounced as the most massive LMC model considered. Contamination from unidentified substructures may bias our wake detections, so wider-area surveys with similar depth are needed for confirmation.

(abridged)Photometric redshifts for AGN (galaxies hosting an accreting supermassive black hole in their center) are notoriously challenging and currently better computed via SED fitting, assuming that deep photometry for many wavelengths is available. However, for AGN detected all-sky, the photometry is limited and provided by different projects. This makes the task of homogenising the data challenging and is a dramatic drawback for the millions of AGN that wide surveys like SRG/eROSITA will detect. This work aims to compute reliable photometric redshifts for X-ray-detected AGN using only one dataset that covers a large area: the 10th Data Release of the Imaging Legacy Survey (LS10) for DESI. LS10 provides deep grizW1-W4 forced photometry within various apertures, thus avoids issues related to the cross-calibration of surveys. We present the results from CircleZ, a machine-learning algorithm based on a Fully Connected Neural Network. CircleZ uses training sample of 14,000 X-ray-detected AGN and utilizes multi-aperture photometry. The accuracy and the fraction of outliers reached in a test sample of 2913 AGN are 0.067 and 11.6%, respectively. The results are comparable to or better than those obtained previously for the same field but with much less effort. We further tested the stability of the results by computing the photometric redshifts for the sources detected in CSC2 and Chandra-COSMOS Legacy, reaching comparable accuracy as in eFEDS when limiting the magnitude of the counterparts with respect to the depth of LS10. The method applies to fainter samples of AGN using deeper optical data from future surveys (e.g., LSST, Euclid), granted LS10-like information on the light distribution beyond a morphological type is provided. With the paper, we release an updated version of the photometric redshifts (including errors and probability distribution function) for eROSITA/eFEDS.

Early galaxy formation, initiated by the dark matter and gas assembly, evolves through frequent mergers and feedback processes into dynamically hot, chaotic structures. In contrast, dynamically cold, smooth rotating disks have been observed in massive evolved galaxies merely 1.4 billion years after the Big Bang, suggesting rapid morphological and dynamical evolution in the early Universe. Probing this evolution mechanism necessitates studies of young galaxies, yet efforts have been hindered by observational limitations in both sensitivity and spatial resolution. Here we report high-resolution observations of a strongly lensed and quintuply imaged, low-luminosity, young galaxy at $z=6.072$ (dubbed the Cosmic Grapes), 930 million years after the Big Bang. Magnified by gravitational lensing, the galaxy is resolved into at least 15 individual star-forming clumps with effective radii of $r_{\rm e}\simeq$ 10--60 parsec (pc), which dominate $\simeq$ 70\% of the galaxy's total flux. The cool gas emission unveils a smooth, underlying rotating disk characterized by a high rotational-to-random motion ratio and a gravitationally unstable state (Toomre $Q \simeq$ 0.2--0.3), with high surface gas densities comparable to local dusty starbursts with $\simeq10^{3-5}$ $M_{\odot}$/pc$^{2}$. These gas properties suggest that the numerous star-forming clumps are formed through disk instabilities with weak feedback effects. The clumpiness of the Cosmic Grapes significantly exceeds that of galaxies at later epochs and the predictions from current simulations for early galaxies. Our findings shed new light on internal galaxy substructures and their relation to the underlying dynamics and feedback mechanisms at play during their early formation phases, potentially explaining the high abundance of bright galaxies observed in the early Universe and the dark matter core-cusp problem.

We used the spectro-photometric information of ~219 million stars from Gaia's DR3 to calculate synthetic, narrow-band, metallicity-sensitive CaHK magnitudes that mimic the observations of the Pristine survey, a survey of photometric metallicities of Milky Way (MW) stars that has been mapping >6,500 deg^2 of the northern sky with CFHT since 2015. These synthetic magnitudes were used for an absolute recalibration of the deeper Pristine photometry and, combined with broadband Gaia information, synthetic and Pristine CaHK magnitudes were used to estimate photometric metallicities over the whole sky. The resulting metallicity catalogue is accurate down to [Fe/H]~-3.5 and is particularly suited for the exploration of the metal-poor MW ([Fe/H]<-1.0). We make available here the catalogue of synthetic CaHK_syn magnitudes for all stars with BP/RP information in Gaia DR3, as well as an associated catalogue of more than ~30 million photometric metallicities for high S/N FGK stars. This paper further provides the first public data release of the Pristine catalogue in the form of higher quality recalibrated Pristine CaHK magnitudes and photometric metallicities for all stars in common with the BP/RP spectro-photometric information in Gaia DR3. When available, the much deeper Pristine data greatly enhance the quality of the derived metallicities, in particular at the faint end of the catalogue (G_BP > 16). Combined, both photometric metallicity catalogues include more than two million metal-poor star candidates ([Fe/H]_phot<-1.0) as well as more than 200,000 and ~8,000 very and extremely metal-poor candidates ([Fe/H]_phot<-2.0 and <-3.0, respectively). Finally, we show that these metallicity catalogues can be used efficiently, among other applications to hunt for the most metal-poor stars, and to study how the structure of the MW varies with metallicity.

Blueshifts - or, more accurately, blue asymmetries - in broad emission lines such as CIV $\lambda$1550 are common in luminous quasars and correlate with fundamental properties such as Eddington ratio and broad absorption line (BAL) characteristics. However, the formation of these blueshifts is still not understood, and neither is their physical connection to the BAL phenomenon or accretion disc. In this work, we present Monte Carlo radiative transfer and photoionization simulations using parametrized biconical disc-wind models. We take advantage of the azimuthal symmetry of a quasar and show that we can reproduce CIV blueshifts provided that (i) the disc-midplane is optically thick out to radii beyond the line formation region, so that the receding wind bicone is obscured; and (ii) the system is viewed from relatively low (that is, more face-on) inclinations ($\lesssim40^\circ$). We show that CIV emission line blueshifts and BALs can form in the same wind structure. The velocity profile of the wind has a significant impact on the location of the line formation region and the resulting line profile, suggesting that the shape of the emission lines can be used as a probe of wind-driving physics. While we are successful at producing blueshifts/blue asymmetries in outflows, we struggle to match the detailed shape or skew of the observed emission line profiles. In addition, our models produce redshifted emission-line asymmetries for certain viewing angles. We discuss our work in the context of the CIV $\lambda$1550 emission blueshift versus equivalent-width space and explore the implications for quasar disc wind physics.

Markov Chain Monte Carlo (MCMC) excels at sampling complex posteriors but traditionally lags behind nested sampling in accurate evidence estimation, which is crucial for model comparison in astrophysical problems. We introduce reddemcee, an Adaptive Parallel Tempering Ensemble Sampler, aiming to close this gap by simultaneously presenting next-generation automated temperature-ladder adaptation techniques and robust, low-bias evidence estimators. reddemcee couples an affine-invariant stretch move with five interchangeable ladder-adaptation objectives, Uniform Swap Acceptance Rate, Swap Mean Distance, Gaussian-Area Overlap, Small Gaussian Gap, and Equalised Thermodynamic Length, implemented through a common differential update rule. Three evidence estimators are provided: Curvature-aware Thermodynamic Integration (TI+), Geometric-Bridge Stepping Stones (SS+), and a novel Hybrid algorithm that blends both approaches (H+). Performance and accuracy are benchmarked on n-dimensional Gaussian Shells, Gaussian Egg-box, Rosenbrock Functions, and exoplanet radial-velocity time-series of HD 20794. Across Shells up to 15 dimensions, reddemcee presents roughly 7 times the effective sampling speed of the best dynamic nested sampling configuration. The TI+, SS+ and H+ estimators recover estimates under 3 percent error and supply realistic uncertainties with as few as six temperatures. In the HD 20794 case study, reddemcee reproduces literature model rankings and yields tighter yet consistent planetary parameters compared with dynesty, with evidence errors that track run-to-run dispersion. By unifying fast ladder adaptation with reliable evidence estimators, reddemcee delivers strong throughput and accurate evidence estimates, often matching, and occasionally surpassing, dynamic nested sampling, while preserving the rich posterior information which makes MCMC indispensable for modern Bayesian inference.