PHOEBE 2.0 Documentation

2.0 Docs

  • 1.0
  • 2.0a
  • ver: 2.0


Prev: Reflection and Heating Next: Eclipse Detection
.

IPython Notebook | Python Script

Beaming and Boosting

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)

Let’s make our system so that the boosting effects will be quite noticeable.

b['rpole@primary'] = 1.8
b['rpole@secondary'] = 0.96

b['teff@primary'] = 10000
b['gravb_bol@primary'] = 1.0
b['teff@secondary'] = 5200
b['gravb_bol@secondary'] = 0.32

b['q@binary'] = 0.96/1.8
b['incl@binary'] = 88

b['period@binary'] = 1.0
b['sma@binary'] = 6.0
Fri, 10 Feb 2017 16:56 PARAMETERS   WARNING 'primary' probably has a radiative atm (teff=10000K>8000K), for which gravb_bol=1.00 might be a better approx than gravb_bol=0.32
Fri, 10 Feb 2017 16:56 PARAMETERS   WARNING misaligned orbits are not currently supported.
Fri, 10 Feb 2017 16:56 PARAMETERS   WARNING misaligned orbits are not currently supported.

We’ll add lc, rv, and mesh datasets so that we can see how they’re each affected by beaming and boosting.

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

Relevant Parameters

b.set_value('irrad_method', 'none')
print b['boosting_method@compute']
Parameter: boosting_method@phoebe01@compute
                       Qualifier: boosting_method
                     Description: Type of boosting method
                           Value: none
                         Choices: none, linear
print b['boosting_method@compute'].choices
['none', 'linear']

Influence on Light Curves (fluxes)

b.run_compute(boosting_method='none', model='boosting_none')
<ParameterSet: 954 parameters | kinds: rv, mesh, lc>
b.run_compute(boosting_method='linear', model='boosting_linear')
<ParameterSet: 954 parameters | kinds: rv, mesh, lc>
axs, artists = b['lc01'].plot()
leg = plt.legend()
../../_images/beaming_boosting_19_0.png
axs, artists = b['lc01'].plot(ylim=(1.01,1.03))
leg = plt.legend()
../../_images/beaming_boosting_20_0.png

Influence on Radial Velocities

fig = plt.figure(figsize=(10,6))
ax1, ax2 = fig.add_subplot(121), fig.add_subplot(122)

axs, artists = b['rv01@boosting_none'].plot(ax=ax1)
axs, artists = b['rv01@boosting_linear'].plot(ax=ax2)
../../_images/beaming_boosting_22_0.png

Influence on Meshes

fig = plt.figure(figsize=(10,6))
ax1, ax2 = fig.add_subplot(121), fig.add_subplot(122)

axs, artists = b['mesh@boosting_none'].plot(time=0.6, facecolor='boost_factors@lc01', edgecolor=None, ax=ax1)
axs, artists = b['mesh@boosting_linear'].plot(time=0.6, facecolor='boost_factors@lc01', edgecolor=None, ax=ax2)
/home/kyle/.local/lib/python2.7/site-packages/phoebe/frontend/plotting.py:167: RuntimeWarning: invalid value encountered in divide
  colorarray /= np.nanmax(colorarray)
/usr/local/lib/python2.7/dist-packages/matplotlib/colors.py:581: RuntimeWarning: invalid value encountered in less
  cbook._putmask(xa, xa < 0.0, -1)
/home/kyle/.local/lib/python2.7/site-packages/phoebe/frontend/plotting.py:257: FutureWarning: elementwise comparison failed; returning scalar instead, but in the future will perform elementwise comparison
  if pckwargs['edgecolors'] in ['none', 'None', None] and pckwargs['facecolors'] not in ['none', 'None', None]:
/home/kyle/.local/lib/python2.7/site-packages/phoebe/frontend/plotting.py:257: FutureWarning: comparison to None will result in an elementwise object comparison in the future.
  if pckwargs['edgecolors'] in ['none', 'None', None] and pckwargs['facecolors'] not in ['none', 'None', None]:
../../_images/beaming_boosting_24_1.png
Prev: Reflection and Heating Next: Eclipse Detection
.
Last update: 06/07/2017 11:30 a.m. (CET)