Dipolar relaxation mechanism of long-lived states of methyl groups.

We analyze the symmetry properties of the dipolar Hamiltonian as the main relaxation mechanism responsible for the observed NMR spectra of long-lived states of methyl groups. Long-lived states exhibit relaxation times that are considerably longer than the spin-lattice relaxation time, T 1 . The analysis is complementary to previous studies and provides insight into the relaxation mechanism of long-lived states by focusing exclusively on the symmetry of the spin Hamiltonian. Our study shows that the dipole-dipole coupling between protons of a methyl group and between the protons and an external spin are both symmetry breaking interactions that can lead to relaxation pathways that transform the polarization from symmetry order to Zeeman order. The net contribution of the internal dipolar interaction to the NMR observation of long-lived states is zero. Our calculation is in good agreement with the reported features of the observed spectra and previous theoretical studies.

Introduction

Implementing states with long relaxation times is an attractive approach to quantum memories, sensors in quantum meteorology and as a tool for studies of slow processes in nuclear magnetic resonance (NMR) [1-4].

NMR is normally concerned with states polarized along the external magnetic field at some temperature. Thus, the spin–lattice interaction is the dominant source of energy relaxation with its characteristic time, . There are, however, alternative ways of polarizing the spin system such that the polarized state is protected from the majority of relaxation mechanisms, allowing the relaxation time to substantially exceed . For instance, a group of indistinguishable spins have symmetry properties that are preserved under collective noise processes. These symmetry degrees of freedom can be used for implementing long-lasting states [5-7].

Long-lived states (LLSs) have been implemented experimentally [6-15] and studied theoretically [16-19]. For the simplest case of two spin particles, the four spin eigenstates can be categorized into two symmetry groups: triplet states and singlet state. The former is invariant under spin exchange and the later is antisymmetric. By creating an imbalance of population between the singlet state and the identity in the triplet subspace, one can prepare a symmetry-polarized state that exhibits long relaxation time due to its resistance to the majority of relaxation mechanisms [8, 9]. One can extend this idea to three identical spins, where the eight eigenstates are categorized into three groups labeled A, and where the last two are degenerate subspaces. Under cyclic permutation of the spins, the four eigenstates in group A remain invariant, and so, totally symmetric states, whereas the two eigenstates in each of degenerate subspaces, , acquire a phase , and so, non-symmetric states [12, 17, 20, 21]. In a similar manner to the two-spin case, one can polarize three identical spins with respect to their symmetry by creating an imbalance of population between the A states and the states. This has been implemented in methyl groups, and the relaxation time has been extended by factor of 7 beyond the single spin [12].

The LLS is immune to collective noise because the symmetry of the states is also collective. Although the symmetry is preserved, the magnetization is not. Consider a state of three identical spins that have both symmetry order and magnetization order.1 In the absence of any symmetry breaking interaction and under a collective depolarizing noise, that system can evolve to a state with no magnetization order. Therefore, no signal of that state can be observed in NMR. However, one can still observe a signature of the symmetry order. In fact, any interaction that does not preserve the permutational symmetry of spins, in principle, can lead to a relaxation process that turns the polarization from symmetry order into Zeeman order, converting an initially prepared LLS to an NMR observable state. For example, in [12], accessing the noise-protected subspaces has been achieved at low temperature where the majority of those interactions that distinguishes spins are frozen or have negligible effect. On the contrary, the NMR observations have been done at room temperature where those spin distinguishing relaxation processes are no longer negligible. A few examples of such interactions are the chemical shift anisotropy (CSA) and the dipole–dipole coupling between spins. The authors of [17] have studied the contribution of both of these interactions on the NMR observation of LLS in spherical top molecules. They derive the NMR spectra in the limit of relatively short correlation time of methyl groups. Here, we describe an approach to calculate the contribution of the dipole–dipole interaction to the NMR observation of long-lived states, which focuses exclusively on the symmetry properties of the dipolar Hamiltonian. Our study is complementary to previous studies [12, 17].

