FullProf User's Guide

Download as txt, pdf, or txt
Download as txt, pdf, or txt
You are on page 1of 76

FullProf user's guide FullProf Rietveld, Profile Matching & Integrated Intensities Refinement of X-ray and/or Neutron Data

(powder and/or single-crystal)

HELP FULLPROF

1.- INTRODUCTION AND GENERAL INFORMATION SHORT REFERENCE GUIDE OF THE PROGRAM FFFFFF F F FFFF F u u F u u F uuuuu l l l l l lll l l l l l lll PPPPPP P P P P PPPPPP P r r r ooo P r o o P r ooo f f f fff f f

*************************************** * Program : FullProf * *************************************** (Version 3.5d Oct98-LLB-JRC) Juan Rodriguez-Carvajal Laboratoire Leon Brillouin (CEA-CNRS) Tel: (33) 1 6908 3343, Fax: (33) 1 6908 8261

Disclaimer: The author is not responsible for erroneous results obtained with FullProf. This guide cannot substitute the lack of knowledge of users on crystallography, magnetism, diffraction physics and data analysis. This short guide is merely a description of the input files with minor explanations on how to proceed. Powder diffraction is becoming more and more powerful but FullProf is not an "automatic" (black-box) program, as is usually found in single crystal structure determination. No attempt has been made in order to predict the behaviour of the program against bad input data. The user must check his(her) data before claiming a misfunction of the program. The author acknowledges all suggestions and notification of possible bugs found in the program. The most recent version of FullProf is either in "pub/divers/fullp" or in "pub/divers/fullp" of the anonymous ftp-area of the LLB unix-cluster.

Users interested in create their own subroutines to link with FULLP-library are asked to read the file "fpreadme" in the above mentioned disk-area. To access this area from Internet, one has to type in the local host the following commands: LocalPrompt> ftp charybde.saclay.cea.fr Answer with the word: anonymous , to the Login request and password. Within the ftp prompt, do: From a UNIX host: ftp>cd pub/divers/fullp ftp>get fpreadme ftp>bye From some VMS-VAX hosts: ftp>set def "pub/divers/fullp" ftp>get "fpreadme" ftp>ex -> Go to FullProf area -> Obtain the document -> Return to host -> Go to FullProf area -> Obtain the document -> Return to host

---------------------------------------------------1.1: Purpose, authors, references and documentation ---------------------------------------------------FullProf is a program for Rietveld analysis (structure profile refinement) of neutron (CW, TOF, nuclear and magnetic scattering) or X-ray powder diffraction data collected as a function of the scattering variable T (2theta or TOF). The program can be also used as a Profile Matching tool, without the knowlegde of the structure. Single Crystal refinements can also be performed alone or in combination with powder data. FullProf has been developed starting from the program of Wiles & Young, J. Applied Cryst.14,149(1981), (DBW3.2S, Versions 8711 and 8804). The modifications of the code are mainly related with the re-organization of the central routines performing the calculation of profile functions, derivatives, structure factors, and the introduction of many other things. The total source is more than 1 Megabyte (more than 28.000 fortran lines). The format of the main control input file (e.g. a control file created for use with DBW-8711 can be used by FullProf with minor modifications). The input file is accepted as "interpreted free format". The source is written in standard FORTRAN 77 language, and is organized as to be easily adapted to different computers. The actual version can be run on VAX, Alpha and Unix computers, MacIntoshes and on PCs (Lahey Computer Systems Inc. FORTRAN-compiler, minimum 386/4Mb with co-processor required). The migration towards genuine Fortran 90 is in progress. --------------------------1.2: Features of FullProf: --------------------------- Choice of line shape (Gaussian, Lorentzian, modified Lorentzians, pseudo-Voigt, Pearson-VII or Thompson-Cox-Hastings) for each phase. - Neutron (constant wavelength and TOF) and X-ray (laboratory and synchrotron sources)

le

One or two wavelengths (Ka1 + Ka2) Background refinement Multi-phase (up to 8 phases) Preferred orientation : two functions available Absorption correction for a cylinder Choice between three weighting schemes: standard least squares, maximum likelihood and unit wheights. Choice between automatic generation of hkl and/or symmetry operators and file given by user. Magnetic structure refinement (crystallographic and spherical representation of the magnetic moments). Two methods: describing the magnetic structure in the magnetic unit cell of making use of the propagation vectors using the crystallographic cell. This second method is necessary for incommesurate magnetic structures. Automatic generation of reflections for an incommensurate structure with up to 24 propagation vectors. Refinement of propagation vectors in reciprocal lattice units. h,k,l dependence FWHM for strain and size effects h,k,l dependence of shift and asymmetry for special kind of defects Profile Matching. The full profile can be fitted without prior knowledge of the structure (needs only good starting cell and profile parameters) Quantitative analysis without need of structure factor calculations. Chemical(distances) and magnetic (magnetic moments) slack constraints Resolution function (for pseudo-Voigt peak shape) may be supplied in a fi

- Structural or magnetic model could be supplied by an external subroutine for special purposes (rigid body, TLS, polymers, form factor refinements, small angle scattering of amphifilic crystals, description of incommensurate structures in real direct space, etc.) - Single crystal data or integrated intensities can be used as observations (alone or in combination with a powder profile) - Neutron (or X-rays) powder patterns can be mixed with integrated intensities of X-rays (or neutron) from single crystal or powder data. -------------------------------------------------1.3: Running the program, input and output files. -------------------------------------------------To run the program the user has to invoke the name of the executable file and press the ENTER/INTRO key. Of course the executable file must be in a directory included in the PATH or an alias should exist. For for doing sequential refinements there is a number of command files that can be used. The command files (scripts) depend on the operating system. A facility is included in FullProf for versions higher that 3.2 that allows sequential of cyclic refinements. Examples: 1: Prompt> FULLPROF 2: Prompt> FullProf , ... ----------Input files : ----------To run the program, you need at least one input file, CODFIL is the code of the control file given by the user. CODFIL.PCR : Input control file

It must be in the current directory to run the program. This file contains the title and crystallographic data and must be prepared by user with a file editor. There are two different formats for this file: the first one is free format and closely related to that of the Young & Wiles's program. The second is based on keywords and commands. (this last format is not available at present) Warning : this file is normally up-dated every time you run the program (see parameter NXT on line 3). In the first stages of a refinement, it is wise to save a copy of this file with a different name. The following files are optional. FILE.DAT : Intensity data file (unit 4) : format depends on instrument. If you do not specify the name FILE, the program takes FILE=CODFIL. Not necessary for pattern calculation modes. FILE.BAC: Background is the following first line : following lines: file (unit 12). The format of this file (as that of FILE.DAT for INSTRM=0): 2theta(initial) step 2theta(final) list of intensities in free format.

CODFILn.HKL : Set of files with the reflections corresponding to phase "n" (n is the ordinal number of a phase). These files are optional and depend on the value of the parameter IRF(n) (see below) If sequential refinements have to be done, this file is called HKLn.HKL MYRESOL.INSTRU: File describing the instrumental resolution function. Any name can be used and its content depend on the value of the paramenter IRESO (see below). GLOBAL.SHP: File providing a numerical table for calculating the or peak shape and its derivative. The peak shape should CODFIL.SHP: be given in a normalized form P(x) where the variable x is chosen to give a FWHM=1 and the area is equal to 1 => Integ{x1,x2}[P(x).dx]= 1. That allows the use of the conventional U,V,W parameters for defining the FWHM as a function of angle. The format of this file is the following: Line1: Any comment Line2: Np8, nupr, (anpr(j),j=1,nupr) Np8 = Number of points nupr= Number of different profiles anpr(j)= Angle to which profile "j" is best adapted The rest of the lines are columns with X, P(X,1), PP(X,1),P(X,2), PP(X,2),...P(X,nupr), PP(X,nupr) in free format. PP(X,j) is the derivative of P(X,j) with respect to X. The profile of a reflection situated between anpr(j) and anpr(j+1) is linearly interpolated between the profiles P(X,j) and P(X,j+1) CODFIL.COR: File with corrections for integrated intensities of profile intensities depending on the value of the variable ICORR. See below. Output files : Except for *.OUT and *.SUM, their creation depends on the value of a flag which is quoted in parenthesis. The ordinal number on the flag list is given in brackets. CODFIL.OUT : This is the main output file (unit 7) which contains

all control variables and structure parameters. CODFIL.PRF : Observed and calculated profile (unit 1) : to be fed into PLOTPOW, PLOTR, ...(if IPL2 different from zero) In the case of ICRYG=1 (Integrated intensity mode) a list of sin(theta)/lambda, Gobs, Gcalc is output after two lines of comments. CODFIL.RPA : Summary of refined parameters (unit 2) : short version of CODFIL.SUM (if JCIL=1) If the file exist the new data are APPENDED at the end. CODFIL.SYM : List of symmetry operators (unit 3) (if IPL1=JSY=1) (The last two files are necessary to run DISTAN or BONDSTR) CODFIL.SUM : Parameter list after last cycle (unit 8) : the summary of the last parameters, their standard deviations and reliability factors. An analysis of the goodness of the refinement is included at the end. CODFIL.FOU : If JFOU=1 H,K,L, Structure Factors in Cambridge format (unit 9) : to be fed into FOURTK (FOURPL) to produce Fourier maps. It corresponds to the file usually called HKLFF.DAT but you must prepare the second file CRYST.cry If JFOU=2 (List of 'observed' structure factors in SHELXS format) H,K,L, Fo, sigma(Fo) (3I4,2F8.2) JFOU=-1 or -2, as above but they are calculated in another way. The Fcalc in JFOU>0 may depend on the peak shape and the integration interval, because they are obtained by integration of the calculated profile in the same way as the 'Fobs' are obtained from 'Iobs'. If JFOU is negative, Fcalc are really the structure factors of the conventional cell in absolute units. JFOU=3 Format suitable for the program FOURIER (3I4,2F10.4,f8.5,f10.4) H,K,L,Freal,Fimag,sintheta/lambda,fobs JFOU=4 Format (3I4,2F10.4,i8) H,K,L,Fobs,Fcalc,nint(10000.* Phrd) For JFOU= 3,4: Phrd is the phase in radians and the observed and calculated structure factors of the conventional cell are in absolute units CODFILn.SHX :If JFOU=2,-2, template of SHELXS *.IN file. CODFILn.INP :If JFOU=3,-3, template of FOURIER *.INP file. CODFILn.HKL :Files that can be input or output files. Depending on the value of IRF(n)

CODFIL.INT :Single integrated intensity file when the program is used for refining with ICRYG=1,2 (see below). CODFIL.HKL : (if JLKH<>0, unit 10). Complete list of reflections of each phase. JLKH=1 --> If JOBTYP less than 2 reflection code, h, k,l, multiplicity, dspacing,twotheta, FWHM, Iobs, Icalc, Iobs-Icalc. --> If JOBTYP >1 h, k, l, multiplicity, Icalc, twotheta, dspacing JLKH=2 --> Output for SIRPOW92 h,k,l,mul,sint/l,2t,Fwhm,F2,sF2 JLKH=-2 --> Output for EXPO h,k,l,Fwhm,F2 JLKH=3,-3 --> Output of real and imaginary part of structure factors (only for crystal structures) h, k, l, mul, Freal, Fimag, 2theta, Intensity If JLKH<0 the structure factors are given for the conventional cell. Otherwise the structure factor corresponds to the non-centrosymmetric part of the primitive cell. (the obtained file can be used as a CODFILn.HKL files for new runs) JLKH=4 --> Output of: h, k, l, F2, sigmaF2. Where F2 is the "observ ed" structure factor squared. The file could be used as inpu t for a "pseudo-single" crystal integrated intensity file using ICRYG=1 and IRF=4 JLKH=5 --> Output of: h, k, l, mult, Fcalc, T, D-spacing, Q. Where Fcalc is the module of the calculated structure fa ctor. This file can be used as input for JBT=-3 and IRF=2 in o rder to perform quantitative analysis without re-calcultating the structure factors for each cycle. The Fcalc are in absolute units for the conventionnel call. CODFIL.SAV : (if JCIL=2) List of reflections between two selected angles h, k, l, multiplicity, Iobs, twotheta, dspacing For build-in sequential refinements (version 3.2 and higher) the user must prepare the data files using names of the form CODnnn.dat. COD stands for the code of these files and can be formed by whatever number of characters (compatible with the actual operating system). nnn stands for a sequence of integers. All CODnnn.dat files must be in the same directory and the numbers nnn should be in between a minimun number (first) and a maximun number (last) that are asked by program. Holes are allowed between first and last. The file CODFIL.PCR can have a different code (CODFILcould be different of COD) and it will be used for refining the whole set of CODnnn.dat files. The final results are contained in the CODFIL.RPA file.

