Bayesian 3D modeling of the orbital dynamics of wide binary stars: gravitational anomaly at low acceleration from a highest-quality sample of 36 wide binaries with accurate 3D velocities

{"Python":[0,83],"modules,":[1],"data":[2,43,67,364],"files,":[3],"and":[4,192,239,248,269,289],"MCMC":[5,41,54,73,208,362,375],"output":[6,42,74,363],"figures":[7,49],"for":[8,101,138,142,156,178,235,243,274,336,346,356],"\\"Detection":[9],"of":[10,20,72,105,126,206,212,220,256,295,303,315,322,326],"Gravitational":[11],"Anomaly":[12],"at":[13],"Low":[14],"Acceleration":[15],"from":[16,182,246,373],"a":[17,102,117,179,183],"Highest-quality":[18],"Sample":[19],"36":[21],"Wide":[22],"Binaries":[23],"with":[24,59,285,360],"Accurate":[25],"3D":[26],"Velocities\\"":[27],"(arXiv:2601.21728v3,":[28],"revised":[29],"version":[30],"submitted":[31],"to":[32,38,88,149,172,196],"the":[33,60,65,80,97,112,123,127,132,135,157,175,207,226,244,251,257,304,323,337,347,357,374],"AAS":[34],"journals).":[35],"[Note:":[36,130],"Due":[37],"upload":[39],"limitations,":[40],"files":[44,75],"other":[45,252],"than":[46],"corner":[47,370],"plot":[48],"are":[50,201],"not":[51,312],"included.":[52],"However,":[53],"outputs":[55],"can":[56,76,108],"be":[57,77,109],"produced":[58],"given":[61,229,258,307],"codes":[62,195],"based":[63],"on":[64],"observational":[66],"file":[68,310,368],"provided":[69],"here.":[70],"Examples":[71],"found":[78],"in":[79,122,225,308],"previous":[81,133],"version.]":[82],"modules:":[84],"'MCMCwb3D_newton.py':":[85],"The":[86,169,194,217,228,279,301,332,342,352],"module":[87],"perform":[89],"Markov":[90],"Chain":[91],"Monte":[92],"Carlo":[93],"(MCMC)":[94],"simulations":[95,376],"using":[96],"public":[98],"package":[99],"'emcee'":[100],"fixed":[103,148],"value":[104],"G=gamma*G_N":[106],"(gamma":[107],"specified":[110],"by":[111,186],"user).":[113],"It":[114,145],"also":[115],"has":[116],"function":[118],"that":[119,200,222],"is":[120,146,166,223],"called":[121],"grid":[124,334,344,354],"search":[125],"parameter":[128,164],"space.":[129],"In":[131],"version,":[134],"prior":[136],"range":[137],"inclination":[139],"was":[140],"(-90deg,90deg)":[141],"this":[143,309],"module.":[144],"now":[147],"(-180deg,180deg).":[150],"'MCMCwb3D_general.py':":[151],"Same":[152],"as":[153,203,282],"'MCMCwb3D_newton.py'":[154,188],"but":[155,284],"generalized":[158],"gravity":[159,339,349,359],"case":[160],"where":[161],"G":[162],"(through":[163],"Gamma)":[165],"free.":[167],"'jobMCMCfit_wb3Dany_basic.py':":[168],"main":[170],"program":[171],"carry":[173],"out":[174],"Bayesian":[176],"fit":[177],"wide":[180,213],"binary":[181],"Gaia":[184,210],"sample":[185,219,281,325],"calling":[187],"or":[189,298],"'MCMCwb3D_general.py'.":[190],"'jobGridSearch_wb3D_newton.py'":[191],"'jobGridSearch_wb3D_general.py':":[193],"obtain":[197],"grid-search":[198],"solutions":[199,335,345,355],"needed":[202],"starting":[204],"values":[205],"procedure.":[209],"samples":[211],"binaries:":[214],"'gaia_purewb_3D_valueadded_basic75_corrected.csv'":[215],":":[216,278,318],"Chae+2026":[218],"binaries":[221,245],"used":[224],"study.":[227],"radial":[230,270,290],"velocities":[231,271,291],"(RVs)":[232],"were":[233,261,272],"corrected":[234,273],"gravitational":[236],"redshift":[237],"(GR)":[238],"convective":[240],"blueshift":[241],"(CB)":[242],"APOGEE":[247],"HARPS.":[249],"For":[250,263],"binaries,":[253,265],"realistic":[254],"errors":[255,302,314],"relative":[259,305],"RVs":[260,306,321],"estimated.":[262],"all":[264],"both":[266],"proper":[267,287],"motions":[268,288],"perspective":[275],"effects.":[276,300],"'gaia_purewb_3D_valueadded_basic75.csv'":[277],"same":[280],"above":[283],"measured":[286,320],"without":[292],"any":[293],"corrections":[294],"GR,":[296],"CB,":[297],"perpective":[299],"do":[311],"include":[313],"GR+CB.":[316],"'gaia_wb_RV_com2026.csv'":[317],"All":[319],"raw":[324],"306":[327],"binaries.":[328],"Auxiliary":[329],"files:":[330,365],"'gaia_purewb_3D_valueadded_basic75_corrected_gridsol_newton.csv':":[331],"pre-calculated":[333,343,353],"Newtonian":[338],"case.":[340,350],"'gaia_purewb_3D_valueadded_basic75_corrected_gridsol_general.csv':":[341],"general":[348],"'gaia_purewb_3D_valueadded_basic75_corrected_gridsol_newton140.csv':":[351],"pseudo-Newtonian":[358],"G=1.4G_N.":[361],"'fig_corner_basic75.tar.gz':":[366],"compressed":[367],"containing":[369],"plots":[371],"derived":[372]}