The spectrum of initially prepared LLS in methyl groups has been reported and has some unique features [12, 21] that differentiate an initial symmetry-polarized state from a thermally polarized state. The scalar coupling between the methyl group with total spin operator and the carbon with spin operator is given by , which leads to splitting the proton and the carbon spectra into multiple peaks. On the proton channel, one would observe two resolved peaks, and , each corresponding to the carbon’s state being aligned or anti-aligned with the field. On the carbon channel, one would observe four resolved peaks, and , each corresponding to the total magnetization of the methyl group being or . The explicit definition of is given in Eq. 33. In case of symmetry-polarized states, the two peaks on the proton channel have equal amplitudes but with opposite phases. The four peaks on the carbon channel have different amplitudes, and they are two-by-two anti-phase with one another [12, 22], as opposed to a thermally polarized state where all peaks on both spectra are in-phase with each other. In this report, we analytically calculate the NMR spectra of an initially prepared LLS and predict the above- described features that have been observed. Our result is consistent with the results in [12, 17]. Our approach relies on the symmetry properties of the dipolar Hamiltonian and explains the underlying physics that leads to extended relaxation time and observed features of the NMR spectra.

In Sect. 2, we review the definition of long-lived states in the case of three indistinguishable spins. In Sect. 3, we expand both the heteronuclear and homonuclear dipolar Hamiltonians in terms of symmetrized spin operators and analyze their properties. The symmetry analysis allows us to solve the master equation analytically in Sect. 4 and to predict the NMR spectrum of LLS in methyl groups. We report a summary of this work in Sect. 5.

Long-lived state in methyl groups

In the presence of a uniform static magnetic field, , the spin Hamiltonian of protons in a methyl group iswhere is the proton frequency, is the scalar coupling constant between any two protons and and account for the chemical shift anisotropy (CSA) and the dipole–dipole (DD) interaction between any two spins. At relatively large field, when the Zeeman interaction is the dominant term in the Hamiltonian, it is a good approximation to treat these protons as three indistinguishable spins. The spin eigenstates arewith . The other four eigenstates are obtained by replacing with . We denote the spin eigenstates with where the first label is the “symmetry label” and the second label is the collective magnetization. The symmetry label reflects that if we cyclically permute the indistinguishable spins, the eigenstates are invariant and the () states acquire a phase (or ). In addition, all are eigenstates of the Zeeman Hamiltonian. So, they are labeled by the eigenvalues of the z component of the total spin angular momentum operator, . The above states are the eigenstates of the cyclic permutation operator denoted by .

As a first approximation, when the chemical shift anisotropy and dipole–dipole interactions are negligible compared to the Zeeman interaction, the Hilbert space of these identical spins can be partitioned as a direct sum of two subspaces, . Each subspace can be further decomposed into a product of the symmetry label s and the magnetization label m, i.e., and . The first subsystem (s) is not affected by collective noise operators, becausewhere and are complex numbers. All components of the total spin angular momentum, with , preserve the symmetry label although they may corrupt the magnetization label. One can take advantage of this symmetry property by polarizing the methyl group in terms of the symmetry order rather than the conventional Zeeman order. This way, the system is protected against collective spin noise and may exhibit long relaxation times. This has been experimentally demonstrated as long-lived states in methyl groups [12] where an imbalance of population is created between the A subspace and the subspaces. In the following section, we provide a mathematical description for these noise- protected states.

We denote a state that is only populated in a particular symmetry subspace by with . These states are totally polarized in terms of the symmetry label, but they are totally mixed in terms of the magnetization label. Explicitly, we define these density states byWith the above definitions, we introduce a polarized long- lived state asThe appearance of the scalar terms, , assures that is not affected by any component of the total spin magnetization, , with . One could extend this idea and define a more general type of protected states aswhere the A subspace is -polarized relative to the E subspace and the subspace is -polarized relative to the subspace. In this report, we are not concerned about the initialization of the system in or . We assume that a long-lived state is given, and we study the dipolar relaxation mechanism that turns this noise-protected state into an NMR observable.

Symmetry of dipolar coupling