For VaX-users using a command file to execute FullProf in cyclic mode: For sequential refinements, *.DAT files will normally be prepared by SEPFIL. In this case CODFIL must have three letters (e.g.XXX) as code followed by a number. The *.PCR file must be named XXXIN.PCR XXXCYC.RES : In the case of sequential refinements, all the above files would rapidly yield a quota exceeded error message. Thus only condensed results are saved in files XXXCYC.RES (similar to CODFIL.RPA) and XXXSUM.RES (similar to CODFIL.SUM). XXXHKL.RES : List of reflections of a selected zone of the diffraction patern (useful with Profile Matching mode) The final version of the file XXX**.PCR (where ** corresponds to the last data set is also saved.

2.- DETAILED DESCRIPTION OF INPUT FILES

- - -CODFIL.PCR This file is free format. That doesn't mean free format in FORTRAN (,*)-sense A routine interprets the items given by the user that must obey the order given below. A space is needed between each item (except when the second is a negative number). When the program is run, messages of error reading a line of this file are normally due to a previous error. For example, the number of atoms you really wrote does not correspond to the number you put in the line following the name of the phase. Empty lines as well as lines starting with the symbol "!" in the first column are considered as comments and are ignored by the program. If the user starts his(her) CODFIL.PCR file with the left-ajusted capital "COMM", the new CODFIL.PCR file has comments with mnemonics for each variable. If the user introduce his(her) own comments, they are not saved in the new version of the file. The unexperienced user can create a template by answering "starting" (without quotes) to the prompt asking for the name of the file. Note that a star after a line number (or a variable) indicates that the line's (or variable) existence depends on the value of a control variable. ============================================================================ LINE 1 : TITLE (any 70 characters to be used to label the printout) If the first four character of TITLE correspond to the word TITL the file is given in "command mode" (not available yet). If the first four characters of TITLE correspond to COMM, comments lines (starting with ! in the first column) are automatically addet to the new CODFIL.PCR (or CODFIL.NEW). The comment lines give a keyword for each variable in order to be easily recognized by the user. This comment line has been included below to ============================================================================ LINE 2 :

JOBTYP, NPROF, NPHASE, NBCKGD, NEXCRG, NSCAT, NORI, IDUM, IWGT, ILOR, IASG, IRESO, ISTEP, NRELL, ICRYG, IXUNIT, ICORR (15 integers) (It is understood that they are separated by a space) --------------Comment line : !Job Npr Nph Nba Nex Nsc Nor Dum Iwg Ilo Ias Res Ste Nre Cry Uni Cor ---------------------------JOBTYP = 0 (Job) 1 2 3 -1 -3 X-ray case Neutron case (constant wavelength, nuclear and magnetic) pattern calculation (X-ray) pattern calculation (Neutron, constant wavelength) Neutron case (T.O.F., nuclear and magnetic) pattern calculation (Neutron, T.O.F.)

If abs(JOBTYP)>1 and IDUM=1 (see below) a calculated pattern is created with the name CODFIL.SIM in format corresponding to INSTRM=0. This pattern corresponds to an "ideal observed" pattern and can be use for simulation purposes in order to investigate the effect of systematic errors on the structural parameters and on the reliability factors. NPROF = Default value for selection of a peak shape. Particular (Npr) values can be given for each phase (see line 11-2) 0 1 2 3 4 5 6 7 8 9 10 11 12 Gaussian Cauchy Modified 1 Lorentzian Modified 2 Lorentzian Tripled pseudo-Voigt pseudo-Voigt Pearson VII Thompson-Cox-Hastings pseudo-Voigt Numerical profile given in CODFIL.SHP or in GLOBAL.SHP T.O.F. Convolution pseudo-Voigt x Double Exponential Not yet used Split pseudo-Voigt function Pseudo-Voigt function convoluted with axial divergence asymmetry function (Finger, Cox & Jephcoat, J. Appl. Cryst. 27, 892, 1994)

NPHASE = number of phases ( max:8) (if NPHASE <0 the number of phases is (Nph) abs(NPHASE) and the asymmetry correction is applied following the app roximation of C.J.Howard, J.Appl.Cryst.15 615-620 (1982) with the Simpson formula f or five points) NBCKGD =0 (Nba) 1 Refine background with polynomial function Read background from file CODFIL.BAC. The format of this file is explained above. Some coefficients are read below. 2,3,.,N linear interpolation between the N given points If NBCKGD<0 but IABS(NBCKGD)>4 the interpolation is performed using cubic splines -1 refine background with Debye-like + polynomial function. -2 Background treated iteratively by using a Fourier filtering technique. An extra parameter is read below. The starting backgroung is read from

-3 NEXCRG = (Nex) NSCAT = (Nsc)

file FILE.BAC as for NBCKGD=1. Read 6 additional polynomial background coefficients number of excluded regions number of scattering sets (zero in most cases) If NSCAT>0, the program performs an internal fit if a table is given in order to get coefficients for the exponential expansion (see below). If NSCAT<0, a linear interpolation is made. preferred orientation function No 1 preferred orientation function No 2 (March)

NORI = 0 (Nor) 1

IDUM =1 If equal to 1 and some of the phases are treated (Dum) with Profile Matching modes, the criterium of convergence when shifts are lower than a fraction of standard deviations is not applied. =2 If equal to 2, the program is stopped in case of local divergence: chi2(icycle+1) > chi2(icycle) =3 If equal to 3 the reflections near excluded regions (excl+/-wdth*2theta) are not taken into account to calculate the Bragg R-factor. These reflections are omitted in the output files with hkl's. If JOBTYP greater than 1 and IDUM is different than zero a file CODFIL.SIM is generated IWGT =0 standard least squares refinement (Iwg) 1 maximum likelihood refinement 2 unit weights ILOR =0 Standard Debye-Scherrer geometry, or Bragg-Brentano if (Ilo) the iluminated area does not exceed the sample surface. If Bragg-Brentano geometry is used but the above condition is not fulfilled, the intensity data must be corrected for the geometric effect before attempting any refinement. (A partial correction can be done by using the parameter SENT0 in line 5) =1 Flat plate PSD geometry =2 Transmission geometry. Flat plate with the scattering vector within the plate (Stoe geometry for X-rays) =3 Polarization correction is applied even if the format of the DATFIL.DAT file does not correspond to one of the synchrotron explicitely given formats (see below). This must be used for synchrotron data given in a X,Y,Sigma format (INSTRM=10). IASG =1 Subroutine ASSIGN is called at each cycle, then reflections (Ias) are re-ordered. =0 Subroutine ASSIGN is called only at the first cycle (If JBT=2 for one phase, IASG must be =1) IRESO=0 Resolution function of the instrument is not given (Res) If IRESO is not zero, the next line contains the name of the file where the instrumental resolution function is given. The profile is assumed to be a Voigt function (NPROF=7). 12 parameters or a table determine the resolution function. Ui,Vi,Wi,Xi,Yi,Zi (i=1,2 for lambda1 and lambda2) The different types of functions are:

=1 HG**2= (Ui*tan(q)+Vi)*tan(q)+Wi HL= Xi*tan(q)+Zi =2 HG**2= (Ui*tan(q)+Vi)*tan(q)+Wi HL= (Xi*(2q)+Yi)*(2q)+Zi =3 HG**2= (Ui*(2q)+Vi)*(2q)+Wi HL= (Xi*(2q)+Yi)*(2q)+Zi =4 List of values 2q, HG(2q), HL(2q) (a linear interpolation is applied for intermediate 2q) ISTEP=1,2,3,.. If ISTEP>1 the number of data points is reduced by (Ste) a factor of ISTEP. Only those points corresponding to the new step size ISTEP*STEP (see Line #3 below) are taken into account in the refinement. Useful for speed-up preliminary refinements. NRELL (Nre) Number of parameters to be constrained within given limits. At the end of the file you must give a list of NRELL lines specifying the number and the limit of each parameter. If not equal to zero, only integrated intensity data will be given. No profile parameters are needed. For ICRYG=2 no least-squares algorithm is applied. Instead a Montecarlo search of the starting configuration is performed. A selected number of parameters NRELL are moved within a box defined by the NRELL relations fixing the allowed values of the parameters. The best (lowest R-factor) NSOLU solutions are printed and the CODFIL.PCR file is updated with the best solution. (See NRELL variable in this line and Line ) IXUNIT Units of the scattering variable (Uni) =0 2theta in degrees =1 T.O.F. in micro-seconds ICORR (Cor) =0 No correction is applied =1 A file with intensity corrections is read. The corrections are applied to the integrated intensities as a multiplicative constant. The file CODFIL.COR starts with a comment and follows with a list of pairs: a simple list of abcisae and correction values. TITLE ...... Scattering variable (T) Value of the correction " " ............. ................ Data are read in free format. For peaks between points provided in the CODFIL.COR file, the correction is linearly interpolated. Example: First line -> This is my correction FILE for Following lines -> 10.0 1.3 20.0 1.1 30.0 1.0 40.0 0.9 80.0 0.8 120.0 0.7 180.0 0.7

ICRYG (Cry)

The intensity of a reflection at scattering variable 40 is assumed to be I(calc)*0.9. =2 A similar file is read but the coefficients of an empirical function and their standard deviations are read instead of directly the corrections. The format is: First line -> TITLE .... Second line -> ITYCORR, ITYFUNC, NPCORR Following lines -> Coefficient Sigma(Coefficient) (NPCORR lines) If ITYCORR = 1 corrections are applied to the integrated intensities. Standard deviations must not be given. ITYCORR = 2 corrections are applied to the observed profile. The corrected observed profile and their variance are obtained as: y(corr) = y(obs)/cor Sigma2(y(corr)) = sigma2(y(obs))/cor^2 + sigma2(cor)/y(obs)^2 NPCORR : Number of coefficients of the empirical function. ITYFUNC=1 ITYFUNC=2 Polynomial function: cor = Sum{i=1,npcorr}{coeff(i)* T**(i-1))} Exponential + Maxwellian for TOF raw data

cor = Coeff(1)+ Coeff(2)*Exp(-Coeff(3)/T^2)/T^5+ + Sum{i=4,NPCORR,2}{Coeff(i)*Exp(-Coeff(i+1)*T^2)} Line 2-1*: FILERES (A16) Name of the file with the instrumental resolution function. To be given only in the case of IRESO<>0. The items in FILERES are read in free format. The first line is considered as a title For IRESO=1,2,3 the 12 parameters Ui, Vi, Wi, Xi, Yi, and Zi are read from lines 2 and 3 (see the above line for the available instrumental functions). Example: Line1: Resolution function of MyXrayDiffractometer Line2: 0.00802 -0.00936 0.01024 0.0029 0.0 0.0 ! U1,V1... Line3: 0.00774 -0.00552 0.00814 0.0000 0.0 0.0 ! U2,V2... For IRESO=4, the file FILERES starts with a line whith the title followed by a line with the number of points (NPOINS) where the instrumental Gaussian and Lorentzian FWHM are given. NPOINS lines follow containing the three items: 2thet, HG and HL. The Bragg peaks of the diffraction pattern must be between 2thet(1) and 2thet(NPOINS). For this case the same resolution function is applied to both wavelengths. The maximum number of NPOINS is 30. ============================================================================ LINE 3 :IOT, IPL, IPC, MAT, NXT, LST1, LST2, LST3, IPL1, IPL2, INSTRM, JCIL, JSY, JLKH, JFOU, ISHOR, IANALY (17 integers)

--------------Comment line : !Ipr Ppl Ioc Mat Pcr Ls1 Ls2 Ls3 Syo Prf Ins Rpa Sym Hkl Fou Sho Ana ---------------------------List of output control flags : normally 0 = off / any value = on IOT = 1 obs. & calc. profile intensities --> CODFIL.OUT (0) (Ipr) 2 The files CODFILn.SUB with the calculated profile of each phase are generated. 3 As 2 but the background is added to each profile. IPL = 1 line printer plot --> CODFIL.OUT (Ppl) 2 Generates the background-file FILE.BAC 3 Puts difference pattern in file FILE.BAC IPC = 1 list of obs. & calc. integr. int. --> CODFIL.OUT (Ioc) = 2 The reflections corresponding to the second wavelength are also writen if different from the first one. MAT = (Mat) correlation matrix --> CODFIL.OUT If MAT=2, the diagonal of LSQ matrix is printed before inversion at every cycle.

NXT = 1 CODFIL.PCR is re-written with updated parameters (Pcr) 2 new input file --> CODFIL.NEW LST1 = reflection list --> CODFIL.OUT (usually 0) (Ls1) LST2 = 1 corrected data list --> CODFIL.OUT (usually 0) (Ls2) LST3 (Ls3) IPL1 (Syo) IPL2 (Prf) = = 4 In some versions of FullProf a plot of the diffraction pattern is diplayed on the screen at each cycle of refinement. merged reflection list --> CODFIL.OUT (usually 0)

symmetry operators --> CODFIL.OUT (+ CODFIL.SYM if JSY=1) plot --> CODFIL.PRF for PLOTPOW, BENSTRAP,PLOTR, etc.. " IGOR (MacIntosh software) " KaleidaGraph (MacIntosh software) & PLOTR (Pc software) " Picsure, Xvgr(Sun-Unix Software)

= output data for = 1 Format suitable = 2 " " = 3 " " = 4 " "

INSTRM= 0 Data supplied in free format (Ins) Up to seven comment lines are accepted. The first three real numbers found at the beginning of a line are interpreted as Ti, step and Tf. The following lines after (Ti, step, Tf) must contain NPTS=(Tf-Ti)/step+1 values of the intensity profile. Data format from Argonne are also interpreted by this value of INSTRM. 1 D1A/D2B format (original Rietveld-Hewat format : the first line must be 2Thetai, step, 2Thetaf, i.e. the first four lines of the POWDER file must be removed. Note however that angles are given in degrees 2Theta, not in hundredths of degree ! ).

2 D1B old format (DEC-10) 3 new format for D1B & D20 (Vax DataBase) +/-4 Brookhaven synchrotron. 4: First line: 2thetamin, step, 2thetamax (free format) Rest of file: pairs of lines with 10 items like Y1 Y2 ......... Y10 -- (10F8) intensities S1 S2 ......... S10 -" sigmas -4: Format given by DBWS program for synchrotron data. (Version DBW3.2S-8711) 5 Data from GENERAL FORMAT for TWO AXIS instrument 3 lines of text followed by two lines with the items: -> NPTS, TSample, Tregul, Ivari, Rmon1, Rmon2 -> Ti, step, Tf Set of lines containing 10 items corresponding to the Intensities in format 10F8.1, up to NPTS points (NPTS=(Tf-Ti)/step+1), followed by the corresponding sigmas in format (10f8.2) if Ivari=1. If Ivari=0 the sigmas are calculated as SQRT(Yi*Rmon1/Rmon2). 6 D1A/D2B standard format for files MYFILE.SUM prepared by D1A(D2B)SUM or equivalent programs. The extension of the data file must be 'dat'. 7 Files from D4 or D20L 8 Data from DMC at Wurenlingen (Paul Scherrer Institut) 9 Data of file CODFIL.UXD generated by the Socabim software on X-Rays diffractometer. 10 X,Y,Sigma format with header lines. In all cases the first 6 lines are considered as comments. If in the first line (left ajusted) appears the keyword XYDATA, then the following 5 lines are considered as the heading of the file. Among these 5 lines the following keywords and values have a meaning to the program: -> INTER fac_x fac_y Interpol Stepin -> TEMP tsamp fac_x internal multiplier of X-values fac_y internal multiplier of Y and Sigma-values Interpol=0 Variable step is used in the program =1 The variable step data are interpolated internally to the constant step Stepin. =2 Data are supplied directly at constant step If no sigma values are provided the program assumes that sigma(Y)=sqrt(Y). You can add comments to the data file if they start with the character ! in the first position of the line. These lines are ignored by the program. 11 Data from variable time X-ray data collection The first four lines are considered as comments The following lines are:

-> 2Thetai, step, 2Thetaf Comment -> (Time, Intensity) in format 5(F6, I10) The program uses the information contained in Time to normalize the observed intensities to the average time and to calculate the variance of the normalized values. 12 The input data file conforms to GSAS standard data file. BINTYP = LOG6, TIME_MAP and LPSD are not yet available JCIL = 1 prepares file CODFIL.RPA (Rpa) If the file exists before running the program the new data are APPENDED. 2 prepares file CODFIL.SAV (sequential refinements) JSY = 1 (Sym) JLKH (Hkl) = 1 = 2 =-2 = 3 prepares CODFIL.SYM (if 1, IPL1 must be set to 1)

prepares CODFIL.HKL as explained above Output for SIRPOW.92 Output for EXPO Output of Real & Imaginary parts of Structure Factors = 4 Output of h, k, l, F2, sF2 (to be re-used with IRF=4) = 5 Output of h, k, l, mult, Fcalc, T, D-spacing, Q. (to be re-used with JBT=-3 and IRF=2)

JFOU prepares CODFIL.FOU (Fou) = 1 Cambridge format = 2 Shelxs format (Prepares also the file CODFILn.SHX) = 3 FOURIER format (Prepares also the file CODFILn.INP) ISHOR = 1 Supress the output from each cycle. Only the information (Sho) from the last cycle is printed. IANALY= 1 Provides an analysis of the refinement at the end of (Ana) the summary file CODFIL.SUM. =2 Prints also the actual dimension of arrays. ============================================================================ LINE 4* : LAMDA1, LANDA2, RATIO, BKPOS, WDT, CTHM, TMR, RLIM, K (9 reals) (not needed if ICRYG is not 0) --------------Comment line : ! lambda1 Lambda2 Ratio BKPOS ---------------------------OR (for T.O.F. data) BKPOS, WDT, IABSCOR (2 real + 1 integer) --------------Comment line : ! BKPOS Wdt Iabscor ---------------------------LAMDA1 = wavelength l1 Wdt Cthm muR AsyLim Rpolarz

LAMDA2

= wavelength l2 (= l1 for monochromatic beam)

ABS(RATIO)=Intensity ratio I(l2)/I(l1) If RATIO<0 the parameters U, V, W (see below) for the second wavelength are read separately. BKPOS WDT CTHM TMR (muR) = origin of polynomial for background (in deg. 2theta or usecs) = width (range) of calc. profile in units of Hk (typically 3.5 for Gaussian and 10 for Lorentzian, 3-3.5 for T.O.F.) = coefficient for monochromator polarization correction (see Mathematical information) = absorption correction coefficient mR, used only for refinement on cylindrical samples and flat samples with symmetrical theta-2theta scanning (the scattering vector lying within the sample plane). m= effective absorption coefficient R= radius or thickness of the sample.

RLIM = peaks below this 2Theta limit are corrected for asymmetry (AsyLim) (see below for more details) K = polarization factor (synchrotron) (Rpolarz) fraction of mosaic-crystal (transmission geometry)

IABSCOR = Type of absorption correction for T.O.F. data (Iabscor) 1: Flat plate perpendicular to the incident beam 2: Cylindrical sample 3: Exponential correction Abs= exp(-c*Lambda**2) ============================================================================ LINE 5 : MCYCLE, EPS,RELAX1,RELAX2,RELAX3,RELAX4, THMIN, STEP,THMAX (required for pattern calculation mode only), ALPSD (required for flat plate PSD geometry only) SENT0 (required for Bragg-Brentano X-ray ILOR=0, when the cross section of the sample is lower than the beam dimensions at low angles) (1 integer and 10 reals) --------------Comment line : !NCY Eps R_at R_an R_pr R_gl ---------------------------MCYCLE = EPS = RELAX = 1 -> Thmin Step Thmax PSD Sent0

number of cycles of refinement forced termination when shifts less than EPS*e.s.d. the four relaxation factors for shifts : Atomic parameters: coordinates, magnetic moments, site occupancies & isotropic displacement (temperature) factors 2 -> anisotropic displacement (temperature) factors 3 -> profile parameters, asymmetry, overall displacement (temperature), cell constants, preferred orientation parameters, strains, size, propagation vectors & user-supplied param eters. 4 -> global parameters, zero-shift T0, background, displacement and transparency. THMIN = starting angle for calculated pattern in degrees T STEP = step size in degrees T THMAX = ending angle for calculated pattern in degrees T

ALPSD = SENT0 =

incident beam angle at sample surface in degrees Theta angle at which the sample intercepts completely the x-ray beam. Below SENT0 part of the beam doesn't touch the sample and the intensity of reflections below SENT0 have to be multiplied by the factor: sclow=sin(theta)/sin(SENT0)

============================================================================ LINE 6* : (2 reals) (not needed if ICRYG is not 0) --------------Comment line : ! 2Theta/TOF Background ---------------------------if NBCKGD > 2 or NBCKGD less than -3, there are iabs(NBCKGD) lines with Pos = position in degrees T Bck = background counts at this position If NBCKGD is positive: linear interpolation If NBCKGD is negative: cubic splines interpolation ============================================================================ - LINE 7* : (2 reals) (not needed if ICRYG is not 0) --------------Comment line : ! Excluded regions (LowT HighT) ---------------------------if NEXCRG > 0, enter limits of excluded regions : ALOW = low scattering variable bound in degrees or microsecs AHIGH = high scattering variable bound in degrees or microsecs ============================================================================ LINE 8* : 2*NSCAT sets of lines (needed only if you wish to enter your own scattering length or form factor instead of using the values stored in internal table; scattering factors and anomalous dispersion corrections incorporated in the program. --------------Comment line : ! Additional scattering factors ---------------------------Line 8-1 : NAM, DFP, DFPP,ITYM (A4, 2 reals and 1 integer) NAM = symbol identifying this set (left justified) This symbol is converted to lower case for X-ray diffraction global data. DFP = Df' or neutron scattering length b DFPP = Df" (ignored in the neutron case) ITYM =1 Indicates that you are giving a magnetic form factor If ITY=0 and and JOBTYP=1 or 3 (neutron case) the next line must not be given. You are just giving the Fermi length of the species NAM in DFP. =2 Indicates that you are giving just Df' and Df" and the program will use tabulated coefficients for the sin(Theta)/l dependent part of f (X-rays). The name NAM must correspond in this case to a

valid tabulated name (See Notes(1, 2) below). At variance with the name used for determining the scattering factor in the description of atoms, the chemical symbol used in NAM must be LOWER CASE. This is the most simple way of giving anomalous dispersion parameters for synchrotron data. Line 8-2 : (9, 7 or 2 reals -see below-) One line of the form A1,B1,A2,B2,A3,B3,A4,B4,C giving the coefficients for the analytic approximation to the X-ray form factor f. or one line of the form A,a,B,b,C,c,D giving the coefficients for the analytic approximation to the magnetic form factor f (P.J. Brown, Vol C new ed. ITC) or a set of lines of the form : sin(Theta)/l - f The set is terminated by a line with -100 in first position. If the first form is desired, A2 must not be zero. Note(1): Scattering length, X-rays and magnetic form factors are stored in internal tables. To use them you must give the "name" of the scatterer using UPPER CASE chemical symbols (scattering length), chemical species (e.g. CU+2, for X-rays) or M followed by the chemical symbol and formal charge state (e.g. MNI2, for magnetic form factor of Ni+2). These names are given in lines 11-4 behind the atom name (see below). In the case of giving user supplied Df' and Df'' the chemical symbol is converted to LOWER CASE. For X-ray diffraction the form factors symbols behind the atom name could be given either in LOWER or UPPER case. If the magnetic form factors of the rare earths are to be used, two options exist. Example: MHO3: magnetic form factor of Ho+3 as JHO3: magnetic form factor of Ho+3 as +c2 where c2 has been calculated using the dipolar approximation. Seven coefficients A,a,B,b,C,c,D are used for approximating +c2. Note(2): If a table is supplied and NSCAT>0 the program performs an internal fit to NINE coefficients and this could fails. If you want a linear interpolation NSCAT must be negative and the list is given as: 1000.0*2sin(Theta)/L- f ============================================================================ LINE 9 : MAXS = number of parameters varied (1 integer) --------------Comment line : !Number of refined parameters (appears in the same line as MAXS) ---------------------------============================================================================ LINE 10* : Global parameters (not needed if ICRYG is not 0) Line 10-1 : ZER, FLGZER, SYCOS, FLCOS, SYSIN, FLSYN, LAMBDA, FLAMBDA, IGLMORE (8 reals and 1 integer) --------------Comment line : ! Zero Code Sycos Code Sysin Code Lambda Code MORE ----------------------------

ZER FLGZER SYCOS FLCOS SYSIN FLSIN LAMBDA FLAMBDA IGLMORE

= zero point for T (in degrees) : Ttrue = Texp.-ZER Note that the shift convention is opposite to that used in Pawley's program. = codeword for zeroshift (codewords are described in the Mathematical section below). = systematic 2Theta shift with cosTheta dependence sample displacement (Theta - 2Theta diffractometers) = codeword for SYCOS = systematic 2Theta shift with sin2Theta dependence sample transparency coefficient = codeword for SYSIN = Wavelength to be refined (only 1-wavelength can be refined) = Codeword for LAMBDA. Cell parameters should be fixed if wavelength is to be refined. if different from zero the following line is read

Line 10-1-1: PO, CPO, Cp, CCp, Tau, CTau (6 reals) --------------Comment line : ! Microabsorption coefficients ! P0 Cod_P0 Cp Cod_Cp ----------------------------

Tau Cod_Tau

Microabsorption coefficients and codes. See mathematical section. (only used if ILOR=0 and JOBTYP=0 or 2) Line 10-1 (bis)* : ZERO, FZERO, DTT1, FDTT1, DTT2, FDTT2, TOFTET (Replaces Line 10-1 for T.O.F. data) --------------Comment line : ! Zero Code Dtt1 Code Dtt2 Code 2sinTh ------------------------------------------------------------------ZERO : Zero point for T (in microseconds) : Ttrue = Texp.-ZER FZERO : Codeword for zeroshift DTT1,DTT2 : The TOF position of a reflection with d-spacing d is calculated usi ng the formula T = ZER + (DTT1 + DTT2 * d) * d FDTT1,FDTT2: Codewords for DTT1,DTT2 TOFTET : value of 2sin(Theta) for the detector bank Used for obtaining the wavelengths and for Lorentz factor. Line 10-2*: BACK1, BACK2, BACK3, BACK4, BACK5, BACK6 FBACK1,FBACK2,FBACK3,FBACK4,FBACK5,FBACK6 (6 reals / 6 reals) If NBCKGD=-3, two lines more with coefficients: BACK7, BACK8, BACK9, BACK10, BACK11, BACK12 FBACK7,FBACK8,FBACK9,FBACK10,FBACK11,FBACK12

--------------Comment line : ! Background coefficients/codes ---------------------------BACK = background coefficients (see Mathematical section) FBACK = codewords for background coefficients If NBCKGD=1 (background read from file), BACK1 cannot be zero Only four coefficients are needed if such a case. The comment

line in this case is: --------------Comment line : ! Background Tranf_coefficients/codes ---------------------------Line 10-3* : BACKs,FBACKs (Only if NBCKGD = -1) (6 reals / 6 reals / 6 reals / 6 reals) --------------Comment line : ! Additional background coefficients/codes ---------------------------Four lines (see Mathematical section): Bc1, Bc2, Bc3, Bc4, Bc5, Bc6 CBc1, CBc2, CBc3, CBc4, CBc5, CBc6 d1, d2, d3, d4, d5, d6 Cd1, Cd2, Cd3, Cd4, Cd5, Cd6 Line 10-4* : FWINDOW (Only if NBCKGD = -2)

--------------Comment line : Window for Fourier filtering (appears in the same line as FWINDOW) ---------------------------Window for Fourier filtering. The value of FWINDOW must be much greater than the number of points subtended by the base of a single Bragg reflections in the widest region (a factor greater than five, at least!). The starting background is read from file FILE.BAC as in the case NBCKGD=1. But, at variance with the case NBCKGD=1, the file FILE.BAC is re-written at the end of the session. ============================================================================ LINE 11 : NPHASE sets of lines ---------------------------------------------------------------------------Line 11-1 : PHSNM = name of phase (A70) --------------Comment line : ! Data for PHASE number: n ==> Current R_Bragg: Rb ------------------------------------------------------------------------------------------------------Line 11-2 : N, NDIST, NMAGC, PREF(1), PREF(2), PREF(3), JBT, IRF, ISYM, ISTR, IFURT, ATZ, NVK, NPRO, IMORE (3 integers, 3 reals, 5 integers, 1 real and 3 integers) --------------Comment line : !Nat Dis Mom Pr1 Pr2 Pr3 Jbt Irf Isy Str Furth ---------------------------N = ATZ Nvk Npr More

number of atoms in asymmetric unit The total number of atoms for all phases cannot be greater than NATS (defined in PARAMETER statement of FUL0.INC) number of distance constraints

NDIST

NMAGC

number of magnetic moment constraints preferred orientation direction (in reciprocal space)

PREF(1,2,3)= JBT

= 0 The phase is treated with the Rietveld Method, then refining a given structural model. = 1 The phase is treated with the Rietveld Method and it is considered as pure magnetic. Only magnetic atoms are required. In order to obtain the correct values of the magnetic moments the scale factor and structural parameters must be constrained to have the same values (except a multiplicative factor defined by the user) that their crystallographic counterpart.(See note on magnetic refinements) The three extra parameters characterizing the atomic magnetic moments corresponds to components (in Bohr magnetons) along the crystallographic axes. =-1 As 1 but the three extra parameters characterizing the atomic magnetic moments corresponds to the value of M (in Bohr magnetons) the spherical Phi angle with X axis and the spherical Theta angle with Z axis. This mode works only if the Z axis is perpendicular to the XY plane. (for monoclinic space groups the the Laue Class 1 1 2/m) is required). = 2 Profile Matching mode with constant scale factor =-2 As 2 but instead of intensity the modulus of the structure factor is given in the CODFILn.HKL file = 3 Profile Matching mode with constant relative intensities for the current phase, but refinable scale factor.In this case IRF must be equal to 2. =-3 As 3 but instead of intensity the modulus of the structure factor in absolute units (effective number of electrons for X-rays/ units of 10(-12) cm for neutrons) is given in the CODFILn.HKL file. This structure factor is given for the non-centrosymetric part of the primitive cell, so for a centrosymmetric space group with a centred lattice the structure factor to be read is: Freduced = Fconventional / (Nlat*Icen) where Nlat is the multiplicity of the conventional cell and Icen=1 for non-centrosymmetric space groups and Icen=2 for centrosymmetric space groups. = 4 The intensities of nuclear reflections are calculated from a routine, supplied by the user, called STRMOD. The default subroutine handles Rigid body groups. = 5 The intensities of magnetic reflections are calculated from a routine, supplied by the user, called MAGMOD. =+10/-10 The phase can contain nuclear and magnetic contributions STFAC is called for reflections with no propagation vector associated and CALMAG is called for satellite reflections. CALMAG is also called for fundamental reflections if there is no propagation vector given but the number of magnetic symmetry matrices is greater than 0. The negative value indicates spherical components for magnetic parameters. For this case the atom parameters are input in a slightly different way. = 0 The list of reflections for this phase is automatically generated from the space group symbol = 1 The list h,k,l,Mult is read from file CODFILn.HKL (where n

IRF

is the ordinal number of the current phase) = 2 The list h,k,l,Mult,Intensity (or Structure Factor if JBT=-3) is read from file CODFILn.HKL =-1 The satellite reflections are generated automatically from the given space group symbol = 3 The list h,k,l,Mult,Freal,Fimag is read from file CODFILn.HKL In this case, the structure factor read is added to that calculated from the supplied atoms. This is useful for simplifying the calculation of structure factors for intercalated compounds (rigid host). =4,-4 A list of integrated intensities is given as observations for the current phase (In the case of ICRYG<>0 this is mandatory) The file CODFILn.HKL can also be named as HKLn.HKL, or CODFIL.INT in the case ICRYG<>0. The format of CODFILn.HKL files is the following: For abs(IRF)<4: The first two lines are read as titles (characters) The rest of the lines consist on: 1) No propagation vectors h h h h k k k k l l l l m (IRF=1) m Coeff (IRF=1+JSOL=1) m Intensity (or F) (IRF=2) m Freal Fimag (IRF=3) (free format) ( " ) ( " ) ( " )

2) NVK propagation vectors In the third line you have to give the number of propagation vectors in format (32x,i2), then you give NVK lines with: Nv K1 K2 K3 , where Nv is the ordinal number of K and Ki are the components of K in free format. h h h h k k k k l l l l nv nv nv nv m (IRF=1) m Coeff (IRF=1+JSOL=1) m Intensity (or F) (IRF=2) m Freal Fimag (IRF=3) (free format) ( " ) ( " ) ( " )

Note: The generated files when JBT=2,3 may content additional items that are not used by Fullprof. These items (sigma,angle,FWHM) can be used by other programs. The case IRF=1+JSOL=1 is to be used when shifts of Bragg reflections are observed and a model for it is known. The user must provide the value of the coefficient COEFF for each reflection. For abs(IRF)=4: - The first line is considered as a TITLE - In the second line the format of the intensity data to be read below is given. Example: (3i4,2f10.2,i4,3f8.4) (don't forget parentheses)

- R_lambda(n),Itypdata,ipow(n) (free format) R_lambda(n) :wavelength for phase n Itypdata = 0 Square of structure factors (F2) and sigma(F2) are input. = 1 Structure factors (F) and Sigma(F) are input. These quantities are transformed internally to case Itypdata=0. Ipow(n) = 0 Single crystal observations. = 1 Twinned single crystal observations. Up to 6 hkl's can contribute to a single observation. = 2 Powder integrated intensities. In this case cluster of peaks can be given. For this case Itypdata is irrelevant. -* cmono, rkks (to be given only for X-rays and Ipow=2) Correspond to variables CTHM and K for monochromator polarization correction.

1) No propagation vectors h k l Gobs Sigma(Gobs) Icode c1 c2 c3

2) Propagation vectors -Nok (number of propagation vectors: must be equal to NVK) NVK lines with: Nv K1 K2 K3 , where Nv is the ordinal number of K and Ki are the components of K in free format. h k l nv Gobs Sigma(Gobs) Icod c1 c2 c3

