Beaming and Boosting

Setup

Let's first make sure we have the latest version of PHOEBE 2.0 installed. (You can comment out this line if you don't use pip for your installation or don't want to update to the latest release).



In [ ]:

    
!pip install -I "phoebe>=2.0,<2.1"

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



In [1]:

    
%matplotlib inline



In [2]:

    
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.



In [3]:

    
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.



In [4]:

    
times = np.linspace(0,1,101)



In [5]:

    
b.add_dataset('lc', times=times, dataset='lc01')









    Out[5]:





<ParameterSet: 15 parameters | contexts: compute, dataset>



In [6]:

    
b.add_dataset('rv', times=times, dataset='rv01')









    Out[6]:





<ParameterSet: 15 parameters | contexts: compute, dataset>



In [7]:

    
b.add_dataset('mesh', times=times[::10], dataset='mesh01')









    Out[7]:





<ParameterSet: 2 parameters | contexts: compute, dataset>

Relevant Parameters



In [8]:

    
b.set_value('irrad_method', 'none')



In [9]:

    
print b['boosting_method@compute']









    



Parameter: boosting_method@phoebe01@compute
                       Qualifier: boosting_method
                     Description: Type of boosting method
                           Value: none
                         Choices: none, linear



In [10]:

    
print b['boosting_method@compute'].choices









    



['none', 'linear']

Influence on Light Curves (fluxes)



In [11]:

    
b.run_compute(boosting_method='none', model='boosting_none')









    Out[11]:





<ParameterSet: 954 parameters | kinds: rv, mesh, lc>



In [12]:

    
b.run_compute(boosting_method='linear', model='boosting_linear')









    Out[12]:





<ParameterSet: 954 parameters | kinds: rv, mesh, lc>



In [13]:

    
axs, artists = b['lc01'].plot()
leg = plt.legend()



In [14]:

    
axs, artists = b['lc01'].plot(ylim=(1.01,1.03))
leg = plt.legend()

Influence on Radial Velocities



In [15]:

    
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)

Influence on Meshes



In [16]:

    
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]:



In [ ]: