Code to plot the Proper Motion anomaly (PMa) of any star in the Kervella+2022 Vizier catalog. The code follows the method detailed in Kervella+2019.
An example of how to use the code is given in example/PMa_example.py with the potential kwargs that can be given to the different functions.
Query can be made:
- star: Name of your target that will be serached through Simbad to get the HIP numbers
- hip: You can bypass the Simbad query by specifying HIP number, e.g. hip= 'HIP XXXX'
Useful features:
- self.data: DataFrame with queried PMa information from Gaia eDR3 and Hipparcos + any geometry information given for the system
- plot_pma(ax, epoch=['eDR3' or 'Hipparcos']): Returns ax where the specifc PMa curve is plotted
- plot_disc_extent(ax, aps, discdata, mstar): Returns ax where the disc extent is plotted as well as a Hill radius argument for the stability of the inner edge location.
- ruwe_cutoff(ax, dpc, mstar, ruwe): Returns ax where specific planet parameters can be ruled out based on RUWE (following Limbach+2024, Kiefer+2024)
Requirements are listed in requirements.txt.