PHOEBE 2.0 Documentation

2.0 Docs

  • 1.0
  • 2.0a
  • ver: 2.0


IPython Notebook | Python Script

Limb Darkening

Setup

As always, let’s do imports and initialize a logger and a new bundle. See Building a System for more details.

%matplotlib inline
import phoebe
from phoebe import u # units
import numpy as np
import matplotlib.pyplot as plt

logger = phoebe.logger()

b = phoebe.default_binary()
WARNING: Constant u'Gravitational constant' is already has a definition in the u'si' system [astropy.constants.constant]
WARNING: Constant u'Solar mass' is already has a definition in the u'si' system [astropy.constants.constant]
WARNING: Constant u'Solar radius' is already has a definition in the u'si' system [astropy.constants.constant]
WARNING: Constant u'Solar luminosity' is already has a definition in the u'si' system [astropy.constants.constant]
/usr/local/lib/python2.7/dist-packages/astropy/units/quantity.py:782: FutureWarning: comparison to None will result in an elementwise object comparison in the future.
  return super(Quantity, self).__eq__(other)

We’ll just add ‘lc’ and ‘mesh’ datasets

b.add_dataset('lc', times=np.linspace(0,1,101), dataset='lc01')
<ParameterSet: 15 parameters | contexts: compute, dataset>
b.add_dataset('mesh', times=[0.75], dataset='mesh01')
<ParameterSet: 2 parameters | contexts: compute, dataset>

Relevant Parameters

print b['ld_func_bol@primary']
Parameter: ld_func_bol@primary@component
                       Qualifier: ld_func_bol
                     Description: Bolometric limb darkening model
                           Value: logarithmic
                         Choices: linear, logarithmic, quadratic, square_root, power
print b['ld_func_bol@primary'].choices
['linear', 'logarithmic', 'quadratic', 'square_root', 'power']
print b['ld_coeffs_bol@primary']
Parameter: ld_coeffs_bol@primary@component
                       Qualifier: ld_coeffs_bol
                     Description: Bolometric limb darkening coefficients
                           Value: [ 0.5  0.5]
                  Constrained by:
                      Constrains: None
                      Related to: None
print b['ld_func@lc01']
ParameterSet: 2 parameters
    ld_func@primary@lc01@dataset: interp
  ld_func@secondary@lc01@dataset: interp
print b['ld_func@lc01@primary'].choices
['interp', 'linear', 'logarithmic', 'quadratic', 'square_root', 'power']

Note that ld_coeffs isn’t visible (relevant) if ld_func==’interp’

b['ld_func@lc01@primary'] = 'logarithmic'
print b['ld_coeffs@lc01@primary']
Parameter: ld_coeffs@primary@lc01@dataset
                       Qualifier: ld_coeffs
                     Description: Limb darkening coefficients
                           Value: [ 0.5  0.5]
                  Constrained by:
                      Constrains: None
                      Related to: None
                 Only visible if: ld_func:!interp

Influence on Light Curves (fluxes)

b.run_compute(model='mymodel')
<ParameterSet: 76 parameters | kinds: mesh, lc>
axs, artists = b['lc01@mymodel'].plot()
../../_images/limb_darkening_19_0.png

Influence on Meshes (Intensities)

COMING SOON

TODO: vary over different ld_func or ld_coeffs and show how they affect values and plots

Prev: Distance Next: Gravitational Redshift (RVs)
.
Last update: 06/07/2017 11:30 a.m. (CET)