format: text

text:
# Méthodes de Kashiwara-Vergne-Rouvière pour certains espaces symétriques

## Introduction

Dans notre article précédent \[To1\], suivant les idées élaborées par Kontsevich \[Ko\] pour la quantification formelle des variétés de Poisson et déjà utilisées dans \[ADS\] et \[AST\], nous avons montré comment résoudre certaines conjectures sur la formule de Campbell-Hausdorff (BCH in English). Plus précisément, nous avons construit une déformation de la formule BCH, qui vérifie une équation différentielle analogue à celle que l’on trouve dans l’article de Kashiwara-Vergne \[KV\]. Ces équations ont la même utilité que celles de l’article de Kashiwara-Vergne. En effet elles permettent de montrer, comme dans \[ADS\] et \[AST\], que l’application exponentielle modifiée par la racine carrée du jacobien, transporte la convolution des distributions invariantes sur les algèbres de Lie (\[Du\],\[KV\], \[Ve\],\[Mo\]).

Dans cet article nous considérons le cas des espaces symétriques.

Dans notre mémoire d’Habilitation \[To2\] nous avons par ailleurs abordé certains problèmes d’analyse harmonique pour les espaces symétriques. Nous avons notamment expliqué en quoi le comportement des espaces symétriques est à la fois proche et éloigné du cas des algèbres de Lie. On peut dire en résumant grossièrement, qu’il existe une catégorie d’espaces symétriques pour laquelle les méthodes des groupes fonctionnent bien (cas nilpotent, résoluble, $`\mathrm{"}G_{\mathrm{I}\mathrm{C}}/G_{\mathrm{I}\mathrm{R}}\mathrm{"}`$, cas ”très symétrique”). Malheureusement le cas général des espaces symétriques ne peut pas être traité par des méthodes s’inspirant simplement du cas des groupes (voir la fin de cet article pour quelques remarques).

Dans cet article nous allons déformer la formule de Campbell-Hausdorff pour les espaces symétriques généraux.

Nous montrerons que cette déformation vérifie, dans les ”bons” cas évoqués plus haut (résoluble ou très-symétrique), des équations analogues au cas des groupes. On en déduira que l’application exponentielle modifiée par la racine carrée du déterminant jacobien, transporte la convolution des germes de distributions invariantes. Nous retrouvons ainsi les résultats de Rouvière \[Rou1\] du cas résoluble. Notre démonstration fournit une preuve dans le cas très symétrique ce qui, à notre connaissance, est un résultat nouveau.

Remerciements: L’auteur remercie M. Pevzner et le département de mathématiques l’Université Libre de Bruxelles (ULB) pour leur invitation (fin janvier 2002). L’auteur y a finalisé les résultats de cet article. Nous remercions aussi les organisateurs du colloque Carmona pour leur invitation. Nous avons eu le plaisir d’y présenter les résultats de cet article.

## 1 Rappels

### 1.1 Rappels sur la formule de Kontsevich

La formule de Kontsevich associe à toute structure de Poisson régulière sur $`\mathrm{I}\mathrm{R}^d`$ un étoile-produit formel associatif. Ce n’est qu’un cas particulier du théorème de formalité démontré par Kontsevich dans \[Ko\] théorème $`6.4`$.

Lorsque $`f`$ et $`g`$ sont deux fonctions régulières sur $`\mathrm{I}\mathrm{R}^d`$ et $`\alpha `$ un deux-vecteur de Poisson régulier sur $`\mathrm{I}\mathrm{R}^d`$, Kontsevich écrit dans \[Ko\] paragraphe $`2`$, la formule suivante :

$$fg=fg+\underset{n=1}{\overset{\mathrm{}}{}}\frac{h^n}{n!}\underset{\begin{array}{c}\mathrm{\Gamma }G_{n\mathrm{,2}}\\ \mathrm{\Gamma }\mathrm{a}dmissible\end{array}}{}w_\mathrm{\Gamma }B_\mathrm{\Gamma }(f,g).$$
(1)
Cette formule munit l’espace des fonctions régulières sur $`\mathrm{I}\mathrm{R}^d`$ d’une structure associative formelle. Dans cette formule $`h`$ est un paramètre formel, $`G_{n\mathrm{,2}}`$ désigne l’ensemble des graphes étiquetés avec $`n`$ points de première espèce et $`2`$ points de seconde espèce, $`\mathrm{\Gamma }`$ est un graphe dit admissible parmi les graphes de $`G_{n\mathrm{,2}}`$, $`w_\mathrm{\Gamma }`$ est un cœfficient obtenu par intégration sur un espace de configurations d’une forme différentielle dépendant de $`\mathrm{\Gamma }`$ et $`B_\mathrm{\Gamma }`$ est un opérateur bidifférentiel construit à partir de $`\mathrm{\Gamma }`$.

Par souci de clarté et d’autonomie du texte nous allons préciser maintenant chacun des termes de cette formule et appliquer dans un second temps cette formule dans le cas des structures de Poisson linéaires c’est-à-dire dans le cas du dual des algèbres de Lie.

#### 1.1.1 Espaces de configurations

On note par $`\mathrm{Conf}_{n,m}`$ l’espace des configurations de $`n`$ points distincts dans le demi-plan de Poincaré (ce sont les points aériens) et de $`m`$ points distincts sur la droite réelle (ce sont les points terrestres). Le groupe:

$$G^{(1)}=\{zaz+b\text{ avec }(a,b)\mathrm{I}\mathrm{R}_+^{}\times \mathrm{I}\mathrm{R}\}$$

agit librement sur $`\mathrm{Conf}_{n,m}`$. Le quotient :

$$C_{n,m}=\mathrm{Conf}_{n,m}/G^{(1)}$$

est une variété de dimension $`2n+m2`$. Compte tenu de l’action de ce groupe sur les points terrestres, on peut identifier deux des points terrestres aux points $`0`$ et $`1`$ (à condition que l’on ait $`m2`$, sinon on peut identifier un des points aériens au complexe $`i`$). Dans \[Ko\] paragraphe $`5.1`$, Kontsevich construit des compactifications de ces variétés notées $`\overline{C}_{n,m}`$. Ce sont des variétés à coins de dimension $`2n+m2`$. Ces variétés ne sont pas connexes pour $`m2`$. On notera par $`\overline{C}_{n,m}^+`$ la composante connexe qui contient les configurations où les points terrestres sont ordonnés dans l’ordre croissant (i.e. on a $`\overline{1}<\overline{2}<\mathrm{}<\overline{m}`$). On introduit de manière analogue des variétés de configurations de $`n`$ points dans le plan complexe modulo l’action du groupe $`G^{(2)}=\{zaz+b\text{ avec }(a,b)\mathrm{I}\mathrm{R}_+^{}\times \mathrm{I}\mathrm{C}\}`$. On les note $`C_n`$, ce sont des variétés de dimension $`2n3`$. On note $`\overline{C_n}`$ les compactifications associées (\[Ko\] paragraphe $`5.1`$).

Ces variétés sont stratifiées et chaque strate est décrite par un arbre. En termes géométriques, les strates sont obtenues par concentrations itérées de points en des amas (\[Ko\] paragraphe 5.2).

#### 1.1.2 Graphes admissibles

La notion de graphes admissibles est maintenant bien établie dans la littérature. On désigne par $`G_{n,m}`$ l’ensemble des graphes étiquetés et orientés (les arêtes sont orientées) ayant $`n`$ sommets du première espèce (sommets aériens) et $`m`$ sommets du deuxième espèce (sommets terrestres). Par graphe étiqueté on entend un graphe $`\mathrm{\Gamma }`$ muni d’un ordre total sur l’ensemble $`E_\mathrm{\Gamma }`$ de ses arêtes, compatible avec l’ordre des sommets.

Les graphes qui vont intervenir dans la formule de Kontsevich sont dans $`G_{n\mathrm{,2}}`$ et vérifie des conditions supplémentaires. On dira qu’ils sont admissibles (\[Ko\] paragraphe $`6.1`$ et paragraphe $`2`$) si

1. Les arêtes partent toutes des sommets de première espèce.
2. Il part deux arêtes de chaque sommet de première espèce.
3. Le but d’une arête est différent de sa source (il n’y a pas de boucle).
4. Il n’y a pas d’arêtes multiples (même source, même but).

#### 1.1.3 Opérateur différentiel associé à un graphe

Soit $`\mathrm{\Gamma }`$ un graphe admissible de $`G_{n\mathrm{,2}}`$. De tout point aérien $`i`$ partent deux arêtes numérotées $`(e_i^a,e_i^b)`$.

Soit $`\alpha `$ un deux-vecteur sur $`\mathrm{I}\mathrm{R}^d`$. On peut alors associer à tout graphe admissible un opérateur bi-différentiel sur $`\mathrm{I}\mathrm{R}^d`$ comme expliqué dans \[Ko\] paragraphe $`2`$. On notera $`B_\mathrm{\Gamma }(f,g)`$ l’opérateur bidifférentiel associé que l’on suppose agir sur les fonctions $`f`$ et $`g`$.

Expliquons heuristiquement la formule. Sur chaque sommet aérien on met le deux-vecteur et sur les sommets terrestres on met les fonctions $`f`$ et $`g`$. Chaque arête arrivant sur un sommet dérive la fonction associée au sommet. On multiplie les fonctions ainsi obtenues et on somme sur toutes les possibilités. Concrètement la formule est la suivante. Pour chaque arête $`e`$, on note par $`s(e)`$ le point aérien source (départ) et par $`b(e)`$ le point but (arrivée). Dans la formule ci-dessous $`I`$ décrit l’ensemble les applications de l’ensemble des arêtes $`E_\mathrm{\Gamma }`$ dans l’ensemble des indices de coordonnées $`\{1,\mathrm{}d\}`$.

$$B_{\mathrm{\Gamma },\alpha }(f,g)=\underset{I}{}\left[\underset{k=1}{\overset{n}{}}\left(\underset{\begin{array}{c}eE_\mathrm{\Gamma }\\ b(e)=k\end{array}}{}_{I(e)}\right)\alpha ^{I(e_k^a)I(e_k^b)}\right]\left(\underset{\begin{array}{c}eE_\mathrm{\Gamma }\\ b(e)=\overline{1}\end{array}}{}_{I(e)}\right)f\left(\underset{\begin{array}{c}eE_\mathrm{\Gamma }\\ b(e)=\overline{2}\end{array}}{}_{I(e)}\right)g.$$

#### 1.1.4 Forme d’angles

Soient deux points distincts $`(p,q)`$ dans le demi-plan de Poincaré muni de la métrique de Lobachevsky. On note

$$\varphi (p,q)=\frac{1}{2i}\mathrm{log}(\frac{(qp)(\overline{q}p)}{(q\overline{p})(\overline{q}\overline{p})}).$$
(2)
C’est l’angle entre la géodésique $`(p,\mathrm{})`$ et $`(p,q)`$ où l’infini peut être vu comme l’infini sur la droite réelle (figure )

La fonction d’angle s’étend à la compactification $`\overline{C}_{\mathrm{2,0}}`$ en une fonction d’angle régulière. La variété $`\overline{C}_{\mathrm{2,0}}`$ est précisément décrite dans l’article de Kontsevich (\[Ko\] paragraphe $`5.2`$), c’est le fameux œil (figure ). On remarquera, mais c’est tautologique vu la construction des compactifications, que lorsque les points $`p,q`$ s’approchent selon un angle $`\theta `$, la fonction d’angle vaut précisément cet angle. Lorsque $`p`$ s’approche de l’axe réel la fonction d’angle est nulle et lorsque $`q`$ s’approche de l’axe réel on obtient deux fois l’angle de demi-droite avec l’axe réel.

Comme la fonction d’angle est régulière sur la compactification, on peut considérer sa différentielle qui est alors une $`1`$-forme sur $`\overline{C}_{\mathrm{2,0}}`$.

#### 1.1.5 Poids associé à un graphe

Si $`\mathrm{\Gamma }`$ est un graphe admissible dans $`G_{n\mathrm{,2}}`$, alors toute arête $`e`$ définit par restriction une fonction d’angle notée $`\varphi _e`$ sur la variété $`\overline{C}_{n\mathrm{,2}}^+`$. Le produit ordonné

