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 [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!
time, phase, density = density_through_orbit_theoretical(resolution=0.01, show_plot=True, load_directly=False)
time, phase, density = density_through_orbit_theoretical(resolution=0.01, show_plot=True, load_directly=False)
In [3]:
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)
