Source code for solrat.atom_model.multi_term_atom_model.data.MnI

from solrat.atom_model.multi_term_atom_model.object.level_registry import LevelRegistry
from solrat.atom_model.multi_term_atom_model.object.multi_term_atom_config import MultiTermAtomConfig
from solrat.atom_model.multi_term_atom_model.object.transition_registry import TransitionRegistry
from solrat.engine.functions.decorators import log_function_experimental




@log_function_experimental
def get_Mn_I_5432_config() -> MultiTermAtomConfig:  # pragma: no cover
    r"""
    Atomic model for Mn I 5432.5 A line, constrained to :math:`J=5/2 \to J=5/2` transition.

    Note that Mn I 5432.5 A transition is :math:`S=3.5 \to S=2.5`, i.e. it is forbidden
    under strict :math:`LS` coupling.
    Here it is modeled by expanding the upper term into two terms of different multiplicities.

    Note that Mn I 5432.5 A line also has significant HFS, which, however, can be modeled
    by artificially increasing the turbulent velocity (provided that the line width is
    large enough to begin with).

    :return: :any:`MultiTermAtomConfig` instance

    Due to rather crude approximations applied, this atomic model is recommended to be used with LTE SEE only.
    """
    # Levels
    level_registry = LevelRegistry()

    # a 6S (L=0, S=2.5) J=2.5
    level_registry.register_level(beta="a6S", L=0, S=2.5, J=2.5, energy_cmm1=0)

    # z 8P (L=1 S=3.5), J = 2.5, 3.5, 4.5
    # Within a rather crude single-J approximation, intercombinational transition can
    # be modeled by pretending that the upper term has only allowed component
    # |upper> = cos theta z6P + sin theta z8P.
    # Here, sin theta is dominant (determines the energy), but cos theta is the only one
    # allowed (determines selection rules).
    # Also need to multiply RTE by cos^2 theta (assuming SEE is LTE), or tau_line needs to be
    # redefined as tau_transition_under_interest
    level_registry.register_level(beta="z6P+z8P", L=1, S=2.5, J=1.5, energy_cmm1=18000)  # Artificial, not coupled
    level_registry.register_level(beta="z6P+z8P", L=1, S=2.5, J=2.5, energy_cmm1=18402.46)  # Energy from 8P
    level_registry.register_level(beta="z6P+z8P", L=1, S=2.5, J=3.5, energy_cmm1=18531.64)  # Energy from 8P

    level_registry.validate()

    # Transitions
    transition_registry = TransitionRegistry()

    transition_registry.register_transition(
        term_lower=level_registry.get_term(beta="a6S", L=0, S=2.5),
        term_upper=level_registry.get_term(beta="z6P+z8P", L=1, S=2.5),
        lower_J_constraint=[2.5],  # used if j_constrained=True
        upper_J_constraint=[2.5],  # used if j_constrained=True
        einstein_a_ul_sm1=6.04e03,
    )

    return MultiTermAtomConfig(
        level_registry=level_registry,
        transition_registry=transition_registry,
        reference_lambda_A_air=5432.55,
        atomic_mass_amu=54.9,
        j_constrained=True,
    )