$$\mathrm{\Omega }_\mathrm{\Gamma }=\underset{eE_\mathrm{\Gamma }}{}d\varphi _e$$
(3)
est donc une $`2n`$-forme sur $`\overline{C}_{n\mathrm{,2}}^+`$ qui est de dimension aussi $`2n`$. Le poids associé sera donc

$$w_\mathrm{\Gamma }=\frac{1}{(2\pi )^{2n}}_{\overline{C}_{n\mathrm{,2}}^+}\mathrm{\Omega }_\mathrm{\Gamma }.$$
(4)
#### 1.1.6 Permutation des arêtes

Soit $`\mathrm{\Gamma }`$ un graphe admissible dans $`G_{n\mathrm{,2}}`$ et $`\alpha `$ un deux-vecteur régulier sur $`\mathrm{I}\mathrm{R}^d`$. Le groupe $`\underset{n}{\underset{}{S_2\times \mathrm{}\times S_2}}`$, produit des groupes de permutations des arêtes attachés à chaque sommet, agit naturellement sur $`\mathrm{\Gamma }`$ par permutation de l’étiquetage des arêtes. On vérifie que l’on a ($`ϵ(\sigma )`$ désigne la signature de la permutation $`\sigma `$):

$`B_{\sigma .\mathrm{\Gamma }}`$ $`=\epsilon (\sigma )B_\mathrm{\Gamma }`$
$`w_{\sigma .\mathrm{\Gamma }}`$ $`=\epsilon (\sigma )w_\mathrm{\Gamma },`$ (5)
de sorte que le produit $`w_\mathrm{\Gamma }.B_\mathrm{\Gamma }`$ ne dépend pas de l’étiquetage. On remarque aussi que le cœfficient $`w_\mathrm{\Gamma }`$ et l’opérateur $`B_\mathrm{\Gamma }`$ ne dépendent pas de l’ordre des sommets.

#### 1.1.7 Principe de la démonstration

L’associativité de l’étoile-produit défini par Kontsevich résulte de deux ingrédients. Le premier ingrédient est la formule de Stokes qui va donner des relations entre les cœfficients $`w_\mathrm{\Gamma }`$. Le deuxième ingrédient est un lemme fondamental (\[Ko\] lemme $`6.6`$) qui assure que les concentrations aériennes de plus de trois points auront des contributions nulles. Au final le défaut d’associativité de cet étoile-produit sera nul car le deux-vecteur $`\alpha `$ est un deux-vecteur de Poisson (\[Ko\] paragraphe $`6.4`$).

### 1.2 Le cas linéaire

Soit $`𝔤`$ une algèbre de Lie de dimension finie sur $`\mathrm{I}\mathrm{R}`$. L’espace vectoriel dual $`𝔤^{}`$ possède une structure de Poisson linéaire donnée par la moitié du crochet de Lie. Ces deux notions sont équivalentes. Si $`e_1,\mathrm{},e_d`$ désigne une base de $`𝔤`$, $`e_1^{},\mathrm{},e_d^{}`$ la base duale et $`_i`$ la dérivée dans la direction du vecteur $`e_i^{}`$, alors le deux-vecteur de Poisson associé est $`\alpha =\frac{1}{2}_{i,j}[e_i,e_j]_i_j`$.

Dans ce cas (\[Ko\] théorème $`\mathrm{8.3.1}`$) l’ordre des opérateurs $`B_\mathrm{\Gamma }`$ est suffisamment croissant de sorte que, lorsque $`f`$ et $`g`$ sont deux fonctions polynomiales, la formule de Kontsevich est en fait une somme finie. Cet étoile-produit vérifie pour $`X,Y`$ dans $`𝔤`$ (considérés comme des fonctions linéaires sur $`𝔤^{}`$) la relation

$$XYYX=h[X,Y].$$

En localisant en $`h=1`$, l’étoile-produit de Kontsevich définit donc sur l’algèbre symétrique $`S[𝔤]`$ une structure isomorphe à l’algèbre enveloppante de $`𝔤`$ (\[Ko\] théorème $`\mathrm{8.3.1}`$).

#### 1.2.1 Géométrie des graphes

Dans le cas linéaire, les graphes qui vont intervenir de manière non triviale dans la formule (1) ont une description relativement simple. En effet chaque sommet de première espèce ne pourra recevoir qu’au plus une arête: ce sont par définition les graphes pertinents (relevant graphs) et on renvoie le lecteur à l’article \[AST\] paragraphe $`3.1`$ pour une description précise des graphes pertinents qui interviennent dans la formule finale.

En fait des arguments élémentaires (détaillés dans \[AST\] paragraphe $`3.1`$) montrent que tout graphe admissible pertinent se décompose en produit de graphes simples. Il y a deux types de graphes simples: ceux qui contiennent une seule roue (et donc pas de racine) comme dans la figure  on les appellera graphes de type roue  (\[AST\] définition $`\mathrm{3.1.1}`$) et ceux qui contiennent une seule racine (et donc pas de roue) comme dans la figure , ce sont les graphes de type Lie (\[Ka\] définition $`3.1`$).

#### 1.2.2 Symbole de $`B_\mathrm{\Gamma }`$

Dans \[AST\] paragraphe $`3`$ on a associé à chaque graphe admissible $`\mathrm{\Gamma }`$ un symbole $`a_\mathrm{\Gamma }`$. C’est une fonction de $`𝔤\times 𝔤`$ à valeurs dans $`S[𝔤]`$ donnée par la formule pour $`X,Y`$ dans $`𝔤`$

$$a_\mathrm{\Gamma }(X,Y)=B_\mathrm{\Gamma }(e^X,e^Y)e^{XY}.$$

Par exemple le symbole associé au graphe de la figure  est

$$a_\mathrm{\Gamma }(X,Y)=\frac{1}{8}[[X,[X,Y]],Y]$$

et le symbole associé au graphe de la figure  est

$$a_\mathrm{\Gamma }(x,y)=\frac{1}{2^5}\mathrm{tr}_𝔤(\mathrm{ad}[X,Y]\mathrm{ad}X\mathrm{ad}Y\mathrm{ad}Y).$$

Lorsque le graphe est simple de type Lie, alors $`a_\mathrm{\Gamma }(X,Y)`$ est naturellement un élément de l’algèbre de Lie engendrée par $`X`$ et $`Y`$. Lorsque $`\mathrm{\Gamma }`$ est produit de graphes simples, le symbole est le produit des symboles associés (\[AST\] lemme $`3.6`$).

#### 1.2.3 Formule de Campbell-Hausdorff en termes de graphes

Comme démontré dans \[Ka\] thérème $`5.1`$, la formule de Campbell-Hausdorff s’écrit alors pour $`X`$ et $`Y`$ dans $`𝔤`$ (la série est convergente pour $`X`$ et $`Y`$ proches de $`0`$)

$$Z(X,Y)=X+Y+\underset{n=1}{\overset{\mathrm{}}{}}\frac{1}{n!}\underset{\mathrm{\Gamma }}{}w_\mathrm{\Gamma }a_\mathrm{\Gamma }(X,Y).$$
(6)
où la somme porte sur les graphes simples de type Lie. Les graphes qui contribuent de manière non triviale dans cette formule n’ont donc qu’une racine et ne possèdent pas de symétries. Par conséquent les graphes étiquetés (numérotés) associés à un graphe géométrique (graphe orienté associé pour lequel on oublie l’étiquetage) sont au nombre de $`n!2^n`$. Le terme $`n!`$ se compense avec le terme $`\frac{1}{n!}`$ et le terme $`2^n`$ disparaît aussi car on a pris le deux-vecteur de Poisson associé à la moitié du crochet de Lie. On notera $`\mathrm{\Gamma }(X,Y)`$ le symbole associé au graphe géométrique lorsqu’on prend le deux-vecteur de Poisson associé au crochet de Lie. On a alors la formule plus synthétique suivante

$$Z(X,Y)=X+Y+\underset{\mathrm{\Gamma }}{}w_\mathrm{\Gamma }\mathrm{\Gamma }(X,Y)$$
(7)
où la somme porte sur les graphes géométriques simples de type Lie. On remarquera toutefois que le symbole $`\mathrm{\Gamma }(X,Y)`$ est mal défini si le graphe n’est pas étiqueté. Pour résoudre ce problème il suffit de remarquer que c’est aussi le cas pour le cœfficient $`w_\mathrm{\Gamma }`$ et que ces deux difficultés se compensent grâce aux équations (5).

#### 1.2.4 Symbole formel de l’étoile-produit dans le cas linéaire et fonction de densité

La propriété multiplicative du symbole et la combinatoire sur les graphes assurent que le symbole formel associé à l’étoile-produit de Kontsevich est, dans le cas des algèbres de Lie, un produit de deux termes de type exponentiel correspondant aux contributions des graphes simples de type Lie (c’est la formule de Campbell-Hausdorff par \[Ka\] théorème $`5.1`$) et des graphes simples de type roue (c’est la fonction de densité par \[AST\] proposition $`3.7`$). Le symbole formel vu dans $`S[𝔤][[h]]`$, vérifie alors (\[AST\] proposition $`3.8`$)

$$\left(\mathrm{exp}(X)\mathrm{exp}(Y)\right)\mathrm{exp}(XY)=D(hX,hY)\mathrm{exp}\left(\frac{1}{h}Z(hX,hY)XY\right)$$
(8)
avec $`D(X,Y)`$ la fonction de densité (série convergente pour $`X`$ et $`Y`$ proches de $`0`$) :

$$D(X,Y)=\mathrm{exp}\left(\underset{\mathrm{\Gamma }}{}w_\mathrm{\Gamma }\mathrm{\Gamma }(X,Y)\right)$$
(9)
où la somme porte sur l’ensemble des graphes géométriques simples et de type roue. On montre (\[AST\] proposition $`3.12`$) que cette fonction vaut

$$D(X,Y)=\frac{j^{1/2}(X)j^{1/2}(Y)}{j^{1/2}(Z(X,Y))}$$
(10)
avec $`j`$ le déterminant jacobien de l’application exponentielle à savoir

$$j(X)=\underset{𝔤}{det}\left(\frac{1\mathrm{exp}(\mathrm{ad}X)}{\mathrm{ad}X}\right).$$
(11)
### 1.3 Rappel sur les espaces symétriques

Soit $`𝔤`$ une algèbre de Lie de dimension finie sur $`\mathrm{I}\mathrm{R}`$, on appelle paire symétrique, tout couple $`(𝔤,\sigma )`$ avec $`\sigma `$ une automorphisme de Lie involutif. Habituellement on note $`𝔤=𝔨𝔭`$ la décomposition en espaces propres relativement à $`\sigma `$.

Lorsque $`𝔤`$ est l’algèbre de Lie d’un groupe de Lie $`G`$ (connexe et simplement connexe) et $`\sigma `$ la différentielle d’un involution de $`G`$ (un automorphisme de groupe vérifiant $`\sigma ^2=1`$), on note $`G^\sigma `$ le sous-groupe des points fixes de $`\sigma `$ dans $`G`$ et $`K`$ la composante connexe de $`G^\sigma `$ qui contient l’élément neutre $`e`$. L’espace $`G/K`$ est appelé espace symétrique. Alors $`𝔨`$ est l’algèbre de Lie de $`K`$, mais en général $`𝔭`$ n’est pas une sous-algèbre de Lie car on a les relations suivantes $`[𝔨,𝔭]𝔭`$ et $`[𝔭,𝔭]𝔨`$. On peut identifier $`𝔭`$ à l’espace tangent $`T_{eK}(G/K)`$. Il existe cependant sur $`𝔭`$ un produit triple donné par la formule pour $`X,Y,Z`$ dans $`𝔭`$

$$(X,Y,Z)[[X,Y],Z].$$
(12)
Les espaces vectoriels munis de produits triples sont pour les espaces symétriques, ce que sont les algèbres de Lie pour les groupes de Lie (\[Lo\] vol I, paragraphe II-$`2`$).

L’application exponentielle définit un difféomorphisme local de $`𝔭`$ sur $`G/K`$. On la note $`\mathrm{Exp}`$ et on notera $`\mathrm{exp}_𝔤`$ l’application exponentielle de $`𝔤`$ dans $`G`$. On en déduit alors l’existence d’une formule de Campbell-Hausdorff pour les espaces symétriques, définie de la manière suivante. Pour $`X`$ et $`Y`$ dans $`𝔭`$ proches de $`0`$, il existe une série convergente $`Z_{sym}(X,Y)`$ à valeurs dans $`𝔭`$ telle que l’on ait

