Articles | Volume 2, issue 2
Magn. Reson., 2, 741–749, 2021

Special issue: Geoffrey Bodenhausen Festschrift

Magn. Reson., 2, 741–749, 2021

Review article 08 Oct 2021

Review article | 08 Oct 2021

Mechanisms of coherent re-arrangement for long-lived spin order

Mechanisms of coherent re-arrangement for long-lived spin order
Florin Teleanu1,2 and Paul R. Vasos1,2 Florin Teleanu and Paul R. Vasos
  • 1Extreme Light Infrastructure Nuclear Physics ELI-NP, Laser Gamma Experiments Department (LGED), Horia Hulubei National Institute for Physics and Nuclear Engineering IFIN-HH, 30 Reactorului Street, 077125 Bucharest-Măgurele, Romania
  • 2Interdisciplinary School of Doctoral Studies, University of Bucharest, Blvd. Regina Elisabeta, 030018 Bucharest, Romania

Correspondence: Paul R. Vasos (


Long-lived spin order-based approaches for magnetic resonance rely on the transition between two magnetic environments of different symmetries, one governed by the magnetic field of the spectrometer and the other where this strong magnetic field is inconsequential. Research on the excitation of magnetic-symmetry transitions in nuclear spins is a scientific field that debuted in Southampton in the year 2000. We advanced in this field carrying the baggage of pre-established directions in NMR spectroscopy. We propose to reveal herein the part of discoveries that may have been obscured by our choice to only look at them through the experience of such pre-established directions at the time. The methodological developments that are emphasised herein are the mechanisms of translation between the symmetric and non-symmetric environments with respect to the main magnetic field B0. More specifically, we look again thoroughly at zero-quantum rotations in the starting blocks of long-lived state populations, magnetisation transfers between hyperpolarised heteronuclei, and protons. These pulse sequences seed subsequent magnetic mechanisms that contribute to further applications. For instance, we show how some of the introduced coherence rotations were combined with classical pulse blocks to obtain two-dimensional correlations between protons and heteronuclei. We hope the pulse sequence building blocks discussed herein will open further perspectives for magnetic resonance experiments with long-lived spin order.

1 Introduction

This paper is an opportunity to present several magnetic resonance concepts free of particular application-specific introductions. This may allow such concepts for what they are worth simply as magnetisation transfer mechanisms and comments on their potential usefulness in further experiments. We point out that all concepts presented herein were already addressed, albeit concisely, in references (Sarkar et al., 2007; Vasos et al., 2009; Ahuja et al., 2010; Sarkar et al., 2010) or in the Supplements of these papers.

Proposing that a presentation free of application-specific introductions may reveal magnetic resonance progress to the fullest implies that the drive for traditional discipline-oriented applications may have obscured part of the concepts in the original papers. A legitimate question is whether these articles would have been accepted by the journal editors without the applications in mind or whether they would have been worth accepting. With hindsight, doubts raised by editors whether our work in Lausanne and Paris brought actual progress for applications were far more severe than any doubts regarding the soundness of the work itself. For instance, in the search for new singlet-based excitation sequences (Carravetta et al., 2004) on the route of hyperpolarised magnetisation to long-lived spin states (LLS), we were never tormented by the question “is transport of hyperpolarisation really long-lived?” (Vasos et al., 2009) (Pileio, 2020). However, “is LLS-based polarisation storage in peptides better than the mere longitudinal relaxation time constant of heteronuclei with which peptides are often isotopically enriched, T1(15N), T1(13C)?” was a harrowing question. Equally present were the doubts “are long-lived states, with their complicated excitation and sustaining mechanism, really a better way of measuring slow diffusion, slow exchange constants than heteronuclei (Ferrage et al., 2003) such as 15N?” or “are long-lived coherences (LLCs) actually a good route to improved spectral resolution in NMR?”.

When we dedicated the first of a series of papers (Sarkar et al., 2007) to Anatole Abragam along with a letter expressing our hopes that the discoveries may be useful for diffusion studies, he seized the essence of our work in his answer (mainly addressed to Geoffrey Bodenhausen): “nice to see a way of skillfully sending spins to sleep in their soft bed”, “envoyer les spins se reposer dans leur lit douillet”. The remark, thus made rhythmic by alliteration, was as concise as it was exact, since the singlet state we were searching for is magnetically inactive, i.e. the spins are “sleeping”. This commentary alone may have replaced the introduction to our original paper.

2 Zero-quantum rotation in the starting block of long-lived states

The structure of singlet–triplet population differences or long-lived state operators, QLLS,

(1) Q LLS = | S 0 S 0 | - 1 3 ( | T - 1 T - 1 | + | T 0 T 0 | + | T 1 T 1 | ) ,

was first discussed in formulae adapted to the zero-field magnetic structure for a two-spin system, as first created (Carravetta et al., 2004) in non-equivalent nuclei. While the preference for writing highly symmetrical long-lived states in spherical tensor operators is natural, we recurred in Lausanne, however, to Cartesian operators (Sørensen et al., 1984) in the Liouvillian space:

(2) Q LLS = - N LLS ( I x S x + I y S y + I z S z ) ,

with NLLS=4/3.

The form of this operator allowed us to understand the structure of coherences prone to evolution:

(3) LLS = - 4 3 Z Q x - 2 3 2 I z S z ,

where ZQx=122IxSx+2IySy is a zero-quantum coherence.

Under such a configuration, the system is immune to the scalar-coupling evolution and also to chemical shift evolution, provided the chemical shift difference between the two spins is eclipsed by ample radio-frequency radiation or by cycling the main field (Cavanagh et al., 1995).

Equation (1) proved to be one the most useful formulae in developing the general theory of long-lived states by pointing out the very nature of their extended lifetime, the population imbalance between states, or manifolds of different symmetries with respect to spin permutations (Stevanato, 2015, 2020) which cannot be interconverted by relaxation mechanisms with certain symmetries. The novelty of Eq. (3) was that it strongly connected the singlet-state explorations to research in the Ionel Solomon consecrated (I, S) homonuclear and heteronuclear magnetisation transfer (Solomon, 1955).

The first method of excitation for long-lived populations, developed by Levitt and collaborators (Carravetta and Levitt, 2004), worked for a pair of spins I and S provided carefully chosen delays dependent on the chemical shifts, νI and νS, were used, making a sweep through frequencies necessary to excite different pairs of coupled spins (I, S), (R, K), etc., in different experiments, just like one-dimensional magnetic resonance spectroscopy necessitated a sweep of the main field through resonance conditions for different chemical environments before the introduction of Fourier transform (Ernst and Anderson, 1966; Ernst, 2021 – Nobel Lecture). The chemical-shift dependency of long-lived states rendered impossible any two-dimensional investigations of phenomena involving two or more spin pairs or several chemical environments of spin pairs with encoded LLS, such as exchange or interaction dynamics, in the same experiment.

The first concept introduced in the Lausanne paper (Sarkar et al., 2007) was the broadband excitation of singlet states. The topic may have deserved, in retrospect, a paper on its own. Our way towards broadband LLS excitation passed through zero-quantum coherences, as explained below. The first attempts to excite QLLS in Lausanne (Fig. 1) posed challenges regarding the evolution and relative orientation of zero-quantum coherences (ZQx, ZQy) and ZZ magnetisation (2IzSz).

Figure 1Pulse sequence adapted from a reference (Sarkar et al., 2007) showing the evolution of the density operator at different stages with an emphasis on its three projections (ZQxZQy, 2IzSz) between time points C and E (figure generated with SpinDynamica; Bengs and Levitt, 2018). Here, τ1=1/(4JIS) and τ2=1/(2ΔνIS), where JIS is the scalar-coupling constant (Hz) between spins I and S and ΔνIS=νI-νS is the chemical shift difference (Hz) between the two spins. The phase cycling is ϕ1=x,-x, ϕ2=x, ϕ3=2(y),2(-y) and ϕrec=x,-x,-x,x.

After the first 45 pulse, at point C in Fig. 1, the density operator takes the following expression:

(4) ρ pre - LLS C = 2 I x S x - 2 I z S z = Z Q x + D Q x - 2 I z S z .

At first, we expected to induce the presence of QLLS at point C in this sequence, due to the presence of projections onto QLLS by both longitudinal two-spin order and 2IxSx. However, we realised these two contributions exactly cancel one another out, leaving us at a loss as to how to excite singlets in a broadband manner. We could have anticipated the mutual cancellation by expressing the operator at point C in the singlet–triplet basis, relevant upon application of a “sustaining” radio-frequency (rf) field:

(5) ρ pre - LLS C = - 1 2 ( | T - 1 T - 1 | + | T 1 T 1 | - 2 | T 0 T 0 | - | T - 1 T 1 | - | T 1 T - 1 | ) ,

indicating there was no singlet order to be found there.

By applying the pulsed field gradient g1, the double-quantum term dissipates and the density operator becomes

(6) ρ pre - LLS D = ( I x S x + I y S y ) - 2 I z S z = Z Q x - 2 I z S z ,

which possesses all the components of the long-lived state (Eq. 1) but displays an opposite orientation of zero-quantum and ZZ components with equal projections onto QLLS. Therefore, these projections cancel each other out. To better understand this apparent conundrum, we can write the operator on the basis of singlet–triplet operators:

(7) ρ pre - LLS D = 1 2 | T 0 T 0 | - | S 0 S 0 | - 1 2 | T - 1 T - 1 | + | T 1 T 1 | - | T 0 T 0 | - | S 0 S 0 | = 1 2 2 | T 0 T 0 | - | T - 1 T - 1 | - | T 1 T 1 | .

The next step was figuring out how to interchange the singlet and central triplet populations in order to get the expression for LLS, a task which is not immediately obvious in this form. To do that, a reconversion in the Cartesian product basis proved fruitful: |T0T0|-|S0S0|=IxSx+IySy=ZQx. After several weeks of calculations, a group seminar was dedicated to the otherwise well-known evolution of ZQx under a scalar-coupling interaction (Cavanagh et al., 1995). The new aspect was that the rotation axis was this time also apparent in the density operator expression, so effectively one of the constituents of spin order was rotating around the other, thus changing the relative sign of the constituent product operators ZQx and 2IzSz and yielding the sought-after LLS (Fig. 2). The atmosphere in the magnetic resonance laboratory in Lausanne should be credited for a substantial contribution to the birth of these concepts. However, as a side note, the physical exercises of magnetisation succeeded in captivating more attention on paperback than in the coffee-table setting around a group-meeting white board.

Figure 2Position of coherence prior to and after zero-quantum rotation with an auto-relaxation rate constant corresponding to the density operator given by the two linear combinations of ZQx and IzSz. Only the dipolar relaxation mechanism was considered for a pair of two spins where bDD is the dipolar coupling constant and J(ω) is the spectral density.


The evolution of the ZQx during free precession is

(8) Z Q x ( 2 π ν I I z + 2 π ν S S z ) τ Z Q x cos 2 π Δ ν I S τ + Z Q y sin 2 π Δ ν I S τ ,

where νI is the Larmor frequency of spin I, νS is the Larmor frequency of spin S, ΔνIS=νI-νS is the chemical shift difference (Hz) between the two spins, and JIS is the scalar-coupling constant (Hz) between spins I and S. Thus, after an evolution period τ2=1/(2Δυ), the ZQx will change sign such that 12(|T0T0|-|S0S0|)(2πνIIz+2πνSSz)τ12(|S0S0|-|T0T0|). During the evolution represented in the coherence clock in Fig. 2, the right-hand side in the first quadrant is transformed into its symmetric image on the other side of the (2IzSz,ZQy) plane. The initial ρpre-LLSD and its plane-symmetric image ρLLSE have quite different relaxation rates triggered by different symmetry rules. The density operator evolves into

(9) ρ LLS E = | S 0 S 0 | - 1 2 ( | T - 1 T - 1 | + | T 1 T 1 | ) .

A note can be made on the time dependence of the relaxation rate constant for the case of free precession of the LLS. Given the evolution of zero-quantum coherences under the chemical shift difference, the density operator, starting from ZQx+IzSz, will have an oscillatory evolution between the two ZQx and ZQy as

(10) Z Q x + I z S z ( 2 π ν I I z + 2 π ν S S z ) τ Z Q x cos 2 π Δ ν I S τ + Z Q y sin 2 π Δ ν I S τ + I z S z .

Considering a relaxation superoperator only for the dipolar interaction between two coupled spins,

(11) Γ ^ = - 6 5 b D D 2 m = - 2 2 - 1 m J m × ω T 2 , m T 2 , - m _ ,

where bDD is the dipolar coupling constant, J(ω) is the spectral density, and T2,m are the spherical tensor spin operators of rank 2, the computed relaxation rate constant for the density operator ZQxcos2πΔνISτ+ZQysin2πΔνISτ+IzSz is

(12) R ( τ ) = - 6 5 b D D 2 ( 1 + cos 2 π Δ ν I S τ J ( 0 ) + 3 J ( ω ) ) ( sin 2 π Δ ν I S τ 2 ) 2 .

Thus, only for τ=2n2πΔνIS in the starting blocks is the relaxation rate constant of the obtained state optimally low.

Juggling with operators in order to drive the spin system in its “soft bed”, we realised we should always look at nature from various perspectives. We learned that if only one of the longitudinal two-spin-order and zero-quantum components could be selected at time point C (Fig. 1), QLLS would have been present already. For instance, ZZ magnetisation alone projects onto the long-lived state, given that 2IzSz=12(|T-1T-1|+|T1T1|-|T0T0|-|S0S0|). Thus, during the sustaining period, only the singlet population will survive for a period much longer than longitudinal magnetisation. In order to do so, we employed a Thrippleton–Keeler (Thrippleton and Keeler, 2003) filter to wipe out the troubling zero- and double-quantum coherences and obtained a singlet state with an amplitude 2 times lower than using both zero-quantum and ZZ magnetisation. Other groups employed the so-called “pseudo singlet order” (Pileio, 2017), which is just ZQx=12|T0T0|-|S0S0|, as the source of a long-lived state, obtaining similar results.

Though broadband excitation of singlet states would have deserved publication as a discovery in its own right, we were cautious to avoid publication of our research in slices of “salami science” (Sweedler, 2019). However, this advance proved relevant for the advancement of long-lived state order (Pileio, 2020; Teleanu et al., 2021; Pileio, 2017; Bengs et al., 2020) and was more challenging to obtain than the two-dimensional spectroscopy application for the study of singlet-state-based exchange we describe in the same paper (SS-EXSY) (Sarkar et al., 2007). In the tradition of finding low-key names for sequences such as “INEPT” (Morris and Freeman, 1979) or “INADEQUATE” (Bax et al., 1980), we could have named the zero-quantum rotation block of the pulse sequence in Fig. 1 a “(ZZ-)ZEROTATION”.

As singlet-based applications aim to store magnetisation for ever-longer time periods, the most adapted systems to this purpose, quasi-equivalent spin pairs with J couplings far overweighing the differences between the chemical shifts of the components, became increasingly studied. Spin dynamics that shift the magnetisation of the two spins differentially to create the singlet state are particularly challenging in such systems.

Pulse sequences of particular interest include the magnetisation to singlet (Tayler and Levitt, 2011) and spin-lock-induced crossover (DeVience et al., 2013). These methods are suited for strongly coupled spins where the scalar coupling is larger than the chemical shift difference, while the pulse sequence described in Fig. 2, which we identify herein as “ZZ+ZQx” (Sarkar et al., 2007), performs better in terms of long-lived state excitation in the weakly coupled regime. Several attempts to efficiently excite singlet order on broader domains of coupling regimes have been recently devised (Bengs et al., 2020; Mamone et al., 2020). Figure 3 depicts numerical simulations performed with SpinDynamica (Bengs and Levitt, 2018) for singlet population excitation using the aforementioned pulse sequences for both weakly and strongly coupled regimes, outlining the difference in excitation efficiency.

Figure 3Numerical simulations using SpinDynamica outlining the efficiency of singlet polarisation for M2S, SLIC and ZZ+ZQx pulse sequences in a two-spin system (by projecting the density operator during the pulse sequence onto the singlet population). The weakly coupled system is described by {ΔνIS=50Hz;JIS=17.4Hz}, while the strongly coupled system features {ΔνIS=2.8Hz;JIS=17.4Hz}. For each scenario, only the coherent evolution was considered in simulations (no relaxation dampening of the amplitude of source operators or LLS is taken into account). The excitation period is followed by a free-precession evolution after the maximum amplitude for the singlet population was reached.


3 Heteronuclei or proton long-lived states for conserving hyperpolarisation

In order to maximise the magnetisation lifetime, heteronuclear longitudinal spin order (mainly on carbon-13) can be excited and used during evolution periods in both room-temperature (Bermel et al., 2005b; Richter et al., 2010) and hyperpolarised NMR. We strove to also preserve hyperpolarisation on a pair of hydrogens entwined in a long-lived state (Vasos et al., 2009). Since the invention of dissolution–dynamic nuclear polarisation (dissolution–DNP) (Ardenkjær-Larsen et al., 2003; Comment et al., 2008; Balzan et al., 2016) and its development in Lausanne by the team of Arnaud Comment, Sami Jannin and Jac van der Klink in the Functional Imaging Laboratory at EPFL (Ardenkjær-Larsen et al., 2003; Comment et al., 2008; Balzan et al., 2016), the topic has been associated with our research due to its conjunction with long-lived spin order (Ardenkjær-Larsen et al., 2003; Comment et al., 2008; Balzan et al., 2016).

The preservation of hyperpolarised magnetisation obtained by dissolution–DNP in long-lived states raised fewer challenges than the comparison of the LLS with heteronuclear lifetimes in terms of performance as polarisation batteries. The hyperpolarised magnetisation in samples stemming from a polariser such as the one developed in Lausanne and the similar one installed as the first dissolution–DNP system in France (Ardenkjær-Larsen et al., 2003; Comment et al., 2008; Balzan et al., 2016) followed intently the rf pulses in our high-resolution magnets. Inhomogeneities due to fast dissolution and injection could be tamed to run pulse sequences for LLS excitation and decoding and to observe the signal. To us, the main issue remained that the benefits of a procedure consisting of QLLS excitation on protons compared to simply preserving hyperpolarised magnetisation in heteronuclei had to be carefully considered. We scrutinised this issue in terms of magnetisation lifetimes in the given conditions (room temperature, molecular size) as well as in the case of extreme molecular sizes or crowding of the environment and came to the conclusion that proton-based LLS were valuable for storing magnetisation even when the molecules contained isotope-enriched heteronuclei like 15N or 13C. In this analysis, we were inspired by similar comparisons between proton and heteronuclear magnetisation carried out for relaxation rates of heteronuclei and protons in large or paramagnetic proteins (Bermel et al., 2005a).

Journal editors understood the interest of the research topic, and the papers we sent for publication encountered no uphill Sisyphean battles (Molinié and Bodenhausen, 2013; The Myth of Sisyphus – Wikiwand,, last access: 21 September 2021, Evslin, 2006). The only significant delays in publication we incurred were for a paper submitted directly to a specialised journal that discussed results obtained using our freshly installed DNP system at the time (Balzan et al., 2017).

Figure 4(a) Pulse sequence designed to generate a two-dimensional correlation spectrum between 1H and 15N via scalar coupling of the type 2JNH and 3JNH starting from a proton long-lived state. (b) Schematic representation of polarisation transfer along a protein's backbone which generates a two-dimensional correlation spectrum via the 2JNH and 3JNH coupling constants. Protons are excited via the first part of the pulse sequence from (a) into a long-lived state spin order. Figure adapted from the reference Fernandes et al. (2013a).

Potential applications of symmetry-adapted states as magnetisation reservoirs for various two-dimensional experiments led us to investigate the polarisation transfer from long-lived states of protons towards heteronuclei across proteins' backbones (Fig. 4) (Fernandes et al., 2013a, b). The particular relaxation rates of LLS reflect conformational exchange and act as probes for unravelling proteins' inner dynamics, while two-dimensional correlations (e.g. 1H−15N) help disentangle complex spectra, which is acutely needed in the analysis of unfolded or intrinsically disordered proteins. The field of two-dimensional heteronuclear experiments using long-lived spin order (Fig. 4) is likely to develop further. As the field advances, long-lived state-based explorations of inter-molecular (Stavarache et al., 2017; Bornet et al., 2011) and intra-molecular (Ahuja et al., 2007; Tayler et al., 2010) interactions become important, and so becomes spectral resolution for the study of systems of increasing complexity.

4 Long-lived coherences, eppur' si muove

Aware of the potential of long magnetisation lifetimes for line narrowing in NMR, we strove to obtain some type of magnetisation akin to LLS to rotate. Long lifetimes of magnetisation had been traded for spectral resolution, e.g. for 1H−15N pairs in cross-correlated relaxation experiments and for 15N heteronuclei for narrowing spectroscopic lines (Goldman, 1984; Pervushin et al., 1997; Vasos et al., 2006).

Our search for “moving” long-lived configurations first involved complicated coherences in alanine, serine, and other molecules with up to five coupled spins (Ahuja et al., 2009). When we finally resorted to our favourite paired Gly aliphatic protons of AlaGly, the IxSx configuration was deduced from the diagonalised Liouvillian (Carravetta and Levitt, 2005). The next hurdle appeared in fitting the exotic long-lived states with products of oscillating and decaying functions and translating them to signals in two-dimensional spectra via Fourier transform-adapted spectroscopy (prior experience of non-conventional heteronuclear two-dimensional experiments (Bertini et al., 2004; Vasos et al., 2005) helped at this point). The simple scheme involving a 180 pulse for creating a QLLC observable in the indirect dimension of a two-dimensional experiment, where

(13) Q LLC = ( I x - S x ) cos ( 2 π J I S τ ) + ( 2 I y S z - 2 I z S y ) sin ( 2 π J I S τ ) ,

was further refined by different approaches (Fig. S1 in the Supplement), of which the most advanced generates coherences with spins pointing in opposite ways in molecules with almost-equivalent nuclei (Sheberstov et al., 2019).

When we sought long-lived coherences in high fields, we were not aware yet that the contemporary work on extremely low-frequency oscillations (Pileio et al., 2009) in low magnetic fields developed at Southampton involved, practically, the same operators.

Just as in the case of LLS for hyperpolarisation safekeeping was compared to heteronuclear storage, once LLCs were developed, comparisons with zero-quantum coherences came to mind. For the case of two J-coupled, non-equivalent spins, I and S, the evolution of ρ1(0)=Ix-Sx and ρZQx(0)=IxSx+IySy during free precession (without any radio-frequency “sustaining” applied) is given by

(14) ρ 1 ( τ ) = [ ( I x - S x ) cos ( π J I S τ ) + ( 2 I y S z - 2 I z S y ) sin ( π J I S τ ) ] cos ( 2 π Δ ν I S τ 2 ) - [ ( I y + S y ) cos ( π J I S τ ) - ( 2 I x S z + 2 I z S x ) sin ( π J I S τ ) ] sin ( 2 π Δ ν I S τ 2 ) ,

(15) ρ Z Q x τ = I x S x + I y S y cos 2 π Δ ν I S τ + I y S x - I x S y sin 2 π Δ ν I S τ = Z Q x cos ( 2 π Δ ν I S τ ) + Z Q y sin ( 2 π Δ ν I S τ ) ,

while in the presence of rf fields with the carrier placed in the middle of their offsets (υ1=-ΔυIS/2 and υ2=ΔυIS/2) and an amplitude υ1≫ΔυIS, the evolutions for the LLC and ZQx are

(16) ρ L L C ( τ ) = ( I x - S x ) cos ( 2 π J I S τ ) + ( I y S z - I z S y ) sin ( 2 π J I S τ ) ,

(17) ρ Z Q x ( τ ) = I x S x + I y S y cos 2 π υ 1 τ 2 + I z S z sin 2 π υ 1 τ 2 + ( I y S z + I z S y ) sin 2 π υ 1 τ cos 2 π υ 1 τ .

In the absence of sustaining rf fields, differences of single-quantum transverse coherences, the source of LLCs, evolve under the chemical shift difference and the scalar coupling, while ZQs are immune to the latter (Cavanagh et al., 1995). In the second scenario, LLCs evolve only under the effect of J coupling (Sarkar et al., 2010), oscillating between in-phase (IxSx) and anti-phase (2IySz−2IzSy), with the coherence order equal to 1. Thus, the scalar-coupling evolution sets LLCs apart from ZQs (Fig. 5).

Figure 5Evolution of zero-quantum coherences during free precession with an oscillating frequency of ΔυIS and evolution of long-lived coherences during a sustaining period with an oscillating frequency equal to the scalar-coupling constant JIS.


Broadband excitation of LLCs in molecules with broadly different J couplings and chemical shifts is still a challenge, despite the progress. We explored part of the territory by exciting with a series of selective 180 pulses (Sarkar et al., 2011) and by sustaining with various pulse trains (Sadet et al., 2014), but we can safely say that LLCs benefit from their simple and parametric-free excitation scheme which consists of a selective π pulse and non-selective π/2 hard pulse followed by spin lock.

5 Conclusions

We present from today's perspective several challenging aspects in the introduction of coherent dynamics designed to render spin order resilient to Chronos' decrees. The focus is placed on zero-quantum inversion with respect to longitudinal two-spin order in homonuclear spin systems. Numerical simulations outlining the efficiency of different pulse sequences in creating long-lived states in different coupling regimes are discussed. The presented methods are delineated in a manner designed to render them useful as building blocks in further applications.

Code availability

The Mathematica notebook used for the simulations is provided as a Supplement.

Data availability

The only data are the notebooks used for simulation which were deposited on ResearchGate as (Teleanu, 2021). The link was also added in the Supplement.


The supplement related to this article is available online at:

Author contributions

FT and PRV designed the study, FT performed the numerical simulations, and FT and PRV wrote the paper.

Competing interests

The authors declare that they have no conflict of interest.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special issue statement

This article is part of the special issue “Geoffrey Bodenhausen Festschrift”. It is not associated with a conference.


The authors thank Diana Serafin and Ioana Fidel for recent discussions in Bucharest as well as Riddhiman Sarkar, Karthik Gopalakrishnan, Simone Cavadini, Jens Dittmer, and Adonis Lupulescu for the discussions in Lausanne. The authors acknowledge the financial support ensured by UEFISCDI, the Romanian Ministry of Research, the Romanian Government and the European Union.

This article is dedicated to Geoffrey Bodenhausen for his 70th birthday.

Financial support

This research has been supported by the UEFISCDI (grant nos. PN-III-P4-ID-PCE-2020-2642, PN-III-P2-2.1-PED-2019-4212 and 545PED/2020), the Romanian Ministry of Research (core project 342 PN 1906 01 05/2019), the Romanian Government through Extreme Light Infrastructure Nuclear Physics (ELI-NP) Phase II project 343, and the European Union through the European Regional 344 Development Fund and the Competitiveness Operational Programme (1/07.07.2016, ID 1334).

Review statement

This paper was edited by Daniel Abergel and reviewed by Malcolm Levitt and Thomas Eykyn.


Ahuja, P., Sarkar, R., Vasos, P. R., and Bodenhausen, G.: Molecular properties determined from the relaxation of long-lived spin states, J. Chem. Phys., 127, 134112,, 2007. 

Ahuja, P., Sarkar, R., Vasos, P. R., and Bodenhausen, G.: Long-lived States in Multiple-Spin Systems, ChemPhysChem, 10, 2217–2220,, 2009. 

Ahuja, P., Sarkar, R., Jannin, S., Vasos, P. R., and Bodenhausen, G.: Proton hyperpolarisation preserved in long-lived states, Chem. Commun., 46, 8192–8194,, 2010. 

Ardenkjær-Larsen, J. H., Fridlund, B., Gram, A., Hansson, G., Hansson, L., Lerche, M. H., Servin, R., Thaning, M., and Golman, K.: Increase in signal-to-noise ratio of >10,000 times in liquid-state NMR, P. Natl. Acad. Sci. USA, 100, 10158–10163,, 2003. 

Balzan, R., Fernandes, L., Comment, A., Pidial, L., Tavitian, B., and Vasos, P. R.: Dissolution Dynamic Nuclear Polarization Instrumentation for Real-time Enzymatic Reaction Rate Measurements by NMR, J. Vis. Exp., 108, e53548,, 2016. 

Balzan, R., Fernandes, L., Pidial, L., Comment, A., Tavitian, B., and Vasos, P. R.: Pyruvate cellular uptake and enzymatic conversion probed by dissolution DNP-NMR: the impact of overexpressed membrane transporters, Magn. Reson. Chem., 55, 579–583,, 2017. 

Bax, A., Freeman, R., and Kempsell, S. P.: Natural abundance carbon-13-carbon-13 coupling observed via double-quantum coherence, J. Am. Chem. Soc., 102, 4849–4851,, 1980. 

Bengs, C. and Levitt, M. H.: SpinDynamica: Symbolic and numerical magnetic resonance in a Mathematica environment, Magn. Reson. Chem., 56, 374–414,, 2018. 

Bengs, C., Sabba, M., Jerschow, A., and Levitt, M. H.: Generalised magnetisation-to-singlet-order transfer in nuclear magnetic resonance, Phys. Chem. Chem. Phys., 22, 9703–9712,, 2020. 

Bermel, W., Bertini, I., Felli, I. C., Pierattelli, R., and Vasos, P. R.: A selective experiment for the sequential protein backbone assignment from 3D heteronuclear spectra, J. Magn. Reson., 172, 324–328, 2005a. 

Bermel, W., Bertini, I., Duma, L., Felli, I. C., Emsley, L., Pierattelli, R., and Vasos, P. R.: Complete assignment of heteronuclear protein resonances by protonless NMR spectroscopy, Angew. Chem. Int. Edit., 44, 3089–3092, 2005b. 

Bertini, I., Duma, L., Felli, I. C., Fey, M., Luchinat, C., Pierattelli, R., and Vasos, P. R.: A heteronuclear direct-detection NMR spectroscopy experiment for protein-backbone assignment, Angew. Chem. Int. Edit., 43, 2257–2259, 2004. 

Bornet, A., Jannin, S., Konter, J. A. (Ton), Hautle, P., van den Brandt, B., and Bodenhausen, G.: Ultra High-Resolution NMR: Sustained Induction Decays of Long-Lived Coherences, J. Am. Chem. Soc., 133, 15644–15649,, 2011. 

Carravetta, M. and Levitt, M. H.: Long-Lived Nuclear Spin States in High-Field Solution NMR, J. Am. Chem. Soc., 126, 6228–6229,, 2004. 

Carravetta, M. and Levitt, M. H.: Theory of long-lived nuclear spin states in solution nuclear magnetic resonance. I. Singlet states in low magnetic field, J. Chem. Phys., 122, 214505,, 2005. 

Carravetta, M., Johannessen, O. G., and Levitt, M. H.: Beyond the T1 Limit: Singlet Nuclear Spin States in Low Magnetic Fields, Phys. Rev. Lett., 92, 153003,, 2004. 

Cavanagh, J., Fairbrother, W. J., Palmer III, A. G., and Skelton, N. J.: Protein NMR Spectroscopy: Principles and Practice, Elsevier, 915 pp., ISBN: 978-0-08-051529-8, 1995. 

Comment, A., van den Brandt, B., Uffmann, K., Kurdzesau, F., Jannin, S., Konter, J. A., Hautle, P., Wenckebach, W. T., Gruetter, R., and van der Klink, J. J.: Principles of Operation of a DNP Prepolarizer Coupled to a Rodent MRI Scanner, Appl. Magn. Reson., 34, 313–319,, 2008. 

DeVience, S. J., Walsworth, R. L., and Rosen, M. S.: Preparation of Nuclear Spin Singlet States Using Spin-Lock Induced Crossing, Phys. Rev. Lett., 111, 173002,, 2013. 

Ernst, R. R.: Nobel Lecture, Nobel Media AB 2021, available at:, last access: 14 March 2021. 

Ernst, R. R. and Anderson, W. A.: Application of Fourier Transform Spectroscopy to Magnetic Resonance, Rev. Sci. Instrum., 37, 93–102,, 1966. 

Evslin, B.: Gods, Demigods and Demons: A Handbook of Greek Mythology, Bloomsbury Academic, ISBN 978-1-84511-321-6, 2006. 

Fernandes, L., Sadet, A., Bolopion, A., Beaune, P., Kateb, F., and Vasos, P. R.: Hauts champs, pour la RMN des protéines désordonnées et la métabolomique de milieux complexes, Spectra Analyse, 295, 42–48, 2013a (in French). 

Fernandes, L., Guerniou, C., Marín-Montesinos, I., Pons, M., Kateb, F., and Vasos, P. R.: Long-lived states in an intrinsically disordered protein domain, Magn. Reson. Chem., 51, 729–733,, 2013b. 

Ferrage, F., Zoonens, M., Warschawski, D. E., Popot, J.-L., and Bodenhausen, G.: Slow Diffusion of Macromolecular Assemblies by a New Pulsed Field Gradient NMR Method, J. Am. Chem. Soc., 125, 2541–2545,, 2003. 

Goldman, M.: Interference effects in the relaxation of a pair of unlike spin-12 nuclei, Journal of Magnetic Resonance (1969), 60, 437–452,, 1984. 

Mamone, S., Rezaei-Ghaleh, N., Opazo, F., Griesinger, C., and Glöggler, S.: Singlet-filtered NMR spectroscopy, Science Advances, 6, eaaz1955,, 2020. 

Molinié, A. and Bodenhausen, G.: On toxic effects of scientific journals, J. Biosci., 38, 189–199,, 2013. 

Morris, G. A. and Freeman, R.: Enhancement of nuclear magnetic resonance signals by polarization transfer, J. Am. Chem. Soc., 101, 760–762,, 1979. 

Pervushin, K., Riek, R., Wider, G., and Wüthrich, K.: Attenuated T2 relaxation by mutual cancellation of dipole–dipole coupling and chemical shift anisotropy indicates an avenue to NMR structures of very large biological macromolecules in solution, P. Natl. Acad. Sci. USA, 94, 12366–12371,, 1997. 

Pileio, G.: Singlet NMR methodology in two-spin-1/2 systems, Prog. Nucl. Mag. Res. Sp., 98–99, 1–19,, 2017. 

Pileio, G.: Long-lived Nuclear Spin Order: Theory and Applications, Royal Society of Chemistry, 462 pp., ISBN: 978-1-78801-568-4, 2020. 

Pileio, G., Carravetta, M., and Levitt, M. H.: Extremely Low-Frequency Spectroscopy in Low-Field Nuclear Magnetic Resonance, Phys. Rev. Lett., 103, 083002,, 2009. 

Richter, C., Kovacs, H., Buck, J., Wacker, A., Fürtig, B., Bermel, W., and Schwalbe, H.: 13C-direct detected NMR experiments for the sequential J-based resonance assignment of RNA oligonucleotides, J. Biomol. NMR, 47, 259–269,, 2010. 

Sadet, A., Fernandes, L., Kateb, F., Balzan, R., and Vasos, P. R.: Long-lived coherences: Improved dispersion in the frequency domain using continuous-wave and reduced-power windowed sustaining irradiation, J. Chem. Phys., 141, 054203,, 2014. 

Sarkar, R., Vasos, P. R., and Bodenhausen, G.: Singlet-State Exchange NMR Spectroscopy for the Study of Very Slow Dynamic Processes, J. Am. Chem. Soc., 129, 328–334,, 2007. 

Sarkar, R., Ahuja, P., Vasos, P. R., and Bodenhausen, G.: Long-Lived Coherences for Homogeneous Line Narrowing in Spectroscopy, Phys. Rev. Lett., 104, 053001,, 2010. 

Sarkar, R., Ahuja, P., Vasos, P. R., Bornet, A., Wagnières, O., and Bodenhausen, G.: Long-lived coherences for line-narrowing in high-field NMR, Prog. Nucl. Mag. Res. Sp., 59, 83–90,, 2011. 

Sheberstov, K. F., Kiryutin, A. S., Bengs, C., Hill-Cousins, J. T., Brown, L. J., Brown, R. C. D., Pileio, G., Levitt, M. H., Yurkovskaya, A. V., and Ivanov, K. L.: Excitation of singlet–triplet coherences in pairs of nearly-equivalent spins, Phys. Chem. Chem. Phys., 21, 6087–6100,, 2019. 

Solomon, I.: Relaxation Processes in a System of Two Spins, Phys. Rev., 99, 559–565,, 1955. 

Sørensen, O. W., Eich, G. W., Levitt, M. H., Bodenhausen, G., and Ernst, R. R.: Product operator formalism for the description of NMR pulse experiments, Prog. Nucl. Mag. Res. Sp., 16, 163–192,, 1984. 

Stavarache, C., Hanganu, A., Paun, A., Paraschivescu, C., Matache, M., and Vasos, P. R.: Long-lived states detect interactions between small molecules and diamagnetic metal ions, J. Magn. Reson., 284, 15–19, 2017. 

Stevanato, G.: Long-lived states in multi-spin systems, PhD thesis, University of Southampton, UK, 188 pp., 2015. 

Stevanato, G.: Chapter 2: Symmetry and Long-lived Spin Order, in: Long-lived Nuclear Spin Order, Royal Society of Chemistry, 33–63,, 2020. 

Sweedler, J. V.: Salami Publications and Duplicate Submissions: Put Them on Your List of Things to Avoid, Anal. Chem., 91, 3177–3178,, 2019. 

Tayler, M. C. D. and Levitt, M. H.: Singlet nuclear magnetic resonance of nearly-equivalent spins, Phys. Chem. Chem. Phys., 13, 5556–5560,, 2011. 

Tayler, M. C. D., Marie, S., Ganesan, A., and Levitt, M. H.: Determination of Molecular Torsion Angles Using Nuclear Singlet Relaxation, J. Am. Chem. Soc., 132, 8225–8227,, 2010. 

Teleanu, F.: M2S SLIC ZZZQx comparison weak strong coupling, ResearchGate [data set],, 2021. 

Teleanu, F., Sadet, A., and Vasos, P. R.: Symmetry versus entropy: Long-lived states and coherences, Prog. Nucl. Mag. Res. Sp., 122, 63–75,, 2021. 

Thrippleton, M. J. and Keeler, J.: Elimination of Zero-Quantum Interference in Two-Dimensional NMR Spectra, Angew. Chem. Int. Edit., 42, 3938–3941,, 2003. 

Vasos, P. R., Hall, J. B., and Fushman, D.: Spin-state selection for increased confidence in cross-correlation rates measurements, J. Biomol. NMR, 31, 149–154, 2005. 

Vasos, P. R., Hall, J. B., Kuemmerle, R., and Fushman, D.: Measurement of 15N relaxation in deuterated amide groups in proteins using direct nitrogen detection, J. Biomol. NMR, 36, 27–36,, 2006. 

Vasos, P. R., Comment, A., Sarkar, R., Ahuja, P., Jannin, S., Ansermet, J.-P., Konter, J. A., Hautle, P., van den Brandt, B., and Bodenhausen, G.: Long-lived states to sustain hyperpolarized magnetization, P. Natl. Acad. Sci. USA, 106, 18469–18473,, 2009. 

Short summary
Long-lived nuclear spin order helps extend the timescale of polarisation storage. This occurs as populations of nuclear spin singlet states feature slow decays compared to longitudinal magnetisation. The excitation of these symmetry-adapted states required spin choreographies that were often overshadowed in scholarly presentations by their applications. We recapitulate some of the milestones achieved in designing methods for long-lived spin order excitation and comment on future perspectives.