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%).
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
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
-
J0404+1112: A 3-Hour Eclipsing White Dwarf-Brown Dwarf Probing Multiple Atmospheric Regimes
J0404+1112 is a 2.93 hr period totally eclipsing WD+BD system with a hot DA white dwarf (T_eff ~28,000 K) and ~40 M_Jup brown dwarf, enabling isolation of nightside emission and serving as a JWST atmospheric benchmark.
-
Strong Constraints on Millisecond Pulsar Injection Spectra from Fermi-LAT Observations of the Galactic Center
Joint prompt and inverse-Compton modeling of MSP-injected e± yields strong upper limits on η_e/η_γ from recent Fermi-LAT GCE spectra that exceed MAGIC globular-cluster bounds.
-
Cosmological constraints from neighbor-density-weighted marked correlation functions
Neighbor-density-weighted marked correlation functions improve FoM for Ω_m–σ_8 by 1.7–2.5× over standard 2PCF using emulators from 129 w0waCDM+∑m_ν simulations.
-
Constraints on Schwarzschild Black Hole in a Generalized Dehnen-Type $(1,4,\gamma)$ Dark Matter Halo via the S2 Star Orbit around Sgr A$^\star$
Derives perihelion shift equations for S2 star in generalized Schwarzschild-Dehnen BH-DM spacetime and constrains gamma, rho_s, rs via MCMC on Do et al. and Gillessen et al. datasets.
-
JWST Advanced Deep Extragalactic Survey (JADES) Data Release 5: stellar population catalogue for galaxies in GOODS-N and GOODS-S
JADES DR5 delivers a public catalog of Bayesian-inferred stellar masses, SFRs, SFHs, dust, metallicities, and AGN contributions for ~500k galaxies via Prospector with an evolving SFMS prior.
-
Euclid preparation: Testing multi-field inflation with galaxy power spectrum and bispectrum
Validates redshift-space power spectrum and bispectrum analysis on Abacus-PNG mocks to recover unbiased f_NL constraints for Euclid spectroscopic sample.
-
Detection of a dark matter subhalo in the strongly lensed system PJ011646
A subhalo of M200 = 2.78e10 solar masses and concentration 30 is detected at 5.8 sigma in PJ011646 via ALMA imaging and grid-based NFW search after fitting an elliptical power-law plus multipole macromodel.
-
Connecting CGM enrichment with Lyman alpha emitters at 2.9 < z < 6.7
Blind MUSE search detects 156 LAEs at 2.9<z<6.7 with 34 C IV and 14 Mg II absorber associations, indicating low-mass galaxies enrich the CGM.
-
Mare versus highland lunar impact flash light curve dichotomy
Highland lunar impact flashes decay more slowly and last longer than mare flashes, indicating that early cratering stages depend on lunar surface lithology.
-
Bulgeless Evolution And the Rise of Discs (BEARD) III. A numerical simulation view of satellites around Milky-Way analogues
Simulation comparison finds bulgeless galaxies host more centrally concentrated, disc-aligned satellites with steeper faint-end luminosity functions than bulge-dominated controls, reflecting co-evolution and quieter merger histories.
-
Modeling the probability distribution for cosmological analysis with photometrically classified samples
A simplified likelihood treating photometric contamination as a redshift-dependent Gaussian mean shift is strongly favored by Bayes factor over the BEAMS two-component model and improves cosmological constraints on the DES-Dovekie sample.
-
Perturbation Dynamics and Structure Formation in Extended Proca-Nuevo Gravity
Extended Proca-Nuevo gravity modifies the background expansion via a vector field algebraic constraint but leaves the matter growth equation identical to general relativity.
-
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.
-
Dynamical dark energy from Kretschmann scalar at low redshifts
Replacing the cosmological constant with the Kretschmann scalar yields a dynamical dark energy model that fits supernova and cosmic chronometer data and produces a phantom-crossing equation-of-state parameter w(z) similar to recent phenomenological fits.
-
Directional Tests of the Cosmic Distance Duality Relation using Pantheon+ and BAO
No evidence for directional anisotropy in the cosmic distance duality relation is found, yielding a robust 95% upper bound of 0.025 on the dipole amplitude after accounting for survey selection effects.
-
Probing soft signals of gravitational-wave memory with space-based interferometers
Space-based detectors can measure soft displacement-memory signals from gravitational waves at SNR greater than or equal to 10.
-
Testing the Spacetime Geometry of Sgr A* with the Relativistic Orbit of S2 star
Multiple black-hole and naked-singularity spacetimes are tested against S2 star astrometric and spectroscopic data, showing statistical degeneracy among several models under current observations.
-
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.
-
Probing the sensitivity of dark energy dynamics to equation of state parametrization flexibility
Flexible dark energy equation of state parametrizations mildly favor dynamical phantom-like behavior at z~1-2, but the preference is only ~2 sigma and sensitive to the model form and extrapolation.
-
GW250114: testing Hawking's area law and the Kerr nature of black holes
GW250114 data confirm the remnant black hole ringdown frequencies lie within 30% of Kerr predictions and that the final horizon area is larger than the sum of the progenitors' areas to high credibility.
-
Combined tracer analysis for DESI 2024 BAO
Combining LRG and ELG tracers with bias weighting improves BAO constraints by 11% on alpha_iso and 7% on alpha_AP in DESI DR1 data for the 0.8<z<1.1 bin.
-
Constraints on primordial non-Gaussianity from Quaia
Reports f_NL = -20.5^{+19.0}_{-18.1} (68% CL) from combined Quaia quasar auto-correlation and CMB lensing cross-correlation assuming p_phi=1, or -28.7^{+26.1}_{-24.6} for p_phi=1.6.
-
RedDots: Magnetic field of the nearby active M dwarf GJ 729, and a search for companions
GJ 729 exhibits a weak, evolving large-scale magnetic field (50-145 G) and a persistent ~7 d radial velocity signal that could be a ~1.5-2 Earth-mass planet or residual stellar activity.
-
Unveiling the Milky Way with a Gaia DR3 census of OB-type stars within 2 kpc. I. Tracing local Galactic structure, massive star-forming regions and core-collapse supernova progenitors
A Gaia DR3-based census of 105,971 OB stars within 2 kpc maps local Galactic structure and identifies over 4,200 core-collapse supernova or black hole progenitor candidates.
-
Colour evolution in the radio afterglow of GRB 241025A
Multi-band data for GRB 241025A require an ad-hoc factor-500 increase in shocked-material optical depth to match the observed radio spectral evolution within a structured-jet forward-shock model.
-
CHEOPS observations of V1298 Tau: updated planetary densities and implications on the early evolution of the young system
CHEOPS photometry improves radii for V1298 Tau b, c, and d, producing revised densities of 0.06-0.23 g cm^{-3} that indicate differential atmospheric evolution and no requirement for past dynamical excitation.
-
Direct VLBI evidence for a buried AGN in the triple-merger LIRG UGC 2369S
VLBI imaging detects a compact radio source with Tb > 10^7 K and flat spectrum in the northern core of UGC 2369S, confirming a buried low-luminosity AGN.
-
EasyScan_HEP 2: Agent-Ready Parameter Scans for High-Energy Physics
EasyScan_HEP 2 adds AI-agent interfaces to a HEP parameter scan framework for natural-language to .ini config translation and new sampler integration.
-
The Companions to B and A Stars Snapshot (C-BASS) Survey: I. Discovery of a Young Brown Dwarf Companion to HIP 17453
Discovery of an L2 brown dwarf companion (mass 53 M_Jup, age ~280 Myr) to HIP 17453 A via proper-motion confirmed imaging and spectroscopy.
-
Probing the Nature of Lyman Continuum Emitting and Low-metallicity Galaxies Using the SKA
SKA-Mid is predicted to yield samples of 10-100 low-metallicity Lyman continuum emitting galaxies per square degree, enabling multi-wavelength studies of feedback processes linked to ionizing photon escape.
-
Fast Simultaneous Surveys with On-the-Fly Mapping
MeerKLASS applies on-the-fly imaging on MeerKAT for a 10,000 sq. deg. UHF-band continuum survey at 14 arcsec resolution and 25 μJy/beam rms, run commensally with HI intensity mapping.
-
Weak Lensing with SKAO: Cosmic Shear Cosmology
Forecasts indicate SKAO cosmic shear can constrain S8 to 5% alone and 3% combined with LSST or Euclid, with added benefits from polarization and kinematics for separating lensing signals.
-
A cross-epoch endpoint-consistency test of a single effective scaling from dark energy to inflation
A single power-law effective scaling anchored at today's Hubble rate is extrapolated to inflation and required to match the Starobinsky plateau within O(1) factors, selecting γ ≈ 0.49 and β ≈ 0.68.
-
Comparative Study of Two Luminous Red Novae I. Progenitor Modeling and Dust Formation
Binary evolution modeling constrains donor masses of 14-23 solar masses for two luminous red novae and shows dust masses are 1-5 orders of magnitude below total ejected envelope masses.
-
Constraints on Cosmic Strings from the Curl-Mode CMB Lensing Power Spectrum measured by ACT DR6
ACT DR6 curl lensing yields 2σ bounds Gμ P^{-1} ≤ 3.5×10^{-5} (small-P) and Gμ ≤ 5.0×10^{-5} (P=1), tightened further with Planck 2013 data.
-
Revisiting the Exo-Mercury Candidate GJ 367 b with ESPRESSO and a Self-Consistent Tidal Distortion Model
Revised mass of 0.503 M_Earth and radius of 0.736 R_Earth for GJ 367 b give a density of 6.9 g cm^{-3} and an iron fraction of 50-70% via new tidal and composition modeling.
-
Detectability of deuterium in spectra of early-type stars
Synthetic spectra and MCMC analysis establish deuterium detection limits in early-type stars as a function of SNR, Teff, and v sin i, with upper limits D/H < -5.5 dex for HD 32115 and < -4.9 dex for 21 Peg.
-
A HINSA view of cosmic-ray ionization in IC 348 and NGC 1333: evidence for a strong low-energy cosmic-ray disparity
HINSA observations show CRIR decreasing with H2 column density in both clouds, with IC 348 values an order of magnitude above NGC 1333, modeled as an order-of-magnitude difference in low-energy CR populations from local acceleration sources.
-
Morphology, sizes, and scatter in a large sample of distant quiescent galaxies
New high-redshift quiescent galaxy sample shows size decreasing with redshift and wavelength, with stellar mass plus redshift sufficient to predict size but large residual scatter.
-
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.
-
The mass of TOI-1883 b: A low density super-Neptune in the ridge regime transiting an early-M dwarf
Mass of 13.7 Earth masses and density 0.4 g cm^{-3} measured for TOI-1883 b, a super-Neptune in the ridge regime around an early-M dwarf, with implications for disk migration and photoevaporation.
-
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.
-
The VMC survey -- LV. The coherent expansion of the SMC
Improved VMC proper motions show SMC expansion along SE-NW axes consistent with LMC tides, radial inward motions with no rotation, and northward motion in older RGB stars from a prior interaction.
-
The Pan-Pacific Planet Search -- IX. A menagerie of companions orbiting evolved stars
Resolves six speculative companions into one giant planet, one eccentric brown dwarf, two low-mass stars, and two stars with no detectable companions.
-
Bayesian estimation of spectral parameters of the 6.7-GHz methanol maser G339.884-1.259 from GRAO observations
Bayesian MCMC decomposition of the G339.884-1.259 methanol maser spectrum identifies seven velocity components and statistically prefers the Voigt profile over Gaussian or Lorentzian based on AIC, BIC, RMSE, and R².
-
It's Not Just Star Formation: A trend of low dark matter densities in the Andromeda dwarf galaxy system
Five of seven modeled M31 dwarf spheroidals show anomalously low central DM densities at 150 pc, with star formation heating disfavored as the sole cause.
-
The energy conditions and model selection in the local Universe
Local HII galaxy and cosmic chronometer data favor the Rh=ct universe over ΛCDM at 92% vs 8% likelihood, with Rh=ct satisfying all energy conditions while ΛCDM violates the strong energy condition at z≲2.
-
Testing cosmic anisotropy with the Combo correlation of gamma-ray bursts
Larger sample of 244 GRBs with Combo correlation shifts best-fit anisotropy longitude by 54° from Pantheon-only result and deviates >1σ in hemisphere method, unlike smaller A118 sample, indicating potential to reduce fake signals from inhomogeneous distributions.