The format of the data corresponds to that given explicitely in line 2 of the CODFILn.HKL file. No data reduction is performed. The program expects to be provided with an independent set of reflections. nv is the ordinal number of the propagation vector corresponding to the current observation (hkl). Gobs and Sigma(Gobs) have different meanings depending on the value of Itypdata and Ipow. Ipow=0 Itypdata= 0 Gobs=F2 = 1 Gobs=F Ipow=1 As above but is Gobs<0 the reflection contributes to the next positive observation. Ipow=2 Gobs= Sums of {jLpF2} Icod: code for reflections indicating the scale factor number to be applied (for twinned crystals or inhomogeneous data). If Ipow=2 and NVK><0 Icod is the multiplicity. If Ipow=2 and NVK=0 the multiplicity is automatically calculated from the symbol of the space group. For IRF=4 c1,c2,c3: Not yet used (coefficients for extinction corrections) For IRF=-4 c1,c2: Real and Imaginary part of the partial calculated structure factor or the reflection. The program will add this

contribution to the structure factor calculated with the given atoms -> Ftot= F + Fp= (A+iB) + (c1+ic2). See comments for IRF=3. Examples: Twinned Orthorhombic crystal with two domains (a,b,c) (b,a,c) h k l Gobs Sigma Icod ............................. 2 0 0 -1.0 0.0 2 0 2 0 3221.0 12.1 1 3 1 1 -1.0 0.0 2 1 3 1 1221.0 8.2 1 Powder cluster of peaks 1 1 1 23.2 0.4 1 ............................. 5 3 1 -1.0 0.0 1 3 4 2 -1.0 0.0 1 4 4 1 -1.0 0.0 1 5 0 3 832.1 9.4 1 ISYM=0 Isolated peak Cluster of peaks: four independent reflections contribute to <- this observation

The symmetry operators are generated automatically from the space group symbol. =+/-1 The symmetry operators are read below. In the case of a pure magnetic phase ISYM must be always equal to 1. For JBT=10 with magnetic contribution ISYM could be 0 but a comment starting with "Mag" should be given after the space group symbol (see below) Note: For Profile Matching mode 2, IRF can be 0 in the first run. In that case, a CODFILn.HKL file is generated and IRF is set to 2 in the new CODFIL.PCR file. The file is updated at each run in the case of JBT=2. Of course ISYM must be 0. If for a phase IRF.LE.0 and ISYM=1, the reflections are generated from the symbol given in the place reserved for the space group. In that case, a file CODFILn.HKL is generated with the relevant (non-zero) reflections and proper multiplicities for the particular model described by user-given symmetry operators. In addition the calculated intensities are given in F2 (corrected for multiplicity, scale and LP-factor) in absolute units. The program doesn't use the intensities in new runs reading this generated file. The contain of this generated file, apart from the features described above, is:

No k-vectors -> h k-vectors -> h

k k

l l

nv

m F2(calc) F2(obs) m F2(calc) F2(obs)

hr

kr

lr

with obvious meaning. ISTR =0 =1 If strain or/and size parameters are used, they are those corresponding to selected models (see below) The generalized formulation of strains parameters will

be used for this phase. See below =2 =3 IFURT The generalized formulation of size parameters will be used for this phase. See below The generalized formulation of strain+size parameters will be used for this phase. See below Number of further parameters defined by user, to be used with user supplied subroutines.

AZT = Z.Mw.f^2/t (useful to calculate the weight percentage of the phase) Z: Number of formula units per cell, Mw=molecular wheight f: Used to transform the site multiplicities used on line 11-41 to their true values. For a stoichiometric phase f=1 if these multiplicities are calculated by dividing the Wyckoff multiplicity m of the site by the general multiplicity M. Otherwise f=Occ.M/m, where Occ. is the occupation number given in line 11-4-1. t: Is the Brindley coefficient that accounts for microabsorption effects. It is required for quantitative phase analysis only. When phases have like absorption (in most neutron uses), this factor is nearly 1. If IMORE=1 (see below) the Brindley-coeff. is directly read in the next line (in such case ATZ=Z.Mw.f^2). NVK = Number of propagation vectors. If NVK<0 the vector -K is added to the list. NPRO = Integer indicating the peak shape for the present phase (see line 2). If NPRO=0, the default value NPROF is taken. IMORE= If different from 0 a new line is read ---------------------------------------------------------------------------Line 11-2-1*: JVIEW, JDIST, JHELIX, JSOL, JMOM, JTER, BRIND, RM1, RM2, RM3, JTYP (6 integers, 4 reals and 1 integer) (Read only if IMORE=1) --------------Comment line : !Jvi Jdi Hel Sol Mom Ter Brind ---------------------------RMua RMub RMuc Jtyp

JVIEW = 1 a file suitable for SCHAKAL is generated = 2 a file suitable for STRUPLO is generated (The extension of the file is in both cases ".sch") =11 If JBT=2 a file CODFILn.INT with a list of overlaped peak clusters is output. JDIST = 1 Creates a file called CODFILn.ATM with all atoms within a primitive unit cell for a magnetic phase. The number "n" corresponds to the number of the current phase.If JBT=10 o nly the list of magnetic atoms is generated. -1 For a magnetic phase creates a file called CODFILn.ATM with a format suitable for further processing with the program MOMENT.

2 As 1 but for a crystal structure, all atoms inside the conventional cell are generated. JHELIX = 1 The real and imaginary components of the Fourier coefficient of a magnetic atom are constrained to be orthogonal. The factor 1/2 is also included (see mathematical section). JSOL JMOM JTER BRIND RM1 = 1 Additional hkl-dependent asymmetry and shifts parameters are read. - unused at present - unused at present - Brindley coefficient (see line 11-2) - Used when IRF=4 and IMORE=1. If RM1=0.0 the program makes RM1=1.0 internally. The meaning RM1 correponds to the global weight of the integrated intensity observations with respect to the global profile. The contribution to the normal equations of the integrated intensity part is multiplied by RM1. - If IRF=4, RM2 is a factor for excluding reflections only the reflections with Gobs>= RM2*Sigma(Gobs) are considered in the refinement. If JVIEW=11 and JBT=2 and IRF<>4 see note below. - If RM3>0.9 the weights are divided by the Chi2 of the precedent cycle (not tested!) for integrated intensity refinements (IRF=4). If JVIEW=11 and JBT=2 and IRF<>4 see note below. - Job type for the phase. Allows the refinement of hetero geneous data (Same values as the global variable JOBTYP in line 2). For the moment is only useful for IRF=4.

RM2

RM3

JTYP

Note:If JVIEW=11 and JBT=2 the parameters RM2 and RM3 are used to control whether two consecutive reflections belongs to a same cluster. This is only for IRF different from 4/-4. The rule is the following: The reflections i and i+1 belong to the same cluster if a) TwTh(i+1)-TwTh(i) < (Fwhm(i)+Fwhm(i+1))*RM2/2 OR b) TwTh(i+1)-TwTh(i) < (Fwhm(i)+Fwhm(i+1))/2 and I(i+1)<Isum*RM3 Isum being the cumulated integrated intensity of the current cluster. If the RM2 and RM3 are given as zeroes, the program uses the values RM2=1.0 and RM3=0.2. ---------------------------------------------------------------------------Line 11-3 : SYMB, Comment (A20,A60) SYMB: Space group symbol e.g. P 63/m for P63/m P 21 21 21 for P212121 Comment: Only needed for JBT=+10/-10 (See below) --------------Comment line : <-- Space group symbol (appears as default Comment) ---------------------------Note that rhombohedral space groups must be given in

the hexagonal description. Warning: don't forget blanks between symmetry operators; it is advisable to check the Laue symmetry and symmetry operators in the output file especially for those space groups for which alternative origins are shown (i.e. use the setting with -1 at the origin). Upper and/or lower case characters can be used. Some space groups are not correctly generated, for those cases you have to change the setting or give your own symmetry operators (see above ISYM). For cubic space group use the old notation, e.g. F d 3 m instead of F d -3 m. The space group symbol must be given even in the case that you are giving your own symmetry operators. The reflections (if they are not read from file) will be generated according to the space group symbol. A comment can be put after column 20. If this comment starts with the keyword "Mag" (without quotes) then the following line is read if JBT=10 ---------------------------------------------------------------------------Line 11-3-0*: Time_rev(i)(i=1,NS+1) (Only if JBT=+10/-10 and Comment=Mag) (up to 25 integers) --------------Comment line : ! Time Reversal Operations on Crystal Space Group ---------------------------NS is the number of independent symmetry operators given in file CODFIL.OUT for the crystallographic space group. Time_rev(i)=-1 if time reversal is associated to operator "i" for magnetic symmetry, otherwise is equal to 1. The order of operators is the same as in CODFIL.OUT, so a first run is needed for knowing the list of crystallographic symmetry operators. For centrosymmetric groups Time_rev(NS+1) tells the program if time reversal is associated (-1) or not (1) to the inversion operator. This last item should be given only for centrosymmetric space groups. This approach assumes that the magnetic symmetry belongs to the family of the crystallographic space group. However the user can treat the problem using subgroups of the space group (making the appropriate constraints in the atomic positions) when needed. Ex of lines 11-3/11-3-0*: P 6/m m m Magnetic symmetry below ! Time Reversal Operations on Crystal Space Group 1 1 1 1 1 1 -1 -1 -1 -1 -1 -1 1 ---------------------------------------------------------------------------Line 11-3-1*: MULT, ICENT, NLAUE, NMAGR (Only if ISYM><0) (4 integers) --------------Comment line : !Nsym Cen Laue MagMat ---------------------------MULT = Number or symmetry operators given below. ICENT= 1 Non centrosymmetric structure

2 Centrosymmetric structure NLAUE= Integer corresponding to the following laue classes: 1:(-1), 2:(2/m), 3:(mmm), 4:(4/m), 5:(4/mmm), 6:(-3,R), 7:(-3m,R), 8:(-3), 9:(-3m1), 10:(-31m), 11:(6/m), 12:(6/mmm), 13:(m3), 14:(m3m) This number is only used for checking the symmetry operators given by users. For a phase described in a hexagonal basis one should put NLAUE=6,7...12, even if the space group symbol used for generating the reflections is of different symmetry. NMAGR= Number of magnetic rotation matrices for each symmetry operator. -------------------------------------------------------------------------Line 11-3-2*: MULT(1+NMAGR) lines of the form: => If ISYM=1 the symmetry operators are given in numeric form: --------------Comment line : !S11 S12 S13 T1 S21 S22 S23 !M11 M12 M13 M21 M22 M23 M31 M32 M33 ---------------------------T2 Ph S31 S32 S33 T3

. S11 S12 S13 T1 S21 S22 S23 T2 S31 S32 S33 T3 (3(3Int,1real)) . R11 R12 R13 R21 R22 R23 R31 R32 R33 .Phase (9Int,1real) MULT . .NMAGR lines blocks . . . => If ISYM=-1 the symmetry operators are given in alpha-numeric form: e.g. ! SYMM MSYM ! SYMM MSYM ! SYMM MSYM

X,Y,Z U,V,W, 0.0 X+1/2,-y, Z -U, V,-W, 0.0 -x,-y,-Z U, V, W, 0.0

The symbols U,V,W are used for the Fourier components of the magnetic moments along X,Y,Z. The numerical value following the MSYM operator is the magnetic phase in units of 2pi. ---------------------------------------------------------------------------Line 11-4: 4 or 2 lines for each of the N atoms For X-ray or nuclear Neutron scattering: ----------> ---------------------------------------------------------------------Line 11-4-1: LABEL, NTYP, X, Y, Z, B, N, IOPIN,IOPFIN,N_Type (2A4, 5 reals and 3 integers) --------------Comment line : !Atom Typ X Y Z Biso Occ In Fin N_t /Codes ----------------------------

If JBT=4,-4 (structural model supplied by user): --> Line 11-4-1: LABEL, NTYP, P1, P2, P3, P4, P5, P6, P7, P8 (2A4, 8 reals) {parameters defined by user in STRMOD} --------------Comment line : !Atom Typ p1 p2 ! p9 p10 p11 ---------------------------p3 p12 p4 p13 p5 p14 p6 p15 p7 p16 p8

----------------------------------------------------------------------

For magnetic Neutron scattering: ------------------> ---------------------------------------------------------------------Line 11-4-1: LABEL, NTYP, IMAGR, IK, X, Y, Z, B, N, Mx, My, Mz (2A4, 2 integers and 8 reals) --------------Comment line : !Atom Typ Mag Vek X Y Z Biso ! Ix Iy Iz beta11 beta22 beta33 ---------------------------Occ MagPh Rx Ry Rz

If JBT=5,-5 (magnetic model supplied by user): --> Line 11-4-1: LABEL, NTYP, IMAGR, IK, P1, P2, P3, P4, P5, P6, P7, P8 (2A4,2 integers and 8 reals) {parameters defined by user in MAGMOD} --------------Comment line : (default JBT=5) !Atom Typ Mag Vek X Y Z Biso ! Phi & Theta of Cone-axis + unused params ---------------------------Occ Mom beta Phase

