FullProf. Rietveld, Profile Matching & Integrated Intensities Refinement of X-ray and/or Neutron Data (powder and/or single-crystal).
SHORT REFERENCE GUIDE OF THE PROGRAM FFFFFF PPPPPP F l l P P f F l l P P f FFFF l l PPPPPP f F u u l l P r r r ooo fff F u u l l P r o o f F uuuuu lll lll P r ooo 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 -> Go to FullProf area ftp>get fpreadme -> Obtain the document ftp>bye -> Return to host From some VMS-VAX hosts: ftp>set def "pub/divers/fullp" -> Go to FullProf area ftp>get "fpreadme" -> Obtain the document ftp>ex -> 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) - 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 file - 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 file (unit 12). The format of this file is the following (as that of FILE.DAT for INSTRM=0): first line : 2theta(initial) step 2theta(final) following lines: 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 "observed" structure factor squared. The file could be used as input 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 factor. This file can be used as input for JBT=-3 and IRF=2 in order 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 X-ray case (Job) 1 Neutron case (constant wavelength, nuclear and magnetic) 2 pattern calculation (X-ray) 3 pattern calculation (Neutron, constant wavelength) -1 Neutron case (T.O.F., nuclear and magnetic) -3 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 Gaussian 1 Cauchy 2 Modified 1 Lorentzian 3 Modified 2 Lorentzian 4 Tripled pseudo-Voigt 5 pseudo-Voigt 6 Pearson VII 7 Thompson-Cox-Hastings pseudo-Voigt 8 Numerical profile given in CODFIL.SHP or in GLOBAL.SHP 9 T.O.F. Convolution pseudo-Voigt x Double Exponential 10 Not yet used 11 Split pseudo-Voigt function 12 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 approximation of C.J.Howard, J.Appl.Cryst.15 615-620 (1982) with the Simpson formula for five points) NBCKGD =0 Refine background with polynomial function (Nba) 1 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 file FILE.BAC as for NBCKGD=1. -3 Read 6 additional polynomial background coefficients NEXCRG = number of excluded regions (Nex) NSCAT = number of scattering sets (zero in most cases) (Nsc) 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. NORI = 0 preferred orientation function No 1 (Nor) 1 preferred orientation function No 2 (March) 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 Number of parameters to be constrained within given (Nre) limits. At the end of the file you must give a list of NRELL lines specifying the number and the limit of each parameter. ICRYG If not equal to zero, only integrated intensity data (Cry) 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 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 Polynomial function: cor = Sum{i=1,npcorr}{coeff(i)* T**(i-1))} ITYFUNC=2 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 = correlation matrix --> CODFIL.OUT (Mat) 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) 4 In some versions of FullProf a plot of the diffraction pattern is diplayed on the screen at each cycle of refinement. LST3 = merged reflection list --> CODFIL.OUT (usually 0) (Ls3) IPL1 = symmetry operators --> CODFIL.OUT (+ CODFIL.SYM if JSY=1) (Syo) IPL2 = output data for plot --> CODFIL.PRF (Prf) = 1 Format suitable for PLOTPOW, BENSTRAP,PLOTR, etc.. = 2 " " " IGOR (MacIntosh software) = 3 " " " KaleidaGraph (MacIntosh software) & PLOTR (Pc software) = 4 " " " Picsure, Xvgr(Sun-Unix Software) 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 prepares CODFIL.SYM (if 1, IPL1 must be set to 1) (Sym) JLKH prepares CODFIL.HKL (Hkl) = 1 as explained above = 2 Output for SIRPOW.92 =-2 Output for EXPO = 3 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 Wdt Cthm muR AsyLim Rpolarz ---------------------------- OR (for T.O.F. data) BKPOS, WDT, IABSCOR (2 real + 1 integer) --------------Comment line : ! BKPOS Wdt Iabscor ---------------------------- LAMDA1 = wavelength l1 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 = origin of polynomial for background (in deg. 2theta or usecs) WDT = 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.) CTHM = coefficient for monochromator polarization correction (see Mathematical information) TMR = absorption correction coefficient mR, used only for (muR) 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 Thmin Step Thmax PSD Sent0 ---------------------------- MCYCLE = number of cycles of refinement EPS = forced termination when shifts less than EPS*e.s.d. RELAX = the four relaxation factors for shifts : 1 -> 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 parameters. 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 = incident beam angle at sample surface in degrees SENT0 = 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 = zero point for T (in degrees) : Ttrue = Texp.-ZER Note that the shift convention is opposite to that used in Pawley's program. FLGZER = codeword for zeroshift (codewords are described in the Mathematical section below). SYCOS = systematic 2Theta shift with cosTheta dependence sample displacement (Theta - 2Theta diffractometers) FLCOS = codeword for SYCOS SYSIN = systematic 2Theta shift with sin2Theta dependence sample transparency coefficient FLSIN = codeword for SYSIN LAMBDA = Wavelength to be refined (only 1-wavelength can be refined) FLAMBDA = Codeword for LAMBDA. Cell parameters should be fixed if wavelength is to be refined. IGLMORE 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 using 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 ATZ Nvk Npr More ---------------------------- N = 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) NDIST = number of distance constraints NMAGC = number of magnetic moment constraints PREF(1,2,3)= preferred orientation direction (in reciprocal space) 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. IRF = 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 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 k l m (IRF=1) (free format) h k l m Coeff (IRF=1+JSOL=1) ( " ) h k l m Intensity (or F) (IRF=2) ( " ) h k l m Freal Fimag (IRF=3) ( " ) 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 k l nv m (IRF=1) (free format) h k l nv m Coeff (IRF=1+JSOL=1) ( " ) h k l nv m Intensity (or F) (IRF=2) ( " ) h k l nv m Freal Fimag (IRF=3) ( " ) 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 Isolated peak ............................. 5 3 1 -1.0 0.0 1 Cluster of peaks: four 3 4 2 -1.0 0.0 1 independent reflections 4 4 1 -1.0 0.0 1 contribute to 5 0 3 832.1 9.4 1 <- this observation ISYM=0 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 l m F2(calc) F2(obs) k-vectors -> h k l nv m F2(calc) F2(obs) hr kr lr with obvious meaning. ISTR =0 If strain or/and size parameters are used, they are those corresponding to selected models (see below) =1 The generalized formulation of strains parameters will be used for this phase. See below =2 The generalized formulation of size parameters will be used for this phase. See below =3 The generalized formulation of strain+size parameters will be used for this phase. See below IFURT 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 only 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 = 1 Additional hkl-dependent asymmetry and shifts parameters are read. JMOM - unused at present JTER - unused at present BRIND - Brindley coefficient (see line 11-2) RM1 - 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. RM2 - 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. RM3 - 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. JTYP - 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. 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 T2 S31 S32 S33 T3 !M11 M12 M13 M21 M22 M23 M31 M32 M33 Ph ---------------------------- . 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 X,Y,Z MSYM U,V,W, 0.0 ! SYMM X+1/2,-y, Z MSYM -U, V,-W, 0.0 ! SYMM -x,-y,-Z MSYM 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 p3 p4 p5 p6 p7 p8 ! p9 p10 p11 p12 p13 p14 p15 p16 ---------------------------- ---------------------------------------------------------------------- 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 Occ Rx Ry Rz ! Ix Iy Iz beta11 beta22 beta33 MagPh ---------------------------- 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 Occ Mom beta Phase ! Phi & Theta of Cone-axis + unused params ---------------------------- ---------------------------------------------------------------------- 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. IMAGR= 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 IK= Number of the propagation vector to which the atom 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 = isotropic displacement (temperature) parameter in angstroms**2 N = occupation number i.e. chemical occupancy x site multiplicity (can be normalized to the multiplicity of the general position of the group). IOPIN,IOPFIN= 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!) 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 = anisotropic displacement (temperature) parameters (betas) OR (For N_typ=4) (AND if N_typ=5, but not available) --------------Comment line : ! Form-factor refinable parameters ---------------------------- Line 11-4-3: f1 f2 f3 f4 f5 f6 f7 (7 reals) Line 11-4-4: Cf1 Cf2 Cf3 Cf4 Cf5 Cf6 Cf7 Line 11-4-5: f8 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, B, N N_type (a) (2a4,2int,5real,1int) Line 11-4-2: CX, CY, CZ, CB, CN (b) (5 reals) Line 11-4-3: Mx My Mz Mxi Myi Mzi MPhas (c) (7 reals) Line 11-4-4: CMx CMy CMz CMxi CMyi CMzi CMPhas (d) (7 reals) Line 11-4-5: b11 b22 b33 b12 b13 b23 (e) (6 reals) Line 11-4-6: Cb11 Cb22 Cb33 Cb12 Cb13 Cb23 (f) (6 reals) Line 11-4-7: f1 f2 f3 f4 f5 f6 f7 (g) (7 reals) Line 11-4-8: Cf1 Cf2 Cf3 Cf4 Cf5 Cf6 Cf7 (h) (7 reals) Line 11-4-9: f8 f9 f10 f11 f12 f13 f14 (i) (7 reals) Line 11-4-10: Cf8 Cf9 Cf10 Cf11 Cf12 Cf13 Cf14 (j) (7 reals) If N_type = 0 Only lines (a) and (b) need to be given If N_type = 1 give the lines (a), (b), (c) and (d) If N_type = 2 give the lines (a), (b), (e) and (f) If N_type = 3 give the lines (a) -> (f) if N_type = 4 give the lines (a), (b) and (g)->(j) (special form-factor) --------------Comment lines : !Atom Typ Mag Vek X Y Z Biso Occ N_type !Line below:Codes ! Rx Ry Rz Ix Iy Iz MagPh !Line below:Codes or... ! beta11 beta22 beta33 beta12 beta13 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 = scale factor GAM1 = profile shape parameter , e.g. : eta0 for NPROF=4,5 but not for NPROF=7 (see line 11-6-1) m0 for NPROF=6 Bov = overall isotropic displacement (temperature) factor in angstroms**2 STR1,STR2,STR3 = strain parameters, defined through the subroutine STRAIN (see additional information) If ISTR=1 set these values to 0.0 IstrainModel = 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) --------------Comment line : ! Scale Factors ! Sc1 Sc2 Sc3 Sc4 Sc5 Sc6 ---------------------------- Line 11-5-2: CSc1, CSc2, CSc3, CSc4, CSc5, CSc6 (6 reals) Sci = scale factor for domain (i) CSci = code of the scale factor for domain (i) 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 LorSiz Size-Model or for NPRO=11 (split pseudo-Voigt) ! Ul Vl Wl Xl Y GauSiz LorSiz Size-Model ---------------------------- (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 =====> 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 =====> 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(Theta) 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) --------------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 Ext7 ---------------------------- Line 11-6-2: CExt1, CExt2, CExt3, CExt4, CExt5, CExt6, CExt7 (7 reals) Exti = Extinction parameter (i) CExti = Code of the extinction parameter (i) 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(2Theta) 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 Xt Yt Z1 Z0 Size-Model -------------------------------------------------------------------------------- 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 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 ---------------------------------------------------------------------------- Line 11-8-3: Abs1, CAbs1, Abs2, CAbs2 (4 reals) --------------Comment line : !Absorption correction parameters ---------------------------------- Abs1, Abs2 : Absorption correction parameters CAbs1, CAbs2 : 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 VALUEPAR CODEPAR --------------Comment line : ! Further parameters: ---------------------------- ---------------------------------------------------------------------------- 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 Symmetry symbol Rotation matrix ( 1) 1 --> ( x, y, z) ( 2) 2 ( 0, 0, z) --> (-x,-y, z) ( 3) 2 ( 0, y, 0) --> (-x, y,-z) ( 4) 2 ( x, 0, 0) --> ( x,-y,-z) ( 5) 3+ ( x, x, x) --> ( z, x, y) ( 6) 3+ (-x, x,-x) --> ( z,-x,-y) ( 7) 3+ ( x,-x,-x) --> (-z,-x, y) ( 8) 3+ (-x,-x, x) --> (-z, x,-y) ( 9) 3- ( x, x, x) --> ( y, z, x) (10) 3- ( x,-x,-x) --> (-y, z,-x) (11) 3- (-x,-x, x) --> ( y,-z,-x) (12) 3- (-x, x,-x) --> (-y,-z, x) (13) 2 ( x, x, 0) --> ( y, x,-z) (14) 2 ( x,-x, 0) --> (-y,-x,-z) (15) 4- ( 0, 0, z) --> ( y,-x, z) (16) 4+ ( 0, 0, z) --> (-y, x, z) (17) 4- ( x, 0, 0) --> ( x, z,-y) (18) 2 ( 0, y, y) --> (-x, z, y) (19) 2 ( 0, y,-y) --> (-x,-z,-y) (20) 4+ ( x, 0, 0) --> ( x,-z, y) (21) 4+ ( 0, y, 0) --> ( z, y,-x) (22) 2 ( x, 0, x) --> ( z,-y, x) (23) 4- ( 0, y, 0) --> (-z, y, x) (24) 2 (-x, 0, x) --> (-z,-y,-x) Hexagonal frames ITnum Symmetry symbol Rotation matrix (25) 1 --> ( x , y, z) (26) 3+ ( 0, 0, z) --> ( -y, x-y, z) (27) 3- ( 0, 0, z) --> (-x+y,-x , z) (28) 2 ( 0, 0, z) --> (-x , -y, z) (29) 6- ( 0, 0, z) --> ( y,-x+y, z) (30) 6+ ( 0, 0, z) --> ( x-y, x , z) (31) 2 ( x, x, 0) --> ( y, x ,-z) (32) 2 ( x, 0, 0) --> ( x-y, -y,-z) (33) 2 ( 0, y, 0) --> (-x ,-x+y,-z) (34) 2 ( x,-x, 0) --> ( -y,-x ,-z) (35) 2 ( x,2x, 0) --> (-x+y, y,-z) (36) 2 (2x, x, 0) --> ( 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 complex 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), scattering 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.
Fullprof Manual Html Español
The Best Free Web Hosting. We offer outstanding free web hosting with Unlimited Resources Supporting PHP, MySQL and Wordpress. Wink Web-Hosting has redefined what.
Document/File: fullprof-tutorial.pdf, filesize: n/a. MANUAL DE BIOINGENIERíA - MÁSTER SERGIO J NAVARRO HUDIEL; FICHA PRODUCTO - COMESTEROGROUP; PDF.
Tutorials and Examples Winplotr Graphics and Fullprof Rietveld (and Single Crystal). WinPlotr Documentation in HTML format. FullProf user's guide in HTML format.
Fullprof Manual Html5
Juan Rodríguez-Carvajal Institute Laue-Langevin Grenoble, FRANCE Introduction to the programs of the FullProf Suite.
FullProf user's guide