Density Through the Orbit and Absorption Column¶

In this Jupyter notebook, we will provide two examples of calculations related to wind density. All calculations assume that the wind is unionized and spherically symmetric, following the CAK model (Castor, Abbott, and Klein) as described in:

Castor, J. I., Abbott, D. C., & Klein, R. I. (1975). Radiation-driven winds in Of stars. Astrophysical Journal, 195, 157-174.

Objectives¶

Calculate the wind density throughout the orbit.
Calculate the variations in the absorption column as a function of orbital phase and time.

Methodology¶

Wind Density Calculation:
- Assume a spherically symmetric, unionized wind distribution.
- Use the CAK model to describe the wind properties.
Orbital Considerations:
- Take into account that the areolar velocity in an eccentric system varies throughout the orbit.
- Provide wind density and absorption variations based on orbital phase and time.

Input Parameters¶

A form will be displayed to input the necessary parameters for the orbital calculations.
An initial set of parameters will be provided as an example to guide the user.

In [5]:

Copied!

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

In [6]:

Copied!

time, phase, density = density_through_orbit_theoretical(resolution=0.01, show_plot=True, load_directly=True)
time, phase, density = density_through_orbit_theoretical(resolution=0.01, show_plot=True, load_directly=True)

semimajor: 1.5
orbitalperiod: 11.5
eccentricity: 0.2
periapsis: 180.0
Rstar: 37.0
wind_infinite_velocity: 500.0
Mass_loss_rate: 2.5e-05
beta: 0.8

No description has been provided for this image

In [7]:

Copied!

time, phase, NH = absorption_column_through_orbit_theoretical(resolution=0.01,show_plot=True, load_directly=False)
time, phase, NH = absorption_column_through_orbit_theoretical(resolution=0.01,show_plot=True, load_directly=False)

semimajor: 1.51
orbitalperiod: 11.5
eccentricity: 0.1
periapsis: 180.0
inclination: 30.0
Rstar: 37.0
wind_infinite_velocity: 500.0
Mass_loss_rate: 2e-06
beta: 0.5

In [ ]: