Travaux Dirigés 2e Ecole du GERM Jean-Marc Nuzillard

Cargèse, 2008

Principes de base

Exercice 1.
Un spectromètre de RMN opère à 500 MHz pour l’observation du proton.
Sachant que γ(1 H) = 26,75.107 s−1 .T−1 et que γ(113 Cd) = -5,96.107 s−1 .T−1 , calculer la
fréquence d’observation du noyau 113 Cd.

La relation qui donne la fréquence de résonance ν en fonction du champ appliqué B0 et

le rapport gyromagnétique γ est

γB0 γ(113 Cd)

ν=− et donc ν(113 Cd) = ν(1 H) × = −111.4 MHz
2π γ(1 H)
−
→
Le signe négatif du résultat indique qu’au même champ B0 la précession de Larmor des
noyaux 1 H et 113 Cd s’effectue dans des sens opposés.

Exercice 2.
Deux noyaux, 1 et 2, ont pour déplacement chimique ν1 = 2.5 ppm et ν2 = 3.0 ppm.
L’intensité de leur couplage scalaire est de 12 Hz.
Sont-ils faiblement ou fortement couplés selon que la fréquence d’observation est 600 MHz
ou 60 MHz ?

Le déplacement chimique δ d’un noyau de fréquence de résonance ν défini par rapport à

une substance de référence est défini par :
ν − ν ref
δ= × 106 (1)
ν ref
Une différence de déplacement chimique ∆δ correspond donc à une différence de fréquence
de résonance ∆ν telle que
106 ν ref
∆δ = ∆ν × ou ∆ν = ∆δ ×
ν ref 106
A 600 MHz, ∆ν = 0,5 × 600 = 300 Hz. En considérant que |J| ≪ |∆ν|, le couplage est faible.
A 60 MHz, ∆ν = 30 Hz = 2.5 J. Le couplage peut alors être considéré comme fort.

Exercice 3.
Un échantillon est excité par une impulsion de fréquence ν RF .
Quelle est la fréquence et le sens de rotation du référentiel tournant par rapport au
référentiel du laboratoire ?
L’intensité du champ magnétique de cette impulsion est de ”25 kHz”. Que cela signifie-t-
il ?

1
 Que vaut B1max , l’amplitude maximale du champ magnétique d’excitation (en RMN du
1 H) ?

Combien de temps faut-il pour réaliser une impulsion d’angle du nutation égal à 90˚ ?
De combien de dB faut-il atténuer la puissance de l’amplificateur d’émission pour ce temps
soit doublé ?

Le réfrentiel tourant tourne autour de l’axe z du référentiel du laboratoire à la fréquence

ν RF , dans le sens de la précession de Larmor.
Tant que dure l’impulsion de RF, l’aimantation macroscopique d’un noyau en résonance
tourne, dans le cas présent, à une fréquence ν1 de 25 kHz autour d’un axe horizontal du
−→
référentiel tournant. Le champ B1 polarisé linéairement a pour expression dans le référentiel
du laboratoire :
−
→
B1 = B1max cos(2πν RF t + ϕ)~ı
qui est la somme vectorielle de deux champs tournants en sens opposés aux fréquences ±ν RF
et d’intensités égales à B1max /2. La résolution des équations de Bloch dans le référentiel
tournant indique que
γB1max /2
ν1 =
2π
et donc que
4πν1
B1max = = 1,17 mT
γ
L’aimantation qui tourne à la fréquence ν1 pendant l’impulsion fait un tour complet en
1/ν1 = 40 µs et donc tourne de 90˚en 10 µs.
Pour multiplier ce temps par un facteur f en gardant l’angle de nutation inchangé, il faut
diviser ν1 par f et donc diviser par f l’intensité du courant électrique qui parcours la bobine
RF. Pour ce faire, il faut aussi diviser par f la tension appliquée à la bobine et donc diviser
par f 2 la puissance fournie par l’amplificateur qui délivre l’impulsion RF. Cela est équivalent
à atténuer cette puissance de

atténuation = 10 log f 2 = 20 log f ≃ 6 dB.

Exercice 4.
On considère une impulsion d’angle de nutation π/2 radians et durée 100 µs.
Pour quelle valeur d’offset (exprimée en Hz) un noyau verra-t-il son aimantation tourner
d’un angle de 2π ?
Quel est l’intérêt d’un tel calcul ?

Pour cette impulsion de durée τ = 100 µs :

2πν1 τ = Ω1 τ = π/2.

Un noyau d’offset ν0 = Ω0 /2π soumis à la même impulsion voit son aimantation tourner d’un
angle q
Ωeff τ = Ω20 + Ω21 τ = 2π
On en déduit
√
√ 15
Ω20 + Ω21 = 16 Ω21 et donc ν0 = 15 ν1 = ≃ 9,7 kHz.
4τ
A 100 MHz en RMN du 13 C, cela représente environ 100 ppm de différence de déplacement
chimique entre le centre du spectre et le noyau considéré. Il est ainsi possible d’exciter les

2
noyaux 13 CO sans toucher les 13 Cα dans les protéines.

Exercice 5.
Quels sont l’angle de rotation et la phase de l’impulsion composite définie par la suite
d’impulsions π/2(x)π(y)π/2(x) ?

Il faut savoir comment les axes X, Y et Z du référentiel tournant sont transformés par
l’impulsion composite :
π/2(x) π(y) π/2(x)
X −→ X −→ −X −→ −X
Y −→ Z −→ −Z −→ Y
Z −→ −Y −→ −Y −→ −Z
L’impulsion composite laisse l’aimantation sur Y invariante et inverse le sens de celle sur X
et sur Z. Il s’agit donc d’une impulsion π(y).
La règle de transformation utilisée ci-dessus est la suivante :
– Si l’axe de rotation associé à l’impulsion est colinéaire avec l’aimantation, rien ne passe.
– Dans le cas contraire
– une impulsion d’angle π et de phase quelconque inverse le sens de l’aimantation.
– une impulsion d’angle π/2 effectue les transformations suivante :

π/2(x) : Y −→ Z Z −→ −Y π/2(y) : X −→ −Z Z −→ X

Le sens de l’axe final est donné par la règle du tire-bouchon. En prenant le premier cas
comme exemple, une rotation d’angle π/2 et d’axe X correspond à une translation
du tire-bouchon dans le sens des X positifs qui amène effectivement l’axe Y sur l’axe
Z. Cela n’est vrai que si les axes X, Y et Z, pris dans cet ordre, forment un trièdre
d’orientation directe et que si le modèle du tire-bouchon est standard.

Exercice 6.
Une impulsion réelle d’angle π/2 radians et de durée τ peut être modélisée au premier
ordre par une impulsion idéale de même angle suivie d’un délai de durée δ.
Exprimer δ en fonction de τ .
Quel est l’effet de ce délai sur le spectre produit par transformation de Fourier du FID
obtenu ?

On considère une impulsion π/2(y) d’intensité Ω1 et de durée τ qui agit sur un noyau
d’offset Ω0 ≪ Ω1 . La pulsation Ω1 est reliée à la durée de l’impulsion par

Ω1 τ = π/2

La fréquence angulaire effective Ωeff de rotation de l’aimantation d’équilibre autour de la

direction du champ effectif est très proche de Ω1 et considérée comme égale à Ω1 au premier
ordre : q q
Ωeff = Ω21 + Ω20 = Ω1 1 + Ω20 /Ω21 ≃ Ω1 (1 + Ω20 /2Ω21 )
L’axe du champ effectif est très proche de l’axe Y et fait avec dernier un angle ϕ tel que

tan ϕ = Ω0 /Ω1 ≃ ϕ
−→
L’angle entre l’axe de rotation et l’aimantation initiale M0 est π/2 − ϕ, qui est aussi (ap-
−→ −→
proximativement) l’angle entre l’aimantation finale M1 et l’axe Y . Au premier ordre M1 est
dans le plan XY et fait un angle ϕ avec l’axe X.

3
 Z Z
−→
M0

90˚
−−
→ −−
→
Ωeff Ωeff
ϕ
Y Y
ϕ
X X −→
M1

−→
Une impulsion parfaite π/2(y) amène M0 exactement sur l’axe X. Pendant un délai δ,
cette aimantation tourne d’un angle ϕ = Ω0 δ dans le plan XY . Ainsi :
ϕ 1 2
δ= = =τ×
Ω0 Ω1 π
L’angle qui repère la position finale de l’aimantation dans le plan XY dépend linéairement
de l’offset Ω0 lorsque Ω0 ≪ Ω1 , ce qui est interprétable par l’existence du délai δ ”inclus”
dans l’impulsion. A ce délai s’ajoute celui qui est inséré entre la fin de l’impulsion et le début
de l’acquisition du signal. Le délai global introduit un déphasage de l’aimantation et donc
du signal qui dépend de l’offset du noyau étudié. Cela conduit à effectuer une correction de
phase dans le spectre qui dépende linéairement du déplacement chimique.

Exercice 7.
L’aimantation d’équilibre d’un échantillon est soumise à l’action d’une impulsion de ra-
diofréquence d’angle de nutation α et de phase ϕ nulle.
Calculer dans le référentiel tournant les coordonnées du vecteur aimantation juste après
l’impulsion.
Pour un noyau d’offset Ω0 (en rad.s−1 ), que deviennent ces coordonnées après t secondes
d’évolution libre ? Les calculs seront effectués soit géométriquement soit à l’aide du formalisme
des opérateurs cartésiens.
Donner l’expression du signal de FID s(t) en tenant compte de la relaxation et des inho-
mogénéités du champ B0 .
Répondre aux mêmes questions pour une impulsion de phase ϕ quelconque.
−→
La phase de l’impusion étant nulle, l’aimantation initiale M0 tourne d’un angle α autour
de l’axe X du référentiel tournant, et se dirige vers l’axe −Y :
−→
M1 (0, −M0 sin α, M0 cos α)

4
 Z Y
−→
M0
−→ −→
M1 M2
M0 cos α −M0 sin α cos(Ω0 t)
α

−Y Y −X X
Ω0 t
X Z
−M0 sin α M0 sin α sin(Ω0 t)
−M0 sin α −→
M1
−Z −Y
Nutation d’angle α Précession d’angle Ω0 t
−→
Pendant la période de précession, M1 tourne d’un angle Ω0 t autour de l’axe Z pour donner
−→ −
→
M2 . La composante de M sur l’axe Z est invariante et celle sur l’axe −Y se dirige vers l’axe
X :
−→
M2 (M0 sin α sin(Ω0 t), −M0 sin α cos(Ω0 t), M0 cos α) (2)
S’agissant de l’aimantation macroscopique de noyaux isolés, le formalisme des opérateurs
cartésiens fournit un résultat rigoureusement identique. L’état initial du système est

σ 0 = Iz

Pendant l’impulsion de durée τ , d’intensité Ω1 et de phase nulle, l’hamiltonien HRF agit sur
le système :
HRF = Ω1 Ix avec Ω1 τ = α.
Les opérateurs HRF et σ0 ne commutent pas et l’état σ1 du système après l’impulsion est
donné par
σ1 = Iz cos α + {Ix , Iz } sin α = Iz cos α − Iy sin α
La notation {A, B} est définie par
1
{A, B} = [A, B]
i
où [A, B] = AB − BA est le commutateur des opérateurs A et B. Cette notation évite
d’introduire des nombres complexes pendant les calculs d’évolution des systèmes, quand ils
sont effectués dans la base des opérateurs cartésiens.
Pendant la période de précession libre

Hevo = Ω0 Iz .

Puisque [Hevo , Iz ] = 0 et [Hevo , Iy ] 6= 0, l’état du système evolue pendant le temps t de l’état

σ1 vers l’état σ2 selon :

σ2 = Iz cos α − sin α(Iy cos(Ω0 t) + {Iz , Iy } sin(Ω0 t))

(3)
= sin α sin(Ω0 t)Ix − sin α cos(Ω0 t)Iy + cos αIz

La similitude formelle des équations 2 et 3 traduit, pour un ensemble de noyaux isolés,

l’équivalence des approches de la dynamique des spins par le modèle vectoriel et par l’opérateur
densité.
Le signal de FID s(t) enregistré au cours de l’évolution libre du système est identifié, à
un facteur multiplicatif près, avec le nombre complexe Mx (t) + iMy (t), où l’aimantation est
mesurée dans le référentiel tournant. Le modèle vectoriel fournit explicitement les valeurs de

5
Mx et de My , alors que ces grandeurs sont respectivement proportionnelles aux coefficients
multiplicatifs des opérateurs Ix et Iy . Ainsi

s(t) = sin α sin(Ω0 t) − i sin α cos(Ω0 t)

= −i sin α(cos(Ω0 t) + i sin(Ω0 t))
= −i sin α exp(iΩ0 t)

Pour tenir compte de la décroissance de l’aimantation transversale sous l’effet conjugué de la

relaxation transversale et de l’inhomogénéité du champ B0 , il faut multiplier l’expression de
∗
s(t) précédemment obtenue par e−t/T2 :
∗
s(t) = −i sin α exp(iΩ0 t) e−t/T2

où T2∗ est le temps de relaxation transversal apparent.

Lorsque la phase de l’impulsion est nulle, l’aimantation transversale fait un angle initial
(à t = 0) égal à −π/2 avec l’axe X, ce qui explique la présence du facteur multiplicatif −i
dans l’expression de s(t). Si la phase de l’impulsion est augmentée de ϕ (et donc égale à
−
→
ϕ), l’angle entre M (t = 0) et l’axe X est augmenté d’autant. Cela n’est vrai que parce que
l’aimantation initiale est alignée sur l’axe Z. Le signal est alors multiplié par exp(iϕ) :
∗
s(t) = −i exp(iϕ) exp(iΩ0 t) e−t/T2
(4)
= ei(ϕ−π/2) exp(iΩ0 t) e−t/T2
∗

Du point de vue des opérateurs, faire tourner l’aimantation d’un angle ϕ autour de l’axe Z
est équivalent à appliquer fictivement un opérateur ΩIz tendant un temps τ tel que Ωτ = ϕ.
D’une manière générale, le résultat de l’action d’un opérateur H = ΩA pendant le temps
τ sur un état B ne dépend que de B et du produit Hτ = Ωτ A. Par abus de langage, l’action
de H pendant le temps τ sera par la suite désignée comme l’action de l’opérateur Hτ :

ΩA Ωτ A
−−−→ ⇐⇒ −−−−→
τ
Ceci se justifie par la relation générale

σ(τ ) = exp(−iHτ ) σ(τ = 0) exp(+iHτ )

Ainsi, faire agir successivement sur σ1 les opérateurs Ω0 tIz et ϕIz est équivalent à faire
agir l’opérateur (Ω0 t + ϕ)Iz :

σ2 = sin α sin(Ω0 t + ϕ)Ix − sin α cos(Ω0 t + ϕ)Iy + cos αIz

qui conduit à l’expression de s(t) donnée par la relation 4.

Exercice 8.
On considère un système de deux noyaux couplés I et S.
Ecrire l’expression de l’opérateur d’évolution libre de ce système. Tracer l’évolution des
états Ix,y,z et 2Ix,y,z Sx,y,z sous l’action de cet opérateur pendant un temps t.
Calculer la FID correspondante et le spectre associé, s’il y a lieu.
Justifier la nomenclature usuelle des états.