Authors

Institutions

Publication Details

Journal
arXiv (Cornell University)
Published
2026-09-18
DOI
https://doi.org/10.5281/zenodo.22803893
Primary Topic
Pulsars and Gravitational Waves Research
Type
preprint
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
preprint

Bayesian 3D modeling of the orbital dynamics of wide binary stars: gravitational anomaly at low acceleration from a highest-quality sample of 36 wide binaries with accurate 3D velocities

B. -C. Lee, K. -H. Chae, Y. -W. Lee, D. A. Turnshek et al.
arXiv (Cornell University)
Pulsars and Gravitational Waves Research
preprint

Bayesian 3D modeling of the orbital dynamics of wide binary stars: gravitational anomaly at low acceleration from a highest-quality sample of 36 wide binaries with accurate 3D velocities

B. -C. Lee, K. -H. Chae, Y. -W. Lee, D. A. Turnshek, X. Hernandez, V. G. Orlov, D. Lim
preprint en

Abstract

Python modules, data files, and MCMC output figures for "Detection of Gravitational Anomaly at Low Acceleration from a Highest-quality Sample of 36 Wide Binaries with Accurate 3D Velocities" (arXiv:2601.21728v3, revised version submitted to the AAS journals). [Note: Due to upload limitations, MCMC output data files other than corner plot figures are not included. However, MCMC outputs can be produced with the given codes based on the observational data file provided here. Examples of MCMC output files can be found in the previous version.] Python modules: 'MCMCwb3D_newton.py': The module to perform Markov Chain Monte Carlo (MCMC) simulations using the public package 'emcee' for a fixed value of G=gamma*G_N (gamma can be specified by the user). It also has a function that is called in the grid search of the parameter space. [Note: In the previous version, the prior range for inclination was (-90deg,90deg) for this module. It is now fixed to (-180deg,180deg). 'MCMCwb3D_general.py': Same as 'MCMCwb3D_newton.py' but for the generalized gravity case where G (through parameter Gamma) is free. 'jobMCMCfit_wb3Dany_basic.py': The main program to carry out the Bayesian fit for a wide binary from a Gaia sample by calling 'MCMCwb3D_newton.py' or 'MCMCwb3D_general.py'. 'jobGridSearch_wb3D_newton.py' and 'jobGridSearch_wb3D_general.py': The codes to obtain grid-search solutions that are needed as starting values of the MCMC procedure. Gaia samples of wide binaries: 'gaia_purewb_3D_valueadded_basic75_corrected.csv' : The Chae+2026 sample of binaries that is used in the study. The given radial velocities (RVs) were corrected for gravitational redshift (GR) and convective blueshift (CB) for the binaries from APOGEE and HARPS. For the other binaries, realistic errors of the given relative RVs were estimated. For all binaries, both proper motions and radial velocities were corrected for perspective effects. 'gaia_purewb_3D_valueadded_basic75.csv' : The same sample as above but with measured proper motions and radial velocities without any corrections of GR, CB, or perpective effects. The errors of the relative RVs given in this file do not include errors of GR+CB. 'gaia_wb_RV_com2026.csv' : All measured RVs of the raw sample of 306 binaries. Auxiliary files: 'gaia_purewb_3D_valueadded_basic75_corrected_gridsol_newton.csv': The pre-calculated grid solutions for the Newtonian gravity case. 'gaia_purewb_3D_valueadded_basic75_corrected_gridsol_general.csv': The pre-calculated grid solutions for the general gravity case. 'gaia_purewb_3D_valueadded_basic75_corrected_gridsol_newton140.csv': The pre-calculated grid solutions for the pseudo-Newtonian gravity with G=1.4G_N. MCMC output data files: 'fig_corner_basic75.tar.gz': compressed file containing corner plots derived from the MCMC simulations

arXiv (Cornell University)
Sejong University (KR)
Pulsars and Gravitational Waves Research
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.