1. Introduction
Investigation of spatial parity (
) violation effects at low energies is an effective tool for developing models of fundamental interactions in physics. Hypothesizing of the neutral weak current led to development of the standard model (SM) [
1,
2,
3]. The claim of the SM stimulated to search for the
-odd effects in different atomic and molecular systems. However, to extract the SM parameters from the atomic (or molecular) experiments, one needs to perform considerably sophisticated theoretical calculations, since all electrons contribute to such effects. There are two main types of experiments to search for the spatial parity nonconserving (PNC) interactions in atomic (or molecular) systems. The first type is the optical dichroism observation, i.e., the observation of asymmetry in the number of emitted (absorbed) left- or right-handed photons. Such experiments were proposed by M. Bouchiat and C. Bouchiat on Cs [
4] and performed by the Wieman group later [
5]. This is the most accurate experiment in a low-energy regime that supports the SM theory and strongly restricts any “new physics” beyond it (see, e.g., the review [
6]). The second type of experiments is the optical rotation of the light polarization plane in atomic (molecular) vapors. It was originally proposed by Khriplovich on Bi [
7] and performed by Barkov and Zolotorev [
8].
Systems with a simpler electronic structure are of interest due to a more direct theoretical interpretation. In neutral atoms, the
-odd electron–nuclear (
e–N) interaction effect increases approximately as
with the charge of the nucleus
Z [
4,
9]. Therefore, a reasonable choice could be few-electron heavy highly charged ions (HCI). They have large nuclear charges and it is possible to perform very accurate calculations of their electronic structure. In [
10], it was proposed for the first time to observe the PNC effects in He-like HCI. Such proposal was stimulated by the crossing of opposite parity levels with certain
Z values. Despite a number of suggestions to search for the PNC effects in HCI, no experimental results have been reported. Recently, a new PNC experiment with HCI in storage rings has been proposed in [
11] where it was suggested to employ a new method of preparation of polarized HCI nuclei—by the polarized electrons capture.
In most of atomic experiments to study PNC effects, the
-odd electron–electron (
e–e) interaction is hidden by the
-odd
e–N one. The latter usually dominates [
4,
9,
12,
13]. The
-odd
e–e interaction was observed exclusively at high energies by the SLAC collaboration [
14]. Up to now, there were no other evidence for the presence of the PNC
e–e interaction. In [
15], it was proposed to measure PNC optical rotation due to the weak
e–e interaction on the oxygen molecule which assumed to be the dominant source of the parity violation due to the formation of a chemical bond. Recent calculation of the PNC effects in the O
molecule [
16] showed that the enhancement of the PNC
e–e exists in this case, compared to the atomic one, but still, the effect is suppressed compared to the dominant PNC
e–N interaction. In [
17], it was shown that the favorable situation for observation of the PNC
e–e interaction effect appears in the para-H
molecule.
In an atomic experiment, [
5], on the nuclear spin-dependent (NSD) parity violating the
e–N interaction effect observation, the value of the nuclear anapole moment was extracted. The effect of the electromagnetic interaction of electrons with the anapole moment of the nucleus prevails over other sources of the PNC-NSD
e–N interaction effects in heavy atomic systems. In that way, nuclear spin-dependent spatial parity nonconserving (NSD-PNC) effects due to the nuclear anapole moment were also observed by the Wieman group [
5]. However, the hydrogen molecule is the unique molecular system in which the dominant source of the NSD parity violation is due to the neutral weak current interaction. Therefore, it is the favorable system to extract the coupling constant for this interaction, the Weinberg angle, in order to test the SM and put constraints on its extensions.
The aim of this paper is to complete the investigation of PNC effects in the H
molecule and to perform ab-initio calculations of the nuclear spin-dependent part of the PNC effects in orthohydrogen. In the ortho-H
molecule, the total nuclear spin is equal to one. The magnetic dipole (M1) transition between different vibrational sublevels with the same rotational numbers of the ground electronic state X
of the orthohydrogen molecule is considered. In [
18] it was demonstrated that such transitions in the H
molecule are of the M1 type. The M1 amplitude is nonzero due to the nonadiabatic corrections. Note that for the observation of the parity violating effects it is favorable to use the M1 transitions [
9]. Here we also demonstrate that using the cavity-enhanced scheme (reported in [
19]) and record achievement of the intracavity absorption spectroscopy (ICAS) (described in [
20,
21]) the experiment on the observation of the PNC-NSD effect due to the neutral weak current in the orthohydrogen molecule looks feasible enough.
2. Theory
Here we study the parity violation mechanism in the transition between and vibrational sublevels of the ground electronic state of the orthohydrogen molecule. The initial and the final states have equal rotational numbers , where is the total angular momentum of the molecule less the electron spin.
Owing to the spatial symmetry violation produced by an effective operator of PNC interaction
the gerade (
g) wave function of the H
ground state gets the admixture of the whole spectrum of the ungerade (
u) ones (i.e., of the opposite parity). As a consequence of such mixing, the pseudoscalar term appears in the probability of the considered process:
Here we introduced the factor
—the PNC, or the circular polarization degree. In our case the pseudoscalar term is
, where
is the spin of photon,
is the direction of photon emission,
is the polarization of photon. In a general way
where
,
are the PNC-induced electric dipole and the ordinary magnetic dipole amplitudes between the states of the same parity (
g), respectively.
In Equation (
3) the electric dipole operator
and in Equation (
4) the magnetic dipole operator are introduced. A spherical components of
and
are denoted by the letters
f and
q, respectively. Since the considered interaction is the
-invariant one, the PNC matrix element is imaginary [
9].
Within the SM framework, the effective operator describing NSD parity violating
e–N interaction for the case of the diatomic molecule is given by
where
is the Fermi constant (
is the proton mass),
is the nuclear spin,
is the Dirac matrix,
is the normalized to 1 nuclear density distribution and
R is the internuclear distance. From Equation (
5) it follows that the effect of the NSD-PNC
e–N interaction is proportional to the dimensionless constant
. It includes three major types of contributions (see the review [
6]). The first one arises from electromagnetic interaction of electrons with the anapole moment of the nucleus [
22,
23]. In a simple nuclear shell model [
23] for heavy nuclei this contribution
can be estimated as
where
(
l is the orbital angular momentum of the single unpaired nucleon),
and
are the magnetic moment in nuclear magnetons and the weak coupling constant of the unpaired nucleon,
A is the number of nucleons. The next contribution arises from the electron–nucleus NSD neutral weak-current interaction [
24]. In the nuclear shell model this contribution
reads
where
is the coupling constant for the valence nucleon. For protons and neutrons it can be expressed as follows:
Here
,
is the Weinberg angle. In our paper, we approximate the Weinberg angle as
. The last contribution arises from the effect of the nuclear spin-independent
interaction combined with the magnetic hyperfine interaction [
25] (i.e., appears as the radiative correction to NSI-PNC effect). This contribution
can be estimated as (see, e.g., [
6])
The
constant thus reads
In heavy atomic and molecular systems usually the anapole moment contribution dominates, i.e., from all previous experiments on the search for the NSD-PNC interaction effects the anapole moment of the nucleus was usually extracted. Let us estimate the contributions to
for the case of the orthohydrogen molecule. In our case Equation (
6) is not valid, since the nuclei are pure protons. In [
26] it was shown that each particle with nonzero spin possesses an intrinsic anapole moment due to radiative corrections to the PNC-NSD
e–N interaction and thus can be attributed to the
constant. In [
9] an upper bound of this contribution was roughly estimated as
. In the present case Equation (
9) also yields in upper bound for the value of
. Since we have unpaired protons in the system considered, we use the value of
. For one nucleus
From Equation (
11) it follows that the orthohydrogen molecule is the unique molecular system where the dominant source of parity violating NSD effects is the
e-N neutral weak-current interaction and where the
coupling constant can be directly observed.
The ground electronic state
of the orthohydrogen contains only odd rotational numbers
[
27]. The NSD-PNC operator (Equation (
5)) is an electronic vector, so the matrix of this operator can contain
entries. However, the major contribution to the effect considered is arising from the mixing of
ground state and
states (since the admixed E1
transition is allowed). Besides,
. In homonuclear orthohydrogen molecule for each rotational quantum number
N only the
state survives.
Spatial parity violation interaction leads to the interference between M1 and PNC-induced E1 amplitudes. The derivation of explicit expressions for them is based on the use of complete wave function of the molecule that takes into account also nuclear spin variables. Applying the Wigner–Eckart theorem and using the expression for the matrix element of a single operator in a coupled scheme (see, e.g., [
28]), we express the M1 amplitude via the reduced matrix element:
Here
, where the nuclear magneton
and the rotational
g-factor
g was introduced. Note that the vibrational wave-functions that correspond to different vibrational sublevels of the same electronic level are orthogonal to each other (
. However, in [
18], it was shown that this
g-factor depends on
R due to the nonadiabatic corrections, therefore such M1 amplitudes are different from zero. In what follows the M1 amplitudes from [
18] are employed. According to [
18],
for
(
e is the charge of the electron and
is the Bohr radius), i.e., they weakly depends on
J.
Now let us write E1
transition amplitudes with intermediate admixed
terms. For simplicity we consider only
term in
. Similar derivations show that the contribution of
term is equal to the
one. Performing angular reduction gives the following result:
Performing now the transformation from the laboratory frame of reference to the molecule-fixed frame according to [
27,
28] results in
Here
(since we consider only singlet molecular states) is the constant that arises from the transformation to the molecule-fixed frame. Note that
,
and
[
27].
Denoting
and summing over the
projections in the M1 and E1 interference contribution it follows that
In Equation (
16)
The minus sign in Formula (
16) is due to the imaginary matrix element
. The choice of the hyperfine sublevel depends on the certain experimental conditions. For different
F and
J, the coefficient
S ranges between
and
. Without restricting the generality we set
. The same holds true for the coefficient at the reduced matrix element of
for M1 amplitude in Equation (
12).
4. Results and Discussion
In the present work we took into account the contribution originating from the mixing between the ground and the first excited C
state. Note also that the energy denominators in Equation (
3) for higher excited terms are larger, thus within the claimed uncertainty (about 10%) their total contribution to the parity violating effects can be neglected. The calculated dependencies of the imaginary component of
and the E1 matrix element on
R are given in
Figure 1a,b.
The final result for the nuclear spin-dependent E1
amplitudes in ortho-H
is as follows:
Then, the parity nonconserving degree for the effect considered now reads
Note that according to [
17] the nuclear spin-independent PNC effects for the molecular system considered are negligible compared to the NSD-PNC one.
Finally, we perform an analytical modeling of the PNC optical rotation experiment on the ortho-H
molecule. In what follows we will use the record achievements of the cavity-enhanced polarimetric scheme introduced in [
19] and intracavity absorption spectroscopy (ICAS) discussed in [
20,
21]: the record optical pathlength
km and the record shot-noise-limited birefringence-phase-shift sensitivity
rad. In our derivations, one should take into account the light absorption while propagating through the medium. Obviously, working at the resonance frequency of the considered transition, the laser light will be immediately absorbed. But in [
19], the possibility to shift off-resonance in such ICAS experiment is discussed. For the off-resonance measurements the PNC optical rotation angle reaches its maximum value and the absorption is too small (see, e.g., [
9]). The rotation angle in optical experiments increases with the optical pathlength
l linearly. The light transmission through a medium is governed by the Beer–Lambert law [
34], that is, falls as
(
is the absorption length). The PNC signal (the product of the rotation angle and the light transmission)
. So the optimal signal-to-noise ratio corresponds to
[
9,
35].
The PNC optical rotation angle
can be expressed as [
9]
Here
is the vapor number density and
resonance frequency of the transition,
is the Doppler width. The latter can be expressed as [
9]
where
. Here
M is the molecule mass,
T is the absolute temperature and
is the Boltzmann constant. The Voigt profile can be introduced as:
Here we introduced two dimensionless variables.
u is determined in the following way:
where
is the detuning of the frequency. Variable
v is determined by
where
is the collisional broadening width (in the case of forbidden M1 transition the contribution originating from the natural line width is negligible) and
where
is the collisional cross-section. The
function in Equation (
20) can be expressed via the Voigt profile in the following way:
The quantity inverse to the absorption length reads
where
For
(i.e., for the off-resonance measurements) the
and
functions have the following asymptotics:
Using the condition
[
9,
35] of the optimal signal-to-noise relation, as well as Equations (
20)–(
29), one obtains
Equation (
27) is valid only for
when the electric quadrupole amplitude is less than the magnetic dipole one (see [
18]). Note that for
there is an additional absorption originating from the E2 contribution, thus the PNC effect is suppressed. In what follows we use
. Since
, cooling of molecular vapor does not improve the situation. We use
s
, the characteristic value of the collision cross section for H
cm
,
K,
km [
20]. Then the optimal PNC optical rotation angle can be estimated as
The value of the PNC rotation angle in Equation (
32) corresponds to the number density
It is worth to mention again that using cavity-enhanced scheme the record shot-noise-limited birefringence phase-shift at the level
rad has been demonstrated [
21].