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%).
137
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.
Derives analytical solutions and fitting formulae for Lyα spectra under cylindrical geometry including recoil and velocity gradients, validated against Monte Carlo simulations.
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.
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.
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.
Simulations forecast that 10 years of Einstein Telescope and Cosmic Explorer data could detect the cosmic dipole magnitude using strongly lensed GW events, with tighter bounds from combining double, triple, and quadruple lensed systems.
High-resolution interferometric imaging of eight post-AGB circumbinary discs reveals diverse inner-rim substructures including azimuthal brightness enhancements and arc-like features not explained by inclination alone.
SPHEREx data confirm 77 new luminous heavily reddened quasars at 1.5<z<3.9 that are hot-dust poor relative to unobscured quasars, supporting a blow-out feedback phase.
An Er3+ spin sensor in CaWO4 measures the hyperfine NMR spectrum of a single 93Nb nuclear spin-9/2 with Hertz resolution, determining its site, position, quadrupolar tensor, and two new terms in the spin Hamiltonian.
Galaxy pairwise peculiar velocities from Cosmicflows-4 yield M_ν = 0.24^{+0.34}_{-0.18} eV and η² = 2.14^{+0.30}_{-0.32} (7σ non-zero asymmetry) in the CMB framework, consistent with prior Planck results.
Non-polynomial quasi-topological gravity models reproduce the standard thermal history, generate dynamical dark energy of geometric origin, and fit supernova, cosmic chronometer, and BAO data competitively with ΛCDM.
DESI DR2 and ACT DR6 data yield 17σ LRG-velocity, 8.3σ ELG-velocity, and 6.8σ QSO-velocity detections plus a 3.1σ velocity-velocity signal, producing f_NL^loc = 15.9_{-34.4}^{+34.6} from the velocity field.
Bayesian MCMC fitting of a phenomenological model reveals that isospin asymmetry strongly suppresses alpha preformation probability in superheavy nuclei while reproducing shell effects and experimental half-lives.
citing papers explorer
-
Detection of a parsec-scale, compact, and fading ejecta from an accreting massive black hole
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.
-
Old and Bright: The Remarkable Radio Brightening of the Engine-driven SN 2012au Several Years After Explosion Signals the Birth of a PWN
Late-time radio observations of SN 2012au show re-brightening best explained by emission from a newborn pulsar wind nebula rather than continued shock interaction with circumstellar material.
-
The properties of tidal disruption event infrared counterparts produced by dust rings and inference of the observing angle
A toy model of dust rings in TDEs predicts brighter IR emission on-axis, explaining X-ray/IR correlations and enabling viewing-angle constraints from observed light curves.
-
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.
-
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.
-
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.
-
The January 2010 flare of Mrk421: Insights from a stochastic acceleration model
The January 2010 flare of Mrk421 shows spectral variability consistent with stochastic acceleration, including development of a Maxwellian pile-up component in the particle distribution on three nights.
-
Revisiting the Reported Period of FRB 20201124A Using MCMC Methods
An MCMC-based period search method recovers previously reported candidate periods in FRB 20201124A observations.
-
Updating the PATH framework with FRB host galaxy models
PATH is extended with three fitted P(m_r|z) prior models combined with P(z|DM), raising host-association confidence for ASKAP FRBs while showing fainter-than-expected host magnitude distribution.
-
Direct simulations of very high energy cosmic ray acceleration in 3D MHD model of a compact star cluster
3D MHD simulations of young massive star clusters find proton acceleration to hundreds of TeV near O-star termination shocks, with even faster acceleration to over 100 TeV in under 100 years when a supernova remnant expands inside the core.
-
Constraining Orbital Eccentricity of a Supermassive Black Hole Binary Candidate PKS 2131-0211
Bayesian fitting of an eccentric Keplerian orbit to the radio light curve of PKS 2131-021 gives e = 0.053 ± 0.015 without red noise but favors a circular orbit plus DRW noise with e < 0.15.
-
Earliest simultaneous multi-color optical observations of GRB 230328B: from 41 seconds to the host-galaxy identification
New early multi-wavelength data on GRB 230328B shows afterglow with early bump and late achromatic rebrightening at ~4000 s, modeled via MCMC as forward shock plus late energy injection in a dusty S0 host with AV~0.8 and no supernova signature.
-
Searching for links between energetic millisecond pulsars and repeating fast radio bursts
Wideband observations show M28A giant pulses differ from FRB 20200120E bursts in duration, luminosity, timing statistics, and spectral structure, yielding no strong evidence for a direct link.
-
Optical Super-orbital Modulation of SMC X-1: Disk Precession and a Revised Pulsar Mass
A modified ellipsoidal modulation model with precessing disk irradiation effects revises the SMC X-1 pulsar mass to approximately 1.35 solar masses.
-
Iron line diagnostics of the stellar wind in X1908+075
Chandra HETGS observations of X1908+075 detect variable Fe Kα fluorescence, model NH orbital changes to derive inclination 46° and wind mass-loss rate 9.1e-7 solar masses per year, confirming a classical wind-fed supergiant X-ray binary.
-
A 14-year-old Mystery: The Peculiar Case of the Engine-driven SN 2012ap
Late-time radio rebrightening in SN 2012ap is consistent with either progenitor mass-loss variation producing a density enhancement or an off-axis energetic jet viewed at large angle, potentially reclassifying it as GRB-like rather than weakly engine-driven.
-
A series of unfortunate events: CHIME/FRB misclassification of a Galactic pulsar as a periodic fast radio burst
A reported periodic fast radio burst is reclassified as Galactic pulsar emission due to CHIME calibration and beam-pointing error.