On désignera par ΩI et ΩS les offset des noyaux I et S dans leurs référentiels tournants
respectifs (qui sont confondus si le système est homonucléaire). La grandeur J désigne l’inten-
sité (en Hz) du couplage scalaire partagé par I et S. L’opérateur d’évolution libre du système
s’écrit
Hevo = ΩI Iz + ΩS Sz + πJ2Iz Sz

6
Les trois termes de cet opérateur commutent entre eux et peuvent donc être appliqués suc-
cessivement à un état initial dans un ordre arbitrairement ou judicieusement choisi.
Les opérateur Iz , Sz et 2Iz Sz sont invariants par action de H car ils commutent avec
chacun de ses termes. Ces trois opérateurs, auxquels il faut adjoindre l’opérateur unité (inva-
riant par toute transformation), décrivent les populations des quatre états de spins associés
aux états quantiques définis par mIs = ±1/2 et mSs = ±1/2. Ils ne peuvent pas évoluer sous
l’action de Hevo vers de l’aimantation transversale détectable.
L’évolution libre de chacun des états Ix , Iy , 2Ix Sz et 2Iy Sz conduit à une combinaison
linéaire de ces mêmes états. Considérons l’évolution de Iz :

Ix
Le calcul est conduit ici graphique-
ΩS t · Sz ment. Une double flèche indique que
Ix l’opérateur d’évolution commute avec
l’état auquel il s’applique. Dans le cas
πJt · 2Iz Sz contraire, un opérateur ΩA qui agit
pendant le temps t sur un opérateur B
Ix 2Iy Sz
conduit à la somme de cos(Ωt)B (flèche
Ω I t · Iz de gauche) et de sin(Ωt){A, B} (flèche
de droite).
Ix Iy 2Iy Sz −2Ix Sz

Au cours de l’évolution, Ix varie comme

1
cos(πJt) cos(ΩI t) = (cos((ΩI − πJ)t) + cos((ΩI + πJ)t))
2
c’est-à-dire aux pulsations ±(ΩI ± πJ) des transitions à ±1 quanta du noyau I, qui corres-
pondent aux fréquences ±(νI ± J/2). Pour cette raison, les états Ix , Iy , 2Ix Sz et 2Iy Sz sont
désignés comme états à ±1 quanta du noyau I.
Si le temps t est celui de l’acquisition du signal à partir de l’état Ix à t = 0, alors

sx (t) = cos(πJt) cos(ΩI t)

sy (t) = cos(πJt) sin(ΩI t)
s(t) = sx (t) + isy (t)
= cos(πJt) exp(iΩI t)
1
= [exp(+iπJt) + exp(−iπJt)] exp(iΩI t)
2
1
= [exp (i(ΩI + πJ)t) + exp (i(ΩI − πJ)t)]
2
La TF de ce signal fait apparaı̂tre un doublet de raies d’intensités relatives 12 aux fréquences
νI ± J/2.
L’évolution libre de l’état Iy est qualitativement comparable à celle de Ix , à ceci près que

1
s(t) = i [exp (i(ΩI + πJ)t) + exp (i(ΩI − πJ)t)]
2
Le déphasage de π/2 de l’aimantation initiale de I dans le plan XY se traduit par un
déphasage identique dans le signal et dans le spectre obtenu par TF.
L’évolution libre de l’état 2Ix Sz est donnée par le graphe suivant :

7
 2IxSz
ΩS t · Sz Les symboles  indiquent qu’il n’y a
2IxSz pas besoin d’écrire les opérateurs en
question si le but du calcul est de
πJt · 2Iz Sz déterminer le signal observable. En ef-
fet, 2Ix Sz et {Iz , 2Ix Sz } ne sont pas
2Ix Sz Iy
détectables. Seuls les opérateurs Ix , Iy ,
Ω I t · Iz Sx et Sy le sont.

  Iy −Ix

ce qui conduit au signal mesuré

s(t) = i sin(πJt) exp(iΩI t)
1
= [exp (i(ΩI + πJ)t) − exp (i(ΩI − πJ)t)]
2
qui donne après TF un doublet de raies d’intensités relatives 12 et − 12 aux fréquences νI ±J/2.
Un tel doublet est dit ”antiphase”, par opposition à celui issu de l’évolution de Ix qui est ”en
phase”.
L’évolution libre de 2Iy Sz est qualitativement comparable à celle de 2Ix Sy et conduit au
même signal, à un facteur i près :

s(t) = − sin(πJt) exp(iΩI t)

L’évolution libre de chacun des états 2Ix Sx , 2Iy Sx , 2Ix Sy et 2Iy Sy conduit à une combi-
naison linéaire de ces mêmes états. Considérons l’évolution de 2Ix Sx :

2IxSx
πJt · 2Iz Sz

2IxSx
L’évolution libre des états Ix,y Sx,y ne
Ω I t · Iz conduit à aucune aimantation obser-
vable.
2IxSx 2Iy Sx
ΩS t · Sz

2IxSx 2Ix Sy 2Iy Sx 2Iy Sy

L’état 2Ix Sx est modulé au cours de son évolution par le facteur

1
cos(ΩI t) cos(ΩS t) = [cos ((ΩI + ΩS )t) + cos ((ΩI − ΩS )t)]
2
qui exprime une variation aux pulsations ±(ΩI ± ΩS ), c’est-à-dire aux fréquences des tran-
sitions à zéro et double quanta, respectivement égales à ±(νI − νS ) et ±(νI + νS ). Les états
Ix,y Sx,y sont dits ”à zero et double quanta”.

Exercice 9.
On considère un système hétéronucléaire de deux noyaux couplés I et S. Tracer les trois
séquences d’écho de spin possibles.
Donner dans chaque cas l’expression de l’opérateur d’évolution réduit qui lui est associé
et indiquer comment il s’utilise.

8
 La séquence de refocalisation peut concerner le noyau I seul (cas A) ou le noyau S seul
(cas B) ou les deux à fois (cas C).

A I : —τ —π—τ — B I : —τ — —τ — C I : —τ —π—τ —
S : —τ — —τ — S : —τ —π—τ — S : —τ —π—τ —
La succession de transformations subie par le système de spins pendant l’écho s’écrit :

Hevo ”π” Hevo

−−−→ −−−→ −−−→
τ τ
où ”π” est la séquence de refocalisation et Hevo est l’hamiltonien d’évolution libre du système,
supposé faiblement couplé :

Hevo = ΩI Iz + ΩS Sz + πJ2Iz Sz

Cette succession d’évènements est équivalente à

”π” Hred
−−−→ −−−→
2τ

où l’hamiltonien réduit Hred se déduit de Hevo en en conservant que les termes qui sont
invariants par action de la séquence ”π”. Dans le cas A :

πIx,y
Iz −−−−→ −Iz
πIx,y
Sz −−−−→ Sz
πIx,y
2Iz Sz −−−−→ −2Iz Sz

et donc Hred = ΩS Sz . Le résultat de la séquence d’écho de spin ne dépend ni du déplacement

chimique de I (via ΩI ) ni du couplage. Ces interactions sont dites ”refocalisées” au cours de
l’écho de spin. De manière analogue, dans le cas B, Hred = ΩI Iz .
Dans le cas C :
πIx,y πSx,y
Iz −−−−→ −Iz −−−−→ −Iz
πIx,y πSx,y
Sz −−−−→ Sz −−−−→ −Sz
πIx,y πSx,y
2Iz Sz −−−−→ −2Iz Sz −−−−→ 2Iz Sz
et donc Hred = 2Iz Sz . Seule l’interaction de couplage n’est pas refocalisée.

Exercice 10.
La séquence d’impulsion suivante est destinée à enregistrer le spectre de RMN 1 H d’une
substance dissoute dans un mélange H2 O / D2 O.

y
π/2(−x) π/2(−x)
Toutes les impulsions sont
∆ sél. sél. ∆ appliquées exactement à
I
la fréquence de l’eau. Les
τ τ
impulsions sélectives sont
supposées n’affecter que le
Gz signal de l’eau.

9
 Montrez qu’une impulsion de gradient de champ statique d’intensité G, appliquée pendant
le temps τ à un noyau de rapport gyromagnétique γ fait tourner l’aimantation autour de l’axe
Z des noyaux à l’altitude z d’une quantité φ = kz où k = −γGτ .
Analysez l’action de la séquence d’impulsions pour le signal de l’eau et pour les autres
signaux.

Pendant une impulsion de gradient de direction z et à l’altitude z de l’échantillon (z = 0

au centre magnétique de la sonde)

B0 (z) = B0 + Gz

La pulsation de Larmor pour un noyau de rapport gyromagnétique γ dépend de z et de G,

elle vaut
ω(z, G) = −γ(B0 + Gz) = −γB0 − γGz = ω(z, G = 0) − γGz
Tout se passe donc comme si à l’altitude z, la pulsation de Larmor des noyaux est augmentée
de −γGz. Le supplément de champ B0 entraı̂ne pendant le temps τ un supplément d’angle
de précession de l’aimantation transversale :

φ = −γGτ z = kz

que ce soit dans le référentiel fixe ou dans le référentiel tournant. Si au début de l’impulsion de
gradient l’aimantation transversale ne dépend pas de z, elle se trouve distribuée de manière
hélicoı̈dale le long de l’axe Z après le temps τ . Le pas de l’hélice λ correspond à l’espacement
∆z tel que |∆φ| = 2π à savoir
2π
2π = |k|∆z soit λ =
|k|

Une telle distribution spatiale d’aimantation transversale ne produit aucun signal si le nombre
de tours d’hélice est très grand devant 1. G est généralement suffisamment intense pour sa-
tisfaire à cette condition. Ainsi, il n’est possible d’enregistre un signal non nul que si l’aiman-
tation transversale est indépendante de z.
Pour le signal de l’eau, les trois impulsions RF du centre de la séquence ont un effet global
nul. L’offset de l’eau étant nul, l’aimantation transversale créée par le première impulsion RF
est déphasée de kz par la première impulsion et par la seconde impulsion de gradient, soit au
total de φ = 2kz. Le signal de l’eau est ainsi supprimé.
Pour les autres signaux que celui de l’eau, tout se passe comme si les impulsions sélectives
n’existent pas. La séquence est alors celle d’un écho de spin homonucléaire. L’action des
impulsions de gradients, identiques entre elles, est compensée en ce qui concerne les offsets
des noyaux. Il faut toutefois maintenir T assez petit pour que les couplages J des noyaux
n’agissent pas manière significative (JT ≪ 1).

Exercice 11.
La séquence INEPT complète (refocalisée) s’écrit :

