Example 14 Electron and nuclear spin mixed state dynamics with RadicalPy

In this tutorial, one constructs MPO of radical pair system (two electron spins and a couple of nuclear spins under magnetic field) by using RadicalPy library.

Import required modules

try:
    import radicalpy
except ModuleNotFoundError:
    # remove uv if you are not using uv
    !uv pip install radicalpy --quiet

!uv pip show radicalpy

import pytdscf

pytdscf.__version__

Using Python 3.13.2 environment at /Users/hinom/GitHub/PyTDSCF-dev/.venv
Name: radicalpy
Version: 0.8.1
Location: /Users/hinom/GitHub/PyTDSCF-dev/.venv/lib/python3.13/site-packages
Requires: dot2tex, graphviz, importlib-resources, matplotlib, numpy, pint, scikit-learn, scipy, seaborn, sympy, tqdm
Required-by:

'1.1.0'

from math import isqrt

import matplotlib.pyplot as plt
import numpy as np
import radicalpy as rp
from IPython.display import HTML
from pympo import (
    AssignManager,
    OpSite,
    SumOfProducts,
)
from radicalpy.simulation import State
from sympy import Symbol

from pytdscf import BasInfo, Exciton, Model, Simulator, units
from pytdscf.dvr_operator_cls import TensorOperator
from pytdscf.hamiltonian_cls import TensorHamiltonian
from pytdscf.util import get_anim, read_nc

Vectorisation of density matrix

Density matrix \(\rho \in \mathbb{C}^{d \times d}\) can be “vectorised” as

\[|\rho\rangle\rangle := \mathrm{vec}(\rho) \in \mathbb{C}^{d^2}\]

where

\[\begin{split}\mathrm{vec}(\rho) = \begin{pmatrix} \rho_{11} \\ \rho_{21} \\ \vdots \\ \rho_{d1} \\ \rho_{12} \\ \vdots \\ \rho_{dd} \end{pmatrix}\end{split}\]

Thus, one-dimensional index (in the sense of column-major order, density matrix: ket->bra) is defined as

\[(ji) := (j-1)d + i\]

and

\[[\mathrm{vec}(\rho)]_{(ji)} = \rho_{ij}\]

If one adopt basis \(\{|i\rangle\}\),

\[|i\rangle\langle j| \to |i\rangle \otimes |j\rangle = |ij\rangle\rangle.\]

Now, dual operation \(A\rho B^\dagger\) can be vectorised as

\[\mathrm{vec}(A\rho B^\dagger) = (B^\ast \otimes A) \mathrm{vec}(\rho)\]

Proof >

\[[A\rho B^\dagger]_{ij} = \sum_{k,l} A_{ik} \rho_{kl} B^\ast_{jl}\]

>

\[\therefore \mathrm{vec}(A\rho B^\dagger)_{(ji)} = \sum_{k,l} (B_{jl}^\ast A_{ik}) \rho_{kl} = \{(B^\ast \otimes A) \mathrm{vec}(\rho) \}_{(ji)}\]

Vectorisation of the von Neumann equation

von Neumann equation is

\[i\hbar \frac{d}{dt} \rho = [H, \rho] = H\rho I^\dagger - I\rho H^\dagger \; (\because H=H^\dagger, I = I^\dagger)\]

Thus, its vectorised form is

\[i\hbar \frac{d}{dt} |\rho\rangle\rangle = (I \otimes H - H^\ast \otimes I) |\rho\rangle\rangle\]

Note!

Since, the density matrix is approximated by MPS, theare no guarantee that

• \(\rho = \rho^\dagger\) (Hermitian)

• \(\rho_{ij} \geq 0\) (positive semi-definite)

• \(\mathrm{Tr}(\rho) = 1\) (normalized)

Total Hamiltonian

\[\hat{H}_{\text{total}} = \hat{H}_{\text{Z}} + \hat{H}_{\text{H}} + \hat{H}_{\text{J}} + \hat{H}_{\text{D}}\]

Define systems

is_small_case = False

if is_small_case:
    # You can use following block instead
    n_nuc_spins = 1

    flavin = rp.simulation.Molecule.fromisotopes(
        isotopes=["1H"] * n_nuc_spins, hfcs=[0.4] * n_nuc_spins
    )
    Z = rp.simulation.Molecule.fromisotopes(
        isotopes=["1H"] * n_nuc_spins, hfcs=[0.5] * n_nuc_spins
    )
    sim = rp.simulation.LiouvilleSimulation([flavin, Z])

    # Parameters
    A = {}  # mT
    isotropic = True

    # Isotropic
    for i in range(len(sim.radicals)):
        for j, nuc in enumerate(sim.molecules[i].nuclei):
            if isotropic:
                A[(i, j)] = np.eye(3) * nuc.hfc.isotropic
            else:
                A[(i, j)] = nuc.hfc.anisotropic
    B0 = 0.2  # 2J
    B = np.array((0.0, 0.0, 1.0)) * B0  # mT
    J = 0.1  # Typically 1.0e+03 scale # mT
    D = 0.1  # mT
    kS = 1.0e06  # Exponential model in s-1
    kT = 1.0e06
    if isinstance(D, float):
        D = 2 / 3 * np.diag((-1.0, -1.0, 2.0)) * D * sim.radicals[0].gamma_mT
else:
    flavin = rp.simulation.Molecule.all_nuclei("flavin_anion")
    trp = rp.simulation.Molecule.all_nuclei("tryptophan_cation")
    sim = rp.simulation.LiouvilleSimulation([flavin, trp])
    A = {}
    isotropic = False
    for i in range(len(sim.radicals)):
        for j, nuc in enumerate(sim.molecules[i].nuclei):
            if isotropic:
                A[(i, j)] = np.eye(3) * nuc.hfc.isotropic
            else:
                A[(i, j)] = nuc.hfc.anisotropic
    B0 = 2.0  # Typically 0.01 mT~10 mT
    B = np.array((0.0, 0.0, 1.0)) * B0
    J = 0.1  # Typically 1.0e+03 scale

    D = (
        np.array(
            [
                [-0.225, 0.156, 0.198],
                [0.156, 0.117, -0.082],
                [0.198, -0.082, 0.107],
            ]
        )
        * sim.radicals[0].gamma_mT
    )
    kS = 1.0e06  # Exponential model in s-1
    kT = 1.0e06
    if isinstance(D, float):
        D = 2 / 3 * np.diag((-1.0, -1.0, 2.0)) * D * sim.radicals[0].gamma_mT
sim