$$\mathrm{exp}_𝔤(X)\mathrm{Exp}(Y)=\mathrm{Exp}\left(Z_{sym}(X,Y)\right).$$

En utilisant l’involution $`\sigma `$ on trouve facilement

$$\mathrm{exp}_𝔤(2Z_{sym}(X,Y))=\mathrm{exp}_𝔤(X)\mathrm{exp}_𝔤(2Y)\mathrm{exp}_𝔤(X).$$
(13)
Cela nous suggère de modifier le crochet de Lie pour $`𝔤`$ et de prendre plutôt deux fois le crochet. Notons $`𝔤_2`$ l’algèbre de Lie ainsi obtenue et $`G_2`$ le groupe de Lie connexe et simplement connexe associé. Notons $`\mathrm{exp}_{𝔤_2}`$ l’application exponentielle associée. On aura alors (le membre de droite est calculé dans $`G_2`$)

$$\mathrm{exp}_{𝔤_2}\left(Z_{sym}(X,Y)\right)=\mathrm{exp}_{𝔤_2}(X/2)\mathrm{exp}_{𝔤_2}(Y)\mathrm{exp}_{𝔤_2}(X/2).$$
(14)
## 2 Méthodes de déformation pour les espaces symétriques

### 2.1 Calcul de la formule de Campbell-Hausdorff

Comment exprimer maintenant $`Z_{sym}(X,Y)`$ en termes de diagrammes? Cela revient à se demander comment on calcule les produits itérés de l’étoile-produit. La proposition suivante résout la question (voir l’article \[Ko\] paragraphe $`\mathrm{8.3.3.2}`$ pour un calcul analogue).

On considère $`G_{n\mathrm{,3}}`$ l’ensemble des graphes étiquetés avec $`n`$ points aériens et $`3`$ points terrestres. On place les points terrestres en $`0,s\mathrm{,1}`$ avec $`s]\mathrm{0,1}[`$ (figure ). On étend la notion de graphes admissibles aux graphes de $`G_{n\mathrm{,3}}`$ (voir section 1.1.2).

Notons $`w_\mathrm{\Gamma }(s)`$ la valeur du cœfficient obtenu par intégration de la forme $`\mathrm{\Omega }_\mathrm{\Gamma }`$ sur les configurations à $`s`$ fixe (en d’autres termes on intègre sur les points aériens).

On considère la structure de Poisson linéaire sur $`𝔤^{}`$, donnée par la moitié du crochet de Lie. On note $`\alpha `$ le deux-vecteur associé.

###### Proposition 1.

Pour $`f,g,h`$ trois fonctions régulières sur $`𝔤^{}`$ on a

$$fgk=fgk+\underset{n=1}{\overset{\mathrm{}}{}}\frac{h^n}{n!}\underset{\begin{array}{c}\mathrm{\Gamma }G_{n\mathrm{,3}}\\ \mathrm{\Gamma }\mathrm{admissible}\end{array}}{}w_\mathrm{\Gamma }(s)B_\mathrm{\Gamma }(f,g,k).$$

Preuve: On doit montrer deux choses, d’une part que l’expression de droite est indépendante de $`s`$ et qu’elle vaut bien le terme de gauche.

On utilise le calcul de la dérivée de $`w_\mathrm{\Gamma }(s)`$ comme dans \[To1\] paragraphe $`4.1`$. Cela revient à utiliser la formule de Stokes de manière infinitésimale qui est à la base de l’associativité (\[Ko\] paragraphe $`6.4`$). En pratique cela veut dire que l’on doit concentrer des points. Comme les points $`0,s\mathrm{,1}`$ sont sur l’axe réel, les concentrations que l’on doit effectuer sont aériennes (sinon pour des raisons de dimension on se retrouve dans le cas \[Ko\] paragraphe $`\mathrm{6.4.2.2}`$). Les concentrations de plus de trois points aériens sont nulles (\[Ko\] lemme $`6.6`$). Les seules expressions non nulles, proviennent alors des concentrations de deux points aériens. Elles se compensent par l’identité de Jacobi, à savoir $`[\alpha ,\alpha ]_S=0`$ (crochet de Schouten). Le terme de droite est donc constant.

Maintenant lorsque $`s`$ tend vers $`0`$, on obtient toutes les concentrations de $`p`$ points aériens avec leurs $`2p`$ arêtes. En effet, il faut tenir compte des composantes de bord dans les compactifications des espaces de configurations (\[Ko\] $`\mathrm{6.4.2.1}`$). Le résultat vaut alors la factorisation

$$(fg)k.$$

On en retrouve l’associativité de l’étoile-produit en considérant la limite quand $`s`$ tend vers $`1`$.$`\mathrm{}`$

Remarque :  Dans la proposition ci-dessus on n’a pas utilisé le fait que la structure de Poisson était linéaire, on en déduit que cette proposition est vraie pour tout deux-vecteur de Poisson.

Plaçons maintenant le point médian en $`s=\frac{1}{2}`$. On place en $`0`$ la fonction $`\mathrm{exp}(X/2)`$, en $`\frac{1}{2}`$ la fonction $`\mathrm{exp}(Y)`$ et en $`1`$ la fonction $`\mathrm{exp}(X/2)`$. On s’intéresse à la contribution des termes dans

$$\mathrm{exp}(X/2)\mathrm{exp}(Y)\mathrm{exp}(X/2)=$$

$$\mathrm{exp}(X+Y)+\underset{n=1}{\overset{\mathrm{}}{}}\frac{h^n}{n!}\underset{\begin{array}{c}\mathrm{\Gamma }G_{n\mathrm{,3}}\\ \mathrm{\Gamma }\mathrm{admissible}\end{array}}{}w_\mathrm{\Gamma }(\frac{1}{2})B_\mathrm{\Gamma }(\mathrm{exp}(X/2),\mathrm{exp}(Y),\mathrm{exp}(X/2)).$$
(15)
Le bivecteur de Poisson intervenant est comme d’habitude la moitié du crochet de Lie (n’oublions pas que l’on travaille dans $`𝔤_2`$). On fait maintenant intervenir la symétrie par rapport à la droite $`x=\frac{1}{2}`$. Si $`\mathrm{\Gamma }`$ est un graphe, on note $`\mathrm{\Gamma }^{}`$ le graphe symétrisé. Clairement, les quantités

$$B_\mathrm{\Gamma }(\mathrm{exp}(X/2),\mathrm{exp}(Y),\mathrm{exp}(X/2))$$

et

$$B_\mathrm{\Gamma }^{}(\mathrm{exp}(X/2),\mathrm{exp}(Y),\mathrm{exp}(X/2))$$

sont égales. Par contre on a $`w_\mathrm{\Gamma }^{}=(1)^nw_\mathrm{\Gamma }`$. On en déduit que les contributions des termes avec $`n`$ impair sont nulles.

La géométrie et la combinatoire des graphes intervenant dans la formule (15) est clairement identique à celle du paragraphe 1.2. On étend donc toutes les notions rappelées au paragraphe 1.2.

On note de la même façon que précédemment le symbole de l’opérateur $`B_\mathrm{\Gamma }`$ lorsqu’on prend pour deux-vecteur de Poisson le crochet de Lie de $`𝔤_2`$. On a donc

$$\mathrm{\Gamma }(X,Y)=B_\mathrm{\Gamma }(\mathrm{exp}(X/2),\mathrm{exp}(Y),\mathrm{exp}(X/2))\mathrm{exp}(XY).$$
(16)
On en déduit, comme dans l’équation (8), que le membre de gauche dans (15) est un produit de deux termes de type exponentiel correspondant aux contributions des graphes simples de type roue et aux contributions des graphes simples de type Lie. La contribution dans (15) des graphes admissibles dont les composantes simples sont de type Lie calcule le terme de droite dans (14) (c’est le théorème $`5.1`$ de \[Ka\]). Par conséquent la contribution des graphes simples de type Lie dans (15) correspond à la formule de Campbell-Hausdorff pour les espaces symétriques par (14).

En considérant les graphes géométriques associés, on en déduit comme pour la formule (7) la proposition suivante :

###### Proposition 2.

La formule de Campbell-Hausdorff pour les espaces symétriques s’écrit pour $`X`$ et $`Y`$ dans $`𝔭`$ :

$$Z_{sym}(X,Y)=X+Y+\underset{\mathrm{\Gamma }}{}w_\mathrm{\Gamma }(\frac{1}{2})\mathrm{\Gamma }(X,Y)$$

avec la somme portant sur les graphes simples géométriques de type Lie avec trois points terrestres. En particulier on a bien $`Z_{sym}(X,Y)𝔭`$.

Pour les mêmes raisons que \[ADS\] lemme $`2.2`$ et $`2.3`$ ou \[AST\] proposition $`3.9`$, les séries associées sont convergentes pour $`X`$ et $`Y`$ dans $`𝔭`$ proches de $`0`$. On peut donc prendre $`h=1`$ et on obtient

$$\mathrm{exp}(X/2)\mathrm{exp}(Y)\mathrm{exp}(X/2)=D_{sym}(X,Y)\mathrm{exp}\left(Z_{sym}(X,Y)\right).$$
(17)
Notons $`J_𝔤(X)`$ le jacobien de l’application $`\mathrm{Exp}`$. On a la formule classique

$$J_𝔤(X)=\underset{𝔭}{det}\left(\frac{\mathrm{sinh}\mathrm{ad}X}{\mathrm{ad}X}\right).$$
(18)
En utilisant la fonction de densité (10) pour les algèbres de Lie et l’associativité de l’étoile-produit, on peut calculer sans difficulté la fonction de densité $`D_{sym}(X,Y)`$. On obtient

$$D_{sym}(X,Y)=\frac{j_{𝔤_2}(X/2)j_{𝔤_2}(Y)^{1/2}}{j_{𝔤_2}(Z_{sym}(X,Y))^{1/2}}$$
(19)
où la fonction $`j_{𝔤_2}`$ (11) est calculée pour le crochet double. Or pour $`X𝔭`$ on a

$$j_{𝔤_2}(X)=\mathrm{exp}(\mathrm{tr}_𝔤(\mathrm{ad}X))J_𝔤(X)^2$$

avec membre de droite calculé pour le crochet simple. On en déduit le lemme

###### Lemme 1.

On a la formule

$$D_{sym}(X,Y)=\frac{J_𝔤(X/2)^2J_𝔤(Y)}{J_𝔤(Z_{sym}(X,Y))}.$$
(20)
### 2.2 Première Déformation

On va effectuer maintenant une première déformation de la formule de Campbell-Hausdorff et de la fonction de densité. Elle consiste à déformer les cœfficients $`w_\mathrm{\Gamma }`$ en déplaçant les points terrestres dans le demi-plan de Poincaré (\[Ko\] $`\mathrm{8.2.3}`$).

La première déformation correspond typiquement au graphe de la figure  où la position du point associé à $`\mathrm{exp}(Y)`$ est repéré par la variable $`t`$. Cette variable désigne la position du point dans $`\overline{C_{\mathrm{1,2}}}^+`$, avec le point aérien sur la droite d’équation $`x=\frac{1}{2}`$. Cela revient à considérer des graphes simples de type Lie de $`G_{n\mathrm{,3}}`$, où on a placé les points de deuxième espèce en $`(0,t\mathrm{,1})`$.

On ajoute un indice $`t`$ pour signifier la dépendance en $`t`$. La règle concernant les symétries dans les cœfficients $`w_\mathrm{\Gamma }(t)`$ joue de la même façon que précédemment. Seuls les termes avec un nombre pair de sommets de première espèce vont intervenir.

On effectue comme dans \[To1\] paragraphe $`4.1`$ le calcul de la dérivée de la fonction

$$Z_{sym,t}(X,Y)=X+Y+\underset{\mathrm{\Gamma }}{}w_\mathrm{\Gamma }(t)\mathrm{\Gamma }(X,Y).$$
(21)
On ne doit considérer que les concentrations d’un point aérien sur $`\mathrm{exp}(Y)`$, les autres concentrations donnant soit directement $`0`$ (\[Ko\] lemme $`6.6`$) soit s’annulant par compensation due à l’identité de Jacobi. On se retrouve avec des cœfficients dérivés correspondant à des graphes du genre de la figure .

