[a E] [b E] [c E] [al E] [be E] [ga E]

Lattice parameters in Angstroms and lattice angles in degrees.

[adps] [u11 E] [u22 E] [u33 E] [u12 E] [u13 E] [u23 E]

adps when used generates the unn atomic displacement parameters with considerations made for special positions. On termination of refinement the adps keyword is replaced with the unn parameters; see example AE1-ADPS.INP. Instead of using the adps keyword the unn parameters can be manually entered.

The unn matrix can be kept positive definite with the site dependent macro ADPs_Keep_PD; this can stabilize refinement.  The ADPs_Keep_PD macro can be used after the unn parameters are created.


Signals that the associated phase is amorphous when calculating degree_of_crystallinity.


Generates the un-normalized and normalized A and correlation matrices. If do_errors is defined then C_matrix_normalized is automatically generated and appeneded to the OUT file.

[append_cartesian] [append_fractional [in_str_format] ]

Appends site fractional coordinates in Cartesian or fractional coordinates respectively to the *.OUT file at the end of a refinement. For the case of append_fractional, the in_str_format keyword formats the output in INP format.

[append_bond_lengths [consider_lattice_parameters] ]

Appends bond lengths to the end of the *.OUT file at the end of refinement. A number corresponding to equivalent positions is appended to site names.

consider_lattice_parameters includes the effects of lattice parameter errors in the calculation of bond length and bond angle errors.

An example of bond lengths output is as follows:

Y1:0    O1:0    2.23143

        O2:0    2.23143    88.083

        O3:0    2.28045    109.799    99.928

The first line gives the distance between the sites Y1 and O2. The first number in the second line gives the distance between sites Y1 and O2. The third number of 88.083 gives the angle between the vectors Y1 to O1 and Y1 to O2. The first number on the third line contains the distance between sites Y1 and O3. The second number in the third line contains the angle between the vectors Y1 to O3 and Y1 to O2. The third number in line three contains the angle between the vectors Y1 to O3 and Y1 to O1. Thus bond lengths correspond to the first number in each line and bond angles start from the second number. The numbers after the site name and after the ‘:’ character corresponds to the site equivalent position as found in the *.SG space group files found in the SG directory


[A_matrix_memory_allowed_in_Mbytes !E]

[A_matrix_elements_tollerance !E]


See section 5.2 for a description of approximate_A.

A_matrix_memory_allowed_in_Mbytes limits the memory used by the A matrix to a maximum of A_matrix_memory_allowed_in_Mbytes. If the matrix requires less than A_matrix_memory_allowed_in_Mbytes then the full matrix is used otherwise the matrix is treated as a sparse matrix.

A_matrix_elements_tollerance removes elements in the A matrix with values less than A_matrix_elements_tollerance. The comparison is made againts normalized elements of A such that the diagonals have a values of 1.  The A matrix is made sparse when A_matrix_elements_tollerance is defined. Typical values of A_matrix_memory_allowed_in_Mbytes range from 0.0001 to 0.01.

A_matrix_memory_allowed_in_Mbytes and A_matrix_elements_tollerance can be used simultanuously.

A_matrix_report_on displays the percentage of non-zero elements in the A matrix.

 [atomic_interaction N E] |  [ai_anti_bump N]…

ai_sites_1 $sites_1 ai_sites_2 $sites_2


[ai_closest_N !E]

[ai_radius  !E]



Defines an atomic interaction with the name N between sites identified by $site_1 and $site_2. For atomic_interaction E is the site interaction equation that can be a function of R and Ri. R returns the distance in Å between two atoms; these distances are updated when dependent fractional atomic coordinates are modified. The name of the atomic_interaction N can be used in equations and in particular penalty equations.

For ai_anti_bump an anti-bump interaction equations is internally generated. For anti-bumping only the ai_anti_bump is faster than using atomic_interaction. The macro AI_Anti_Bump usesai_anti_bump.

no_self_interaction prevents any interactions between equivalent positions of a site. This is useful when a general position is used to describe a special position.

ai_closest_N: interactions between $sites_1 and $sites_2 are sorted by distance and only the first ai_closest_N number of interactions are considered.

ai_radius: only the interactions between $sites_1 and $sites_2 that are within the distance ai_radius are considered.

When ai_radius and ai_closest_N are both defined then interactions from both sets of corresponding interaction are considered.

ai_exclude_eq_0: only interactions that is not the first equivalent positions in $sites_2 are considered. For example, in the following:



ai_sites_1 Pb

ai_sites_1 “O1 O2”

the following interactions are considered:

Pb:0 and O1:n      (n ¹ 0)

Pb:0 and O2:n      (n ¹ 0)

where the number after the ‘:’ character corresponds to the equivalent positions of the sites.

ai_only_eq_0:only interactions between equivalent positions 0 are considered.


The atomic_interaction’ equation can be a function of the following functions:

AI_R(#ri): Returns the distance between the current site and the atom defined with Ri = #ri.

AI_R_CM: A function of no arguments that returns the geometric center of the current atom and the atoms defined in $sites_2.

AI_Flatten(#toll): A function that returns the sum of distances of the current atom and those defined in $sites_2 to an approximate plane of best fit. The plane of best fit is constructed such that the sum of the perpendicular distances to the current atom plus those defined in $sites_2 are a minimum

AI_Cos_Angle(#ri1, #ri2): Returns the Cos of the angle between the atom define as Ri=#ri1, the current atom and the atom defined as Ri=#ri2.

AI_Angle(#ri1, #ri2) : Similar to AI_Cos_Angle except that the value returned is the angle in degrees.

Examples BENZENE_AI1.INP, BENZENE_AI2.INP and BENZENE_AI3.INP demonstrates the use of the atomic_interation functions.

atomic_interaction’s can be used to apply geometric restraints. For example, an anti-bump interaction between symmetry related molecules can be formulated as follows:

atomic_interaction ai1 =








ai_sites_1 C*

ai_sites_1 C*

ai_radius 3

penalty = If(Cycle_Iter < 10, ai1, 0);

This example demonstrates anti-bumping between molecules for the first ten iterations of a refinement cycle.

[atom_out $file [append] ]…

Used for writing site dependent details to file. See the keyword out for a description of out_record. The Out_CIF_STR uses atom_out


filament_length E sample_length E receiving_slit_length E

[primary_soller_angle E]

[secondary_soller_angle E]

[axial_n_beta !E]

Defines the full axial divergence model using the method of Cheary & Coelho (1998b).

filament_length, sample_length and receiving_slit_length define the length of the tube filament, sample and receiving slit in the axial plane in mm.

primary_soller_angle: primary Soller slit angle in degrees.

secondary_soller_angle:  secondary Soller slit angle in degrees.

axial_n_beta: Define the number of rays emanating from a point X-ray source in the axial plane. Larger values for axial_n_beta increases both accuracy and calculation time.

The macro Full_Axial_Model simplifies the use of axial_conv.

Personal Tools