ϕ1
Les impulsions étroites et
déc.
larges sont respectivement
I
τ τ τ τ d’angle π/2 et π. Leur
phase vaut x sauf indica-
tion contraire.
S

10
avec τ = 1/4J, où J est la constante de couplage des noyaux I et S.
Quelle doit être la valeur de la phase ϕ pour que le transfert d’aimantation de I vers S
soit le plus efficace possible ?
Quels sont les rôles du premier et du second écho de spin ?
Quel est l’intérêt de cette séquence ?

On considère que l’état initial du système s’écrit

σ0 = Iz + aSz où a = γS /γI

Avec I ≡1 H, a = 1/4 si S ≡13 C et a = −1/10 si S ≡15 N. L’évolution de l’état initial

commence par

π/2Ix
1ère impulsion : Iz −−−−→ −Iy
πIx πSx πJ2I Sz
1er écho : −Iy −−−→ Iy −−−→ Iy −−−−−z−−→ −2Ix Sz
2 × 1/4J

L’état obtenu évoluerait vers un doublet antiphase si le signal était détecté à partir de cet
instant. Un tel état est appelé état couplé d’aimantation transversle du noyau I. Le but des
deux impulsions simultanées d’angle π/2 suivantes est de réaliser un transfert d’aimantation
de I vers S via un état de type ”population” 2Iz Sz . Il faut pour le créer que ϕ = ±y. Avec
ϕ=y :
π/2Iy π/2Sx
transfert : −2Ix Sz −−−−→ 2Iz Sz −−−−−→ −2Iz Sy
un état couplé d’aimantation transversale du noyau S est créé. Il est facile de vérifier que
ϕ = ±x produit un état à zéro et double quanta, qui ne fournira pas d’aimantation détectable
pendant la période d’acquisition du signal.
Finalement :
πIx πSx π/2 · 2Iz Sz
2ième écho : −2Iz Sy −−−→ 2Iz Sy −−−→ −2Iz Sy −−−−−−−−→ Sx .

La séquence de découplage exercée sur les noyaux I pendant l’acquisition du signal provenant
des noyaux S peut être assimilée à une sucession d’écho de spins où seuls les noyaux I
sont soumis aux impulsions d’angle π. En conséquence l’hamiltonien d’évolution pendant
l’acquisition se réduit à Hdec = ΩS Sz . L’état Sx évolue alors pendant l’acquisition pour
donner un singulet d’offset ΩS .
L’aimantation initiale de S, décrite par aSz , subit les transformations :

π/2Ix
1ère impulsion : aSz −−−−→ aSz
πIx πSx π/2 · 2Iz Sz
1er écho : aSz −−−→ aSz −−−→ −aSz −−−−−−−−→ −aSz
π/2Iy π/2Sx
transfert : −aSz −−−−→ −aSz −−−−−→ aSy
πIx πSx π/2 · 2Iz Sz
2ième écho : aSy −−−→ aSy −−−→ −aSy −−−−−−−−→ a · 2Iz Sx .
ΩS tSz
détection : a · 2Iz Sx −−−−−→ rien de détectable

Seule l’aimantation initiale de I conduit à un signal détectable de S, sous forme d’un singulet
d’offset ΩS et d’intensité relative 1. L’expérience ”impulsion–détection” appliquée aux noyaux
S donnerait un signal d’intensité relative a, c’est-à-dire plus faible que celui fourni par la
séquence INEPT. Un transfert d’aimantation par effet Overhauser est toutefois possible et
fournirait une amplification du signal de S par un facteur au plus égal à 1 + 1/2a, soit 3 si
S ≡13 C et -4 si S ≡15 N. La séquence INEPT procure des amplifications par un facteur 4 et
-10, respectivement pour les noyaux 13 C et 15 N.

11
 Le premier écho de spin sert à créer de l’aimantation tranversale couplée de I convertible
en aimantation transversale couplée de S. Le second écho de spin est nécessaire pour obtenir
de l’aimantation transversale de S non couplée, susceptible de fournir un signal pendant la
période d’acquisition. La séquence complète (aimantation en phase I vers aimantation en
phase de S) nécessite une durée de l’ordre de 1/J.

Exercice 12.
La séquence suivante est appliquée à système hétéronucléaire de deux spins I et S couplés.
Pour simplifier, on considèra d’abord que ΩI = 0.
Comment le signal observé varie-t-il en fonction de l’angle α ?
Quelle est l’application pratique de cette séquence ?

π/2(y)

I
1 α(y)
2J
S

Le calcul de l’évolution de l’aimantation initiale de I donne

π/2Iy 0Iz Ω Sz π/2Iz Sz αSy
Iz −−−−→ Ix −−→ Ix −−S−−→ Ix −−−−−−→ 2Iy Sz −−−→ cos α2Iy Sz + sin α2Iy Sx

Lorsque α vaut π/2, il n’y a aucune aimantation détectable pour le noyau I car l’aimantation
initiale a été intégralement convertie en état à 0 et 2 quanta. En augmentant progressivement
α depuis la valeur nulle, l’intensité du doublet antiphase produit par l’évolution de 2Iy Sz
diminue. Il est facile de vérifier que le résultat est inchangé si ΩI 6= 0 et si on tient compte
de l’évolution de l’aimantation initiale du noyau S.
La recherche du doublet d’intensité nulle permet de calibrer une impulsion appliquée
aux noyaux S en détectant le noyau I. Cette situation est celle des expériences 2D dites de
”détection inverse” (HSQC, HMBC).

Exercice 13.
La séquence BIRD

I
1 1
2J 2J

est appliquée à l’aimantation d’équilibre d’un système IS couplé.

Comment est transformée l’aimantation initiale du noyau I ?
Répondre à la même question pour un noyau I isolé soumis à la même séquence.

12
 Quel est l’intérêt de cette séquence ?

Pour un noyau I couplé :

π/2Ix
1ère impulsion : Iz −−−−→ −Iy
πIy πSx π/2 · 2Iz Sz
écho : −Iy −−→ −Iy −−−→ −Iy −−−−−−−−→ Iy
π/2Ix
dernière impulsion : Iy −−−−→ Iz

et pour un noyau I non couplé

π/2Ix
1ère impulsion :
Iz −−−−→ −Iy
πIy
écho : −Iy −−→ −Iy
π/2Ix
dernière impulsion : −Iy −−−−→ −Iz

Dans le cas où I ≡1 H et S ≡13 C, pour un échantillon en abondance naturelle du 13 C, il y

a environ 99 noyaux I isolés pour 1 noyau I couplé. Dans l’expérience HMQC, ce sont ces
derniers qu’on souhaite observer. Les signaux des noyaux isolés sont en théorie éliminables
par le programme de phase mais l’efficacité de cette seule méthode est limitée. La séquence
BIRD laisse l’aimantation d’intérêt à sa position d’équilibre alors que celle de noyaux isolés est
inversée le long de l’axe Z. En attendant le délai nécessaire pour que l’aimantation inversée
passe par zéro (par relaxation longitudinale) avant de commencer la séquence HMQC, on
réalise une suppression substantielle de signaux qui conduisent à la création du bruit ”t1 ”.

Exercice 14.
Un système de spins homonucléaire IS est soumis à une impulsion d’excitation sélective
d’angle π/2 sur le noyau I. Le système est ensuite soumis à un train d’impulsions de re-
focalisation de durée t. On considère que l’opérateur qui agit alors sur le système de spins
s’écrit
Hiso = πJ(2Ix Sx + 2Iy Sy + 2Iz Sz ).
Déterminer l’état du système à la fin de cette séquence de mélange isotrope.
Combien faut-il de temps pour transférer l’aimantation de I vers un maximum d’aiman-
tation en phase de S ?

En considérant que l’impulsion initiale a créé un état Ix , son évolution sous l’action de
Hisopendant le temps t est fournie par le schéma suivant :

Ix
πJt · 2Ix Sx

Ix Les trois termes de Hiso com-

mutent entre eux, ils peuvent
πJt · 2Iy Sy
donc être appliqués succes-
sivement et dans n’importe
Ix −2Iz Sy
quel ordre.
πJt · 2Iz Sz

Ix 2Iy Sz −2Iz Sy Sx

Hiso crée de l’aimantation transversale non couplée de S à partir d’aimantation trans-

versale non couplée de I, en quantité proportionnelle à sin2 (πJt). Cette quantité atteint son

13
premier maximum quand πJt = π/2, c’est-à-dire quand t = 1/2J. A cet instant tous les
termes de l’état du système autres que Sz ont une valeur nulle. La séquence de mélange
isotrope est à la base de la séquence d’implsions 2D TOCSY. Elle est utilisée ici avec une
impulsion sélective d’excitation dans la séquence 1D TOCSY.

Exercice 15.
Ecrire l’opérateur d’évolution libre d’un système de trois spins I, S et L couplés deux à
deux.
Comment évoluent les états Ix , 2Ix Sz et 4Ix Sz Lz pendant le temps t ?
Quel est l’opérateur réduit correspondant à une séquence d’écho de spin où seule l’aiman-
tation du noyau I est refocalisée ?

L’opérateur hamiltonien d’évolution libre d’un système de trois spins I, S et L faiblement

couplés s’écrit

Hevo = ΩI Iz + ΩS Sz + ΩL Lz + πJIS 2Iz Sz + πJIL 2Iz Lz + πJSL 2Sz Lz

Les états initiaux Ix , 2Ix Sz et 4Ix Sz Lz correspondent à de l’aimantation transversale du

noyau I, respectivement non couplée, couplée avec S, et couplée avec S et L. Chacun de ces
trois états commute avec les opérateurs d’évolution ΩS Sz , ΩL Lz et πJSL 2Sz Lz . L’évolution
libre de Iz est déterminée par :

Ix
πJIS t · 2Iz Sz

Ix
JIL
πJIL t · 2Iz Lz
JIS
Ix
Ω I t · Iz
ΩI
Ix Iy

Seules les branches de l’arbre qui contribuent au signal sont représentées. Le signal enre-
gistré est ainsi déterminé par

s(t) = cos(πJIS t) cos(πJIL t) exp(iΩI t)

1
= (exp(iπJIS t) + exp(−iπJIS t)) (exp(iπJIL t) + exp(−iπJIL t)) exp(ΩI t)
4
1 1
= exp (i (ΩI + πJIS + πJIL ) t) + exp (i (ΩI + πJIS − πJIL ) t)
4 4
1 1
+ exp (i (ΩI − πJIS + πJIL ) t) + exp (i (ΩI − πJIS − πJIL ) t)
4 4