Le cœfficient $`w_{A,B}(t)`$ correspondant à un graphe comme dans la figure  contient un terme en $`dt`$. Pour les mêmes raisons que celles expliquées dans \[To1\] proposition $`4.1`$, le cœfficient $`w_{A,B}(t)`$ se factorise en un produit $`\stackrel{~}{w_A}(t)w_B(t)`$$`w_B(t)`$ se calcule à $`t`$ constant et $`\stackrel{~}{w_A}(t)`$ contient toute la contribution de la différentielle en $`t`$. En regroupant les graphes qui dégénèrent de la même manière, on fait apparaître un terme, au niveau des symboles, du genre

$$_YB(X,Y)[Y,A(X,Y)]$$
(22)
$`A,B`$ désignent par abus les symboles des graphes $`A,B`$. Ces graphes sont forcément des graphes simples admissibles de type Lie (avec points de deuxième espèce placés en $`(0,t\mathrm{,1})`$). Remarquons par ailleurs qu’il s’agit bien du crochet de Lie, car on a pris le deux-vecteur de Poisson associé au crochet de Lie. Pour des raisons de symétrie, on aura une contribution non triviale des cœfficients lorsque le graphe $`B`$ a un nombre pair de sommets de première espèce et le graphe $`A`$ un nombre impair de sommets de première espèce (à cause de l’arête qui joint le point maintenant aérien $`t`$ à la racine de $`A`$). La factorisation du cœfficient permet alors de tout factoriser comme dans \[To1\] théorème $`4.2`$ :

$$\underset{A,B}{}w_{A,B}(t)[Y,A(X,Y)]_YB(X,Y)=\underset{A,B}{}w_B(t)\stackrel{~}{w_A}(t)[Y,A(X,Y)]_YB(X,Y)$$
(23)
$$[Y,\underset{A}{}\stackrel{~}{w_A}(t)A(X,Y)]_Y\left(\underset{B}{}w_B(t)B(X,Y)\right)=[Y,F_t]_YZ_{sym,t}(X,Y).$$
(24)
La convention des signes de \[AMM\] $`I\mathrm{.2.1}`$, pour l’orientation des strates donne le signe final. Par ailleurs l’étiquetage des graphes n’a plus d’importance une fois qu’on a considéré le produit du cœfficient et du symbole (car il y a compensation des signes). On peut donc travailler directement avec les graphes géométriques.

Pour la fonction de densité on procède comme dans \[To1\] paragraphe $`5.2`$ et on obtient alors les mêmes équations que celles que nous avons établies dans \[To1\] théorème $`4.2`$ et théorème $`5.2`$. Par suite on peut énoncer le théorème suivant :

###### Théorème 1.

Soit $`(𝔤,\sigma )`$ une paire symétrique réelle. Il existe une série convergente en $`X,Y`$ dans un voisinage de $`0`$ dans $`𝔭`$, on la note $`F_t(X,Y)`$. Elle est à valeurs dans $`𝔨`$ et c’est une $`1`$-forme en $`t`$. Il existe une déformation de la formule de Campbell-Hausdorff pour les espaces symétriques, on la note $`Z_{sym,t}(X,Y)`$ et il existe enfin une déformation de la fonction de densité, on la note $`D_{sym,t}(X,Y)`$ telles que l’on ait

$$\mathrm{d}Z_{sym,t}(X,Y)=[Y,F_t(X,Y)]_YZ_{sym,t}(X,Y)$$
(25)
et

$$\mathrm{d}D_{sym,t}(X,Y)=[Y,F_t(X,Y)]_YD_{sym,t}(X,Y)+\mathrm{tr}_𝔤\left(D_YF_t\mathrm{ad}(Y)\right)D_{sym,t}(X,Y).$$
(26)
La $`1`$-forme $`F_t(X,Y)`$, correspond à

$$\underset{A}{}\stackrel{~}{w_A}(t)A(X,Y)$$
(27)
avec $`A`$ un graphe simple géométrique de type Lie (où les points de deuxième espèce sont placés en $`(0,t\mathrm{,1})`$) et $`w_A`$ le cœfficient calculé en ajoutant une arête du point $`t`$ vers le sommet de $`A`$ comme dans la figure . La notation $`D_YF_t(X,Y)`$ signifie que l’on a pris la différentielle partielle de $`F_t`$ en $`Y`$.

### 2.3 Modification de la fonction de densité pour les espaces <br>symétriques: cas résoluble et très symétrique

Les équations obtenues dans le théorème précédent ne sont pas satisfaisantes du point de vue des espaces symétriques. En effet au lieu de trouver une $`\mathrm{tr}_𝔤`$ on aurait aimé trouver une $`\mathrm{tr}_𝔨`$ dans (26).

Comme $`D_YF_t`$ et $`\mathrm{ad}(Y)`$ sont des endomorphismes de $`𝔤`$ qui envoient $`𝔨`$ sur $`𝔭`$ et $`𝔭`$ sur $`𝔨`$, on en déduit que l’on a

$$\mathrm{tr}_𝔭\left(\mathrm{ad}(Y)D_YF_t\right)=\mathrm{tr}_𝔨\left(D_YF_t\mathrm{ad}(Y)\right)$$
(28)
puis

$$\mathrm{tr}_𝔤\left(D_YF_t\mathrm{ad}(Y)\right)=2\mathrm{t}\mathrm{r}_𝔭\left(\mathrm{ad}(Y)D_YF_t\right)\mathrm{tr}_𝔨\left([D_YF_t,\mathrm{ad}(Y)]\right)$$
(29)
###### Lemme 2 (\[Rou1\] page $`573`$).

Pour les espaces symétriques résolubles et les paires très symétriques le commutateur $`[D_YF_t,\mathrm{ad}(Y)]`$ est à trace nulle sur $`𝔨`$.

On en déduit alors, que l’on peut transformer facilement les traces sur $`𝔤`$ en des traces sur $`𝔭`$ à condition que le terme $`\mathrm{tr}_𝔨([D_YF_t,\mathrm{ad}Y])`$ soit nul. C’est exactement l’argument de Rouvière. En particulier on trouve que ce terme est nul dans le cas des paires symétriques résolubles et des paires très symétriques (\[Rou1\] paragraphe $`5.3`$). Une paire très symétrique est la donnée d’une paire symétrique et d’un automorphisme commutant aux $`\mathrm{ad}X`$ et qui envoie $`𝔭`$ sur $`𝔨`$ et $`𝔨`$ sur $`𝔭`$ (\[Rou1\] paragraphe $`5`$). Pour les paires très symétriques on a facilement (\[Rou1\] paragraphe $`5.3`$ page $`573`$)

$$\mathrm{tr}_𝔨\left(D_YF_t\mathrm{ad}Y\right)=\mathrm{tr}_𝔭\left(D_YF_t\mathrm{ad}Y\right)=\mathrm{tr}_𝔨\left(\mathrm{ad}YD_YF_t\right).$$
(30)
Des exemples connus de paires très symétriques, sont les paires $`𝔤\times 𝔤/𝔤`$ et $`𝔤_{\mathrm{I}\mathrm{C}}/𝔤_{\mathrm{I}\mathrm{R}}`$. Dans le premier cas, on prend $`(X,Y)(X,Y)`$ et dans le second cas on prend la multiplication par $`i`$.

Conclusion:  Dans le cas des espaces symétriques il est utile de prendre la racine carrée de la fonction de densité pour obtenir une trace sur $`𝔨`$ (afin de compenser le cœfficient $`2`$ dans (29)).

### 2.4 Calcul de la fonction de densité à l’infini

On fait tendre $`t`$ vers l’infini, ce qui revient au même de fixer le point aérien en $`\frac{1}{2}+i`$ et de faire tendre les deux points terrestres vers $`\frac{1}{2}`$. On détermine facilement la valeur de la fonction de densité à l’infini.

Avant cela introduisons une autre fonction de densité que l’on note $`E_{1/2}(X,Y)`$, correspondant aux contributions des graphes admissibles dont les composantes simples sont de type roue et dont les points de deuxième espèce sont placés en $`\frac{1}{2}`$ (correspondant à $`\mathrm{exp}(X)`$) et en $`\frac{1}{2}+i`$ (correspondant à $`\mathrm{exp}(Y)`$) (voir figure ). Par ailleurs les points de première espèce sont associés au deux-vecteur de Poisson associé au crochet de Lie.

De manière analogue on notera $`\mathrm{exp}(Z_{1/2}(X,Y))`$ la contribution des graphes admissibles dont les composantes simples sont de type Lie et dont les points de deuxième espèce sont placés en $`\frac{1}{2}`$ et en $`\frac{1}{2}+i`$ (voir figure ).

###### Proposition 3.

Lorsque le point $`t`$ tend vers l’infini dans la direction imaginaire, on obtient

$$D_{sym,\mathrm{}}(X,Y)=E_{1/2}(X,Y)\frac{J_𝔤^2(X/2)}{J_𝔤(X)}$$
(31)
Preuve:  En approchant les deux points terrestres vers $`\frac{1}{2}`$, on doit tenir compte de toutes les concentrations de points aériens sur ces points terrestres. Cela ajoute un terme correspondant à la fonction de densité pour les algèbres de Lie (10), où l’on place $`\mathrm{exp}(X/2)`$ et $`\mathrm{exp}(X/2)`$ aux points terrestres, ce qui donne

$$\frac{j_{𝔤_2}(X/2)}{j_{𝔤_2}(X)^{1/2}}=\frac{J_𝔤^2(X/2)}{J_𝔤(X)}.\mathrm{}$$

### 2.5 Interprétation dans l’œil de Kontsevich

Lorsqu’on a rapproché nos deux points terrestres en $`\frac{1}{2}`$, on retrouve des graphes admissibles où les points de seconde espèce sont placés en $`(\frac{1}{2}+i,\frac{1}{2})`$. C’est donc un point dans $`\overline{C_{\mathrm{2,0}}}`$ positionné sur la paupière de l’œil et sur l’axe de symétrie voir figure (). Cela explique notre choix de notation : on est à la moitié sur la paupière.

### 2.6 Application à la convolution des distributions $`K`$-invariantes sur les espaces symétriques résolubles ou très symétriques

En reprenant d’une part les arguments énoncés dans \[AST\] paragraphe $`4`$ et \[KV\] paragraphe $`2`$ sur la convolution des germes de distributions invariantes et d’autre part les définitions de \[Rou1\] paragraphe $`3`$ pour la convolution sur $`𝔭`$, on déduit le théorème suivant (pour simplifier les notations on a noté $`J`$ pour $`J_𝔤`$) :

###### Théorème 2.

Soit $`(𝔤,\sigma )`$ une paire symétrique résoluble ou très symétrique. Soient $`u(x)`$ et $`v(y)`$ des germes de distributions $`\mathrm{ad}𝔨`$-invariantes en $`0`$ dans $`𝔭`$ à supports compatibles. On a alors la formule suivante

$`u(x)v(y),{\displaystyle \frac{J^{1/2}(x)J^{1/2}(y)}{J^{1/2}(z_{sym}(x,y))}}\mathrm{\Phi }(z_{sym}(x,y))=u(x)v(y),E_{1/2}^{1/2}(x,y)\mathrm{\Phi }(z_{1/2}(x,y))`$
pour $`\mathrm{\Phi }`$ une fonction test à support dans un voisinage de $`0`$ dans $`𝔭`$.

Preuve: Les problèmes de convergence et la définition de supports compatibles sont étudiés dans \[Rou1\] paragraphe $`3.3`$. La méthode est par ailleurs strictement identique. On se contente donc de la partie formelle.

On a introduit une déformation de la formule de Campbell-Hausdorff et de la fonction de densité. On a donc une dépendance en $`t`$

$$\mathrm{\Psi }(t)=u(x)v(y),D_{sym,t}^{1/2}(x,y)\mathrm{\Phi }(z_{sym,t}(x,y)).$$
(32)
On calcule facilement la dérivée en $`t`$, compte tenu des équations (25) et (26) du théorème 1. L’invariance du germe $`v(y)`$ signifie que l’on a

$$v(y)[y,F_t(x,y)]_y=\mathrm{tr}_𝔭\left(\mathrm{ad}(y)D_YF_t\right)v(y),$$
(33)
On obtient alors

