barycorrpy (v0.4.4)
Please join the google group for updates regarding bug reports, new versions etc:
To sign up for updates, please join the Google Group linked here -
https://groups.google.com/g/barycorrpy
Barycorrpy is the Python version of Jason Eastman and Jason Wright's IDL code BaryCorr based on Wright and Eastman (2014). BCPy can be used to calculate the barycentric velocity correction for a star with an accuracy well below 1 cm/s .
To do this, it takes into consideration the following effects-
- Revolution of the Earth to consider position and velocity of the geocenter with respect to the Solar System barycenter
- Rotation of the Earth
- Precession, nutation and polar motion of the Earth, along with the above to calculate the position and velocity of the observatory with respect to the geocenter
- Gravitational time dilation due to objects of the Solar System
- Leap second offset
- Proper motion and systemic radial velocity of the star
- Parallax
- Shapiro delay
The installation instructions and the guide on how to run and use the code are explained in the [wiki](https://github.com/shbhuk/barycorrpy/wiki
There is also a video tutorial (courtesy of the Exoplanet Modeling and Analysis Center ) available here, describing how to get started with barycorrpy.
It is important to note that the code does not just output a barycentric velocity which must be subtracted from the raw Radial Velocity. It outputs the net radial velocity after correcting for barycentric correction. This is because the correction involves a cross term due to the relativistic addition. Therefore include the zmeas in the input parameters.
JDUTC to BJDTDB converter
As explained in Eastman et al. 2010, we also include a JDUTC to BJDTDB time converter.
We include the following corrections -
- Clock Correction - To correct for difference between UTC and TDB time scales.
- Geometric Correction - Light travel time from observatory to Solar System Barycenter.
- Einstein Correction - Relativistic correction since Earth is not an inertial frame.
The output of our function - utc_tdb.JDUTC_to_BJDTDB() matches the web applet converter to about 0.2ms . Therefore for applications requiring higher precision it should not be used.
We advocate not using the standard Astropy converter this due to the same leap second staleness issue. Therefore for cases requiring such high precision, the leap second should be checked for and be updated as is done by this code.
INSTALLATION INSTRUCTIONS
The instructions for installation and getting started for this package are detailed in the wiki.
Citation
To cite using this code you can cite this paper - RNAAS.
Guidelines to cite this package can be found here.
To sign up for updates, please join the Google Group linked here -
https://groups.google.com/forum/#!forum/barycorrpy