So far, we have neglected that not all interaction terms in the Hamiltonian preserve the cyclic permutation symmetry. Ideally, if the first two terms of Eq. 1 are the only terms in the Hamiltonian, spins are completely indistinguishable, and therefore, the long-lived states can never be converted into an NMR observable. If they are protected against the collective noise, they are protected against the measurement as well. Quantitatively, for all . The chemical shift anisotropy, the internal dipolar coupling among protons in a methyl group and the external dipolar coupling between the group of protons and an external spin (such as C) are all possible interactions that do not preserve the spin symmetry. These physical interactions that distinguish spins can potentially transform the symmetry polarization into Zeeman polarization, leading to NMR observation of what were once protected states. The focus of our study is the contribution of dipole–dipole interaction in the spin read out of an initially prepared protected state. The CSA is normally small for protons in methyl groups in comparison with the dipolar coupling and hence is not considered in this study, although the symmetry arguments would be similar.

Dipolar Hamiltonian

Two magnetic dipole moments and , which are at a distance apart, interact through space via the dipolar Hamiltonian [23],Here, and with and for the case of spin half particles where are the Pauli operators. The use of different symbols S and I emphasizes that these operators act on different Hilbert spaces and may represent two different spin species. The dipolar Hamiltonian is commonly written in a product form of spatial functions and irreducible rank-2 tensors [23],where is a constant, a function of space parameters and a normalized bilinear spin operator. Explicitly,Given these definitions, we haveIn case of the homonuclear interaction, since , each term changes the z component of the total spin magnetization from m to . So, we may refer to q as the order number.

In case of heteronuclear interactions, each is further decomposed to based on its commutation with the Zeeman Hamiltonian, i.e., . The explicit form of these bilinear spin operators and their corresponding frequencies are listed in Table 1. One can write them in a closed form by noting and , where or represents , and or represents .

Table 1

Components of the dipolar Hamiltonian

q	T^_(q,p)			ω_(q,p)
	p=1	p=0	p=-1	p=1	p=0	p=-1
0	-16S_+I_-	83S_zI_z	-16S_-I_+	ω_s-ω_I	0	-ω_s+ω_I
+ 1	S_+I_z	S_zI_+	–	ω_s	ω_I	–
+ 2	S_+I_+	–	–	ω_s+ω_I	–	–

Given the above definitions, the homonuclear coupling between protons in a methyl group is where represents the dipole–dipole coupling between the ith proton and the jth proton and is given in Eq. 7 with and . In addition, is the relative distance between the two protons and is a constant.

Similarly, the heteronuclear coupling between the collective spins, , and the test spin, I, consists of three terms, where represents the dipole–dipole coupling between the jth proton, and the test spin, I, and is the relative distance between them in the Zeeman frame.

The heteronuclear dipolar Hamiltonian does not commute with the spin cyclic permutation operator, i.e., . Similarly, the homonuclear dipolar Hamiltonian does not commute either. This non-commutativity implies that the dipolar Hamiltonian and the cyclic permutation operator, , do not have a common basis, and hence, the dipolar interaction can induce transitions between the eigenstates with different symmetry labels. Consequently, a dipolar relaxation might lead to a process that takes polarization from symmetry order to Zeeman order, which is measurable.

Symmetrized Operators

We introduce the symmetrized bilinear spin operator and rewrite the dipolar Hamiltonian in terms of them. These symmetrized operators provide insight into the key components of the dipolar coupling that leads to observing a protected state. We start with the external dipolar interaction and then elaborate on the internal dipolar coupling and discuss their differences and similarities.

For heteronuclear coupling, we rewrite the Hamiltonian in Eq. 12 by substituting [24, 25] and obtain,where . We emphasize that in the last line, the sum over different spin indices is replaced by a sum over the symmetry label , corresponding to symmetries. The symmetrized complex functions and the symmetrized bilinear spin operator are given byIn closed form, , where the bold notation reminds us that acts on all three spins.