$$\mathrm{\Psi }^{}(t)=u(x)v(y),m_t(x,y)D_{sym,t}^{1/2}(x,y)\mathrm{\Phi }(z_{sym,t}(x,y))$$
(34)
avec $`m_t(x,y)=\frac{1}{2}\mathrm{tr}_𝔨([D_YF_t,\mathrm{ad}y])`$. En effet l’invariance fait disparaître le terme en $`\mathrm{tr}_𝔭`$ dans (29) (n’oublions pas que l’on a pris la racine carrée). Lorsque $`(𝔤,\sigma )`$ est une paire très symétrique ou une paire résoluble, la trace sur $`𝔨`$ du crochet est nulle par le lemme 2. On en déduit que la fonction $`\mathrm{\Psi }`$ est constante.

Il vient $`\mathrm{\Psi }(0)=\mathrm{\Psi }(\mathrm{})`$, puis en utilisant (20) et (31)

$`u(x)v(y),{\displaystyle \frac{J(x/2)J(y)^{1/2}}{J(z_{sym}(x,y))^{1/2}}}\mathrm{\Phi }(z_{sym}(x,y))=`$
$`u(x)v(y),E_{1/2}^{1/2}(x,y){\displaystyle \frac{J(x/2)}{J(x)^{1/2}}}\mathrm{\Phi }(z_{1/2}(x,y)).`$ (35)
En changeant un peu l’ordre des termes et en simplifiant par $`J(x/2)`$ (la fonction $`J`$ est invariante par $`K`$ et l’énoncé vaut pour tous germes de distributions invariantes) on trouve

$`u(x)v(y),{\displaystyle \frac{J(x)^{1/2}J(y)^{1/2}}{J(z_{sym}(x,y))^{1/2}}}\mathrm{\Phi }(z_{sym}(x,y))=`$
$`u(x)v(y),E_{1/2}^{1/2}(x,y)\mathrm{\Phi }(z_{1/2}(x,y)).\mathrm{}`$ (36)
Conclusion: dans les cas de paires symétriques étudiées, on a trouvé une déformation adéquate vis à vis de la convolution des germes de distributions invariantes sur les espaces symétriques.

### 2.7 Deuxième déformation

On utilise maintenant la deuxième déformation, qui consiste à déplacer dans $`\overline{C_{\mathrm{2,0}}}`$ notre point de la paupière vers l’iris de l’œil tout en restant sur l’axe de symétrie vertical de l’œil. En d’autres termes, on fait se rapprocher les deux points de seconde espèce sur l’axe $`x=\frac{1}{2}`$.

C’est la déformation que l’on a utilisé dans \[To1\]. La situation est alors stricto sensu identique. Les équations que l’on obtient sont celles de \[To1\] théorème $`4.2`$ et théorème $`5.2`$. Compte tenu des symétries pour les cœfficients et le fait que $`X`$ et $`Y`$ sont dans $`𝔭`$, il est certain que nos champs adjoints sont dans $`𝔨`$. Comme on restreint notre étude au cas des espaces symétriques résolubles ou très symétriques, on peut grâce au lemme 2, transformer dans l’équation différentielle (26) la trace sur $`𝔤`$ en une trace sur $`𝔨`$ (ou sur $`𝔭`$) à condition de prendre la racine carrée de la fonction de densité.

Sur l’iris de l’œil cette fonction de densité vaut $`1`$ d’après \[Sh\] et la déformation de CBH vérifie $`Z_{1/2,iris}(X,Y)=X+Y`$. On a donc établi le théorème suivant (on note $`J`$ pour $`J_𝔤`$):

###### Théorème 3.

Soit $`(𝔤,\sigma )`$ une paire symétrique résoluble ou très symétrique. Soient $`u(x)`$ et $`v(y)`$ des germes de distributions $`\mathrm{ad}𝔨`$-invariantes en $`0`$ dans $`𝔭`$ à supports compatibles. Pour $`\mathrm{\Phi }`$ une fonction test à support dans un voisinage de $`0`$ dans $`𝔭`$, on a alors la formule suivante

$$u(x)v(y),\frac{J^{1/2}(x)J^{1/2}(y)}{J^{1/2}(z_{sym}(x,y))}\mathrm{\Phi }(z_{sym}(x,y))dxdy=u(x)v(y)\mathrm{\Phi }(x+y).$$

Remarques finales:

1. Quand $`(𝔤,\sigma )`$ est une paire symétrique résoluble, le théorème ci-dessus est démontré par Rouvière \[Rou1\] théorème $`3.6`$ en utilisant une déformation homogène de la formule de Campbell-Hausdorff comme dans \[KV\].
2. Les déformations utilisées dans cet article sont régulières mais leur jonction n’est pas différentiable. Cela explique notre choix d’exposition en deux parties et nous permet heureusement de retrouver les résultats de \[Rou1\].
3. Lorsque $`u(x)`$ et $`v(y)`$ correspondent à des distributions de support $`0`$, on retrouve que la symétrisation modifiée par la racine carrée du jacobien, constitue un isomorphisme d’algèbre de $`S[𝔭]^𝔨`$ (les symboles $`K`$-invariants) sur $`(U(𝔤)/U(𝔤)𝔨)^𝔨`$, l’algèbre des opérateurs différentiels invariants sur l’espace symétrique $`G/K`$. C’est l’isomorphisme de Rouvière dans le cas des espaces symétriques résolubles (voir \[To2\] pour une synthèse sur ce problème).
4. On sait qu’en général on ne peut pas obtenir un tel transport des distributions invariantes en utilisant une modification aussi simple de l’application exponentielle. On peut s’en convaincre par exemple en examinant ce qu’est l’homomorphisme d’Harish-Chandra du cas semi-simple (non très symétrique). Pour pallier cette difficulté Rouvière a introduit une fonction auxiliaire dans \[Rou2\], \[Rou3\], \[Rou4\] pour les fibrés en droites sur $`G/K`$, notée $`e_\lambda (x,y)`$. Il semble opportun de relier ces résultats avec les méthodes décrites dans cet article. Ce programme de recherche est annoncé dans \[To2\].

C. Torossian
UMR 8553 du CNRS, DMA, École Normale Supérieure
45 rue d’Ulm 75230 Paris cedex 05


text:
# Spin orientation of two-dimensional electron gas under intraband optical pumping

## I Introduction

In recent years spin-dependent phenomena in low-dimensional structures attract a great attention. One of the problems of common interest is a possibility and methods of the spin polarization in systems with a two-dimensional electron gas. The conventional way to achieve the spin polarization experimentally is the optical orientation of electron spins . Under optical excitation with the circularly polarized light the direct transitions from the valence to the conduction band are allowed only if the angular momentum is changed by $`\pm 1`$. These selection rules lead to the spin orientation of optically excited electrons, with the sign and the degree of polarization depending on the light helicity.

Up to now the theoretical consideration of the optical spin orientation was focused on direct interband transitions  and, partially, on direct intersubband transitions in the complicated valence band $`\mathrm{\Gamma }_8`$ of a zinc-blende lattice semiconductor . Here we report for the first time a theory of optical spin orientation under the indirect intrasubband excitation of electron gas by circularly polarized light and give an estimate of the spin generation rate. Note that the “monopolar” optical orientation under consideration can serve as a model of spin injection because the only type of carriers, electrons, is involved.

In addition, the short-range Elliot-Yafet mechanism of electron spin relaxation considered previously only for bulk semiconductors (, ch. 3) is extended to quantum-well structures. This spin-relaxation mechanism is shown to be governed by the same interband matrix elements of the electron-phonon or electron-defect interaction as those which govern the indirect intrasubband optical orientation.

## II Monopolar optical orientation of electron gas

The light absorption under intrasubband excitation in quantum wells (Drude-like absorption) is possible only if it is assisted by a phonon or a static imperfection in order to satisfy simultaneously both the energy and momentum conservation laws. Theoretically, these indirect optical transitions with initial and final states in the same conduction subband $`n`$ are described by second-order processes with virtual intermediate states. The compound matrix element for the indirect optical transition has the standard form

$$M_{ns^{}𝒌^{}ns𝒌}=\underset{\nu }{}\left(\frac{V_{ns^{}𝒌^{},\nu 𝒌}R_{\nu 𝒌,ns𝒌}}{E_{\nu 𝒌}E_{n𝒌}\mathrm{}\omega }+\frac{R_{ns^{}𝒌^{},\nu 𝒌^{}}V_{\nu 𝒌^{},ns𝒌}}{E_{\nu 𝒌^{}}E_{n𝒌}}\right),$$
(1)
where $`E_{n,𝒌}`$, $`E_{n,𝒌^{}}`$ and $`E_\nu `$ are the electron energies in the initial $`|n,s,𝒌`$, final $`|n,s^{},𝒌^{}`$ and intermediate $`|\nu `$ states, $`s`$ is the spin index, $`𝒌`$ is the electron wavevector, $`R`$ is the electron-photon matrix element, $`V`$ denotes the matrix element of interaction which allows the momentum transfer, e.g. scattering from impurities or phonon-assisted scattering. In the latter case the photon energy $`\mathrm{}\omega `$ is assumed to exceed the energy of an involved phonon.

The most probable processes are the transitions with intermediate states in the same subband (Fig. 1). This is the channel that determines the light absorption. However such transitions conserve the electron spin and, hence, do not contribute to the optical orientation.

In order to obtain the spin orientation of electron gas under intrasubband optical transitions one should take into account the virtual processes with intermediate states in the valence band, namely in the heavy, light and spin-split hole subbands. We assume that the carriers occupy the lowest electron subband $`e1`$. Then in the case of the flat QW with infinite barriers, the interband optical transitions between the $`e1`$ and $`\nu `$ states are allowed only for the hole subbands $`\nu =hh1,lh1`$ and $`so1`$. Fig. 2 demonstrates schematically the spin orientation under the $`\sigma ^+`$ normal-incidence excitation. In this case only the light and spin-split hole subbands contribute to the optical orientation. Because of the selection rules for interband matrix element, the electron transitions with the spin reversal $`|e1,1/2|e1,+1/2`$ occur via the intermediate states $`|lh1,\pm 1/2`$ and $`|so1,\pm 1/2`$, while the inverse processes $`|e1,+1/2|e1,1/2`$ are forbidden. While examining the optical orientation under oblique incidence, the heavy hole subband should also be taken into account as an intermediate state.

We consider the conduction electrons to form a nondegenerate 2D gas at the temperature $`T`$ or a degenerate 2D gas with the Fermi energy $`\epsilon _F`$, the interaction $`V`$ to describe the electron scattering by bulk acoustic phonons. Then an expression for the spin generation rate at $`\mathrm{}\omega k_BT`$ or $`\epsilon _F`$ ($`k_B`$ is Boltzmann’s constant) has the form

$$\dot{𝑺}=\frac{1}{6}\frac{\mathrm{\Xi }_{cv}^2}{\mathrm{\Xi }_c^2}\frac{\mathrm{\Delta }_{so}^2}{E_g(E_g+\mathrm{\Delta }_{so})(3E_g+2\mathrm{\Delta }_{so})}\left[𝒐_{}+\frac{a}{3}\sqrt{\frac{2m^{}\omega }{\mathrm{}}}𝒐_z\right]\eta I_0P_{circ}.$$
(2)
Here $`\mathrm{\Xi }_{cv}`$ and $`\mathrm{\Xi }_c`$ denote the interband and intraband deformation potential constants, $`\mathrm{\Delta }_{so}`$ and $`E_g`$ stand for the energies of the spin-splitting and the band gap, $`a`$ is the width of the quantum well, $`m^{}`$ is the electron effective mass, $`I_0`$ and $`P_{circ}`$ are the light intensity and the degree of circularly polarization inside the sample, $`𝒐_{}`$ and $`𝒐_z`$ are the in-plane and $`z`$ components of the unit vector o in the light propagation direction. The factor $`\eta `$ in Eq. (2) is the fraction of the light energy flux absorbed in the QW under phonon-assisted intrasubband optical transitions calculated for intermediate states in the same conduction subband, it is given by

$$\eta =\frac{3\pi \alpha }{n_\omega }\left(\frac{\mathrm{\Xi }_c}{\mathrm{}\omega }\right)^2\frac{k_BT}{\rho av_s^2}N_e,$$
(3)
where $`\alpha `$ and $`n_\omega `$ are the fine structure constant ($`1/137`$) and the refraction index of the medium, $`\rho `$ is the crystal density, $`v_s`$ is the sound velocity, and $`N_e`$ is the 2D carrier concentration.