Number of electrons: 2
Number of nuclei: 27
Number of particles: 29
Multiplicities: [2, 2, 3, 3, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 3, 3, 2, 2, 2, 2, 2, 2, 2, 2]
Magnetogyric ratios (mT): [-176085963.023, -176085963.023, 19337.792, 19337.792, 19337.792, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 19337.792, 267522.18744, 19337.792, 19337.792, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 267522.18744, 267522.18744]
Nuclei: [14N(19337792.0, 3, 0.5141 <anisotropic available>), 14N(19337792.0, 3, -0.001275 <anisotropic available>), 14N(19337792.0, 3, -0.03654 <anisotropic available>), 1H(267522187.44, 2, 0.05075 <anisotropic available>), 1H(267522187.44, 2, -0.1371 <anisotropic available>), 1H(267522187.44, 2, -0.1371 <anisotropic available>), 1H(267522187.44, 2, -0.1371 <anisotropic available>), 1H(267522187.44, 2, -0.4403 <anisotropic available>), 1H(267522187.44, 2, 0.4546 <anisotropic available>), 1H(267522187.44, 2, 0.4546 <anisotropic available>), 1H(267522187.44, 2, 0.4546 <anisotropic available>), 1H(267522187.44, 2, 0.009597 <anisotropic available>), 1H(267522187.44, 2, 0.4263 <anisotropic available>), 1H(267522187.44, 2, 0.4233 <anisotropic available>), 1H(267522187.44, 2, -0.02004 <anisotropic available>), 14N(19337792.0, 3, 0.1784 <anisotropic available>), 1H(267522187.44, 2, 1.605 <anisotropic available>), 14N(19337792.0, 3, 0.3216 <anisotropic available>), 14N(19337792.0, 3, 0.1465 <anisotropic available>), 1H(267522187.44, 2, 0.04567 <anisotropic available>), 1H(267522187.44, 2, -0.03213 <anisotropic available>), 1H(267522187.44, 2, -0.04 <anisotropic available>), 1H(267522187.44, 2, -0.09307 <anisotropic available>), 1H(267522187.44, 2, -0.278 <anisotropic available>), 1H(267522187.44, 2, -0.3634 <anisotropic available>), 1H(267522187.44, 2, -0.4879 <anisotropic available>), 1H(267522187.44, 2, -0.5983 <anisotropic available>)]
Couplings: [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
HFCs (mT): [0.5141 <anisotropic available>, -0.001275 <anisotropic available>, -0.03654 <anisotropic available>, 0.05075 <anisotropic available>, -0.1371 <anisotropic available>, -0.1371 <anisotropic available>, -0.1371 <anisotropic available>, -0.4403 <anisotropic available>, 0.4546 <anisotropic available>, 0.4546 <anisotropic available>, 0.4546 <anisotropic available>, 0.009597 <anisotropic available>, 0.4263 <anisotropic available>, 0.4233 <anisotropic available>, -0.02004 <anisotropic available>, 0.1784 <anisotropic available>, 1.605 <anisotropic available>, 0.3216 <anisotropic available>, 0.1465 <anisotropic available>, 0.04567 <anisotropic available>, -0.03213 <anisotropic available>, -0.04 <anisotropic available>, -0.09307 <anisotropic available>, -0.278 <anisotropic available>, -0.3634 <anisotropic available>, -0.4879 <anisotropic available>, -0.5983 <anisotropic available>]

Now, one defines matrix product state (MPS) in the following order

(nuclei in flavin) \(\to\) (electronic states \(\{|T_{+}\rangle, |T_{0}\rangle, |S\rangle, |T_{-}\rangle\}\)) \(\to\) (neclei in Z)

Extract one particle operator

RadicalPy provides variety of spin operators such as

• \(\hat{s}_x, \hat{s}_y, \hat{s}_z\) for radical singlet-triplet basis

• \(\hat{I}_x, \hat{I}_y, \hat{I}_z\) for nuclear Zeeman basis

# Clear nuclei temporally
_nuclei_tmp0 = sim.molecules[0].nuclei
_nuclei_tmp1 = sim.molecules[1].nuclei
sim.molecules[0].nuclei = []
sim.molecules[1].nuclei = []

# for Singlet-Triplet basis
sx_1 = sim.spin_operator(0, "x")
sy_1 = sim.spin_operator(0, "y")
sz_1 = sim.spin_operator(0, "z")
sx_2 = sim.spin_operator(1, "x")
sy_2 = sim.spin_operator(1, "y")
sz_2 = sim.spin_operator(1, "z")

Qs = sim.projection_operator(rp.simulation.State.SINGLET)
Qt = sim.projection_operator(rp.simulation.State.TRIPLET)

# Revert nuclei
sim.molecules[0].nuclei = _nuclei_tmp0
sim.molecules[1].nuclei = _nuclei_tmp1

plt.title("Projection operator $Q_S$")
plt.imshow(Qs.real)
plt.xticks([0, 1, 2, 3])
plt.yticks([0, 1, 2, 3])
plt.show()

Define OpSite and coefficients

• RadicalPy uses Hz in energy unit but it is too large to keep numerical stabiltiy

• Thus, one will use GHz in energy unit

• For some reasons, extraction of imaginary unit of pauli y operator leads stable simulation for degenerate system (why??)

SCALE = 1.0e-09
gamma = [p.gamma_mT for p in sim.particles]

g_ele_sym = [
    Symbol(r"\gamma_e^{(" + f"{i + 1}" + ")}") for i in range(len(sim.radicals))
]
g_nuc_sym = {}
for i in range(len(sim.radicals)):
    for j in range(len(sim.molecules[i].nuclei)):
        g_nuc_sym[(i, j)] = Symbol(r"\gamma_n^{" + f"{(i + 1, j + 1)}" + "}")

subs = {}
for i, ge in enumerate(g_ele_sym):
    subs[ge] = sim.radicals[i].gamma_mT
for (i, j), gn in g_nuc_sym.items():
    subs[gn] = sim.molecules[i].nuclei[j].gamma_mT

Define radical spin operators

def get_OE(op):
    """
    O* ⊗ 𝟙
    """
    return np.kron(op.conjugate(), np.eye(op.shape[0], dtype=op.dtype))


def get_EO(op):
    """
    𝟙 ⊗ O
    """
    return np.kron(np.eye(op.shape[0], dtype=op.dtype), op)


ele_site = len(sim.molecules[0].nuclei)

SxE_ops = []
SyE_ops = []
SzE_ops = []
ESx_ops = []
ESy_ops = []
ESz_ops = []

S1S2E_op = OpSite(
    r"(\hat{S}_1\cdot\hat{S}_2)^\ast ⊗ 𝟙",
    ele_site,
    value=get_OE(sx_1 @ sx_2 + sy_1 @ sy_2 + sz_1 @ sz_2),
)
ES1S2_op = OpSite(
    r"𝟙 ⊗ (\hat{S}_1\cdot\hat{S}_2)",
    ele_site,
    value=get_EO(sx_1 @ sx_2 + sy_1 @ sy_2 + sz_1 @ sz_2),
)
EE_op = OpSite(
    r"\hat{E} ⊗ \hat{E}", ele_site, value=get_OE(np.eye(sx_1.shape[0]))
)

QsE_op = OpSite(r"\hat{Q}_S ⊗ 𝟙", ele_site, value=get_OE(Qs))
EQs_op = OpSite(r"𝟙 ⊗ \hat{Q}_S", ele_site, value=get_EO(Qs))
QtE_op = OpSite(r"\hat{Q}_T ⊗ 𝟙", ele_site, value=get_OE(Qt))
EQt_op = OpSite(r"𝟙 ⊗ \hat{Q}_T", ele_site, value=get_EO(Qt))


SxE_ops.append(OpSite(r"\hat{S}_x^{(1)\ast} ⊗ 𝟙", ele_site, value=get_OE(sx_1)))
SxE_ops.append(OpSite(r"\hat{S}_x^{(2)\ast} ⊗ 𝟙", ele_site, value=get_OE(sx_2)))
SyE_ops.append(OpSite(r"\hat{S}_y^{(1)\ast} ⊗ 𝟙", ele_site, value=get_OE(sy_1)))
SyE_ops.append(OpSite(r"\hat{S}_y^{(2)\ast} ⊗ 𝟙", ele_site, value=get_OE(sy_2)))
SzE_ops.append(OpSite(r"\hat{S}_z^{(1)\ast} ⊗ 𝟙", ele_site, value=get_OE(sz_1)))
SzE_ops.append(OpSite(r"\hat{S}_z^{(2)\ast} ⊗ 𝟙", ele_site, value=get_OE(sz_2)))

ESx_ops.append(OpSite(r"𝟙 ⊗ \hat{S}_x^{(1)}", ele_site, value=get_EO(sx_1)))
ESx_ops.append(OpSite(r"𝟙 ⊗ \hat{S}_x^{(2)}", ele_site, value=get_EO(sx_2)))
ESy_ops.append(OpSite(r"𝟙 ⊗ \hat{S}_y^{(1)}", ele_site, value=get_EO(sy_1)))
ESy_ops.append(OpSite(r"𝟙 ⊗ \hat{S}_y^{(2)}", ele_site, value=get_EO(sy_2)))
ESz_ops.append(OpSite(r"𝟙 ⊗ \hat{S}_z^{(1)}", ele_site, value=get_EO(sz_1)))
ESz_ops.append(OpSite(r"𝟙 ⊗ \hat{S}_z^{(2)}", ele_site, value=get_EO(sz_2)))

SrE_ops = [SxE_ops, SyE_ops, SzE_ops]
ESr_ops = [ESx_ops, ESy_ops, ESz_ops]

Define nuclear spin operators

IxE_ops = {}
IyE_ops = {}
IzE_ops = {}
EIx_ops = {}
EIy_ops = {}
EIz_ops = {}

for j, nuc in enumerate(sim.molecules[0].nuclei):
    val = nuc.pauli["x"]
    IxE_ops[(0, j)] = OpSite(
        r"\hat{I}_x^{" + f"{(1, j + 1)}" + r"\ast} ⊗ 𝟙",
        j,
        value=get_OE(val),
    )
    EIx_ops[(0, j)] = OpSite(
        r"𝟙 ⊗ \hat{I}_x^{" + f"{(1, j + 1)}" + "}",
        j,
        value=get_EO(val),
    )
    val = nuc.pauli["y"]
    IyE_ops[(0, j)] = OpSite(
        r"\hat{I}_y^{" + f"{(1, j + 1)}" + r"\ast} ⊗ 𝟙",
        j,
        value=get_OE(val),
    )
    EIy_ops[(0, j)] = OpSite(
        r"𝟙 ⊗ \hat{I}_y^{" + f"{(1, j + 1)}" + "}",
        j,
        value=get_EO(val),
    )
    val = nuc.pauli["z"]
    IzE_ops[(0, j)] = OpSite(
        r"\hat{I}_z^{" + f"{(1, j + 1)}" + r"\ast} ⊗ 𝟙",
        j,
        value=get_OE(val),
    )
    EIz_ops[(0, j)] = OpSite(
        r"𝟙 ⊗ \hat{I}_z^{" + f"{(1, j + 1)}" + "}",
        j,
        value=get_EO(val),
    )

for j, nuc in enumerate(sim.molecules[1].nuclei):
    site = ele_site + 1 + j
    val = nuc.pauli["x"]
    IxE_ops[(1, j)] = OpSite(
        r"\hat{I}_x^{" + f"{(2, j + 1)}" + r"\ast} ⊗ 𝟙",
        site,
        value=get_OE(val),
    )
    EIx_ops[(1, j)] = OpSite(
        r"𝟙 ⊗ \hat{I}_x^{" + f"{(2, j + 1)}" + "}",
        site,
        value=get_EO(val),
    )
    val = nuc.pauli["y"]
    IyE_ops[(1, j)] = OpSite(
        r"\hat{I}_y^{" + f"{(2, j + 1)}" + r"\ast} ⊗ 𝟙",
        site,
        value=get_OE(val),
    )
    EIy_ops[(1, j)] = OpSite(
        r"𝟙 ⊗ \hat{I}_y^{" + f"{(2, j + 1)}" + "}",
        ele_site + 1 + j,
        value=get_EO(val),
    )
    val = nuc.pauli["z"]
    IzE_ops[(1, j)] = OpSite(
        r"\hat{I}_z^{" + f"{(2, j + 1)}" + r"\ast} ⊗ 𝟙",
        site,
        value=get_OE(val),
    )
    EIz_ops[(1, j)] = OpSite(
        r"𝟙 ⊗ \hat{I}_z^{" + f"{(2, j + 1)}" + "}",
        site,
        value=get_EO(val),
    )

IrE_ops = [IxE_ops, IyE_ops, IzE_ops]
EIr_ops = [EIx_ops, EIy_ops, EIz_ops]

Hyperfine coupling Hamiltonian

\[\hat{H}_{\text{H}} = \sum_i \sum_j \hat{S}_i\cdot A_{ij}\cdot \hat{I}_{ij} = \sum_i \sum_j \sum_{r\in\{x, y, z\}} A_{ij}\hat{S}_{r}^{(i)}\hat{I}_{r}^{(ij)}\]

hyperfine = SumOfProducts()

xyz = "xyz"

for i in range(len(sim.radicals)):
    for j in range(len(sim.molecules[i].nuclei)):
        for k, (SrE_op, ESr_op) in enumerate(
            zip(SrE_ops, ESr_ops, strict=True)
        ):
            for l, (IrE_op, EIr_op) in enumerate(
                zip(IrE_ops, EIr_ops, strict=True)
            ):
                if A[(i, j)][k, l] == 0.0:
                    continue
                Asym = Symbol(
                    "A^{"
                    + f"{(i + 1, j + 1)}"
                    + "}_{"
                    + f"{xyz[k]}"
                    + f"{xyz[l]}"
                    + "}"
                )
                subs[Asym] = A[(i, j)][k, l].item() * SCALE
                hyperfine -= Asym * g_ele_sym[i] * SrE_op[i] * IrE_op[(i, j)]
                hyperfine += Asym * g_ele_sym[i] * ESr_op[i] * EIr_op[(i, j)]

hyperfine = hyperfine.simplify()
hyperfine.symbol

$\displaystyle - A^{(1, 1)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 1)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 1)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 1)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 1)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 1)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 1)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 1)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 1)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 1)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 1)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 1)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 1)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 1)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 1)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 1)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 1)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 1)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 10)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 10)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 10)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 10)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 10)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 10)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 10)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 10)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 10)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 10)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 10)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 10)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 10)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 10)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 10)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 10)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 10)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 10)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 11)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 11)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 11)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 11)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 11)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 11)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 11)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 11)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 11)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 11)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 11)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 11)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 11)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 11)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 11)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 11)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 11)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 11)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 12)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 12)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 12)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 12)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 12)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 12)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 12)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 12)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 12)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 12)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 12)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 12)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 12)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 12)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 12)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 12)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 12)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 12)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 13)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 13)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 13)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 13)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 13)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 13)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 13)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 13)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 13)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 13)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 13)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 13)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 13)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 13)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 13)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 13)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 13)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 13)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 14)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 14)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 14)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 14)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 14)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 14)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 14)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 14)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 14)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 14)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 14)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 14)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 14)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 14)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 14)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 14)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 14)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 14)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 15)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 15)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 15)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 15)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 15)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 15)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 15)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 15)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 15)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 15)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 15)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 15)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 15)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 15)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 15)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 15)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 15)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 15)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 16)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 16)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 16)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 16)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 16)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 16)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 16)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 16)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 16)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 16)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 16)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 16)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 16)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 16)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 16)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 16)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 16)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 16)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 2)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 2)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 2)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 2)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 2)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 2)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 2)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 2)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 2)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 2)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 2)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 2)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 2)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 2)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 2)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 2)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 2)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 2)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 3)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 3)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 3)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 3)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 3)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 3)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 3)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 3)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 3)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 3)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 3)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 3)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 3)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 3)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 3)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 3)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 3)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 3)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 4)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 4)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 4)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 4)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 4)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 4)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 4)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 4)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 4)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 4)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 4)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 4)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 4)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 4)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 4)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 4)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 4)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 4)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 5)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 5)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 5)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 5)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 5)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 5)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 5)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 5)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 5)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 5)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 5)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 5)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 5)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 5)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 5)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 5)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 5)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 5)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 6)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 6)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 6)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 6)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 6)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 6)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 6)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 6)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 6)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 6)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 6)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 6)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 6)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 6)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 6)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 6)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 6)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 6)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 7)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 7)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 7)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 7)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 7)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 7)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 7)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 7)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 7)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 7)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 7)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 7)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 7)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 7)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 7)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 7)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 7)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 7)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 8)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 8)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 8)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 8)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 8)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 8)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 8)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 8)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 8)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 8)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 8)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 8)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 8)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 8)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 8)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 8)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 8)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 8)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 9)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 9)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 9)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 9)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 9)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 9)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 9)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 9)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 9)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 9)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 9)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 9)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 9)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 9)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 9)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 9)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 9)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 9)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(2, 1)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 1)} - A^{(2, 1)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 1)} - A^{(2, 1)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 1)} - A^{(2, 1)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 1)} - A^{(2, 1)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 1)} - A^{(2, 1)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 1)} - A^{(2, 1)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 1)} - A^{(2, 1)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 1)} - A^{(2, 1)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 1)} - A^{(2, 10)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 10)} - A^{(2, 10)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 10)} - A^{(2, 10)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 10)} - A^{(2, 10)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 10)} - A^{(2, 10)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 10)} - A^{(2, 10)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 10)} - A^{(2, 10)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 10)} - A^{(2, 10)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 10)} - A^{(2, 10)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 10)} - A^{(2, 11)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 11)} - A^{(2, 11)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 11)} - A^{(2, 11)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 11)} - A^{(2, 11)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 11)} - A^{(2, 11)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 11)} - A^{(2, 11)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 11)} - A^{(2, 11)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 11)} - A^{(2, 11)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 11)} - A^{(2, 11)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 11)} - A^{(2, 2)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 2)} - A^{(2, 2)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 2)} - A^{(2, 2)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 2)} - A^{(2, 2)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 2)} - A^{(2, 2)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 2)} - A^{(2, 2)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 2)} - A^{(2, 2)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 2)} - A^{(2, 2)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 2)} - A^{(2, 2)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 2)} - A^{(2, 3)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 3)} - A^{(2, 3)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 3)} - A^{(2, 3)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 3)} - A^{(2, 3)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 3)} - A^{(2, 3)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 3)} - A^{(2, 3)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 3)} - A^{(2, 3)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 3)} - A^{(2, 3)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 3)} - A^{(2, 3)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 3)} - A^{(2, 4)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 4)} - A^{(2, 4)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 4)} - A^{(2, 4)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 4)} - A^{(2, 4)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 4)} - A^{(2, 4)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 4)} - A^{(2, 4)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 4)} - A^{(2, 4)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 4)} - A^{(2, 4)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 4)} - A^{(2, 4)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 4)} - A^{(2, 5)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 5)} - A^{(2, 5)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 5)} - A^{(2, 5)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 5)} - A^{(2, 5)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 5)} - A^{(2, 5)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 5)} - A^{(2, 5)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 5)} - A^{(2, 5)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 5)} - A^{(2, 5)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 5)} - A^{(2, 5)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 5)} - A^{(2, 6)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 6)} - A^{(2, 6)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 6)} - A^{(2, 6)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 6)} - A^{(2, 6)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 6)} - A^{(2, 6)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 6)} - A^{(2, 6)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 6)} - A^{(2, 6)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 6)} - A^{(2, 6)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 6)} - A^{(2, 6)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 6)} - A^{(2, 7)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 7)} - A^{(2, 7)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 7)} - A^{(2, 7)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 7)} - A^{(2, 7)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 7)} - A^{(2, 7)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 7)} - A^{(2, 7)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 7)} - A^{(2, 7)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 7)} - A^{(2, 7)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 7)} - A^{(2, 7)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 7)} - A^{(2, 8)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 8)} - A^{(2, 8)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 8)} - A^{(2, 8)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 8)} - A^{(2, 8)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 8)} - A^{(2, 8)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 8)} - A^{(2, 8)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 8)} - A^{(2, 8)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 8)} - A^{(2, 8)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 8)} - A^{(2, 8)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 8)} - A^{(2, 9)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 9)} - A^{(2, 9)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 9)} - A^{(2, 9)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 9)} - A^{(2, 9)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 9)} - A^{(2, 9)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 9)} - A^{(2, 9)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 9)} - A^{(2, 9)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 9)} - A^{(2, 9)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 9)} - A^{(2, 9)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 9)}$

