EL Arbaoui Mohamed

Mini projet sur la double provskite

1- Nom de compos mon compos est le Ba2VNbO6

2- La structure cristalline

Ba2VNbO6 se cristallise dans la structure cristallin cubique, o le vanadium occupe le milieu

des arrtes et niobium occupe les sommet et les centre des faces de cube , et loxygne faire les
coins sur V et Nb.

Paramtre de la maille a=7.918 angchtrom et = = = 90

Je ne reprsente pas le Ba sur la structure, car les atomes qui plus intressent pour mois cest V et
Nb ce sont des atomes qui peuvent ton possde un moment magntique

La structure cristalline de Ba2VNbO6

Les atomes blancs represente le Nb

Les atomes noirs reprsentent le Vanadium

Les atomes en gris reprsentent loxygne

Pour le Ba2VNbO6 les atomes magntiques sont V et Nb, alors pour dtermin les moment on a
besoin de connaitre la valence de chaque atome.

3- Configuration lectronique de Ba et V et Nb et O

Ba : [Xe] 6 2 donc 2+ est ltat ionique stable

O :1 2 2 2 24, ltat ionique la plus stable de O est 2
V : [Ar] 33 4 2 les tats doxydation sont 2+ 3+ 4+ 5+
Nb : [Kr] 44 51 les tats doxydation sont 3+ 5+

La neutralit de notre compos (6*-2)+2*2+x+y=0 donc on a : -8+x+y=0 partir des

possibilit de V et Nb il prendre 3+ et 5+
4- La valance des atomes

Pour vanadium : on a 3+ sa configuration est: [Ar] 33 4 2 donc leur valeur de spin est
2*1/2=1 le spin prendre donc 3 valeur soit +1 ou -1 ou bien 0.

Pour niobium : on a 5+ sa configuration est [Kr] 44 51 donc leur valeur de spin est S=0, le
niobium ne possde donc dans ce cas un moment magntique.
5- Hamiltonien de Ba2VNbO6

Puisque on a un seul atome magntique, on peut liminer dans notre structure tout
les autres atomes en garde seulement les atomes qui sont magntique c..d. le vanadium,
la structure devient cubique face centr o les sommets et les centre des faces occup
par le vanadium.

-Dans ce cas on a un seul couplage entre les atomes de vanadium qui sont spar par la
distance a2 o a est le paramtre de la maille.

-le systme est ferromagntique

-puisque la valeur de spin est dfrent de 1/2 on peut tien compte un champ cristallin
qui va minimiser lnergie.

Alors Hamiltonien de Ba2VNbO6 est :

 = 2

J : le couplage vanadium-vanadium

: : 1,0,1

6-cherchons ltat fondamental

Pour chercher le fondamental il faut calculer les nergies des tats possible,
lexpression de lnergie est donne par

= ( ) 2
2

a-lnergie de la phase (1,1) ou encore de (-1,-1) car les deux phase sont quivalant :

(1,1) = 6

b-la phase (1,-1) est quivalant (-1,1), leur nergie :

(1,1) = 6

c- la phase (0,0)

(0,0) = 0

7-diagramme de phase
Pour J>0 et > 0 la phase qui minimise lnergie est la phase (1,1)

Pour > 0 et J<0 la phase (1,-1) qui est minimise lnergie

Pour < 0 et J<0 ou J>0 la phase qui minimise lnergie est la phase (0,0)

Affin de tracer le diagramme de phase on prendre les galits suivant :

E(1,1)=E(0,0) donc 6 = 0 on peut tirer donc lquation de la droite = 6

De mme E (-1,1)=E(0,0) donc 6 = 0 on peut tir donc lquation de la droite =

6

Et pour E(1,1)=E (-1,1) donc on a 6 = 6 = 0

8-dterminations des positions des atomes
Si je prendre latome qui se trouve lorigine c..d. la position (0,0 ,0) il possde 12
proches voisins leur positions sont :

(1/2,1/2,0) ; (-1/2,1/2,0) ; (1/2,-1/2,0) ; (-1/2,-1/2,0) ; (0,1/2,1/2) ; (0,-1/2,1/2) ; (0,1/2,-1/2) ;

(0,-1/2,-1/2) ; (1/2,0, 1/2) ; (-1/2,0, 1/2) ; (1/2,0, -1/2) ; (-1/2,0,- 1/2)

De mme pour chaque atome on peut dterminer les positions de leur voisin.