An estimation for typical GaAs/AlGaAs structures shows that, at the comparable light intensities, the spin generation rate under intrasubband phonon-assisted optical transitions is by a factor of $`10^4÷10^5`$ smaller than that under interband excitation. However, consideration of other mechanisms of electron scattering, e.g. by impurity-assisted coupling, can essentially increase this ratio.

## III Spin relaxation mechanism induced by interband scattering

Besides the monopolar optical orientation, the virtual interband scattering described by the constant $`\mathrm{\Xi }_{cv}`$ can lead to a short-range mechanism of spin relaxation of the 2D electron gas. Microscopically this mechanism is connected with the $`𝒌p`$-induced admixture of the valence band states, $`\mathrm{\Gamma }_8`$ and $`\mathrm{\Gamma }_7`$, into the wave function of the conduction band $`\mathrm{\Gamma }_6`$ and the phonon- or defect-assisted interband coupling of these states. For bulk semiconductors this short-range Elliot-Yafet mechanism of electron spin relaxation was considered by Pikus and Titkov (see , ch. 3).

The spin relaxation time due to the mechanism under consideration can be calculated using the spin-flip matrix element (1) where $`\nu `$ are the valence band states $`\nu =hh1,lh1`$ and $`so1`$ and, in the matrix element $`R_{\nu 𝒌,ns𝒌}=(e/cm_0)𝑨𝒑_{\nu ,ns}`$, the vector $`(e/c)𝑨`$ is replaced by the in-plane momentum $`\mathrm{}𝒌`$. Here $`e`$ is the electron charge, $`c`$ is the light velocity in vacuum and A is the amplitude of the vector potential of the electro-magnetic wave.

Assuming the interband coupling $`V`$ to be caused by scattering on bulk acoustic phonons one can derive the relaxation rates for the in-plane and $`z`$ electron spin components which, after the averaging over the Boltzmann distribution, take the form

$`{\displaystyle \frac{1}{\tau _s}}`$ $`=`$ $`{\displaystyle \frac{1}{6}}{\displaystyle \frac{\mathrm{\Xi }_{cv}^2}{\mathrm{\Xi }_c^2}}{\displaystyle \frac{\mathrm{\Delta }_{so}^2k_BT}{E_g(E_g+\mathrm{\Delta }_{so})(3E_g+2\mathrm{\Delta }_{so})}}{\displaystyle \frac{1}{\tau _p}},`$ (4)
$`{\displaystyle \frac{1}{\tau _{sz}}}`$ $`=`$ $`{\displaystyle \frac{2\sqrt{\pi }a}{3}}\sqrt{{\displaystyle \frac{2m^{}}{\mathrm{}^2}}k_BT}{\displaystyle \frac{1}{\tau _s}},`$ (5)
where the momentum relaxation time determined by the acoustic phonon-assisted scattering is given by

$$\frac{1}{\tau _p}=\frac{3}{2}\frac{m^{}\mathrm{\Xi }_c^2}{\rho av_s^2\mathrm{}^3}k_BT.$$
(6)
One can see that, in case of the scattering by acoustic phonons, the spin relaxation of the 2D electron gas governed by the short-range Elliot-Yafet mechanism is very anisotropic, $`\tau _s\tau _{sz}`$.

### Acknowledgements

This work was supported by the Russian Ministry of Science, the RFFI (projects 00-02-16997 and 01-02-17528), the INTAS (project 99-00015), and by the Science Programme of the Presidium of RAS ”Quantum low-dimensional structures”.


text:
# Is there a true Model-D critical dynamics?

## Abstract

We show that non-locality in the conservation of both the order parameter and a noncritical density (model D dynamics) leads to new fixed points for critical dynamics. Depending upon the parameters characterizing the non-locality in the two fields, we find four regions: (i) model-A like where both the conservations are irrelevant (ii) model B-like with the conservation in the order parameter field relevant and the conservation in the coupling field irrelevant (iii) model C like where the conservation in the order parameter field is irrelevant but the conservation in the coupling field is relevant, and (iv) model D-like where both the conservations are relevant. While the first three behaviours are already known in dynamical critical phenomena, the last one is a novel phenomena due entirely to the non-locality in the two fields.

Based on the long time behaviour, the dynamics of different critical systems are classified HH  in various dynamic universality classes. These classes are characterized by several factors like the number of conserved field, non-vanishing Poisson brackets etc, in addition to the usual quantities necessary for equilibrium. Among the models which have no propagating modes, Model A is a system in which there are no conserved fields, Model B is a system with a conserved order parameter while Model C is the class of systems with a non-conserved order parameter field coupled to a conservative field. A fourth classification was also made - Model D, in which both the order parameter field and the coupled field are conserved. This model, however, was shown to be equivalent to Model B as the coupling to the field becomes irrelevant in the large length long time scale limit.

The dynamical behaviour in systems with one or more conserved field like Models B, C and D in the above classification is usually studied by considering local conservation. However, conservation is by no means necessarily local and in fact, non-local conservation is a more general condition. A few recent studies on critical or near-critical dynamics show that non-locality in the conserved field may change the dynamical class, besides the obvious modifications in nonuniversal details. This kind of effect of non-locality in the conservation has been noted in the study of critical dynamics of model B and C and in the numerical study of the early time effect of a model believed to belong to the class of model C and on the zero temperature dynamics (phase-ordering kinetics or Ostwald ripening) bray1 ; bray2 ; ps ; satya ; zheng ; conti . Besides the theoretical interest on its own right, one of the important applications of nonlocal dynamics is in speeding up numerical simulations, prompting implementation of various types of nonlocal movesmoves . Hence the need of a general classification of effects of nonlocality in dynamics.

In a sense, model D is the most generalised class of dynamical critical phenomena with dissipative dynamics. However, with local conservation, it is equivalent to Model B showing no novel features. Motivated by the results of model Cps , we investigate the problem of non-local conservation in Model D in the most general way by keeping both the conservations in the order parameter field and the coupled field non-local. One of our main results is to show the emergence of a true model D behaviour by nonlocal dynamics, unlike the local conservation case.

We introduce non-locality by using a non-local kernel bray1 ; bray2 ; ps  in the chemical potential keeping the continuity equation intact. Using $`\varphi `$ for the order parameter and $`m`$ for the secondary coupling parameter, the Hamiltonian of model D is given by

$`H={\displaystyle }d^d𝐱[{\displaystyle \frac{1}{2}}r\varphi ^2(𝐱)+{\displaystyle \frac{1}{2}}(\varphi (𝐱))^2+\stackrel{~}{u}\varphi ^4(𝐱)`$
$`+\gamma \varphi ^2(𝐱)m(𝐱)+{\displaystyle \frac{1}{2}}C^1m^2(𝐱)],`$ (1)
where $`\varphi ^2(𝐱)=_{i=1}^n\varphi _i^2(𝐱),\varphi ^4(𝐱)=[_{i=1}^n\varphi _i^2(𝐱)]^2,r=0`$ is the mean field critical point, and $`C>0`$ ensures noncriticality of $`m`$. There is a short distance cutoff which in momentum space is an ultraviolet cutoff $`\mathrm{\Lambda }`$. We shall set $`\mathrm{\Lambda }=1`$ in the calculation. The above Hamiltonian, Eq. ( Is there a true Model-D critical dynamics?) is the same as in Model Cps ; HH . Making $`\gamma =0`$ effectively reduces the above model to model A.

The dynamics 777 Additional torques for $`n=3`$ (Heisenberg model) may lead to new dynamic universality classesjdas . We do not consider such effects. is given by the following equations:

$`\left({\displaystyle \frac{1}{\mathrm{\Gamma }_\rho k^\rho }}+{\displaystyle \frac{1}{\mathrm{\Gamma }_0}}\right){\displaystyle \frac{\varphi _𝐤}{t}}`$ $`=`$ $`{\displaystyle \frac{\delta H}{\delta \varphi _𝐤}}+\eta (𝐤,t)`$ (2a)
$`\left({\displaystyle \frac{1}{\lambda _\sigma k^\sigma }}+{\displaystyle \frac{1}{\lambda _0}}\right){\displaystyle \frac{m_𝐤}{t}}`$ $`=`$ $`{\displaystyle \frac{\delta H}{\delta m_𝐤}}+\zeta (𝐤,t).`$ (2b)
The noises $`\eta `$ and $`\zeta `$ obey

$`\eta =0,\zeta =0,`$ (3a)
$`\eta (𝐤,t)\eta (𝐤^{},t^{})=2({\displaystyle \frac{1}{\mathrm{\Gamma }_\rho k^\rho }}+{\displaystyle \frac{1}{\mathrm{\Gamma }_0}})\delta (𝐤+𝐤^{})\delta (tt^{}),`$ (3b)
$`\zeta (𝐤,t)\zeta (𝐤^{},t^{})=2({\displaystyle \frac{1}{\lambda _\sigma k^\rho }}+{\displaystyle \frac{1}{\lambda _0}})\delta (𝐤+𝐤^{})\delta (tt^{}).`$ (3c)
Here, $`k=𝐤`$, and we have introduced two different parameters $`\rho `$ and $`\sigma `$ to denote the non-localities in the order parameter field and the coupling field respectively. The dissipative or the kinetic coefficients ($`\mathrm{\Gamma }_0`$ and $`\lambda _0`$) are denoted by the subscript $`0`$ while the transport coefficients ($`\mathrm{\Gamma }_\rho `$ and $`\lambda _\sigma `$) by the corresponding powers of $`k`$. For the local case, both $`\rho `$ and $`\sigma `$ assume a value equal to 2. Results are known for some limiting values of $`\rho `$ and $`\sigma `$: a) ($`\rho =0,\gamma =0`$) corresponds to model A HH , b) ($`\rho =2,\gamma =0`$) and ($`\rho =2,\sigma =2`$) give model B like behaviour HH , c) any value of $`\rho >0`$ and $`\gamma =0`$ is a non-local Model B bray1 ; satya , and d) $`\rho =0`$ with any $`\sigma >0`$ is non-local model C ps .

A Hamiltonian of the type of Eq.  Is there a true Model-D critical dynamics? occurs near the critical wings of a tricritical system or in critical systems with additional constraints. Model D dynamics are relevant in such cases, for example in three or four component mixtures or in spin-1 or Potts modelslawrie . A spin-$`\frac{1}{2}`$ example, though a bit unrealistic, would be an antiferromagnetic system with conserved magnetisation (as in model C ps2 ; zheng ) as well as conserved staggered magnetisation. The Kawasaki dynamics can be used to conserve the magnetisation maintaining at the same time that the exchanges between opposite spins have to be done in such a way that the order parameter is unchanged.

A momentum-shell renormalization group approach is used to study the long distance long time behaviour. In this approach, small scale fluctuations in space (between $`\mathrm{\Lambda }e^{\delta l}`$ and $`\mathrm{\Lambda }`$ in momentum space) are integrated out, and the cutoff rescaled to $`\mathrm{\Lambda }`$. One then obtains an effective Hamiltonian and effective equations for the dynamics valid for longer scales. The universal features are obtained from the fixed points of the flow equations for the various parameters of the problem. As per the usual RG scheme, any parameter that vanishes (grows) in the long scale limit is called an irrelevant (relevant) variable. A fixed point is then characterized by the set of relevant parameters. In this scheme the four dynamics classes correspond to the relevance and irrelevance of $`\lambda _\sigma `$ and $`\mathrm{\Gamma }_\rho `$ as mentioned in the abstract.

The statics of this system is the same as in Model C and the fixed point values of the parameters $`r,u`$ and $`\gamma ^2C`$ are as given in ref ps . The dynamical equations, on the other hand, take the forms:

$`{\displaystyle \frac{\lambda _\sigma ^1}{l}}`$ $`=`$ $`(z+\sigma +{\displaystyle \frac{\alpha }{\nu }})\lambda _\sigma ^1`$ (4a)
$`{\displaystyle \frac{\lambda _0^1}{l}}`$ $`=`$ $`(z+{\displaystyle \frac{\alpha }{\nu }})\lambda _0^1+(n\gamma ^2K_d)Q`$ (4b)
$`{\displaystyle \frac{\mathrm{\Gamma }_\rho ^1}{l}}`$ $`=`$ $`(z+\rho +2)\mathrm{\Gamma }_\rho ^1`$ (4c)
$`{\displaystyle \frac{\mathrm{\Gamma }_0^1}{l}}`$ $`=`$ $`(z+2)\mathrm{\Gamma }_0^1+(4C\gamma ^2K_d)Q{\displaystyle \frac{1}{1+\frac{Q}{CP}}}`$ (4d)
where $`Q=\mathrm{\Gamma }_\rho ^1+\mathrm{\Gamma }_0^1`$, $`P=\lambda _\sigma ^1+\lambda _0^1`$ and $`K_d=2^{1d}\pi ^{d/2}\mathrm{\Gamma }(d/2)`$. The parameters $`\lambda _\sigma ,\lambda _0,\mathrm{\Gamma }_\rho ,\mathrm{\Gamma }_0,\gamma `$ and $`C`$ appearing in the above equations are all $`l`$ dependent. To simplify notation, this $`l`$-dependence has been suppressed.

