TOI-837 b has a true obliquity of 25.9+7.5-6.3 deg, the first planet younger than 100 Myr with accessible ψ incompatible with an aligned orbit, favoring primordial disc torque followed by disc-driven migration.
super hub Mixed citations
emcee: The MCMC Hammer
Mixed citation behavior. Most common role is method (48%).
180
Pith papers citing it
Method
48% of classified citations
abstract
We introduce a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010). The code is open source and has already been used in several published projects in the astrophysics literature. The algorithm behind emcee has several advantages over traditional MCMC sampling methods and it has excellent performance as measured by the autocorrelation time (or function calls per independent sample). One major advantage of the algorithm is that it requires hand-tuning of only 1 or 2 parameters compared to $\sim N^2$ for a traditional algorithm in an N-dimensional parameter space. In this document, we describe the algorithm and the details of our implementation and API. Exploiting the parallelism of the ensemble method, emcee permits any user to take advantage of multiple CPU cores without extra effort. The code is available online at http://dan.iel.fm/emcee under the MIT License.
hub tools
citation-role summary
method 16
background 13
dataset 1
other 1
citation-polarity summary
claims ledger
- abstract We introduce a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010). The code is open source and has already been used in several published projects in the astrophysics literature. The algorithm behind emcee has several advantages over traditional MCMC sampling methods and it has excellent performance as measured by the autocorrelation time (or function calls per independent sample). One major advantage of the algorithm is that it requires hand-tuning of only 1 or 2 parameters compared to $\sim
- method The resultingP 1D mocks are assigned the same covariance matrix as the one em- ployed in the DESI DR1 analysis. These mocks do not include contaminants or systematics; however, we analyze them with the model presented in the previous section, which accounts for both. We extract cosmological constraints using the publicly available affine-invariant Markov chain Monte Carlo (MCMC) ensemble sampleremcee 7 [118]. At each step,emceeproposes 7https://emcee.readthedocs.io/en/stable/ - 23 - 0.9 1.0 1.1
- method For the BAO sector, we utilize the dataset of [89], which incorporates fiducial cosmology corrections through the ratior d/rfid, allowing for a consistent com- parison between theoretical predictions and observations. The exploration of the parameter space is carried out using the Markov Chain Monte Carlo (MCMC) method, implemented through the publicly available Python pack- ageemcee[90]. In this context, the SNIa nuisance pa- rameterMis treated as a free parameter of the Pan- theon dataset (see
- method RPG-2019-350, and Royal Society Grant No. RGS-R2-202004. LIGO was con- structed by the California Institute of Technology and Massachusetts Institute of Technology with funding from the National Science Foundation and operates un- der cooperative agreement PHY-0757058. We acknowledge use of iPython [86], Matplotlib [87, 88], NumPy [89], SciPy [90], emcee [91] and SeaBorn [92]. This is LIGO Document P2000231 Glossary and main symbols α(⃗λ ) The fraction of physical sources that are detectable by
- method PLANCK + PP, and (ii) CC + DESBAO + PLANCK + DES, as −2 lnLtot = 𝜒2 tot. (6.8) We estimate the likelihood by implementing the model in Python and use the publicly available affine-invariant Markov Chain Monte Carlo (MCMC) ensemble sampler emcee [145] to obtain the posterior distributions of the model parameters. The resulting samples are then analyzed usingGetDist [146] to obtain marginalized 1D and 2D posterior distributions. Finally, we compare the statistical preference of the current model r
- method consistent solutions across cosmic epochs, and adopt uniform priors as listed in Tab. II. For sampling and likelihood eval- uation, we employ the nested samplerPolyChord, which is well suited for high-dimensional parameter spaces compared to affine-invariant MCMC samplers such asemcee[99-101]. The chains are analyzed usingGetDistto extract the best-fit values [102, 103]. Finally, we compare our models against flat ΛCDM using information criteria such as the Akaike Informa- tion Criterion (AIC) a
- method package that provides a common framework for writing DESI likelihoods. The BAO the- ory and likelihood is implemented in JAX [126] 20. Even though gradient-based sampling methods were implemented, we found that with analytic marginalization over broadband parameters that leaves a few sampled parameters, and using Jax just-in-time compilation and parallelization capabilities, the ensemble sampler emcee [127]21 provided well-sampled posterior estimates in a just a few minutes. In addition to MCMC
authors
co-cited works
representative citing papers
The DM distribution of unlocalized FRBs yields H0 = 73.8 +14.0/-12.3 km/s/Mpc with 18% uncertainty.
A new histogram-free likelihood method applied to simulated JWST observations of brown dwarfs shows that globular cluster ages can be determined with formal errors under 0.2 Gyr.
Correcting NIRISS charge migration reveals KELT-7 b as a metal-enriched (~92x solar) ultra-hot Jupiter with H2O, CO2 and TiO but no H- or clouds, orbiting a young metal-rich star.
Pseudo double-source plane lenses enable statistical strong lensing tomography that forecasts σ(w0) ~ 0.45 from the LSST 10-year photometric sample in flat w0waCDM cosmology.
Native-resolution retrievals on 0.97-5.27 micron JWST spectra of VHS 1256 b find a ~79% cloud deck, solar-like metallicity and C/O, and depleted 18O in a partially cloudy disequilibrium atmosphere.
LRD host galaxies show average metallicity 0.08 Z_sun with narrow stable range, challenging pristine-gas formation models while ruling out typical local AGN.
One acceleration measurement equals ~10^5 phase-space measurements for local dark matter density estimation, with acceleration outperforming Jeans modeling in both equilibrium and perturbed Milky Way simulations.
A hierarchical Bayesian framework that uses the empirical anti-correlation between AGN variability amplitude and luminosity to infer cosmological parameters from moderate-baseline light curves via importance reweighting.
GraphNPE recovers a significantly lower central density for Boötes I consistent with a core while Draco remains marginally cuspy, and demonstrates that higher-order velocity moments reduce bias in dynamical modeling.
Asteroseismic masses average 1.29 Msun for Ba dwarfs versus 1.96 Msun for Ba giants, supporting main-sequence accretion evolution from dwarfs to giants, though models fail to match the observed [hs/ls] ratio.
First narrow-line flux ratios for lensed quasar RXJ1131-1231 measured with JWST/NIRSpec IFS at ~5% precision, detecting cusp anomaly consistent with prior work.
76 N/O-enhanced galaxies at 4<z<8.5 are observed shortly after starbursts, either in the WR enrichment phase within 10 Myr or the AGB phase after 30-40 Myr following outflows.
JWST MIRI observations of 634 galaxies at 0.2<z<2 yield IR luminosity functions with faint-end slope α≈0.147, implying lower dust-obscured SFRD than previous ALMA/Herschel/Spitzer studies.
Direct detection of CO(3-2) at z=7.31 in REBELS-25 gives M_mol ~ 10^11 M_sun with f_gas ~0.95, confirming a massive molecular reservoir and showing low-J CO remains detectable in the Epoch of Reionization.
VLBI at 4.9 GHz plus multi-epoch survey data reveal a fading parsec-scale radio component near the optical center of a dwarf galaxy, interpreted as transient ejecta from IMBH accretion.
DIXE proposes CSS-based survey strategies covering 72.5% of the sky with 26-68 ks exposures and shows via simulation that MCMC collimator demodulation yields 1° point-source localization and 3° extended-source resolution.
A 6D kinematic census identifies 18 anomalous Cepheids with extreme orbits, including one possibly scattered by globular cluster E3, and finds consistency between dynamical and stellar ages.
Simultaneous measurement of low- and high-mass IMF slopes in 214 star-forming galaxies reveals diversity, weak correlation between ends, and links to stellar mass, star formation rate, and metallicity.
Four parameters suffice to describe dust attenuation curve diversity in TNG simulations, yielding a new symbolic-regression model that recovers curves and fluxes better than existing parameterizations while linking parameters to SFR surface density, metallicity, and geometry.
Derives idealized P-L-T relation from hydrostatic envelope model for LRDs and uses sparse data to obtain preliminary H0 = 120.7 with large errors as proof-of-concept.
Introduces MGIC_rv, an information criterion that combines conditional RV likelihood with an effective parameter count for selecting multi-GP models focused on radial velocities.
High-resolution M-band spectroscopy detects super-stellar SiO in TWA 5 B, implying no significant magnesium-silicate clouds and formation consistent with core accretion beyond the CO snowline or gravitational instability with solid enrichment.
Deep interferometric observations of a z≈1.12 barred spiral reveal bar-driven molecular inflows at a rate matching the galaxy's star formation rate of ~36 M⊙/yr.
citing papers explorer
-
The Hubble Missing Globular Clusters Survey IV. Ultra-faint compact satellites of the Milky Way. The case of Koposov 2
Koposov 2 is an ancient, extremely metal-poor, compact stellar system best explained as a dissolving globular cluster, not a dwarf galaxy.
-
Fabricating fiber cavity mirror substrates compatible with high coupling efficiency
A reflectometry-based pre-selection method maintains 95.3-99.2% mode matching after single-pulse CO2 ablation on 138 fiber tips pre-screened for 96.5-99.5% compatibility.
-
JADES: the mass-metallicity relation at $z=1-10$. New calibrations, extremely metal-poor galaxies, and chemical diversity
New stack-based strong-line calibrations from ~1500 spectra yield mass-metallicity relations at z=1-10 with decreasing metallicity toward higher redshift and no slope change, plus 50 EMPG candidates at 1-4% solar metallicity showing large scatter and opposite sSFR trends.
-
The colour variability of low-z SNe Ia is entirely explained by dust
Bayesian hierarchical modeling of ZTF DR2 and Foundation DR1 datasets shows dust explains all low-z SN Ia color variability after correcting for color-cut selection bias, with no residual intrinsic color term needed.
-
A Disappearing Act: Constraints From "Missing" Flares of Repeating Partial TDE Candidates
Non-detections of expected third flares in TDE 2022dbl and TDE 2020vdq support rpTDE interpretation over independent events, with modeling favoring bound main-sequence star orbits and deep initial encounters.
-
Steep Redshift Evolution of the Ionizing Escape Fraction at $z = 5$--$12$: Empirical Constraints and Comparison with Simulations
Empirical three-parameter fit to f_esc(M_h,z) yields steep redshift evolution with population-averaged escape fraction rising from ~2% at z=5 to ~9% at z=12.
-
Yoctosecond imaging of the ground state of $^{129}$Xe at the Large Hadron Collider
Bayesian global fit to Xe-Xe and Pb-Pb LHC data infers nearly maximal triaxiality for the 129Xe ground state and extracts two- and three-particle correlations.
-
21cmEMUv3: a hybrid diffusion-LSTM emulator of 21cmFAST summary observables
21cmEMUv3 emulates the cylindrical 21cm power spectrum via score-based diffusion and six other 21cmFAST observables via LSTM networks at sub-percent accuracy, then uses the emulator to infer a lower limit on soft-band X-ray luminosity from HERA data.
-
Subarcsecond Multi-line Observations of NH$_3$ with VLA toward the Class 0 Source IRAS 16293-2422
New VLA multi-line NH3 observations of IRAS 16293-2422 reveal hot gas at 200-300 K near the protostars, attributed to shock and accretion heating.
-
JWST Observations of Asteroid 2024 YR4 Rule Out a 2032 Lunar Impact and Demonstrate a New Regime for Planetary Defense Follow-up
JWST/NIRCam observations of 2024 YR4 extend the orbital arc by eight months, reduce 2032 lunar encounter uncertainty by >30x, and rule out impact with a 22,900 ± 800 km miss distance.
-
Latitudinal variations in Neptune's temperature profile observed with ALMA
ALMA CO(3-2) observations with MCMC retrievals show latitudinal temperature trends on Neptune, including a cold southern polar layer at 300-600 mbar, consistent with Voyager at 10-100 mbar but new at deeper levels.
-
Prospect of Measuring the Cosmic Dipole by Associating Strongly Lensed Gravitational Waves with Galaxy Surveys
Simulations for Einstein Telescope and Cosmic Explorer forecast that decade-scale samples of strongly lensed GWs can constrain the cosmic dipole magnitude at the few×10^{-3} level when combined with multi-image events.
-
Quenching of X-ray emission in little red dots by both Compton-thick gas and high accretion rates
X-ray non-detection of little red dots requires both Compton-thick gas columns (~1e25 cm^-2) and intrinsically weak X-ray emission with bolometric correction k_bol,X ≳ 30.
-
Milky Way Dynamics Favor Dark Matter over Modified Gravity Models
Milky Way radial and vertical dynamics are inconsistent with MOND and STVG but consistent with dark matter halos, disfavoring the former at high significance.
-
A single power law for the TRAPPIST-1 flare distribution across four orders of magnitude in energy
TRAPPIST-1 flares follow a single power law N(≥E_TESS) ∝ E_TESS^{-0.753} from 10^{29} to 10^{33} erg after sensitivity corrections and bandpass conversion.
-
Confronting Color Glass Condensate at next-to-leading order with HERA data
A Bayesian global fit at full NLO+NLL accuracy extracts the posterior distribution for the non-perturbative initial condition of the NLO Balitsky-Kovchegov equation from HERA inclusive and charm data.
-
If at First You Don't Succeed, Trispectrum: I. Estimating the Matter Power Spectrum Covariance with Higher-Order Statistics
Estimators from squeezed bispectrum and collapsed trispectrum recover unbiased small-scale matter power spectrum covariance at the percent level using 25 Quijote simulations.
-
AI-assisted modeling and Bayesian inference of unpolarized quark transverse momentum distributions from Drell-Yan data
An AI-assisted Bayesian framework extracts TMD PDFs from global Drell-Yan data using surrogate models for scalable MCMC sampling.
-
Probing Collapsed Dark Matter Halos with Fast Radio Bursts
Core-collapsed SIDM halos produce longer FRB image time delays than CDM halos, enabling future surveys to constrain self-interaction cross sections above roughly 18-40 cm²/g depending on collapse timing.
-
Posterior Predictive Checks for Gravitational-wave Populations: Limitations and Improvements
Maximum-likelihood-based posterior predictive checks detect model misspecification better than event-level versions for uncertain spin tilts, but current detector sensitivity limits their power; the Gaussian Component Spins model underpredicts high spin magnitudes and overpredicts anti-aligned tilts
-
Model-independent test of the cosmic distance duality relation with recent observational data
Two model-independent methods applied to latest SN and BAO data find the cosmic distance duality relation consistent with observations within 1 sigma and no evidence of violation.
-
Simultaneous Estimation of Seabed and Its Roughness With Longitudinal Waves
An infinite-dimensional Bayesian framework estimates seabed topography and roughness simultaneously from acoustic data by assuming statistical isotropy and using fractional differentiability.
-
Impact of projection-induced optical selection bias on the weak lensing mass calibration of galaxy clusters
Projection-induced selection bias causes 20-50% overestimation of weak lensing masses for optically selected galaxy clusters, larger on scales >3 Mpc.
-
The universal growth of magnetic energy during the nonlinear phase of subsonic and supersonic small-scale dynamos
Ensemble simulations find linear magnetic energy growth (p_nl=1) in subsonic small-scale dynamos and quadratic growth (p_nl=2) in supersonic ones, with universal conversion efficiency ~1/100 of turbulent kinetic energy flux and nonlinear phase duration ~20 t0.
-
Bring the noise: exact inference from noisy simulations in collider physics
Introduces pseudo-marginal MCMC with unbiased Poisson likelihood estimator for exact inference despite noisy collider Monte Carlo simulations.
-
Reconstructing the Thermal Sunyaev Zeldovich Power Spectrum from Planck using the ABS Method
The ABS method applied to Planck PR3 data yields a tSZ power spectrum amplitude 34% lower than the Planck 2015 best-fit when trispectrum is included.
-
Modeling and measuring the anisotropic halo 3-point correlation function: a coordinated study
Novel implementation of anisotropic 3PCF model and estimator tested on 298 halo catalogs shows degeneracy breaking between f and b1 in 3PCF-only analysis but limited added value in joint 2PCF+3PCF due to tree-level model shortcomings on small scales.
-
Gapless neutron superfluidity can explain the late time cooling of transiently accreting neutron stars
Gapless neutron superfluidity from vortex pinning explains late-time cooling of KS 1731-260 and MXB 1659-29 without requiring suppression of superfluidity.
-
A closer look at the electromagnetic signatures of Bethe-Heitler pair production process in blazars
An empirical kernel for Bethe-Heitler pair spectra is introduced and applied to argue that synchrotron emission from these pairs can explain gamma-ray emission in low- and intermediate-peaked blazars for jet regions of size ~10^15 cm and B = 5-500 G.
-
Dark Energy Survey Year 3 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing
DES Y3 3x2pt analysis constrains S8=0.776±0.017 and Ωm=0.339±0.032 in flat ΛCDM, consistent with Planck CMB results at p=0.13-0.48.
-
Gravitational-wave parameter estimation with gaps in LISA: a Bayesian data augmentation method
Bayesian data augmentation reintroduces missing LISA data segments as auxiliary variables during posterior sampling to enable consistent parameter estimation for galactic binaries despite gaps.
-
Constraints on Horndeski Gravity with Phantom Crossing
Asymptotic Cubic Galileon models can fit Planck+DESI+DES expansion history with phantom crossing and are tightly constrained by galaxy-ISW and void-force observables.
-
The 20 GeV Galactic Halo Excess: Pixel-Level Confirmation and Consistency with Sub-TeV WIMP Annihilation
Pixel-level Fermi-LAT analysis confirms a high-latitude 20 GeV Galactic halo excess consistent with ~0.55–0.72 TeV WIMP annihilation at ⟨σv⟩~10⁻²⁴ cm³ s⁻¹, with dwarf tension reduced to R≃1.6–9.3 by systematics.
-
The DESI View of the Faint Radio Source Population in LoTSS DR2
Probabilistic spectroscopic classification of 251k LoTSS radio sources yields the largest high-confidence sample and confirms LERGs accrete below ~1% Eddington while HERGs accrete above it.
-
An equal mass ratio supermassive binary black holes in Q J0158-4325 with periodic microlensing signature?
A triple-disk SMBBH model with mass ratio q~1 reproduces the observed 173-day microlensing periodicity and SED in Q J0158-4325 via Bayesian fitting.
-
The Impact of Non-Gaussian Line Spread Functions on Stellar Kinematic Recovery: Consequences for Dynamical Models
Non-Gaussian LSF shapes bias kinematic extraction from spectra; matching the LSF of templates to the target reduces dispersion bias below 1%.
-
Late-time evolution of the interacting stripped-envelope supernova 2017dio
Late-time data on SN 2017dio yield mass-loss rates of ~0.2 M_sun/yr peak and ~0.06 typical, with H-rich CSM from a companion and dust masses 0.001-0.02 M_sun, indicating sudden mass-loss increase.
-
TOA_SP: A Multi-Strategy Framework for Single-Pulse Timing
toa_sp applies multiple single-pulse timing strategies to 688 pulses from RRAT J1913+1330, yielding 1.33 ms weighted RMS residuals (24% better than PSRCHIVE) while keeping all pulses.
-
A hidden reionization prior biases cosmological inference
No monotonic reionization history fits Planck τ, patchy kSZ from SPT/ACT, and Lyα endpoint; an early ionization phase at z≳12 relaxes ∑mν<0.39 eV and shifts σ8 via As-τ degeneracy.
-
MSA-3D: Rotation Curves and Dark Matter Fractions at z~0.5-1.7 with JWST/NIRSpec
New JWST data on 23 galaxies at 0.5<z<1.7 show median dark matter fraction of 0.63 at effective radius with 0.2 dex scatter, and a mix of rising, flat, and falling rotation curves.
-
Dust destruction signals shock-accelerated outflows in the nearby active galaxy NGC 1068
Line ratio diagnostics in NGC 1068 indicate AGN outflows are shock-accelerated, with outflowing gas dust-free and 19-110 times denser than disk gas.
-
Revisiting the 'Lensing is Low' Problem with UNIONS
New UNIONS galaxy-galaxy lensing data around CMASS galaxies indicates no significant lensing is low problem, with joint HOD fits to GGL and GC favoring a slightly lower matter power spectrum amplitude than Planck.
-
Observing a 542-day transiting giant with large TTVs: The 2025 transit of HIP 41378 f and new constraints on the outer system
New 2025 transit timing for HIP 41378 f confirms large TTVs and is combined with prior data on planets d and e in an N-body model to update ephemerides and predict future transits.
-
When the black holes align: a subpopulation of aligned massive binary black holes observed via gravitational waves
A non-parametric reconstruction of the GWTC-5 black-hole population finds a separate high-mass subpopulation (M1 ≳ 50 M☉) with positively aligned spins (χ_eff ≈ 0.4).
-
Weak-lensing mass calibration of \emph{Planck} Sunyaev--Zel'dovich clusters with HSC-SSP Year~3
Weak-lensing calibration of 19 Planck SZ clusters with HSC Year 3 yields 1-b = 0.73^{+0.10}_{-0.11} at z_eff ~0.24 after forward-modeling selection, Eddington bias, and miscentering.
-
MINDS survey of silicates in T Tauri disks: Correlation between dust and gas
Spectral decomposition of 26 T Tauri disks finds Mg-rich silicate dust with 5-24% crystallinity and correlations linking annealed silica to stronger CO2 emission and forsterite to stronger H2O emission.
-
JWST Absorption-Line Analysis of UV-Bright Galaxies at $z=7.2-10.6$: Early Chemical Enrichment Traced by C, O, Mg, Al, Si, and Fe
JWST spectra reveal that two z~7 galaxies already show near-solar iron-to-silicon ratios with no strong odd-even effect, favoring early Type Ia supernovae over pair-instability supernovae as the source of iron enrichment.
-
Probing a new subclass of llGRB-SN transients: Insights from EP250304a and its associated supernova
EP250304a/SN 2025fhm is presented as a member of an emerging subclass of shocked cocoon-dominated low-luminosity GRB-SNe based on spectral, photometric, and light-curve modeling comparisons to prior events.
-
Early Multiwavelength Observations of AT 2026fgk: The Luminous Afterglow to Sub-luminous GRB 260310A, Identified Independently of a Gamma-ray Trigger
First blind optical identification of a z=0.153 sub-luminous GRB afterglow with Ic-BL SN, yielding a volumetric rate consistent with on-axis high-luminosity long GRBs.
-
Posterior sampling in the Age of Emulators
Compares MH, MALA, HMC, NUTS, and AIES on differentiable likelihood emulators for ΛCDM and sterile-neutrino models, finding MALA and MH competitive in wall time despite NUTS needing fewer samples.