Zeeman Hamiltonian

\[\hat{H}_{\text{Z}} = - \vec{B} \cdot \sum_i \left(\gamma_{\text{e}}\vec{S}^{(i)} + \left(\sum_j\gamma_{\text{n}}^{(ij)}\vec{I}^{(ij)}\right)\right)\]

zeeman = SumOfProducts()

xyz = "xyz"

for k, (SrE_op, ESr_op, IrE_op, EIr_op) in enumerate(
    zip(SrE_ops, ESr_ops, IrE_ops, EIr_ops, strict=True)
):
    if B[k] == 0.0:
        continue
    r = xyz[k]
    Br = Symbol(f"B_{r}")
    subs[Br] = B[k] * SCALE
    for i in range(len(sim.radicals)):
        zeeman += -Br * g_ele_sym[i] * ESr_op[i]
        zeeman += Br * g_ele_sym[i] * SrE_op[i]
        for j in range(len(sim.molecules[i].nuclei)):
            zeeman += -Br * g_nuc_sym[(i, j)] * EIr_op[(i, j)]
            zeeman += Br * g_nuc_sym[(i, j)] * IrE_op[(i, j)]

zeeman = zeeman.simplify()
zeeman.symbol

$\displaystyle B_{z} \gamma_e^{(1)} \hat{S}_z^{(1)\ast} ⊗ 𝟙 - B_{z} \gamma_e^{(1)} 𝟙 ⊗ \hat{S}_z^{(1)} + B_{z} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 - B_{z} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} + B_{z} \gamma_n^{(1, 1)} \hat{I}_z^{(1, 1)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 1)} 𝟙 ⊗ \hat{I}_z^{(1, 1)} + B_{z} \gamma_n^{(1, 10)} \hat{I}_z^{(1, 10)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 10)} 𝟙 ⊗ \hat{I}_z^{(1, 10)} + B_{z} \gamma_n^{(1, 11)} \hat{I}_z^{(1, 11)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 11)} 𝟙 ⊗ \hat{I}_z^{(1, 11)} + B_{z} \gamma_n^{(1, 12)} \hat{I}_z^{(1, 12)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 12)} 𝟙 ⊗ \hat{I}_z^{(1, 12)} + B_{z} \gamma_n^{(1, 13)} \hat{I}_z^{(1, 13)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 13)} 𝟙 ⊗ \hat{I}_z^{(1, 13)} + B_{z} \gamma_n^{(1, 14)} \hat{I}_z^{(1, 14)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 14)} 𝟙 ⊗ \hat{I}_z^{(1, 14)} + B_{z} \gamma_n^{(1, 15)} \hat{I}_z^{(1, 15)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 15)} 𝟙 ⊗ \hat{I}_z^{(1, 15)} + B_{z} \gamma_n^{(1, 16)} \hat{I}_z^{(1, 16)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 16)} 𝟙 ⊗ \hat{I}_z^{(1, 16)} + B_{z} \gamma_n^{(1, 2)} \hat{I}_z^{(1, 2)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 2)} 𝟙 ⊗ \hat{I}_z^{(1, 2)} + B_{z} \gamma_n^{(1, 3)} \hat{I}_z^{(1, 3)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 3)} 𝟙 ⊗ \hat{I}_z^{(1, 3)} + B_{z} \gamma_n^{(1, 4)} \hat{I}_z^{(1, 4)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 4)} 𝟙 ⊗ \hat{I}_z^{(1, 4)} + B_{z} \gamma_n^{(1, 5)} \hat{I}_z^{(1, 5)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 5)} 𝟙 ⊗ \hat{I}_z^{(1, 5)} + B_{z} \gamma_n^{(1, 6)} \hat{I}_z^{(1, 6)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 6)} 𝟙 ⊗ \hat{I}_z^{(1, 6)} + B_{z} \gamma_n^{(1, 7)} \hat{I}_z^{(1, 7)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 7)} 𝟙 ⊗ \hat{I}_z^{(1, 7)} + B_{z} \gamma_n^{(1, 8)} \hat{I}_z^{(1, 8)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 8)} 𝟙 ⊗ \hat{I}_z^{(1, 8)} + B_{z} \gamma_n^{(1, 9)} \hat{I}_z^{(1, 9)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 9)} 𝟙 ⊗ \hat{I}_z^{(1, 9)} + B_{z} \gamma_n^{(2, 1)} \hat{I}_z^{(2, 1)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 1)} 𝟙 ⊗ \hat{I}_z^{(2, 1)} + B_{z} \gamma_n^{(2, 10)} \hat{I}_z^{(2, 10)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 10)} 𝟙 ⊗ \hat{I}_z^{(2, 10)} + B_{z} \gamma_n^{(2, 11)} \hat{I}_z^{(2, 11)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 11)} 𝟙 ⊗ \hat{I}_z^{(2, 11)} + B_{z} \gamma_n^{(2, 2)} \hat{I}_z^{(2, 2)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 2)} 𝟙 ⊗ \hat{I}_z^{(2, 2)} + B_{z} \gamma_n^{(2, 3)} \hat{I}_z^{(2, 3)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 3)} 𝟙 ⊗ \hat{I}_z^{(2, 3)} + B_{z} \gamma_n^{(2, 4)} \hat{I}_z^{(2, 4)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 4)} 𝟙 ⊗ \hat{I}_z^{(2, 4)} + B_{z} \gamma_n^{(2, 5)} \hat{I}_z^{(2, 5)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 5)} 𝟙 ⊗ \hat{I}_z^{(2, 5)} + B_{z} \gamma_n^{(2, 6)} \hat{I}_z^{(2, 6)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 6)} 𝟙 ⊗ \hat{I}_z^{(2, 6)} + B_{z} \gamma_n^{(2, 7)} \hat{I}_z^{(2, 7)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 7)} 𝟙 ⊗ \hat{I}_z^{(2, 7)} + B_{z} \gamma_n^{(2, 8)} \hat{I}_z^{(2, 8)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 8)} 𝟙 ⊗ \hat{I}_z^{(2, 8)} + B_{z} \gamma_n^{(2, 9)} \hat{I}_z^{(2, 9)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 9)} 𝟙 ⊗ \hat{I}_z^{(2, 9)}$