To clarify the effect of on the collective spin system, consider a more simple example of two spins, where and (or ). Thus, is either symmetric or antisymmetric. By looking at the nonzero matrix elements of in the triplet–singlet basis, we conclude that the symmetric tensors cause transitions only within each symmetry subspace and the antisymmetric tensors cause transitions between two subspaces of the triplet states (A) and the singlet state (E). Indeed, for two spins we obtainwhere is the triplet subspace and is the singlet subspace.

Extending this to three spins is straightforward. The of three spins is a totally symmetric operator, and it has nonzero matrix elements only within each symmetry subspaces A, or . Consequently, it is block diagonal in the eigenbasis of the cyclic permutation operator. In a similar manner to the two-spin case, connects two different symmetry subspaces. Moreover, in case of three spins, is neither symmetric nor antisymmetric, and one needs to be careful about the direction of transitions. We expand the in the basis and conclude that transforms the symmetry of the eigenstates in the right cycle as and the transforms the symmetry of the eigenstates in left cycle as . The upper index of determines whether the transformation is within each symmetry subspace or between them in a right/left cyclic direction (. The lower index determines the change in the magnetization. Explicitly, when , takes and when the takes for all values of . In other words,where the sum in is mod 3. The final remark is that the reverse transition process occurs via , whereFor example, the transition from to occurs through but the reverse process, from to occurs through . The allowed transitions due to are demonstrated in Fig. 1.

Fig. 1

Allowed transitions due to : the nonzero transitions between the eigenstates of are shown that are due to the nonzero matrix elements of the symmetrized collective spin operators . The lower index acts on the magnetization label where it takes m to and the upper index acts on the symmetry label in a cyclic manner. The blue/red arrows indicate the nonzero transitions between different symmetry spaces in the right/left cyclic order, and the green arrows refer to nonzero transition within each symmetry space

The homonuclear dipolar Hamiltonian can also be written as a sum over symmetrized operators. We substitute in Eq. 10 and rewrite the Hamiltonian as,whereWe observe that the symmetrized spin operators for internal dipolar coupling, , exhibit similar symmetry properties to that of the external dipolar coupling, . Similarly, the fully symmetrized operators, or A, have nonzero matrix elements only within each symmetry subspace and not between them, and hence, they are not of interest. Thus, the important terms are , and Fig. 2 shows all the nonzero transitions induced by for . Note that in the figure the terms are not considered because they contribute to shifting the energy levels not to the transitions.

Fig. 2

Allowed transitions due to : the nonzero transitions between the eigenstates of are shown that are due to the nonzero matrix elements of the symmetrized bilinear spin operators . From left to right, the figure shows allowed transitions due to the SQ (), the ZQ () and the DQ () spin operators. The blue/red arrows indicate the nonzero transitions between different symmetry spaces in the right/left cyclic. The green arrows indicate those transitions that do not effectively contribute in the NMR signal

The above transitions can be categorized in three groups: the zero-quantum transitions (ZQ) also known as flip–flop terms when and , the single-quantum transitions (SQ) or single spin flip terms when and , and the double-quantum transitions (DQ) or the flop–flop terms when and . Similar to the operators in case of heteronuclear Hamiltonian, the internal flip–flop interaction among protons also leads to transitions between states with different symmetry in a cyclic order (left /right for ). More precisely,These transitions are passive in the sense that even though they lead to an exchange of population between different symmetry subspaces, they will not effectively contribute in transferring polarization from the symmetry order to Zeeman order, as it is explained in the next section.

The and the transitions, the and the operators, are slightly different. They do induce transitions between the A and the subspaces but do not allow transition between the and the subspaces as opposed to the external dipolar coupling. For we haveNow that we visualized the effect of these symmetrized irreducible tensors on the symmetrized eigenbasis, and we proceed to solve the master equation in the next section.

Master equation

At high temperature and in liquid phase, the space coordinates randomly fluctuate in time, as does the dipolar Hamiltonian. In the absence of CSA, the total Hamiltonian consists of a Zeeman term , a scalar coupling , and a fluctuating dipolar term and . The scalar coupling is totally symmetric and does not contribute to breaking the symmetry of the protected states, so, it is neglected in the following discussion, unless otherwise stated. In the following discussion, we analyze one of the two types of dipolar relaxation process (homo vs. hetero) in the absence of the other. Of course, in experiment, these two processes are not independent of each other and the evolution of one affects the other, because their corresponding Hamiltonians do not commute. For a short time period, , one may consider them as two independent processes. We first analyze the pure contribution of heteronuclear coupling in the NMR observation of initially prepared protected states and then analyze the contribution of the homonuclear coupling and discuss their differences and similarities.

Relaxation induced by heteronuclear coupling

At high field, we treat as a perturbed term and apply Redfield’s semiclassical theory [23, 26] to study the collective spin dynamics. Then, we solve the master equation analytically and predict the NMR signal.

The Lindbladian form of the semiclassical master equation gives us [26]where is the density matrix in the rotating frame of the Zeeman interaction. The coefficient is the real part of the symmetrized spectral density of noise and is defined as:Here, denotes the symmetrized autocorrelation function, and the overbar notation refers to averaging over the random variables. The imaginary part of the spectral density of noise leads into the dynamical shift and can be absorbed in the coherence evolution part [27].

For the purpose of the following discussion, the explicit form of the is not required, because we are not interested in the exact dynamics of the system. Rather we would like to know which components of the heteronuclear dipolar coupling convert the symmetry polarization of the into the measurable Zeeman polarization.

We define the Lindbladian map with in which L is the Lindblad operator and rewrite the master equation aswhere , and for all cases. Considering the classical noise and neglecting the term which just shifts the energy, one can break the above master equation into three parts: the zero-quantum transitions (ZQ), the double-quantum transitions (DQ) and the single-quantum transitions (SQ),Here, we replace . The ZQ and the DQ terms exchange energy between the collective spin and the test spin, and the SQ terms change either the collective spin states or the test spin states. As was mentioned before, the totally symmetric Lindblad operators do not cause a transition between the two different symmetry subspaces of the collective spin. Therefore, if it happens that the spectral density of noise is very well approximated with just the totally symmetric component, i.e., , one can conclude that the system is very robust against noise and exhibits very long relaxation time. This is in agreement with the result in [12, 17] where in the limit of very fast rotations of the methyl group, the LLS becomes the eigenstate of the relaxation superoperator, preserving the population imbalance. For the sake of simplicity in the following discussion, we ignore all of the totally symmetric Lindblad operators, since they do not play a critical role in observing the protected state. We also neglect the terms, because we are interested in those transitions that the collective spin exchanges the energy with the test spin. The only important terms in the dissipator are and .

To obtain the allowed transitions due to ZQ and DQ terms, we need to calculate the effect of on the eigenbasis . The nonzero components are,Similarly,The proportionality constant is 1 for transitions from or to the and is for all other levels.

Based on the above relations, the allowed transition due to and terms is indicated in Fig. 3 in which transitions are color coded with red and blue, respectively.

Fig. 3

Left segment: the selection rule due to transitions. Right segment: the selection rule due to transitions. Red versus blue refers to nonzero transitions due to versus

The calculation in Eq. 25 convinces us that if the initial state is a probabilistic mixture of different energy levels, the above master equation reduces to the classical rate equations to:where [x](t) is the population of the xth energy level at time t with , and is the transition rate between two energy levels x and y. To solve the above differential equation, the evolution time is discretized into N steps where and is small compared to the energy scales of the system. Given the population distribution at time , the change in the population of xth level at time is approximated by . Therefore, if the initial conditions are known, one should be able to calculate the NMR signal by solving the rate equations recursively. Note that at any instance of time, the transition occurs only between those energy levels that first are allowed due to and terms, and second, there is an imbalance of population at .

The initial state of interest iswhere is the polarization of the test spin and is replaced from Eq. 6. Since is diagonal, we alternatively represent it by a population vector whereThe represents the vector of population of three identical spins at in which all energy levels with the same symmetry label are equally populated, but there is an imbalance of population between energy levels with different symmetry labels. For use in future analysis, we computewhere is the imbalance of population between the states and states and is the imbalance of population between states and states for all possible values of m and and at .

Considering an unpolarized test spin, , the solutions to the rate equations for short evolution time, , arewhere we denoted with and with for short. Since all energy levels in A subspace are equally populated at , in the absence of and terms, we obtainThe expressions in Eq. 31 may appear very complicated and may sound difficult to get an insight into the relaxation. If we pay attention to the symmetry, there is a delicate and simple relation between the population of different energy levels. For all values of (s, m), the change in the population of level is the same as that of with just the difference of replacing . Second, the change of population in each level is the same as that of with just the difference of replacing . It will be shown that the anti-phase feature of the NMR peaks observed by others arises from these two properties, which are also visually captured from Fig. 3.

For the subspace, after summing over the symmetry labels and doing some algebra, we obtainwhere . We see in the following that these very neat relations between different energy levels enable us to anticipate the NMR signal analytically.

To be more specific, we start from the NMR signal on the test spin channel. The scalar coupling that was neglected so far shifts the frequency of the test spin condition on the total spin magnetization of protons. Therefore, it is expected to observe four distinguishable peaks on the I channel corresponding to and . We define a set of operators that project the collective spins into these magnetization subspaces with:For a short evolution time, the expected NMR peaks at the test spin channel arein which . To compute the above NMR signals, we used and replaced the expressions from Eq. 31 into it. On the carbon channel, both and parts of the density matrix become an NMR observable, leading to an anti-phase term that is proportional to and an in-phase term that is proportional to . Similarly, the anticipated NMR peaks on the proton channel arein which . On the proton channel, becomes an NMR observable, leading to two equal peaks with opposite phases.

It remains to solve the rate equations in case of . When the test spin has some initial polarization, , in Eq. 30 additional terms show up. Precisely, the imbalance of population difference between different energy levels is now replaced by and . By replacing these initial imbalance of population in the rate equation (Eq. 26) and after doing some tedious calculations, we obtainInterestingly, all terms that have dependency are proportional to . This means that the polarization difference between the and the subspaces become observable if . This also means that if one is able to create an imbalance of population between the and subspaces, that wavefunction leads to an extra feature in the observed NMR spectra which reveals information about the symmetry of the spectral density of noise.

In case of , Eq. 35 reduces toThe anti-phase contribution in the NMR spectra is proportional to , the initial imbalance of population between the A and E symmetry subspaces. The in-phase contribution is proportional to , which is the polarization of the test spin that is dipolarly coupled with the three protons. We do not know the explicit value of different orders of the spectral density of noise and therefore cannot make any conclusion about the relative magnitude of these two contributions. Nevertheless, we can conclude that the peaks in pair of or have the same phase and the same amplitude. Now, if , even for the case of , we can no longer make any judgment about the relative amplitude and phase of these four peaks. Because, according to Eq. 35, even when an additional anti-phase term survives for the peaks. In “Appendix,” we compute the NMR peaks for one step further when and conclude that for a longer evolution time, additional terms show up that differentiates the amplitude of the from that of . Thus, four peaks with unequal amplitudes are expected. This argument holds true without any assumptions about the magnitude of the anti-phase contribution relative to the in-phase contribution. Now, we assume those terms with a factor are larger than the others and conclude that the m peaks are overall anti-phase with the peaks but with unequal amplitudes. A large corresponds to the long-lived state that was experimentally demonstrated in [12], and surprisingly the above model, which considers only the and the transitions of the DD coupling between the collective spin and an external spin, results in an analytic solution that predicts most features of the NMR peaks that have been experimentally observed. Here, we are not interested in the exact dynamics of the spins system and focus on the permutation symmetry properties of spin operators of the Hamiltonian to obtain insightful understanding of the symmetry breaking mechanism. This is complimentary to the previous study [17] where the authors calculate the correlation function of the spherical top molecules explicitly and solve the master equation with the zeroth-order approximation numerically. In their approach, the nonzero matrix elements of the relaxation operator in the fully symmetric subspace are computed where in the limit of very fast rotation, the symmetry breaking components vanish. Our analytic study avoids the complexity of computing the explicit form of the spectral density of noise and avoids calculating large relaxation matrices and instead finds certain symmetry breaking patterns that lead to converting LLS into an NMR observable state. In both approaches, in a limit where the spectral density of noise is fully symmetric, LLS exhibits extremely long relaxation time.

Relaxation induced by homonuclear coupling

In case of internal dipolar coupling, the analogy to Eq. 24 iswhereis the spectral density of noise due to the homonuclear dipolar Hamiltonian. We simply denote by .

We solve the master equation using the same approach we used in the previous section for heteronuclear coupling. In order to calculate the change of population of each energy level during the short evolution time , we need to know which transitions are allowed (the nonzero matrix elements of as shown in Fig. 2) and the initial imbalance of population between the allowed transitions. Thus, for each energy level , we should calculate in which are the transition rates and is the population of the xth level at .

For the initial state of interest (Eq. 27), at , all the energy levels have the same population as the energy levels. On the other hand, the transition rates induced by terms of the internal dipolar coupling are identical for both states and . As a result, the change of population of is equal to that of f , leading to a zero polarization in the Zeeman basis. Therefore, even though there is an exchange of population between different symmetry subspaces, as shown in the middle section of Fig. 2, the net contribution of flip–flop interaction among protons to the NMR spectra of a protected state is zero.

Considering the above point, the solution of the master equation in Eq. 36 for a short evolution time, , isGiven the above relations and for all , the expected NMR spectrum of protons for short evolution time isWe repeat the above calculation for the second time step, , using the method explained in “Appendix” and concluded that the signal is proportional to . This means that if the and the subspaces are equally populated at , the higher-order contribution of the homonuclear dipolar relaxation to the NMR spectra is also zero. We justify this counterintuitive result by comparing the spin operators of the heteronculear coupling to that of the homonuclear couplings. The and the transitions of the external dipolar interaction occur when the protons of methyl group exchange energy with an external spin. This effectively appears as a single spin flip (or a operator) in the space of protons. In that sense the of the homonuclear coupling is similar to the of the heteronuclear coupling. There are two important points that differentiate the two relaxation processes. The first point is that the is block diagonal in the E subspace, meaning that there is no exchange of population between the and the subspaces as opposed to the operator. As a result, when these two spaces are equally populated at , after an evolution time , the net polarization in the Zeeman basis remains zero, because both of these subspaces exchange population with the A subspace with the same rate but in opposite direction. Thus, the homonuclear coupling passively breaks the symmetry of an initially prepared protected state, meaning that even though it induces nonzero transition between different symmetry subspaces, its effective contribution to the NMR observation of LLS is zero. The second point is that, in case of heteronuclear interaction, when a level experiences a transition with the rate for instance, its counterpart experiences a transition with the rate . The differences between the two rates, , appear as a nonzero factor in the NMR signal (Eq. ). This holds true even if the and the subspaces have identical initial population. Thus, the heterouclear coupling actively breaks the symmetry of an initially prepared protected state because it leads to a nonzero signal.

Conclusion

We analyzed the contribution of both the internal and the external dipolar interactions of methyl groups in NMR observation of long- lived states that are initially prepared as symmetry-polarized states. The symmetry properties of spin operators and the density matrix allowed us to first show why these states exhibit long relaxation times and second to solve the master equation analytically and obtain the NMR spectrum. Our study provides insight into the nature of symmetry breaking interactions that lead to converting LLS into an NMR observable state and is in agreement with the reported experimental observations [12]. As complimentary to the previous study [17] where the authors solved the master equation by considering a particular correlation function, one can use the symmetry breaking patterns that we derived here to simplify the calculation of the master equation for any arbitrary spectral density of noise. We find that even though both the heteronuclear and homonuclear dipolar coupling break the symmetry of LLS, only the former leads to a relaxation pathway with nonzero polarization transfer from symmetry order to Zeeman order. This approach might find application in other areas in NMR where the underlying physics of the relaxation mechanisms is of interest.