---------------------------------------------------------------------The variables defined either for nuclear or magnetic scattering have the following meaning: LABEL= Identification characters for atom or object. NTYP= Link to scattering data for atom : either NAM from 8.1 or chemical symbol and valence to access internal table (use only upper case letters). See note above in line 8. Also a series of special form factors (under test!) are available with refinable parameters. For using this option NTYP should be equal to one of the following words: SPHS (available), SPHE, ELLI, DISK, TORE SASH (not available yet), FUD1,FUD2,FUD3,FUD4 (to be supplied by user in Fdum1.for (see Form_Factor subroutine). See mathematical section for details. Ordinal number of the magnetic rotation matrices applied to the magnetic moment of the atom. To be given only in the case of a magnetic phase Number of the propagation vector to which the atom

IMAGR=

IK=

contributes. If IK=0 the atom is used for all the propagation vectors in the calculation of structure factor. If IK<0 the atom contributes to VK(abs(IK)) and to the vector VK(abs(IK)+NVK/2) X, Y, Z =fractional atomic coordinates B = N = isotropic displacement (temperature) parameter in angstroms**2 occupation number i.e. chemical occupancy x site multiplicity (can be normalized to the multiplicity of the general position of the group). Ordinal number of first and last symmetry operator applied to the atom, apart from the identity which must always be the first one. Useful to describe pseudosymmetries. This option is normally used when the user supply their own list of symmetry operators (ISYM=1).Be careful with multiplicity of reflections!. It is suggested that the users supply also their list of reflections. If IOPIN=IOPFIN=0 all the symmetry operators are applied. Only used for crystallographic structures. N_Type = 0 -> Isotropic atom (no anisotropic temperature factors are given) 2 -> Anisotropic atom. The temperature factors should be given below. 4 -> The form-factor of this atom is calculated using a special subroutine and refinable parameters should be given below (under test!)

IOPIN,IOPFIN=

Mx,My,Mz =Components along the crystallographic axis of the magnetic moments (Bohr magnetons), if JBT=1. In the case JBT=-1 these three parameters correspond to the spherical components of the magnetic moment, in the following order: M, Phi and Theta. M: magnitude of the magnetic moment Phi and Theta are spherical angles of vector M (see note on JBT=-1) If the magnetic phase is incommensurate or described in the crystallographic cell with the help of a propagation vector, these component are actually the real part of the Fourier components of the magnetic moment of the atom. ---------------------------------------------------------------------------Line 11-4-2: CX, CY, CZ, CB, CN, CMx , CMy , CMz (8 reals) CX, CY, CZ =codewords for fractional atomic coordinates (see below) CB = codeword for isotropic displacement (temperature) parameter CN = codeword for occupation number CMx,CMy,CMz =codewords for magnetic moment components

If JBT=4,-4/5,-5 CP1, CP2, CP3, CP4, CP5, CP6, CP7, CP8 For X-ray or nuclear Neutron scattering: ----------> ---------------------------------------------------------------------------STANDARD: Line 11-4-3: b11, b22, b33, b12, b13, b23 (For N_typ=2) (7 reals) --------------Comment line : ! beta11 beta22 beta33 beta12 beta13 beta23 /Codes ---------------------------bij OR = anisotropic displacement (temperature) parameters (betas) (AND if N_typ=5, but not available)

(For N_typ=4)

--------------Comment line : ! Form-factor refinable parameters ---------------------------Line 11-4-3: Line 11-4-4: Line 11-4-5: f2 f3 (7 reals) Cf1 Cf2 Cf3 f8 f1 f4 Cf4 f5 Cf5 f6 f7

Cf6 Cf7

f9 f10 f11 f12 f13 f14 (7 reals) Line 11-4-6: Cf8 Cf9 Cf10 Cf11 Cf12 Cf13 Cf14 (7 reals) The parameters f1 to f14 are used for describing the form-factor of the current object RIGID BODY OR NON-STANDARD: If JBT=4,-4 ---> P9, P10, P11, P12, P13, P14, P15 {parameters defined by user in STRMOD} ---------------------------------------------------------------------For magnetic Neutron scattering: ------------------> ---------------------------------------------------------------------------Line 11-4-3: Mxi, Myi, Mzi, b11, b22, b33, MPhas (7 reals) If JBT=5,-5 ---> P9, P10, P11, P13, P14, P15, P12 {parameters defined by user in MAGMOD} ---------------------------------------------------------------------Mxi,..= Imaginary components of the Fourier coefficient in Bohr magnetons (If JBT<0, spherical components as for real components Mx,My,Mz, see line 11-4-1) bii = diagonal part of anisotropic temperature factors MPhas = Magnetic phase of the atom (see mathematical section) If JHELIX=1 (see line 11-2-1) the third component Mz is calculated by the program in order to have an imaginary vector orthogonal to the real vector. If JBT<0, then the phi-angle of the imaginary part is calculated by the program for keeping the orthogonal constraint.

---------------------------------------------------------------------------Line 11-4-4: CB11, CB22, CB33, CB12, CB13, CB23, CMPhas (7 reals) CBIJ = codeword for anisotropic displacement (temperature) parameters, or for imaginary components of the magnetic Fourier coefficient. CMPhas= codeword for magnetic phase. If JBT=4,-4 CP9, CP10, CP11, CP12, CP13, CP14, CP15 ************************************************************************ New input format for atom parameters The lines 11-4-1 to 11-4-4 should be changed for the case JBT=10,-10 ************************************************************************ For Xrays or nuclear+magnetic Neutron scattering: ----------> ---------------------------------------------------------------------Line 11-4-1: LABEL, NTYP, IMAGR, IK, X, Y, Z, (2a4,2int,5real,1int) Line 11-4-2: CX, CY, CZ, (5 reals) Line 11-4-3: Mx My Mz Mxi Myi (7 reals) Line 11-4-4: CMx CMy CMz CMxi CMyi (7 reals) Line 11-4-5: b11 b22 b33 b12 b13 (6 reals) Line 11-4-6: Cb11 Cb22 Cb33 Cb12 Cb13 (6 reals) Line 11-4-7: f1 f2 f3 f4 f5 (7 reals) Line 11-4-8: Cf1 Cf2 Cf3 Cf4 Cf5 (7 reals) Line 11-4-9: f8 f9 f10 f11 f12 (7 reals) Line 11-4-10: Cf8 Cf9 Cf10 Cf11 Cf12 (7 reals) If If If If if N_type N_type N_type N_type N_type = = = = = 0 1 2 3 4 Only give give give give lines (a) the lines the lines the lines the lines B, N N_type (a) CB, CN Mzi (b) MPhas (c)

CMzi CMPhas (d) b23 Cb23 f6 f7 (e) (f) (g) (h) (i) (j)

Cf6 Cf7 f13 f14 Cf13 Cf14

and (b) need to be given (a), (b), (c) and (d) (a), (b), (e) and (f) (a) -> (f) (a), (b) and (g)->(j) (special form-factor)

--------------Comment lines : !Atom Typ Mag Vek X Y !Line below:Codes ! Rx Ry Rz Ix !Line below:Codes or... ! beta11 beta22 beta33 beta12 ----------------------------

Z Iy beta13

Biso Iz

Occ MagPh

N_type

beta23 /Line below:Codes

This input could be also used for Xrays, in such case IMAGR and IK should be zero for all the atoms and Jobtyp or Jtyp(n)=0. In such case the space group symbol can be used for generation of reflections and symmetry operators.

For a phase with magnetic contributions NTYP should be equal to the magnetic form factor symbol. The program extracts internally the fermi length symbol from NTYP. If there are magnetic contributions the symmetry should be controlled by the user (ISYM=1) and the magnetic part should be described with the formalism of propagation vectors, the magnetic contribution is calculated only for the satellite reflections. If fundamental reflections have magnetic contribution the propagation vector k=(0,0,0) must be included explicitely if there are other propagation vectors. If the magnetic cell is the same as the chemical cell propagation vectors are not needed. For the moment, the symmetry operators must belong to the group of the propagation vector Gk, so some atoms need,in general, to be repeated for the rest of positions not generated by Gk. ------------------------------------------------------------------------------------------------------------------------------------------------====> For iabs(IRF)<>4 (no integrated intensity data) and ConstWavelength ---------------------------------------------------------------------------Line 11-5-1: S, GAM1, Bov, STR1, STR2, STR3, IstrainModel (6 reals and 1 integer) --------------Comment line : ! Scale Shape1 Bov Str1 Str2 Str3 Strain-Model ---------------------------S GAM1 = scale factor = profile shape parameter , e.g. : eta0 for NPROF=4,5 but not for NPROF=7 (see line 11-6-1) m0 for NPROF=6 = overall isotropic displacement (temperature) factor in angstroms**2

Bov

STR1,STR2,STR3 = strain parameters, defined through the subroutine STRAIN (see additional information) If ISTR=1 set these values to 0.0 = Integer to select a particular model for strains in subroutine STRAIN. ----------------------------------------------------------------------------------------------------------------------------------------------------Line 11-5-2: CS, FLGAM1, CBov, CSTR1,CSTR2,CSTR3 (6 reals) CS = codeword for scale factor FLGAM1 = codeword for GAM1 CBov = codeword for overall isotropic displacement (temperature) factor CSTR1,CSTR2,CSTR3= codeword for strain parameters If ISTR=1 set these values to 0.0 -----------------------------------------------------------------------------------------------------------------------------------------------------====> For iabs(IRF)=4 (integrated intensity data) ---------------------------------------------------------------------------Line 11-5-1: Sc1, Sc2, Sc3, Sc4, Sc5, Sc6 (6 reals) IstrainModel

--------------Comment line : ! Scale Factors ! Sc1 Sc2 Sc3 ---------------------------Line 11-5-2:

Sc4

Sc5

Sc6

CSc1, CSc2, CSc3, CSc4, CSc5, CSc6 (6 reals) = = scale factor for domain (i) code of the scale factor for domain (i)

Sci CSci

For powder data only the first scale factor is used -------------------------------------------------------------------------===> End iabs(IRF)=4 ---------------------------------------------------------------------------====> For iabs(IRF)<>4 (no integrated intensity data) and ConstWavelength ---------------------------------------------------------------------------Line 11-6-1 U, V, W, X, Y, IG, SZ, IsizeModel FWHM (or shape) parameters : --------------Comment line : ! U V W X Y GauSiz or for NPRO=11 (split pseudo-Voigt) ! Ul Vl Wl Xl Y GauSiz ---------------------------(7 reals and 1 integer) For profiles 0 to 6 and 12: FWHM^2 = (U+DST^2)*Tan^2(Theta) + V*Tan(Theta) + W + IG/cos^2(Theta) =====> For NPROF = 4 (tripled pseudo-Voigt), the three components are assumed to have the same eta0 and FWHM, so the effective total width depends on the additional shape parameter Shp1 (see line 11-8-3). The profile function is given by the formula: p4(x) = X*pV(x-D) + (1-X-Y)*pV(x) + Y*pV(x+D) where D = Shp1/d.costheta pV(x) = Eta0*L(x)+(1-Eta0)*G(x) So, apart from the FWHM that is calculated from U,V,W DST and IG parameters for a single component, the profile function has FOUR shape parameters Eta0, X, Y and Shp1. This function is adapted for medium resolution neutron powder diffractomers having defects on the monochromator and/or the guide spacial spectral distribution giving rise to a non-gaussian distribution of wavelengths. =====> For NPROF = 5 and 12 (pseudo-Voigt) the eta parameter can be dependent on X through the formula: pV(x) = Eta*L(x)+(1-Eta)*G(x) Eta = Eta0 + X * 2Theta LorSiz Size-Model LorSiz Size-Model

=====>

For NPROF = 11 (split pseudo-Voigt) DST and IG are common to the left and right parts of the profile. Moreover additional FWHM parameters are used as new shape parameters, so the expression of the left FWHM(L) for NPROF=11 is

FWHM^2(L) = (Ul+DST^2)*Tan^2(Theta) + Vl*Tan(Theta) + Wl+IG/cos^2(Theta)+ Shp1/tan^2(2Theta) Shp1 is applied only for 2Theta<=90. The value of ETA for the left part is given by: Etal = Etal0 + Xl * 2Theta The expression of the right part is: FWHM^2(R) = (Ur+DST^2)*Tan^2(Theta) + Vr*Tan(Theta) + Wr+IG/cos^2(Theta) + Shp2/tan^2(2Theta) Shp2 is applied only for 2Theta >90. The value of ETA for the right part is given by: Etar = Etar0 + Xr * 2Theta The FWHM and shape parameters for the right part are read in next lines =====> For NPROF = 6 (Pearson-VII) the m parameter can be dependent on X and Y through the formula: m = m0 + 100* X / 2Theta + 10000* Y / (2Theta)**2 =====> ta) FWHM(lorentzian) = X tan(Theta) + (Y+ F(SZ))/cos(Theta) All expressions are in (degrees 2Theta)^2 U,V,W = Half-width parameters (normally characterizing the instrumental resolution function). X = Lorentzian isotropic strain parameter. DST(STR) = Anisotropic gaussian contribution of microstrain. It is calculated in subroutine STRAIN as a function of IstrainModel or ISTR. If IstrainModel <>0 then ISTR must be zero. DST depends on STR1,STR2,...parameters and hkl. IG= Isotropic size parameter of gaussian character F(SZ)= anisotropic lorentzian contribution of particle size. It is calculated in subroutine SIZE and depend on parameter SZ and hkl. IsizeModel= Integer to select a particular model for F(SZ) in subroutine SIZE. ---------------------------------------------------------------------------Line 11-6-2: CU, CV, CW, CX, CY CIG CSZ : codewords for the FWHM (or shape) parameters (7 reals) ---------------------------------------------------------------------------Line 11-6-3: Ur, Vr, Wr, Etar0, Xr (5 reals) Read only if NPRO=11 (split pseudo-Voigt) For NPROF = 7, the FWHM of the two components is calculated as

