Doppler Theoretical in four scenarios

Doppler Evolution in Various Orbital Scenarios¶

In this notebook, we will provide theoretical examples of Doppler evolution through different orbital scenarios:

Conic orbit
Spiral orbit
Disc in a main orbit (or ballistic motion, referred to as disc for simplicity)
Spiral in a main orbit

The input in these scenarios will be a time array, which we will generate artificially in the next cell.

In [1]:

Copied!

import xraybinaryorbit
from xraybinaryorbit import *
import xraybinaryorbit
from xraybinaryorbit import *


HELLO, nice to see you! :)

PLEASE READ THIS, IT'S VERY IMPORTANT:

These are the units that must be used within this package:

- Rstar: Solar radius
- Mstar: Solar masses
- Inclination: Sexagesimal degrees
- Periapsis: Sexagesimal degrees
- Semimajor: Stellar radius
- Periods: Days (Periods in the case of the period_sliding_window function will support any units)
- Iphase: Radians

A list of the functions contained in this package will be displayed by runing the function list_functions().

As these functions use a lot of parameters, which can sometimes be difficult to handle, we have implemented a user-friendly method for parameter input:
A form will be displayed, and the parameters will be saved in the directory for further interactions. These saved parameters will be used if new parameters are not provided.
For the function to work, the submit button must be pressed.

Please, take into account that fits in general will take A LOT of time to complete.

If you need help, contact graciela.sanjurjo@ua.es.

In [2]:

Copied!

#In first place we define a time array
time_array=np.arange(0,2*24*60*60,1)
#In first place we define a time array
time_array=np.arange(0,2*24*60*60,1)

In [3]:

Copied!

time_orbit, phase_orbit, ev_orbit = doppler_orbit_theoretical(time_array, units="keV", show_plot=True, load_directly=False)
time_orbit, phase_orbit, ev_orbit = doppler_orbit_theoretical(time_array, units="keV", show_plot=True, load_directly=False)

No description has been provided for this image

In [4]:

Copied!

time_spiral, phase_spiral, ev_spiral = doppler_spiral_theoretical(time_array,  units="keV", show_plot=True, load_directly=True)
time_spiral, phase_spiral, ev_spiral = doppler_spiral_theoretical(time_array,  units="keV", show_plot=True, load_directly=True)

In [5]:

Copied!

time_disc_in_orbit, phase1_disc_in_orbit, phase2_disc_in_orbit, ev_disc_in_orbit = doppler_disc_theoretical(time_array, units="keV",show_plot=True, load_directly=True)
time_disc_in_orbit, phase1_disc_in_orbit, phase2_disc_in_orbit, ev_disc_in_orbit = doppler_disc_theoretical(time_array, units="keV",show_plot=True, load_directly=True)

In [6]:

Copied!

time_spiral_in_orbit, phase1_spiral_in_orbit, phase2_spiral_in_orbit, ev_spiral_in_orbit = doppler_spiral_in_orbit_theoretical(time_array,units="keV",show_plot=True, load_directly=True)
time_spiral_in_orbit, phase1_spiral_in_orbit, phase2_spiral_in_orbit, ev_spiral_in_orbit = doppler_spiral_in_orbit_theoretical(time_array,units="keV",show_plot=True, load_directly=True)