9-dcomposition de rseau en sous rseau

Affin de faciliter la simulation de notre systme, on va le dcompos en deux
systme, on se basant sur la valeur de moment de spin comme critre de dcomposition,
autrement dit on prendre les atomes qui possdant la mme valeur de spin dans mme
sous rseau on obtient donc deux sous rseau. Pour une maille lmentaire
les position des atomes de sous rseau 1 est : (0,0,0) ; (1,0,0) ; (1,1,0) ; (0,0,1) ; (0,1,0) ;
(1,0,1) ; (0,1,1) ; (1,1,1)

Les positions des atomes de sous rseau 2 est : (1/2,1/2,0) ; (1,1/2,1/2) ; (1/2,1/2,1) ;
(0,1/2,1/2) ; (1/2, 0 ,1/2) ; (1/2, 1 ,1/2)

Les deux sous rseau sont des structures cubiques simples, on peut donc calculer
laimantation facilement, et pour calculer laimantation globale de notre systme il suffit de
prendre

m=m1+m2 cas de ferromagntique

m=m1-m2 cas de antiferromagntique

10- gnralisation des positions des atomes et condition priodique

Pour calculer les alimentations et tout autre paramtre, on a besoin de connaitre les
positions des atomes sur notre systme, ce l nous permettre de simplifier la simulation de
notre systme.

Alors pour calculer laimantation on va prendre lorigine de repre sur un atome de

rseau 1, et on ajoute chaque fois un pas de 1 suivant x et y et z.

Si je prendre par exemple un rseau 4*4*4 et si je travaille seulement suivant laxe des x ( la
mme chose se rpte pour les direction y et z) alors la premire atome possde comme
coordonnes (0,0,0) la 2 (1,0,0) la 3 (2,0,0) , le 4 (3,0,0) et la dernire (4,0,0) .
pour dtermin les positions des atomes il faut dajouter un pas de 1 suivant x et suivant y et
suivant z
 Alors pour dtermin les positions des atomes de sous rseau 2 on fait la mme chose
sauf que dans ce cas on prendre latome (1/2,1/2,1/2) comme origine.

Comment on peut faire celle l pratiquement ?

Les tapes expliqu au niveau de cette paragraphe peut se faire pratiquement, il
suffit de boucler sur tout les i et j et les k on commence par i=j=k=0 pour le rseau 1 et
i=j=k=1/2 pour le sous rseau 2 jusqu N avec un pas de 1.

Comment ce faire les conditions priodique ?

Les systmes physiques sont gnralement infinis, alors pour les modalis et simul, il
faut compte cette caractristique dans le programme de simulation pour obtenir des
comportement semblable celle rel.

La faon de reprsenter cette caractristique est de mettre des conditions priodiques c..d.
de boucler le systme.

Exemple :

Si on prendre par exemple le cas dun rseau 4*4*4

Latome qui occupe la position (4, 4,4) en le boucl avec (0, 0, 0)

Latome qui occupe la position (4, j, k) en le boucl avec (0, j, k)

Latome qui occupe la position (i, 4, k) en le boucl avec (i, 0, k)

Latome qui occupe la position (i, j, 4) en le boucl avec (i, j, 0)

11- le programme
parameter (N=100,Jc=0.89, Delta=10)! Jc : couplage d'interaction entre
vandium_vandium , N : largeur du rseau carr , Netap : tapes monte Carlo

real T1(N,N,N),T2(N,N,N) !T(N,N,N): matrice qui donne l'tat du spin dans chaque site

real Mi1,Mi2,SM1,SM2,Ei,SE,Mt!Mf,Mf Mi,Mf : aimantation initiale,finale , Ei,E : nergie

initiale,finale

open(1,file='aimantation initiale.dat')

open(2,file='energie initiale.dat')
 ! Initialisation de la matrice qui donne l'etat de spin

do i=0,N,1

do j=0,N,1

do k=0,N,1

T1(i,j,k)= 0.

enddo

enddo

enddo

do i=1./.2,N,1

do j=1./.2,N,1

do k=1./.2,N,1

T2(i,j,k)= 0.

enddo

enddo

enddo

Mi1=0

SM1=0

do i=0,N,1

do j=0,N,1

do k=0,N,1

x=rand(0)

if(x .le. 1./3. ) then

T1(i,j,k)=-1

endif
 if (x.gt.1./3. .and. x .le. 2./3. ) then