FWHM^2(gaussian) = (U+DST^2)*Tan^2(Theta) + V*Tan(Theta) + W + IG/cos^2(The

--------------Comment line : ! Ur Vr Wr Etar0 Xr ---------------------------FWHM and shape parameters for the right part of the split pseudo-Voigt function. This function is similar to NPROF=5 but the left (x<0) and right (x>0) parts of the profile have different U,V,W,eta0 and X parameters. Additional shape parameters are also read. Line 11-6-4*: CUr, CVr, CWr, CEtar0, CXr: codewords for the FWHM and shape parameters (5 reals) Read only if NPRO=11 ---------------------------------------------------------------------------===> End iabs(IRF)<>4 ---------------------------------------------------------------------------====> For iabs(IRF)=4 (integrated intensity data) ---------------------------------------------------------------------------Line 11-6-1: Ext1, Ext2, Ext3, Ext4, Ext5, Ext6, Ext7 (7 reals) --------------Comment line : ! Extinction Parameters ! Ext1 Ext2 Ext3 Ext4 Ext5 Ext6 ---------------------------CExt1, CExt2, CExt3, CExt4, CExt5, CExt6, CExt7 (7 reals) = = Extinction parameter (i) Code of the extinction parameter (i)

Ext7

Line 11-6-2:

Exti CExti

At present only the first extinction parameter is used. ===> End iabs(IRF)=4 ---------------------------------------------------------------------------Line 11-7-1: a, b, c, a, b, g cell parameters in A and degrees (6 reals) --------------Comment line : ! a b ---------------------------c alpha beta gamma

---------------------------------------------------------------------------Line 11-7-2: CA, CB, CC, CD, CE, CF : (6 reals) codewords for cell constants A, B, C, D, E & F defined by : 1/d2 = Ah2 + Bk2 + Cl2 + Dkl + Ehl + Fhk Note that these codewords do not refer directly to the cell parameters; for instance, in the hexagonal system, the last codeword CF must be the same as CA and CB. ---------------------------------------------------------------------------Line 11-8-1: G1, G2, Pas1, Pas2, Pas3, Pas4 (6 reals) --------------Comment line :

! Pref1 Pref2 Asy1 Asy2 ----------------------------

Asy3

Asy4

G1, G2 =preferred orientation parameters (see Math. section) when NORI = 0, G1 = 0 means no preferred orientation when NORI = 1, G1 = 1 means no preferred orientation Pa1,..Pas4 =asymmetry parameters applied to angles below RLIM (given on line 4 : see Mathematical section) If NPHASE is negative only the first parameter is relevant. ---------------------------------------------------------------------------Line 11-8-2: CG1, CG2, CPas1, CPas2, CPas3, CPas4 (6 reals) CG1, CG2 = codewords for preferred orientation parameters CPas1,...CPas4 = codewords for asymmetry parameters ---------------------------------------------------------------------------Line 11-8-3*: Shp1, CShp1, Shp2, CShp2 (4 reals) --------------Comment line : !Additional Shape parameters ---------------------------Additional shape parameters and corresponding codewords. Read only if NPRO=4 or if NPRO >8 For NPRO=11 (split pseudo-Voigt) they correspond to the additional contribution to the FWHM for the Left (L) and Right (R) part of the profile for 2Theta<90 and 2Theta>90 respectively. addFWHM^2(L) = Shp1/tan^2(2Theta) , addFWHM^2(R) = Shp2/tan^2(2Thet a) For NPRO=12 (Convoluted pseudo-Voigt with axial divergence asymmetry) Shp1= S_L is source width/detector distance Shp2= D_L is detector width/detector distance These parameters play the role of asymmetry parameters, they are used only for reflections below 2Theta = RLIM. ---------------------------------------------------------------------------Line 11-8-4*: U2, V2, W2 (3 reals) U,V,W parameters for the second wavelength present in the diffraction pattern. Read only if RATIO is negative. --------------Comment line : !Additional U,V,W parameters for Lambda2 ------------------------------------------------------------------------------------------------------Line 11-8-5*: CU2, CV2, CW2 (3 reals) Codewords of the additional U,V,W parameters Read only if RATIO is negative. ---------------------------------------------------------------------------TIME OF FLIGHT DATA =============================================================================== ====> LINES 11-5-1 to 11-8-5 are substituted by the following lines for TIME OF FLIGHT DATA

=============================================================================== Line 11-5-1: S, Ext, Bov, STR1, STR2, STR3, IstrainModel (6 reals + 1 integer) --------------Comment line : ! Scale Extinc Bov Str1 Str2 Str3 Strain-Model --------------------------------------------------------------------------Same parameters as in 11-5-1 for CW, except that GAM1 is replaced by the extinction parameter EXT. Line 11-5-2: CS, FLEXT, CBov, CSTR1, CSTR2, CSTR3 (6 reals) Codewords of the above parameters

---------------------------------------------------------------------------Line 11-6-1: Sig2, Sig1, Sig0, Xt, Yt, Z1, Z0, IsizeModel (7 reals and 1 integer) --------------Comment line : ! Sig-2 Sig-1 Sig-0 l ------------------------------------------------------------------------------Gaussian FWHM parameters : (d2=d*d,d4=d2*d2, d: d-spacing) Sigma^2 = (sig2 + GSIZ ) d4 + (sig1 + DST ) * d2 + sig0 Xt,Yt : Not used at present Z1 = GSIZ : Gaussian isotropic size component Z0 : Not used at present Units :: sig2,GSIZ: (microsecs/Angstrom^2)^2 sig1, DST: (microsecs/Angstrom)^2 sig0 : (microsecs)^2 Xt Yt Z1 Z0 Size-Mode

DST depends on STR1, STR2, ... through the selected strain model. The Gaussian FWHM is Sigma*sqrt(8 Ln2) Line 11-6-2 CSig2, CSig1, CSig0, CXt, CYt, CZ1, CZ0 (7 reals) Codewords of the above parameters

---------------------------------------------------------------------------Line 11-6-3: Gam2, Gam1, Gam0, LStr, LSiz (5 reals) --------------Comment line : ! Gam-2 Gam-1 Gam-0 LStr LSiz ---------------------------------------------------Lorentzian FWHM parameters : gamma = (gam2 + DSIZ ) d2 + (gam1 + LStr) * d + gam0 gamma: Lorentzian FWHM LStr : Lorentzian isotropic strain LSiz : Lorentzian isotropic strain

DSIZ=F(LSiz) : F depends on LSiz (and eventually on more size parameters) through the selected size model Units :: gam2,DSIZ: microsecs/Angstrom^2 gam1,LStr: microsecs/Angstrom gam0 : microsecs

Line 11-6-4

CGam2, CGam1, CGam0, CLStr, CLSiz (5 reals) Codewords of the above parameters

---------------------------------------------------------------------------Line 11-7-1: a, b, c, alpha, beta, gamma (6 reals) --------------Comment line : ! a b c alpha beta gamma ---------------------------------------------------------------------------------------------------------------------------------------------Line 11-7-2: CA, CB, CC, CD, CE, CF : (6 reals) Cell parameters and codewords as above. ---------------------------------------------------------------------------Line 11-8-1: G1, G2, alph0, beta0, alph1, beta1 (6 reals) --------------Comment line : ! Pref1 Pref2 alph0 beta0 alph1 beta1 -------------------------------------------------G1, G2 : Preferred orientation parameters (as above) alph0, beta0, alph1, beta1 : Parameters defining the variation of the exponential decay constants with d-spacing. Fast decay: alpha = alpha0 + alpha1/d Slow decay: beta = beta0 + beta1/d4 alpha and beta are in reciprocal microseconds and d in angstroms. Line 11-8-2: CG1, CG2, Calph0, Cbeta0, Calph1, Cbeta1 (6 reals) Codewords of the above parameters ---------------------------------------------------------------------------Abs1, CAbs1, Abs2, CAbs2 (4 reals) --------------Comment line : !Absorption correction parameters ---------------------------------Abs1, Abs2 CAbs1, CAbs2

Line 11-8-3:

: Absorption correction parameters : Codewords The physical meaning of these parameters depend on the function selected by IABSCOR (see Line 4) ----------------------------------------------------------------------------

=====================END=OF=SPECIFIC=INPUT=FOR=T.O.F.=========================== Line 11-8-6*: Ahkl, Shf1, Shf2, IASV,ISHIF (3 reals and 2 integers) --------------Comment line : ! AsyP Shift1 Shift2 ModA ModS ---------------------------Read only if JSOL=1 Ahkl = HKL-dependent asymmetry parameter. Shf1,Shf2 = HKL-dependent shift parameters. The three last parameters are defined by the user through the subroutines ASYMHKL & SHIFHKL, where a particular model for displacement and asymmetry of Bragg reflections is built. IASV= Model for asymmetry. ISHIF= Model for shifts.

---------------------------------------------------------------------------Line 11-8-7*: CAhkl, CShf1, CShf2 Read only if JSOL=1 (3 reals) CAhkl = codeword for hkl-dependent asymmetry parameter CShf1, CShf2 = codewords for hkl-dependent shift parameters. ---------------------------------------------------------------------------Line 11-8-8* : Sh1,Sh2,Sh3 --------------Comment line : 'Shift-cos(1) or Shift-sin(-1) axis' in the same line as the numbers ---------------------------(3 reals) If Ishift= +/- 1, [Sh1,Sh2,Sh3] is the vector defining the axial "shift-platelets". ---------------------------------------------------------------------------Line 11-9* : Sz1,Sz2,Sz3 (3 reals) If IsizeModel= +/- 1, [Sz1,Sz2,Sz3] is the vector defining the platelets. ---------------------------------------------------------------------------Line 11-10-1* : St1,St2,St3 (3 reals) If IstrainModel=7, [St1,St2,St3] is the vector defining the axial microstrain. ---------------------------------------------------------------------------Line 11-10-2* : Str4, Str5, Str6, Str7, Str8 (5 reals/5 reals) : CStr4,CStr5,CStr6,CStr7,CStr8 If IstrainModel > 8, 5 additional strain parameters and codes. --------------Comment line : ! 5 additional strain parameters (IstrainModel>8) ------------------------------------------------------------------------------------------------------Line 11-11*: NVK pair of lines with propagation vectors and codes Components of K in reciprocal lattice units. (3 reals) Kx, Ky, Kz (3 reals) CKx, CKy, CKz --------------Comment line :

! Propagation vectors: ------------------------------------------------------------------------------------------------------Line 11-12*: IFURT lines with further parameters introduced by users, each line contains: (A4,2 reals): NAMEPAR --------------Comment line : ! Further parameters: ---------------------------VALUEPAR CODEPAR

---------------------------------------------------------------------------Line 11-13*: Generalized strain parameters and codewords (5 reals/5 reals) If ISTR=1,3 (STR(j),j=1,5),(CODSTR(j),j=1,5) (STR(j),j=6,10),(CODSTR(j),j=6,10) --------------Comment line : ! Additional strain parameters: ---------------------------Two sets of five strain parameters and their codes Eg. STR1 STR2 ... STR5 CSTR1 CSTR2 ... CSTR5 STR6 STR7 ... STR10 CSTR6 CSTR7 ...CSTR10 ---------------------------------------------------------------------------Line 11-14*: Generalized size parameters and codewords (6 reals/6 reals) If ISTR=2,3 --------------Comment line : ! Generalized strain parameters: ---------------------------(Siz(j),j=1,6),(CODSiz(j),j=1,6) A set of six size parameters and their codes Eg. Siz1 Siz2 ... Siz6 CSiz1 CSiz2 ... CSiz6 ---------------------------------------------------------------------------Line 11-15*: NDIST number of lines of distance constraints CATOD1, CATOD2, ITnum, T1, T2, T3, Dist, Sigma (2A4, 1 integer and 5 reals) --------------Comment line : ! Soft distance constraints: ---------------------------CATOD1 and CATOD2: Names of the atoms to be constrained They must coincide with labels in the asymmetric unit. ITnum: Integer for selecting the rotation part of the symmetry operator to be applied to the coordinates of the atom CATOD2.

(T1, T2, T3): Translation part of the above symmetry operator Dist : Value of the required distance. Sigma : Standard deviation of the distance. The numbering of symmetry operators to be given in distance constraints conditions. The integer number to be given is ITnum. If combination with a center of symmetry is needed the value must be entered as negative. Non-hexagonal frames ITnum ( 1) ( 2) ( 3) ( 4) ( 5) ( 6) ( 7) ( 8) ( 9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) Symmetry symbol Rotation matrix 1 --> ( x, y, z) 2 ( 0, 0, z) --> (-x,-y, z) 2 ( 0, y, 0) --> (-x, y,-z) 2 ( x, 0, 0) --> ( x,-y,-z) 3+ ( x, x, x) --> ( z, x, y) 3+ (-x, x,-x) --> ( z,-x,-y) 3+ ( x,-x,-x) --> (-z,-x, y) 3+ (-x,-x, x) --> (-z, x,-y) 3- ( x, x, x) --> ( y, z, x) 3- ( x,-x,-x) --> (-y, z,-x) 3- (-x,-x, x) --> ( y,-z,-x) 3- (-x, x,-x) --> (-y,-z, x) 2 ( x, x, 0) --> ( y, x,-z) 2 ( x,-x, 0) --> (-y,-x,-z) 4- ( 0, 0, z) --> ( y,-x, z) 4+ ( 0, 0, z) --> (-y, x, z) 4- ( x, 0, 0) --> ( x, z,-y) 2 ( 0, y, y) --> (-x, z, y) 2 ( 0, y,-y) --> (-x,-z,-y) 4+ ( x, 0, 0) --> ( x,-z, y) 4+ ( 0, y, 0) --> ( z, y,-x) 2 ( x, 0, x) --> ( z,-y, x) 4- ( 0, y, 0) --> (-z, y, x) 2 (-x, 0, x) --> (-z,-y,-x) Hexagonal frames ITnum (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36) Symmetry symbol 1 --> 3+ ( 0, 0, z) --> 3- ( 0, 0, z) --> 2 ( 0, 0, z) --> 6- ( 0, 0, z) --> 6+ ( 0, 0, z) --> 2 ( x, x, 0) --> 2 ( x, 0, 0) --> 2 ( 0, y, 0) --> 2 ( x,-x, 0) --> 2 ( x,2x, 0) --> 2 (2x, x, 0) --> Rotation matrix ( x , y, z) ( -y, x-y, z) (-x+y,-x , z) (-x , -y, z) ( y,-x+y, z) ( x-y, x , z) ( y, x ,-z) ( x-y, -y,-z) (-x ,-x+y,-z) ( -y,-x ,-z) (-x+y, y,-z) ( x , x-y,-z)

---------------------------------------------------------------------------Line 11-16*: NMAGC number of lines of magnetic moment constraints --------------Comment line : ! Soft moment constraints:

---------------------------(A2,2 reals) : CATOM, Moment, Sigma CATOM: Two letters equal to the two first character of the label of atoms in asymmetric unit which are constrained. Moment: Value of the required magnetic moment. Sigma : Standard deviation of Moment. (It doesn't work with incommensurate magnetic structures) ---------------------------------------------------------------------------LINE 12* : NRELL lines containing the following items (1 integer and 2 reals) --------------Comment line : ! Hard limits for selected parameters: ---------------------------NUMPAR, LowLIMIT, HighLIMIT Where NUMPAR is the "number" of the parameter (as given by the parameter code) to be constrained within the limits specified by the interval [LowLIMIT, HighLIMIT] For a proper use of this option one has to put limits to the variable appearing for the first time with the wished parameter code. This should have a positive sign and a unit multiplier. ---------------------------------------------------------------------------Line 12-1*: NCONFG, NSOLU, NREFLEX,NSCALEF (3 integers) --------------Comment line : ! Nconfg Nsolu Num_Ref Nscalef -----------------------------------Read only if ICRYG=2. The program tries NCONFG configurations and select the best NSOLU solutions (lower R-factors) using the first NREFLEX reflections of the file CODFILn.HKL. If NSCALEF is different from zero, then the scale factor used in the program is obtained from the relation: Sum{Iobs}= Scale*Sum{Icalc} A configuration means a set of NRELL values of the selected parameters within the box defined in LINE 12. The constraints established with the coding of parameters have no the same meaning as with least-squares(LS). In LS refinement for two variables having codes xx1.00 and xx0.5 the shift applied to the second variable is half the shift applied to the first one irrespective of their initial values. In Montecarlo search, the value of the second variable is just half the value of the first variable. Only one parameter, numbered here as "xx", controls the value of the two variables either in LS or in Montecarlo searh.

---------------------------------------------------------------------------LINE 13* : ISCALE, IDIF used only if IPL.ne.0 (line 3) (2 integers) --------------Comment line : ! Iscale Idif ----------------------------

ISCALE =counts per character position for observed and calculated curves on line print plot IDIF = counts per character for difference curve ---------------------------------------------------------------------------LINE 14* : THET1, THET2 (2 reals) The reflection list between these angles is saved in the file CODFILHKL.SAV --------------Comment line : ! 2Th1 2Th2 ---------------------------==================END OF CODFIL.PCR'S DESCRIPTION=============================

3.- MATHEMATICAL INFORMATION

--------------------------------------------3.1: Calculated profile. Structure factors. --------------------------------------------Calculated counts yci at the ith step are determined by summing the contribution from neighbouring Bragg reflections plus the background : yci = S SUM(Lh.Fh^2 .Omeg(Ti - Th).Ah.Th. Ph + ybi where S is the scale factor Lh contains the Lorentz, polarization and multiplicity factors Fh is the structure factor. The ratio of the intensities for the two wavelengths is absorbed in the calculation of Fh^2, so that only a single scale factor is required. Ah is the asymmetry function Th is the transmission factor (line 4) Ph describes the preferred orientation of the sample Omeg is the reflection profile function which approximates the effects of both instrumental and, possibly, specimen parameters. ybi is the background intensity -----------------------------------------Form-factor calculations and Refinements -----------------------------------------Apart from the standard scattering factor for indivual atoms existing in an internal library of FullProf. The version 3.2 and higher can handle comple

x form-factors as a standard option. In the general expression of the nuclear structure factor: F(H)=Sum(s) {ns.f(H)s. Sum(j)[Tjs.exp(2pi( H Gj rs +tj)]} the form factor f(H) is normally dependent on the module of H. For molecular plastic crystals the treatment of rotating molecules cannot be done using an atomic description. The approach of a molecular form-factor that takes into account the particular dynamics of the object is more reliable. f(H) depends on a series of parameters for different types of objects. Coefficients of Symmetry Adapted Spherical Harmonics, geometrical parameters (radius of a sphere, length and radius of a cylinder or disk, etc), scatterin g density, etc could serve for describing the scattering factor of a complex object. In the FullProf version 3.2 or higher the available (or projected) objects are the following: Sphere: Elipsoidal: Cylinder of elliptical section: ---------------------------------------Magnetic scattering calculations ---------------------------------------For a magnetic phase F(Q)^2 is calculated using the general expresion of Halpern and Johnson: F(Q)^2 = {Fm(Q)^2 - (e.Fm(Q)^2}

where Fm(Q) is the magnetic structure factor, Q=H+k and e is the unit vector along the scattering vector Q. The magnetic moment is considered as a Fourier superposition of type: m(l,j) = Sum(k) { S(k,j) exp[-2.pi.i.k.R(l)]} In such a case the magnetic structure factor is given by: Fm(H+k) = Sum(j){S(k,j) fj(H+k) exp [2.pi.i.(H+k).r(j)]}

The Sum(j) is over all the atoms in the crystallographic cell If symmetry relations are established for coupling the different Fourier components S(k,j), phase factors are added to the exponential and the sum is only for the asymmetric unit: Fm(Q=H+k)=0.2695 Sum(s){ns.f(Q)s. Sum(j)[(Rj(s).S(k,s)Tjs.exp(2pi(Q Mj rs - Psik(j,s))]} The sum over (s) concerns the magnetic atoms of the asymmetric unit for the wavevector k (the Fourier component k contributes only to k-satellite),

ns is the occupation factor and f(Q)s is the form factor of atom s. The sum over (j) concerns the different symmetry operators of the crystal space group Mj={g,t}j (g and t are the rotational and translational part of Mj. The matrix Rj(s) transform the components of the Fourier term S(k,s) of the starting atom s to that numbered as "j" in the orbit of s. Tjs is the temperature factor. The phase factor Psik(j,s) has two components: Psik(j,s) = Mphas(s) + Phase(j) Mphas(s) is a phase factor which is not determined by symmetry. It is a refinable parameter and it is significant only for an independent set of magnetic atoms which respect to another one. Phase(j) is a phase factor determined by symmetry. The Fourier component k of the magnetic moment of atom s, S(k,s) is transformed to S(k,sj) = Rj(s) S(k,s) exp {-2pi i Phase(j)} The sign of Psik(j,s) changes for -k. The reflection H+k has the "negative" sign indicated in the above formulas and the reflection H-k has the positive sign. In the general case S(k,s) is a complex vector (in general there are six components. In old versions of FullProf one could simulate this complex vector by splitting the atom contributing to propagation vector k in two parts by putting a "magnetic phase" of pi/2 with respect to the real part of S(k,s). The magnetic phases are given in fractions of 2pi, then for the above purpose one can use a fixed Phase=0.25. In the new version the above splitting is not needed. The imaginary components are read in the place which was reserved for anisotropic temperature factors (see lines 11-4-3). For the scattering vector H-k the Fourier component is the complex conjugate of the Fourier component used for calculating the structure factor for H+k. The program takes into account automatically this fact. If k is at the interior of the Brillouin Zone a factor 1/2 is applied to the Fourier coefficient. Let us consider a single index "j" for the sublattice j of the site "s". The Fourier coefficient for the sublattice is given by: S(k,j)= { 1/2 [MRxj e1 + MRyj e2 + MRzj e3] + 1/2 i [MIxj e1 + MIyj e2 + MIzj e3] } exp(-2.pi.i.Psik(j)) The vector -k must also be given either explicitly or implicitly by giving NVK<0 (see NVK on line 11-2 ). If NVK<0 the program applies the factor 1/2 because it is supposed that k is non equivalent to -k even if k belong to the surface of the Brillouin zone. If the option JHELIX=1 is used, the number of free parameters per magnetic atom is reduced. The Fourier coefficients are considered of the form: S(k,j)= 1/2 [m1s uj + i m2j vj] exp(-2.pi.i.Psik(j)) where uj and vj are orthogonal unit vectors. If m1j=m2j=m0 the magnetic structure for the sublattice j corresponds to a classical helix (or spiral) of cylindrical envelope. All j atoms have a magnetic moment equal to m0. If m1j/=m2j the helix has an elliptical envelope and the moments have values between min(m1j,m2j) and max(m1j,m2j). If m2j=0 the magnetic structure corresponds to a modulated sinusoid of amplitude A=m1j. In general, the user has to calculate the real magnetic moments from the

refined values of the Fourier components: the phrase "Magnetic Moment" in the output file means the modulus of the corresponding Fourier component. The program MOMENT has been written in order to help the user with these calculations. In any case the calculation of the magnetic moment of the atom "j" in the unit cell of index "l" should be done by using the formula: m(l,j) = Sum(k) { S(k,j) exp[-2.pi.i.k.R(l)]}= = Sum[k] {[MRxj e1+ MRyj e2+ MRzj e3] cos2pi[k.R(l)+Psik(j)]+ [MIxj e1+ MIyj e2+ MIzj e3] sin2pi[k.R(l)+Psik(j)]} where Sum[k] is the sum extended for half the number of propagation vectors, i.e. over the number of pairs (k,-k). If the propagation vector k is commensurate (rational components) one can use the magnetic unit cell and m(k,j) can be identified with the magnetic moment at site j. In this case one can describe the magnetic structure with Psik(j)=0 and Q=H, being H an integer vector of the reciprocal lattice of the magnetic cell. If k=1/2H, one can use the chemical unit cell and real magnetic moments. In such a case only one propagation vector is needed: if NVK is given as negative the generation of magnetic reflections could be in error. For centered crystallographic unit cells one can use only the content of a primitive cell and generate the satellites from the symbol of the centering followed by -1 (e.g. I -1 for a I-centered cell). In order to take the advantage of the crystallographic conventions (propagation vector given with respect to the reciprocal basis of the conventional cell) one can use the dimensions and the metrics of the convencional cell provided that, putting the content of a primitive cell in the conventional cell frame, the ocupation factors are multiplied by the number of centering vectors. See the two files HOBK1.PCR and HOBK2.PCR in the anonymous FTP area. -----------------3.2: Background -----------------Background intensity ybi at the ith step is obtained (line 2 and 6* or 10-2*) either from an user-supplied table of background intensities (optional lines 6*),or from a refinable background function ybi = SUM(Bm.(Ti/BKPOS - 1)^m) + SUM(Bcj sin[Qidj]/Qidj) with 0<=m<=5 for NBCKGD=0, or 0<=m<=11 for NBCKGD=-3. The origin of the background polynomial is given by a selectable input parameter BKPOS (line 4) and should be supplied by the user. The second sum (six terms) is used only if NBCKGR=-1. The parameter to be refined are: Bo,B1,...B5, Bc1,...Bc6, d1,...d6 Qi is given by Qi=4pi.sin(Thetai)/Lambda(1) The parameters dj are distances in Angstroms. The background can be also read from a supplied file FILE.BAC. The actual background is calculated from the read background applying the following formula: Backg(2theta) = a * BackgRead [ (1+c) 2theta + d ] + b

The background parameters and codes, given in line 10-2*, correspond to the coefficients of the above formula in the following order: a,b,c,d. If a is given as zero, the program puts a=1. Limits against divergence are fixed by program. The parameter c is allowed to vary up to a maximum value abs(c)=0.1 and abs(d) is kept below 3 degrees (2theta). The user can chek the excursion of those parameters out of the allowed range when they are strictly zero and their standard deviation is (fixed arbitrarily to) 0.99999. This option is useful when complicated background shapes are present due to sample environment. A new procedure for treating the background has been introduced. The background is ajusted iteratively at each cycle by using a Fourier filtering technique. The starting background is read from a file. At cycle "n" the new background is calculated from the old one, cycle "n-1", with the formula: Back(n) = Back(n-1) + Filtered[yobs-ycal](n-1) where Filtered[yobs-ycal] is a strongly smoothed version of yobs-ycal. The parameter controlling the smoothing is FWINDOW which is equivalent to PST described in subroutine SMOOFT described in Numerical Recipes (see ref 19). The implementation of SMOOFT in FullProf is not the same as in ref 19. When using this method it is wise to draw the final background to see if it is really a smooth curve. This option is only justifiable in cases of a very wavy background. The starting background should be close to the real one. ---------------------------3.3: Peak-shape functions ---------------------------The profile function is selected by the control variable NPROF (line 2). The currently available functions are given above and their particular definitions can be found in the literature. --------------------------------------------------------3.4: Monochromator, Lorentz and geometrical corrections --------------------------------------------------------Monochromator polarization correction (CTHM and K in line 4); the Lorentz polarization factor LP is calculated as : LP = [1 - K + K.CTHM.cos(2Theta)^2]/2(sinTheta)^2.cosTheta with CTHM = cos(2Thetam)^2. For instance with a Graphite monochromator, CTHM = 0.8351 and 0.7998 for CuKb and CuKa respectively. Note that K is used only for synchrotron data (INSTRM=4, or ILOR=3)and does not have to be input for other kind of data: - for neutrons, LP = 1/ 2sin2Theta.cosTheta i.e. K is ignored - for characteristic X-Ray radiation (unpolarized beam) the formula used is: LP = [1 + CTHM.cos(2Theta)^2]/2(sinTheta)^2.cosTheta that corresponds to the general formula above for K=0.5 multiplied by 2.

- for synchrotron radiation : K must be given (K ~ 0.1) In the transmission geometry, flat plate with the scattering vector within the plate (Stoe geometry for X-rays,ILOR=2) the lorentz factor is: L = 1/sin2Theta and the polarization correction is given by: P = K (1+ c1 cos2Theta^2)/(1+c1) + (1-K) (1+ c2 cos2Theta^2)/(1+c2) where c1= sqrt(CTHM) = abs(cos(2Thetam)) c2= CTHM = (cos(2Thetam))**2 K = fraction of perfect crystal contribution Note that the program, internally, calculates the Lorentz factor combined with the absortion correction. A non-zero absorption coefficient (mt) must be given to get correct results (see section 3.7 to know the total correction for this case). ----------------3.5: Asymmetry ----------------Asymmetry correction (RLIM in line 4) : the shape of the peaks below RLIM is corrected using the semi-empirical function : A = 1 - P.sign(2Thetai-2Thetah). (2Thetai-2Thetah)2 /tan2Thetak where P is a refinable parameter. The above formula for asymmetry correction has been removed from the program. It was really very bad!. The expression given by Berar & Baldinozzi (J.Appl.Cryst, 26, 128 (1993)) has been introduced The asymmetry correction adopts now the form: A = 1 + {P1 Fa(z)+ P2 Fb(z)}/tan Theta + {P3 Fa(z)+ P4 Fb(z)}/tan 2Theta Where z = (2Thetai-2Theta-shift)/FWHM, (Shift includes the zero-point and other shifting terms) Fa(z)=2z exp(-z^2) , Fb(z) = 2(2z^2-3)Fa(z) and has four independent parameters, which are read at the same place as before. If NPHASE<0 the method proposed by C.J.Howard, in J.Appl.Cryst.15 615-620 (1982), is used. The approximation of the convolution integral is performed using the Bode's rule (Simpson formula for five points). The profile is calculated as a superposition of five profile functions (only pseudo-Voigt #5 and #7 are implemented for this correction) calculated at displaced points: Omega(x)={7g(x) + 32g(x+c2.P) + 12g(x+c3.P) + 32 g(x+c4.P) + 7 g(x+c5.P)}/90 Where x= 2Thetai-2Theta-shift, P is the asymmetry parameter and cotan(2Theta) = 16.c2 = 4.c3 = 16/9.c4 = c5

---------------------------3.6: Preferred orientation ---------------------------Two functions are currently implemented in the program : NORI = 0 --> The usual Rietveld function : Ph = (G2+(1-G2)*exp(G1*a^2) ) where G1 and G2 are refinable parameters and a is the acute angle between the scattering vector and the normal to the crystallites (platey habit). Note that setting G1 to any number > 99.0 for a phase causes the program to generate for that phase only those reflections for which d* is parallel to the preferred orientation vector PREF specified in line 11.2. NORI = 1 --> March's function Ph = G2 + (1-G2) * ((G1cosa)^2 + (sina)^2/G1)^-1.5 where G1 is a refinable parameter. This expression is adapted both to fiber and platey habits : G1 < 1 platey habit ( a is the acute angle between the scattering vector and the normal to the crystallites) G1 = 1 no preferred orientation G1 > 1 needle-like habit (a is the acute angle between the scattering vector and the fiber axis direction) Note that these values of G1 correspond to the Bragg geometry of usual X-ray powder diffractometer, for the Debye-Scherrer geometry of most neutron powder diffractometers the opposite holds. The parameter G2 represents the fraction of the sample that is not textured. The program put its value between 0 and 1 in case of divergence. ---------------------------------------3.7: Absorption (and microabsoprtion) ---------------------------------------Absorption correction : for Debye-Scherrer data, intensities may be corrected for the effects of sample absorption by applying the following transmission factor: Th = exp[-(1.7133 - 0.0368sinTheta^2)mR + (0.0927 + 0.375sinTheta^2)(mR)^2] where m is the linear absorption coefficient and R the radius of the cylindrical sample. For a flat plate in transmission geometry for X-rays (ILOR=2) the absorption correction is implemented as:

Th = exp(-mt/cosTheta)/cosTheta where t is the "effective" thickness of the sample. The product mt is given in the same place as mR. For Bragg-Brentano geometry in X-rays an angular-dependent microabsorption correction has been introduced following reference 20. The factor Th becomes: Th ---> Th Sr = Th (1-P) where Sr is given by the formula 6 of the first reference 20 (or P is given by formula 17 of the second reference 20). P = P0 * Cp * (Tau/sinTheta) * (1-Tau/sinTheta) The limitations and degree of applicability of the above formula are explained in reference 20. ----------------------------3.8: Systematic line shifts ----------------------------Systematic line position errors D2Theta: powder diffraction data are sometimes affected by systematic aberrations arising either from the sample itself or from an improper setting of the sample or diffractometer. FullProf gives the possibility to correct for two of the most commonly occurring errors by refining the parameters called SYCOS and SYSIN (line 10-1). These parameters relate to errors having a cosTheta and sin2Theta dependency, respectively. The corresponding errors originate from a different physical or/and geometrical problem depending on the diffraction geometry. They are summarized below: a) Bragg-Brentano parafocusing arrangement: the two largest systematic aberrations of Theta-2Theta powder diffractometers operating in this geometr y arise from specimen displacement and transparency; the sample displacement error is one of the largest systematic error affecting line positions in this geometry. It is given by: D2Theta = -2s/R cosTheta [in radians] where s is the displacement of the sample surface with respect to the axis of the goniometer and R the radius of the goniometer circle. The negative sign means that a displacement away from the center of the focusing circle moves the diffraction lines to lower 2Theta angle. The refinable parameter is SYCOS = -2s/R. This is by far the largest systematic aberration in this geometry. As the angles are expressed in degrees in FullProf, the sample offset can be calculated as: s = pi.R.SYCOS/180 The transparency correction is given by the relation: D2Theta = 1/(2mR) sin2Theta [in radians] where m is the linear absorption coefficient of the sample. This relation holds in the case of thick absorbing samples and the refinable parameter is SYSIN = 1/(2mR). For thin transparent samples, the correction would write:

D2Theta = t/R cosTheta

[in radians]

where t is the sample thickness; this (less usual) correction is not explicitly included in the code but can be accounted for by the displacement correction which turns out to show the same 2Theta dependency. Note however that samples requiring that kind of correction would also give biased integrated intensities; correction for this effect is not implemented in FullProf. For further details on systematic aberrations in Bragg-Brentano geometry, see [13]. b) Debye-Scherrer: the largest shifts of Debye-Scherrer rings result from sample off-centering and absorption. Eccentricity perpendicular to the incident beam direction is normally a second order effect if both sides of the Debye-Scherrer ring are measured. If only one side of the cone is measured, the line shift takes the form: D2Theta = e/R cosTheta [in radians] where e denotes the eccentricity, that is the refined parameter is SYCOS = e/R. Eccentricity in the incident beam direction is observed in both cases and takes the form: D2Theta = e/R sin2Theta [in radians] i.e., the refined parameter is SYSIN = e/R. The correction is negative for a shift along the beam direction towards the detector. For highly absorbing specimen with radius r, diffraction is limited to a cylindrical surface layer resulting in a maximum peak shift: D2Theta = r/R cosTheta [in radians] i.e., in this case, SYCOS = r/R. The latter effect also leads to disymmetrical line profiles. c) Curved position sensitive detector with flat plate sample: for the asymmetric geometry of diffractometers using a curved position sensitive detector (CPSD) with a flat-plate sample, the displacement correction takes the form: D2Theta = -s/(R.sina) sin2Theta [in radians] where R is the radius of the CPS and a the incident beam angle (in degrees) at sample surface. Thus, the parameter refined by FullProf is SYSIN = -s/(R.sina). The negative sign means, as in the case of Bragg-Brentano geometry, that a displacement away from the center of the focusing circle moves the diffraction lines to lower 2Theta angle; the value of the sample offset is given by: s = pi.R.SYSIN.sina/180 -----------------3.9: Codewords ------------------

Codewords Cx they are entered for each refined parameter. A zero codeword means that the parameter is not being refined. For each refined parameter, the codeword is formed as : Cx = sign(a).(10p + |a|) where p specifies the ordinal number of the parameter x (i.e. p runs from 1 to MAXS) and a(multiplier) is the factor by which the computed shift will be multiplied before use. The calculated shifts are also multiplied by a relaxation factor (line 5) before being applied to the parameters. ---------------------------3.10: Standard deviations ---------------------------Standard deviations are estimated from the formula : si = |a|.[Chi2*Mii]1/2 where a is the coefficient of the codeword for the parameter Mii is the corresponding diagonal element in the inverted matrix The Chi2 index used in the above formula is always calculated for the points in the pattern having Bragg contributions, thus si could be greater than the corresponding value calculated with other programs. As Chi2 is calculated in two ways (see below) the user can easily calculate the other value. ---------------------------3.11: Method of refinement ---------------------------Least squares refinement (the standard method) : Mp = Sum(i) { wi (yoi-yic)^2} the weights of the observations are calculated as : wi = 1/variance(obs)i Maximum likelihood refinement : the weights of the observations are calculated at each cycle as : wi = 1/variance(calc)i -------------------------3.12: Agreement factors -------------------------The quality of the agreement between observed and calculated profiles is measured by a set of nowadays conventional factors. In FullProf two sets of indices are calculated, according to the meaning of the integer N. In the first set N is the total number of points used in the refinement (N=NPTS-NEXC=total number of points in the patterntotal number of excluded points). In the second set only those points where there are Bragg contributions are taken into account. The definition of the indices is as follows: The profile Rp = 100Si|yoi-yci| / Si|yoi|

The weighted profile Rwp = 100[Siw|yoi-yci|2 / Siw|yoi|2]1/2 The Bragg RB = 100Sk|Ik-Ick| / Sk|Ik| The expected Rexpected = 100[(N-P+C)/Si(wi yoi2)]1/2 The goodness of fit Chi2 = [Rwp / Rexpected]^2

where N-P+C is the number of degrees of freedom (The meaning of N has been given above,P the number of refined parameters and C the number of strict constraint functions). Excluded regions are always excluded from the calculation of all agreement factors. Conventional Rietveld R-Factors: cRp, cRwp are calculated as above but using background corrected counts. The magnetic R-factor is defined as the Bragg RB-factor but is applied to magnetic intensities. The "observed" integrated intensity Ik is in fact calculated from the Rietveld formula: Ik = Ick Si { Omeg(Ti-Tk) (yoi-Bi)/(yci-Bi) } This formula is equivalent to a "proportional sharing" of the integrated intensity of a cluster between its components according to the actual model. Then, if the model contains a strictly zero integrated intensity for the component "k" (Ick=0), the observed integrated intensity is also zero: Ik=0, even if it is obvious that Ik is not zero from the experimental pattern. This has as a consequence that the reflections with Ic=0 do not contribute to the Bragg R-factor. Although commonly used in crystallography, the Rp,Rw, RB agreement factors are not satisfactory from a statistical point of view. Therefore, a number of statistically more significant parameters are calculated by FullProf: a) the deviance [15] defined as: D = 2 Sum(i){ yio ln(yio/yic) - (yio-yic)} b) from the deviance, one can derive two other measures of discrepancy which are useful as model selection criteria (somewhat analogous to Hamilton's criterion). These criteria take account of both the goodness of fit of a model and of the number of parameters used to achieve that fit. They take the form: Q = D + a.maxs where maxs is the number of refined parameters and a represents the "cost" of fitting an additional parameter. Akaike's information criterion uses a=2 while Schwarz's criterion has a=ln(maxs). c) the Durbin-Watson statistic parameters: d and Q. quantities to assess the quality of the refinement Hill and Flack [16]. This statistic which measures adjacent residuals (serial correlation) is defined The use of these two has be advocated by the correlation between as:

d= Sum(i=2toN){[wi(yi-yic)-w(i-1)(y(i-1)-y(i-1)c)]^2}/Sum(i){[wi (yi-yic)]^2} Serial correlation is tested (at the 99.9% confidence level) by comparing

the value of d to that of Q which is given by the relation: Q = 2 { (N-1)/(N-P) - 3.0902/Sqrt(N+2) } Three cases may occur: - if d < Q, there is positive serial correlation: successive values of the residuals tend to have the same sign. This is the most common situation in profile refinement. - if Q < d < 4-Q, there is no correlation - if d > 4-Q, there is negative serial correlation: successive values of the residuals tend to have opposite sign. --------------------------------3.13: Analysis of the refinement --------------------------------An emprical analysis of the refinement has been introduced at the end of the file CODFIL.SUM. A part from some subjective comments that could appear in that analysis there are some important quantities that have to be known by the users because they have not been published yet. - Expected Rp factors are calculated supposing the best possible model. Rp = Sum abs(Yiobs-Yical)/Sum(Yiobs), where Yical is calculated with the help of a Poissonian function, Genpoi(x), from Yiobs. The argument of Genpoi is an integer value representing an observation which is equal to their variance, Genpoi returns another possible value compatible with the deterministic value x of variance x. The value of x is calculated form Yiobs and Variance(Yiobs) as follows: c1=Yiobs/Var(Yiobs) x=Yiobs*c1 Yical=Genpoi(x)/c1 - The percentage of the contribution to the total integrated intensity (Sum(Iobs) for all phases) of each phase is now written. - The number of refined parameters distributed in three classes are written: Nglb : Number of global parameters (not depending of the phase index) Nprofp: Number of profile parameters. Nintdp: Number of intensity-dependent parameters (x,y,z,B,occ, Mx...) The preferred orientantion parameters are included in this class. - An effective number of reflections is calculated in order to get the ratio Refni=(Effective Number of Reflections)/Nintdp = Enref/Nintdp. - The concept of effective number of reflections is introduced in order to take into account the effect of the resolution in the refinement. It is clear that well separated independent reflections give better results that when the reflections are overlaped. A global effective number of reflections is calculated by the program. For each phase, a similar indicator is written. A reflection contributes as x/(x+nearest), where "x" is the fraction of the total area of the current phase and "nearest" is the number of adjacent reflections verifying the formula:',

2theta-p*FWHM <= 2theta(adjacent) <= 2theta+p*FWHM "nearest" is weighted by the corresponding "x(s)", and p is a parameter lower than unity. The general formula for calculating the global effective number of reflections is: Enref = Sum(i){ x[ph:i]/Sum(ni){x[ph:ni]} } where Sum(i) is a sum over all the reflections contributing to the allowed areas of the diffractogram. Sum(ni) is the sum extended to the reflections near to the reflection "i" (including this reflection). The symbol x[ph:i] is the "x"-value of the phase to which the reflection "i" belongs. The same formula restricted to reflections of a single phase is applied to calculate Enref(Iphase): Effective number of reflections of the phase Iphase. The program calculates these numbers (Enref's) and the ratios (Refni's) for three values of the parameter p (1, 1/2 and 1/4).

4.- ADDITIONAL NOTES

--------------------------4.1: Magnetic Refinements --------------------------For a conmensurate structure two descriptions of the magnetic structure are possible: the general formalism using fourier components of magnetic moments through the propagation vectors or a crystallographic-like description using the magnetic unit cell. The most simple one, for a non-expert user, is to describe the magnetic structure in the magnetic unit cell. In that case two important points must be taken into account: a) For magnetic structures described in a magnetic unit cell larger than the crystallographic cell, the coordinates of atoms must be changed consequently, as well as the unit cell parameters. The codes for common (or related) parameter with the crystallographic counterpart, must be changed applying a the correct multiplicative factor. It is worth stressing that codes for cell parameters are actually applied to the "cell constants" defined by: 1/d2 = Ah2 + Bk2 + Cl2 + Dkl + Ehl + Fhk Therefore, if you are dealing with an orthorhombic structure which has a magnetic structure with propagation vector k=[1/2 0 0] you have to use the magnetic unit cell 2xa, b, c and if the codes of the crystallographic unit cell a,b,c are for example: Cell CodeCell a 81.00 b 91.00 c 101.00 90 0 90 0 90 0

the corresponding values for the magnetic counterpart are:

Cell CodeCell

2a 80.25

b 91.00

c 101.00

90 0

90 0

90 0

the reason is that the crystallographic cell constant A is in this case: Ac=1/a**2 and the magnetic cell constant is Am=1/(2a)**2= 0.25*Ac. As explained above, the multiplicative factor is applied to shifts of parameters. b) The scale factors of the crystallographic and magnetic parts have to be related in some way in order to get good values of magnetic moments. This relation depends on the way the users describes the magnetic structure, however several rules can be useful to avoid bulk errors: i) Use the correct occupation numbers in the crystallographic part (=multiplicity of special position/general multiplicity) ii) The number of magnetic atoms in the chemical unit cell must coincide with the description above, therefore the occupation numbers in the magnetic part are related to the number of symmetry operators given, the centrosymmetry (or not) of the magnetic structure. If these requirements are satisfied and the magnetic unit cell is the same as the chemical one, the scale factors are strictly the same numbers; and, therefore represents the same parameter with a shift equal to unity. If the magnetic unit cell is a multiple one, and the above requirements are satisfied, the relation between the scale factors is a multiplicative factor given by: Sc: crystallographic scale factor Sm: magnetic scale factor Sm = 1.0/(Volm/Volc)**2 Sc In the case of large magnetic cells it can be more convenient to modify the occupation numbers of magnetic atoms in such a way that the two scale factors coincide. For incommensurate magnetic structures the general formalism must be applied. When the magnetic structure is described using the formalism of propagation vectors, the components Mx,My,Mz no longer represent true magnetic moments. The user should be cautious in interpreting the output files. The modulus of the "magnetic moment" represent the Fourier component modulus of an atomic magnetic moment which have to be calculated externally. The calculation of the intensity is based on the expression of magnetic structure factor given in mathematical section, therefore the user knows how to play with his input items in order to obtain physically sound results. For spherical description of the magnetic moments the following must be taken into account: The orthonormal system with respect to which are defined the spherical angles verifies: X axis coincides with the crystallographic A Y axis belongs to the plane A,B Z axis is perpendicular to the plane A,B

The particular implementation of spherical components in magnetic structure refinements is that Z axis must coincide with C. That works in all crystallographic systems except for triclinic. The monoclinic setting must be changed to -> 1 1 2/m to satisfy the above prescription. -------------------------4.2: Propagation vectors -------------------------A complete list of reflections can be generated when propagation vectors of an incommensurate structure are present. To each fundamental reflection it is added the corresponding satellites. For n propagation vectors k1, k2,...kn, there are n satellites obtained from each fundamental reciprocal lattice vector h: Fundamental reflection: h=[h,k,l] Satellites h1=h+k1, h2=h+k2,...hn=h+kn In the present version of the program no symmetry analysis is performed. We recommend to use the triclinic space group L -1 (where L=P,A,B,C,F,I,R) in order to have a full set of reflections with the proper multiplicity when the true magnetic symmetry is not known. The program generates first a list of unique reflections corresponding to the required space group and then adds the satellites. This method had to be modified for reflections belonging to the boundary planes and lines of the asymmetric region of the reciprocal space in order to obtain the correct number of reflections and not miss (or repeat) some of them. Be careful with propagation vectors k equivalent to -k!. Two vectors k1 and k2 are "equivalent" if k1-k2 is a vector of the reciprocal lattice. So, for k2=-k1, if K=2*k1 belongs to the reciprocal lattice, k1 is "single" and belongs to a point of high symmetry of the Brillouin Zone. In such cases only ONE propagation vectors should be introduced NVK=1, if the user puts NVK=-1, the satellite reflections are not correctly generated. For centered cells a propagation vector k having components +- 1/2, verifies that 2k has integer components, but that does not mean that k and -k are equivalent, because 2k could not belong to the reciprocal lattice. For a C lattice the propagation vector k1=(1/2 0 0) is not equivalent to k2=(-1/2 0 0) because K=2*k1=(k1-k2) K=(1 0 0) does not belong to the reciprocal lattice: h+k=2n is the lattice C condition for components (hkl). On the contrary, the vector (0 0 1/2) is "single" because (001) is a reciprocal lattice point of the C lattice. ------------------------------------------------------------------------------4.3: Microstrains and domain size effects. HKL-dependent shifts and asymmetry. ------------------------------------------------------------------------------In real materials it is frequent to find more or less strong microstructure effects in the shape and width of Bragg reflections. The capabilities of the program in order to handle these effects is reflected in two subroutines (STRAIN and SIZE) which calculate the broadening of Bragg reflections as a function of hkl and a set of microstructural parameters depending defining a model. In general the broadening of reflections due to microstrains has an angular

dependence of the form: DFWst = e tan(Theta). The corresponding dependence for the size effect is given by the Scherrer formula DFWsi= c/cos(Theta). The reader i s referred to the literature for further details. An increase of the U halfwidth-parameter with respect to the instrumental value is indicative of the presence of microstrains. The value (U-Uins) is an estimate of the isotropic broadening. The parameter IG (line 11-6-1) is a measure of the isotropic size effect. strain(gaussian) = (pi/1.8).Sqrt(U - Uins) (in %)

size(gaussian) = 180 K Lambda / (pi.sqrt(IG)) (Lambda & size in A) where K denotes the Scherrer crystal shape constant near to unity. If the user tries to obtain the maximum physically significant information on the microstructure of his(her) sample, the profile function NPROF=7 must be used. Only in that case can be performed a correct convolution of the instrumental function with the intrinsic profile function (considered both as Voigtians). In the case NPROF=7, isotropic lorentzian and gaussian contributions of size and microstrains are represented by the following parameters: Gaussian component Size: Strain sqrt(IG) sqrt( U) Lorentzian component Y X

The FWHM of the gaussian and lorentzian components of the Voigt function for NPROF=7 is calculated using the expressions: FWHM^2(gaussian) = [U + c.DST(STR)^2] * Tan^2(Theta) + V * Tan(Theta) + W + IG/cos^2(Theta) FWHM(lorentzian) = X tan(Theta) + (Y + F(SZ))/cos(Theta) The isotropic particle size broadening Y provides a size value given by: size(lorentzian) = 180 K Lambda /(pi.Y) (Lambda & size in A)

K = FWHM/Integral-Breadth = 2/pi for a lorentzian (with such value for K, size represent the volume averaged diameter of crystallites in all directions) For anisotropic contributions {DST(STR) and F(SZ)} the actual version of FullProf has the following set of microstructural models: Size effect:{only lorentzian component is taken into account for anisotropic broadening of size origin through F(SZ)} ISTR=0 and IsizeModel /= 0 {Integer to select a particular model for F(SZ) in subroutine SIZE} IsizeModel = 1 Platelet coherent domains. F(SZ) is assumed to be of the form F(SZ)=SZ*cos(phi), where SZ is the refined parameter and phi is the acute angle between the scattering vector (h,k,l) and the vector defining the platelet shape of domains {[Sz1,Sz2,Sz3] in line 11-9*} Needle-like coherent domains. F(SZ) is assumed to be of the form F(SZ)=SZ*sin(phi), where SZ is the refined parameter and phi is the acute angle between the scattering

IsizeModel =-1

vector (h,k,l) and the vector defining the needle shape of domains {[Sz1,Sz2,Sz3] in line 11-9*} IsizeModel = 2 to 7, in these cases the broadening is considered only for reflections of the form: (00l),(0k0),(h00), (hk0), (h0l) and (0kl). The function is: F(SZ)=SZ. IsizeModel =8 Only satellites reflections are considered to be broadened. In this case also F(SZ)=SZ. IsizeModel =9 IsizeModel =10 IsizeModel =11 IsizeModel =12 IsizeModel =13 Reflections (HKL) with H=2n+1 and K=2m+1 are broadened. In this case also F(SZ)=SZ. Reflections (H0L) with H+L=2n are broadened. In this case also F(SZ)=SZ. Reflections (HKL) except (HHL) are broadened. In this case also F(SZ)=SZ. Reflections (HKL) with H=2n+1 are broadened. In this case also F(SZ)=SZ. Reflections (HKL) with H=2n and K=2m+1 are broadened. In this case also F(SZ)=SZ.

ISTR > 1 : A generalized formulation of the size effect is applied and the size parameters (up to six) are given in line 11-14*. In that case the F(SZ) function is given by: F(p1...p6) = 180 Lambda/(pi^2) * d(hkl)^2* Quad(p1...p6) Quad = p1 h^2 + p2 k^2 + p3 l^2 + 2 (p4 hk + p5 hl +p6 kl) where the size parameters are p1,p2,...p6. This formulation assumes that the average diameter can be expressed as an ellipsoid of refinable parameters p1..p6. Therefore = d(hkl) * Quad Strain effects:{only gaussian component is taken into account for anisotropic broadening due to strains through DST(STR)} ISTR=0 and IstrainModel /=0 {Integer to select a particular model for DST(STR ) in subroutine STRAIN} FWHM(strain) = DST tanTheta = c Sigma(1/d^2)/(1/d^2) tanTheta the constant c is equal to 2 sqrt(Ln2) (180/pi). In the program the effective value used is c*10e-3, therefore it is necessary to divide the refined strain parameters by 1000. The user can consult the modules FDUM's.FOR for details and reference 17. IstrainModel = 1 IstrainModel = 2 IstrainModel = 3 Orthorhombic strain in a tetragonal lattice. Refined parameter STR1. Strain along "a" in an orthorhombic lattice. Refined parameter STR1. Strain along a diagonal of the ab plane in an orthorhombic lattice for which an equivalent monoclinic double cell has

been used. Refined parameter STR1. IstrainModel = 4 IstrainModel = 5 Same as 3 but strains along the two diagonals. Refined parameters STR1, STR2. Strain in an orthorhombic lattice with fluctuations along "a" and "b" and correlation between "a" and "b". Refined parameters STR1, STR2, STR3. Strain in an orthorhombic lattice with fluctuations along "a", "b" and "c" and corr(a,b)=-1, corr(a,c)=1, and corr(b,c)=0. Refined parameters STR1, STR2, STR3. Uniaxial microstrain. DST(STR) is assumed to be of the form DST(STR)=STR1*cos(phi), where STR1 is the refined parameter and phi is the angle between the scattering vector (h,k,l) and the vector defining the axial direction of strains {[St1,St2,St3] in line 11-10*} General anisotropic strain of hexagonal symmetry. (See ref. 17)

IstrainModel = 6

IstrainModel = 7

IstrainModel = 8

STR1= Saa=Sigma(A), STR2= Scc=Sigma(C) STR3= Cac=Correl(A,C) 1/d^2 = SQ = A (h^2+k^2+hk) + C l^2 = A m^2 + C l^2 Var(SQ)= Saa^2 m^2 + Scc^2 l^4 + 2 Saa Scc Cac m l^2 FWHM^2 (strain)= Var(SQ) (P/SQ)^2 {P=180/pi*sqrt(8Ln2) 1.E-03 tan(theta)} IstrainModel = 9 Strain adapted for compounds LaNi5-like. Formula by: P.Thompson et al, J.Less Comm Met 129, 105-114 (1987)

Four param.: Gaussian FWHM=Sqrt(Sumj[S(j)]) tan(theta) STR1 STR2 STR3 STR4 -> -> -> -> S(1) S(2) S(3) S(4) = = l^4 /(h^2+k^2+l^2) = (h^2 k^2+k^2 l^2)/(h^2+k^2+l^2) = h^2 k^2 /(h^2+k^2+l^2) x x x x Par(1) Par(2) Par(3) Par(4)

IstrainModel =10

General anisotropic strain for orthorhombic symmetry. (See ref. 17) 6 strain parameters corresponding to :

STR1=Saa=Sigma(A), STR2=Sbb=Sigma(B), STR3=Scc=Sigma(C), STR4=Cab=Corr(A,B), STR5=Cac=Corr(A,C), STR6=Cbc=Corr(B,C) Where: SQ = A h^2+ B k^2 + C l^2 Var(SQ)= Saa^2 h^4 + Sbb^2 k^4 + Scc^2 l^4 + 2 Saa Sbb Cab h^2 k^2 + 2 Saa Scc Cac h^2 l^2 +

2 Sbb Scc Cbc k^2 l^2 FWHM^2= Var(SQ)(P/SQ)^2 {P=180/pi*sqrt(8Ln2) 1.E-03 tan(theta)} IstrainModel =11 General anisotropic strain for monoclinic symmetry being the setting gamma/=0 (2-fold axis along 001) (See ref. 17) {limited to 8 parameter by putting arbitrarily zero two correlation values (Cad=Cbd=0)} 8 strain parameters corresponding to : STR1=Saa=Sigma(A), STR2=Sbb=Sigma(B), STR3=Scc=Sigma(C), STR4=Sdd=Sigma(D), STR5=Cab=Corr(A,B), STR6=Cac=Corr(A,C) ,STR7=Cbc=Corr(B,C), STR8=Ccd=Corr(C,D) Where: SQ = A h^2+ B k^2 + C l^2 + D hk 0 } 0 } Var(SQ)=( h^2, k^2,l^2,hk){ } { SaaSccCac SbbSccCbc Scc^2 SccSddCcd} { 0 0 SccSddCcd Sdd^2} FWHM^2= Var(SQ)(P/SQ)^2 {P=180/pi*sqrt(8Ln2) 1.E-03 tan(theta)} ISTR=1,3 The generalized formulation of the strain broadening is used and 10 parameter are read in line 11-13*. The formalism is that described in reference 17. The strain parameters correspond to the fluctuations and correlations of direct cell parameters. The lower symmetry is monoclinic with the conventional setting. The hexagonal lattice is treated as a special case. Strain parameters: Px = sigma(x), corr(x,y) = sin(Pxy) x,y = a,b,c, beta { Saa^2 SaaSbbCab SaaSccCac { SaaSbbCab Sbb^2 SbbSccCbc (h^2) (k^2) (l^2) (hk)

The parameter Pxy is an angle in degrees. This particular description is for taking into account that corr(x,y)=covar(x,y)/(sigma(x).sigma(y)) and must belong to the interval [-1,1]. The strain parameters are read in the following order: Pa, Pac, Pb, Pabeta, Pc, Pbc, Pbeta, Pbbeta, Pab Pcbeta

Note: ----In some cases the use of size-strain options of the program, the computing time is greatly increased. This is related to the fact that the number of reflections contributing to each particular point of the diagram increases because the broadening is quite important. An excessive computing time may indicate a divergence of the refinement that has attributed a too much large FWHM to some reflections.

-----------------------------------HKL-dependent shifts and asymmetry. -----------------------------------As a part of microstructural effects, the program can handle some cases of peak-shifts and asymmetry effects. No hkl-dependent asymmetry model is currently available. Two models are available for peak-shifts: ISHIF = 1/-1 Uniaxial shifts along a direction [S1,S2,S3] Shift of Bragg reflections of the form: Shift(hkl) = Shift1 * cos(Phi(hkl)) (Ishif= 1) or Shift(hkl) = Shift1 * sin(Phi(hkl)) (Ishif=-1) Where Phi is the angle: between [hkl] and [S1,S2,S3] ISHIF = 2 Shift(hkl) = Shift1 * Coeff(hkl) * tan(Theta(hkl)) The parameter Coeff(hkl) must be given by user in the file from which Bragg reflection indices and multiplicities are read. This option cannot be used simultaneously with the output of a Fourier file, so JFOU is set to zero y the program. -----------------------------------4.4: Quantitative phase analysis -----------------------------------For quantitative analysis it is essential that two conditions are fulfilled: - sample must be carefully prepared to comply to the definition of a powder: homogeneity, sufficient number of particles with random orientation - structure factors must be correctly calculated According to Brindley, it is convenient to classify mixed powders according to the value of the product mD where m is the linear absorption coefficient and D a measure of the linear size of a particle. Four cases must be considered: - fine powders: mD< 0.01 The individual particles of the powder have negligible absorption and no correction has to be applied to the data - medium powders: 0.01 < mD < 0.1 - coarse powders: 0.1 < mD < 1 - very coarse powders: mD > 1 In a mixture of N crystalline phases the weight fraction Wj of phase j is given by: Wj ={ Sj Zj Mj Vj / tj} / Sum(i)[Si Zi Mi Vi /ti] where Sj Zj Mj Vj tj is the scale factor of phase j is the number of formula units per unit cell for phase j is the mass of the formula unit is the unit cell volume is the Brindley particle absorption contrast factor for phase j defined as:

tj = (1/Vj) Integ[ exp{-(mj-mu)x} dVj] where Vj is the volume of a particle of phase j mj is the particle linear absorption coefficient mu is the mean linear absorption coefficient of the solid material of the powder x is the path of the radiation in the particle of phase j when reflected by the volume element dVj

The latter parameter accounts for microabsorption effects that become important when the compounds of the powder have rather different linear absorption coefficients. Its calculation requires only the knowledge of the particle radius R and linear absorption coefficient m. Values of t as a function of the product (mj-mu)R have been tabulated by Brindley and are reproduced below: R t t t R t t t (mj-mu) (2Th=0) (2Th=90) (2Th=180) (mj-mu) (2Th=0) (2Th=90) (2Th=180) -----------------------------------------------------------------------------0.50 2.068 2.036 2.029 -0.40 1.813 1.807 1.827 -0.30 1.508 1.573 1.585 -0.20 1.352 1.353 1.362 -0.10 1.159 1.162 1.163 -0.09 1.142 1.143 1.144 -0.08 1.124 1.125 1.125 -0.07 1.107 1.108 1.108 -0.06 1.090 1.091 1.091 -0.05 1.074 1.073 1.074 -0.04 1.059 1.058 1.059 -0.03 1.043 1.042 1.042 -0.02 1.028 1.027 1.027 -0.01 1.014 1.014 1.014 0.00 1.000 1.000 1.000 0.01 0.986 0.986 0.986 0.02 0.972 0.973 0.973 0.03 0.959 0.960 0.960 0.04 0.945 0.946 0.947 0.05 0.932 0.933 0.934 0.06 0.918 0.919 0.921 0.07 0.905 0.906 0.908 0.08 0.892 0.893 0.895 0.09 0.878 0.879 0.882 0.10 0.865 0.866 0.870 0.20 0.742 0.753 0.760 0.30 0.640 0.653 0.671 0.40 0.545 0.569 0.587 0.50 0.468 0.496 0.529 R.J. Hill & C.J. Howard, J. Appl. Cryst. 20, 467-476 (1987) G.W. Brindley, Phil. Mag. 36, 347-369 (1945) ------------------------------------------------------------------------------4.5: User-supplied parameters & subroutines for structure-factor calculation. ------------------------------------------------------------------------------Some powder diffraction applications need user-supplied parameters and subroutines calculating the structure factors. Examples of this kind of applications are the following: Form factor determination. Rigid body-generalized coordinates refinements. TLS refinements. Description of incommensurate phases in real space (not Fourier components) by some simple parameters. - Anharmonic displacement parameters. To do that, the user can prepare his(her) own subroutines and link the object FDUM1,FDUM2,FDUM3,... with the rest of the modules and the FullProf library. For crystal structures the subroutine to calculate the structure factor

squared for a reflection (hkl) must be called STRMOD. For a magnetic phase the name is MAGMOD. Examples of these are contained in the public source code FDUM's.FOR. The user can also modify or introduce new models for strains and size effects in the subroutines SIZE and STRAIN which are also included in FDUM's.FOR Useful information for users concerning the internal parameters: The parameters (refinable variables) used in FullProf are divided in three categories: Global parameters, Atom parameters and Phase parameters. For each category there is an array to store the values: GLB, XL and PAR. The global parameters are the following: GLB(J) : J=1,2....40 J=1 : T0 J=2,...7: b1,b2,...b6 J=8,9 : SYCOS,SYSIN J=10..15: Bc1,.....Bc6 J=16..21: d1,..... d6 J=22,23,24,25: a,b,c,d J=26 : Flambda OF J=27,28,29: P0,Cp,Tau J=30..35: J=35..40: Microabsorption coefficients Additional background parameters Not used at present Zero Background parameters Systematic 2Theta-shifts/ DTT1,DTT2 for TOF Background coefficients for Debyelike function: DBback distances in Angstroms in DBback Background transforming coefficients Wavelength (for CW neutron)/ 2SinTheta for T

XL(I,J): Atom parameters, the index I runs over atoms I=1,2,..N1,N1+1,N1+2,....N1+N2,N1+N2+1,....N1+N2+N3....N N1: number of atoms of phase 1 N2: number of atoms of phase 2 ............................................... N : total number of atoms : N1+N2+N3+...+NNPHASE The index J represents different atom parameters J=1,2,3 : J=4 : J=5 : J=6....8: 9...11: J=12 : J=13..18: J=19..32: (x,y,z) B Occ. (Mx,My,Mz) (Mxi,Myi,Mzi) MPhas (betas11,22..) fractional coordinates isotropic temperature factor Occupation factor. Magnetic moment components Imaginary Magnetic moment components Magnetic phase of the atom Anisotropic thermal factors Coefficients of generalized form-factors

These atomic parameters can be used by the subroutines STRMOD and MAGMOD for other purposes. Moreover, there are other non-atomic parameters which can be used also for defining generalized coordinates. They are internally stored in the array PAR(K,J) where the index K run over the phases and J represents different parameters: PAR(K,J):K=1,2,....NPHASE

J=1 J=2 J=3,4,5 J=6,..11 J=12,13 J=14 J=15,16 J=17 J=18..27 J=28 J=29..34 J=35,36,37 J=38 J=39,40 J=41,42 J=43,44,45 J=46,47,48 J=49,50 J=51,..59 oad.)

: : : : : : : : : : : : : : : : : : :

S Bov U,V,W A,B,C,D,E,F G1,G2 As1 X,Y eta(m) Str1...Str10 IG Siz1...Siz6 As2,As3,As4 Asymhkl Shf1,Shf2 Shap1,Shap2 U2,V2,W2.. Ud,Vd,Wd Eta0d, Xd CYi

Scale factor Overall temperature factor Halfwidth parameters Cell constants Preferred orientation parameters First Asymmetry parameter Shape parameters (depend on NPROF) Shape parameter of p-Voigt(Pearson) Strain parameters Isotropic gaussian size parameter Anisotropic size parameters Additional asymmetry parameters HKL-dependent asymmetry parameter HKL-dependent shift parameters Additional shape parameters U,V,W parameters for lambda(2) U,V,W for right part of SpV(x) Eta and X for right part of SpV(x) Coeffs. of Spherical Harmonics (size br User-supplied parameters Number of the Last parameter used internally. To be used as offset for user-supplied subroutines.

J=60..MPAR : P1,P2,.... LASTMPAR=59

---------------------------------------------------------------------------4.6: Examples of user-supplied subroutines for structure-factor calculation ---------------------------------------------------------------------------CRYSTAL STRUCTURE REFINEMENTS In FullProf the default JBT=4 option is an updated subroutine written by V.Rodriguez which handles rigid body objects.

Guide for the General Rigid-Body-Constraints/TLS Subroutine for Rietveld Refinements ********************** V1.4 (January 1997) Dr V.Rodriguez Laboratoire de Spectroscopie Moleculaire et cristalline * (U.A.124 CNRS) Universite de Bordeaux I, 351 cours de la Liberation F-33405 TALENCE Cedex Tel : (33) 56.84.63.61 Fax : (33) 56.84.84.02 * UMR Physico-Chimie Moleculaire since January 1st 1997 For further details see "A routine for crystal-structure refinements based on rigid-body model with constrained generalized coordinates and mean thermal

displacements" Rodriguez V. and Rodriguez-Carvajal J., J. Applied Cryst., (to be published) ********************* In the present version of Fullprof the default JBT=4 option refers to a subroutine which handles rigid body objects. Here we give a short guide in order to facilitate the use of this subroutine. This routine should be used with caution and the user should be familiar enough with "conventional" Rietveld refinements. The subroutine is still under improvement. At this time, all the options have been well tested except option Nr 4 (see below). I acknowledge all suggestions and notifications of possible bugs found in the program. It should be pointed out that option JBT=4 is not restricted to "perfect" rigid-body-groups. In the standard definition of Rietveld atomic positions, reduced coordinates are used. Within this subroutine, molecules or atomic groups are defined by spherical internal coordinates and six additionnal parameters, which define the groups position and orientation in the crystal. In addition, this extended version of Rietveld refinement allows to distort groups in a final refinement step if required. Note: Each parameter refering to recognized by the appended Each parameter refering to recognized by the appended A-Description of the parameters ******************************* - Any kind of rigid body groups built from atoms can be generated. - All parameters described below are refinable except optional parameters P6 (always true) and P16 (specific to "satellite" rigid body groups). - Each group of atoms (rigid body group, RBG) is identified by 1 or 2 letters and a number from 1 to 99 (both compatible with the standard DBWS format and the free format of Fullprof). The number indicates that a RBG is defined. The label(two letters max.) and the number of the atom must not be separated. - Representation of RBG is performed as follow : 1) Definition of the isolated group in an orthonormal molecular system. Each atomic position (xm,ym,zm) is defined through spherical coordinates (dm=distance, thetam, phim). The centre of this molecular system may or may not coincide with an atom. 2) Definition of the absolute position (Xc,Yc,Zc) and orientation (angles THETAc, PHIc and CHIc) of the molecular orthonormal system in the crystallographic orthonormal system. These are standard EULER angles definitions. In this way, each group is entirely defined within a crystal. - Anti-clockwise rotations are always applied to spherical angles. These standard spherical angles are defined as follow : THETAc / thetam : inclination with respect to the Zc/zm axis of the the orthonormal crystal system is "c". the orthonormal molecular system is "m".

corresponding orthonormal system [-Pi, Pi]. PHIc / CHIc and phim: rotation around the Zc/zm axis of the corresponding orthonormal system [0, 2*Pi]. - The orthonormal crystallographic system, with respect to which the spherical angles THETAc, PHIc and CHIc of any RBG are defined, fulfills the following conditions: Xc axis coincides with the crystallographic direction a Yc axis belongs to the plane (a, b) Zc axis is perpendicular to the plane (a, b) (parallel to c*) One can imagine how to place a rigid molecule in the correct position in the unit cell by making the following operations: - The RBG has been completely defined by the spherical coordinates (dm,thetam,phim) of each atom in the internal orthogonal system, that coincides at the beginning with the orthogonal crystallographic system. - Perform a rotation CHIc of the whole RBG around the zm axis of matrix: ( COS(CHIc) M1= ( SIN(CHIc) ( 0 -SIN(CHIc) COS(CHIc) 0 0 0 1 ) ) )

- Perform a rotation of the wole RBG inclining their z-axis to the angles THETAc and PHIc of matrix: ( COS(THETAc)*COS(PHIc) M2= ( COS(THETAc)*SIN(PHIc) ( -SIN(THETAc) -SIN(PHIc) COS(PHIc) 0 SIN(THETAc)*COS(PHIc) ) SIN(THETAc)*SIN(PHIc) ) COS(THETAc) )

- Translate the origin of the RBG to the position (Xc,Yc,Zc) which are transformed in reduced coordinates (xo, yo, zo) with xo=Xc/a, yo=Yc/b, zo=Zc/c - "Free atoms" (unconstrained or isolated) can be added with number 0 or without number. Each isolated atom has the same definition of parameters as described in the Fullprof manual (see lines 11.4.1 and 11.4.3 of the Fullprof guide). Current information about "free atoms" is printed to the screen when parameter P6 (see below) is set to any negative value. Note that the following parameters: P6, P7, P8, P15 and P16, are obsolete in this case. - Rigid body satellite groups (RBSG) can be also included in this version, for example a methyl group within a rigid group such as [N(CH3)4]+(tetramethyl ammonium). The definition and the structure of the parameters are almost the same as those for a main RBG. The coordinates of the centre of the satellite group should not be specified since they are specified through the knowledge of its absolute position in the input file "x.pcr". The orientation THETAc,PHIc of satellites groups is also defined following the value of parameter P16 of the 1st satellite atom (see below, option abs(P6)=2). One easy external degree of freedom of the RBSG is the rotation around the z-molecular axis i.e. the N-C bond. This degree of freedom is accesible through the EULER angle PHIc that can be refined. - All output is identical for isolated and/or constrained atoms.

- Every atom has 16 items stored in (X(Iphase,j),j=1,16) and each parameter is denoted here as P1, P2, ...P16. The reading of these parameters is described in lines 11-4-1 and 11-4-3 of the Fullprof guide. Example: P1 P8-PHIc TD1 N P2 P3 P4 P5 P6 P7-THETAc -.290 71.00 .142 61.00

.38900 .54090 .28796 4.00000 0.16667 1.000 .00 .00 .00 .00 .00 .00 1.49300 1.57080 -.95532 .33333 .66667 .25000 1.04100 0 .00 .00 .00 .00 .00 .00 51.00

P9-dm P10-thetam P11-phim P12-xo P13-yo P14-zo P15-CHIc P16 x y z TD2 C .47744 .84054 .28796 4.00000 0.16667 .000 .000 .00 .00 .00 .00 .00 .00 .00 .75000 1.57080 1.57080 .00000 .00000 .00000 .00000 .00 .00 .00 .00 .00 .00 .00 P9-dm P10-thetam P11-phim O O .03957 .00 beta11 x y .28146 -.04125 .00 .00 .07078 .09069 .00 .00 beta22 beta33 z .41456 .00 .00314 .00 beta12 .00000 1.00000 -1.000 .00 .00 .00 .01571 -.01501 0 .00 .00 beta13 beta23 .000 .00

.000 .00

.000 .00

-> In a RBG (here, the generic name is TD) each atom (TD1 and TD2) has its internal coordinates stored in the following items: distance to the centre dm : P9 Spherical angle thetam : P10 Spherical angle phim : P11 The atomic reduced cordinates are only for information with option 1 (RBG), see below for details. Reduced coordinate x : P1 Reduced coordinate y : P2 Reduced coordinate z : P3 The parameters P4 and P5 have their usual meaning for each atom: Isotropic temperature factor: P4 Ocupation number : P5 -> The first atom of each RBG (atom No 1, here TD1) contains the fractional coordinates (xo, yo, zo) of the centre of the RBG (P12, P13, P14) and the three EULER orientation angles THETAc, PHIc and CHIc (in radians, respectively P7=-.290, P8=.142 and P15=1.041) of the whole group. Here, the centre of the RBG does not coincide with an atom (N or C in this example) since the first atom has a non zero value dm=P9=1.493. Of course, it is possible to build a RBG with an atom coinciding with the centre of the RBG. -> The third atom (oxygen, Number 0 (or none which is equivalent to 0)) is

unconstrained ("free atom") and its items are the standard ones : here, reduced coordinates (x, y, z) + atomic temperature parameters (betaij). -> No terminal output is required for the whole RBG TD (P6>0). Terminal output at each cycle is required for "free atom" oxygen (P6<0). B-Practical details ******************* 1- The item abs(P6) of the first atom indicates the option selected. If P6<0, the current information is printed on the screen and in the file CODFIL.OUT. (file x.out). If the value of P6 is 0 for a RBG, the parameter defaults to "1", corresponding to the standard RBG option. 2- The distance are expressed in the same unit as wavelength and cell parameters (usually in angstroms) and the angles are always expressed in radians: 3- The listing of the different options is a function of the main optional parameter P6. In the following, the temperature parameters are exclusively Boveral or Biso (isotropic Debye-Waller factors) except for TLS option with abs(P6)=5. abs(P6)= 1.0x : Normal rigid body option. If x=1 the fractional coordinates of the center of mass is output. That supposes that every atom of the Molecule has been given explicitely in the asymmetric unit. 2.xx : Satellite group (RBSG) option (Int(abs(P6))=2). The integer value xx=100*(int(abs(P6))-2) gives the absolute number of the reference group (as they appear following the writing order, whatever the number of phase). --> The parameter P15 is assigned to the rotation CHIc of RBSG as for RBG. --> The parameter P16 of the first atom of the RBSG is defined as follows: P16 = N1.N2 with N1= Int(P16) : Nr of the first reference atom of the reference RBG. N2=100*(P16-Int(P16)): Nr of the second reference atom of the reference RBG. N1: Centre of the RBSG N2: Optional (if N2 is given then N1-N2 defines the internal z-axis) Example: 3.02 1st reference atom = Nr 3 2nd reference atom = Nr 2 Atom Nr 3 is the centre of the RBSG and the z axis of the RBSG is oriented in the direction: atom Nr 2-> atom Nr 3. Of course, the (xo, yo, zo) fractional coordinates of the RBSG are not needed. The program calculates automatically the corresponding values.

a) If the second atom is defined (N2#0), the spherical angles of the RBSG are calculated from these two reference atoms. b) If the second atom is not defined (N2=0), the centre of the main RBG is taken as the second reference atom of the RBSG. c) If N2=N1, the spherical orientation angles of the RBSG are those of the main reference RBG. Apart from these constraints, a RBSG is treated as a normal RBG with as many atoms as desired. The refinement of the centre of the RBSG as well as the orientation angles THETAc and PHIc are performed in the main RBG. 3.0x : The spherical coordinates (parameters P9, P10 and P11) are generated from the reduced coordinates (P1, P2 and P3) at the first cycle. Then the parameter P6 is automatically set to 1 and the sign of the option is kept i.e. the next option is RBG. For selecting the internal orthogonal system the user has to give in the parameters (P12,P13,P14) the origin of the intenal orthogonal frame and in (P9,P10,P11) the coordinates of an atom (that may be fictitious) for defining the plane "xz". The (P1,P2,P3) coordinates of the first atom defines the z-axis of the internal frame, which is in the direction V3=(P1-P12,P2-P13,P3-P14) in the conventional crystal frame (fractional components). the "y" axis in perpendicular to the plane "xz" defined by the vectors: V3=(P1,P2,P3)-(P12,P13,P14) -> z-axis V1=(P9,P10,P11)-(P12,P13,P14) -> within the xz plane y-axis is in the direction V3 x V1 This option is very usefull as it facilitates the use of standard input file with JBT#4, the conversion of published structures into spherical internal coordinate systems, etc... If x=1 the fractional coordinates of the center of mass is output as in option 1. 4.xx : Option to generate an idealized molecule like aliphatic chains, planar or helicoidal molecules..., where xx is the number of atomic planes along the RBG Zc axis which is the reference axis. The spherical coordinates are calculated from three parameters P9,P10 and P11. These parameters are lost after processing. They are defined as: P9 - distance of one atom to the center of the axis P10 - angle between two atoms lying in consecutive planes and the corresponding mid-point lying along the Zc-axis. P11 - order of the generating axis (integer) Examples: 1-For an aliphatic chain Cn along Zc in the

(Xc, Zc) plane (conformation all trans). distance = 1.54/2. (angstroms) angle = 114.*pi/180. (radians) axis order = 2 2-For benzene lying in the plane (Xc, Yc) distance # 1.40*2. (angstroms) angle = 0.*pi/180. (radians) axis order = 6 The parameter P15 is assigned to CHIz as in the RBG option and the parameter abs(P6) is set to 1 in the 1st cycle, the sign of the option is kept. Therefore, the generated group is treated as a normal RBG in the next cycles. Remark : Further improvements are schedduled for this option which is not very trustly at the present time. 5.0x : TLS option for RBG including satellite groups, if any. This TLS version is based on the formalism of V.Schomaker and K.N.Trueblood, Acta Cryst.(1968), B24, p-63. In this case, the refinement is performed with the so called one step process i.e. atomic positions and temperature factors TLS are refined together. The origin of the main RBG MUST be the centre of mass of the entire group concerned !!! The subroutine calculates the centre of mass of the molecule if x=1. This option MUST be set when the centre of the group does not coincide with the centre of mass. In this case, make sure that all atoms constituting the RBG are defined in the input file since the RBG option does not generate atomic positions. Such a situation can occur in non-centrosymmetric space group, as in urea SPG: P -4 21 M , since the molecule is located on a C2v site (mm.). Here, one must enter the entire molecule in the asymmetric unit. The elements of the T(6), L(6) and S(9) matrices are read in "further parameters IFURT" (see Fullprof guide, line 11.12) following the usual order ie. 11 22 33 12 13 and 23 for T and L. As the S matrix is not symmetric, 9 elements can be required in the general case. The six first elements are as defined previously, and the three last ones are respectively S21, S31 and S32. Similarly to the atomic temperature parameters betaij, the components of the TLS matrices have symmetry constraints, which are imposed by the symmetry of the crystallographic site (not the molecular symmetry !). These symmetry relations can be found in the paper by Schomaker and Trueblood (1968). The T elements are expressed in angstroms^2, the L components in radians^2 and the S ones in angstroms* radians. For convenience, the output to the screen

and to unit 7 (x.out) of the TLS components are expressed in the following units: angstroms^2 for T, in degrees^2 for L and degrees*angstroms for S. Remark : When using this option, one refines the observed structure factors by assuming that the atomic temperature factors are constrained ab initio to satisfy the rigid-body hypothesis. It is wellknown that the RBG/TLS can greatly reduce the number of atomic and thermal parameters, especially when the molecule is located at a site of high symmetry, but the user should be familiar enough with the TLS hypothesis not to perform inconsistent refinements. 4- It is advised to use specific rigid body refinement codes as damping parameters,especially for the rotational parameter which need to be currently three time higher than the other ones. It is highly recommended to use the options with screen output (P6 with a negative value) to control each RBG refinement steps since good constrained refinements are not so easy to perform. 5- Concerning the refinements: when a main RBG has satellites groups, the derivatives of the general RBG parameters, i.e. the 3 orientation angles THETAc, PHIc and CHIc and the coordinates (xo, yo, zo) of the centre of the group, contain the contributions of the satellites groups. In contrast, the derivatives of the internal spheric parameters do not account for these contributions. 6- When this subroutine is used, three additionnal files of fixed names are created at the end: x#.m, x#.bs and x#.ortep (where # stands for the number of the phase), containing atomic parameters in a format suitable to be used with well-known software packages (Molview, Balls &Sticks and ORTEP). Note that files for Molview are directely executable. The output for the two other software are incomplete in the sense that they contain either cartesian coordinates (Balls& Sticks) or reduced coordinates + atomic temperature parameters (ORTEP) of atoms located in the asymmetric unit. MAGNETIC STRUCTURE REFINEMENTS In the present version of FullProf the default JBT=5 option corresponds to the calculation of the magnetic structure factor for a conical magnetic structure. The magnetic intensities are calculated following standard formula as given, for instance, in the paper by J.M.Hasting & L.M.Corliss published is PhysRev 126(2),556 (1962). No symmetry operation can be introduced: all the magnetic atoms within a primitive unit cell must be given (constraints have to be introduced through the codes of the parameters). The subroutine MAGMOD, in this case, does not take into account symmetry for calculations (see the modules FDUM's.FOR for details). However ISYM must be set to 1, the value of the four parameters MULT, ICENT, NLAUE, NMAGR should be: 0 1 1 0. The atom parameters correspond to the following variables: -> P1,P2,P3,P4 and P5 correspond to x,y,z,B,occ of the atom -> P6 is the magnetic moment of the atom (in Bohr magnetons) -> P7 is the half-angle cone of the atom (degrees)

-> P8 is the magnetic -> P9 of the first -> P10 of the first (Phi,Theta) are the

phase of the atom (in fractions of 2 pi) atom correspond to Phi (degrees) atom correspond to Theta (degrees) spherical angles(degrees) of the cone axis

The orthonormal system with respect to which are defined the spherical angles verifies: X axis coincides with the crystallographic A Y axis belongs to the plane A,B Z axis is perpendicular to the plane A,B The particular implementation of spherical components in magnetic structure refinements is that Z axis must coincide with C. That works in all crystallographic systems except for triclinic. The monoclinic setting must be changed to -> 1 1 2/m It is recommended to generate the reflections using P-1 as Space group symbol and use for conical structures the value IRF=0 (Generates satellites+ fundamentals). For a pure helix IRF=-1 (only satellites are generated). The number of propagations vectors must be set to -1 (to generate + and - satellites). For calculating the magnetic moments in different cells the following formula should be used: The magnetic intensity is given by the following formula: Fm**2 =(p.sinw)**2.Sum(j){m(j) f(H,j) cos(bj) exp(2pii.H rj)} Fm**2 =p**2.(1+cosw**2)/4.Sum(j){m(j) f(H+-k,j) sin(bj) exp(2pii.Hrj-+fj)} The correspondence with the parameters P1,..P10 is the following: rj = ( P1, P2, P3)j f(H+-k,j) = P5j. FormFactor(H+-k,j) exp (-P4j(sintheta/lamda)**2) m(j) = P6j bj = P7j fj = P8j The unitary vector defining the axis of the cone is given by: n = (cos P9 sin P10, sin P9 sin P10, cos P10) cos w = n . Q/mod(Q) , where Q is the scattering vector ------------------------------------------------------------------5.- SINCLE CRYSTAL AND INTEGRATED INTENSITY REFINEMENTS 5.1: General comments. The version 3.0.0 (or higher) of FullProf permits the refinement of integrated intensity data. Single crystal and/or powder integrated intensities can be included (or used alone) as observations for refining a structural model. This possibility increases the capabilities of FullProf for handling other kind of data, but it cannot compete (or substitute) more elaborated programs specialized in single crystal work.

The structure factor calculation is exactly the same as in powder diffraction, so all the available features can be directly used with integrated intensity data. The R-factors in single crystal work are calculated according to the following formulae: Optimized function: M = Sum{n}[w(n)(F2obs(n)-Sum{k}[F2cal(k)])^2] The index "n" runs over the observations (1,..Nobs) The index "k" runs over the reflections contributing to the observation "n" F2 is the square of the structure factor (intensity corrected for Lp-factor) RF2 -factor=100*Sum{n}[|F2obs(n)-Sum{k}[F2cal(k)]|]/Sum{n}[F2obs(n)] RF2w-factor=100*Sqrt(M/Sum{n}[w(n)F2obs(n)^2]) RF -factor=100*Sum{n}[|Fobs(n)-sqrt(Sum{k}[F2cal(k)])|]/Sum{n}[Fobs(n)] Fobs=Sqrt(F2obs) Chi2(Intens)=M/(Nobs-Npar) Notice that RF2w for w(n)=1.0 is not the same value as RF2. For integrated intensity powder data, F2obs is in fact Sum{k}[jLpF2(k)]. When the option IWGT=2 (unit weights) is used. The inverse of the normal matrix is not multiplied by Chi2 as is usual in weighted refinements. 5.2: The extinction correction. At present only the "empirical correction" (a single parameter), as used in SHELX, has been implemented. I hope to have time for putting the anisotropic Becker-Coppens extinction coefficients into the program. 5.3: Mixed refinements. This option is still in an exploring stage. For the moment only a single powder diffraction profile can be used with different sets of integrated intensity data that are related to each phase. A global powder diffraction pattern can be given as primary information, and some "phases" can be given in addition with their own integrated intensity data. For mixing X-ray and neutron refinement, the best thing to do is to provide a "powder" neutron diffraction pattern as the global data information and a set of integrated intensities for the crystallographic phase to be refined. A present it is necessary to "repeat" the same phase, as in that modelling the neutron pattern, using the appropriate constraints in the refined parameters (using the same code for physically identical parameters). ------------------------------------------------------------------6.- REFERENCES

1.- H.M. Rietveld, Acta Cryst. 22, 151-1152 (1967) 2.- H.M. Rietveld, J. Applied Cryst. 2, 65-71 (1969) 3.- A.W. Hewat, Harwell Report No. 73/239, ILL Report No. 74/H62S

4.- G. Malmros & J.O. Thomas, J. Applied Cryst. 10, 7-11 (1977) 5.- C.P. Khattak & D.E. Cox, J. Applied Cryst. 10, 405-411 (1977) 6.- D.B. Wiles & R.A. Young, J. Applied Cryst. 14, 149-151 (1981)/ 15, 430-438 (1982) 7.- G.S. Pawley, J. Applied Cryst. 14, 357-361 (1981) 8.- Halpern and Johnson 9.- E. Prince, J. Appl. Cryst. 16, 508 (1983) 10.- W.A. Dollase, J. Applied Cryst. 19, 267- (1986) 11.- K.D. Rouse, M.J. Cooper, E.J. York & A. Chakera, Acta Cryst A26, 682 (1970) 12.- A.W. Hewat, Acta Cryst. A35, 248 (1979) 13.- A.J.C. Wilson, Mathematical theory of X-ray powder diffraction Philips Technical Library, Eindhoven (1963) 14.- J.F. Berar & *. Lelann, J. Appl. Cryst. 24, 1-5 (1991) 15.- A. Antoniadis, J. Berruyer & A. Filhol, Acta Cryst. A46, 692-711 (1990) 16.- R.J. Hill & H.D. Flack, J. Applied Cryst. 20, 356-361 (1987) 17.- J. Rodriguez-Carvajal, M.T. Fernandez-Diaz and J.L. Martinez, Journal of Physics: Condensed Matter 3, 3215-3234 (1991) 18.- J. Rodriguez-Carvajal, Physica B 192, 55-69 (1993) 19.- Numerical Recipes, by W.H. Press, B.P. Flanery, S.A. Teukolsky and W.T. Vetterling (Fortran version), Cambridge University Press, 1990. 20.- W. Pitschke, N. Mattern and H. Hermann, Powder Diffraction 8(4), 223-228 (1993). (In this paper the authors refine simultaneously all occupation numbers. This practice should be avoided. If the calculation were exact that gives rise to a singular matrix when the scale factor is also refined) See also W. Pitschke, H. Hermann and N. Mattern, Powder Diffraction 8(2), 74-83 (1993).

7.- DIMENSIONS OF ARRAYS OF THE STANDARD VERSION

Dimensions can be adapted to the available memory by changing the values in PARAMETER statements IDSZ: IRS: NATS: Maximum number of points in the diffraction pattern Maximum number of reflections Maximum number of atoms in asymmetric unit (all phases included)

MSZ: NOV: NMAGM: NATM: NEQV: NPR: INPR:

Maximum Maximum Maximum Maximum Maximum

refinable parameters overlapping reflections rotation-matrices sets for magnetic structure magnetic atoms (NATM+nonmag-atoms <=NATS) user-supplied symmetry operators/ /propagation vectors Maximum number of points for defining a numerical profile Maximum number of different numerical profiles (constraint: IDSZ >= MSZ*MSZ) (constraint: IDSZ >= 4*IRS )

number number number number number

of of of of of

IEXCL: IBACP: NPHT: MPAR: NGL : NCONST:

Maximum Maximum Maximum Maximum Maximum Maximum

number number number number number number

of of of of of of

excluded regions background points phases (<=8) non atomic parameters for each phase global parameters slack constraints per phase

Integer*4 + IDSZ,IRS,NATS,MPAR,MSZ,NOV,IEXCL,IBACP,NPHT,NMAGM,NATM,NEQV,NGL, + i_Dat, i_Sym, i_Rpa, i_Pcr, i_Out, i_Sum, i_hkl, + i_Fou, i_Ghkl, i_Sav, i_Shkl, i_Prf, i_Bac, i_Res, i_Sim, + i_Sch, i_New, i_Shx, i_Atm, i_Sub Usual Parameters for PC, 4Mbytes core memory PARAMETER PARAMETER PARAMETER PARAMETER (IDSZ=8600,IRS=2600,NATS=100,MPAR=90,MSZ=90,NOV=256) (IEXCL=30,IBACP=100,NPHT=6,NMAGM=6,NATM=46,NEQV=24) (NAT_P=20,NGL=40,NCONST=16,N_spe=8,N_form=60,NPR=150) (INPR=4)

Usual Parameters for PC, 8Mbytes core memory PARAMETER PARAMETER PARAMETER PARAMETER (IDSZ=24000,IRS=5300,NATS=230,MPAR=90,MSZ=150,NOV=600) (IEXCL=30,IBACP=100,NPHT=8,NMAGM=8,NATM=56,NEQV=48) (NAT_P=25,NGL=40,NCONST=16,N_spe=8,N_form=60,NPR=150) (INPR=4)

Logical units used in the program Parameter + (i_Pcr=1,i_Rpa=2,i_Sym=3,i_Dat=4,i_Out=7,i_Sum=8) Parameter + (i_Fou=9,i_Ghkl=10,i_Sav=11,i_Bac=12,i_Res=13,i_Sim=14) Parameter + (i_hkl=15,i_Shkl=16,i_Prf=17,i_Sch=18,i_New=19,i_Shx=21) Parameter + (i_Atm=22,i_Sub=23,i_shp=24) LLB, jan. 98

You might also like