Skip to main content

Package for multiple exoplanet system modelling

Project description

ExoplanetPy

ExoplanetPy is a python package for modelling the transit light curves of systems with multiple exoplanets orbiting around their host stars. To obtain transit curves, set up the Keplerian orbital elements for each planet in the system.

The Planet() class handles each planet's orbital elements, and the System() class can be used to plot the final transit curves.

The following orbital elements are set up:

  • e: eccentricity
  • a: semi-major axis
  • omega: argument of periapsis (ω)
  • Omega: longitude of ascending node (Ω)
  • i: orbital inclination
  • r_p: planet:star radius

Each Planet() can have different initial true anomaly (ν) values, varied by the first_periastron time argument.

Installation

Installation is recommended via pip for Python 3.

pip install exoplanetpy

The package can then be imported using:

import ExoplanetPy

Usage

Access the modules using the following statements.

from ExoplanetPy import Planet
from ExoplanetPy import System

Transit Curves

Define a single Planet() and input as planet_list argument in System().
Limb darkening models are chosen in the plot() method.

p1 = Planet(e=0.0, a=8, omega=0, Omega=0, i=89.9, r_p=0.1, first_periastron=0.0)
sys = System(star_prop={'Mass': 4}, planet_list=[p1], sort=True)
sys.plot(model='Quadratic', normalise=True)

Additional Planet() objects are inputted as planet_list argument in System().

p1 = Planet(e=0.0, a=8, omega=0, Omega=0, i=89.9, r_p=0.1, first_periastron=0.0)
p2 = Planet(e=0.0, a=2, omega=0, Omega=0, i=89.9, r_p=0.05, first_periastron=0.0)
sys = System(star_prop={'Mass': 4}, planet_list=[p1,p2], sort=True)
sys.plot(model='Quadratic', normalise=True)

Direct Imaging

Take the following planetary system:

p1 = Planet(e=0.0, a=8, omega=0, Omega=0, i=89, r_p=0.1, first_periastron=0.03)
p2 = Planet(e=0.0, a=4, omega=0, Omega=0, i=89, r_p=0.07, first_periastron=0.52)
p3 = Planet(e=0.0, a=2, omega=0, Omega=0, i=87, r_p=0.05, first_periastron=0.0)
sys = System(star_prop={'Mass': 4}, planet_list=[p1, p3, p2], sort=True)

The visualize() method allows the user to obtain visual images of the actual transit.

sys.visualize(time=0.26, model='Quadratic')

Changing the fractional time parameter allows the user to obtain images at any point during transit.

sys.visualize(time=0.27, model='Quadratic')

Similarly, for time = 0.28:

sys.visualize(time=0.28, model='Quadratic')

Dependencies

ExoplanetPy has the following dependencies:

License

MIT License
© 2020 ExoplanetPy

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Files for ExoplanetPy, version 0.1.1
Filename, size File type Python version Upload date Hashes
Filename, size ExoplanetPy-0.1.1-py3-none-any.whl (7.0 kB) File type Wheel Python version py3 Upload date Hashes View
Filename, size ExoplanetPy-0.1.1.tar.gz (5.3 kB) File type Source Python version None Upload date Hashes View

Supported by

AWS AWS Cloud computing Datadog Datadog Monitoring Facebook / Instagram Facebook / Instagram PSF Sponsor Fastly Fastly CDN Google Google Object Storage and Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Salesforce Salesforce PSF Sponsor Sentry Sentry Error logging StatusPage StatusPage Status page