Going back to basics : Accelerating exoplanet transit modelling using Taylor-series expansion of the orbital motion
(2020) In Monthly Notices of the Royal Astronomical Society 499(3). p.3356-3361- Abstract
A significant fraction of an exoplanet transit model evaluation time is spent calculating projected distances between the planet and its host star. This is a relatively fast operation for a circular orbit, but slower for an eccentric one. However, because the planet's position and its time derivatives are constant for any specific point in orbital phase, the projected distance can be calculated rapidly and accurately in the vicinity of the transit by expanding the planet's x and y positions in the sky plane into a Taylor series at mid-transit. Calculating the projected distance for an elliptical orbit using the four first time derivatives of the position vector (velocity, acceleration, jerk, and snap) is ∼100 times faster than... (More)
A significant fraction of an exoplanet transit model evaluation time is spent calculating projected distances between the planet and its host star. This is a relatively fast operation for a circular orbit, but slower for an eccentric one. However, because the planet's position and its time derivatives are constant for any specific point in orbital phase, the projected distance can be calculated rapidly and accurately in the vicinity of the transit by expanding the planet's x and y positions in the sky plane into a Taylor series at mid-transit. Calculating the projected distance for an elliptical orbit using the four first time derivatives of the position vector (velocity, acceleration, jerk, and snap) is ∼100 times faster than calculating it using the Newton's method, and also significantly faster than calculating z for a circular orbit because the approach does not use numerically expensive trigonometric functions. The speed gain in the projected distance calculation leads to 2-25 times faster transit model evaluation speed, depending on the transit model complexity and orbital eccentricity. Calculation of the four position derivatives using numerical differentiation takes ∼ 1 μs with a modern laptop and needs to be done only once for a given orbit, and the maximum error the approximation introduces to a transit light curve is below 1 ppm for the major part of the physically plausible orbital parameter space.
(Less)
- author
- Parviainen, Hannu and Korth, J. LU
- publishing date
- 2020-12-01
- type
- Contribution to journal
- publication status
- published
- keywords
- Methods: numerical, Planets, Satellites: general, Techniques: photometric
- in
- Monthly Notices of the Royal Astronomical Society
- volume
- 499
- issue
- 3
- pages
- 6 pages
- publisher
- Oxford University Press
- external identifiers
-
- scopus:85097018794
- ISSN
- 0035-8711
- DOI
- 10.1093/mnras/staa2953
- language
- English
- LU publication?
- no
- additional info
- Publisher Copyright: © 2020 The Author(s)
- id
- 552cb1f7-14f0-4974-b17b-0030b2b760d6
- date added to LUP
- 2023-02-01 10:41:57
- date last changed
- 2023-02-20 13:29:50
@article{552cb1f7-14f0-4974-b17b-0030b2b760d6,
abstract = {{<p>A significant fraction of an exoplanet transit model evaluation time is spent calculating projected distances between the planet and its host star. This is a relatively fast operation for a circular orbit, but slower for an eccentric one. However, because the planet's position and its time derivatives are constant for any specific point in orbital phase, the projected distance can be calculated rapidly and accurately in the vicinity of the transit by expanding the planet's x and y positions in the sky plane into a Taylor series at mid-transit. Calculating the projected distance for an elliptical orbit using the four first time derivatives of the position vector (velocity, acceleration, jerk, and snap) is ∼100 times faster than calculating it using the Newton's method, and also significantly faster than calculating z for a circular orbit because the approach does not use numerically expensive trigonometric functions. The speed gain in the projected distance calculation leads to 2-25 times faster transit model evaluation speed, depending on the transit model complexity and orbital eccentricity. Calculation of the four position derivatives using numerical differentiation takes ∼ 1 μs with a modern laptop and needs to be done only once for a given orbit, and the maximum error the approximation introduces to a transit light curve is below 1 ppm for the major part of the physically plausible orbital parameter space.</p>}},
author = {{Parviainen, Hannu and Korth, J.}},
issn = {{0035-8711}},
keywords = {{Methods: numerical; Planets; Satellites: general; Techniques: photometric}},
language = {{eng}},
month = {{12}},
number = {{3}},
pages = {{3356--3361}},
publisher = {{Oxford University Press}},
series = {{Monthly Notices of the Royal Astronomical Society}},
title = {{Going back to basics : Accelerating exoplanet transit modelling using Taylor-series expansion of the orbital motion}},
url = {{http://dx.doi.org/10.1093/mnras/staa2953}},
doi = {{10.1093/mnras/staa2953}},
volume = {{499}},
year = {{2020}},
}