Exchange Hamiltonian

\[\hat{H}_{\text{J}} = J \gamma_e \left(\vec{S}_1\cdot\vec{S_2}+\frac{1}{2}\hat{E}\right)\]

exchange = SumOfProducts()
if J != 0.0:
    Jsym = Symbol("J")
    subs[Jsym] = J * SCALE
    exchange += Jsym * g_ele_sym[0] * (2 * ES1S2_op + 0.5 * EE_op)
    exchange -= Jsym * g_ele_sym[1] * (2 * S1S2E_op + 0.5 * EE_op)
    exchange = exchange.simplify()
exchange.symbol

$\displaystyle 0.5 J \gamma_e^{(1)} \hat{E} ⊗ \hat{E} + 2 J \gamma_e^{(1)} 𝟙 ⊗ (\hat{S}_1\cdot\hat{S}_2) - 2 J \gamma_e^{(2)} (\hat{S}_1\cdot\hat{S}_2)^\ast ⊗ 𝟙 - 0.5 J \gamma_e^{(2)} \hat{E} ⊗ \hat{E}$

Dipolar Hamiltonian

\[\hat{H}_{\text{J}} = \vec{S}_1^\top \cdot D \cdot \vec{S_2}\]

dipolar = SumOfProducts()
for k in range(3):
    for l in range(3):
        if D[k, l] == 0.0:
            continue
        else:
            Dsym = Symbol("D_{" + f"{xyz[k]}" + f"{xyz[l]}" + "}")
            subs[Dsym] = D[k, l] * SCALE
            dipolar += Dsym * ESr_ops[k][0] * ESr_ops[l][1]
            dipolar -= Dsym * SrE_ops[k][0] * SrE_ops[l][1]
dipolar = dipolar.simplify()
dipolar.symbol

$\displaystyle - D_{xx} \hat{S}_x^{(1)\ast} ⊗ 𝟙 \hat{S}_x^{(2)\ast} ⊗ 𝟙 + D_{xx} 𝟙 ⊗ \hat{S}_x^{(1)} 𝟙 ⊗ \hat{S}_x^{(2)} - D_{xy} \hat{S}_x^{(1)\ast} ⊗ 𝟙 \hat{S}_y^{(2)\ast} ⊗ 𝟙 + D_{xy} 𝟙 ⊗ \hat{S}_x^{(1)} 𝟙 ⊗ \hat{S}_y^{(2)} - D_{xz} \hat{S}_x^{(1)\ast} ⊗ 𝟙 \hat{S}_z^{(2)\ast} ⊗ 𝟙 + D_{xz} 𝟙 ⊗ \hat{S}_x^{(1)} 𝟙 ⊗ \hat{S}_z^{(2)} - D_{yx} \hat{S}_y^{(1)\ast} ⊗ 𝟙 \hat{S}_x^{(2)\ast} ⊗ 𝟙 + D_{yx} 𝟙 ⊗ \hat{S}_y^{(1)} 𝟙 ⊗ \hat{S}_x^{(2)} - D_{yy} \hat{S}_y^{(1)\ast} ⊗ 𝟙 \hat{S}_y^{(2)\ast} ⊗ 𝟙 + D_{yy} 𝟙 ⊗ \hat{S}_y^{(1)} 𝟙 ⊗ \hat{S}_y^{(2)} - D_{yz} \hat{S}_y^{(1)\ast} ⊗ 𝟙 \hat{S}_z^{(2)\ast} ⊗ 𝟙 + D_{yz} 𝟙 ⊗ \hat{S}_y^{(1)} 𝟙 ⊗ \hat{S}_z^{(2)} - D_{zx} \hat{S}_z^{(1)\ast} ⊗ 𝟙 \hat{S}_x^{(2)\ast} ⊗ 𝟙 + D_{zx} 𝟙 ⊗ \hat{S}_z^{(1)} 𝟙 ⊗ \hat{S}_x^{(2)} - D_{zy} \hat{S}_z^{(1)\ast} ⊗ 𝟙 \hat{S}_y^{(2)\ast} ⊗ 𝟙 + D_{zy} 𝟙 ⊗ \hat{S}_z^{(1)} 𝟙 ⊗ \hat{S}_y^{(2)} - D_{zz} \hat{S}_z^{(1)\ast} ⊗ 𝟙 \hat{S}_z^{(2)\ast} ⊗ 𝟙 + D_{zz} 𝟙 ⊗ \hat{S}_z^{(1)} 𝟙 ⊗ \hat{S}_z^{(2)}$

Haberkorn relaxation

\[i\hbar \frac{d}{dt} \rho = [H, \rho] - \sum_{i\in\{S, T\}}\frac{i\hbar k_i}{2}\{Q_i, \rho\}\]

The superoperator for vectorised form density matrix is

\[Q_i \otimes \mathbb{1} + \mathbb{1} \otimes Q_i\]

haberkorn = SumOfProducts()
if kS != 0.0:
    kSsym = Symbol("k_S")
    subs[kSsym] = kS * SCALE
    haberkorn -= 1.0j * kSsym / 2 * (QsE_op + EQs_op)
if kT != 0.0:
    kTsym = Symbol("k_T")
    subs[kTsym] = kT * SCALE
    haberkorn -= 1.0j * kTsym / 2 * (QtE_op + EQt_op)
haberkorn = haberkorn.simplify()
haberkorn.symbol

$\displaystyle - 0.5 i k_{S} \left(\hat{Q}_S ⊗ 𝟙 + 𝟙 ⊗ \hat{Q}_S\right) - 0.5 i k_{T} \left(\hat{Q}_T ⊗ 𝟙 + 𝟙 ⊗ \hat{Q}_T\right)$

Construct matrix product operator (MPO) for MPS propagation

L = hyperfine + zeeman + exchange + dipolar + haberkorn
L = L.simplify()
am = AssignManager(L)
_ = am.assign()
# display(*am.Wsym)
mpo = am.numerical_mpo(subs=subs)
L.symbol

$\displaystyle - A^{(1, 1)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 1)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 1)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 1)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 1)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 1)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 1)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 1)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 1)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 1)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 1)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 1)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 1)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 1)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 1)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 1)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 1)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 1)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 1)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 1)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 10)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 10)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 10)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 10)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 10)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 10)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 10)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 10)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 10)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 10)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 10)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 10)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 10)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 10)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 10)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 10)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 10)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 10)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 10)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 10)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 11)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 11)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 11)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 11)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 11)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 11)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 11)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 11)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 11)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 11)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 11)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 11)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 11)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 11)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 11)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 11)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 11)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 11)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 11)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 11)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 12)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 12)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 12)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 12)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 12)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 12)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 12)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 12)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 12)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 12)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 12)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 12)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 12)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 12)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 12)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 12)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 12)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 12)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 12)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 12)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 13)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 13)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 13)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 13)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 13)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 13)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 13)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 13)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 13)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 13)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 13)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 13)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 13)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 13)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 13)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 13)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 13)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 13)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 13)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 13)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 14)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 14)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 14)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 14)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 14)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 14)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 14)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 14)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 14)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 14)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 14)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 14)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 14)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 14)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 14)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 14)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 14)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 14)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 14)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 14)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 15)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 15)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 15)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 15)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 15)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 15)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 15)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 15)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 15)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 15)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 15)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 15)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 15)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 15)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 15)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 15)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 15)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 15)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 15)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 15)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 16)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 16)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 16)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 16)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 16)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 16)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 16)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 16)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 16)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 16)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 16)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 16)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 16)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 16)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 16)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 16)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 16)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 16)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 16)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 16)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 2)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 2)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 2)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 2)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 2)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 2)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 2)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 2)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 2)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 2)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 2)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 2)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 2)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 2)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 2)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 2)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 2)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 2)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 2)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 2)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 3)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 3)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 3)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 3)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 3)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 3)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 3)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 3)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 3)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 3)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 3)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 3)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 3)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 3)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 3)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 3)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 3)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 3)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 3)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 3)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 4)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 4)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 4)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 4)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 4)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 4)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 4)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 4)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 4)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 4)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 4)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 4)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 4)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 4)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 4)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 4)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 4)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 4)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 4)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 4)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 5)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 5)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 5)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 5)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 5)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 5)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 5)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 5)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 5)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 5)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 5)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 5)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 5)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 5)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 5)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 5)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 5)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 5)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 5)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 5)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 6)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 6)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 6)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 6)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 6)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 6)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 6)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 6)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 6)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 6)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 6)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 6)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 6)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 6)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 6)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 6)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 6)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 6)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 6)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 6)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 7)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 7)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 7)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 7)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 7)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 7)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 7)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 7)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 7)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 7)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 7)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 7)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 7)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 7)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 7)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 7)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 7)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 7)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 7)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 7)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 8)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 8)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 8)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 8)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 8)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 8)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 8)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 8)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 8)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 8)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 8)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 8)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 8)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 8)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 8)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 8)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 8)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 8)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 8)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 8)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 9)}_{xx} \gamma_e^{(1)} \hat{I}_x^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{xx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 9)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 9)}_{xy} \gamma_e^{(1)} \hat{I}_y^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{xy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 9)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 9)}_{xz} \gamma_e^{(1)} \hat{I}_z^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_x^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{xz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 9)} 𝟙 ⊗ \hat{S}_x^{(1)} - A^{(1, 9)}_{yx} \gamma_e^{(1)} \hat{I}_x^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{yx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 9)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 9)}_{yy} \gamma_e^{(1)} \hat{I}_y^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{yy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 9)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 9)}_{yz} \gamma_e^{(1)} \hat{I}_z^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_y^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{yz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 9)} 𝟙 ⊗ \hat{S}_y^{(1)} - A^{(1, 9)}_{zx} \gamma_e^{(1)} \hat{I}_x^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{zx} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_x^{(1, 9)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 9)}_{zy} \gamma_e^{(1)} \hat{I}_y^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{zy} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_y^{(1, 9)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(1, 9)}_{zz} \gamma_e^{(1)} \hat{I}_z^{(1, 9)\ast} ⊗ 𝟙 \hat{S}_z^{(1)\ast} ⊗ 𝟙 + A^{(1, 9)}_{zz} \gamma_e^{(1)} 𝟙 ⊗ \hat{I}_z^{(1, 9)} 𝟙 ⊗ \hat{S}_z^{(1)} - A^{(2, 1)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 1)} - A^{(2, 1)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 1)} - A^{(2, 1)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 1)} - A^{(2, 1)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 1)} - A^{(2, 1)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 1)} - A^{(2, 1)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 1)} - A^{(2, 1)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 1)} - A^{(2, 1)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 1)} - A^{(2, 1)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 1)\ast} ⊗ 𝟙 + A^{(2, 1)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 1)} - A^{(2, 10)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 10)} - A^{(2, 10)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 10)} - A^{(2, 10)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 10)} - A^{(2, 10)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 10)} - A^{(2, 10)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 10)} - A^{(2, 10)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 10)} - A^{(2, 10)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 10)} - A^{(2, 10)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 10)} - A^{(2, 10)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 10)\ast} ⊗ 𝟙 + A^{(2, 10)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 10)} - A^{(2, 11)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 11)} - A^{(2, 11)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 11)} - A^{(2, 11)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 11)} - A^{(2, 11)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 11)} - A^{(2, 11)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 11)} - A^{(2, 11)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 11)} - A^{(2, 11)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 11)} - A^{(2, 11)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 11)} - A^{(2, 11)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 11)\ast} ⊗ 𝟙 + A^{(2, 11)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 11)} - A^{(2, 2)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 2)} - A^{(2, 2)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 2)} - A^{(2, 2)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 2)} - A^{(2, 2)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 2)} - A^{(2, 2)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 2)} - A^{(2, 2)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 2)} - A^{(2, 2)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 2)} - A^{(2, 2)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 2)} - A^{(2, 2)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 2)\ast} ⊗ 𝟙 + A^{(2, 2)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 2)} - A^{(2, 3)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 3)} - A^{(2, 3)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 3)} - A^{(2, 3)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 3)} - A^{(2, 3)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 3)} - A^{(2, 3)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 3)} - A^{(2, 3)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 3)} - A^{(2, 3)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 3)} - A^{(2, 3)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 3)} - A^{(2, 3)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 3)\ast} ⊗ 𝟙 + A^{(2, 3)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 3)} - A^{(2, 4)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 4)} - A^{(2, 4)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 4)} - A^{(2, 4)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 4)} - A^{(2, 4)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 4)} - A^{(2, 4)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 4)} - A^{(2, 4)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 4)} - A^{(2, 4)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 4)} - A^{(2, 4)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 4)} - A^{(2, 4)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 4)\ast} ⊗ 𝟙 + A^{(2, 4)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 4)} - A^{(2, 5)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 5)} - A^{(2, 5)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 5)} - A^{(2, 5)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 5)} - A^{(2, 5)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 5)} - A^{(2, 5)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 5)} - A^{(2, 5)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 5)} - A^{(2, 5)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 5)} - A^{(2, 5)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 5)} - A^{(2, 5)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 5)\ast} ⊗ 𝟙 + A^{(2, 5)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 5)} - A^{(2, 6)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 6)} - A^{(2, 6)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 6)} - A^{(2, 6)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 6)} - A^{(2, 6)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 6)} - A^{(2, 6)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 6)} - A^{(2, 6)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 6)} - A^{(2, 6)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 6)} - A^{(2, 6)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 6)} - A^{(2, 6)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 6)\ast} ⊗ 𝟙 + A^{(2, 6)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 6)} - A^{(2, 7)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 7)} - A^{(2, 7)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 7)} - A^{(2, 7)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 7)} - A^{(2, 7)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 7)} - A^{(2, 7)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 7)} - A^{(2, 7)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 7)} - A^{(2, 7)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 7)} - A^{(2, 7)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 7)} - A^{(2, 7)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 7)\ast} ⊗ 𝟙 + A^{(2, 7)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 7)} - A^{(2, 8)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 8)} - A^{(2, 8)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 8)} - A^{(2, 8)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 8)} - A^{(2, 8)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 8)} - A^{(2, 8)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 8)} - A^{(2, 8)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 8)} - A^{(2, 8)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 8)} - A^{(2, 8)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 8)} - A^{(2, 8)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 8)\ast} ⊗ 𝟙 + A^{(2, 8)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 8)} - A^{(2, 9)}_{xx} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{xx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 9)} - A^{(2, 9)}_{xy} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{xy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 9)} - A^{(2, 9)}_{xz} \gamma_e^{(2)} \hat{S}_x^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{xz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_x^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 9)} - A^{(2, 9)}_{yx} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{yx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 9)} - A^{(2, 9)}_{yy} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{yy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 9)} - A^{(2, 9)}_{yz} \gamma_e^{(2)} \hat{S}_y^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{yz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_y^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 9)} - A^{(2, 9)}_{zx} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_x^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{zx} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_x^{(2, 9)} - A^{(2, 9)}_{zy} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_y^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{zy} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_y^{(2, 9)} - A^{(2, 9)}_{zz} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 \hat{I}_z^{(2, 9)\ast} ⊗ 𝟙 + A^{(2, 9)}_{zz} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} 𝟙 ⊗ \hat{I}_z^{(2, 9)} + B_{z} \gamma_e^{(1)} \hat{S}_z^{(1)\ast} ⊗ 𝟙 - B_{z} \gamma_e^{(1)} 𝟙 ⊗ \hat{S}_z^{(1)} + B_{z} \gamma_e^{(2)} \hat{S}_z^{(2)\ast} ⊗ 𝟙 - B_{z} \gamma_e^{(2)} 𝟙 ⊗ \hat{S}_z^{(2)} + B_{z} \gamma_n^{(1, 1)} \hat{I}_z^{(1, 1)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 1)} 𝟙 ⊗ \hat{I}_z^{(1, 1)} + B_{z} \gamma_n^{(1, 10)} \hat{I}_z^{(1, 10)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 10)} 𝟙 ⊗ \hat{I}_z^{(1, 10)} + B_{z} \gamma_n^{(1, 11)} \hat{I}_z^{(1, 11)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 11)} 𝟙 ⊗ \hat{I}_z^{(1, 11)} + B_{z} \gamma_n^{(1, 12)} \hat{I}_z^{(1, 12)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 12)} 𝟙 ⊗ \hat{I}_z^{(1, 12)} + B_{z} \gamma_n^{(1, 13)} \hat{I}_z^{(1, 13)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 13)} 𝟙 ⊗ \hat{I}_z^{(1, 13)} + B_{z} \gamma_n^{(1, 14)} \hat{I}_z^{(1, 14)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 14)} 𝟙 ⊗ \hat{I}_z^{(1, 14)} + B_{z} \gamma_n^{(1, 15)} \hat{I}_z^{(1, 15)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 15)} 𝟙 ⊗ \hat{I}_z^{(1, 15)} + B_{z} \gamma_n^{(1, 16)} \hat{I}_z^{(1, 16)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 16)} 𝟙 ⊗ \hat{I}_z^{(1, 16)} + B_{z} \gamma_n^{(1, 2)} \hat{I}_z^{(1, 2)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 2)} 𝟙 ⊗ \hat{I}_z^{(1, 2)} + B_{z} \gamma_n^{(1, 3)} \hat{I}_z^{(1, 3)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 3)} 𝟙 ⊗ \hat{I}_z^{(1, 3)} + B_{z} \gamma_n^{(1, 4)} \hat{I}_z^{(1, 4)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 4)} 𝟙 ⊗ \hat{I}_z^{(1, 4)} + B_{z} \gamma_n^{(1, 5)} \hat{I}_z^{(1, 5)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 5)} 𝟙 ⊗ \hat{I}_z^{(1, 5)} + B_{z} \gamma_n^{(1, 6)} \hat{I}_z^{(1, 6)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 6)} 𝟙 ⊗ \hat{I}_z^{(1, 6)} + B_{z} \gamma_n^{(1, 7)} \hat{I}_z^{(1, 7)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 7)} 𝟙 ⊗ \hat{I}_z^{(1, 7)} + B_{z} \gamma_n^{(1, 8)} \hat{I}_z^{(1, 8)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 8)} 𝟙 ⊗ \hat{I}_z^{(1, 8)} + B_{z} \gamma_n^{(1, 9)} \hat{I}_z^{(1, 9)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(1, 9)} 𝟙 ⊗ \hat{I}_z^{(1, 9)} + B_{z} \gamma_n^{(2, 1)} \hat{I}_z^{(2, 1)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 1)} 𝟙 ⊗ \hat{I}_z^{(2, 1)} + B_{z} \gamma_n^{(2, 10)} \hat{I}_z^{(2, 10)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 10)} 𝟙 ⊗ \hat{I}_z^{(2, 10)} + B_{z} \gamma_n^{(2, 11)} \hat{I}_z^{(2, 11)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 11)} 𝟙 ⊗ \hat{I}_z^{(2, 11)} + B_{z} \gamma_n^{(2, 2)} \hat{I}_z^{(2, 2)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 2)} 𝟙 ⊗ \hat{I}_z^{(2, 2)} + B_{z} \gamma_n^{(2, 3)} \hat{I}_z^{(2, 3)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 3)} 𝟙 ⊗ \hat{I}_z^{(2, 3)} + B_{z} \gamma_n^{(2, 4)} \hat{I}_z^{(2, 4)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 4)} 𝟙 ⊗ \hat{I}_z^{(2, 4)} + B_{z} \gamma_n^{(2, 5)} \hat{I}_z^{(2, 5)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 5)} 𝟙 ⊗ \hat{I}_z^{(2, 5)} + B_{z} \gamma_n^{(2, 6)} \hat{I}_z^{(2, 6)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 6)} 𝟙 ⊗ \hat{I}_z^{(2, 6)} + B_{z} \gamma_n^{(2, 7)} \hat{I}_z^{(2, 7)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 7)} 𝟙 ⊗ \hat{I}_z^{(2, 7)} + B_{z} \gamma_n^{(2, 8)} \hat{I}_z^{(2, 8)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 8)} 𝟙 ⊗ \hat{I}_z^{(2, 8)} + B_{z} \gamma_n^{(2, 9)} \hat{I}_z^{(2, 9)\ast} ⊗ 𝟙 - B_{z} \gamma_n^{(2, 9)} 𝟙 ⊗ \hat{I}_z^{(2, 9)} - D_{xx} \hat{S}_x^{(1)\ast} ⊗ 𝟙 \hat{S}_x^{(2)\ast} ⊗ 𝟙 + D_{xx} 𝟙 ⊗ \hat{S}_x^{(1)} 𝟙 ⊗ \hat{S}_x^{(2)} - D_{xy} \hat{S}_x^{(1)\ast} ⊗ 𝟙 \hat{S}_y^{(2)\ast} ⊗ 𝟙 + D_{xy} 𝟙 ⊗ \hat{S}_x^{(1)} 𝟙 ⊗ \hat{S}_y^{(2)} - D_{xz} \hat{S}_x^{(1)\ast} ⊗ 𝟙 \hat{S}_z^{(2)\ast} ⊗ 𝟙 + D_{xz} 𝟙 ⊗ \hat{S}_x^{(1)} 𝟙 ⊗ \hat{S}_z^{(2)} - D_{yx} \hat{S}_y^{(1)\ast} ⊗ 𝟙 \hat{S}_x^{(2)\ast} ⊗ 𝟙 + D_{yx} 𝟙 ⊗ \hat{S}_y^{(1)} 𝟙 ⊗ \hat{S}_x^{(2)} - D_{yy} \hat{S}_y^{(1)\ast} ⊗ 𝟙 \hat{S}_y^{(2)\ast} ⊗ 𝟙 + D_{yy} 𝟙 ⊗ \hat{S}_y^{(1)} 𝟙 ⊗ \hat{S}_y^{(2)} - D_{yz} \hat{S}_y^{(1)\ast} ⊗ 𝟙 \hat{S}_z^{(2)\ast} ⊗ 𝟙 + D_{yz} 𝟙 ⊗ \hat{S}_y^{(1)} 𝟙 ⊗ \hat{S}_z^{(2)} - D_{zx} \hat{S}_z^{(1)\ast} ⊗ 𝟙 \hat{S}_x^{(2)\ast} ⊗ 𝟙 + D_{zx} 𝟙 ⊗ \hat{S}_z^{(1)} 𝟙 ⊗ \hat{S}_x^{(2)} - D_{zy} \hat{S}_z^{(1)\ast} ⊗ 𝟙 \hat{S}_y^{(2)\ast} ⊗ 𝟙 + D_{zy} 𝟙 ⊗ \hat{S}_z^{(1)} 𝟙 ⊗ \hat{S}_y^{(2)} - D_{zz} \hat{S}_z^{(1)\ast} ⊗ 𝟙 \hat{S}_z^{(2)\ast} ⊗ 𝟙 + D_{zz} 𝟙 ⊗ \hat{S}_z^{(1)} 𝟙 ⊗ \hat{S}_z^{(2)} + 0.5 J \gamma_e^{(1)} \hat{E} ⊗ \hat{E} + 2 J \gamma_e^{(1)} 𝟙 ⊗ (\hat{S}_1\cdot\hat{S}_2) - 2 J \gamma_e^{(2)} (\hat{S}_1\cdot\hat{S}_2)^\ast ⊗ 𝟙 - 0.5 J \gamma_e^{(2)} \hat{E} ⊗ \hat{E} - 0.5 i k_{S} \left(\hat{Q}_S ⊗ 𝟙 + 𝟙 ⊗ \hat{Q}_S\right) - 0.5 i k_{T} \left(\hat{Q}_T ⊗ 𝟙 + 𝟙 ⊗ \hat{Q}_T\right)$

Liouville simulation

backend = "numpy"  # or "jax"  for large bond dimension and GPU
m = 48  # Bond dimension
Δt = 5.0e-10 / SCALE * units.au_in_fs  # dt = 0.5 ns

basis = []
for nuc in sim.molecules[0].nuclei:
    basis.append(Exciton(nstate=nuc.multiplicity**2))
basis.append(Exciton(nstate=4**2))
for nuc in sim.molecules[1].nuclei:
    basis.append(Exciton(nstate=nuc.multiplicity**2))
basinfo = BasInfo([basis], spf_info=None)

nsite = len(basis)

op_dict = {
    tuple([(isite, isite) for isite in range(nsite)]): TensorOperator(mpo=mpo)
}
H = TensorHamiltonian(
    nsite, potential=[[op_dict]], kinetic=None, backend=backend
)
op_dict = {
    ((ele_site, ele_site),) : TensorOperator(mpo=[Qs[None, :, : , None]],
                                         legs=(ele_site, ele_site))
}
singlet = TensorHamiltonian(
    nsite, potential=[[op_dict]], kinetic=None, backend=backend
)

def singlet_state():
    hp = []
    for isite in range(nsite):
        if isite == ele_site:
            op = Qs
        else:
            # Mixed states !
            op = np.eye(isqrt(basis[isite].nstate))
        # Automatically nomarized so that trace=1 in internal code
        hp.append(op.reshape(-1).tolist())
    return hp


def process(H):
    operators = {"hamiltonian": H, "S": singlet}
    model = Model(
        basinfo=basinfo, operators=operators, space="Liouville"
    )  # <- changed!
    model.m_aux_max = m
    hp = singlet_state()
    model.init_HartreeProduct = [hp]

    jobname = "radicalpair_liouville"
    simulator = Simulator(
        jobname=jobname,
        model=model,
        backend=backend,
        verbose=0,
    )
    # Initiate the propagation setting with maxstep=0
    nstep = 100
    ener, wf = simulator.propagate(
        reduced_density=(
            [(ele_site, ele_site)],
            2,
        ),
        maxstep=nstep,
        stepsize=Δt,
        autocorr=False,
        energy=False,
        norm=False,
        populations=False,
        observables=True,
        integrator="arnoldi",
    )
    data = read_nc(f"{jobname}_prop/reduced_density.nc", [(ele_site, ele_site)])
    time_data = data["time"]
    density_data = data[(ele_site, ele_site)]
    return density_data, time_data

Since nuclear spins are mixed state, one needs to consider ensemble of nuclear spin initial states.

dm, time_data = process(H)
time_data_μs = time_data * SCALE * 1e06 / units.au_in_fs

14:28:32 | INFO | Wave function is saved in wf_radicalpair_liouville.pkl
14:28:32 | INFO | Start initial step    0.000 [fs]

14:28:33 | INFO | End     0 step; propagated    0.012 [fs]; AVG Krylov iteration: 6.04
14:30:20 | INFO | End    99 step; propagated    1.197 [fs]; AVG Krylov iteration: 5.86
14:30:20 | INFO | End simulation and save wavefunction
14:30:20 | INFO | Wave function is saved in wf_radicalpair_liouville.pkl

plt.plot(time_data_μs, dm[:, 0, 0].real, linewidth=2, label="T+")
plt.plot(time_data_μs, dm[:, 1, 1].real, linewidth=2, label="T0")
plt.plot(time_data_μs, dm[:, 2, 2].real, linewidth=2, label="S")
plt.plot(time_data_μs, dm[:, 3, 3].real, linewidth=2, label="T-")
plt.plot(time_data_μs, np.einsum("ijj->i", dm.real), linewidth=2, label="trace")
plt.title(f"MPS: {m=} B = {B0} mT J = {J} mT")
plt.legend()
plt.xlabel("Time [μs]")
plt.ylabel("Popultation")
plt.grid()
# plt.ylim([0, 1])
plt.savefig("population.png")
plt.show()

from scipy.linalg import logm
plt.plot(time_data_μs, [-np.trace(rho * logm(rho)).real for rho in dm], linewidth=2, label="")

plt.title(f"MPS: {m=} B = {B0} mT J = {J} mT")
plt.legend()
plt.xlabel("Time [μs]")
plt.ylabel(r"Entanglement entropy -Tr($\rho \ln \rho$)")
plt.grid()
plt.show()

/Users/hinom/GitHub/PyTDSCF-dev/.venv/lib/python3.13/site-packages/scipy/linalg/_matfuncs.py:202: LogmExactlySingularWarning: The logm input matrix is exactly singular.
  F = scipy.linalg._matfuncs_inv_ssq._logm(A)
/var/folders/sj/nm7r7k8969b31b686cw52trh0000gn/T/ipykernel_1570/1584447722.py:5: UserWarning: No artists with labels found to put in legend.  Note that artists whose label start with an underscore are ignored when legend() is called with no argument.
  plt.legend()

dt = time_data_μs[1] - time_data[0]
phi_s = np.cumsum(dm[:, 2, 2].real * dt * kS * np.exp(-kS * time_data_μs))

plt.plot(time_data_μs, phi_s)
plt.xlabel("Time [μs]")
plt.ylabel(r"$\Phi_S(t)=\int_0^t dt \langle P^S \rangle(t) ke^{-kt} $")
plt.grid()
plt.savefig("population.png")
plt.show()

# Note: spin coherent state sampling has slower convergence for off-diagnal term than projection state sampling
fig, anim = get_anim(
    dm[::2, :, :],
    time_data_μs[::2],
    title="Reduced density matrix",
    time_unit="μs",
    save_gif=True,
    dpi=30,
    gif_filename="rdm-radical.gif",
    row_names=[
        r"$|T+\rangle$",
        r"$|T0\rangle$",
        r"$|S\rangle$",
        r"$|T-\rangle$",
    ],
    col_names=[
        r"$\langle T+|$",
        r"$\langle T0|$",
        r"$\langle S|$",
        r"$\langle T-|$",
    ],
)
plt.show()
HTML(anim.to_jshtml())

Saving animation to rdm-radical.gif...
Animation saved successfully!

Once Loop Reflect

Compare with RadicalPy simulation

if len(basis) < 6:
    assert isinstance(D, np.ndarray)
    H = sim.total_hamiltonian(B0=B0, D=D, J=J)
    time = np.arange(0, 5.0e-8, 5e-10)
    kinetics = [
        rp.kinetics.Haberkorn(kS, rp.simulation.State.SINGLET),
        rp.kinetics.Haberkorn(kT, rp.simulation.State.TRIPLET),
    ]
    sim.apply_liouville_hamiltonian_modifiers(H, kinetics)
    rhos = sim.time_evolution(State.SINGLET, time, H)

    time_evol_s = sim.product_probability(State.SINGLET, rhos)
    time_evol_tp = sim.product_probability(State.TRIPLET_PLUS, rhos)
    time_evol_tz = sim.product_probability(State.TRIPLET_ZERO, rhos)
    time_evol_tm = sim.product_probability(State.TRIPLET_MINUS, rhos)
    x = time * 1e6

    plt.plot(x, time_evol_tp, linewidth=2, label="T+")
    plt.plot(x, time_evol_tz, linewidth=2, label="T0")
    plt.plot(x, time_evol_s, linewidth=2, label="S")
    plt.plot(x, time_evol_tm, linewidth=2, label="T-")
    plt.plot(
        x,
        time_evol_tp + time_evol_tz + time_evol_s + time_evol_tm,
        linewidth=2,
        label="trace",
    )
    plt.legend()
    plt.title(f"RadicalPy: Density Matrix Approach {B0=}, {J=}")
    plt.xlabel(r"Time ($\mu s$)")
    plt.ylabel("Probability")
    # plt.ylim([0, 1])
    plt.grid()
    plt.show()

    rhos = rhos.reshape(
        rhos.shape[0], isqrt(rhos.shape[1]), isqrt(rhos.shape[1])
    )

if len(basis) < 4: fig, anim = get_anim( rhos[::2, :, :], x[::2], title="Full density matrix", time_unit="μs", dpi=30, ) plt.show() display(HTML(anim.to_jshtml())) elif len(basis) < 6: red_rhos = np.einsum("tikjk->tij", rhos.reshape(rhos.shape[0], 4, rhos.shape[1] // 4, 4, rhos.shape[2] // 4)) fig, anim = get_anim( red_rhos[::20, :, :], x[::20], title="Reduced density matrix", time_unit="μs", dpi=30, ) plt.show() display(HTML(anim.to_jshtml()))

# Clean files
!rm -rf radicalpair_*_prop
!rm -f wf_radicalpair_*.pkl