From (4b) we get

$$\frac{\mathrm{\Gamma }_0}{\lambda _0}=\frac{(1+x)n\gamma ^2K_d}{z\alpha /\nu }$$
(5)
where $`x=\frac{\mathrm{\Gamma }_0}{\mathrm{\Gamma }_\rho }`$ is the dimensionless ratio of the two coefficients for the primary order-parameter (note: $`\mathrm{\Lambda }=1`$). This quantity $`x`$ satisfies

$$\frac{x}{l}=(\rho Y)x,\mathrm{where}Y=\frac{2\alpha }{\nu n}\frac{1}{\mu ^1+(x+1)^1}.$$
(6)
Here we have introduced the ratios of the transport and kinetic coefficients of $`m`$ with the kinetic coefficient of $`\varphi `$, namely,

$$\mu _\sigma =\frac{\mathrm{\Gamma }_0C}{\lambda _\sigma },\mu _0=\frac{\mathrm{\Gamma }_0C}{\lambda _0}\mathrm{and}\mu =\mu _\sigma +\mu _0.$$

By using the fixed point value of $`C`$ and Eq. (5), $`\mu _0`$ can be rewritten as

$$\mu _0=\frac{(x+1)\alpha /\nu }{2(z\alpha /\nu )}.$$

The flow equation for $`\mu _\sigma `$ is

$$\frac{\mu _\sigma }{l}=\mu _\sigma [\sigma 2+\alpha /\nu Y].$$
(7)
We consider the different solutions of the equations for $`x`$ and $`\mu _\sigma `$. Corresponding to the relevance of the conserved quantities, one gets different regions in the $`\rho \sigma `$ plane with different dynamic critical behaviour. It should be mentioned here that for $`n4`$, the coupling $`\gamma `$ scales to zero and one effectively gets a non-local Model B. This has already been considered and shown to be model A-like when globally conserved bray1 . Our discussions are for the $`n<4`$ region where the coupling survives.

There are three fixed points of $`x`$: $`x=0`$, $`x=\mathrm{}`$ and a non-zero finite fixed point value for $`x`$ from Eq. (6):

Case 1. $`x=0`$ in the long length scale limit:

This corresponds to non-local Model C: the conservation of the order parameter is irrelevant. The fixed points and their stabilities of non-local model C, which were obtained in Ref. ps  from the solutions of $`\mu _\sigma `$, now depends on the value of $`\rho `$. In general, $`x=0`$ will be a stable fixed point as long as $`\rho <Y`$. We focus on the stability of the various regions for non-zero values of $`\rho `$.

(a) $`\mu _\sigma `$ has a finite fixed point value and $`z=\sigma +\alpha /\nu `$. This is a stable solution in the region $`(2/n)1>p>1`$ (where $`p=\frac{\sigma 2}{\alpha /\nu }`$) for $`\rho =0`$. For non-zero $`\rho `$ this remains stable for $`\rho <p+1`$.

(b) $`\mu _\sigma =\mathrm{}`$ : this is valid for $`p>(2/n)1`$. Here $`z=2+2\alpha /n\nu `$. This solution $`x=0`$, $`\mu _\sigma =\mathrm{}`$ is stable as long as $`\rho <2\alpha /n\nu `$.

(c) For $`\mu _\sigma =0`$, the stability of $`x=0`$ is valid only for $`p<1`$ and $`\rho <O(ϵ^2)`$. Here $`z=2+O(ϵ^2)`$ and the behaviour is model A like. Any non-zero value of $`\rho `$ thus destroys this model A like region. Therefore, a region where both the conservations are irrelevant is restricted strictly to $`\rho =0`$.

Case 2. A finite non-zero value of $`x`$ in the long length scale limit: This occurs when $`\rho =Y`$. In order that $`\mathrm{\Gamma }_0`$ reaches a finite fixed point value, we must have $`z=2+\rho `$. In this case, $`\mu _\sigma `$ again has the following fixed points:
(a) $`\mu _\sigma =\mathrm{}`$ for $`\sigma 2+\alpha /\nu >\rho `$
(b) $`\mu _\sigma =0`$ for $`\sigma 2+\alpha /\nu <\rho `$
(c) $`\mu _\sigma `$ nonzero finite for $`\sigma 2+\alpha /\nu =\rho `$

Physically a finite valued fixed point of $`x`$ means that the conservation in the order parameter $`\varphi `$ is relevant. If at the same time $`\mu _\sigma `$ attains a non-zero fixed point value, the conservation of the coupling field $`m`$ is also relevant. This occurs if the fixed points (a) and (c) are stable.

For case (a), $`x=\frac{\rho 2\alpha /\nu n}{2\alpha /\nu n}`$ and for non-negative value of $`x`$, $`\rho >2\alpha /n\nu `$. Hence in between $`\rho =2\alpha /n\nu `$ and $`\rho =\sigma 2+\alpha /\nu `$, a stable region is obtained where both the conservations are relevant.

For case (b), $`x+1=\frac{\rho (2z\alpha /\nu )}{2\alpha ^2/n\nu ^2}`$. This has a finite valued solution for $`x`$ only for $`\rho O(ϵ^2)`$. Hence for $`\rho >\sigma 2+\alpha /\nu `$, $`x\mathrm{}`$ is the only solution for $`\mu =0`$. It is a conventional model B fixed point.

Case (c): Here $`z=\sigma +\alpha /\nu `$ and in principle along the entire line $`\sigma 2+\alpha /\nu =\rho `$ a solution exists where both the conservations are relevant. However, there is a restriction on the value of $`\sigma `$ from the condition of existence of a non-zero finite value of $`\mu _\sigma `$. We find that

$$\mu _\sigma +\mu _0=\frac{(\sigma 2+\alpha /\nu )(x+1)}{2\sigma \alpha /\nu +\frac{2\alpha }{n\nu }(x+1)}.$$
(8)
Hence with $`\sigma =2+p\alpha /\nu `$; we get the condition $`(1+p)<2(x+1)/n`$. Since $`x`$ is also non-zero, this implies that this solution with both $`x`$ and $`\mu _\sigma `$ finite is valid for $`p+1=\rho /(\alpha /\nu )`$ with at least $`p<(2/n)1`$.

In the rest of the $`\rho \sigma `$ plane, the only solution is $`x=\mathrm{}`$ and $`1/\lambda _\sigma =0`$ for which only the conservation of the order parameter is relevant. Here $`z=\rho +2`$.

All the above possibilities of different dynamic behaviours are summarized in Fig. 1 in the $`\rho \sigma `$ plane (for $`n=1`$). The possibility of seeing a region where both the conservation conditions are relevant raises new issues like early time effect and boundary or surface effectsearly  in the new regime. These remain to be studied.

To summarize, Model D can be viewed as the most general dynamic model with purely dissipative dynamics and we indeed obtain all the four types of behaviours as soon as non-locality in the conservations (as in Eq. (2)) are introduced. Since our general model subsumes the previously studied nonlocal model C and local model D, it is natural to expect A-like, B-like and C-like regimes. The model A like region, however, turns out to be unstable for any finite value of $`\rho `$. The model C like regions are stable for small nonzero values of $`\rho `$ but ultimately disappear for larger values of $`\rho `$. The most non-trivial result, for general values of $`\rho `$ and $`\sigma `$, is the appearance (Region IV in Fig. 1) of a region with a new dynamical behaviour where both the conservations are relevant. This we call a true model D like region – a consequence of nonlocal conservation laws.

Acknowledgments: PS acknowledges financial support from DST project SP/S2/M-11/99.


text:
# Chandra Observation of SNR G54.1+0.3 — a Close Cousin of the Crab Nebula

## 1 Introduction

The Crab Nebula has played an essential role in our understanding of pulsar wind nebulae (PWNe). The basic observational X-ray characteristics of the Crab Nebula are a central pulsar, an ellipse-shaped torus, two jets and the absence of an outer thermal shell. The X-ray torus is explained as shocked equatorial wind material from the central pulsar (e.g., Aschenbach & Brinkmann 1975; Pelling et al. 1987; Weisskopf et al. 2000), while the jets correspond to the outflows from the two rotation poles of the pulsar (e.g., Aschenbach 1992). The observations and modeling of these characteristics have advanced our understanding of the energetics and geometry of pulsar wind, as well as the physics of ultra-relativistic shocks and particle acceleration (e.g., Rees & Gunn 1974; Aschenbach & Brinkmann 1975; Kennel & Coroniti 1984a, 1984b).

However, morphologically the Crab Nebula is unique. No other source has to date been found to mimic all the basic X-ray characteristics of the Crab Nebula (Gaensler 2001). The PWNe surrounding the Vela Pulsar and PSR B1509-58 have X-ray arcs, yet these X-ray arcs are not complete X-ray rings as observed in the Crab Nebula, and furthermore, these two PWNe are embedded inside bright thermal diffuse emission (Helfand, Gotthelf, & Halpern 2001; Gaensler et al. 2001).

In order to deepen our understanding of the properties of the pulsar wind and its interaction with the environments, it is important to find a PWN which closely resembles the Crab Nebula. Supernova remnant (SNR) G54.1+0.3 is such a candidate. Radio emission from G54.1+0.3 is centrally-filled, shows strong polarization, and has a flat spectrum, typical for a Crab-like SNR (Reich et al. 1985; Velusamy & Becker 1988). A nonthermal X-ray spectrum from ASCA observations confirmed its Crab-like nature (Lu, Aschenbach, & Song 2001). As we will present in this letter, our recent Chandra X-ray Observatory observation of this remnant shows that it is the closest “cousin” of the Crab Nebula.

## 2 Observation and Data Reduction

SNR G54.1+0.3 was observed with Chandra on 6-7 June 2001. The remnant was positioned at the aimpoint on the back-illuminated CCD chip S3 of the Advanced CCD Imaging Spectrometer (ACIS) in “faint” mode and at a working temperature of -120$`\mathrm{°}`$ C. The spatial resolution is $``$0$`\stackrel{}{\mathrm{.}}`$5, and the spectral resolution is $``$120 eV at 2.0 keV. The CCD frame read-out time for this observation is 3.2 seconds.

We calibrated the data using CIAO (2.1.3). After excluding two time intervals of high background we obtained a net exposure of 30.9 ks. Based on a comparison of the positions of the X-ray point source CXOU J193034.62+185110.8 (southeast of the remnant, see Fig. 1) and its optical counterpart, we find that the absolute positional accuracy of the X-ray observations is about 0$`\stackrel{}{\mathrm{.}}`$6, consistent with the position accuracy claimed in the Chandra Proposer’s Observatory Guide.

## 3 Analysis and Results

While Fig. 1 presents an overview of the Chandra data, Fig. 2 is a close-up of G54.1+0.3, with the image data adaptively smoothed using a Gaussian to achieve a signal-to-noise ratio of $``$5. The structure of the remnant is a complex, but we decompose it into several components (Fig. 1): a central bright point-like source (CXOU J193030.13+185214.1), a surrounding ring, elongations to the east and west, and extended low surface brightness emission in between and beyond the central structures. We find no evidence for a shell-like structure that would correspond to the thermal emission usually expected for an SNR. Therefore the remnant resembles the Crab Nebula in morphology quite closely.