et sa TF fourni quatre signaux d’intensités relatives 41 aux fréquences νI ± JIS /2 ± JIL /2.
Le même type de calcul, mené pour un état initial 2Ix Sz conduit à

14
 2Ix Sz
JIL
πJIS t · 2Iz Sz
JIS
Iy
πJIL t · 2Iz Lz
ΩI
Iy
Ω I t · Iz

Iy −Ix

Ainsi :
s(t) = i sin(πJIS t) cos(πJIL t) exp(iΩI t)
1
= (exp(iπJIS t) − exp(−iπJIS t)) (exp(iπJIL t) + exp(−iπJIL t)) exp(ΩI t)
4
1 1
= exp (i (ΩI + πJIS + πJIL ) t) + exp (i (ΩI + πJIS − πJIL ) t)
4 4
1 1
− exp (i (ΩI − πJIS + πJIL ) t) − exp (i (ΩI − πJIS − πJIL ) t)
4 4
dont la TF fournit quatre raies d’intensités relatives 41 , − 41 , 14 et − 14 aux mêmes fréquences que
précédemment. La constante de couplage JIS est ici dite ”active” car deux raies distantes de
JIS Hz sont de signe opposés. La constante de couplage JIL est en conséquence dite ”passive”.
Ces termes apparaissent dans l’analyse des spectres COSY des systèmes a 3 spins.
L’évolution du terme 4Ix Sz Lz a lieu comme suit :

4Ix Sz Lz
JIL
πJIS t · 2Iz Sz
JIS
2Iy Lz
πJIL t · 2Iz Lz
ΩI
−Ix
Ω I t · Iz

−Ix −Iy

Ainsi :
s(t) = − sin(πJIS t) sin(πJIL t) exp(iΩI t).
Le développement en exponentielles complexes des deux fonctions sinus fait apparaı̂tre un
doublet de doublet de pics doublement antiphase, d’intensités relatives 41 , − 41 , − 41 et 14 où les
deux constantes de couplage du noyau I sont actives. Chaque opérateur de couplage intervient
par la branche droite (sinus) de l’arbre puisqu’il y a nécessité de faire disparaı̂tre le ”Sz ” et
le ”Lz ” de 4Ix Sz Lz pour obtenir un terme Ix,y détectable.
Les seuls termes de Hevo qui sont invariants par action de πIx,y sont Sz , Lz , 2Sz Lz . En
conséquence l’hamiltonien réduit Hred qui correspond à un écho de spin où seule l’aimantation
de I est inversée est donné par

Hred = ΩS Sz + ΩL Lz + πJSL 2Sz Lz

15
Exercice 16.
La séquence d’impulsions de l’expérience ”J-résolue”

sat. déc.
I
1 1
J J

est appliquée à un système de spins In S (n = 1, 2 ou 3) où les n spins I sont équivalents et

sont couplés à S avec la constante J.
Ecrire la succession des transformations qu’il faut appliquer à l’état initial du système
pour trouver son état juste avant l’acquisition.
Simplifier cette liste en utilisant le concept d’opérateur réduit.
En déduire pourquoi le signe des pics spectraux dépend de la parité de n.

On considère ici un système de trois noyaux (1 H) de type ”I” notés I, J et K et d’un noyau
S (13 C), le tout représentant par exemple un groupe méthyle. Les noyaux I, J et K couplent
identiquement au noyau S avec la constante J. Le champ RF appliqué aux 1 H pendant la
période de relaxation a pour but d’égaliser les populations des états mI,J,K s = ±1/2 et de
produire un acroissement de l’aimantation initiale de S par effet Overhauser. La première
impulsion crée un état Sx . En ne gardant de l’hamiltonien d’évolution libre que les termes
qui ne commutent pas avec Sx , la succession des évènement qu’il faut lui appliquer est la
suivante :
π2Iz Sz π2Jz Sz π2Kz Sz ΩS Sz πSx ΩS Sz
−−−− −→−−−−−→−−−−−−→−−−−→−−−→−−−−→
1/J 1/J
Sachant que
Ω Sz πSx ΩS Sz πSx
−−S−−→ −−−→−−−−→ est équivalent à −−−→
1/J 1/J
puisque l’action de πSx ne laisse pas Sz invariant (règle de fabrication de Hred pour l’écho de
spin), le calcul à effectuer est :
π2Iz Sz π2Jz Sz π2Kz Sz πSx
Sx −−−−−→−−−−−→−−−−−−→−−−→
Chacun des trois premiers opérateurs ne couplant pas avec Sx , ce calcul fournit le résultat
π2Iz Sz π2Jz Sz π2Kz Sz πSx
Sx −−−−−→ −Sx −−−−−→ Sx −−−−− −→ −Sx −−−→ −Sx
Si S couple avec un nombre n impair (ou pair) de noyaux I, il y a un nombre impair
(ou pair) de signes négatifs appliqués à l’état initial, d’où la modulation de signe des pics
spectraux observés en fonction de (−1)n et donc en fonction de la parité de n.

Voir aussi : www.univ-reims.fr/LSD/JmnSoft/livre.pdf

16