T1(i,j,k)=0

endif

if (x .gt. 2./3.)then

T1(i,j,k)=1

endif

SM1=SM1+T1(i,j,k)

enddo

enddo

enddo

Mi1= SM1/(N*N*N)

write(*,*)'aimentation initial 1',Mi1

Mi2=0

SM2=0

do i=1/2,N,1

do j =1/2,N,1

do k =1/2,N,1

x=rand(0)

if(x .LE. 1./3. ) then

T2(i,j,k)=-1

else if (x.LE. 2./3. ) then

 T2(i,j,k)=0

else if (x.LE. 2./3. ) then

T2(i,j,k)=1

endif

SM2=SM2+T2(i,j,k)

enddo

enddo

enddo

Mi2= SM2/(N*N*N)

write(*,*)'aimentation initial 2',Mi2

Mt=Mi1+Mi2

write(1,*)'aimentation initial ',Mt

SE1=0

do i=0,N,1

do j=0,N,1

do k=0,N,1

SE1=SE1+ T(i,j,k)* T(i+1/2,j+1/2,k) + T(i,j,k)* T(i-1/2,j+1/2,k)&

T(i,j,k)* T(i+1/2,j-1/2,k)+ T(i,j,k)* T(i-1/2,j-1/2,k) &

+T(i,j,k)*T(i,j+1/2,k+1/2)++T(i,j,k)*T(i,j-1/2,k+1/2)&

& +T(i,j,k)*T(i,j+1/2,k-1/2)++T(i,j,k)*T(i,j+1/2,k+1/2)&
& +T(i,j,k)*T(i,j-1/2,k-1/2)+T(i,j,k)*T(i+1/2,j,k+1/2)&

& +T(i,j,k)*T(i-1/2,j,k+1/2)+T(i,j,k)*T(i+1/2,j,k-1/2)&

& +T(i,j,k)*T(i-1/2,j,k-1/2) - Delta/Jc * T(i,j,k)*T(i,j,k)

enddo

enddo

enddo

Ei1 = - Jc * SE1 ! E = -J SOMME(Si*Sj)avec <i,j> proches voisins

write(*,*)'\n l''nergie initiale = ',Ei1

SE2=0

do i=1/2,N,1

do j=1/2,N,1

do k=1/2,N,1

SE2=SE2+ T(i,j,k)*T(i+1/2,j+1/2,k) + T(i,j,k)*T(i-1/2,j+1/2,k)&

T(i,j,k)*T(i+1/2,j-1/2,k)+T(i,j,k)*T(i-1/2,j-1/2,k)

+T(i,j,k)*T(i,j+1/2,k+1/2)++T(i,j,k)*T(i,j-1/2,k+1/2)&

& +T(i,j,k)*T(i,j+1/2,k-1/2)++T(i,j,k)*T(i,j+1/2,k+1/2)&

& +T(i,j,k)*T(i,j-1/2,k-1/2)+T(i,j,k)*T(i+1/2,j,k+1/2)&

& +T(i,j,k)*T(i-1/2,j,k+1/2)+T(i,j,k)*T(i+1/2,j,k-1/2)&

& +T(i,j,k)*T(i-1/2,j,k-1/2) - Delta/Jc * T(i,j,k)*T(i,j,k)

enddo

enddo

enddo

Ei1 = - Jc * SE2 ! E = -J SOMME(Si*Sj)avec <i,j> proches voisins

 Ei=Ei1+Ei2

write(*,*)'\n l''nergie initiale = ',Ei2

write(2,*)'energie total=',Ei

do Istep=1,Nstep ! boucle des tapes Monte Carlo

l=int(rand(0)*(N)) !on tire un site de faon alatoire

m=int(rand(0)*(N))

n =int(rand(0)*(N))

if(l.ne.1.and.l.ne.N.and.m.ne.1.and.m.ne.N.and.n .ne.1.and.n.ne.N)then

som= T(l-1/2,m-1/2,n)+ T(l-1/2,m+1/2,n)+T(l+1/2,m-1/2,n)+

T(l+1/2,m+1/2,n)+T(l,m-1/2,n-1/2)+T(l,m+1/2,n-1/2)+

T(l,m-1/2,n+1/2)+T(l,m-1/2,n-1/2)+T(l+1/2,m,n-1/2)+

T(l-1/2,m,n+1/2)+ T(l+1/2,m,n+1/2)

end