Extinction (ebv, Av, & Rv)

Setup

Let's first make sure we have the latest version of PHOEBE 2.3 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.3,<2.4"

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()

Relevant Parameters

Extinction is parameterized by 3 parameters: ebv (E(B-V)), Av, and Rv. Of these three, two can be provided and the other must be constrained. By default, ebv is the constrained parameter. To change this, see the tutorial on constraints and the b.flip_constraint API docs.


In [3]:
print(b.filter(qualifier='ebv'))


ParameterSet: 2 parameters
C                      ebv@system: 0.0
                   ebv@constraint: {Av@system} / {Rv@system}

In [4]:
print(b.get_parameter(qualifier='ebv', context='system'))


Parameter: ebv@system
                       Qualifier: ebv
                     Description: Extinction E(B-V)
                           Value: 0.0
                  Constrained by: Av@system, Rv@system
                      Constrains: None
                      Related to: Av@system, Rv@system


In [5]:
print(b.get_parameter(qualifier='ebv', context='constraint'))


Constrains (qualifier): ebv
Expression in solar units (value): {Av@system} / {Rv@system}
Current Result (result): 0.0

In [6]:
print(b.get_parameter(qualifier='Av'))


Parameter: Av@system
                       Qualifier: Av
                     Description: Extinction Av
                           Value: 0.0
                  Constrained by: 
                      Constrains: ebv@system
                      Related to: Rv@system, ebv@system


In [7]:
print(b.get_parameter(qualifier='Rv'))


Parameter: Rv@system
                       Qualifier: Rv
                     Description: Extinction law parameter
                           Value: 3.1
                  Constrained by: 
                      Constrains: ebv@system
                      Related to: Av@system, ebv@system

For more details on the contribution of extinction to observables, see the following example scripts:


In [ ]: