The topic was discussed in a paper “Dynamic nuclear polarization by photoexcited-triplet electron spins in polycrystalline samples” by K. Takeda, K. Takegoshi, and T. Terao. I am going to give some supplementary on how the calculation can be done.
In a spin-1 system, the spin operator and the rotation operator are easy to calculate.
The zero-field splitting Hamiltonian in the ordinary basis of is
The eigen vectors in the paper are
with eigen values , , and . The convention that the has an imaginary number is not necessary. In the paper, there are typos on the eigen values.
In other way, the zero-field splitting Hamiltonian can be diagonalized into
where is the transformation changing from ordinary basis to eigen basis.
Given that the population of the eigen basis is , the density matrix is
The corresponding population in ordinary basis is the diagonal elements of
In the paper, it gives . In order to get the population when the molecule X-axis is aligned with the Z-axis Lab frame, the has to be rotated by 90 degree, and we have the population .
The Zeeman Hamiltonian is , where . In the paper, the author fix the crystal and rotating the external magnetic field. To express the into eigen basis. We simply apply a transformation as the density matrix. But first we need to rotate the magnetic field. The rotation is simple, as the direction of the magnetic field is . Thus, the Zeeman Hamiltonian is .
The transformation is
again, the total Hamiltonian is difference from the author.
The difference is because the normal vector is used in the paper.
After that, the calculation for the power sample is the same.
Instead of using integration, we can generated the power using isotropic distribution. I generated 10k molecule, the result of 9.65 GHz microwave is
where the y axis is the magnetic field in Tesla. The result is similar to the paper. I think that their calculation somehow changed the eigen value, so that the eigen vectors were matching with the eigen energy. And the notation of the normal vector and the constant factor of the eigen vectors does not matter.