xnd: phases description in the input file.

Introduction

The phase block has to be repeated for each phase according to the nPha value declared in the program header. After the header the variables are encoded using the standard method seen in encoding variables in xnd.

Phase Header and common blocks

Phase Header

#optional comments
NbHkld 0classical phase, the hkl lines are generated
n > 0classical phase but the hkl lines are read
-1 reserved for amorphous sample( not yet implemanted)
nTransitdnumber of phase transitions to take into account in the scale factor
Temp_OrigdOrigin of Temp in the phase variable expansion
TitlegsTitle of the phase

Scale factor block

#optional comments
Orderdkind of block
Scale@Scale factor for the phase
B_global@Overall isotropic Thermal factor

Phase transition block

This block is read only if nTransit in the phase header has a positive value. In this case the effective scale factor is : Scale_eff = Scale * exp(u) /(exp(u)+exp(-u)) with u = signe(T-T_c) * T_a * abs(T-T_C)^expo
#optional comments
Orderdkind of block
the following group of lines is repeated nTransit times
T_c@Origin Temp
T_a@Activation Temp
Exponant@Overall isotropic Thermal factor
6.2 / # commentaire facultatif si nb_hkl >= 0 ou -4 ou -5 s nom_sym identificateur declare au &2 ou dans le fichier sym d n_orien nb fonctions pour orientation preferentielle d n_prof nb fonctions pour profil et assymetrie d nb_atom nombre d'atomes independants d nb_bloc nombre de blocs rigides d incohrent code pour diffus (reserve pour) d nb_bond nombre de temps pour calcul des contraintes # commentaire facultatif d o_maille ordre du developpement *@* centrage *@* a *@* b *@* c *@* alpha *@* beta *@* gamma 6.21/ ..... (projet) si la phase a ete declaree magnetique ou modulee on intercale a ce niveau les declarations qui permettent de construire cette phase modulee: 6.211/ # commentaire facultatif d dim nb de vecteurs de modulation d nature 0 : incommensurable (integration par defaut) >=1 : comensurable d'ordre nature <=-1 : incommensurable avec -nature pts d'integration d nb_sat nb de satellites a generer y compris la structure de base f t_org pour les struct comensurables, t de la section d'origine 6.212/ les groupes d'indices des satellites dim * n_sat valeurs d n_hkl 0 (genere) ou nb de hkl a lire d hkl_sat indices complementaires de la famille de satellites 6.213/ nom du groupe complementaire de symetrie s nom_sym identificateur declare au &2 ou dans le fichier sym 6.214/ # commentaire facultatif @ les dim vecteurs de modulation d o_vect *@* x_vect *@* y_vect *@* z_vect 6.22 / les orientations preferentielles, deux cas n_orien positif, les fonctions sont communes a touts les experiences n_orien negatif, il y a n_orien fonction par experiences, il faut donc en lire n_orien * n_manip 6.221 / les fonctions de definitions # commentaire facultatif @ si n_orien > 0 d c_orien les n_orien numero des fonctions utilises d ... @ si n_orien < 0 d c_orien les n_orien*n_manip numero des fonctions utilises d ... 6.222 / les coefficients # commentaire facultatif d o_orien ordre du developpement n_orien(*n_manip) fois le groupe suivant *@* coef coefficient de la premiere fonction *@* teta orientation de la premiere fonction *@* phi orientation de la premiere fonction *@* coef.. ... 6.23 / les fonctions de profils 6.231 / les fonctions de definitions # commentaire facultatif @ n_prof > 0 d c_prof les n_prof numero des fonctions utilises d ... 6.232 / les coefficients pour le profil # commentaire facultatif d o_prof ordre du developpement n_prof fois le groupe suivant ou chaque terme est multiplie par la valeur de la fonction *@* Wl_C largeur composante lorentz : / cos(theta) *@* Wg_C largeur composante gauss : / cos(theta) *@* WlT largeur composante lorentz : * tan(theta) *@* WgT largeur composante gauss : * tan(theta) *@* teta orientation de l'axe *@* phi orientation de l'axe 6.234 / les coefficients pour l'assymetrie # commentaire facultatif d o_prof ordre du developpement n_prof*max_assym fois le groupe suivant ou chaque terme est multiplie par la valeur de la fonction la boucle sur max_assym est interne vis a vis de n_prof *@* A0 assymetrie (constante) *@* AT assymetrie : * tan(theta) 6.24 / les atomes # commentaire facultatif @ si nb_atom > 0 c7 nom identificateur c7 espece nom longueur ou facteur dif d nature 0 : atome simple B_isotrope +1 : B_anisotrope +2 : Avec moment magnetique (reserve pour) +4, 8, 12, 16 : Avec modulation ordre 1, 2, 3, 4... # commentaire facultatif d o_atom ordre du developpement *@* X *@* Y... Z, T, puis selon nature dans l'ordre des codes de nature *@* B ou par exemple *@* B11... B22, B33, B12, B13, B23 *@* B... Kx, Ky, Kz *@* B11... B22, B33, B12, B13, B23, Kx, Ky, Kz *@* B... A_x, A_y, A_z, P_x, P_y, P_z, A_t, P_t modulation A sin (2 pi ( t + P )) 6.25 / les blocs rigides # commentaire facultatif @ les nb_blocs d o_bloc *@* x *@* y *@* z *@* teta *@* phi *@* psi 6.251/ puis pour chaque bloc d n_atom les n_atom du bloc 6.252/ # commentaire facultatif c7 nom identificateur c7 espece nom longueur ou facteur dif d nature 0 : B_isotrope 1 : B_anisotrope 2 : B_isotrope magnetique 3 : B_anisotrope magnetique # commentaire facultatif d o_atom ordre du developpement *@* X dans le repere ortho du bloc *@* Y... Z, T, puis selon nature *@* B *@* B11... B22, B33, B12, B13, B23 *@* B... K1, K2, K3 *@* B11... B22, B33, B12, B13, B23, K1, K2, K3 6.252/ bloc suivant .... 6.26/ # commentaire facultatif @ si nb_hkl > 0 et n=atom!=0 : nb_hkl groupes d o_int 0 si n_atom == 0 ou ordre des intensites d h d k d l *@* int intensites : F^2, ces valeurs ne sont pas corrigees d h k l de Lorentz-Polar et de la multiplicite : un groupe est cependant necessaire (P1 possible) on peut introduire directement les valeurs du fichier hkl il y a cependant une legere difference dans les deux modes de calcul. Cette operation peut etre automatisee avec le programme utilitaire xhkl2k ou ce bloc est identifie par les commentaires #BEGINHKL2H et #ENDHKL2K. En cas de phase modulee, on doit relire un bloc 6.26 pour chaque groupe de satellite ou n_hkl est non nul. 6.27 / ..... (projet) # commentaire facultatif @ si nb_dif > 0 : nb_dif groupes d o_int ordre d h d k d l *@* int intensites hors corrections instrumentales *@* wl largeur lorentz d h k l .... 6.28/ les distances interatomiques et leur contraintes # commentaire facultatif @ si nb_bond > 0 d nb_dist nb de parametres de distances d nb_angl nb de parametres angulaires f temps ou est calcule le systeme de contrainte @f temps nb_bond fois 6.281/ # commentaire facultatif @ nb_dist fois d d_coord f sigma ponderation du systeme de contrainte # commentaire facultatif d o_dist ordre du developpement *@* dist @ d_coord fois c7 nom_at0 distance entre At0 et sym (At1) c7 nom_at1 s% sym chaine du type 1+x,y-1/2,-z 6.282/ .... (projet) @ nb_angl fois d a_coord f sigma ponderation du systeme de contrainte d o_angl ordre du developpement *@* angl @ a_coord fois c7 nom_at0 sommet c7 nom_at1 extremite = sym (At1) s% sym1 c7 nom_at2 extremite = sym (At2) s% sym2 6.3 / si nb_hkl == -1 ( amorphe, projet non affinable) d nb_dist nombre de diffusion 6.31 / nb_dist fois d nb_pos nombre de points c7 espece origine (nom du facteur de diffusion) c7 espece extremite d o_dist ordre du developpement f dif_pos distance du point ou le fond est defini @* dif_val amplitude 6.4 / nb_hkl == -2 ( raies parasites ) d n_prof nb de fonction de profil 0 ou 1 d nb_dist nb de raies de diffraction 6.41 / @ n_prof > 0 d c_prof numero de la fonction utilisee, 0 d ... d o_prof ordre du developpement *@* Wl_C largeur composante lorentz : / cos(theta) *@* WlT largeur composante lorentz : * tan(theta) *@* Wg_C largeur composante gauss : / cos(theta) *@* WgT largeur composante gauss : * tan(theta) *@* A0 assymetrie atan : constante *@* B0 assymetrie exp : constante *@* AT assymetrie atan : * tan(theta) *@* BT assymetrie exp : * tan(theta) *@* teta orientation de l'axe *@* phi orientation de l'axe 6.31 / nb_dist groupes d o_int int *@* dist pour chaqur raie distance et *@* int intensite .... 6.4 / nb_hkl == -3 ( quasicubique ) d n_prof nb de fonction de profil 0 ou 1 d nb_dist nb de raies de diffraction d o_cell ordre du developpement *@* centrage *@* a parametre de la maille 6.41 / @ n_prof > 0 d c_prof les n_prof numero des fonctions utilises d ... d o_prof ordre du developpement *@* Wl_C largeur composante lorentz : / cos(theta) *@* WlT largeur composante lorentz : * tan(theta) *@* Wg_C largeur composante gauss : / cos(theta) *@* WgT largeur composante gauss : * tan(theta) *@* A0 assymetrie degre1 : constante *@* B0 assymetrie degre3/2 : constante *@* AT assymetrie degre1 : * tan(theta) *@* BT assymetrie degre3/2 : * tan(theta) *@* teta orientation de l'axe *@* phi orientation de l'axe 6.42 / nb_dist groupes d o_int d m d n d 0 *@* int intensites d ml ....