tardis.montecarlo.montecarlo_numba.opacities module

tardis.montecarlo.montecarlo_numba.opacities.calculate_tau_electron(electron_density, distance)

Calculate tau for Thomson scattering

Parameters

electron_densityfloat

distancefloat

Returns

tau_electronfloat

tau for thomson scattering

tardis.montecarlo.montecarlo_numba.opacities.chi_bf_interpolator(opacity_state, nu, shell)

Interpolate the bound-free opacity.

This function interpolates the tabulated bound-free opacities and cross-sections to new frequency values nu.

Parameters

opacity_stateOpacityState

nufloat, dtype float

Comoving frequency of the r-packet.

shellint, dtype float

Current computational shell.

Returns

chi_bf_totfloat

Total bound-free opacity at frequency nu.

chi_bf_contributionsnumpy.ndarray, dtype float

Cumulative distribution function of the contributions of the individual bound free continua to the total bound-free opacity.

current_continuanumpy.ndarray, dtype int

Continuum ids for which absorption is possible for frequency nu.

x_sect_bfsnumpy.ndarray, dtype float

Photoionization cross-sections of all bound-free continua for which absorption is possible for frequency nu.

tardis.montecarlo.montecarlo_numba.opacities.chi_continuum_calculator(opacity_state, nu, shell)

Returns

chi_bf_totfloat

Total bound-free opacity at frequency nu.

chi_bf_contributionsnumpy.ndarray, dtype float

Cumulative distribution function of the contributions of the individual bound free continua to the total bound-free opacity.

current_continuanumpy.ndarray, dtype int

Continuum ids for which absorption is possible for frequency nu.

x_sect_bfsnumpy.ndarray, dtype float

Photoionization cross-sections of all bound-free continua for which absorption is possible for frequency nu.

chi_fffloat64

Free Free opacity at frequency nu

Attributes

opacity_stateOpacityState

nufloat64

Comoving frequency of the r_packet

shellint64

Current shell id of the r_packet

tardis.montecarlo.montecarlo_numba.opacities.chi_electron_calculator(opacity_state, nu, shell)

Calculate chi for Thomson scattering

Parameters

opacity_stateOpacityState

nufloat

Comoving frequency of the r-packet.

shellint

Shell of current r_packet

Returns

numpy.ndarray, dtype int

Continuum ids for which absorption is possible for frequency nu.

tardis.montecarlo.montecarlo_numba.opacities.chi_ff_calculator(opacity_state, nu, shell)

Returns

chi_fffloat64

Free Free opacity

Attributes

opacity_stateOpacityState

nufloat64

Comoving frequency of the r_packet

shellint64

Current shell id of the r_packet

tardis.montecarlo.montecarlo_numba.opacities.compton_opacity_calculation(energy, electron_density)

Calculate the Compton scattering opacity for a given energy (Rybicki & Lightman, 1979)

$ rho / 2 p_m times 3/4 sigma_T ((1 + kappa) / kappa^3 ((2kappa(1 + kappa)) / (1 + 2kappa) - ln(1 + 2kappa) + 1/(2kappa) ln(1 + 2kappa) - (1 + 3kappa) / (1 + 2kappa)^2) $

Parameters

energyfloat

The energy of the photon

ejecta_densityfloat

The density of the ejecta

Returns

float

The Compton scattering opacity

tardis.montecarlo.montecarlo_numba.opacities.compton_opacity_partial(energy, fraction)

Partial Compton scattering opacity, from artis file gamma.cc

Parameters

energyfloat

packet energy in terms of electron rest energy

fractionfloat

1 + 2 * packet energy

Returns

np.float64

Compton scattering opacity

tardis.montecarlo.montecarlo_numba.opacities.get_current_bound_free_continua(opacity_state, nu)

Determine bound-free continua for which absorption is possible.

Parameters

opacity_stateOpacityState

nufloat

Comoving frequency of the r-packet.

Returns

numpy.ndarray, dtype int

Continuum ids for which absorption is possible for frequency nu.

tardis.montecarlo.montecarlo_numba.opacities.kappa_calculation(energy)

Calculates kappa for various other calculations i.e. energy normalized to electron rest energy 511.0 KeV

Parameters

energyfloat

Returns

kappafloat

tardis.montecarlo.montecarlo_numba.opacities.pair_creation_opacity_artis(energy, ejecta_density, iron_group_fraction)

Calculates pair creation opacity for a given energy Approximate treatment from Ambwani & Sutherland (1988) as implemented in ARTIS

Parameters

energyfloat

Photon energy

ejecta_densityfloat

The density of the ejecta

iron_group_fractionfloat

Fraction of iron group elements in the shell

Returns

float

Pair creation opacity

tardis.montecarlo.montecarlo_numba.opacities.pair_creation_opacity_calculation(energy, ejecta_density, iron_group_fraction)

Calculates pair creation opacity for a given energy Approximate treatment from Ambwani & Sutherland (1988)

Parameters

energyfloat

Photon energy

ejecta_densityfloat

The density of the ejecta

iron_group_fractionfloat

Fraction of iron group elements in the shell

Returns

float

Pair creation opacity

tardis.montecarlo.montecarlo_numba.opacities.photoabsorption_opacity_calculation(energy, ejecta_density, iron_group_fraction)

Calculates photoabsorption opacity for a given energy Approximate treatment from Ambwani & Sutherland (1988) Magic numbers are from the approximate treatment

Parameters

energyfloat

Photon energy

ejecta_densityfloat

The density of the ejecta

iron_group_fractionfloat

Fraction of iron group elements in the shell

Returns

float

Photoabsorption opacity

tardis.montecarlo.montecarlo_numba.opacities.photoabsorption_opacity_calculation_kasen(energy, number_density, proton_count)

Calculates photoabsorption opacity for a given energy Approximate treatment from Kasen et al. (2006)

Parameters

energyfloat

Photon energy

number_densityfloat

The number density of the ejecta for each atom

proton_countfloat

Number of protons for each atom in the ejecta

Returns

float

Photoabsorption opacity