The data for "HD 164604 c: a second giant planet on a 15-yr orbit and the constraint of the planet-planet mutual inclination"
收藏资源简介：
This file contains radial-velocity (RV) data as well as astrometric measurements from both Hipparcos and Gaia for the star HD 164604. It also provides the RV + astrometry code needed to reproduce the orbital solution of the HD 164604 planetary system. File description: RVAst_mcmc_v18.1.py: the source code to perform joint fitting of radial velocity (RV), eclipse timing variation(ETV) and Hipparcos-Gaia astrometry data. The detailed methodology are presented in Feng et.al (2023) and Xiao et.al (2024). *_posterior_*.txt: the posterior distribution of fitting parameters generated by the code. The basic fitting parameters include the orbital period, RV semi-amplitude, eccentricity, argument of periastron of stellar reflex motion, orbital inclination, longitude of ascending node, mean anomaly at the minimum epoch of RV data and five astrometric offsets of barycenter relative to GDR3. *png/pdf: resulting figures. *vels/dat: RV data from PFS and MIKE. *gost.csv: the synthetic Gaia epoch data generated from Gaia Observation Forecast Tool *hip2.abs: Hipparcos Intermediate Astrometric Data (IAD) *hg123: Astrometric measurements of HIpparcos and Gaia DR2 & DR3 Additional auxiliary files are provided in the HD164604/data directory. RVAst_mcmc_v18.1.py: ############################## Global configuration ############################target = 'HD164604'nplanet = 2RV_Enable = True # if True, then RV model will be usedTTV_Enable = False # Light travel time effect (LTTE) model,see https://zenodo.org/records/17986248GDR4_Enable = False # Only for systems with Gaia DR4 epoch astrometry, e.g. GaiaBH3, see https://zenodo.org/records/17986248HGAst_Enable = True # hiparcos-Gaia astrometry ####### do not recommmend ######corr_GDR23 = False # if True, then use GDR2-DR3 correlation of Thompson et al. 2025 (doi: 10.3847/1538-3881/ae0cbd) calibrate_rui = False # if True, you should provide Gaia Gmag and color for the target. The frame rotation between Hipparcos and Gaia will be corrected using the recommendation of Feng et al. 2024 (doi: 10.3847/1538-4365/ad27d2) ############################## PTmcmc ######################################: ndim = len(new_index) ntemps = 30 nwalkers = max(100, 2*(ndim+1)) # The number of walkers must be greater than ``2*dimension`` nsteps = 80000 # steps for each chain thin = 200 # reserved steps for 1 walker -> nsteps/thin buin = min(200, int(nsteps/thin/2)) # should < nsteps/thin, total saved sample: nwalkers*(nsteps/thin-buin) nthreads = 32 # about 2 hours for a target Please feel free to contact us (gyxiao_tdli@sjtu.edu.cn or ffeng@sjtu.edu.cn) if you have any questions about the use of this script.