Fig. 3 displays spectral model fits of the individual components of G54.1+0.3 and Table 1 lists the corresponding model parameters. The background used in the spectral analysis is extracted from an annulus centered on CXOU J193030.13+185214.1 and with the inner and outer radii equal to 70$`\mathrm{}`$ and 130$`\mathrm{}`$ (with CXOU J193034.62+185110.8 removed). The spectrum of each of these components can be well fitted with a power-law model. There is no evidence for any emission line, and fitting the spectra with thermal models yields temperatures around 10 keV, which is much higher than the plasma temperature ($``$3 keV) typically observed for young SNRs. Therefore we detect no significant thermal emission. We assume that the X-ray absorption column density is the same across the field as the absorption column densities of all these components derived from individual spectral fitting are similar. The jointly fitted value (1.6$`\pm 0.1\times 10^{22}`$ cm-2) is about half of the total Galactic X-ray absorption ($``$3.6$`\times `$1022 cm-2) derived from the IRAS 100 $`\mu `$m emission in this direction (Boulanger & Perault 1988; Wheelock et al. 1994). Assuming that the absorbing gas is relatively uniformly distributed along the line of sight, the distance to G54.1+0.3 should be approximately 5 kpc, about half way to the edge of the Galaxy from the Sun.

The source CXOU J193030.13+185214.1 most likely represents the pulsar that powers G54.1+0.3. The radial surface brightness profile of this source is consistent with the point spread function of the instrument. The slow CCD frame read-out time (3.2 s) of the observation, however, prevents a search for any periodicity typical of a young pulsar. We find no significant signal with a period $`>`$6 s. The count-rate of the source, 0.064$`\pm `$0.002 cts s-1, indicates a pile-up fraction of only $``$6$`\%`$ with a negligible effect on the spectrum. The relation between the photon spectral indices of the source and the nebula is well consistent with the one derived by Gotthelf & Olbert (2001) for pulsars and their respective nebulae. We find no emission peaks at the source position at radio (with 4.8 GHz VLA data, Velusamy & Becker 1988) and infrared (with the online 2MASS data) wavelengths. From our 0.1-2.4 keV X-ray luminosity (1.8$`\times `$1033 $`d_5^2`$ erg s-1,where $`d_5`$ is the distance in units of 5 kpc) of this point source, we infer the pulsar spin-down luminosity $`\dot{E}`$ as $``$2$`\times `$1036 $`d_5^2`$ erg s-1 (Becker & Trümper 1997).

An eyeball fit to the ring with an ellipse centered on CXOU J193030.13+185214.1 gives the north-south semi-major axis as 5$`\stackrel{}{\mathrm{.}}`$7 and the east-west semi-minor axis as 3$`\stackrel{}{\mathrm{.}}`$7. The eastern part of the ellipse is much brighter than the western part. The spectrum of the ring is the hardest among the diffuse components.

The two elongations are not as well defined as the ring (Figs. 1-2). The western elongation is oriented perpendicular to the ring and extends $``$32$`\mathrm{}`$ horizontally from near the central point source to (RA, Dec; J2000) $`19^\mathrm{h}30^\mathrm{m}28^\mathrm{s}`$, $`18\mathrm{°}52^{}14\mathrm{}`$. The eastern elongation extends $``$14$`\mathrm{}`$ ending at about 19$`{}_{}{}^{\mathrm{h}}30_{}^{\mathrm{m}}31^\mathrm{s}`$, $`18\mathrm{°}52^{}10^{\prime \prime }`$. The orientation has an average inclination of $``$18$`\mathrm{°}`$ to the south.

The outer low surface brightness diffuse emission marks the accumulated pulsar wind material. The size of this low surface brightness emission (1$`\stackrel{}{\mathrm{.}}5\times 1\stackrel{}{\mathrm{.}}2`$) is much larger than that of the central bright region (0$`\stackrel{}{\mathrm{.}}7\times 0\stackrel{}{\mathrm{.}}3`$) composed of the point-like source, the ring and the elongations. The low surface brightness component shows some extensions to the west, north and northeast and compares well with the radio map both in size and overall shape (Velusamy & Becker 1988). The spectrum is significantly softer than those of the ring and western elongation, suggesting that, if we assume synchrotron radiation, the particles in this region are more evolved. A detailed comparison and explanation of the X-ray and radio maps and a more extensive study of the low surface brightness diffuse emission will be presented in a separate paper.

## 4 Comparison with the Crab Nebula

The X-ray morphology of G54.1+0.3 is strikingly similar to that of the Crab Nebula although some differences exist. Both remnants show a complete X-ray ring (torus) around the corresponding pulsar, plus two opposite elongations (jets), as well as the lack of a thermal outer shell. But the elongations of G54.1+0.3 are more diffuse than the Crab jets, that we will discuss more specifically in $`\mathrm{\S }`$5.2. G54.1+0.3 shows large scale low surface brightness X-ray emission beyond the ring and the elongations. This emission coincides well with the radio morphology of the remnant (Velusamy & Becker 1988). In contrast, little X-ray emission is presented beyond the Crab torus (Weisskopf et al. 2000). The radio size of the Crab Nebula is $``$3 times bigger than the X-ray size (Bietenholz & Kronberg 1990; Weisskopf et al. 2000). The morphological difference of the outer low surface brightness emission suggests that the synchrotron cooling efficiency in the G54.1+0.3 ring is much lower than in the Crab torus.

G54.1+0.3 and the Crab Nebula have quite similar spectral characteristics. The spectra of both remnants are characterized by power-laws which steepen progressively from pulsar, through torus (ring) and jets (elongations), to the outer diffuse emission (Table 1; Willingale et al. 2001). The similar spectral steepening trend indicates that particles follow similar tracks in the two remnants. However, the spectra of G54.1+0.3 are flatter than those of the corresponding components in the Crab Nebula, consistent with lower efficiency of the synchrotron cooling in G54.1+0.3.

Both the total X-ray luminosity and the luminosities of the corresponding components in G54.1+0.3 differ substantially from those of the Crab Nebula (see Table 1 and Weisskopf et al. 2000). The putative pulsar and the X-ray nebula of G54.1+0.3 are about two orders of magnitude less luminous than those of the Crab. The ring of G54.1+0.3 accounts for only about 10$`\%`$ of the total extended emission, but the torus of the Crab Nebula dominates the X-ray emission of the whole nebula. In G54.1+0.3, the luminosities of the elongations are comparable with that of the ring, whereas in the Crab Nebula, the jets are obviously much weaker than the torus. These differences are likely due to the particular pulsar wind geometry and energy input to the various components in the two remnants.

## 5 Pulsar Wind Properties

We suggest that the X-ray ring and elongations are due to the equatorial wind and the two polar outflows from the pulsar, respectively, as illustrated schematically in Fig. 4.

### 5.1 The X-ray Ring

The elliptical ring in G54.1+0.3 is apparently the projection of an inclined circular ring. Similar to the X-ray torus in the Crab Nebula, this ring might be due to the shocked pulsar equatorial wind (Aschenbach & Brinkmann 1987). The ratio of observed semi-major and semi-minor axes suggests an inclination angle $`\theta `$=41$`\mathrm{°}`$.

Following Pelling et al. (1987), we calculate whether the Doppler boosting of the downstream motion can explain the X-ray brightness variation across the ring (Fig. 5). The intensity ($`I`$) of the ring at position angle $`\varphi `$ measured counter-clockwise from the north can be expressed as

$$I=I_0[\frac{\sqrt{1v^2}}{1\mathrm{cos}(\theta )v\mathrm{sin}(\varphi )}]^{\alpha +1}+I_b,$$
(1)
where $`I_0`$ is the unshifted intensity, $`v`$ is the bulk motion velocity in units of $`c`$ — the speed of light, $`\alpha `$ is the photon spectral index, and $`I_b`$ is the intensity of the diffuse background emission. Fig. 5 shows that the Doppler boosting model fits the measured intensity variation well. The fitted parameters and their 1$`\sigma `$ uncertainties are: $`I_0`$=7.6$`\pm `$2.4 cts arcsec-2, $`v`$=0.40$`\pm `$0.12 c, and $`I_b`$=5.7$`\pm `$2.5 cts arcsec-2. Using the $`v`$ value and following Kennel & Coroniti (1984b) we infer a ratio of the electromagnetic energy flux to the particle energy flux in the wind of $`\sigma `$$``$0.06. Therefore, the surface brightness variation across the ring is consistent with the assumption that the enhanced X-ray emission arises from the freshly shocked pulsar wind material.

### 5.2 The Elongations

We propose that the elongations are associated with the pulsar polar outflows. This association is suggested by the morphology of the western elongation (Figs. 1-2). The spectrum of the elongation is flatter than that of the surrounding low surface brightness emission, suggesting that the particles in the western elongation have undergone fewer synchrotron losses. The eastern elongation is not well resolved, but appears bent and shorter than the western elongation. These two elongations, however, are broader and more irregular than the relatively well confined jets observed in Crab (Weisskopf et al. 2000), Vela (Helfand et al. 2001), and PSR B1509-58 (Gaensler et al. 2001), indicating that the pulsar polar flows in G54.1+0.3 are less collimated.

### 5.3 Geometry

The X-ray morphology of the ring and the two elongations in G54.1+0.3 suggest that the pulsar wind is concentrated in the equatorial plane and the polar outflow directions. The width of the ring is about 2$`\stackrel{}{\mathrm{.}}`$5 in the minor axis direction. As the radius of the deprojected ring is 5$`\stackrel{}{\mathrm{.}}`$7, the width indicates that the equatorial wind particles fan out in a half opening angle of $``$16$`\mathrm{°}`$, corresponding to $``$1.1$`\pi `$ steradian. The width of the head of the western elongation is measured to 6′′, implying that the half opening angle of the western polar outflow is $``$4$`\mathrm{°}`$, assuming that the outflow is perpendicular to the equatorial plane as implied by Fig. 2. The corresponding solid angle is 5$`\times 10^3`$$`\pi `$ steradian. The open angle of the eastern outflow remains unclear because of the insufficient knowledge about its direction and thus length.

The observed pulsar wind geometry reflects the fractional wind power injected to the ring and elongations. Because the reverse shock occurs when the wind ram pressure balances the diffuse nebula pressure, $`p_n`$, we have

$$r_s=(\frac{\delta \dot{E}}{\mathrm{\Omega }cp_n})^{1/2},$$
(2)
where $`r_s`$ is the distance from the pulsar, $`\delta `$ is the fraction of the wind power $`\dot{E}`$ in one particular wind component, and $`\mathrm{\Omega }`$ is the solid open angle of that component. We assume that $`\dot{E}`$ is split between $`\dot{E}_e`$ (injected to the equatorial ring) and $`\dot{E}_p`$ (released in the polar outflows) that is equally divided between the two polar outflows. We also assume that $`p_n`$ does not change substantially with position in the diffuse nebula, because the sound speed in a relativistic plasma is $`c/\sqrt{3}`$ and so no pressure gradient in the bubble can be maintained on time scales that are dynamically important (e.g., Rees & Gunn 1974). From Equation (2), with the geometric parameters and assumptions, we infer that $`1/3`$ of $`\dot{E}`$ is released to the polar outflows.

## 6 Summary and Conclusions

Based on a Chandra ACIS-S observation, we have resolved the spatial and spectral structures of SNR G54.1+0.3. We find that both its morphology and the spectra closely resemble those of the Crab Nebula, though the synchrotron cooling efficiency in G54.1+0.3 may be substantially lower than in the Crab Nebula. In particular, we have discovered a point-like source, CXOU J193030.13+185214.1, at a location close to the morphological center of G54.1+0.3. The source is most likely the pulsar powering the SNR, although a spin period of the pulsar is yet to be determined. We estimate the pulsar spin-down luminosity to be $``$2$`\times `$1036 $`d_5^2`$ erg s-1 from the X-ray luminosity of the point source. We suggest that the diffuse X-ray emission around the pulsar originates from the shocked pulsar wind. The presence of the X-ray ring indicates a strong equatorial wind with an inclination angle of 41$`\mathrm{°}`$. From the surface brightness variation across the ring we infer the downstream flow speed of the shocked pulsar wind material to be 0.40$`\pm `$0.12 c, which indicates that the wind is particle dominated. The two elongations are oriented nearly perpendicular to the ring and are likely to be associated with the polar outflows from the pulsar. The geometry of the X-ray ring and the elongations suggest at least one-third of the spin-down energy is injected into the two polar outflows. These characteristics of SNR G54.1+0.3 provide new constraints for modeling of the elusive pulsar wind and its interaction with the environment.

We thank the referee D. J. Helfand for his comments. This work is supported partially by the NASA-grant SAO GO-12068X and NASA LTSA grant NAG5-7935. FJL and LMS also appreciate support of the Special Funds for Major State Basic Research Projects and the National Natural Science Foundation of China.