Jana2006 Cookbook
Jana2006 Cookbook
Authors
Michal Dušek, Václav Petříček, Lukáš Palatinus, Jan Rohlíček, Karla Fejfarová,
Václav Eigner, Monika Kučeráková, Morgane Poupon and Margarida Henriques
Diamond can be replaced by any other plotting software capable to read CIF in command line. For the
workshop we recommend Diamond because the connection with Jana2006 is very well tested.
Installation of Diamond
Install demo version of Diamond from the directory Setup\Diamond installation in the flash disk.
Copy the licence file “yourlicence.lic” to the same directory where Diamond.exe has been installed.
Usually it will be c:\Program files\Diamond 4
This time limited license has been kindly provided by Crystal Impact for this workshop.
Note: Diamond started from Jana2006 may behave differently than when started separately. For
instance, it doesn’t offer automatic wizards.
Installation of Vesta
Unpack the archive “Vesta_2.1.2.zip” from the directory Setup\Vesta to a directory of your choice
Installation of Jana2006
• Don’t use your “old” preinstalled
Execute janainst.msi from the directory versions of Jana2006
Setup\Jana2006 and follow instructions.
The installation creates shortcut to Jana2006 in the • Don’t use flash disk for calculation.
desktop and in Start → Programs Copy the example files to hard disk.
Start Jana2006, go to “Tools → About Jana2006” and
verify the version string. • Don’t install Jana2006 to desktop or
The version date should be close to the workshop date flash disk
Configuration of Jana2006
Start Jana2006
Go to “Tools → Preferences”, adjust size of the window (for instance 60% of the screen height) and font
height (about 11 points)
Go to “Tools → Programs”; focus with mouse the textbox “Graphic viewer”; press “Browse” and localize
Diamond.exe (usually in Program Files x86)
Focus "3d visualization of maps"; press “Browse” and localize VESTA.exe
(VESTA is needed for examples 9.x, 10.x, 12.x and 13.x)
For using a different text editor than Notepad define its pathname in the textbox “Editor name”.
Launching PETS:
PETS must be called with the input file as a command-line argument. There are several ways of launching
PETS:
Method 1 – the most convenient:
Locate the file .pts (input to PETS) on your system
Right-click to open context menu
“Open → Select program from a list → Browse…”
Locate the file pets.exe on your system; Open
Make sure the checkbox next to “Always use the selected program to open this kind of file” is checked
OK
PETS starts. Next time you want to run PETS, it is sufficient to double click the icon of the pts file you
want to run. Unfortunately, this method fails on some systems, notably if the extension pts is already
associated with another software.
Method 2:
Place the shortcut to PETS on Desktop:
Locate the file pets.exe on your system
Select it and copy it to clipboard (Ctrl+C)
Place cursor on the place on the desktop where you want to place the icon of PETS
Right-click → Paste shortcut
Now you can run PETS by drag-and-dropping the file .pts on the icon on the desktop. Unfortunately, even
this method does not work on all systems.
Method 3 – works always:
Open DOS terminal (Start → All programs → Accessories → Command Prompt)
Move to the directory with the file .pts you want to use (use command cd)
Run PETS by typing “pets filename.pts”, where filename.pts is the file you want to use.
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction → CCD”; NEXT
Select “Input from “sum” file”; OK
Leave all settings unchanged; NEXT
Leave all settings unchanged; NEXT
The program reads 24240 reflections from hkl file
For absorption correction select “None or done before importing”; NEXT
FINISH
3. Data Repository
Each line in data repository window corresponds to one data set converted to Jana format.
In our case, we use only one data set. Data repository is saved in file Zn.m95.
OK; YES to accept the data set
4. Symmetry Wizard
[On the screen: “Tolerances for crystal system recognition”.]
Symmetry Wizard can be started separately by “File → Reflection file → Make space group
test.
Leave all settings default; NEXT
[On the screen: Select Laue symmetry]
Select Orthorhombic Laue symmetry; NEXT
Select primitive unit cell; NEXT
[On the screen: Select space group]
The window shows possible space groups.
“Obs/all” column shows for each space group number of reflections, which should be
9. Refinement
[On the screen: basic window of Jana]
Right-click the icon of Refine.
[On the screen: refinement options]
Define 100 of cycles; OK
Choose “YES+START”
Refinement converges with R value about 6%
Start “Edit atoms”
[On the screen: list of atoms]
Press “Select all” and “Action → Edit/define”
Choose “harmonic” for ADP parameters. This defines anisotropic temperature parameters
for all atoms. OK, OK, Yes to save changes
Double-click the icon Refine.
Refinement converges with R value about 4 % for 73 refined parameters.
Standard refinement calculates with F2 while Jana2006 is still using F by default. In order to
get a standard result, we will re-refine the structure on F2:
Right-click the icon of Refine.
[On the screen: refinement options]
Check “Refinements on R(obs)**2”; OK; YES to start the program
Refinement converges with R factor 2.19%, goodness of fit 1.61, 74 parameters.
For good data, refinement on F and F2 usually gives very similar results.
Jana2006 does not refine the weighting scheme. The goodness of fit therefore corresponds
to the experimental sigma of reflections corrected by the instability factor. For the moment,
we are using the default instability factor 0.01. In the next step, we will use the instability
factor calculated during merging of reflections (see above):
PbSO4
Powder data measured with laboratory diffractometer
Input files: PbSO4.mac (powder profile data)
PbSO4.txt (additional information)
2. Import Wizard
Select “Various CW formats”, NEXT
CW stands for constant wavelength
Select “MAC format” and “Another/unknown method”; NEXT
File name is automatically set to”PbSO4.mac”.
Fill cell parameters 8.48 5.4 6.96 90 90 90
For wavelength type first select Kalpha1/Kalpha2 doublet and then select X-ray tube Cu
Select “Parallel setting”, “Glancing angle” 13.28815, “Perfectness 0.5”; NEXT; FINISH
Accept the data in Data repository
Select "Yes, I would like to continue with the wizard"; OK
At this point Jana2006 offers wizard for determining powder profile parameters, space
group and structure. We shall use the wizard but all its tools could be also started
separately.
The observed peaks are systematically broader which may indicate that Lorentzian
broadening is to be combined with the Gaussian profiles.
“Edit profile parameters”
[On the screen: Powder options]
In page “Profile” select Pseudo-Voigt and activate refinement of LY
“Run Refine”
After10 refinement cycles Rp~8%
“Show powder profile”
“Edit profile parameters”
[On the screen: Powder options]
In page “Profile” activate refinement of LX.
“Run Refine”
After 10 refinement cycles Rp~6.2%
“Show powder profile”
5. Structure solution
[On the screen: Structure solution]
This part of wizard can be also started separately by “Run → Solution”
Type formula Pb S O4; Formula units 4
Select “Use Superflip”
Check “Repeat Superflip until the convergence is detected”
Set “Maxcycles” to 1000
Press “Accept last solution”. This action will also close the structure solution wizard.
At the beginning of “Changes overview” there will be probably list of spurious atoms
automatically excluded from the refinement as their ADP is larger than the allowed limit. In
this example this is O4. The excluded atoms are still present in M40 file with zero
occupation. In order to delete them:
Run “Edit atoms”, press “Select rejected”
Press the button Action and delete them; OK
Start EditM50, select page "Composition" and press "Formula from M40". It should display
"Pb S O4".
Powder data measured with laboratory diffractometer experimental setup with Johannson-
type monochromator and transmission mode in asymmetric Guinier type arrangement.
2. Import procedure
Select “Structure: from CIF”; NEXT
[On the screen: Select input CIF file]
Select “Y2O3-icsd.cif”; OK
[On the screen: “CIF file does not contain any reflection block”]
Answer “Yes” to the question: “Do you want to import data from file?”
[On the screen: Specify type of the file to be imported]
Select “Powder data: various CW formats” (CW = constant wavelength); NEXT
[On the screen: Powder data from:]
Select “Riet7 formats”; NEXT
[On the screen: Complete/correct experimental parameters]
For wavelength type use “X-ray tube” Cu, deselected Kalpha1/Kalpha2 doublet and
wavelength 1.54051
Select “Parallel setting”, “Glancing angle” 13.65232 (CuKα, Ge[111]) “Perfectness 1”; NEXT;
FINISH
Accept the data in Data repository
The observed peaks are systematically broader which can indicate that Lorentzian
broadening is to be combined with the Gaussian profiles.
Start “Edit profile”
[On the screen: Powder options]
In page “Profile” select Pseudo-Voigt and activate refinement of LY; OK
Yes to rewrite changes
Start refinement (icon Refine)
Refinement converges in 10 cycles to Rp=7.1%
Start “Edit profile”
[On the screen: Powder options]
In page “Profile” activate refinement of GU; OK
Yes to rewrite changes
Start refinement (icon Refine)
Refinement converges in 10 cycles to Rp~7.0%
Start “Profile viewer”
[On the screen: Jana 2006 profile viewer]
Press “X exactly” and select the 2θ interval from 19 to 35.
5. Asymmetry by Berar-Baldinozzi
[On the screen: Basic Jana window]
Start “File → Structure → History”
[On the screen: History]
Select the original structure Y2O3; OK
[On the screen: Basic Jana window]
Start “File → Structure → Save As”
[On the screen: File manager for saving structure]
In File manager focus the right-bottom textbox and define a name for testing of the Berar-
Baldinozzi model i.e. “Y2O3-BB”; OK
Yes to continue with the new structure
[On the screen: Basic Jana window]
Start “Edit profile”
[On the screen: Powder options]
In page “Asymmetry” select “Berar-Baldinozzi”, use the default value and activate
refinement of the asym1, asym2, asym3 and asym4 parameters; OK;
Yes to rewrite changed files
[On the screen: Basic Jana window]
Start "Refine commands" (right-click on the icon "Refine")
In page "Basic" set 100 refinement cycles and Damping factor 0.3
OK; YES+START to save commands and start refinement.
Refinement converges after many cycles to Rp~5.8%
Start “Profile viewer”
[On the screen: Jana 2006 profile viewer]
Press “X exactly” and select the 2θ interval from 19 to 35.
30 Jana2006 Cookbook, version October 2017
The fit is even better than for the previous method but the number of asymmetry
parameters is now four.
The parameter HmS/L oscillates , within three times the standard uncertainty, about zero
which induces instability in the refinement. For this reason, we can change it to zero and
keep it fixed.
Start “Edit profile”
In page “Asymmetry” set HmS/L to zero and unselect it; OK;
Refinement reaches after 20 cycles Rp~4.9%
Start “Profile viewer”
[On the screen: Jana 2006 profile viewer]
Press “X exactly” and select the 2θ interval from 19 to 35.
With only
one additional free parameter the fit is significantly better than for the previous methods.
In Jana2006 there is also implemented the fundamental approach following closely the
approach described by Cheary and Coelho [J. Appl. Cryst. (1998). 31, 851-861; J. Appl. Cryst.
(1998). 31, 862-868]. This method was not applied in our example because it was developed
for Bragg-Brentano geometry.
2. Import procedure
Select “Powder data – various CW formats”; NEXT
[On the screen: Powder data form]
Select "PANalytical XRDML"
For method select “Debye-Scherrer”
[On the screen: Complete/correct experimental parameters]
Fill cell parameters 9.19 21.38 8.12 90 113.78 90
Check the radiation (should be interpreted from the xrdml file as Copper with doublet):
“Run Refine”
Refinement of profile slowly converges to the profile R value Rp=1.4 and GOF ~3.0.
“Show Powder profile”
5. Structure solution
[On the screen: information that the solution wizard will start]
OK
Type formula Mg O10 H40 Br2 N8 C14; Formula units 2
Select “Use Superflip”
Check "Repeat Superflip: Number of runs" and choose 10 runs
With this setting Superflip will run ten times and a solution with the best fit will be used
OK
The agreement factors of the symmetry printed at the bottom of the Superflip listing will be
between 10-20, indicating rather worse fit for the screw axis.
Press the button “Draw structure”
The initial structure model should display octahedrally coordinated Mg located in one of
inversion centers, isolated Br, two oxygens of the water molecules, a fragment of the N4C7
cage and probably some false maxima. For instance:
Rietveld refinement with initial structure model could destroy profile parameters:
Start Parameters → Powder
Clear refinement keys for all refined powder parameters in page "Cell", "Profile",
"Asymmetry/Diffractometer" and "Corrections"
The stability of the structure model refined with powder data requires binding together ADP
parameters of some atoms and keeping them isotropic. For instance, carbon and nitrogen
atoms will all have the same isotropic ADP.
Plot the structure and note atoms around Mg
Start “Edit atoms”;
change the atoms around Mg to oxygen O1, O2, O3;
change all other atoms except Br and Mg to carbon.
Start Refinement options, page “Various”
Start “Restrictions” and define:
The same ADP for O*
The same ADP for C*
OK;
In page Basic change number of refinement cycles to 100 and damping factor to 0.5
Clear “Disable atoms having too large isotropic ADP parameter”
Clear “Make only profile matching”
OK; Yes+Start
Rietveld refinement slowly converges with Robs=9%, Rp=6%, and GOF=13.
The numbers are only indicative!
In case of oscillations, note the label of the oscillating atom. Plot the structure and try to
identify and delete false maxima. In case of troubles, repeat charge flipping.
Change ADP of Br to anisotropic and repeat refinement
Br is slightly disordered and it can be described with an ellipsoid even from powder data
For isolated carbon atoms out of the fragment of the cage change their chemical type to
oxygen
Repeat refinement
Two of the isolated positions will be oxygens of water. The others will probably move to the
cage during the refinement.
Start Fourier option
In the page “Basic” define “F(obs)-F(calc) – difference Fourier)
OK; YES+START
OK to start procedure for including new atoms
The charge of difference maxima will be extremely low comparing with single crystal
Run refinement
Rietveld refinement converges with Robs~3.2%, Rp~2.4%, and GOF~5.2.
At this point ADP of nitrogens should be positive while ADP of carbons may still be close to
zero
Start “Edit atoms”
Add hydrogen atoms to carbons of the cage (automatic mode)
Run refinement
Rietveld refinement converges with Robs~3.2%, Rp~2.5%, and GOF~5.4.These value are
slightly worse than in the previous step but ADP of carbon atoms become positive
Start File → Structure → Save As and make a backup copy
Start Parameters → Powder
Activate refinement of all powder parameters which were fixed before starting with the
Rietveld refinement
Run refinement
Rietveld refinement converges with almost the same Robs, Rp and GOF.
Run refinement
Rietveld refinement converges with Robs~2.5%, Rp~2.0%, and GOF~4.0.These value are
significantly better.
Oxygen atoms of water, which is not coordinated to Mg, may have different ADP
Start refinement options, page Various
Start "Restrictions"
Delete the restriction for O*
Define separately the same ADP for O1, O2, O3 and for O4, O5:
Run refinement
The results are almost the same. Visual inspection of the powder profile shows good fit
between the calculated and observed lines.
The preparation of M45 has no special support in the current version of Jana2006. Part of the
job must be done with a text editor.
In the following steps we shall use the CSD structure HICTOD and we shall prepare two versions
of the cage: with C1 symmetry and Cs symmetry.
4. Transformations
M45 templates contain information about the point group symmetry. The axis c is used as
the unique axis of the point group. For this reason the unit cell taken from CSD must be
transformed in order to get the mirror plane of the cage perpendicular to c and located at
z=0. Moreover, all atoms of the cage (not only symmetry independent atoms) must be
present in the template.
Start "Tools → Transformations → Origin shift"
Move the origin to point 0, 0.25, 0
Now the mirror plane of the cage is located at y=0 (can be verified in M40)
Start "Edit Atoms"
Press "Select all"
Start "Action → Expand selected atoms"
Press "Next matrix" ; "Fill by a symmetry operation"
Select "x, -y, z"; OK;
OK; OK; YES to rewrite changes
Start "EditM50"
Change symmetry to P1; OK; YES to rewrite M50
Start "Edit Atoms"
Press "Select all"
Start "Action → Merge selected atoms"; use default merging distance;
The program merges 7 atoms, resulting in 25 atoms of the cage
Start "Tools → Transformations → Cell transformation"
Press "by a matrix"
Define transformation matrix -1 0 0 / 0 0 1 / 0 1 0
OK; OK; select "rewrite the old structure"; OK
Start "Edit Atoms"
Press "Select all"
Start "Action → Transform selected atoms"
Use unit transformation matrix and translation vector 1, 0, 0
This transforms the x coordinates back to the unit cell
42 Jana2006 Cookbook, version October 2017
Start "Plot structure"; "Draw + continue"
Use , and to plot atoms from M40
Use Picture → Viewing direction → b
Use and to display the unit cell
This shows the resulting cage of 25 atoms with mirror plane located at z=0
5. Labels
Quit Diamond
Start "Edit Atoms"
Double-click the first atom
Using "Define" and "Edit" page change the suffix to "m" if z=0, "1" if z<0 and "2" if z>0
Arrows close to the "List" button can be used to move between atoms in the editing mode
Renaming can also be done in a text editor (Edit/View → Edit of M40 file).
The resulting labels will have "m" for atoms at the mirror plane, "1" for atoms below the
mirror plane and "2" for atoms above the mirror plane. The suffixes "1" and "2" resulted
from the previous expansion by the mirror plane.
The keyword scdist is used for scaling of the complete object using arbitrary number
described with arbitrary identifier (without spaces). In our case the identifier is “C-C” and the
number is a typical C-C distance in the cage. But the identifier could be also “scale” and the
number could be 100.
Each line except header lines contains the following items: label, number of chemical type
(indicates that chemical types are different - they will be specified during application of the
rigid body in Jana2006), unused, site occupancy, x, y, z, symmetry identifier (0 means
symmetry independent with respect to the used point group).
2. Import procedure
Select “Powder data – various CW formats”; NEXT
[On the screen: Powder data form]
Select "PANalytical XRDML"
For method select “Debye-Scherrer”
[On the screen: Complete/correct experimental parameters]
Fill cell parameters
Check the radiation (should be interpreted from the xrdml file as Copper with doublet):
5. Structure solution
[On the screen: information that the solution wizard will start]
OK
Type formula Mg O10 H40 Br2 N8 C14; Formula units 2
Select “Use Superflip”
Check "Repeat Superflip: Number of runs" and choose 10 runs
With this setting Superflip will run ten times and a solution with best R value will be used
“Run Solution”
In case that Superflip does not converge interrupt the calculation, discard the results and
repeat it. It will start from a new set of random phases and the result may be different. The
agreement factors of the symmetry printed at the bottom of the Superflip listing will be
between 10-20, indicating rather worse fit for the screw axis.
Plot the cage and compare with the template prepared in the previous example. Probably
C9 corresponds to C1_m; C3 corresponds to C2_m; C8 corresponds to C4_m (your labels
can be different!)
Return to Jana2006
Start “Parameters → Molecules → New molecule”
[On the screen: Atoms of the new molecule from:”]
Select “Model file”
For name of the molecule type “kage”
We need a name not starting with C, N, O, B, H and M
For “Model name” browse for “Magcag-symmC1.m45” created in the previous example
Keep the scaling distance unchanged
“Show the model molecule” is disabled because we don’t have Magcag-symmC1.jpg
NEXT
[On the screen: Specify atomic types in the model molecule]
Do not change the suggested atomic types; NEXT
[On the screen: Define the molecular reference point]
Select “Explicit”
In “Reference point” textbox type N1_m
The reference point should be some atom from the mirror plane in order to reflect the
symmetry of the molecule. Such atoms in the model molecule have suffix _m (see Example
2.3.2)
NEXT
[On the screen: Define and complete molecular position #1]
For the pairs “model atoms” – “actual positions” choose c1_m and C9; C2_m and C3; C4_m
and C8 (with actual names taken previously from Diamond)
In this example one water oxygen is too close to the cage but it may change during the
refinement
Return to Jana2006
Start “Refine commands”; go to page “Various”
Press “Fixed commands”
Fix “all parameters” for C*, H*, N* (i.e. for atoms of the cage)
Because the molecule is called “kage”, the wild character C* does not include the molecule
and therefore molecular parameters (rotations and translation of the cage), will not be
fixed.
OK; OK; YES+START to save changes and start refinement
During refinement you may change damping from 0.1 to 0.5
After many cycles refinement converges with R~4.5%, Rp~4%, GOF~9
Start “Plot structure”
The cage should be now in better position, i.e. without short distances to lattice water
All oxygens belong to water; refinement of their occupancy mimics the missing hydrogens
OK
Go to page “Various”, press “restrictions”
Define the same ADP for all oxygen atoms (O*); OK
Be careful not to make identical coordinates!
OK; YES to save commands without starting Refine
Start “Edit atoms”
Set harmonic ADP for Br1
This works because Br is slightly disordered
Activate refinement of occupancy (ai) of O1
OK, YES to save changes
Start “Refine”
The refinement converges with R~3%, Rp~2.3%, GOF~5
Final plot:
Start "Plot structure" and try to orient your cage similarly. Make use of the fact that there is
only one C-C bond in the cage.
Press “Apply+End”
Start "Edit atoms"
The atomic part should contain only Br1, Mg1 and five oxygens
Start "Plot structure"
The cage should be correctly placed – verify by plotting
In a plot slash "/" means the atoms is transformed by the point group symmetry of the
molecule, while "a" indicates the first actual position of the molecule.
Input files:
pfphenyl.xrdml with powder profile
pfphenyl.cif with solution obtained by the program FOX
Additional information:
unit cell 10.77024 9.957075 26.70771 90 91.683 90
expected symmetry P21/n
composition C30 N3 H29 P O2 F3 (Z=4)
Chemical scheme:
2. Import procedure
[On the screen: Specify type of the file to be imported]
Select “Structure – from CIF”; NEXT
[On the screen: select input CIF file]
Select pfphenyl.cif; OK
[On the screen: CIF does not contain any reflection block]
YES to import data from file
[On the screen: Specify type of the file to be imported]
Select "various CW formats"; NEXT
[On the screen: Powder data from:]
Select "PANalytical XRDML"
For method select “Debye-Scherrer”; NEXT
[On the screen: Complete/correct experimental parameters]
NEXT; FINISH; OK to close data repository; YES to accept the data set
4. Verification of symmetry
Start File → Reflection file → Make space group test
YES to continue
[On the screen: Tolerances for crystal system recognition]
Use default; NEXT
[On the screen: Select Laue symmetry]
Select the highest one – monoclinic; NEXT
[On the screen: Select cell centering]
(calculation takes time! It is finished when the program prints warning about cell centering
at the end of the list of centering symbols).
The table shows Rp values calculated for various centering vectors. For instance, for A
centering the program discards corresponding extinct Bragg reflections and calculates Rp. If
A centering would really exist, the new Rp would be similar like the one for primitive unit cell
because the discarded forbidden reflections would not be vital for the profile fit. In our case,
A centering is not possible because the corresponding Rp is much worse than for the
primitive unit cell.
Select P centering as it gives the best profile fit; NEXT
[On the screen: Information about progress of space group tests. It takes some time…]
[On the screen: Select space group]
The space group P21/n, which is the space group used for the FOX solution has good FOM
(figure of merit). This space group has been already entered with the CIF structure so we can
cancel the process at this point. We could also consider P2/n but we will rely on the FOX
symmetry.
CANCEL
Before refinement we have to fix geometry for hydrogen atoms. Hydrogen positions cannot
be determined from this powder data set.
Return to Jana
Start "Edit atoms"
Select all carbon atoms and N9
Start "Action → Adding of hydrogen atoms"
[On the screen: Adding of "hydrogen" atoms"]
Use default settings; OK
Hydrogen atoms of carbon have been fixed automatically while for N9 the program offers a
dialogue
[On the screen: Adding "hydrogen" atoms for "N9"]
Select "Trigonal"; C28 and P1 for neighbors; H54 for hydrogen
Press "Apply → Go to next"
OK; YES to save changes
Start Refinement options (Right-click on the Refine icon)
Select the page "Basic"
Set 0 refinement cycles
Uncheck "Make only profile matching"
This changes le Bail fitting into Rietveld refinement, where intensities of the profile points
are calculated from the structure model
We can see that the PN3O group has been refined correctly but phenyl groups are not
planar and the CF3 group has been destroyed
6. Fixing geometry with keep commands
Start "Tools → Recover files"
This makes one step back. Otherwise "File → Structure → Copy in" would be used for
importing the previously created backup file.
Start "Plot structure" in order to see atom labels.
With "Draw+continue" Jana2006 and Diamond can be used together and atom labels seen
The old keep commands acting for atoms replaced by the actual molecular positions, except
of those selected as a model molecule, were automatically deleted. The atoms will be still
kept in ideal position as before but now through the keep commands applied to the model
molecule.
However, we have also used keep commands involving atoms from the phenyl groups and
N-CH2- group. The command defining positions of two hydrogen atoms bonded to C19 has
not been deleted, and it still contains the atom C16, which belongs to the model molecule
and not to some actual position. Therefore, C16 must be changed to C16a:
The commands for keeping the geometry of hydrogen atoms in N-CH2- groups for remaining
three atoms (C55, C15 and C13) have to be re-defined:
OK
Open the page "Basic"
Set number of cycles to 0; OK; YES+START to save commands and start refinement
This will apply keep commands without any further changes
Start "Plot structure"
Verify the structure after application of the new keep commands
There should be only small changes
Start "Refine Commands"
Change number of refinement cycles back to 100
OK; YES+START to save commands and start refinement
The refinement converges with Rp~2.6%, R~7.8% and GOF~2.4.
Start "Plot structure"
The resulting structure has much better geometry of phenyl groups
However, the distance C13-N8 is too large
Room temperature (RT) phase of this sample is known from single crystal data. At cooling, a
phase transition occured which repeatedly destroyed the single crystals. Therefore, the low
temperature (LT) phase was measured at 150K as a powder. The powder sample at 150K still
contained an admixture of the room temperature phase.
Unit cell of the LT phase was determined with CRYSFIRE 2004. Its structure was solved by direct
methods implemented in the program FOX, using the RT structure as a starting template.
Input files:
twoph.xrdml with powder profile
RT.cif crystal structure of the RT phase solved and refined from single-crystal
LT.cif crystal structure of the LT phase solved from powder data by FOX. Not refined yet.
Additional information:
Powder measurement was done with Empyrean (PANalytical) using CuKα radiation, focusing
mirror and Debye-Scherrer geometry (sample in a capillary)
RT phase unit cell: 6.81 8.56 17.78 90 96.8 90
RT phase symmetry: P21/n
LT phase unit cell: 8.54 17.68 7.15 90 111.25 90
LT phase symmetry: P21/a
2. Import procedure
[On the screen: Specify type of the file to be imported]
Select “Powder data – various CW formats”; NEXT
6. Le Bail refinement
Start "Edit profile"
[On the screen: Powder options for phase: RT]
Right-Click on “Refine”
[On the screen: Refine commands]
Set "Number of cycles" to 30
OK; Yes+start
Refinement converges with GOF~3.3%, Rp~1.27%, 37 profile parameters
8. Backup copy
The Le Bail fitting is now complete, create a backup copy!
Start “File → Structure → Save as"
Make a backup copy "twoph_lebail"
OK; NO (Do not continue with the new structure)
Right-Click on Refine
[On the screen: Refine commands]
Set "Number of cycles" to 100
Uncheck "Make only profile matching"
Go to "Various" and open "Fixed commands"
Fix all structure parameters (“*” for Atoms/parameters; “Add”)
OK; OK; Yes+start
Refinement converges with Rp ~ 6%
OK; Yes;
Make a backup copy (File->Structure->Save as)
Start refinement
Press "Use all" as a reply to the message about significant dicsrepancy
Refinement converges with R~4.7% (RT), 4.3% (LT), Rp ~ 3%
The CIF for the RT contains unit cell determined at room temperatire. However, our powder
pattern was measured at 150 K. The unit cell parameters are changing with temperature and it
may happen that the change it so large that refinement from the room temperature unit cell
does not converge.
We can estimate the deviation of the RT unit cell from the powder profile before its refinement:
Right-click on Refine
[On the screen: Refine commands]
Set Number of cycles to 0;
Check "Make only profile matching" to perform only Le Bail refinement.
OK; Yes+start; (This calculate theoretical pattern)
Go to the "Profile viewer"
Display the range from 9° to 19° 2Θ.
We can see that while for LT Bragg positions fit well, RT Bragg positions (second line in the plot)
are deviated, e.g., 0 1 1, -1 0 1, 0 1 3.
For rough estimation of the proper starting point for refinement of the RT unit cell, we can
assume that beta does not change with temperature and we can calculate a, b, c from
observed positions of d011= 7.7 Å, d013 = 4.85 Å and d-101 = 6.58 Å: a = 6.84 Å, b = 8.55 Å, c =
17.77 Å.
For obtaining d, use “Draw in d’s”, find reflection hkl by “Go to” button and find its real position
by “Info” line (it must stay a while at the peak position).
Bis[N-(2-benzylidenepropylidene)phenyl]ether
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction → CCD”; NEXT
Select “Input from “sum” file”; OK
Leave all settings unchanged; NEXT
Leave all settings unchanged; NEXT
For absorption correction select “None or done before importing”; NEXT
FINISH
3. Data Repository
The line in data repository window corresponds to one data set.
OK; YES to accept the data set
The results may be slightly different for each run of Superflip. This is because Superflip
starts from random phases and the resulting electron density map differs for each run.
Compare the plot with the chemical scheme and note labels of atoms, which must be
oxygen or nitrogen. In our example, C2 and C3 should be nitrogen and C1 should be
oxygen.
Quit Diamond
9. Refinement
[On the screen: basic window of Jana]
Right-click the icon of Refine.
[On the screen: refinement options]
Define 100 of cycles, Instability factor 0.02; OK
Choose “YES+START”
Refinement converges with R value about 16%, GOF about 8, 141 refined parameters
Right-click Fourier; select “F(obs)-F(calc) - difference Fourier”; make sure that “Weighting
of reflections” is checked and “Apply sin(th)/lambda limits” unchecked; in the page
“Scope” press “Reset to default”; in the page “Peaks” select “in a fixed sphere of radius
0.8”
OK; YES+START; NO to include new atoms; NO to see the listing
This figure shows diffraction pattern of pyNinit – peak positions projected along b* with unit
cells of both twins related by the above-mentioned twinning matrix
91 Jana2006 Cookbook, version October 2017
1. Creating new jobname
Start Jana2006
“File → Structure → New” opens a file manager
Left pane: locate directory with input files
Right pane: double-click pyNinit_twin1
3. Data Repository
The line in data repository window corresponds to one data set.
OK; YES to accept the data set
4. Symmetry Wizard
[On the screen: an information that symmetry wizard follows]
NEXT
[On the screen: “Tolerances for crystal system recognition”.]
Leave all settings default; NEXT
[On the screen: Select supercell]
Select continue with the basic cell; NEXT
[On the screen: Select Laue symmetry]
Select Monoclicnic; NEXT
[On the screen: Select cell centering]
Select primitive unit cell; NEXT
[On the screen: Select space group]
Using the criteria described in the first example Select P21/c; NEXT
[On the screen: Final step of the space group test]
Accept the space group transformed into the original cell; FINISH
Follow the wizard for creation of the refinement reflection file
5. Structure solution
[On the screen: window of Structure solution]
In “Formula” textbox check list of chemical elements for PyNinit: Ni C N O H
Select “Superflip”, “Peaks from Jana2006”; leave other settings default; “Run Solution”
Superflip converges (after noise suppression) with R value ≈ 23%.
Fit of symmetry elements is very good for c, slightly worse for 21
Press CLOSE to leave the listing
Press “Draw structure” to see the result
Press “Accept last solution” to close the Structure solution wizard
7. Twinning
In “EditM50” check “Twinning” and define twinning matrix (-1 0 -0.733| 0 -1 0| 0 0 1)
The twinning matrix is known from the diffractometer software and represents a rotation by
180º around c. The program offers re-creation of the refinement reflection file because
twinning may change rules for merging of symmetrically equivalent reflections. In this case,
it is not necessary.
In “Parameters → Twin fractions” enable refinement of “twvol2”
Repeat refinement
Refinement converges with R value around 13% and GOF 8 for 317 parameters
Open listing of Refine
ADP tensors are still not positive definite for some atoms.
This twin has a diffraction pattern where some reflections are partially overlapped. In
refinement options there are two limits available for treating partial overlaps. The
reflections with difference in theta smaller than “Maximal angular difference for twin
overlap” are treated as fully overlapped; the reflections with theta difference larger than
“Minimal angular difference for full separation” are treated as fully separated; the
reflections between the limits are discarded. The limits should be adjusted carefully not to
93 Jana2006 Cookbook, version October 2017
discard too many reflections.
In Refinement options change values of “Maximal angular difference for twin overlap” (L1 in
the table) and “Minimal angular difference for full separation” (L2 in the table) and watch
their impact on R factors, number of excluded reflections and ADP tensors. Find an optimal
combination.
With maximal angular difference for twin overlap 0.1 and minimal angular difference for full
separation 0.4, refinement converges with R value 3.9% and GOF 1.4, all ADP tensors are
positive definite, the fraction of the second twin domain refines as 0.468(2) (see
“Parameters → Twin fractions)”.
Using “Structure → Save as”, create a backup copy, for instance pyNinit_1domain
In this step the program offers sorting of the data sets to data blocks. Each data block has
an independent scale factor. In our case, both twin domains come from the same
measurement and they should be used as one data block.
Set Data collection block #1 for both datasets
Follow the wizard for creation of the refinement reflection file; select “on the same scale”
for both data sets
Refine
Refinement converges with R 4.1%, GOF 1.36, 3800 excluded reflections, fraction of the
second domain is 0.4641(6) (the number with s.u. is printed in the refinement listing).
Create a backup copy, e.g., pyNinit_2domains
The scale is calculated from 4334 pairs of equivent reflections, which should have the same
intensity in both data sets after the scaling.
Refine
Refinement converges with R about 4.1%, GOF 1.39, 3800 excluded reflections, fraction of
the second twin domain is 0.4716(6). However, in this refinement we have unrealistic weight
of the duplicate reflections because they were entered two times with the same sigma.
NEXT;
OK to confirm information about 68907 imported reflections
FINISH
OK; YES to close and save Data repository
Follow the wizard for averaging of reflections and creation of the refinement reflection file
Jana2006 merges only fully separated reflections from both domains
Refine
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction-CCD”; NEXT
Select “Input from “sum” file”; OK
Leave all settings unchanged; NEXT
Leave all settings unchanged (leave Twinning checkbox clear); NEXT
For absorption correction select “None or done before importing”; NEXT
FINISH
3. Data Repository
The line in data repository window corresponds to one data set.
OK; YES to accept the data set
4. Symmetry Wizard
[On the screen: an information that the next tool is the symmetry wizard]
NEXT
[On the screen: “Tolerances for crystal system recognition”.]
Leave all settings default; NEXT
[On the screen: Select Laue symmetry]
For exactly hexagonal cell parameters, merohedric twinning is highly probable. Another
indication of twinning is the fact that the three equivalent orthorhombic settings have quite
different Rint. It might be caused by merging reflections belonging to different unequally
populated twin domains.
In the next step, we will repeat test of Laue symmetry but this time considering twinning.
That means for any subgroup of 6/mmm the program will generate corresponding pseudo-
merohedric twinning matrices and it will use only symmetry elements which are present in a
twinned structure for general case of not equally populated twin domains.
Press BACK
Check „Introduce twin laws in case of subgroups“
NEXT
This time, the three orthorhombic cases gives the same internal R factor.
Select the highlighted orthorhombic Laue symmetry; NEXT
The three orthorhombic cases are related the three twin domains. The selected one gave us
in the previous test the lowest Rint. This means that this domain is more populated than the
others.
Select C centering; NEXT
Hexagonal unit cell transformed to an orthorhombic one must be C centered
[On the screen: Select space group]
Select Ccce space group; NEXT
There are several possibilities for selecting space group and at this stage it is unclear which
is the correct one.
Accept the space group in the standard setting; FINISH
Follow wizard for creation of the refinement reflection file
QUIT the wizard for structure solution without solving the structure
eformat91 0
datblock Block1
nref 2021 obslim 3 sctom90 1 norefitems 1 diffscales 1 dataave 1 scalelim
10 hklf5 0
indslowest 3 indfastest 1 addcentrsymm 0 sigimethod 1 multave 1
flimprint 5 flimcull -1 avesigwt 0 unstab 0.022791
radtype 1 dattype 1
lambda 0.71073 lpfactor 1 monangle 6.0821 perfmono 0.5
datcolltemp 293
end
Data Block1
2 0 0 166.8 3.4 1 0 1 0.0000
4 0 0 115151.2 357.3 1 0 1 0.0000
6 0 0 13.8 9.2 1 0 1 0.0000
8 0 0 56835.4 139.6 1 0 1 0.0000
10 0 0 -8.5 16.5 1 0 1 0.0000
The numbers indicated in bold are twin flags. All flags except several reflections at the end
of the file have the twin flags equal to “1” because most reflections can be indexed in the
first domain. The reflections at the end of the file with the twin flag “2” and “3” are the
ones, which cannot be transformed to the first domain because after such transformation
they would break extinction rules for the selected space group.
7. What is wrong?
“File → Transformation → Cell transformation”
Press “to those from data collection”; OK
Select “rewrite the old structure”; OK
Follow steps of the wizard for creation of the refinement reflection file
We are now back at the original hexagonal unit cell. In EditM50, you can see there is no
symbol for the transformed symmetry operation. However, we will check for the symmetry
again.
“File → Reflection file → Make space group test”
This manually starts the Symmetry wizard.
[On the screen: Tolerances for crystal system recognition]
Check „Introduce twin laws in case of subgroups“
Select “look for centering vectors composed from 0 and ½”
This last setting is different compared with the previous symmetry check. Normally, one
looks for vectors 0, 1/3, 2/3 when the hexagonal unit cell is present, but we already know
that we will use an orthorhombic symmetry where 0 and ½ centering vectors are
reasonable.
NEXT
Select the highlighted orthorhombic Lause symmetry (like before); NEXT
[On the screen: Select cell centering]
The transformation suggested at the final step of the space group test produces an
orthorhombic unit cell with four time smaller volume compared with the previously tested
orthorhombic unit cell. Such unit cell alone could not index all diffraction spots, but
combined with a three-fold twinning it can. The indexing of the hexagonal diffraction
pattern with a three-fold twin of this small orthorhombic unit cell is indicated in the next
figure.
Run refinement
Start “Edit atoms”; change ADP of all atoms to harmonic
Run refinement
Refinement converges with GOF ~ 4.4 and R(obs) ~ 4.4% for 36 parameters
Right-click Fourier
Select “Map type” = “F(obs)-F(calc) – difference Fourier”
OK; YES+START; YES to include new atoms
[On the screen: Inserting/replacing of atoms]
OK
[On the screen: List of peaks]
The strongest difference peak (Max1) has reasonable Li-O distances to oxygen atoms. We
can add it as Lithium in the first highlighted symmetry position.
Press “Include selected peak”
For “Name of atom” type Li1; TAB; OK; FINISH; YES to include the new atom.
Run refinement
Refinement converges with GOF ~ 3 and R(obs) ~ 2.7% for 39 parameters
106 Jana2006 Cookbook, version October 2017
YES (or Edit/View → View of Refine) to open the listing
Go to → Statistics Fo, sin(th)/lambda after refinement
The Statistics show that the strongest low-angle reflections are systematically weaker. This
is an indication that extinction correction is needed.
Close the refinement listing
Using “Parameters → Extinction” introduce isotropic Type 1 Gaussian or Lorentzian
extinction correction
Run Refine
Refinement converges with GOF ~ 2.7 and R(obs) ~ 2.2% for 40 parameters
Start “Edit atoms”; change ADP of Lithium to harmonic
Refinement converges with GOF ~ 2.6 and R(obs) ~ 2.2% for 43 parameters
Right-click Refine; go to the page Basic
Check “Instability factor from reflection statistics”
OK; YES+START
Refinement converges with GOF ~ 1.4 and R(obs) ~ 2.3% for 43 parameters
This last step has changed the weighting scheme – see Example 1.1 for details
H. Kruger, V. Kahlenberg, V. Petříček, F. Phillipp & W. Wertl, J. Solid State Chem., 182, (2009),
1515-1523.
Main features of the solution for the HT phase: The atoms Fe2 and O3 are occupying the
positions close to the Wyckoff position 4(e), point group mm2. The mirror m, with normal along
a axis, having half shift along the internal x4 axis, transforms the “left” oriented ribbon of Fe
tetrahedrons to the “right” ones and the occupational modulations of both atoms have a crenel
like shape.
During the first order phase transition both phases (LT and HT) are coexisting in a wide region
of temperatures and with x-ray diffraction they both contribute to the same diffraction pattern.
In this example we shall start from the known structures of both phases and we shall
demonstrate that the generalized twinning option in Jana2006 can also be applied to single
crystals composed from two or more phases.
As a starting point you have at disposal the HT structure refined against the data set having the
highest possible fraction of the HT phase from all measured data sets. In the first part of the
example it is demonstrated main features of the solution. In the second part the LT phase will
be read in from the supplied CIF file and then the final refinement will be performed.
There are 381 observed reflections, some of them relatively strong, which do not follow the
extinction rules for Imma(00γ)s00. This means that the reflections of the second
“P-centered” phase are really present.
110 Jana2006 Cookbook, version October 2017
3. Run refinement
[On the screen: basic window of Jana]
Right-click the icon of Refine.
[On the screen: refinement options]
Define 100 of cycles, Damping factor 1; OK
Choose “YES+START”
The refinement converges with R(obs) about 5.98%, 4.48% and 9.12%.
5. Run dist
[On the screen: basic Jana window]
Right-click the icon of Dist.
[On the screen: distance commands]
Use the default parameters except for the listing form which should be changed to “With
symmetry code”.
OK; Start the program;
Open the listing and look for Fe2 coordination;
The atom Fe2 is coordinated by the atom O2, O2#s2, O3 and O3#s4c2t0,0,1.
Repeat it for the center at -0.5 0.25 0.625 (a position close to O3);
Press the button “Atoms edit” to show the modulation curves for O3 and the symmetry
related atom O3#s8t-1,0,0. The translation part -1,0,0 returns the symmetry related atom
The crenel function for atoms Fe2 and O3 assures the alternation between the “Left” and
“Right” configurations of Fe tetrahedrons. During the refinement the “delta” values were
fixed to 0.5 for both atoms. But there were no restrictions applied to x40. Now we have to
check:
(a) if there is no overlap between two symmetry related positions of Fe2 and Fe2#s8 and O3
and O3#s8t-1,0,0.
The symmetry operation #8 has the following form: -x1, x2, x3, x4+1/2
for the atom Fe2 we have: t min = x40 − ∆ 2 − q ⋅ r t max = x40 + ∆ 2 − q ⋅ r and for the
symmetry related Fe2#s8: t min = x40 + 1 2 − ∆ 2 − q ⋅ r ′ t max = x40 + 1 2 + ∆ 2 − q ⋅ r ′
From the fact that the modulation vector is (0,0, γ ) and that the third component of the
symmetry related atoms remains the same for both positions, it follows that the t intervals
are fully separated with no overlaps.
The same conclusion is valid also for the atom O3.
(b) if the t intervals in which exist atom Fe2, O3 and O3#s4c2t0,0,1 are identical.
This means that the refined value was really close to the ideal value
For O3#s4c2t0,0,1 the symmetry operation s 4c 2t 0,0,1 is x1+1/2 -x2+1/2 -x3+3/2 -x4+1/2
modifies tmin to:
z (O 3# s 4c 2t 0,0,1) = 3 2 − z (O 3) = 0.825
Numerically:
x40 (O 3# s 4c 2t 0,0,1) = 1 2 − x40 (O 3) = −0.1755
This means that: t min = −0.1755 − 0.6101 ⋅ 0.875 − 0.25 = −0.959 and therefore we have to
modify the equation by adding one:
and by applying the symmetry operation s 4c 2t 0,0,1 ~ x1+1/2 -x2+1/2 -x3+3/2 -x4+1/2:
Combining (2) with the previous one we have for x40 (Fe 2 ) :
x40 (O 3) = 3 4 ⋅ (1 − γ ) + γ ⋅ z (O 3) ⇒ t 40 (O 3) = 3 4 ⋅ (1 − γ )
x40[Fe2]=0.2924+0.6101*z[Fe2]
x40[O3]=0.2924+0.6101*z[O3]
[Cu(CH3-NH-CH2-CH2-NH2)2] [Pt(CN)4]
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction-CCD”; NEXT
[On the screen: “Define basic input file”]
Use "input from sum file"; OK
[On the screen: “Your sum file consists of 6 appended data reduction runs”]
OK;
[On the screen: “Complete/correct experimental parameters”]
Leave all settings unchanged; NEXT
[On the screen: “Define the reference cell/split by twinning”]
Leave all settings unchanged; NEXT
For absorption correction select “None or done before importing”; NEXT
FINISH
3. Data Repository
The line in data repository window corresponds to one data set.
OK; YES to accept the data set
4. Symmetry Wizard
[On the screen: “Tolerances for crystal system recognition”]
Leave all settings default; NEXT
[On the screen: Select Laue symmetry]
Select Tetragonal Laue symmetry 4/m; NEXT
Select primitive unit cell; NEXT
[On the screen: Select space group]
Using the criteria described in the first example Select P42/m; NEXT
Select "Accept the space group transformed into the original cell"
Clear “Accept twinning matrices induced by the space group test”
117 Jana2006 Cookbook, version October 2017
The highest Laue symmetry consistent with the unit cell parameters is 4/mmm while the
selected Lause symmetry is lower – 4/m. This opens possibility for twinning, but in this case
twinning does not occur (as one would find during refinement of twin volume fractions).
FINISH
The original cell is used because the suggested transformation does not change the unit cell
and the space group symbol
YES to create refinement reflection file M90
OK
[On the screen: list of peaks]
In the left panel, Max1 is focused. Its integrated charge is displayed at the bottom. In the
right panel, we can see distances of symmetry equivalent positions of Max1 to already
existing atoms of the structure model, up to the previously defined limit 4Å. We can expect
distances ~2 Å for Pt-C, ~1.5 Å for C-C, ~1 Å for C-N of cyano group and ~3 Å for Pt-N.
Therefore, the maximum Max1 is probably nitrogen.
In the left panel, focus Max1 (or other peak if in your case Max1 has no good distance)
In the right panel, focus its reasonable position having good distance to some existing atom
Press “Include selected peak”
[On the screen: Complete information for the new atom]
For “Name of atom” type C* (even for nitrogen); TAB; OK
We will intentionally add all atoms like carbon, it will help later on to analyse disorder.
The chemical type was assigned to carbon, the name will be C+sequence number.
Repeat the procedure for other maxima. Skip distances below 1 Å as well as maxima with
too small charge or strange distances. At this point, all atoms should be added as carbon,
we will distinguish between C and N later on.
FINISH; YES to include the new atoms
Finally the structure model should contain 9 symmetry-independent atoms. In case some
improper maxima have been added, they will be probably discarded during refinement
based on large isotropic ADP.
Run Refine
Refinement converges with R value about 4%
Start “Edit atoms” and check if the button “Select rejected” is enabled.
If “Select rejected” is enabled some atoms have zero occupancy. It was set by Refine because
their Uiso was refined above the limit defined in the Refinement options (default limit is 0.2).
Delete rejected atoms.
Note names of atoms (C1, C5 in our example) which should be nitrogen of cyano groups
Leave Diamond and change chemical type and names of C1, C5
Run Refine
8. Solution of disorder
Open “Edit atoms” and define harmonic ADP for all atoms
Double-click is necessary in the locked yellow-indicated options
Run Refine
Refinement converges with R value ~3% for 60 parameters
Plot the disordered part of the structure: press ; click on Pt; Ctrl-M selects the Pt
“molecule”; Delete deletes it – remaining atoms are coodinated to Cu
Using “Action → Rename selected atoms manually”, rename atoms according to this plot
This is not necessary but will help to follow easily the cookbook
Start “File → Structure → Save As” and create a backup copy, for instance PtCu_completed
Open “Edit atoms”, select C6
“Action” → “Split atomic position”
[On the screen: split atomic position window]
Leave all settings default; Apply
[On the screen: list of atoms in “Edit atoms”]
Select C6, go to “Edit/define atoms”, switch to the page “Edit”, activate refinement of
occupancy (the parameter ai)
120 Jana2006 Cookbook, version October 2017
Close "Edit atoms", YES to save changes
The next steps show automatically created restriction for C6, C6’
Open Refinement options, go to the page “Various”, press “Restrictions”
Focus the line “restric C6 2 C6’
Press Load
In the lower part of the window there are listed properties of the restriction. The atom C6
was split to two positions, C6 and C6’. Both atoms will have the same ADP and their
occupancy will be complementary.
Leave Refinement options
Run Refine
Refinement converges with R=~2% for 59 parameters
Plot the same part of the structure
C2 belongs to methyl
C2’ and C6 belong to CH2
N4 belongs to NH2
N6 belongs to NH
123 Jana2006 Cookbook, version October 2017
Open “Edit atoms”
Select C2, C2’, C6, N4 and N6
Start “Action → Adding of hydrogen atoms”
Clear “Try automatic run”; OK
[On the screen: Adding hydrogen atoms for C2]
C2 belongs to methyl, therefore it has one neighbor and three hydrogen atoms
Select “Tetrahedral”
Change number of neighbors to 1
Click in the textbox for the neighbor: this activates the button “Select neighbors” in the
lower right corner
Press “Select neighbors”
Select N6; OK
In the textboxes for hydrogens type H1C2, H2C2 and H3C2
Press “Apply → Go to next”
[On the screen: Adding hydrogen atoms for C2’]
C2’ belongs to CH2, therefore it has two neighbors and two hydrogen atoms
Select “Tetrahedral”
Change number of neighbors to 2
Click in the textbox of the first neighbor; button “Select neighbors” will be activated
Press “Select neighbors”
Select N4 and C6; OK
In the textboxes for hydrogens type H1C2’ and H2C2’
Press “Apply → Go to next”
Similarly, continue for C6. For C6 two atoms N6 will be shown as neighbors but the one with
too short distance belongs to the second chain
Similarly, continue for N4 belonging to NH2. Here one of the neighbors will be copper.
Similarly, continue for N6 belonging to NH. Geometry for N6 is trigonal, not tetrahedral.
Leave “Edit atoms”
Check in Diamond geometry of hydrogens (needs manually deleting of unwanted distances
using Table of created bonds)
Run Refine
Refinement converges with R value 1.6% for 67 parameters
At this structure model, hydrogen atoms attached to nitrogen are strictly kept in idealized
geometry. This is correct for carbon but not for nitrogen. In the next step, we will redefine
constraints (keep commands) to restraints.
Start refinement commands; select the page Various; open “Keep commands”
Disable the command “keep hydro tetrahed N4”
Disable the command “keep hydro triang N6”
The commands “keep ADP riding N4” and “keep ADP riding N6” will be retained
OK; open “distance restraints”
Select “Fix to the specified value” and for the “Value” type 0.87
In the textbox “Atoms”, type: N4 H1N4; N4 H2N4; N6 H1N6
Press “Add”; OK; OK; “YES+START”
In case of bad convergence, change the dampling factor to 0.5.
Refinement converges with R value ~1.5% for 76 parameters
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction → CCD”; NEXT
Select “Input from cif_od file”; OK
Leave all settings unchanged; NEXT
Leave all settings unchanged; NEXT
7999 reflections has been imported
For absorption correction select “None or done before importing”; NEXT
FINISH
[On the screen: Data repository]
OK; YES to accept the data set
Start refinement
Refinement converges with R=3.4%
Occupancy of Na1 refines to 0.62
Occupancy of O1 refines to 0.69
From refinement listing we can see there are no warnings about negative ADP
Na1 must be mixed with Ni or Fe because the refined occupancy exceeds the site occupancy.
The BVS larger than 1 support this idea.
For oxygen 0.69*2*8=11; 2 is the site multiplicity, 8 is number of electrons. The result is not
too far of 10 which is sum of electrons for water. Therefore, O1 could be water.
Using “Edit atoms” change chemical type of O1 to Na without renaming
Run refinement
Calculate distances
R value is the same. Occupancy of O1 (now chemical type Na) refines to 0.46. BVS of O1 is
By this way we have created one possible model. Unfortunately there is a lot of another
possibilities which should be checked, for instance:
The water could be partially occupied sodium (it slightly decreases R value to 3.07%) and the
valence of the iron mixed with Ni could be 2+. The composition would be Na1.769 Ni1.27
Fe1.844 P3 O12.
1.844*3+1.27*2+1.769=9.84, i.e. exceed of positive charge
For the above case we could also consider mixed valence for Fe. To get easily the formula we
can temporarily change chemical type of Fe1 to something else, for instance aluminium. But
first it must be added to the end of the formula (to become available in “Edit atoms”).
EditM50 then calculates the formula Na1.769 Ni1.269 Fe1.324 P3 O12 Al0.52.
For Al representing trivalent iron the charge would be 0.52*3*1.32*2+1.27*2+1.769=8.51,
i.e. lack of positive charge
For Al representing divalent iron the charge would be 0.52*2*1.32*3+1.27*2+1.769=9.31,
131 Jana2006 Cookbook, version October 2017
i.e. exceed of positive charge
For the above case with lack of positive charge (sum 8.51) we can also consider mixing of
Na1 not with nickel but with trivalent iron. It yields: Na1.748 Ni1.153 Fe1.325 P3 O12
Al0.656 with trivalent iron represented by aluminium. The charge is
0.66*3+1.33*2+1.15*2+1.75=8.69, i.e. still lack of the positive charge
The Ni1 site is maybe not mixed (BVS is smaller than for Ni2 and mixing has small impact on
R). But iron alone in this site increases R value. It can only be nickel alone yielding the
formula: Na1.716 Ni1.713 Fe1.437 P3 O12. The charge would be 9.46 (trivalent iron), 8.02
(divalent iron), 8.2 or 9.3 (for mixed valence of iron).
It seems that another model with proper charge is not easy to find (maybe it will be
discovered during the workshop). Good proof would be to find hydrogen atoms of water in
the difference Fourier map but they are not there. Or we could find hydrogen in the P2O12
group which would allow the more logical solution with positive charge 8.02, but it is also
not present in difference map.
Conclusions:
or
The latter yields slightly better R value and corresponds better to found bond valences sums.
However, hydrogens cannot be determined from the data.
(Fast microprobe analysis yields Na1.56 Ni1.69 Fe0.78 P3 O6 which is insufficient precision
for making a decision)
Y(PO3)3
Single crystal data measured with Oxford Diffraction four-circle diffractometer
Input files: YPO.hkl, YPO_red.sum
Frame scaling, absorption correction: done
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction → CCD”; NEXT
[On the screen: “Define basic input file:”]
Select “Input from sum file”; file name “ypo_red.sum”; OK
[On the screen: “Warning about appended data reduction runs”]
OK
[On the screen: “Complete/Correct experimental parameters”]
Leave all settings unchanged; NEXT
[On the screen: “Define the reference cell/split by twinning”]
Leave all settings unchanged; NEXT
For absorption correction select “None or done before importing”; NEXT
FINISH
3. Data Repository
The line in data repository window corresponds to one data set.
OK; YES to accept the data set
[On the screen: “Information about starting symmetry wizard”]
NEXT
4. Symmetry Wizard
[On the screen: “Tolerances for crystal system recognition”.]
Symmetry wizard can be started separately by “File → Reflection file → Make space group
test”.
Leave all settings default; NEXT
[On the screen: Select supercell]
Select "Continue with the basic cell"; NEXT
[On the screen: Select Laue symmetry]
Select Monoclinic Laue symmetry; NEXT
[On the screen: Select cell centering]
The plot illustrates that the electron density from Superflip contains reliable information
about position modulation. The modulation curve results by interpretation of the electron
density map by Jana.
Quit Contour
Analogously (starting with button “New plot” followed by “Calculate new ones”), calculate
and plot the section x2-x4 and the section x3-x4:
As the difference Fourier is more sensitive to light oxygen atom we should save the
temporary result for a final check of the refined modulation of that oxygen atom.
Close “Contour”
“File → Structure → Save as” and use a new name for example “YPO-for-maps” but do not
continue with the new structure.
Switch to the page “Edit” and press the button “Occupancy”, which is close to “Edit
modulation parameters”
define “1” for “delta” and “0” for “x40”; OK
This is width and centre of the crenel function
Go to the page "Define"
Select "Legendre polynomials in crenel interval"
OK; OK; Yes to rewrite
Draw Contour section x4-x2 for (-1/4 1/4 -1/2) - like previously
Press "Atoms edit" and add O5; OK
Run “Edit atoms”, select all atoms, press “Action” and select “Edit/Define atoms”
Modify all ADP parameters to harmonic ones by double-clicking the corresponding radio
button
Yellow indicates locked objects. Locking occurs when the corresponding parameters are not
the same (or of the same type) for the selected atoms. In our case ADP is harmonic for most
of atoms but isotropic for the newly added O5, this is the reason for locking of radio buttons.
The text box for site occupancy is also locked because some atoms are in special positions.
Drawing programs such as Diamond cannot draw a modulated structure directly. Jana2006
creates an approximant based on the refined modulated functions. During this process the
superspace symmetry is used to generate all atomic positions and the symmetry is reduced
to P1. The size of approximant and the starting t value are to chosen to visualize different
portions of the non-periodic matter. In our case the most convenient supercell is 1x8x1 as
follows from the modulation vector being close to (0,3/8,0).
“Plot structure”
In Diamond symmetry operations must not be used for the approximant (no automatic
structure, no supercell, no coordinations …. ). This is because Diamond does not know
modulation displacements of newly created positions.
Start “Grapht”
Press New/Edit
For “Parameter to draw” select “ADP parameter”
For “Central atom” use “O2” (the label may differ)
For “Principal” select “U(min)”; OK
The following plot shows that U(min) is negative for a small region of t
Na2CO3
Single crystal data measured with KUMA four-circle diffractometer
(KUMA has changed to Oxford Diffraction, later on Varian, later on Agilent, now Rigaku Oxford
Diffraction)
Input files: naco.hkl, naco.cif_od
Frame scaling, absorption correction: done
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction → CCD”; NEXT
[On the screen: “Define basic input file:"]
Select “Input from cif_od file”; file name “naco.cif_od”; OK
[On the screen: “Complete/Correct experimental parameters"]
Leave all settings unchanged; NEXT
q-vector read from the cif_od file should be: 0.1776 0 0.3252
3. Data Repository
The line in data repository window corresponds to one data set.
OK; YES to accept the data set
[On the screen: “Information about starting symmetry wizard”]
NEXT
4. Symmetry Wizard
[On the screen: “Tolerances for crystal system recognition”.]
Symmetry wizard can be started separately by “File → Reflection file → Make space group
test”.
Change “Maximal deviation for cell angles” to 0.3
Without this change unit cell angles will not be recognized as monoclinic
Leave all other settings default; NEXT
[On the screen: “Select supercell”.]
Leave all settings default; NEXT
[On the screen: Select Laue symmetry]
Select Monoclinic Laue symmetry; NEXT
[On the screen: Select cell centering]
Select C centered unit cell; NEXT
[On the screen: Select space group]
Select C2/m; NEXT
[On the screen: Select superspace group]
C2/m(a0g)0s
[On the screen: Final step of the superspace test]
Accept the superspace group transformed into the original cell: C2/m(a0g)0s; FINISH
Superflip converges (after noise suppression) with R value 21%. It confirms the symmetry
C2/m(a0g)0s.
7. Reflection report
Start “Edit/View → View of Reflection report”
Reflection report shows which symmetry operations were used for merging of reflections
and which extinction rules were applied.
Note that there is no oxygen in M40 – later on some carbon atoms will be changed to
oxygen. The coordinates may differ from the ones shown here dependent on which
equivalent origin position was used by Superflip.
Close the editor
Start “Edit atoms”. It shows the same information in a graphical interface; parameters of
atoms can be accessed by double clicking on atom or by “Action → Edit/Define atoms”;
OK
Start "Plot structure"
Select “Draw average structure”
Start "Draw+continue"
149 Jana2006 Cookbook, version October 2017
With this option the basic 3d structure will be plotted
In Diamond use Ctrl-A and Delete to remove an automatically created plot (if present)
In Diamond use to include all symmetry independent atoms from M40
Delete Na atoms
In Diamond use to add coordinated atoms
Using mouse cursor and rotations identify atoms which should be oxygens of CO3 group:
C1 and C2
Quit Diamond (without saving the plot)
Start "Edit atoms"
Change chemical type of C1, C2 to oxygen
Use “Action → Rename selected atoms to atom type + number” for oxygens and carbon
The structure now contains: Na1, Na2, Na3, O1, O2, C1
Return to the basic window
Repeat refinement
Robs factors: 12%, 7%, 10%, 16%, 23%, 48%
Start “Edit atoms” and add another position modulation wave for all atoms
Repeat refinement
Robs factors: 8%, 8%, 7%, 9%, 15%, 25%
Start “Edit atoms” and add another position modulation wave for all atoms
Repeat refinement
Robs factors: 8%, 6%, 7%, 9%, 14%, 21%
Refinement with four position modulation waves is reasonable because there are satellites
up to the 4th order in data. However, it is not sure at this point that all atoms really need so
many waves.
OK; Yes+start
Fourier calculates a four-dimensional parallelepiped x1–x4–x2–x3 with dimensions 2-?-1-1
Å. In x4, the calculation will run over two periods: from 0 to 2 with step 0.02. The step in
other directions will be 0.05 Å. The two-dimensional sections x1-x4 will be stacked along x2;
the three-dimensional "sections" x1-x4-x2 will be stacked along x3.
Start Contour; press “New plot”
Select “General Contour procedure”, “Use old maps” and “Draw maps as calculated”; OK
[On the screen: the first x1-x4 section]
Press “Sum on” and select summation over x2 and x3
This merges all x1-x4 sections stacked along x2 and x3
[On the screen: merged section]
Press “Atom edit”; for “Atom name” type Na1; Add; OK
[On the screen: merged section with refined modulation function of Na1]
Primary scope means the scope of the coordinate displayed in the plot; secondary scope
means the scope of the other two coordinates. Sections will be automatically summed.
Press “Define atoms for de Wolff’s sections”
Run refinement
Robs factors: 6%, 5%, 5%, 7%, 12%, 20%
Start “Edit/View → editing of M40 file”.
This is the structure parameter file. Each atom has four modulation waves (number of waves
and 24 parameters indicated in red) and – except carbon – one ADP modulation wave
(number of ADP modulation waves and 12 parameters indicated in blue)
According to de Wolff sections, O1 has the second smallest position modulations. Because
ADP modulation usually corresponds to changes in environment of an atom due to position
modulation, we should check whether ADP modulation of O1 is reasonable.
For the remaining atoms we shall add ADP waves one-by-one to see where to stop.
Add the 3rd ADP modulation wave for Na1
Run Refinement
Robs factors: 4.5%, 3.8%, 3.5%, 5.3%, 8.8%, 14.2%
Add the 4th ADP modulation wave for Na1
Run Refinement
Robs factors: 4.5%, 3.8%, 3.5%, 5.2%, 8.6%, 14.3%
Remove the 4th ADP modulation wave for Na1
Run Refinement
Na1 4 3
Na2 4 2
Na3 4 2
O1 4 2
O2 4 4
C1 3 0
The shell of the strongest reflections has R factor 2.66% while the next strongest shell has
2.32%.
Close the refinement listing
Start “Parameters → Extinction”
Select Isotropic, Type1, Gaussian; OK; Yes to rewrite changes
Run refinement
Robs factors: 4.1%, 3.6%, 3.2%, 4.7%, 8.1%, 12.0% (for 193 parameters)
Check refinement listing
The partial R factor for the strongest reflections is better (2.48%) but still slightly worse than
the neighboring shell.
O1 O2 C1
Use “File → Structure → Save As” and create backup copy “naco_final”
This plot allows for estimation which supercell we shall need for plotting of approximant
structure. We shall use supercell 2 x 3 x 10
Close Diamond without saving
[On the screen: basic window of Jana2006]
Start “Plot structure”
Select “Draw approximant structure” and define supercell 2 x 3 x 10
(In case bonds obscure atoms too much, use Ctrl-A + Right click on a selected object + Edit
→ Atom design and change bond style from “thick” to “thin”)
Note starting letters (Na1, Na2 or Na3) of labels of sodium atoms which are not in the
middle of the plot
On the right sidebar switch to the table “Table of created atoms”
Order atoms by name and select the noted sodium atoms
“Build → Polyhedra → Add coordination polyhedra”
Select oxygen for “Ligand atom groups”
OK
With make sure that Rendering and Lightning is activated
The resulting plot will be:
OK; OK
[On the screen: Plot of three C-O distances as a function o t]
167 Jana2006 Cookbook, version October 2017
As we have already seen in Dist listing, C-O distances are almost constant over the whole
interval of t
Press “New/Edit”
Change “Parameter to draw” to “Bond valence sum”
For “p” change lower limit to 3.5 and upper limit to 4.2
Press “Define bond valence coefficients”
The dialogue shows the coefficients defined in Dist commands are still active
OK; OK
[On the screen: Plot of BVS of carbon as a function of t]
Using the mouse cursor we can find the fluctuation of BVS is between 3.94 and 4.04
Both Grapht and Dist confirm the C-O bonds are constant. Now we shall check planarity of
the CO3 group:
Start "Edit Atoms"
Select atoms O1, O2 and C1
Choose Action "Make symmetrically contiguous motif"
OK; YES to save changes
With this procedure we transformed positions of C1, O1 and O2 to be connected without
need of symmetry transformation
Start commands for Dist
In page “Basic” for “Listing form” check “With symmetry code”
In Page “Select atoms” make sure the central atoms is still carbon and the neighbor atoms
are oxygens
OK; YES+Start; YES to open the listing
“Goto → List of distances”
There are three distances: C1-O1, C1-O2 and C1-O2. For the third distance ,C1-O2, oxygen
O2 must be transformed from the position given in M40 by the symmetry operation x,-y, z,
which is encoded to the Jana symmetry code #s4 (i.e. the fourth symmetry operation). Other
168 Jana2006 Cookbook, version October 2017
distances do not need a transformation because "Symmetry contiguous motif" has been
created in the previous step.
Symmetry codes of may differ for various Superflip runs!
Close the listing of Dist
Start commands for Dist
In page “Basic” press the button “Planes”
For “Atoms defining the best plane” type “O1 O2 O2#s4”
For “Additional atoms” type C1
Add; OK; OK; YES+Start; YES to open listing
“Goto → List of best planes”
The listing shows that the distance of C1 from the plane of oxygens does not change with
the t coordinate which means the planarity of CO3 is not influenced by modulation.
At the next step we will test whether the atoms of CO3 group move during modulation in the
phase, i.e. like a rigid body. This can be visualized by “animated” plot in Contour.
Run refinement (at least zero cycles) to get input M80 for Fourier calculation
Start Fourier commands, select “Map type” “F(obs) – Fourier”
In the page “Scope” select “automatically”, “Use default map orientation”
OK; YES+Start; don’t include new atoms; don’t open listing
Start Contour; New Plot
Select “General Contour procedure”, “Use old map”, “Draw a general section”; OK
[On the screen: form for definition of the general section]
For the first atom type O1
For the second atom type O2
For the third atom type O2#s4
The three atoms define the plane of the section through the electron density map. The
dimensions of the plot are defined in Scope
For “Scope” type “4 4 0”
For “Interpolation step” type “0.02”
This means the plot will be 4x4 Å. Number of parallel sections calculated by the program is
1+s3/is, where “s3” is the width of the section and “is” is the interpolation step. Here it
means one section, because the third parameter is zero.
For “First point put to” type “1 1 0”
The first point of the plot is O1. The horizontal axis of the plot will be O1-O2. The last atom
O2#s4 will be above the horizontal axis.
Using Grapht, plot Na-O distances for Na1, Na2 and Na3. In the Grapht dialogue select
“Draw non-modulated curves” – they will appear in the plots like dashed lines indicating
distances in the basic structure.
Distances Na1-O
As it has been previously indicated by plot of the structure, Na1 which alternates with
carbon in a column along c, has different coordination comparing with octahedrally
coordinated Na2, Na3. Na1 has more neighbors with strongly varying distances which are all
longer comparing with Na2 and Na3. Good test if this complicated model is reasonable
provides plot of the bond valence sum.
Using Grapht plot graphs of bond valence sums for Na1, Na2 and Na3 (do not forget to
define coefficients for Na-O bond)
Graphs of BVS show that Na1, Na2 are oversaturated while Na3 is slightly undersaturated.
However, the fluctuation of BVS is in all cases small and similar. This means the modulation
of Na3 and its coordination is properly described with approaching and departing oxygen
atoms in balance.
Cr2P2O7
Single crystal data measured with Oxford Diffraction four-cycle diffractometer
Input files: directory CrPO_RT available at http://crysa.fzu.cz/jana/workshops/CrPO_RT.zip
This example is based on the Crysalis 171.38.43g. Crysalis is a program for controll and data
processing of Rigaku Oxford Diffraction (ROD) diffractometers. It has very friendly licencing
policy: everybody can install Crysalis for processing data acquired with a ROD instrument
(which is the case of Cr2P2O7).
2. Starting Crysalis
Start Crysalis
In a dialogue that follows locate the file CrPO_228_4.par
Crysalis opens the experiment.
Use to decrease the image size
[On the screen: the basic window of Crysalis]
In this wizard, each button has a default function started by pressing the centre of the
button and various optional functionalities, which are offered after pressing “>” part of the
button
Start “Peak hunting” with default options
About 6500 peaks have been found
Start “Unit cell finding” with default options
The program finds a unit cell 7.03, 8.42, 4.64, 90, 108.60, 90 and indexes about 30% of
peaks
From this projection, it is clear that the additional spots cannot be indexed by any
reasonable supercell.
Set projection along a*
The peaks between the equidistant planes are evidently some spurious peaks arising e.g.
from a small crystal attached on the main sample. We will now remove these peaks
In the right panel, switch to the page “Selection/Lattice-It”
Activate the radio button “Lattice-It”
[On the screen: Ewald Explorer with activated function “Lattice-it”]
At the bottom of the right panel, clear the checkbox of the group 2.
The spurious peaks now become invisible
At the right panel, switch to the page “Crystal”
Select vie along b*
Now we have a clear figure suitable for determination of the q-vector
Play with the arrows close to “Modify” button until the satellite position fit with the
projections
Using the button “Modify” at the right panel, increase the satellite intex to four.
The reason for four is that when we can see the third order satellites by an eye there might
be the forth order satellites not found by the default peak search
Check views along a* and c* to be sure that the satellites are described
Close the window of the Ewald explorer
OK
The program refines the unit cell and q-vector an indexes 93% of the peaks. This is quite
acceptable number because the refinement included also the previously suppressed spurious
peaks.
Press the “>” part of the button “Save information” and select “Save specific peak table
only”. Make backup copy of the peak table under some name which cannot be overwritten
by Crysalis, e.g. “incommensurate”
6. Data reduction
[On the screen: Basic window of Crysalis]
In the right panel, press the button “Start-Stop”
In the next window, select “Data reduction with options”
[On the screen: Step 1: Orientation matrix for data reduction]
Uncheck “Follow significant sample wobbling” (we know that the sample was not wobbling)
Note: In this case, no special parameters are necessary. For more complicated modulated
structures, e.g. with short modulation vector, there are two typical parameters one should
consider: (1) Disabling the option “Refine q-vectors(s)”; (2) “Overriding integration mask
size” available through the button “Edit special parameters”; (2) “Skiping model refinement”
available through the button “Edit special parameters” followed by Alt-E keyboard shortcut.
7. Output
1) The resulting data are placed to the directory struct\tmp. For Jana2006 we need the files
CrPO_228_4.cif_od and CrPO_228_4.hkl
2) In the left toolbar of Crysalis, there is a button Jana. The button triggers the following
actions: directory struct\ jana_CrPO_228_4 is created; files CrPO_228_4.cif_od and
CrPO_228_4.hkl are copied to this directory; Jana2006 is launched
The Rint is calculated together for observed and unobserved reflections. In order to see Rint
separately for observed and unobserved reflections, Jana2006 must be used
Cr2P2O7
Single crystal data measured with Oxford Diffraction four-cycle diffractometer
Input files: crpo.hkl, crop.cif_od
Frame scaling, absorption correction: done with software of the diffractometer
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction – CCD”; NEXT
In “Define basic input file“ select “Input from cif_od file“; OK
Leave all settings unchanged; NEXT
Leave all settings unchanged; NEXT
For absorption correction select “None or done before importing”; NEXT
FINISH
[On the screen: Data repository]
OK; YES to accept the data set
3. Symmetry Wizard
Use default tolerances for crystal system recognition
Select monoclinic Laue symmetry
Select C centering
Select space group C2/m
Select superspace group C2/m(a0g)0s
Accept the superspace group in the original cell
The symmetry codes will be different for different Superflip solutions, note your codes!
Start Contour and New plot
Select “General Contour procedure”
Select “Calculate new ones”; OK (if the option is not active, Contour has no old map and it
will calculate a new map by default)
In Fourier option, page Basic, select “Map type” = “F(obs) – Fourier”
In page Scope set calculation of the x1-x4 section with the central point Cr1
OK;
[On the screen: Contour window with plot of x1-x4 section through position of Cr1]
Press “Atoms edit” and add Cr1 to the list of atoms to be plotted
For making de Wolff sections of other atoms, we will use the automatic tool:
O1: section x1-x4 O2: section x1-x4 O3: section x2-x4, scope along x2 is 4 Å
For O3 we also display the modulation function related by symmetry operation O3#s3t0,1,3
(or your actual symmetry operation). To display such plot, use the following steps:
Press the Contour button “From list”
Select the section O3_y; OK
Press Atom edit and add O3#s3t0,1,3; OK
Run refinement
Refinement converges with Robs factors 21%, 10%, 19%, 54%, 65% and 129%
Run Contour and plot again the three sections:
Cr1 O2 O3
Start Edit atoms and add another one position modulation wave for all atoms
Repeat refinement
Refinement converges with Robs factors 11%, 4%, 13%, 20%, 57% and 97% for 68
parameters
Make backup copy “crpo_initial_crenel”
The resulting plot shows distances P1-O2, P1-O2#s4t0,1,0, P1-O3 and P1-O3#s3t0,1,3 as a
function of t and it is obvious that discontinuities occur for slightly different values of t. This
191 Jana2006 Cookbook, version October 2017
is however not possible because PO4 cannot be so much distorted. The discontinuous
changes of oxygen positions must therefore occur at the same t value.
The centre x40 of the crenel function of O3 must be fixed at 0.75 (or equivalent value
following from your structure solution), because this atom is in a special position; for
different center, there would be an overlap with symmetry related O3#s3t0,1,3 leading to a
disordered structure. The t-coordinate of this centre is t40[O3].
For O3 delta=0.5, therefore one of the two discontinuity points is t40[O3]+0.25
For O2 delta=1.0 and one of the two discontinuity points is t40[O2]-0.50
The following equation will align the two discontinuity points:
t40[O3]+0.25 = t40[O2]-0.50 => t40[O2]= t40[O3]+0.75 => t40[O2]= t40[O3]-0.25
(The symbol t40 means x40 transformed to t coordinate.) In the last step the equivalent
translation has been applied to get t40[O2] from the interval (0,1).
The discontinuity points for O2 and O3 now occur at exactly the same value of t
Make backup copy “crpo_aligned”
Quit Contour
Start Refinement option, page “Various”, button “Equations”
Define equation t40[P1] = t40[O3]
This bounds discontinuity points of P1 with O3
Run refinement
Refinement converges with Robs factors 10%, 5%, 10%, 14%, 56% and 96% (77 parameters),
i.e. better fit for 1st and 2nd order satellites
Make backup copy (rewrite crpo_aligned)
Oxygen O1 is only slightly modulated and does not need more waves.
Add another one position modulation wave for Cr1, P1, O2 and O3
Refine
Refinement converges with Robs factors 9%, 5%, 10%, 7%, 58% and 85%
Add the first ADP modulation wave for Cr1, P1, O2 and O3
Start Refinement options, change number of refinement cycles to 100 and set damping
factor 0.4
Refine
Refinement converges with Robs factors 7%, 4%, 6%, 12%, 38% and 57%
Add another ADP modulation wave for Cr1, P1, O2 and O3
Refine
194 Jana2006 Cookbook, version October 2017
Refinement converges with Robs factors 5%, 3%, 5%, 10%, 24% and 39%
Add another position modulation wave only for Cr1
Refine
Refinement converges with Robs factors 5%, 3%, 5%, 7%, 18% and 38%
Add another ADP modulation wave only for Cr1
Refine
Refinement converges with Robs factors 4%, 3%, 4%, 6%, 18% and 34% with 176 refined
parameters.
Make backup copy: crpo_simple
The result obtained here is not bad but a question remains if the poor fit of the 3rd and 4th
order satellites is caused by structure model or weak satellite intensity. For 3rd order
satellites there is still 50 % of reflections above the observability limit. GOF 3.42 also
suggests the model could be improved.
Start Contour and plot the following difference Fourier sections: x1-x4 for Cr1, x2-x4 for P1,
x1-x4 for O2 and x2-x4 for O3. For plotting use always the contour step 0.2.
For each of the atoms delta will be shortened and additional atoms with small
complementary delta will be placed to the difference maxima found in the Contour sections.
The resulting structure model must keep formula unchanged. The delta of all disordered
atoms will be kept the same.
The atoms Cr1 and O3 have to be duplicated and the atoms O2 triplicated. Then the atomic
parameters of the original and split atoms must be modified to account for the disordered
model as found from difference Fourier maps.
Start “Edit atoms”
Select Cr1, O2 and O3; Press “Action” → Select “Split atom”;
Choose 0 for relative occupancy of split atom;
Check the manually defined shift of atoms equal 0 0 0;
Uncheck “Reset ADP to isotropic”;
Uncheck “Generate restrict command for refine;
Press “Apply->Go to next”;
Repeat this procedure for O2 and O3;
Select the atoms O2’ and repeat the split procedure once more to create O2’’;
OK; Yes to the question about rewriting;
The new positions will describe disorder at the discontinuity area. For Cr1’ and O3’ two
positions are generated by symmetry; for O2 we have introduced both positions.
We want to assign some non-zero delta for additional atoms and view modulation functions
but with the original electron density
Start “Edit atoms” and change delta of Cr1 from 1 to 0.9
In refinement we have still zero cycles
Run refinement
This recalculates delta based on equations and creates new (unwanted) input for Fourier
Start “File → Structure → Save As” and create copy “tmp”
Do not continue with the new structure
Change delta of Cr1 back to 1
Run refinement (with zero cycles)
This creates good input for Fourier
Start “File → Structure → Copy in” and copy to the current structure only M40 from the
previously created backup copy “tmp”
Now we have old input for Fourier and new input for plotting modulation functions
Start Contour
For Cr1 plot x1-x4 difference section and visualize Cr1, Cr1’, Cr1’#s2t0,0,2
For O2 plot x1-x4 difference section and visualize O2, O2’ and O2’’
For O3 plot x3-x4 difference section and visualize O3, O3’ and O3’#s2t0,1,3
13. Conclusions
The disordered model does not provide proper geometry of disordered atoms. However, we
can draw an approximant structure in large area and try to determine in which places of the
real space the structure becomes disordered.
In view along c we can distinguish columns without disorder and columns with disorder
Press “New/Edit”
Change Minimum and Maximum for t: from 0.21 to 0.7
OK
From refinement of disordered structure we have found delta of the disordered area 0.105.
In order to get reliable distances from the simple model we limit the interval where distances
are plotted to reliable interval.
Press New/Edit
Change Minimum and Maximum for t: from 0.21+0.05=0.26 to 0.7-0.05=0.65
OK
Ba2TiSi2O8
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction – CCD”; NEXT
In “Define basic input file“ select “Input from cif_od file“; OK
Change all q-vector components to 0.302 keeping the original sign; NEXT
Leave all settings unchanged; NEXT
For absorption correction select “None or done before importing”; NEXT
FINISH
[On the screen: Data repository]
OK; YES to accept the data set
3. Symmetry Wizard
Use default tolerances for crystal system recognition
Select tetragonal Laue symmetry 4/mmm
Press Show/modify X centering
[On the screen: select individual centering vectors]
Select centering vectors: (0 0 0 0 0) and (0 0 1/2 1/2 1/2); OK; Next
Select space group X4/mbm and note the non-centrosymmetric possibilities
Accept the superspace group in the original cell
For five-dimensional structures the wizard does not offer selection of superspace group. It
will be done later manually.
The figures indicate a disorder of SiO4 groups by the mirror perpendicular to the c axis. The
disorder can be described alternatively by the two fold axis perpendicular to c and a+c.
The disorder in the average structure can be produced either by the fact the structure is not
centrosymmetric or by modulation. First we shall test non-centrosymmetric description.
The highest subgroups not giving a split model for SiO4 group are P-421m and P4bm.
Now we have to check if basic positions (x,y,z) of atoms Si1, O1, O2 and O4 belong to a
unique tetrahedron.
Return to the structure “ave-centro”; Set distance commands: in the “Basic” define
calculation of angles and the listing form “with symmetry codes”.
The blue lines indicate the bond distances and red lines bond angles of the basic SiO4 group.
From the symmetric codes we can see that only the atom O2 should be transformed by the
symmetry operation #s10 i.e. by x,y,-z
Start “Edit atoms”: select the atom O2 and transform it by x,y,-z
Before testing of non-centrosymmetric space groups we have to reset the occupancies to
their the values corresponding the SiO4 split. This means that their full values are to be
reduced by 1/2. The atom O3 is coordinating Ti and it can also be shown that it has to
represent a split position.
Start “Edit atoms” set “ai” for O1, O2, O3 and O4 to 0.25, 0.5, 0.25 and 1, respectively.
9. Superspace symmetry
In the previous paragraph the program asked for completion of the translation part of
superspace symmetry operators. For the present case two symmetry operators should be
completed:
and
Each translation part has two elements which can be either 0 or 1/2. Thus we have sixteen
combinations: 0 0 + 0 0, 0 0 + 1/2 0, 0 0 + 0 1/2, 0 0 + 1/2 1/2 … However, the translation
part can be combined with the X centering (0 0 1/2 1/2 1/2) which makes equivalent for
instance the pairs 0 0 + 1/2 1/2 and 1/2 1/2 + 0 0. Moreover, some combinations are
forbidden because they lead to inconsistent translation (0 0 0 1/2 1/2).
Finally we have only two independent possibilities:
0 0 + 1/2 1/2
00+00
208 Jana2006 Cookbook, version October 2017
Start EditM50
Go to page “Symmetry”
Focus the textbox “Space group” and press TAB. The program offers the first symmetry
operator needing to complete.
At the first attempt try combination 0 0 for the first operator and 1/2 0 for the second one.
The program will not allow it.
Then try the combination 0 0 + 0 0. Close EditM50 and re-create refinement reflection file.
From the screen output of the symmetry wizard note number of rejected reflections: 54372
rejected reflections from which 455 is observed with the largest I/sigma(I) around 15
Return to EditM50
TAB in the textbox “Space group”
For translation parts use combination 0 0 + 1/2 1/2
Re-create refinement reflection file
From the screen output of the symmetry wizard note number of rejected reflections: 55353
rejected reflections from which 463 is observed with the largest I/sigma(I) around 15
This means: the translation part 1/2 1/2 of the second operator leads to 1000 more forbiden
reflections but the number and sigma level of the observed forbidden reflections remains
almost the same like with the translation part 0 0. This confirms the presence of the
translation part 1/2 1/2 for the second symmetry operator.
This will plot section x1-x4 but the shape of atomic domain will be function of x5
OK to calculate Fourier map
In Contour press “Atoms edit” and add Ba1
The difference maxima around Ba1 are very strong indicating ADP modulation
Refinement of fresnoite tends to be trapped in a false minimum because of very good fit of
main reflections and relatively weak satellites. Inversion of the position of barium helps to
leave the false minimum.
Ca2MgSi2O7
Single crystal data measured with Oxford Diffraction four-cycle diffractometer
Input files: directory Melilite available at http://crysa.fzu.cz/jana/workshops/Melilite.zip
This example is based on the Crysalis 171.38.43. Crysalis is a program for controll and data
processing of Agilent diffractometers. It has very friendly licencing policy: everybody can install
Crysalis for processing data acquired with an Agilent instrument (which is also the case of
Melilite).
2. Starting Crysalis
Start Crysalis
In a dialogue that follows locate the file run.par
Crysalis opens the experiment.
Use to decrease the image size
The projection shows main reflections with radial distribution of spurious peaks. These peaks
are tails of very strong reflection and their radial districution allows quick conclusion that
these are not satellites.
Majority of satellite peaks can be described by two q-vectors: q1 = (α, α,0), q2 = (-α, α,0). The
214 Jana2006 Cookbook, version October 2017
positive direction of a* and b* is indicated in the shrinked Ewald view by very small labels.
The peaks indicated by white circles are λ/2 reflections. The experiment was done with
classical Mo x-ray tube and graphite monochromator and with this equipment such
reflections may exist between very strong regular reflections.
The most important are weak reflections indicated by red circles. These reflection are
described by combination of q1 and q2, which proves that our structure is (3+2)-dimensional.
Without the combinations there would also be possibility of a (4+1)-dimensional crystal
twinned by 90° rotation around c*.
The projections also indicate that the combined satellites may be overlapped due to short
distance between them
In the next step, we will visualize the reciprocal plane 0,k,l and 1,k,l to see better the
combined satellites
4. Unwarping
Close the Ewald explorer
In Lattice wizard, press “Unwarping – precession images”
NEXT; NEXT
[On the screen: Step 3: Creation of a layer list”]
Press “Generate layers”
Type 1 for “max order”, 1.2 for “resolution”
Delete all layers except 1_k_l
NEXT
[On the screen: Step 4: Use background subtraction”]
Check “Use background subtraction”
NEXT
[On the screen: Step 5: Data corrections”]
Press “Edit user scale”
Define a user scale 40 (to see weak features)
Check “Don’t apply final convolution filter” (to see better separation of spots)
FINISH
Check the radio button “Layer” to see the progress of the unwarp procedure.
Wait until the layer is finished
Open the menu in the lower right cornes of Crysalis and select “Unwarp reciprocal lattice
grid”
The weak spots indicated by yellow rectangle are the λ/2 reflections. The ones indicated by
red circles are the combined satellites. Again, the possibility of overlaps between the
combined satellites is evident.
The resulting image projects all peaks positions into one unit cell. We can recognize q1 and
q2 satellites, combined satellites (the double-points in the middle of the unit cell axes) and
the λ/2 reflections in the middle of the plot.
Use to determine satellites coordinates
Use Refine cell to obtain refined q-vectors
Quit Jana2006
7. Data reduction
In the right panel of Crysalis, press “Start/Stop” and select “Data reduction with options”
218 Jana2006 Cookbook, version October 2017
[On the screen: Step1: Orientation matrix for data reduction]
Select “Other (reduction list)”
OK to confirm a message about non-existing reduction list
Press “Generate”
[On the screen: Proffit reduction list generation]
This tool generates list of reflections, which will be used for the data reduction.
Crysalis could also read an existing list of reflections (the “Load” button). This list can be
produced by Jana2006 and it is very useful for cases when the list cannot be generated by
Crysalis, e.g. for composite structures or for structures with X centering vectors. In this case,
however, we can generate the list by Crysalis.
In the case of Melilite, the combined (cross) satellites need a special treatment. The previous
attempts showed that these satellites are partially overlapped and that the best way how to
process them is making the data reduction with a large integration mask. With such mask,
the pairs of close combined satellites are integrated as a one spot and this can be taken into
the account in the Jana2006 refinement. However, this large mask is not suitable for main
reflections and for the q1/q2 satellites. For this reason, the combined satellites must be
processed separately.
The Rint of satellites is poor but we should consider that the program does not show
separately the observed satellites. Rint of observed satellites will be shown in the Reflection
report of Jana2006.
Select the page Graphs
The graph of frame scaling is almost flat, indicating sufficient absorption correction
In the next step, we will process the combined satellites. We want to ignore the real position
of the cross satellites, use instead the position between the pairs of close satellites and large
integration mask to cover the complete pairs. For this, we have to disable refinement of
instrument parameters and profile analysis.
Quit the tool “Inspect data reduction results”
Start Lattice wizard
Use the “>” button at “Incommensurate/qusicrystal” and set q-vectors
q1 = ¼ ¼ 0
q2 = -¼ ¼ 0
CANCEL; CANCEL
In the main window of Crysalis - the left toolbar, press the button CMD
Type command “dc extinct hklmn 0 0 0 1 1 2 1”; ENTER
This defines that reflections with m+n odd will be skipped
In the right panel of Crysalis, press “Start/Stop” and select “Data reduction with options”
[On the screen: Step1: Orientation matrix for data reduction]
Select “Other (reduction list)”
Press “Generate”
In the line for q1, press “Edit mmax” and define mmax=1
In the line for q1, press “Edit mmax” and define mmax=1
Check “Use cross q-vector reflections”
Press “Edit max index” for the cross q-vector reflections and define mmax=1
At the bottom of the “Generate hkl list” window we can see that our custom extinction rule
is active. It can be checked by “Show rules” button
OK;
Crysalis displays warning about missing “nada/um q-vector refinement” but this is OK
because or q-vetors should not be refined
Out of Crysalis, go to the experimental directory and open in a text editor the file
run.dcredhklmn
After finishing the data reduction, start “Inspect data reduction” by pressing
The results for combined satellites are poor but this might be caused by their weakness. We
cannot make any conclusion without refinement of a structure model.
The example 5.5.2. where structure of melilite is solved uses different input files: 4d data for
q1, 4d data for q2 and 5d data for combined satellites. The reason is that we need to have
the user-independent results. Data created in this example are tested at the end of the
Example 5.5.2.
Ca2MgSi2O7
Revised: 2 December 2016
Single crystal data measured with Oxford Diffraction four-cycle diffractometer
The structure has two q-vectors: q1=(α,α,0) and q2=(-α,α,0). Satellite reflections h,k,l,m,0 and
h,k,l,0,n are always well separated in the CCD frames. Combined satellites hklmn, m≠0, n≠0, are
sometimes separated and sometimes overlapped, depending on their geometry in particular
CDD frame (i.e. angle ψ). The combined satellites needed different data processing strategy
using large integration mask because we want to treat them like “always fully overlapped“.
Schematic view of the diffraction pattern in melilite. Main reflections are indicated by large
dark circles, small circles symbolize first-order satellites (grey for m x n = 0 and white for m x n =
±1 satellites, respectively).
Input files:
q1.hkl, q1.cif_od (main reflections + satellites described by q1)
q2.hkl, q2.cif_od (main reflections + satellites described by q2)
q1q2.hkl, q1q2.cif_od (main reflections + combined overlapped satellites)
Frame scaling, absorption correction: done with software of the diffractometer
2. Import Wizard
Select “Known diffractometer formats”; NEXT
8. Superspace symmetry
We will need to complete the translation part of superspace symmetry operators.
“EditM50 → Symmetry”: focus the textbox with the space group symbol and press TAB.
The program shows two symmetry operators where the translation part in x4 and x5 needs
to be completed:
and
Each translation part has two elements which can be either 0 or 1/2. Thus we have sixteen
combinations: 0 0 + 0 0, 0 0 + 1/2 0, 0 0 + 0 1/2, 0 0 + 1/2 1/2 … However, some of them are
forbidden because they lead to a non-consistent translation vector, which cannot ne used for
generation of the superspace group. Finally, we have only eight combinations. In the
following table we show number of extinct reflections (observed/unobserved) which
Jana2006 reports for each of these possibilities:
Only four possibilities lead to additional extinctions. They are all equivalent as we can judge
from the exactly same number of extinct reflections as well as from the fact that EditM50
does not offer creation of new refinement reflection file when we switch between them.
Another four possibilities without additional extinctions are equivalent, too.
231 Jana2006 Cookbook, version October 2017
Thus we have two possibilities for the superspace group: with and without additional
extinctions. The choice is simple because the additional extinctions mostly concern
unobserved reflections and therefore the symmetry leading to additional extinctions is
confirmed by the experiment.
In order to be consistent with the published structure we will choose 0 0 for the first
symmetry operator and 1/2 1/2 for the second one.
Select 0 0 and 1/2 1/2 for the first and the second symmetry operator, respectively
Close EditM50, save the result
Create new refinement reflection file
Press “Movie” to see the changes of this figure with x5. ESC to stop.
Run Grapht and investigate distances, bond valence sums etc. All obtained graphs will
depend on “u” coordinate, which can be changed by “Movie” button.
Use the old reference parameters - this causes that 0.25 atrufucial components of q-vectors
will not be used
Refinement converges with the following Robs values for all reflections, main reflections, 1st
and 2nd order satellites: 4.4%, 2.3%, 4.1%, 7.9% for 238 parameters
Single crystal data measured with a four-circle diffractometer and CCD detector.
Input files: iodomol.p4p, iodomol.ram
absorption correction: done
2. Import Wizard
Select “known diffractometer formats”; NEXT
Select “Bruker CCD”;
Click the button “Browse” and double click “Iodomol.ram”; NEXT
[On the screen: complete/correct experimental parameters]
Leave all settings default; NEXT;
[On the screen: Define the reference cell/split by twinning]
Set target dimension to 4; NEXT; OK TO confirm the import of 52313 reflections
In the text window, leave all settings default; NEXT; Finish
3. Data Repository
The line in data repository window corresponds to one data set.
OK; YES to accept the data set
[On the screen: “Information about starting symmetry wizard”]
NEXT
4. Symmetry Wizard
[On the screen: “Tolerances for crystal system recognition”.]
Symmetry wizard can be started separately by “File → Reflection file → Make space group
test”.
Leave all settings of the wizard default; NEXT
[On the screen: “Select supercell”.]
Confirm “continue with the basic unit cell”; NEXT
[On the screen: Select Laue symmetry – the Rint values are quite elevated]
Confirm monoclinic Laue symmetry; NEXT
[On the screen: Select cell centering]
Confirm primitive unit cell; NEXT
[On the screen: Select space group]
Confirm P21/m; NEXT
[On the screen: Select superspace group]
Select P21/m(a0g)0s – out of 1038 potentially absent reflection none is observed, this space
237 Jana2006 Cookbook, version October 2017
group is thus the most likely.
[On the screen: Final step of the superspace test]
Accept the superspace group transformed into the original cell: P21/m(a0g)0s; FINISH
Close the viewer of the listing, wait until the structure solution is interpreted
Click “Draw structure
Choose the radio button “Draw average structure”
If you have followed this tutorial exactly, you should see a picture identical to the following
(to obtain the picture in Diamond, click the icon “Get molecules directly” - )
If the structure looks slightly different (e.g. shifted by ½ along a or b), you may still continue
with the tutorial, but the details (atom names, atom positions) may be different.
Alternatively you may rerun Superflip until you obtain the solution in the correct setting.
Close the structure viewer and exit the solution wizard by clicking “Accept last solution”
Important: Make sure you copy the job to a new job name to preserve the original and to
have the chance to return to this point.
Go to “File->Structure->Save as
In the lower right text box type “Iodomol-ref”; OK
Yes to continue with the new structure
Right-click the icon Refine
Set the number of refinement cycles to 100:
9. Refinement of ADP
So far the displacement parameters are refined as isotropic. We will now set them for
anisotropic refinement:
The R-value on the satellite reflections decreased by almost 2%, the modulation of ADPs
seems significant.
Save the structure as a new structure under the name iodomol-mol and continue with the
new structure;
Draw the average structure and check the geometry. Plot the structure with thermal
ellipsoids. The molecule looks like this:
Ring 2
Ring 1
The result is not a disaster, but the displacement parameters are uneven. The data quality
may not be sufficient to refine all atoms independently. We can define the three phenyl rings
as three positions of the same molecule.
Note the atom names for the free phenyl rings. Marking the phenyl rings as in the image
above, you should get the following lists:
Ring 1: C3, C6, C9, C13, C18, C12
Ring 2: C5, C15, C20, C22, C10, C8
Ring 3: C16, C17, C19, C21, C14, C11
Close the structure viewer.
Go to Parameter -> Molecules -> New molecule;
Type “phenyl” in the text box “Name of the molecule”; Next
Select atoms C3, C6, C9, C13, C18, C12 from the list; OK;
Choose gravity center as the molecular reference point; Next;
Note that the occupancy of this position is set to 0.5, because Ring 3 is located on the mirror
plane.
Click “Apply + Next position”;
Add the next position by identifying the model atoms with the atoms of the next phenyl
ring. Fill the table according to the following picture:
Check again that all six atoms of the model molecule coincide with the actual atoms in the
atomic part of the model.
Click “Apply + End”
Now we defined all three phenyl rings as instances of a single molecule. The positions of
atoms in the molecule will still be refined, but one molecule will be placed in three positions
in the crystal. We also have to redefine the modulation parameters – remove the
OK; Yes;
Note that the R-values are only marginally higher than in the previous refinement, but the
number of parameters decreased by a factor of two, from 344 to 176.
For future convenience we will now rename the atoms in the molecule sequentially as
C1…C6:
Note: the following part includes some slightly more advanced techniques, and it turns
out the alternative model is not better than the previous one in this particular case. It is a
very useful exercise though. You are encouraged to follow it, but if you are short of time
or need to focus on something else, you may skip this part and continue with the next
section (12. – viewing the structure).
If you check carefully the de Wolff’s sections of some of the carbon atoms, you will see that
it seems more appropriate to describe the modulation not as a harmonic oscillation of the
atom, but as a discontinuous switching between two positions. A good example is the x
coordinate of atom C2b or C3b.
In the following section we will set up this modulation for the phenyl rings to check, if this
model may be more appropriate than the harmonic model.
Setting up the crenel occupational modulation is in general more difficult than the harmonic
modulation, but often the model fits better and makes more sense.
General idea: Rings 2 and 3 need to be split to two copies, each will be present only in one
half of the complete domain along x4 and the domains will be complementary along x4.
Ring 1 lies on a mirror plane and it does not need to be split. If its occupation domain is
First save the structure as iodomol-crenel and continue with the new structure
Check the de Wolff’s sections for the atoms of the ring 3 (position 3 of the phenyl
molecule):
Go through the de Wolff’s sections of atoms C1a, C2a, C3a, C4a, C5a, C6a and try to make a
conclusion about the behavior of the ring as it is modulated. Note the x4 coordinate of the
largest deviations of the atoms from the average position.
Conclusion: the largest modulation is along the y (x2) coordinate. The smallest modulation is
on the C1c atom, which is bonded to the phosphorus. C4c has also a relatively small
modulation amplitude compared to the rest of the atoms. This indicates that the phenyl ring
is pivoting around the axis C1c-C4c. This is also confirmed by the shape of the modulation
functions of C2c, C3c, C5c and C6c. The first two have very similar modulation functions and
it is just the opposite of the latter two:
The occupancy domain should have a length of 0.5 and the center approximately 0.4. To
generate a perfect model, we should also rotate the molecule from its average position.
However, this is difficult to do, as the rotation angles of the molecule are, in general, not
simply related to the crystallographic axes. We will therefore leave the molecule as it is and
let its position refine.
Close Contour
Go to “Parameters -> Molecules -> Edit”
Set number of positional modulation waves for the Position 1 to 0
Check the check box “Use: crenel”
Set the parameter aimol from 0.5 to 1
The previous step is done because now,due to the shortening of the domain, the molecule
does not sit on a special position anymore, and the occupancy need not be reduced.
Click “Edit modulation parameters -> Occupancy”
Set delta to 0.5 and x40 to 0.4
Check the checkbox next to x40 to allow its refinement
OK; OK; Yes to confirm changes
Go back to the de Wolff’s sections and check that the occupation domain of all atoms of the
Ring 1 reduced to one half of the full length.
Make the analysis of the modulation functions of the atoms from Ring 2 and Ring 3. Can you
estimate the modulation characteristics of these rings similarly to Ring 1?
247 Jana2006 Cookbook, version October 2017
Conclusion for Ring 2: largest modulation along x2, all atoms have a similar modulation
phase (maxima and minima at the similar x4). Minima for x4 between -0.05 and 0.2,
average around 0.07. Modulation amplitude along x1 about 0.01 in fractional coordinates.
Conclusion – the main feature of the modulation of Ring 2 is the shifting as a rigid body
along x2
Conclusion for Ring 3: the largest modulation is along x1, all atoms have approximately the
same amplitude of the function along x1 and also the minimum and maximum of the
modulation function of all atoms is located at about the same x4 – minimum at ~0.15,
maximum at ~0.65.The amplitude of the modulation along x1 is about 0.04 in fractional
coordinates. Conclusion – the main feature of the modulation of Ring 3 is the shifting as a
rigid body along x1.
Rings 2 and 3 are not on a special position, and therefore we first need to split each of them
to two independent molecules
Repeat the procedure for positions #3, #4 and #5 with these settings:
#3: x40=0.15, delta 0.5, xtrans -0.345 (i.e. shift by -0.04)
#4: x40=0.57, delta=0.5, ytrans=0.137
#5: x40=0.65, delta 0.5, xtrans -0.265
Do not to forget to reset the number of positional modulations to 0 for all positions!
OK; Yes to confirm changes
Check the de Wolff’s sections again to check that the new parameters indeed generated
crenel-like modulations and shifted the atoms to approximately the right positions. You
should see images similar to this (example for atom C2e, coordinate x1):
The R-value on the satellites is higher than in the previous refinement, but we have not yet
introduced positional modulation of the newly generated positions.
The result is slightly better than the previous refinements with harmonic modulation.
However, the number of parameters is much larger. Checking the de Wolff’s sections shows
that the refined positional modulations are not very realistic, the amplitudes are large and
the shapes tend to converge to the harmonic shape. We can thus conclude that in this
particular case the harmonic modulation is probably a better model.
We have so far not looked at the modulated structure in detail. We will now try to interpret
the modulation.
Traditionally the modulation functions are interpreted with the help of plots of distances,
angles, bond-valence sums etc. as a function of the modulation phase. This may be difficult
for molecular crystals, where the modulation is usually quite complex. It is more useful to
view the structure either in a supercell or as an animation of the modulation.
Open the structure iodomol-mol with the harmonic modulation of the molecules.
Double click the icon “Plot structure”
Select “Draw approximant structure”
Set the range of x from 0 to 5 and for y from 0 to 6
5x6 supercell closely corresponds to the reciprocal values of the modulation vector
components of 0.19 and 0.15. We should see a good representation of various phases of the
modulation.
Click “Draw+continue”
In Diamond fill the unit cell and draw bonds (icons and ).
If Diamond is too slow, you can speed it up by switching off the rendering – icon
The best view of the structure is along a. You can see how the rows of the hydrogen iodide
molecule (only iodine present) form a wave, followed by the oxygen, to which the hydrogen
iodide molecule forms a hydrogen bond. Therefore the distance iodine-oxygen remains quite
constant. Also the large changes of the torsion angle of Ring 3 are visible throughout the
supercell. However, it is still difficult to make sense of the interactions between the molecule
and the overall picture of the modulation.
The program MoleCoolQt can read and display modulated structure refined by Jana.
Unfortunately, so far it cannot read modulated molecules. To see the movie, we first need to
transform the molecular model to purely atomic model.
First save the structure as iodomol-display and continue with the new structure
Open “Edit atoms”
Select all carbons and hydrogens from the phenyl molecule (in blue)
Click “Action -> Atoms from molecule to atomic part
OK; Yes
Start the program MoleCoolQt
Go to “File -> Open…”, locate the file iodomol-display.m50 and open it
You may also drag an icon of M50 displayed in your file manager on the icon of MoleCoolQt
Keep the selection “Load as modulated structure”; Open
Switch off the atom labels – icon in the top toolbar
You should obtain this image:
Modulation movie starts. It shows the molecule as a function of modulation phase. Keep in
mind that this is NOT a movie of the actual molecular motion, although it is tempting to
interpret it like this! These are snapshots of the unit cell configurations at various modulated
phases, which are distributed across the whole crystal.
You may use the options in the menu “Pack” to fill larger part of the structure to see the
intermolecular interactions and their changes with the modulation phase.
Single crystal data measured with KUMA four-circle diffractometer and point detector.
Reference: Petříček, Dušek and Černák (2005), Acta Cryst. B61, 280-286
Input files: cacni.hkl
absorption correction: done
2. Import Wizard
Select “reflection file corrected for LP and abrorption”; NEXT
Select “General file on I”;
In the text box “Input format” Type “*” (single star, without quotes); NEXT
[On the screen: complete/correct experimental parameters]
Fill unit cell parameters: 8.4867 15.9511 7.6038 90 90.48 90
Type the modulation vector on the “1st modulation vector” box: 0.4365 0 0.1307; NEXT;
In the text window, leave all settings default; NEXT; OK; Finish
3. Data Repository
The line in data repository window corresponds to one data set.
OK; YES to accept the data set
[On the screen: “Information about starting symmetry wizard”]
NEXT
4. Symmetry Wizard
[On the screen: “Tolerances for crystal system recognition”.]
Symmetry wizard can be started separately by “File → Reflection file → Make space group
test”.
Leave all settings of the wizard default; NEXT
[On the screen: “Select supercell”.]
Confirm “continue with the basic unit cell”; NEXT
[On the screen: Select Laue symmetry]
Confirm monoclinic Laue symmetry; NEXT
[On the screen: Select cell centering]
The data were measured with C-centering assumed. Confirm C centered unit cell; NEXT
[On the screen: Select space group]
Confirm C2/c; NEXT
[On the screen: Select superspace group]
Confirm C2/c(a0g)0s suggested by Jana based on systematic absences
253 Jana2006 Cookbook, version October 2017
[On the screen: Final step of the superspace test]
Accept the superspace group transformed into the original cell: C2/c(a0g)0s; FINISH
If the structure looks slightly different (e.g. shifted by ½ along a or b), you may still continue
with the tutorial, but the details (atom names, atom positions) may be different.
Alternatively you may rerun Superflip until you obtain the solution in the correct setting.
Obviously, Jana placed a single carbon atom in place of every cyano group. There is also a
split Cd atom of which only one position is marked. However, this picture is just one of many
unit cells in the modulated structure. To see the complete picture, we will run a movie of all
configurations:
Click the button “Movie”, set the start and end of x3 to 0.25, set delay time to 0.5 s and run
the movie by clicking OK. The movie can be interrupted by pressing Esc. Watch the movie
carefully and try to understand, what happens in the structure.
- Atoms C2 and C3 should be deleted and replaced by properly placed C and N atoms of
the cyano group
- Atom Cd is modelled by a continuous modulation function, but in fact it is jumping
between two positions and it is better modelled with a split position.
We shall now make these modifications by first deleting atoms Cd1, C2 and C3 and then
placing them in the correct positions:
In the form that opens, click “Include the selected point as a new atoms in M40”
Fill the next form according to this picture (i.e. atoms Cd1, atomic type Cd):
N3
C3
C2
Cd2
N2
Cd1
Please make sure to follow the scheme, so that the atomic names further in the tutorial
correspond to your structure model.
If you now click on the Cd atoms and let Jana show the position, you will probably see
something like this:
The Cd atom sits on the special position (two-fold axis) and the clicked position was slightly
away. We will correct this offset later.
Click the button Quit and confirm adding new atoms in the structure model.
Now we will similarly check the other section of the structure, the plane b-c at x=0.5:
Run movie (button movie, set x1 range to 0.5-0.5, delay time to 0.5) and investigate the
behavior of the modulation in this plane. Do we need to modify the structure model?
Your conclusion should be that the model is essentially OK. Atom C4 is modulated between two
positions and is always present close to the neighboring Cd atoms. The Cd atom oscillates
between four-fold and two-fold coordination (in this plane!). Atom C4 should be nitrogen
(ammonium). Atom C1 should also be a nitrogen atom. We can now make the final corrections
to the initial model:
The name changes to the first free number for nitrogen, N1 in this case
Now we have checked and updated the basic structure – atomic types and positions. Now we
can check the modulation functions of individual atoms. Jana tries to assign harmonic
modulation functions to each atom after the solution.
We will now check the modulation function for all atoms and modify them where necessary.
Double click the icon “Contour”
Click the button “New plot”
Select the radio button “Show de Wolff’s sections”
Set the scope of the sections: primary scope to 3, secondary to 0.5
Click the button “Define atoms for de Wolff’s sections
In the text box “Atoms” type * and click the button “Add”; OK
OK
The program calculates the Fourier maps for modulation functions of all atoms (x, y and z
direction for each atom). It may take some time.
Go through the de Wolff’s sections for all atoms (using the buttons “W-“ and “W+” or by
selecting from the list and make notes about the type of modulation for each atom and for
special function note their centers (denoted x40, it is the x4 coordinate of the center of the
domain) and widths (length of the domain along x4, denoted Δ).
Quit Contour
Double click “Edit atoms”
Double click atom Ni1
Click the check box “Use crenel”
In the form that opens set delta to 1, x40 to 0.75; OK
260 Jana2006 Cookbook, version October 2017
Select the radio button “Legendre polynomials in crenel interval”
Note: we defined the function as crenel, although it looks like a saw tooth. Legendre
polynomials contain a linearly sloping term, and a crenel function can turn to a saw tooth
during the refinement.
Repeat the procedure for N1, Cd1, Cd2, set the parameters according to the list above
Right click atom N2, select “Split atomic position”
Set the atom name to N2b; Uncheck the option “Generate restrict commands for REFINE”;
Apply
Double click atom N2 and set its Crenel function according to the list above. On the “Edit”
tab set its x-coordinate to 0.21. Repeat for N2b, set parameters according to the list above.
Repeat the whole splitting procedure, setting of Crenels and change of x coordinates for the
atom C2.
As a last step we will release the width of the crenel domain of Cd1 for later refinement.
Only this parameter needs to be explicitly released, other domain widths and centers will be
restricted before the refinement – see below.
Double click the atom Cd1
On Card “Edit” click “Occupancy” next to the text “Edit modulation parameters”
Check the check box next to “delta”
OK; OK; OK; Yes to quit “Edit Atoms”
The difficult part of the work is over. We have set up the structure model, we have
understood the structure in quite a detail and now we can start the refinement from a good,
elaborated structure model.
Quit Contour and return to the main window.
Important: what you have done up to now is very precious. Make sure you copy the job to
a new job name to preserve the original and to have the chance to return to this point.
Go to “File->Structure->Save as
In the lower tight text box type “CaCNi-ref”; OK
Yes to continue with the new structure
Plot the average structure. Use “Draw and Continue”, keep the structure view open and
return to Jana2006
Right-click the icon Refine
Set the number of refinement cycles to 100 and damping factor to 0.5:
The last step we need to make before the refinement is setting up restrictions on the
modulation parameters. The logic of the structure requires that the groups Cd1+N2b+C2b
have the same “domain of existence”, i.e. are present together in the structure. The group
Cd2+N2+C2 should also be present together and the two groups are mutually exclusive.
Furthermore, no vacancy is expected, i.e. either Cd1 or Cd2 is always present. Thus, the
centers of the crenel domains within each triplet should be equal and between the triplet
shifted by 0.5. The length of the domains is equal within each triplet, and the domain
lengths of one triplet are complementary to 1 with the domain lengths of the second triplet.
We need to define these relationships by means of equations between the atomic
parameters:
Note: There is no equation on t40 of Cd2, because this parameter is fixed by symmetry.
262 Jana2006 Cookbook, version October 2017
OK; OK; Yes + Start
The refinement converges to the following R-values:
It is a promising start. The high R-values for the 2nd and 3rd-order satellites are not
surprising, we have not added higher-order positional modulation yet.
Satellites up to third order are observed. We may add up to three harmonic waves, but first
we will add only one waves for each light atom and two waves for the heavy atoms.
Open “Edit atoms”
Double click Ni1
Increase the number of positional waves to 2
Repeat for other atoms, use 2 waves also for Cd1 and Cd2 and 1 wave for all other atoms.
Close Edit atoms and run refinement again.
The refinement converges to the following R-values:
Recalculate de Wolff’s sections and check, if the newly refined modulation functions are still
follow the traces in the Fourier maps.
The sections look OK, the modulation functions are appropriately describing the electron
density. The only exception may be the Cd1 and Cd2 atoms, where the functions are too
curved towards the end of the occupancy domain:
Checking the shapes of the modulation waves is an important check to be made after every
significant change of the structure model to avoid divergence of the model to unrealistic
parameters.
So far the displacement parameters are refined as isotropic. We will now set them for
anisotropic refinement:
Main reflections and first-order satellites now have nice R-factors, but the fit to the higher
satellite orders is not yet good. We need to add higher harmonic positional waves to
properly describe these parameters. There is not universal recipe on the number of harmonic
waves to be set except that their number should almost never be larger than the number of
satellite orders. At the same time, the refined waves should not be too wavy and
unreasonable.
A little bit of testing (try!) shows that an acceptable setting which lowers the R-value and
does not lead to unphysical modulation function is this:
Ni1: 3 waves
Cd1, Cd2, C3, N3, N4: 2 waves
N1, N2, N2b, C2, C2b: 1 wave
The fit is already quite good, but still not perfect. Further testing and checking of the
difference Fourier maps shows that the modulation of ADPs of Cd atoms is required to
improve the fit to the high-order satellites:
2. Import Wizard
Select “Various CW format”, NEXT
CW stands for constant wavelength
Type “File name” KSm.dat
Select “Free format of 2th, I, [sig(I)]” and “Another/unknown method”; NEXT
Fill cell parameters 5.5304 5.3019 11.7893 90 90 91.138
The string can be copied from KSm.txt by CtrlC + CtrlV
Change “Target dimension” to 4 and fill q vector 0.56883 -0.12885 0
The string can be copied from KSm.txt by CtrlC + CtrlV
For wavelength type 0.7114 and select “Linearly polarized beam”; NEXT; FINISH
Accept the data in Data repository
Select “Yes, I would like to continue with the wizard”; OK
3. Checking data
[On the screen: Refinement of the powder profile by the le Bail algorithm]
“Show powder profile”
Profile viewer shows experimental powder profile
The powder profile is stored in KSm.m95 and in KSm.m90 in analogy to single crystal data
Leave profile viewer
“Edit refinement commands”; 0 refinement cycles; OK
“Run Refine”
“Show powder profile”
Now the profile viewer shows experimental and calculated powder profile. The calculated
profile is based on default profile parameters. Bragg positions of satellites are green. By
default satellites up to the order 1 are used.
The calculated profile is stored in KSm.prf
As in the space group I2/b the minimal site multiplicity is 4 it is clear the Sm and K are either
alternatively occupying the same position or two positions with half occupancy. The first
possibility would change the average density at Sm/K position to a level not to different
from that at Mo position. The average atomic number of Sm/K 40.5 comparing with 42 for
Mo atoms. This means that then we could have problem to distinguish Mo and Sm atom.
Select “Repeat Superflip: Number of runs”, type 10 runs and 3000 cycles
Select “Starting model” = “Patterson superposition map”
Peak from Jana2006 but first run Fourier
“Run solution”
Start “Edit atoms” or “Edit/View → Editing of M40” and make sure the atomic positions
returned by Superflip are as follows (or symmetry transformed):
The position of heavy atoms makes situation even more complicated and it induces
additional translation symmetry (1/2,1/2,0) amd (0,0,1/2)
This will unify the origin selection with the selection used in the manual.
The positive maximum in the difference map indicates region with lack of electron density.
This is the area where fully occupied samarium should be present (instead of the partially
occupied one). From this map we can conclude that the definition interval of potassium has
centre in x4=0 and width 0.5. Samarium would occupy a complementary interval with centre
at x4=0.5 and width 0.5.
Leave Contour and start “Edit atoms”
Right-click Sm1 and select “Split atomic position”
In the dialogue that follows:
Type “K1” for “Name of the split atom”
Select “K” for “Atomic type”
Type “0.5” for “Relative occupancy of the split atom”
Type “0 0 0” for “Shift of the split atom”
Clear “Generate restrict command for Refine”
Apply
Switch to the page “Edit”, ensure that “ai” = 0.5 and its refinement is disabled
Press “Occupancy”
Type 0.5 for “delta” and 0.5 for x40, OK, OK
By this way we have defined two identical crenel functions for both atoms Sm1 and K1. Note
that the occupancy “ai” was recalculated back from 0.25 to 0.5
[On the screen: list of atoms in “Edit atoms”]
Unselect Sm1
Right-click K1, go to “Action → Edit/define atoms”
Switch to page “Edit”; press “Occupancy”, change x40 from 0.5 to 0; OK; OK
Close “Edit atoms” and save changes
Contour step in the difference map is 0.1. The residua in the difference Fourier map are smaller
and in different places than before. The position modulation now occurs in x2-x4 section. The
previously observed modulation in the x1-x4 section was caused by pseudosymmetry.
The profile width can be probably better described with anisotropic strain broadening
Start “Parameters → Powder”, go to the page “Profile”
Select “Tensor method”
Press “Edit tensor parameters”
Press “Refine all”
Leave “Powder options”
In Refinement options go to page “Basic” and select “Use Marquart technique” with the
default “Fudge factor” 0.001
Run Refine
and fix it to zero in “Edit tensor parameters”.
The refinement converges with Rp=3.8%, R(main)=1.7% and R(sat)=2.9%
The section now displays very small residual density. Depending on the fixed St components
the resulting section may look differently.
C10H16NOCl
2. Import Wizard
Select “reflection file corrected for LP and absorption”; NEXT
Select “SHELX on I”; File “Ephedrine.hkl” NEXT
The measurement was done in a supercell using three indices
[On the screen: complete/correct experimental parameters]
Fill supercell parameters as used in data collection. They can be copied (Ctrl-C) from the file
Ephedrine.ins and pasted (Ctrl-V) to the “Cell parameters” textbox; NEXT
[On the screen: Define the reference cell/split by twinning]
Change target dimension for 4
Define modulation vector 0 0 1/4
Press “Define transformation matrix” and fill the matrix array by
100
010
0 0 1/4
NEXT; FINISH
3. Data Repository
The line in data repository window corresponds to one data set.
OK; YES to accept the data set
Run EditM50 and go to page Composition, press Formula from M40. The found formula
should be C10 H15 N O Cl which means that one hydrogen atom bounded to O is missing.
This atom will be localized later from the Fourier map.
Run refinement
Refinement converges with R value about 6.2%, 3.4%, 4.9% and 31.6%.
“DO you want to re-create refinement reflection file now?” NO in order to skip re-creation
of the reflection file
Go to “Parameters/Twin fractions” and allow refinement of twvol2
Run Refine.
Refinement converges with R value about 3.8%, 2.9%, 4.1 and 8.0%.
Save the structure with a new name e.g. “Ephedrine-comm-1” but continue with the old
one.
This structure will be used for a final comparison with two commensurate cases.
On the other hand, in incommensurate structure model overlaps must be determined based
on the proximity of reflection angles
Refinement converges with R value about 3.8%, 2.9%, 4.1 and 8.0%.
Conclusion: Modulations are not very strong and therefore all three models give very similar
results
- for selected model second modulation wave can be added, but it has almost no influence on
results
- we can deactivate twinning and verify by refinement if it was really needed (it was!)
- hydrogen of OH group probably should not have individual modulation waves. They can be
made the same like modulations of the parent atom in “Refinement options → page Various →
Restrictions”
- In the incommensurate refinement the overlaps detected by the program can be printed in the
refinement listing if we activate the option in “Refinement options → page Select/Listing →
Print twin/overlap details”
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction – CCD”; NEXT
Define basic input file: select Input from “sum” file; OK
Leave all settings unchanged; NEXT
Leave all settings unchanged; NEXT
For absorption correction select “None or done before importing”; NEXT
FINISH
3. Data Repository
The line in data repository window corresponds to one data set.
OK; YES to accept the data set
4. Symmetry Wizard
Use default tolerances for crystal system recognition
Select monoclinic Laue symmetry
Select C centering
Select space group C2/m
Select superspace group C2/m(a0g)0s
Accept the superspace group in the original cell
From the plot, it is clear that the two positions of the bridging oxygen O2 are probably
separated in the modulated structure with a crenel function
Oxygen O2 has x2=-0.07, its symmetry equivalent position has x2=+0.07. Therefore we can
define a crenel function of the width 0.5 and the center 0.75 which will be transformed by
the symmetry to the crenel function with the same width, the center 0.25 and x2=+0.07. This
will describe “skipping” of the bridging oxygen between two positions separated by the shift
in x4.
Similarly like in example 5.3.2 (incommensurate phase of Cr2P2O7), we will combine crenel
functions with position modulations expressed by Legendre polynomials.
The Contour plot now shows the two crenel functions related by symmetry operation
x1 -x2 x3 x4+1/2
Using the button “From list” in Contour, plot x1-x4 sections for Cr1
Cr1 should be described with by crenel function with center 0.25 and width 1 combined with
positional modulation to create discontinuity at x4=0.75.
Similarly like for O2 define crenel function for Cr1: x40 = 0.25, delta = 1
Select “legendre polynomials in crenel interval”
Using the button “New plot” in Contour plot x2-x4 sections for O1
9. Refinement
With commensurate structures, crenel function can be used but its center and width cannot
be refined, as their changes give zero derivations, which would lead to singularity problems
during refinement.
Start “Editm50”: check “Commensurate case”, “Simple supercell”, supercell 3 1 2, Tzero = 0
Follow wizard for creating the refinement reflection file
Start Refinement commands: 100 cycles, damping 0.5
291 Jana2006 Cookbook, version October 2017
Refine
Refinement converges with 18%, 8%, 13%, 46% and 56% for all, main, 1st order, 2nd order
and 3rd order (31 refined parameters).
Add the second position modulation wave for all atoms except O2
Refine
Refinement converges with 8%, 4%, 6%, 15% and 48% for all, main, 1st order, 2nd order and
3rd order (48 refined parameters).
Define harmonic ADP for all atoms
Refine
Contour: plot relevant sections for O2, O3 and Cr1
O2 O3 Cr1
OK to plot
Options for Dist: in “Basic” page set listing “With symmetry code”; in “Select atoms” page
set central atom P1, neighbor atoms O2, O3.
Start Dist
Open the listing of Dist
Atom P1 is coordinated with O2 and O3, both in the basic symmetry position x,y,z
The upper limit of the crenel interval should be the same when expressed in t:
tmax(O3) = tmax(O2)
t40(O3) + delta(O3)/2 = t40(O2) + delta(O2)/2
t40(O3) + 0.5 = t40(O2) + 0.25
t40(O3) = t40(O2) - 0.25
Start Grapht and verify that the discontinuity in O2 nad O3 now occurs at the same t
Chemical diagram:
Introduction
Structure of this compound has been published in Journal of Science and Technology in the
Tropics (2005), 1, 157-163, as a commensurately modulated structure. The data from Nonius
KappaCCD diffractometer contained strong main reflections and weak first and second order
satellites with ratio of intensities 1000 : 0.157 : 0.075. Later on the compound was remeasured
with much more intensive source of the Rigaku Oxford Diffraction SuperNova. In the new
experiment, the satellites were strong and the structure could be solved in a five-fold supercell
as a classical three-dimensional structure. This example uses data of SuperNova and compares
possibilities of the commensurate and supercell structure refinement.
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction – CCD”; NEXT
Define basic input file: select Input from “cif_od” file; OK
Leave all settings unchanged; NEXT
Leave all settings unchanged; NEXT
For absorption correction select “None or done before importing”; NEXT
FINISH
3. Data Repository
The line in data repository window corresponds to one data set.
OK; YES to accept the data set
4. Symmetry Wizard
Use default tolerances for crystal system recognition
Select monoclinic Laue symmetry
Select P centering
The structure model returned by charge flipping contains most of the expected structure
and it could be finished by adding maxima from difference Fourier and possibly by applying
geometry constraints imposed during the refinement. However, we would like to try
commensurate refinement in order to see its pros and cons.
NEXT
The wizard transforms indices according to the defined transformation matrix. Resulting
integer indices are sorted as main reflections while the others as satellites defined by the q-
vector.
For absorption correction select “None or done before importing”; NEXT
FINISH
OK; YES to accept the data set in the Data repository
8. Symmetry Wizard
Use default tolerances for crystal system recognition
Select monoclinic Laue symmetry
Select P centering
Select space group P21/c
Although it is not first in the list, FOM is almost the same like for P2/n and the button
"Details" shows that the discarded reflections are very weak
NEXT
Select superspace group P21/c(a0g)0s; NEXT
Accept the superspace group transformed into the original cell; NEXT
The structure model obtained from charge flipping contains fragments of phenyls and
phenantroline, Sn, Cl and water oxygen assigned as carbon C1. In commensurate description
we have only one symmetry independent Sn.
Start "Edit atoms"; change chemical type of C1 (may be also N1) to oxygen and its label to
O1;
OK; YES to save changes
Start "EditM50"
Check "Commensurate case"
Clear "Use simple supercell"
Press "Define supercell matrix" and type 1 0 1 ; 0 -1 0 ; 3 0 -2
301 Jana2006 Cookbook, version October 2017
This is transformation from commensurate to supercell parameters
OK
Press "Select supercell group"
The program lists all supercell space groups, which can be obtained by changing the tzero
value. As we know from the preliminary calculations that the symmetry in the supercell
should be P21/n, we should select one of the two equivalent possibilities.
Select Tzero 0+n/5; OK
OK; YES to rewrite M50; YES to re-create refinement reflection file
Follow steps for merging the reflections; FINISH
Now we have set up the commensurate refinement. In our case, it means that modulation
functions will only have 5 definition points.
Start refinement
The refinement converges with Rall~38%, Rmain~18% (your numbers can be different)
Start "Edit atoms"
Double-click Sn1 and change number of position modulation waves to 2
OK; OK; YES to save changes
Start Refinement options; set 100 refinement cycles; OK; YES+START
The refinement converges with Rall~24%, Rmain~14%, Rsat1~27% and Rsat2~31%
These numbers can be different because Superflip may have returned different number of
carbon atoms. However, if R values at this point are significantly larger, for instance Rall >
26% or Rsat2 ~ 50%, please do not continue and repeat structure solution (return to the
point 10.).
Quit Diamond or let it open for a later reference (depends on style of work)
Press "Apply+End"
Start "Edit atoms"
Delete all carbon atoms, which are not blue (i.e. carbons out of molecules).
Be careful not to delete chlorine
Press "Select all"
The structure should contain 36 atoms
Press "Refresh"; OK; YES to save changes
Plot the structure using "Draw average structure
Each "t" value indicates one of the defined commensurate sections, which is realized as one
configuration in the five-fold superstructure. We can see that changes of angles between the
planes #1, #2 and #2, #3 reach 30 degrees, which cannot be reliably described by rigid body
modulation function as used in Jana2006.
308 Jana2006 Cookbook, version October 2017
Start “Grapht”
Press “New/edit”
For “Parameter to draw” set “plane”
For “p” set limits from 0 to 90
For “Plane #1” type “C1pa C2pa C3pa C4pa C5pa C6pa”
For Plane #2” type “0 0 1”
Grapht now displays angle between the plane through the first position of Ph and the axis c.
Press “Options”; check “Indicate commensurate t sections”; OK
Grapht now also displays the sections, for which the structure is defined. The orientation of
the single plane changes by ~40 degrees. Again, this cannot be reliably described by rigid
body modulation function as used in Jana2006.
Actual value of t-sections is connected with the origin selection made automatically during
solution by Superflip. Also the phenyl molecules can be placed in different order compared
with this manual. Thus the actual plot will differ.
The plot indicates that for one of the sections (here for t=0; its actual value is connected with
the origin selection made automatically during solution by Superflip) some of the
phenantroline atoms are strongly disordered, while the other not. This could not be
described with the rigid body model.
Quit Grapht
Draw the average structure, select only the phenantroline molecule
For this Ctrl-M keyboard shortcut of Diamond can be used, which selects molecule belonging
to a selected atom
Note atoms N1t, C18t, C17t and N2t.
Phenantroline is often twisted and we will calculate torsion angle of these atoms
21. Water
The last missing atoms in our structure model are hydrogens of the water molecule
Start Fourier options
For "Map type" select difference Fourier
Check "Apply sin(th)/lambda limits
define "max" as 0.5
OK; YES+START
YES to start procedure for including new atoms
Add two difference maxima (probably the first two ones) 0.82Å close to O1 and name them
H1O1 and H2O1
Start "Parameters → Molecules → New molecule"
Select "Atomic part" and "water" as a name; Next
Select O1, H1O1 and H2O1; OK
For the reference point select "Explicit" and use O1; NEXT
Apply+End
Start "Edit atoms"
Delete all modulation waves for O1, H1O1 and H2O1; OK; YES to save changes
Start "Parameters → Molecules → Edit"
Set one position modulation wave for "water"; OK; YES to save changes
Start Refinement options
In page "Various" define distance restraint 0.82Å for O1-H1O1 and O1-H2O2
In page "Various" define keep command for ADP of hydrogens H1O1, H2O1, to be 1.2
multiple of Ueq(O1).
24. Conclusions
Commensurate refinement can help to find solution of very large structures, it allows to save
parameters by using less than the full set of modulation parameters, and it facilitates
interpretation of a superstructure through calculation of various parameters as a function of
t. However, commensurate refinement may be difficult when satellites are too strong,
because strong satellites indicate large displacements, which must be described by
modulation functions defined only in few points. Despite of this, commensurate refinement
was possible for our structure and it helped us to understand its geometry. On the other
hand, the phenyl rigid body could be only possible in the supercell refinement, as well as the
full description of the disordered phenantroline molecule.
N-(3-nitrobenzoyl)-N’,N’’-bis(tert-butyl)phosphoric triamide
Chemical formula: C90H150N24O24P6
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction – CCD”; NEXT
Leave all settings unchanged; NEXT
Leave all settings unchanged; NEXT
For absorption correction select “None or done before importing”; NEXT
FINISH
3. Data Repository
Each line in data repository window corresponds to one data set.
OK; YES to accept the data set
4. Symmetry Wizard
[On the screen: “Information” about symmetry wizard]
This is unrealistic value which indicates that the structure can contain more independent
structural motifs. As the density should be larger than 1 the most probable number of
independent structural motifs is 6.
Set "Formula units" to 6;
Select “Superflip”, “Peaks from Jana2006”; leave other settings default; “Run solution”
Close listing of Superflip
Press Accept last solution” to leave the structure solution wizard
8. Refinement
[On the screen: basic window of Jana]
Right-click the icon of Refine.
316 Jana2006 Cookbook, version October 2017
[On the screen: Refinement commands]
Define 100 of cycles; OK
Choose “YES+START”
Refinement converges with R(obs) value about 20-22%
12. Refinement
[On the screen: basic window of Jana]
Double-click the icon Refine.
Refinement converges with R(obs) value about 16%
Start “Edit atoms”
[On the screen: list of atoms]
Press “Select all” and “Action → Edit/define”
Choose “harmonic” for ADP parameters. This defines anisotropic temperature parameters
for all atoms. OK, OK, YES to save changes
Double-click the icon Refine.
Refinement converges with R(obs) value about 9%
e
2
2
.
5
1
.
5
0
.
5
-
0
.
5
-
1
.
5
e
2
5
0
.
5
-
0
.
5
-
1
.
5
Locate the split positions by clicking at corresponding points of the section and save them
for future use by “Save the selected point for future use” button.
Make the same for the second molecule (the section C132, C142, C152). Here the splitting
of atoms C142 and C152 is less pronounced, however we only need some points for initial
orientation of rigid bodies.
Make the same for the third molecule (the section C133, C143, C153). Here all positions are
clearly split but we only need two pairs of points..
The split positions localized during preparation of this example were also written to the file
“TerBut-points.txt” and can be later to pasted to Jana2006 textboxes.
17. Verification
Use Diamond for verification of proper geometry: each molecule should have one single
and one disordered tertiary group.
Use Parameters → Molecules → Edit for verification that positions #4 – #9 have occupancy
(i.e. the parameter “aimol”) 0.5
Use “EditM50”, page “Composition”, button “Formula from M40”
If occupancy of molecular positions has been properly defined and the coincidence ratio has
been always 5/5 the formula should be C15 N4 O4 P H0 (for Z=6)
Use “Edit atoms” to see changes in the atomic part of M40
Black atoms are the ones not used in molecular description
Blue atoms are the ones defining the pattern molecule
Start “File → Structure → Save As” and create a backup copy “rigid_ini” but do not
continue with the new structure.
18. Refinement
[On the screen: basic window of Jana]
Right-click the icon of Refine.
[On the screen: Refinement commands]
Change the damping factor to 0.5;
OK, “YES+START”
Refinement converges with R(obs) value about 10.4%
Resulting number of atoms: 42 in atomic part (here hydrogens are still missing); 14 in the
pattern molecule.
26. Refining of the structure with anisotropic ADP for tertiary butyls
Start “Edit Atoms”;
[On the screen: tool for selection of atoms]
Select non-hydrogen atoms of the pattern molecule.
Double-click; Select "Harmonic ADP";
OK; OK; Rewrite changes
Start refinement;
Refinement converges with R(obs) value about 9.4% for 466 parameters
27. Refining of the structure with TLS tensors for tertiary butyls
Start “Edit Atoms”;
[On the screen: tool for selection of atoms]
Select all atoms of the pattern molecule.
Select “Use TLS” using double-click;
OK; OK; Rewrite changes
Start refinement;
Refinement converges with R(obs) value about 7.5% for 556 parameters
31. Refinement of the structure with local symmetry 31m – atomic positions of the tertiary
butyl fixed
Start the refinement;
Refinement converges with R(obs) value about 6.6% for 556 parameters, i.e. the same value
like before, because the shape of the pattern molecule is still fixed.
(3aRS,7aRS )-1,3-dinitrosooctahydro-1H-benzimidazole
Chemical formula: C7H12N4O2
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction – CCD”; NEXT
Leave all settings unchanged; NEXT
Leave all settings unchanged; NEXT
For absorption correction select “None or done before importing”; NEXT
FINISH
3. Data Repository
Each line in data repository window corresponds to one data set.
OK; YES to accept the data set
4. Symmetry Wizard
[On the screen: “Tolerances for crystal system recognition”.]
Leave all settings default; NEXT
[On the screen: Select Laue symmetry]
Select Monoclinic; NEXT
Select primitive unit cell; NEXT
[On the screen: Select space group]
Select P2/c; NEXT
The results may be slightly different for each run of Superflip, because it starts from
random phases.
Draw rectangle to select all atoms of the molecule or press Ctrl-A
Right-click on one of the selected atoms, choose “Add → Atom labels”
[On the screen: dialogue of Diamond for atom labeling]
Note incorrect chemical types (using the chemical diagram above);
Quit Diamond;
Run “Edit atoms”;
Correct the improper chemical types;
Rename atoms according to chemical types;
Make symmetry contiguous motifs;
Quit “Edit atoms”, save the changes
For this structure making the symmetry contiguous motifs is very important. To make sure
the motifs have been correctly created:
If these motifs are not created the atoms must be put to the motif manually using
symmetry transformations in "Edit atoms". Without proper motifs this example will not
work!
Quit Diamond
8. Refinement
[On the screen: basic window of Jana]
Right-click the icon of Refine.
[On the screen: Refinement commands]
Define 100 of cycles; OK
Choose “YES+START”
Refinement converges with R(obs) value about 13%
Plot the structure
It should be the same like before
Start “Edit atoms” and change ADP of all atoms to “harmonic”
Repeat refinement
Refinement converges with R(obs) value about 9.4%
Plot the structure
Activate ADP ellipsoids using “Picture → Model and Radii”
One of the molecules has very large ellipsoids for carbon atoms:
9. Labeling
At this point we would need some more practical labeling of atoms. Jana2006 doesn’t
include tools for graphical renaming. The labels and their changes must be planned in
Diamond, noted and realized in “Edit atoms” tool or by direct editing of M40
The desired labeling for the molecule with large ADP’s (starts from C1) and the one with
smaller ADP’s (starts from C5):
The splitting procedure created pairs C2-C2’ and C4-C4’, their ADP are isotropic and the
same for each pair, partial occupations are ½ and the sum of occupancies for each pair is
kept by a Restriction for the case it would be refined
Run refinement
Refinement converges with R(obs) value about 10%
Start “Edit atoms” and activate refinement of occupancy for C2 and C4
Run refinement
Refinement converges with R(obs) value about 8.4%
The refined occupancies in both pairs are in the ratio 1:3, i.e. one of the split atom in each
pair is much weaker.
Plot the structure; distinguish by colors the weakly occupied carbon atoms
Activate plotting of ADP ellipsoids and plotting of labels
Warning: splitting of atomic positions was automatic process. In particular solution it may
happen that the labels of the strongly and weakly occupied atoms (i.e. atoms with and
without comma) will be exchanged. In such a case we recommend relabeling the atoms
according to this plot in order to follow easily the cookbook.
’
For “Occupancy” type 0.75;
“Apply + Next position”
[On the screen: Define and complete molecular position #2]
In this dialogue we will define rotations and a translation vector for placing the model
molecule to the second actual position. This is done by defining three atoms of the model
molecule and three atoms in the atomic part of the structure model, which should
coincide. In order to get correct placement of the second position the split atom positions
must be used.
Use C1, C2’, C4’ as model atoms and C1, C2 and C4 for the actual position
Press “Show coinciding atoms”
We expected that the 7 atoms in the actual position will coincide with 7 atoms in the
atomic part. However, we obtained only 3 coinciding atoms .
OK to quit information window
Select “Apply inversion”; Press “Show coinciding atoms”
Note that all positions except C1 are now more or less split. Close positions in the rigid
body model are not problematic because the individual atomic parameters will be refined
only once for the model molecule.
In our case, if we define the local symmetry 1, still keeping only one half of atoms in the
model molecule, the molecule would be completed after correct placement on a two-fold
crystallographic axis. However, we could not keep hydrogens of the carbon C1 in
geometrically expected positions because the keep command applies to the model
molecule and hydrogen positions need to be generated by symmetry. For this reason, we
have to define the local symmetry.
In our case, the model molecule has a two-fold axis along b and its point group will be 2.
Because the internal definition of the point group expects rotation along c we have to
define the point group 2 and change the local coordinate system.
Start EditM50
Go to the page “Composition”
Press “Formula from M40”
It should print C7 H12 N4 O2 for Z=4
Start Refinement options, page Basic
Change Instability factor from 0.01 to 0.02
Run refinement
Refinement converges with R(obs)=3.9%, GOF=2.0.
Ideally, the instability factor should be determined for a given instrument. In the practice
we know that it ranges between 0.01 and 0.02.
GOF about 2.0 may indicate that disorder is not described exactly. It may also indicate that
the data contains information about bonding effects.
Start Fourier options
In page Basic define “Map type” “Fobs-Fcalc difference Fourier”
Run difference Fourier
In the listing of Fourier check list of maxima and minima
The residua are small like for fully solved structure
Last step, which is skipped in this example would be comparison of the molecular model
with a model built from split atom positions. The resulting R value would be very similar to
that of the molecular model.
The purpose of this exercise is to understand the log-file of Superflip, point out the most
important information to check, and to learn some basic tricks for the case Superflip does not
converge with default parameters. The structure used for the exercise is an organic structure
with symmetry Pna21.
Plot the structure, if you wish, check EditM50 to see the composition and symmetry.
2. Run Superflip.
The first part is the header of the file, with information about the version. If you encounter
strange behavior of Superflip, make sure in http://superflip.fzu.cz you have the newest
version.The executable of Superflip downloaded from this web page should be copied to the
SUPERFLIP directory of Jana2006 installation.
**************************************************************
Superflip
Version: 12/04/13 7:57
Palatinus, L. & Chapuis, G.(2007): J. Appl. Cryst. 40, 786-790
http://superspace.epfl.ch/superflip
**************************************************************
-----------------------------------------------
Start of the calculation: 28.JAN 2014, 14:29:49
-----------------------------------------------
###################################################################
# Following data were read from the input file or set as default: #
###################################################################
------------------------
-----------------------------
Crystallographic information:
-----------------------------
Superspace dimension: 3
Dimension of the physical space: 3
Direct cell parameters: 15.3079 9.2915 13.8059 90.0000 90.0000 90.0000 Volume: 1963.65942
Reciprocal cell parameters: 0.0653 0.1076 0.0724 90.0000 90.0000 90.0000 Volume: 0.00051
Wilson plot will be based on this composition: C84H104N8O16
4 symmetry operations found, their list follows:
4: -1 0 0 0.5000
0 1 0 0.5000
0 0 1 0.5000
--------------------------
Settings of the algorithm:
--------------------------
Number of voxels: 40 24 36 Total: 34560
Number of voxels for polishing: 80 48 72 Total: 276480
The density modification method will be charge flipping.
Delta will be determined automatically.
The iteration will be stopped when the convergence is detected or after ***** cycles.
After the regular iteration 5 cycles of "density polishing" will be applied.
The random number generator will be initialized automatically.
The resulting density will be shifted and averaged according to the symmetry operations given above.
Following symmetry operations will be used to locate the origin of symmetry: 3 4
Most of the information above is just diagnostic and is not important for ordinary user.
However, the number of reflections treated as weak might be changed, if Superflip does not
converge with default parameters. See below for more discussion on this parameter.
-------------------------
Averaging of reflections:
-------------------------
Redundancy : 1.000
Rint : -----
Wilson plot is always a useful indicator of the data quality and consistency. If you encounter
problems with your structure, do not forget to check the Wilson plot and make sure it does
not show anomalies. The value Biso is the Wilson-plot based estimate of the average Debye-
Waller factor. If it is larger than 2, a correction of the input data might be needed. We will
learn how to make this correction later in this exercise.
------------
Wilson Plot:
------------
circles: <I>/SUM(f^2)
dots : linear fit y = 1.5450 + -3.480 * [sin(theta)/lambda]^2
4.15377-| o o
-| . . o
Scale^2 -| o o . o o
-| o o . o
-| . o o o
2.99073-| o . . o o
-| o . . o
-| . . o
-| o . .
-| o o . .
2.15334-| . . o
-| . o
-| o o
-| .
-| o . .
1.55041-| . .
-| o . . o
-| o o .
-| o
-| o
1.11630-| o
I-------------------I-------------------I-------------------I----------
0.00 0.10 0.20 0.30
[sin(theta)/lambda]^2
Coverage statistics shows the coverage as a function of resolution. The coverage should be
high (optimally more than 90%) up to a reasonable resolution (dmin<1.1 or better). If one of
these conditions is not fulfilled, Superflip may fail. Very low coverage at low resolution is
especially dangerous.
----------------------------------------------------------
Coverage statistics of the expanded reflections by shells:
----------------------------------------------------------
Resolution (sin(th)/l): 0.050 0.100 0.150 0.200 0.250 0.300 0.350 0.400
Resolution (d_min): 10.000 5.000 3.333 2.500 2.000 1.667 1.429 1.250
Obs. refl. in shell: 0 25 69 139 238 349 515 655
Total refl. in shell: 0 27 73 139 238 349 515 655
Coverage in shell: 0.0% 92.6% 94.5% 100.0% 100.0% 100.0% 100.0% 100.0%
Commulative coverage: 0.0% 92.6% 94.0% 97.5% 98.7% 99.3% 99.6% 99.7%
At this moment the calculation is converged (or the maximum number of iteration cycles
was reached). Superflip has the capacity to check the resulting density for symmetry. It is
advisable to always check this part of the log-file, and to make sure that the detected
symmetry corresponds to your expectations. Derived space group P1 in cases where higher
symmetry is expected usually indicates that the structure is not solved.
Superflip lists all symmetry operations compatible with the lattice. Each line in the listing
corresponds to one symmetry operation. The symbol of the operation is followed by a one-
line representation as used in International Tables for crystallography. Next is the
“symmetry agreement factor”. This number indicates quantitatively, how well the operation
is present in the electron density. Numbers below 10 indicate a very good match, number
between 10 and 25 are still probable symmetry operations of the density, and numbers
above 25 usually mean that the corresponding symmetry operation is not present. In the
343 Jana2006 Cookbook, version October 2017
example below the n-glide, a-glide, and two-fold screw axis have agreement factors below
0.1, indicating an excellent agreement. All the rest of the list has a far too high value. The
rows of XXXXXX after the symmetry agreement factor are just a graphical representation of
the agreement factor, and serve as a quick guide for eye.
#####################################
# Checking the density for symmetry #
#####################################
An experienced crystallographer can “see” the space group already in the list of the
symmetry operations, but Superflip can analyze the list automatically. The default threshold
for accepting a symmetry operation is the symmetry agreement factor below 25. All such
operations are combined together, and a final space group is printed as Herman-Mauguin
symbol, Hall symbol, Fingerprint, and a list of symmetry operations. The “Fingerprint” is an
internal space-group symbol used by Superflip for a database search, and currently does not
have any other use.
-------------------------------------------------
Space group derived from the symmetry operations:
-------------------------------------------------
HM symbol: Pna21
Hall symbol: P 2c -2n
Fingerprint: 3300223}240qY3 (3/4,3/4,0)
Symmetry operations:
1: x1 x2 x3
2_1(0,0,1): -x1 -x2 1/2+x3
n(1,0,0): 1/2-x1 1/2+x2 1/2+x3
a(0,1,0): 1/2+x1 1/2-x2 x3
As the density is reconstructed in P1, the origin is arbitrary. To allow for further processing,
the density must be shifted to the standard origin of the space group, and averaged to have
the desired symmetry exactly. The lines below show the record of this process.
Important note: The “Check for symmetry” described above is just indicative, and does not
have any impact on further symmetry processing (unless you force it explicitly in the input
file of Superflip – see manual, keyword derivesymmetry for details). Thus, the output
symmetry will be the symmetry indicated in the input file. If you call Superflip from Jana with
default parameters, this is the symmetry recorded in m50.
“Agreement factors of individual generators” show agreement factors for the generators of
the input space groups. These generators are used to locate the origin of the space group. If
the input space group is incorrect, or the solution did not converge, these agreement factors
will be high (more than 25). Low values of the agreement factors are another reliable
indicator that the structure was correctly solved.
These are the final notes in the log file with information about the file with the output:
You can obtain more information about the reconstructed reflection phases by using 'expandedlog yes'.
---------------------------------------------
End of the calculation: 18.JAN 2010, 13:17:29
---------------------------------------------
Superflip version: 07/16/09 09:12
Press “Close” to close the listing and return to the Structure solution dialogue
If everything goes well, the structure model correctly describes the density produced by
Superflip. However, in some cases this can lead you to the wrong impression that this is
what Superflip produces. In case of problems (typically disorder or wrong chemical types in
the formula) viewing the density may be very helpful.
Press “Draw 3d map” to view the 3D isosurface representation of the raw density produced
by Superflip
The crystallographic (not only) viewer VESTA must be installed on your computer and
enabled it in Jana. If not, consult the installation notes and install the program first. Jana
passes the volumetric data to VESTA along with the structure model. VESTA shows an
overlay of both, with the density represented by isosurfaces, and atoms by spheres. In this
case the density shows nice spherical maxima at the positions of atoms. All maxima are
“occupied” by an atom. Thus, the structure model is likely to be complete and ready to be
refined.
By default, the structure factors from Superflip are saved in the file jobname.m80. That
means the Fourier map can be recalculated by Jana2006 as long as m80 is not overwritten
by refinement.
You may select a different number of repetitions than the default 10.
If necessary, modify the maximum number of cycles per run, for example to 5000:
This option may help if the solution appears unstable or if it gets trapped in a local
minimum. Some of the repeated runs may escape from the local minimum. Superflip selects
the best of all attempts for the output.
With this option Superflip first calculates a Patterson superposition map and then uses this
map as a starting model for the iteration. Very often this option speeds up the convergence,
and sometimes it makes the difference between solved and non-solved structure. However,
sometimes – especially for simpler structures – the convergence is so rapid that it cannot be
detected and Superflip keeps running. This option saves time only for larger structures,
ideally with at leas one non-light atom.
The AAR algorithm appears to perform better than the standard algorithm for very
incomplete data or for low-resolution data. However, it is less stable than the standard
algorithm and it is more difficult to detect its convergence. AAR is the default algorithm for
electron diffraction data.
346 Jana2006 Cookbook, version October 2017
Try one or more of the above hints and see how they affect the speed of convergence. You
should come to the conclusion that for this simple structure there is hardly any difference
between the options, but most of them speed up the convergence., However, combination
AAR + local normalization is too aggressive and does not lead to the solution.
Superflip offers a large number of other options, which are described in detail in the manual.
If you want to use any of these advanced options, then:
Check the checkbox “Allow manual editing of the command file before start:
Before Superflip starts, the input file opens in a text editor and you may edit the file and
include whatever additional options you wish.
The price to pay for faster convergence is the increased noise in the result, and sometimes
more difficult detection of convergence. Because the default parameters work usually quite
well and yield the best densities, we recommend using them for first attempts. Only if the
iteration does not converge, try modifying the parameters according to the hints described
above.
Data Import
Output: M95
Symmetry wizard
Updates symmetry in M50
Structure solution
(Superflip calculates density, Jana interprets it)
Draw structure
Reads M50,M40
Draw 3d map
Reads M81,M40
*) Superflip expands M90 to P1 symmetry. However, its results are to some extent
influenced by the symmetry set by the user with the Symmetry wizard, because M90
contains merged reflection and systematically extinct reflections are omitted
This relatively simple structure illustrates the advantages of using Superflip to determine the
space group.
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction – CCD”; NEXT
Leave all settings unchanged; NEXT
OK to confirm the name of the basic input file
Leave all settings unchanged (note the cell parameters and thus the probable crystal
system); NEXT
Finish the import wizard with all settings default
3. Symmetry Wizard
[On screen: Information about the symmetry wizzard.]
NEXT to start the wizard;
Leave all settings default; NEXT; OK;
[On screen: Select Laue symmetry]
Jana2006 suggests orthorhombic symmetry based on Rint. In this case the choice appears
clear.
Keep the suggested selection; NEXT
Select primitive unit cell; NEXT
[On screen: Select space group]
You are now offered a list of candidate space groups sorted by their likelihood. Note that
the first two space groups Pnma and Pn21a differ only by the presence of inversion center
and they have the same FOM:
Experimental Theoretical
Value All data hkl 0kl h0l hk0 Acentric Centric Hypercen.
This table shows the numbers typically used to estimate the centrosymmetry. If you
compare the column All data with the theoretical values, you can see that all numbers
indicate a centric or hypercentric distribution. Our guess Pnma thus appears to be strongly
supported by the statistics
Click the button “Close” to close the listing and return to the main window.
6. Structure Solution
#####################################
# Checking the density for symmetry #
#####################################
-------------------------------------------------
Space group derived from the symmetry operations:
-------------------------------------------------
HM symbol: Pn21a
Hall symbol: P -2ac -2n
Fingerprint: 33002241Y40qY3 (3/4,0,3/4)
Symmetry operations:
1: x1 x2 x3
a(0,0,1): 1/2+x1 x2 1/2-x3
n(1,0,0): 1/2-x1 1/2+x2 1/2+x3
2_1(0,1,0): -x1 1/2+x2 -x3
Look at the table. You can see that the first three operators belonging to the Pn21a space
group have very low agreement factors (less than 1%) indicating a very good match. All
other operators have agreement factors above 50%, which is very poor. In particular the
mirror plane perpendicular to b has agreement factor 71.7%. It is thus clearly absent in the
structure. The space group is without doubt Pn21m.
Click “Run solution” to interpret the previous solution with the new setting.
A new screen appers with the run-time information from EDMA: When EDMA finishes, a
listing of Superflip opens.
Close the listing and return to the structure solution dialog
Double click “Draw 3D map”
This time you sould see all maxima in the map occupied by an atom. Most atoms should
have a correctly assigned atomic type, although carbon, nitrogen and oxygen may be
occasionally swapped.
Here is the correct molecule for reference (hydrogens omitted for clarity):
Input files: flo19.cif, flo19.hkl. The structure is an organic molecule built by four aromatic
rings bridged by one or two nitrogen atoms.
1. 1If you accidentally close the wizard or if you want to repeat the test, you can access the
symmetry determination procedure through the EditM50 icon, card Symmetry.
355 Jana2006 Cookbook, version October 2017
Select the space group P21221 and click “Details”;
In the listing you can find basic information about the space group and the list of observed
reflections that violate the rules for systematic absences. Check the details for several space
groups. Which space groups would you consider possible and which would be the most
likely?
Because the solution in Superflip is not sensitive to the choice of the space group symmetry,
we will choose just one of the space groups without any systematic absences to keep the
solution unbiased.
Select the space group Pmmm; NEXT; FINISH; OK; YES to finish the space group selection
We selected the space group Pmmm not because it is the most likely choice, but because it
does not induce any systematic absences. We will now process the reflections, and ask
Superflip to estimate the space group. If we selected a space group with systematic
absences, some reflections would be excluded from the reflection file, and the result from
Superflip could be biased.
7. Run Superflip.
Next Jana offers you the dialog for structure solution. If you want to open this dialog later,
just double click the icon “Structure solution” in the main window of Jana.
Make sure Superflip is selected;
Leave other settings default; Click Run solution;
Superflip starts and runs for many thousand cycles. No convergence is detected.
Press “Interrupt” to stop Superflip;
The log file of Superflip opens. Check the symmetry determination part. You will see that no
symmetry operation has agreement factor lower than 25% and thus the derived space group
is P1. Because we expect higher symmery (although we do not know the exact space group),
it is a sign that the calculation did not converge.
356 Jana2006 Cookbook, version October 2017
Inspection of the log file does not reveal any serious problem, but a series of hints shows
that the structure can be difficult to solve. The potential problems are that the coverage
starts dropping already at dmin≈1.1 Å, and that the coverage at low resolution is very low
(only 21% at dmin≈ 5 Å). It can be assumed that quite a few strong reflections are missing in
the dataset. The first change to try is to use normalized intensities. When using normalized
intensities, Superflip also automatically treats the missing reflections in a special way that
improves the performance in case of data with missing low-order reflections.
Close the view of the log file;
Check the checkbox “local normalization”
Click “Run solution”
Superflip starts and probably runs again for many thousand cycles without detecting
convergence 2.
If you look at the evolution of the convergence indictors, you will most probably see behavior
similar to this2,3:
20 R: 56.175 Charge: 20.84 Peaks: 1.11 Score: ---
30 R: 57.537 Charge: 20.13 Peaks: 1.16 Score: ---
40 R: 58.291 Charge: 19.03 Peaks: 1.15 Score: ---
50 R: 58.075 Charge: 19.08 Peaks: 1.22 Score: ---
60 R: 57.460 Charge: 18.65 Peaks: 1.22 Score: ---
70 R: 59.161 Charge: 18.79 Peaks: 1.32 Score: ---
80 R: 57.279 Charge: 18.88 Peaks: 1.33 Score: ---
90 R: 59.624 Charge: 18.89 Peaks: 1.39 Score: ---
100 R: 56.954 Charge: 18.29 Peaks: 1.41 Score: ---
200 R: 55.838 Charge: 18.37 Peaks: 1.29 Score: 0.15
300 R: 57.292 Charge: 18.12 Peaks: 1.35 Score: 0.15
400 R: 55.889 Charge: 17.89 Peaks: 1.32 Score: 0.15
500 R: 55.456 Charge: 18.50 Peaks: 1.32 Score: 0.15
600 R: 55.902 Charge: 18.88 Peaks: 1.22 Score: 0.15
700 R: 53.192 Charge: 18.50 Peaks: 1.24 Score: 0.15
800 R: 53.102 Charge: 18.23 Peaks: 1.27 Score: 0.15
900 R: 50.902 Charge: 17.41 Peaks: 1.23 Score: 0.15
1000 R: 45.968 Charge: 18.16 Peaks: 1.19 Score: 0.15
2000 R: 43.139 Charge: 18.59 Peaks: 1.12 Score: 0.15
3000 R: 42.466 Charge: 18.12 Peaks: 1.11 Score: 0.15
4000 R: 43.158 Charge: 19.33 Peaks: 1.16 Score: 0.15
5000 R: 43.369 Charge: 18.80 Peaks: 1.11 Score: 0.15
6000 R: 44.689 Charge: 18.11 Peaks: 1.12 Score: 0.15
Note the bold part. We can see a clear drop in R-value accompanied by no change in the
charge and a small decrease (instead of expected increase) of the Peaks indicator. Although
the change in R-value is clearly visible, together with the decrease in Peaks it is not enough
to trigger automatic convergence detection. However, it is likely that the iteration
converged.
Press “Interrupt” to stop Superflip;
In the log file locate the symmetry determination part.
You will see a listing similar to this:
2
Occasionally it may detect convergence, but in most cases it does not.
3
Rarely the calculation does not converge even after several thousand cycles. In that case interrupt Superflip and
run it again.
357 Jana2006 Cookbook, version October 2017
Go carefully through the list, and make sure you understand, why Superflip suggests the
space group P21221. Note also that the space group has nothing to do with the Pmmm space
group used as input. This confirms that the symmetry derivation by Superflip is indeed
independent from the input space group.
Note also that the space group found by Superflip is the same as the most likely space group
suggested by Jana2006 on the basis of systematic absences. Thus, in this case the two
indications agree. However, the analysis of systematic absences was not conclusive, while
the result of Superflip leaves no doubt about the correct space group.
Close Diamond;
Double click the icon “Refine”
The refinement will converge to an R-value around 10% or less depending on the number of
incorrectly identified atoms. The structure is solved. Now the atoms must be correctly
assigned to chemical types, and the structure refined in the standard manner.
(Mn(H2O)4)(VO)2(PO4)2
Single crystal data measured with Oxford Diffraction four-circle diffractometer
Input files: MnVPO.m50, MnVPO.m95
You opened a structure which contains crystal information about the compound and
hkl-data. There is no structural information present yet. To preserve the original files, save
the structure under another name:
2. Symmetry Wizard
“File → Reflection file → Make space group test”
Proceed with the space group test. Check the table of the Laue groups, and note Rint for Laue
groups 4/m and 4/mmm. Continue until the screen with the title “Select space group”. You
should see this output:
In this case it is impossible to determine the space group only from systematic absences.
Select the space group I4/m, NEXT, FINISH
Superflip suggests the space group I4/mmm. This space group is incompatible with the
detected Laue group (Check your notes with Rint of the Laue groups 4/m and 4/mmm). It is
a supergroup of all three possible space groups suggested by the symmetry wizard. This is
an indication that the structure is peculiar in some respect.
Note: Superflip suggests a different space group, but it does not use this space group to
produce the final density. The final result has always the symmetry defined in EditM50. If
you want to accept the symmetry suggested by Superflip, you have to change the symmetry
and rerun Superflip to get new solution.
The structure looks correct. To convince yourself that the cations have correct coordinations,
rotate the structure and increase the accepted V-O bond length to 2 (“Build →
Connectivity”). Try to understand, why Superflip suggested a wrong space group before you
go to the next page and read the solution.
If you look closely at the list of symmetry operators and their agreement factors, you will see
indeed that the symmetry operators of the space group I4/m are all at the beginning of the
list with the lowest agreement factors. But the difference in agreement factors between the
true opeators and the remaining operators of the I4/mmm space group is too small to be
detected by the program.
This exercise demonstrates that the best approach to determining the space group is to
combine the information from diffraction pattern and the symmetry determination of
Superflip. In most cases the two approaches yield the same result. If the results differ, the
structure or data usually exhibit some non-standard features. It can be a strong
pseudosymmetry (as in this example), problem in the evaluation of systematic absences, or
twinning. The two approaches are complementary, and have different weaknesses. Their
combination gives the best result.
7. Additional exercises
a) Try to solve and refine the structure in I4/mmm. Apart from the high Rint, would you be
able to recognize that something is wrong with the structure model? What R-value do
you get?
b) For training purposes you can try to finish off the structure – check for the presence of
hydrogen atoms, define them and refine.
Input files: 8rings.hkl, 8rings.p4p. This is a huge organometallic structure with three eight-
membered and one three-membered ring containing chromium, nickel, iron and cobalt. Cell
volume 67312 Å3.
4. Structure Solution
As the next step the wizard suggests starting the structure solution. When you accept the
suggestion, a form with options for structure solution opens.
This is a large structure, and normalizing the intensities is most likely necessary.
Check the checkbox “Use local normalization”.
It is also necessary to give some composition. It is not used for solution but for interpretation. It
is sufficient to enter just a list of present atoms. Iron, cobalt and nickel will most probably be
indistinguishable in the solution, but we may enter them anyway, so that we do not have to add
them later during the refinement process. The same holds for hydrogen.
In the files Composition enter “Cr Fe Co Ni C H F N O”
Leave other settings default.
Click “Run solution”
Superflip starts and tries to solve the structure with default options for the charge flipping
algorithm. You may let the iteration run for a few hundred cycles and you will see that there are
no signs of convergence in the three indicators (R-value, total charge and peakiness).
There is still a chance that the structure is solved – sometimes the convergence is just not
visible. We will interrupt the iteration and check the result.
In the window of Jana click the button “Interrupt”
Superflip finishes the calculation and the log-file opens in the main window of Jana. Check the
symmetry analysis. The final part will look like this:
-------------------------------------------------
Space group derived from the symmetry operations:
-------------------------------------------------
HM symbol: P1
Hall symbol: P 1
Fingerprint: 33001041a0 (0,0,0)
Symmetry operations:
1: x1 x2 x3
It is most unlikely that the structure really has symmetry P1. Together with the lack of indicators
of convergence this means that the structure is not solved.
The structure contains moderately heavy atoms. It is thus sensible to start the iteration with
Patterson superposition map instead of the random map.
Click the button “Close” in the upper row to close the viewer of the log file;
In the selection list labeled “Starting model” select the radiobutton “Patterson superposition
map”:
This time the solution has the correct symmetry and it is a strong indication that the structure is
solved.
Close the log-file viewer;
Wait until Jana performs the peak search (be patient, it is a large structure…)
Click the button “Draw 3D map” to view the resulting electron density
The density and the interpreted model are shown overlaid in VESTA. The structure is large and it
is not easy to understand the density. Nevertheless, after a while it becomes clear that the
structure is correctly solved and the main expected motifs are well visible. It is also visible that
some of the peaks present in the electron density were not interpreted as atoms by the peak-
seaching routine. These peaks are saved as “Fourier maxima” in Jana and you may change them
to regular atoms when you go to Parameters → Atoms → New → Peaks from the last Fourier
calculation. To see which maxima are interesting, we will plot the strongest peaks together with
the structure:
Close VESTA;
Click the button “Draw structure”;
In the next dialog click “Add Fourier peaks”;
A dialog with a long list of maxima opens. First about 100 maxima are deselected, because they
represent existing atoms. Then a large number of maxima are selected (typed in red).
OK to include all maxima;
The structure-plotting program (typically Diamond) opens and you may view the structure.The
added maxima are shown as the lightest element from the element list. In this case it is
hydrogen. You can see that a large number of these maxima are just artifacts close to the heavy
atoms, but among the added maxima you can also find most of the missing atoms. This
illustrates the fact that often the raw solution – the electron density – is more complete and
informative than the structure interpretation, and it is worth checking the electron density after
the solution, if the interpreted structure appears incomplete or otherwise wrong.
Here is for reference the image of the two chemically different entities found in the structure:
In this example we shall need installation of VESTA for visualization of electron density
associated with anharmonic ADP. VESTA is supplied as a zip archive which can be unpacked to
arbitrary location (with simple path – not to Desktop). Then the path to VESTA.exe must be
defined in “Parameters → Programs”.
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction – CCD”; NEXT
Select “Input from “cif_od” file”; OK
The used diffractometer software splits reflection file and another information (cell
parameters, wavelength …). There are several output formats for both cases.
Leave all settings unchanged; NEXT
Leave all settings unchanged; NEXT
The program reads 11948 reflections from hkl file
For absorption correction select “None or done before importing”; NEXT
FINISH
3. Data RepositoryOK;
Each line in data repository window corresponds to one data set converted to Jana format.
In our case we use only one data set.
OK; YES to accept the data set
4. Symmetry Wizard
[On the screen: “Tolerances for crystal system recognition”.]
Symmetry wizard can be started separately by “File → Reflection file → Make space group
test.
Leave all settings default; NEXT
[On the screen: Select Laue symmetry]
8. Refinement
[On the screen: basic window of Jana]
Right-click the icon of Refine.
[On the screen: refinement options]
Define 100 of cycles, damping factor 0.5; OK
Choose “YES+START”
Refinement converges with R value about 12%
Start “Edit atoms”
[On the screen: list of atoms]
Press “Select all” and “Action → Edit/define”
Choose “harmonic” for ADP parameters. This defines anisotropic temperature parameters
for all atoms. OK, OK, Yes to save changes
Double-click the icon Refine.
Refinement converges with R value about 4%.
9. Extinction correction
The results seem to be OK, but:
[On the screen: basic window of Jana2006]
Open listing of Refine (by “Edit/View → View of Refine”)
Press “Go to” and select “Statistics Fo, sin(th)/lambda after refinement”
R value should decrease with increasing intensity or with lowering the diffraction angle.
Here the R value for the strongest reflections as well as for the low angle reflections is higher
than in the previous shell. This suggests that extinction correction is required.
Close the refinement listing
“Parameters → Extinction”
Select “Isotropic”; OK
Refinement of extinction coefficient is enabled by default
Double-click Refine and refine the final structure
Refinement converges with R factor 3.67 %.
Check the extinction coefficient through “Parameters → Extinction”
View listing of Refine; Go to “Statistics”; check improved R value for strongest reflections.
Quit Contour
In “Atoms Edit” of Contour define the Ag1 atom and all four S1 atoms using their symmetry
codes.
Press “Run 3d maps”
Jana2006 will call an external program for drawing isosurfaces, previously defined in
“Tools → Programs”, in our case VESTA
Setting drawing options:
In “Edit → Bonds” define bond between Ag and S; maximal distance 3.15 Å,
In “Objects → Structural model” select “Stick”,
In “Properties” check “Do not show unit cell”.
In “Curves” “Use mouse to define polyline” and draw the line Ag1(a)—P—Ag1(b)—P—Ag1(c)—
P—Ag1(a) by mouse (2nd mouse button quits the selection). You can also fill in these points
directly in the textboxes of the “Curves” dialogue. Fractional coordinates of the point P are
approximately 0.267 -0.63 0.235. Ag atoms should be written with corresponding symmetry
codes, i.e. Ag1#s3c2.
Press “Continue”
Adjust “y-Scale”: Min and Max in meV: 0 600;
The resulting curve is one particle potential path. The peaks represent the potential barrier
between Ag1 atoms when the exchange path runs through the point P. These values suggest
that the silver transfer within and between neighboring Ag–S–Ag spirals is possible.
X-ray data from an oxalic acid sample measured at 110K to a high resolution limit. The crystal
contains also crystal water: C2O4H2.2H2O
Input files: oxa.m40, oxa.m50, oxa.m90 and oxa.m95 – Jana files from a conventional
refinement
For defining of multipole parameters it is convenient to select a local coordinate systems for
each atom. The following figure shows a recommended choice for atoms C1, O1, O2 and H1:
Using the list box on the top left corner go to the next atom – O1.
Change multipole parameters to order 3 and define the coordinate system x → C1 and
y → O2. Go to the next atom – O2.
Change multipole parameters to order 3 and define the coordinate system x → C1 and
y → O1. Go to the next atom – O3.
Change multipole parameters to order 3 and define the coordinate system x → H2 and
y → H3. Go to the next atom – H1.
Change multipole parameters to order 1 and define the coordinate system x → O1 and
y → C1. Go to the next atom – H2.
Change multipole parameters to order 1 and define the coordinate system x → O3 and
y → H3. Go to the next atom – H3.
Change multipole parameters to order 1 and define the coordinate system x → O3 and
y → H2.
OK; OK; Yes to the question “Do you want rewrite changed files?”
[On the screen: Basic Jana window]
Right click on the icon “Refine”;
In the page “Basic” select the number of cycles 100 and damping factor 0.1;
In the page “Various” press “Fixed commands” and fix the following parameters:
385 Jana2006 Cookbook, version October 2017
“Set individual”, P00[*], “Set to” 0;
“Individual”, kappa'[*];
[*] stands for all atoms. P000 wil be zero for all atoms while kappa’ will be fixed to the
current value.
OK; OK; YES+Start
Refinement converges to R(obs)=1.75%.
The map does not show any clear positive or negative maxima which means that the model
explain the experimental density adequately.
Fractional coordinates
Bond ρ Laplacian Hess eigenvalues Ellipsity
x y z
Number of
Atom Charge
electrons
C1 5.147 0.853
O1 9.124 -1.124
O2 8.773 -0.773
O3 9.378 -1.378
H1 0.396 0.604
H2 0.407 0.593
H3 0.419 0.581
Sr1.272CoO3
Revised: 28 April 2014
Single crystal data measured with Oxford Diffraction four-circle diffractometer
Input files: bezcox.hkl, bezcox.cif_od
Frame scaling, absorption correction: done by diffractometer software
2. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Oxford Diffraction – CCD”; NEXT
[On the screen: “Define basic input file:"]
Select “Input from cif_od file”; file name “bezcox.cif_od”; OK
[On the screen: “Complete/Correct experimental parameters”]
Leave all settings unchanged; NEXT
[On the screen: “Define the reference cell/split by twinning”]
Leave all settings unchanged; NEXT
For absorption correction select “None or done before importing”; NEXT
FINISH
3. Data Repository
OK; YES to accept the data set
[On the screen: “Information about starting symmetry wizard”]
NEXT
4. Symmetry Wizard
[On the screen: “Tolerances for crystal system recognition”.]
Symmetry wizard can be started separately by “File → Reflection file → Make space group
test”.
Leave all settings default; NEXT
[On the screen: Select Laue symmetry]
Select Trigonal -3m1 Laue symmetry; NEXT
[On the screen: Select cell centering]
In the list we can see that R-reverse centering is almost fulfilled.
Click on “Details” and examine strongest reflections contradicting R-reverse centering
The most of printed reflections follow the rule -h+k+l=3n. Similarly, reflections contradicting
R-obverse centering follow the rule h-k+l=3n. This can be a sign of obverse/reverse twinning.
Return back to the page with tolerances and check there “Introduce twin laws in case of
subgroup”
5. Define W-matrix
Using Tools → Reciprocal space viewer we can clearly distinguish the composite subsystems.
This compound has the simplest W matrix, which only exchanges the “l” and “m” index.
Than means that “h k l m” are reflections of the first subsystem while “h k m l” are
reflections of the second subsystem.
Open “Edit M50”, change number of composite parts to 2 and define W-matrix for the
second subsystem: (1 0 0 0 | 0 1 0 0 | 0 0 0 1 | 0 0 1 0 )
“Induced cell paramemeters” are parameters of the second subsystem, as follow from the W
matrix
6. Structure Solution
Double-click the icon “Structure solution”.
[On the screen: window of Structure solution]
In “Formula” textbox type list of chemical elements for Sr1.272CoO3 : Sr Co O
Select “use Superflip”; “Peaks from Jana2006 but first run Fourier”
Select “Use a specific random seed” and type value “111”
Press “Run solution”;
The option “Peak from Jana2006 but first run Fourier” in necessary for proper interpretation
of Fourier map for composite structure because we need to look for bot subsystems
Superflip converges (after noise suppression) with R value 20%. It confirms the symmetry
R-3m(00g)0s.
Fixing the random seed to 111 guarantees for most processors that the Superflip returns
exactly the same results like described in this cookbook. For this example finding proper
solution would require many attempts with Superflip.
Press “Quit” to leave the structure solution wizard
Start “Parameters → Twin fractions” and change the twin fraction to 0.5
Activate refinement of twin fraction
In Refinement options change number of cycles to 100
Run refinement
Refinement converges with R values 34%, 24%, 41%, 58%, 67% and 97%
Start “Edit/View → View of Refine” and read information following the header “Composite
part #2”.
Here we can see unit cell and symmetry of the second composite subsystem. Note that the
second subsystem has also different centering because centering vectors are also
transformed by W matrix.
9. Oxygen
Run refinement to obtain input for Fourier
Run difference Fourier for composite part No 1, with Scope=”automatically”, default map
orientation, independent parallelepiped
We use composite part 1, because we expect oxygen belonging to Co coordination
YES to start procedure for including atoms
Locate the strongest maximum with reasonable distance to Co (probably Max2)
Include this maximum as O1 in position about 0.17, 0.04, 0.5
Start Contour and plot section x3-x4 through O1, scope 1 1 2, difference Fourier, using
summation of x1, x2 over 1 Å;
394 Jana2006 Cookbook, version October 2017
section x3-x4 of O1
The section reveals O1 should be described with crenel function, x40=0.75, delta=0.5
Plot section x1-x2, difference Fourier, x4 from 0 to 1, scope 5 5 0, step 0.5A, center 0 0 0.5
The section x1-x2 with center at (0,0,1/2) with scope about 5 Å as a function x4 shows the
character of the modulation
The largest different peaks are around Sr1 position. Now we shall attempt to describe
modulation function of Sr1 with Lagrange polynomials, with crenel defined at x40=0 with
delta=1
The neighborhood of Co1 is changing from octahedron to trigonal prisma (check it with
drawing program). The trigonal coordination induces a slightly triangular shape of ADP.
However, the region in x4 in which the effect is detectable is very narrow and the number of
modulation waves for ADP’s of the order 3 would be very large to describe such abrupt
change.
The atom Co1 should therefore be split into two positions by crenel – the first one Co1
(x40=0, delta=0.45) and the second one Co1’ (x40=0.25, delta 0.05). After adjusting of
occupancies to ai=1/6 for both atoms the two Legendre waves for position and harmonic
ADP Co1 should be used.
For the Co1’ only average parameters for position and ADP up to the third order.
Ba5Co5ClO13
Powder data measured with the ILL neutron source at 1.4 K
For more details: Mentré, O., Kauffmann, M., Ehora, G., Daviero-Minaud, S., Abraham F. &
Roussel, P. (2008). Solid State Sciences 10, 471-475.
In order to facilitate the work with magnetic structures, Jana2006 offers special wizards and
protects the basic file - the parent structure. However, for experienced users of Jana, all
operations offered by wizards can be also done manually by starting the tools directly from the
basic user interface. The standard work with magnetic structures based on powder data follows
this scheme:
Close;
Select one of the hexagonal magnetic space groups; NEXT
This will pass all present hexagonal magnetic space groups to further testing.
[On the screen: Select Shubnikov space group]
For an individual Shubnikov space group the tool shows components of the magnetic vector
of all magnetic atoms in the structure. The button “Show details” gives a more detailed
406 Jana2006 Cookbook, version October 2017
analysis. In our case only two Shubnikov space groups – P63/mm'c' and P63' /m'm'c – allow
magnetic ordering for all Co atoms. Note that for the first of them the global magnetic
moment for is generally non-zero, which means that this arrangement allows a
ferromagnetic ordering. As ferromagnetism has not been detected for this compound, the
first model can be excluded from the analysis. However, we shall test all models to check if
such a decision can also be based on the diffraction experiment.
Select the Shubnikov space group P63/mm'c' and press “Start profile simulation”. Use "Go
to 2th/hkl" textbox and find position of the reflection 2 -1 3 . Click on the tick
corresponding to the reflection 2 -1 3:
Although the peak is clearly visible, the calculated structure factor is zero.
ESC
Press “Quit” and the select Shubnikov space groups P63'/m'm'c and investigate again the
reflection 2 -1 3:
For this Shubnikov space group the intensity of (2 -1 3) is non-zero in accordance with the
diffraction pattern. As this reflection is systematically absent for the nuclear structure (c
glide plane), it is clear that it originates from magnetic scattering. As apparent from the
Shubnikov space group test all magnetic moments are directed along the z-axis. This leads
to the systematic extinctions l = 2n + 1 for reflections (2h,− h,− h, l ) but only in the case that
the glide plane ”c” is combined with the time inversion operation. The fact that the peak 2 -1
3 is clearly detectable gives us another support for the Shubnikov space group P63'/m'm'c.
Select again the Shubnikov space group P63/mm'c' and press “Start graphic simulation”;
Select “Use Vesta to draw the magnetic structure”;
Press "Draw+return" or "Draw+continue";
In the case that bond between atoms are created go to “Edit→Bonds…” and delete all
automatically generated bond limits
If necessary, you can modify the length of the displayed vectors by “Edit→Vectors…”
Click “Restrictions” and make sure ADP parameters all atoms are restricted to be the same:
The “Fixed commands” and “Restrictions” are created automatically by the program when
preparing the testing structure.
OK; select the page "Basic"
Check if the “Number of cycles” is 100 and the “Damping factor” is 0.1;
“Use Marquart technique” check that the Fudge factor is 0.001;
410 Jana2006 Cookbook, version October 2017
OK; Start refinement;
Refinement converges with Rp~3.3%, R(obs)~8.5%.
Click “Restrictions” and make sure ADP parameters all atoms are restricted to be the same:
Ba6Co6ClO16
Powder data measured with the ILL neutron source at 1.4 K
For more details: Mentré, O., Kauffmann, M., Ehora, G., Daviero-Minaud, S., Abraham F. &
Roussel, P. (2008). Solid State Sciences 10, 471-475.
2. Import Wizard
Select “Magnetic: nuclear model from CIF”; NEXT
Right pane: select the input file Ba6Co6_nuclear.cif; OK
Press “Quit” and repeat this procedure for the next magnetic Shubnikov group P[c]-6c2 with
origin shifted by (0,0,1/4):
The first structure model has better R values. Moreover, the presence of magnetic moments
at the tetrahedral site corresponds to the result in Example 12.1
2. Import Wizard
Select “Magnetic: nuclear model from CIF”; NEXT
[On the screen: Select input CIF file]
Right pane: select the input file PrSrMnO_basic.cif; OK
The first peak has a strong asymmetry. For neutron and synchrotron cases, such asymmetry
can be described by the correction for divergence as introduced by Finger-Cox-Jephcoat, J.
Appl. Cryst. (1994). 27, 892-900.
“Quit”
“Edit profile parameters”
[On the screen: Powder options]
In the page “Asymmetry” select “correction by divergence”
NEXT
[On the screen: Select Shubnikov space group]
OK; “Yes+start”;
425 Jana2006 Cookbook, version October 2017
Refinement converges with Rp~10.1%, R(obs)~3.7%.
This refinement is very unstable.
Double click on the icon “Plot structure” and select Vesta to draw the structure:
Use “Tools → Powder → Profile viewer” and check fit of the calculated powder profile:
8. Refinement in the Shubnikov superspace group C[B]mcm with the origin shift (0,1/4,1/4)
(№ 4 in the list)
[On screen: Basic window]
“File → Structure → History” and select the parent structure PrSrMnO; OK;
[On screen: Basic window]
“Tools → Special tools → Representation analysis for magnetic structures”;
[On screen: Representation analysis]
NEXT; NEXT; Select the fourth Shubnikov superspace group “C[B]mcm+(0,1/4,1/4)”; FINISH
YES to the question “Do you want to continue with the last kernel/epikernel?”;
[On screen: Select structure name]
Use the default name i.e. PrSrMnO_02; OK
No to the question “Do you want to apply this transformation”.
Right click on the icon “Refine”; Select the “Basic” and change the “Number of cycles” to
100 and the “Damping factor” to 0.1;
Select the page “Various” and use the button “Restrictions”
Check that all atoms are restricted to have identical ADP parameters and not restricted
occupancies;
Restrict atoms Pr1 and Sr1, occupying the same position in the nuclear structure: identical
Coordinates ADP parameters and “keep overall sum” for occupancies; Add;
OK; OK; “Yes+start”;
Refinement converges with Rp~6.0%, R(obs)~3.2%.
Double click on the icon “Plot structure” and select Vesta to draw the structure::
Use “Tools → Powder → Profile” viewer and check fit of the calculated powder profile:
9. Refinement in the Shubnikov superspace group C[B]mcm with the origin shift (0,1/4,1/4)
(№ 6 in the list)
The symbol in the table is the same as for the № 4 but both epikernels are generally
different. Compare details for both. The first one has in the original setting the symbol
A[C]mma and the second one A[C]mam. Both give the same symbol when transformed to
the standard setting.
Make the same steps and calculation as in the previous test.
2. Import Wizard
Select “Magnetic: nuclear model from CIF”, NEXT
Right pane: select the input file MnWO4_nuclear.cif
NEXT
[On the screen: List of kernels and epikernels]
Notice that for both SSG’s there are non-zero amplitudes for sine and cosine waves.
However only three amplitudes are independent as can be found from “Details”:
The cosine amplitudes have generally non-zero values but they are related with the sine
amplitudes. However the coefficients are functions of y coordinate of the Mn1 atom, which
means that these restrictions are non-linear. Later we shall see how these restrictions can be
simplified.
In this case, both possibilities yield similar diffraction patterns and we cannot decide at this
stage which solution makes a better fit.
Note: The presence of a rational non-zero modulation vector component qy = 0.5 leads to a
specific problem of the origin selection. While for the space group P 2 c the c glide planes
The second solution gives significantly better profile fit and also better Bragg R values.
Quit
Now the y component of the sine wave and x and z components of the cos wave are fixed to
zero values. This setting makes a full separation of internal and external space (before the
transformation some wave components depended on the y coordinate).
2. Import Wizard
Select “Magnetic: nuclear model from CIF”, NEXT
Right pane: select the input file HoNi_nuclear.cif
Select the only space group and press ‘Continue with the selected Shubnikov space group’;
[On the screen: Select structure name]
Use the default name "HoNi_02"; OK
442 Jana2006 Cookbook, version October 2017
NO to the question “Do you want to apply this transformation?”;
[On the screen: Basic window]
Start refinement;
Refinement converges with R(obs)~ 12.8%, 12.3%, 14.3% for all, nuclear and magnetic
reflections, respectively.
[On the screen: Basic window]
“Parameters → Extinction”
[On the screen: Extinction model]
Start “Parameters → Extinction”
Select Isotropic, Type1, Gaussian; OK; Yes to rewrite changes
[On the screen: Basic window]
Run refinement
Refinement converges with R(obs)~ 9.6%, 7.8%, 14.6% for all, nuclear and magnetic
reflections, respectively.
Right click on the icon “Refine”;
Select the page “Various” and use the button “Restrictions”;
Disable the command “restrict * 12” which keeps ADP's of all atoms identical. OK;
Use the button “Fixed” and disable the command “fixed xyz *” which fixes coordinates of all
atoms. OK; OK; Yes+start;
Refinement converges with R(obs)~ 8.4%, 6.7%, 13.0% for all, nuclear and magnetic
reflections, respectively.
Start “Parameters → Twin fractions”;
[On the screen: Edit twin fraction]
Selection of the orthorhombic Shubnikov space group represents a lowering symmetry from
the original one (parent symmetry. This can create different magnetic domains related by
the lost symmetry operations. The lost symmetry operation has been automatically created
by the program as a twinning operation. The default value of the relative fraction was 0.5
which means 1:1 ratio. In the final stage this value can be refined.
Press the button “Show twinning matrix”.
[On the screen: Twinning matrix #2]
As expected the twinning operation is the four-fold axis, which was present in the tetragonal
description but it is missing in the orthorhombic description.
OK;
[On the screen: Edit twin fraction]
Change the volume fraction to zero and activate the refinement key for the volume fraction.
OK; Yes to the question “Do you want to rewrite changed files?”.
[On the screen: Basic window]
Run refinement
Refinement converges with R(obs)~ 8.3%, 6.7%, 12.6% for all, nuclear and magnetic
reflections, respectively.
[On the screen: Basic window]
Start “Parameters → Twin fractions”;
[On the screen: Edit twin fraction]
Check that the volume fraction return to the value close to 0.41.
[On the screen: Basic window]
Double click on the icon “Plot structure”; Use the Vesta program to draw the structure;
Press "Draw+return" or "Draw+continue";
443 Jana2006 Cookbook, version October 2017
In Vesta enlarge the boundaries in z direction (like in the previous plot):
Press Change;
Press OK;
Select irreducible representation
[On the screen: ISODISTORT irreducible representation ]
Here we can see the kernel and the set of epikernels (isotropy subgroups) associated with
mGM5+.
Return one step back in your Internet browser and select the representation mGM2+; OK;
[On the screen: ISODISTORT: order parameter direction]
Return back to the page ISODISTORT irreducible representation and select now mGM4+
Repeat the work as for GM2+ and save the CIF file as a Text Document with name GM4+
into the working directory.
Return to Jana2006
DyMn6Ge6
Neutron powder data measured at 11K on DMC at PSI (Switzerland)
Available at the Bilbao crystallographic server, database magndata, #1.1.10
References:
J. Rodríguez-Carvajal & F. Bourée, EPJ Web of Conferences 22, 00010 (2012)
P. Schobinger-Papamantellos et al, Journal of Alloys and Compounds, 203 (1994) 243
2. Import Wizard
Select “Magnetic parent structure: nuclear model from CIF”; NEXT
Right pane: select the input file dymn6ge6_nuclear.cif; OK
This window has an information character. You can see here all irreps leading to magnetic
ordering. The “details” buttons provide more information about an individual irrep and its
connection to the kernel magnetic symmetry.
NEXT
Select one of the hexagonal Shubnikov superspace groups based on average space group
P6/mmm; NEXT;
This will pass all present hexagonal Shubnikov superspace groups of type P6/mmm to
further testing. We assume that both magnetic atoms should have non-zero magnetic
moment.
[On the screen: Select Shubnikov space group]
Only the second Shubnikov superspace group allows for non-zero magnetic moments for Dy
atom.
Select the second Shubnikov superspace group: P6/mmm1’(00g)00sss and press “Start
profile simulation”;
456 Jana2006 Cookbook, version October 2017
Press "X exactly" and focus the region 8-16 degrees; Press “Fit Y”:
Two dominant magnetic peaks in this region have zero intensity. This means that the tested
model does not describe the magnetic ordering correctly. However, in making such
conclusion we have to take into account that the simulation procedure uses randomly
generated magnetic moments and that the result can be affected by an actual random
choice. This is why we recommend repeating the simulation when some doubt exists.
“Quit”; “Back” and return to the “List of kernels and epikernels” and select one of
hexagonal superspace groups based on P622.
While the first one P622.1’(00g)t00s restricts the magnetic moments of the Dy atom to zero,
the second one P622.1’(00g)-h00s allows non-zero magnetic moments for both magnetic
atoms. The simulation profiles demonstrate that the intensity of both magnetic peaks in the
interval (8-16) degs are generally non-zero:
There are several maximal subgroups without this operation. Then can be generated by the
Jana tool “Tools→Transformations→Go to subgroup structure”. As generators for the
working group until now can be used:
{6+|0 0 0 5/6}: x1-x2, x1, x3, x4+5/6
{2_xy| 0 0 0 0}: x2, x1, -x3, -x4
{1’|0001/2} : x1, x2, x3 , x4+1/2
459 Jana2006 Cookbook, version October 2017
The possible maximal subgroups will then be generated by:
P622(00g)-h00
{6+|0 0 0 5/6}: x1-x2, x1, x3, x4+5/6, m
{2_xy| 0 0 0 0}: x2, x1, -x3, -x4, m
P6’22’(00g)t00
{6+’|0 0 0 1/3}: x1-x2, x1, x3, x4+1/3, -m
{2_xy| 0 0 0 0}: x2, x1, -x3, -x4, m
P62’2’(00g)-h00
{6+|0 0 0 5/6}: x1-x2, x1, x3, x4+5/6, m
{2’_xy| 0 0 0 1/2}: x2, x1, -x3, -x4+1/2, -m
P6’2’2(00g)t00
{6+’|0 0 0 1/3}: x1-x2, x1, x3, x4+1/3, -m
{2’_xy| 0 0 0 1/2}: x2, x1, -x3, -x4+1/2, -m
Make sure that your actual structure is the final refined structure from the previous step.
x1-x2 x1 x3 x4+5/6 m
x2 x1 -x3 -x4 m
NEXT; NEXT
[On the screen: Select one representative of the 2nd coset]
Select the first operation.
Generally, all selections would give the same result.
NEXT;
[On the screen: Specify the output structure]
Select “dymn6ge6_01_subgroup_1” or any other
NEXT;FINISH;”Yes” to the question “Do you want to continue with the new structure?”
[On the screen: Basic window]
Run “Edit Atoms”; Select atoms Dy1 and Mn1; Press the button “Action”
Select “Reset modulation/magnetic parameters”
Deselect “Magnetic parameters k<>0”
Choose “Reset parameters randomly in the interval (-pmax,pmax)
OK;OK;”Yes” to the question “Do you want to rewrite changed files?”
Run “Refine”;
Repeat the procedure for all subgroups listed above
Use the “File→ Structure → History” to return to the final refined structure in
P622.1'(00g)-h00s.
P622(00g)-h00
Rp=4.54%, R(obs)=4.75%,4.35%,5.47% for all, nuclear and magnetic reflections, respectively.
P6’22’(00g)t00
Rp=4.42%, R(obs)=4.46%,3.96%,5.38% for all, nuclear and magnetic reflections, respectively.
461 Jana2006 Cookbook, version October 2017
P62’2’(00g)-h00
Rp=3.84%, R(obs)=3.11%,2.44%,4.35% for all, nuclear and magnetic reflections, respectively.
P6’2’2(00g)t00
Rp=4.45%, R(obs)=3.91%,4.35%,5.40% for all, nuclear and magnetic reflections, respectively.
One of the symmetries clearly fits the data better. Observe that the strong misfit of some of
the main reflections has disappeared.
OK
Press button “Properties…”; Go the page “Atoms” and modify “Radius” for all atoms
The result as derived step by step from the representation theory is expressed in a different
chiral domain than used in MAGNDATA. Therefore, we should transfer it to the equivalent
chiral domain.
[On the screen: Basic window]
“Tools → Transformation → Change enantiomorph”
Create a new structure with name e.g. “dymn6ge6_01_subgroup_3_chiral”
OK; Yes to continue with the new structure
Run EdiM50; Go the page “Symmetry”
Now the superspace group is the same as used in MAGNDATA however there is a non-zero
origin shift in x4 equal 1/6.
ESC
“Tools → Transformation → Origin shift”
Use shift to the point 0 0 0 -1/6
OK
The last point that has to be checked for a full comparison is if we are using in the same
symmetry positions for the magnetic atoms Dy1 and Mn1. The positions of these atoms
from 1.1.10 are:
466 Jana2006 Cookbook, version October 2017
Dy1 (0,0,0) and Mn1 (0.5,0,0.251).
[On the screen: Basic window]
Run “Edit atoms”; Double click on Dy1 atom; Go to the page “Edit”
The coordinates are (1,1,1) and this means that we can just rewrite them to (0,0,0) as such
transformation does not change the modulation functions.
Move to the atom Mn1
Its coordinates are (0.5,1,0.749). From the list of symmetry operations we can find that the
symmetry operation: x1-x2,-x2,-x3,-x4+1/2 combined with the lattice translation (1,1,1,0)
makes what is needed.
ESC;ESC
[On the screen: Select atoms to be used]
Select the atom Mn1; Press button “Action”; Select “Transform selected atoms”
Use the button “Fill by a symmetry operation”; Select x1-x2 -x2 -x3 -x4+1/2 –m; Modify the
translation vector to 1 1 1 1/2
OK; OK; Yes to rewrite files
[On the screen: Basic window]
Run “Refine” to check that all transformation did not change the solution
[On the screen: Basic window]
“Edit/View → View of Refine”; Press “Go to” and select “Changes overview”
13. Use of the option of polar coordinates for the moments to constraint further the Mn spin
to have a circular modulation.
[On the screen: Basic window]
“File → Structure → Save as”; Use e.g. “dymn6ge6_01_subgroup_3_final”
OK; No to the question “Do you want to continue with the new structure?”
[On the screen: Basic window]
“Parameters → Magnetic moments in polar coordinates”;
Select both atoms
OK; Yes to rewrite changed files
Run “Edit atoms”
[On the screen: Select atoms to be used]
Double click on the atom Dy1; Go the the “Edit” page;
The windows for the magnetic moment parameters have changed to polar coordinates, but
the control of the parameters to be refined does not change automatically, and must be
done manually: for the Mo parameters of Dy1, fix the two polar angles, and leave to refine
the modulus (the first one). For the modulation parameters, modify the refinement controls
if necessary. Only Mrsin1 parameters should be free, the rest are forced by symmetry
namely Mrcos1=Mrsin1
For Mn1, the modulation parameters except the theta angle are symmetry free. In order to
force a circular modulation, the constraints must be used: Mrcos1=Mrsin1, Mpcos1=Mpsin1-
90. These equations are only roughly satisfied by the previous refined model, as shown after
the transformation.
Right click on the “Refine” icon; Select the page Various; Click on the button “Equations”;
Define the equations mentioned above:
Mp0[Dy1]=0
Mt0[Dy1]=0
Mrcos1[Dy1]=Mrsin1[Dy1]
Mpsin1[Dy1]=90
Mpcos1[Dy1]=0
Mtsin1[Dy1]=90
Mtcos1[Dy1]=90
Mp0[Mn1]=0
Mt0[Mn1]=180
Mrcos1[Mn1]=Mrsin1[Mn1]
Mpcos1[Mn1]= Mpsin1[Mn1]-90
Mtsin1[Mn1]=90
OK;OK; Yes+start
Rp=3.83%, R(obs)=2.93%,2.43%,3.82% for all, nuclear and magnetic reflections, respectively.
Which means that it is not significantly worse than the full refinement – see point 9.
cos sin
Mx My Mz Mx My Mz
MAGNDATA 6.20 0.0 0.0 3.58 7.15 0.0
Dy1
Jana 6.12 0.0 0.0 3.53(4) 7.07 0.0
MAGNDATA -1.88 -0.21 0.0 -0.85 -2.05 0.0
Mn1
Jana -1.76(10) -0.19(9) 0.0 -0.79(9) -1.92(8) 0.0
General note: Two different strategies for handling refinement keys in Jana2006
From Petříček, V., Dušek, M. & Palatinus, L. (2014) Z.Kristallogr. 229, 345-352.
“All refinable parameters are divided into two groups: automatic and user-controlled ones. The
automatic parameters are those which can be usually refined without numerical problems: scale
factors and extinction parameters; parameters describing position, ADP and their modulation
for individual atoms as well as for molecules; background parameters of the powder profile;
magnetic parameters. The user-controlled parameters are twin/phase volume fractions, atomic
or molecular site occupancies, parameters of crenel modulation functions, parameters of
multipole refinement and powder profile parameters for individual Bragg reflections.
The automatic parameters are activated by default at the beginning of the refinement process
while the user-controlled ones must be always explicitly activated by the user. The program
always fixes the parameters which cannot be refined due to symmetry restrictions. In order to
fix other automatic parameters the user is provided with a range of special fixing commands.
This way is very useful for most cases because it minimizes mistakes in setting refinement flags.
However, the program also allows switching to a completely manual mode, where all
parameters are to be activated explicitly by the user, analogically to other refinement programs
such as FullProf [29].”
References:
Zhu et al., PRL 113, 076406 (2014)
PHYSICAL REVIEW B 92, 094419 (2015)
Propagation vector:
k= (0,0,0)
Magnetic site:
Cr 4e (0,0,z)
2. Import Wizard
Select “Magnetic parent structure: nuclear model from CIF”; NEXT
Right pane: select the input file cr2wo6_nuclear.cif; OK
This window has an information character. You can see here all irreps leading to magnetic
ordering. The “details” buttons provide more information about an individual irrep and its
connection to the kernel magnetic symmetry.
NEXT
Select the second Shubnikov superspace group: P42’/mnm’ and press “Start profile
simulation”;
Press "X exactly" and focus the region 12-33 degrees; Press “Fit Y”:
The second Shubnikov space group has the global magnetic moment zero.
Select the first Shubnikov superspace group: Cm’mm’ and press “Start profile simulation”;
Press "X exactly" and focus the region 12-33 degrees; Press “Fit Y”;
Both reflections (1,0,0) and (0,0,1) have generally non-zero intensity and this means that we
cannot reject this space groups only from violation of systematically absent reflection as for
the tetragonal case. This is also true for Cmm’m.
While the P centered cell are identical with that used for the parent tetragonal cell the C
centered cell is related to the original cell by equation:
The orthogonal ordering of magnetic moments can have some influence to the cell
parameters a and b which need not be same anymore. In order to test orthorhombic strain
without mixing it with the refined magnetic structure we can use le Bail profile fitting.
The cell parameters are: a=6.4849(23) and b=6.4850(23) which means that there is no strain
within the accuracy following from the refinement.
For good quality data set, the white curve contains distinct steps, and, consequently, the grey
curve has distinct peaks.
Press Enter
In the next step an autoconvolution of the diffraction pattern (difference space) is analyzed,
and the groups of peaks in the autoconvolution (clusters) are replaced by the cluster centers.
Again, a distance distribution is displayed, and clear-cut steps on the curve indicate a well-
defined lattice.
Press Enter
Now the first step is finished and PETS prepared files needed for unit cell indexing. Keep the
window of PETS open, you will need it later!
IMPORTANT: before you continue, make sure you have the same orientation of the lattice on
your computer as displayed on the picture above (check the direction of the axes on the right).
If not, you may obtain slightly different results in the rest of the tutorial. While not necessarily
wrong, it makes following this tutorial more difficult.
Click the button in the lower-left corner
Click, hold the mouse button, and drag the marker to obtain the result on this image:
7. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Pets electron diffractometer”
The file name automatically changes to ttb.cif_pets, which is a file produced by PETS
containing all important information, including the list of intensities.
NEXT
Leave all settings unchanged (note the wavelength of 120kV electrons!); NEXT
Leave all settings unchanged; NEXT
The program reads 3777 reflections from the file
OK
For absorption correction select “None or done before importing”; NEXT
FINISH; OK; YES to accept the data set
You just read in the cell parameters and the intensity list. If you do not have the data in the
cif_pets format, you can read them in from a general hkl file using the option “reflection file
corrected for LP and absorption” at the beginning of the import wizard. You then have to input
the radiation type, wavelength and cell parameters by hand.
8. Symmetry Wizard
NEXT to close the information window and start the wizard
Symmetry Wizard can be started separately by “File → Reflection file → Make space group
The direction and length of the b-axis is difficult to determine. The lattice is obscured by quite
a lot of noise. This noise comes from the diffuse nature of the reflections along b* and from
the mosaicity of the crystal. If such noise is present, it is amplified by the data processing
preceding the creation of the file .clust. In such cases it is useful to check the correctness of
Quit the indexing plugin to create file kankite.smr with the orienting matrix, and continue
data processing in PETS by typing u[Enter] to read in the orienting matrix, and i[Enter] to
integrate the intensities.
You will note significant overlap of integration boxes (white circles) in the images displayed
by ImageJ. This is, however, not a real problem, because the lattice is C-centered, and every
second reflection is systematically absent.
PETS produced file kankite.cif_pets, which will be imported in Jana2006 in the next step.
Type h[Enter] and y[Enter] ]to produce standard sections through reciprocal space.
With this step the data processing in PETS is finished.
Type q[Enter] to close PETS. You may also close ImageJ.
3. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Pets electron diffractometer”
The file name automatically changes to kankite.cif_pets, which is a file produced by PETS
containing all important information, including the list of intensities.
NEXT
Leave all settings unchanged; NEXT
Leave all settings unchanged; NEXT
The program reads 5921 reflections from the file
OK
For absorption correction select “None or done before importing”; NEXT FINISH; OK; YES to
accept the data set
You just read in the experimental parameters, cell parameters and the intensity list.
4. Symmetry Wizard
NEXT to close the information window and start the wizard
9. Structure refinement
[On the screen: basic window of Jana]
Right-click the icon of Refine.
[On the screen: refinement options]
Set Number of cycles to 100; OK Choose “Yes+start”
Refinement converges to an R value about 19%
10. Adding the missing atom[On the screen: basic window of Jana]
Right click the icon Fourier
On the first card, select F(obs)-F(calc) – difference Fourier
Modify the setting on the second card to agree with this screenshot:
Background information:
Orthopyroxene (opx) is a rockforming mineral – iron-magnesium silicate (Fe,Mg)2Si2O6,
Pbca, a=18.268, b= 8.868, c=5.202. It contains two symmetry-independent octahedral sites
with Mg and Fe. The distribution of Mg and Fe in these two sites depends on the
crystallization conditions, and it can therefore be used to infer information about the
formation of the rock.
Orthopyroxene often forms only microcrystalline aggregates or intergrowths with other
minerals. It is thus desirable to determine the occupancies of Fe and Mg on the two sites
from individual micro- or nano-crystals of opx.
Electron diffraction data were measured on FEI Tecnai G2 20 operated at 200 kV with
precession device DigiStar.
Input data: opx_001.tif: the experimental diffraction pattern, zone axis [001] of
orthopyroxene.
opx_001.pts: input file for PETS with experimental information
opx_xray.cif: cif file with the model of the opx structure refined against x-ray
data set.
opx_001_dyn_start.cif_pets: result of data processing with PETS. This file is
provided so that this example can be solved without the data processing part.
If you get lost, open the original diffraction pattern in ImageJ and compare the distribution
of reflections with the peak positions in Jana2006.
Once the two vectors are defined, the selected 2D cell can be refined:
“Refine cell → Refine UB+cell”
You should obtain a result very close to the following 2D cell:
OK; Quit
Now you have an orientation matrix that you can read back to PETS, and extract intensities.
[On screen: window with PETS, main menu]
Type u[Enter]
PETS read in the orientation matrix calculated by Jana2006.
Type d[Enter]
D stands for “integration for dynamical refinement”. PETS extracts all intensities up to
d*=1.4, and produces file opx_001_dyn.cif_pets.
This dialog contains experimental information relevant for the dynamical refinement, and
information about the zone axis (or several zone axes, if there are more zone axes
imported). Most options are now blank, because the reflection file for dynamical refinement
has not yet been prepared.
Check the box next to “Use dynamic approach”
OK to create refinement reflection file
Follow the wizard and accept default choices
[On screen: Electron diffraction dialog]
Now we have to find initial values for the sample thickness, scale factor and orientation. As
a first quick check, we will calculate only the thickness curve, using fixed orientation. This
This plot shows a clearly defined minimum at the thickness 400 Å. This together with the low
R-value shows that the starting point for the refinement is reasonable.
Quit the display of the thickness plot
Click “except of scale optimize also: orientation Run optimization”
The orientation of the experimental pattern is never perfect. Dyngo finds the parameters of
the misalignment of the zone axis from the perfect orientation predicted by the orientation
matrix. This calculation requires repeated evaluations of reflection intensities and may take
some time. You may follow the progress of the simplex algorithm in the window of Dyngo.
After the optimization, the zone axis parameters look like this:
The R-value decreased only slightly to 7.9963%. Value 0.15 of EDtheta means that the
diffraction pattern was recorded at an angle of 0.15 degrees away from the perfect
504 Jana2006 Cookbook, version October 2017
orientation. EDPhi is the azimuth of the tilt axis. Given the large precession angle, this value
is within the precision of the method, and means that the pattern was essentially perfectly
oriented.
Quit the Electron diffraction dialog and confirm changes in the m42 file.
Input data: si_ZAPs.cif_pets: a cif_pets file containing the orientation matrix, positional angles
and intensity information from three zone-axis patterns collected on the same sample – a
wedge-shaped silicon single crystal. The file was obtained from raw diffraction patterns using
the procedure analogical to example 13.1.
1. Import data to Jana2006
Start Jana2006
“File → Structure → New”
In the left pane, locate the example directory
In the right pane, double click green label si_ZAPs
Select “Single crystal: known diffractometer format”; NEXT
Select “Pets – electron diffractometer”; NEXT; NEXT; NEXT; OK; NEXT; FINISH; OK; Yes
You have just imported the data from cif_pets using default options.
As a next step Jana starts the symmetry determination wizard. The data set we have is not
really complete, but we may nevertheless attempt to proceed with the wizard:
Next to start the wizard
On the next screen, set up the tolerances according to this screenshot:
Next; OK;
You now look at a table of Rint values for various Laue classes. The cubic Laue classes have
acceptable Rint values, and there is little doubt the structure is indeed cubic.
Select Laue class m-3m; Next
The data were already integrated assuming F-centering, therefore there are zero reflections
matching the F-centering absences in the data set.
Select centering F; Next
On the next screen a list of possible space groups is presented. For the space group Fd-3m, 7
out of 24 forbidden reflections have I>3σ(I). However, looking at the details (click the button
“Details”) you will find that these reflections are in fact only two - 0 0 2 and 0 0 6. The d-
glide is thus likely to be present. Indeed, the space group of silicon is Fd-3m.
Select “Fd-3m”; Next; Finish
In the next step we prepare a “standard” reflection file suitable for structure solution.
Next; OK; OK
[On screen: the window with three radio buttons for the choice of the handling of scales]
Jana treats can data from each zone axis on a separate scale. For the structure solution and
kinematical refinement, we will, however, assume that the data are on the same scale. The
number of common reflections is insufficient to determine a reliable common scale.
Open the refinement listing (Edit → View of Refine) and check the statistics.
The displacement parameter of silicon refines to -0.0106.
3. Origin shift
The default origin of the space group Fd-3m is at the point -43m, which does not contain
inversion center. For the dynamical refinement it is desirable to have the origin in the
inversion center, if possible. We will now shift the origin of the unit cell to the inversion
center:
[On screen: Main window of Jana2006]
“Tools → Transformations → Origin shift”
Type “1/8” in all three fields; OK
4. Initiation of dynamical refinement
[On screen: Main window of Jana2006]
The displacement parameter of silicon did not refine to physically meaningful value. Before
initiating the dynamical refinement it is good to reset it to a more correct value. It is not
strictly necessary in this example, but it is generally a good practice.
Double click the icon Edit atoms
Double click Si1
On the card Edit, change Uiso to 0.005
OK; OK; Yes
[On screen: Main window of Jana2006]
“Parameters → Electron diffraction”
Default excitation error is too small for the current data – it must be increased:
Change Maximal excitation error to 0.03
Check the box next to “Use dynamic approach”
OK to create refinement reflection file
Follow the wizard and accept default choices
[On screen: Electron diffraction dialog]
Now we have to find initial values for the sample thickness, scale factors and orientations of
all three zones. As a first quick check, if the data quality is acceptable, we will calculate only
the thickness curve, using fixed orientation. This calculation is much faster that the
optimization of orientation, and it provides a first estimate of the data quality.
Click “Optimize thickness and scale”
Jana2006 generates reflections and calls Dyngo.exe for each zone separately. Dyngo.exe
calculates the reflection intensities using dynamical diffraction theory. Wait until Dyngo.exe
terminates.
[On screen: Electron diffraction dialog]
Click “Show thickness plots”
Go through the thickness plots. You should see three curves with more-or-less well defined
minima between 500 and 1000 Å. The best R-values are around 15%. This is a good start,
but it is still not satisfactory. It is necessary to optimize the orientations:
Quit the display of the thickness plot
Click “optimize thickness, scale and orientation”
This option optimizes all three parameters simultaneously, using a simplex algorithm. This
calculation requires repeated evaluations of reflection intensities and is thus much longer
that just the thickness optimization. Wait unil all three runs of Dyngo finish.
[On screen: Electron diffraction dialog]
Check the refined tilt angles. EDtheta is the azimuth of the tilt, EDphi is the amplitude. You
will see that the refined misorientations range between 0.13 and 0.39 degrees – very
reasonable values.
Click “Optimize thickness and scale”
There is a clear and reasonable minimum at t=350 Å, but the lowest R-value is found for the
thickness of 2000 Å. It is necessary to manually reset the thickness to 350 Å before the
refinement.
Quit; Locate zone 3 and set the thickness according to the following screenshot:
The data are now checked and ready for least-squares refinement.
OK; Yes to exit the Electron diffraction dialog
5. Structure refinement
[On screen: Main window of Jana2006]
The dynamical refinement tends to be unstable with default damping factor, it is thus better
to decrease it:
Right click “Refine”
Set damping factor to 0.2
OK; Yes+Start
The dynamical refinement starts. You may follow the progress in the window of Jana2006.
You may also interrupt the refinement or change the parameters the same way as for the
standard refinement.
The refinement converges after 21 cycles to R(obs)=6.8%. (Compare with 23% obtained with
kinematical refinement).
By repeating the optimization of orientation and rerefining the structure, the R-value can be
decreased even more to 5.0%. However, Uiso(Si1) remains negative.
Electron diffraction data was measured on a transmission electron microscope Philips CM120
with a precession divice DigiStar. Accelerating voltage of 120 kV.
Input data:
directory dp-150 with the set of experimental data (reflection images)
directory img with the set of experimental images from the crystal tilt
Ni2Si.pts: input file for PETS
Ni2Si_start.cif_pets: result of the data processing. This file is provided so that this
example can be solved without the data processing part, and to warrant
reproducibility of the solution and refinement part. File Ni2Si.cif_pets will be the
output of the data reduction.
Ni2Si_start_dyn.cif_pets: result of the data processing for dynamical refinement. This
file is provided so that the dynamical approach can be used without the processing
part, and to warrant the reproducibility of the refinement.
type a[Enter]
PETS starts the analysis of the peak positions obtained in the peak hunting procedure. In
the first step, distance distribution between peaks in the image plane is analyzed. A sorted
plot of inter-peak distances is displayed (white curve on screen, black on the image below)
with its derivative (grey curve):
For good quality data set, the white curve contains distinct steps, and, consequently, the
grey curve has distinct peaks.
Press Enter
In the next step an autoconvolution of the diffraction pattern (difference space) is
analyzed, and the groups of peaks in the autoconvolution (clusters) are replaced by the
cluster centers. Again, a distance distribution is displayed, and clear-cut steps on the curve
indicate a well-defined lattice.
Press Enter
Now the first step is finished and PETS prepared the files needed for unit cell indexing.
Keep the window of PETS open, you will need it later!
IMPORTANT: before you continue, make sure you have the same orientation of the lattice
on your computer as displayed on the picture above (check the direction of the axes on the
right). If not, you may obtain slightly different results in the rest of the tutorial. While not
necessarily wrong, it makes following this tutorial more difficult.
Click the button in the lower-left corner and choose the radio button for a vector
setting (right to the button ).
Click, hold the mouse button, and drag the marker to obtain the result on this image:
7. Import Wizard
Select “Known diffractometer formats”; NEXT
Select “Pets electron diffractometer”
The file name automatically changes to Ni2Si.cif_pets, which is a file produced by PETS
containing all important information, including the list of intensities.
NEXT
Leave all settings unchanged (note the wavelength of 120kV electrons!); NEXT
Leave all settings unchanged; NEXT
The program reads 614 reflections from the file
OK
For absorption correction select “None or done before importing”; NEXT
FINISH; OK; YES to accept the data set
You just read in the cell parameters and the intensity list. If you do not have the data in
the cif_pets format, you can read them in from a general hkl file using the option
“reflection file corrected for LP and absorption” at the beginning of the import wizard. You
then have to input the radiation type, wavelength and cell parameters by hand.
8. Symmetry Wizard
518 Jana2006 Cookbook, version October 2017
NEXT to close the information window and start the symmetry wizard.
Symmetry Wizard can be started separately by “File → Reflection file → Make space group
Test”. The default settings are prepared for x-ray diffraction. For electron diffraction we
have to increase the tolerances. Change the tolerances according to this screenshot:
“Select zones for editing Select all OK” and type 150 for “EDThick”. OK, Yes.
Right-click on icon “Refine”, change the number of cycles to 100 and the damping factor
to 0.3 in order to avoid instability of the refinement.
The refinement runs 23 cycles, each cycle takes less than 2 minutes (depending on your
computer). You can see the details of the refinement in the file Ni2Si.ref (Edit/View View
of Refine). The refinement converges with R=7.79%. However, the ratio reflections per
parameter refined is around 5 (454/84) and it is advisable to increase it. It can be done by
changing RSg (selecting more reflections for the refinement).
Another important thing to check is the number of observed reflections per zone. Also in
the Ni2Si.ref file, you can see “R-factors for individual zones”, where it is shown the list of
zones with the respective R-factors and number of reflections.
Now the refinement converges with R= 6.52% after 15 cycles. If you check the ellipsoid
model, you will see that Si is elongated along a*, which was the noisiest direction when
defining the unit cell.
Warning: this is a complex and time-consuming example. The computations in the last stages
may take up to several hours.
Unit cell is orthorhombic, primitive. The following screenshot is the view of the reciprocal
lattice along the axis c (a* horizontal, b* vertical):
Quit the indexing plugin to create file opx.smr with the orienting matrix, and continue data
processing in PETS by typing u[Enter] to read in the orienting matrix, and i[Enter] to
integrate the intensities.
PETS produced file opx.cif_pets, which contains intensities suitable for structure solution and
kinematical refinement. We have to produce another intensity set for dynamical refinement:
In PETS type d[Enter]
PETS will make another integration, this time producing file opx_dyn.cif_pets.
Type h[Enter] and y[Enter] ]to produce a standard set of sections through reciprocal space.
The sections will be useful for confirming the space group.
With this step the data processing in PETS is finished.
Type q[Enter] to close PETS. You may also close ImageJ.
4. Symmetry Wizard
NEXT to close the information window and start the wizard
Symmetry Wizard can be started separately by “File → Reflection file → Make space group
test.
The default settings are prepared from x-ray diffraction. For electron diffraction the
tolerances should be increased so that Jana properly detect that the cell is orthorhombic.
Change the tolerances according to this screenshot:
9. Structure refinement
[On the screen: basic window of Jana]
Double click the icon Refine.
Refinement converges smoothly to an R value about 25%.
The octahedral positions are not exclusively occupied by iron or magnesium, but by both
elements. We have to refine the mixed occupancy in these sites.
[On the screen: basic window of Jana]
Double click icon Edit atoms
Right-click atom Fe1 → Split atomic position
Set up the dialog according to the following screen shot:
The procedure generated a pair of atoms with restricted ADP parameters and coordinates
and with a restricted sum of occupancies. The only small problem (although rather esthetic
than real) is that in one pair the Fe atom is the “Master” and Mg the “slave”, while for the
other pair it is vice versa. We will modify the restrictions to be both the same.
[On the screen: basic window of Jana]
Right-click Refine → Various → Restrictions
Load the restriction on atoms Mg2 and Fe2
Exchange the order of atoms in the restriction (type Fe 2 Mg2 in the line
“Atoms/molecules), so that Fe is the “master” atom. Rewrite the restriction. When finished,
the dialog should look like this:
At this moment there is very little we can do with the kinematical refinement. We have a
satisfactory refinement, the structure makes sense, the occupancies have been refined, but
their accuracy cannot be estimated from the kinematical refinement. Furthermore, some
ADPs are negative. We will now proceed to the structure refinement using dynamical
diffraction theory.
File → Structure → Save as
Save the structure as opx_dyn
Answer Yes to the question if you want to continue with the new structure
File → Reflection file → Import/modify the reflection file
Reimport; Browse for the file opx_dyn.cif_pets; NEXT; NEXT; NEXT; OK; NEXT; Finish; OK;
Yes
You have just replaced the previous reflection data with new data set suitable for dynamical
refinement.
Cancel to avoid the production of refinement reflection file
If we continued, we would have produced a symmetry-averaged file for kinematical
refinement. For dynamical refinement the intensities on each frame must be kept separate.
We will now enable the dynamical refinement option:
Go to “Parameters → Electron diffraction”
Check the checkbox “Use dynamical approach”
Proceed with creation of the reflection file, leave all settings default.
At this point we are ready to start dynamical calculations. Because the calculations are very
time consuming, we will not, in this example, use all 90 zones, but only nine selected zones.
Click “Select zones for refinement”; In the line below the table, set the following:
The data collection was performed on one crystal and there is also no reason to expect that
the thickness varies for each zone. We will therefore restrict the thickness to the sam value
for all zones:
Click “Define zones of equal thicknesses; Select all; OK
Dyngo – the external executable that performs the heavy calculations – can run in a parallel
mode. In the field “Number of threads” you can decide, how many parallel threads Dyngo
may use. For the fastest performance set this number to the number of cores of your
processor. Important! Most Intel processors indicate twice as many threads available than
As a first step we have to check, if our model fits the data sufficiently well to start the
refinement. For that, we will run thickness optimization. This procedure is relatively fast, and
for each selected zone it calculates an R-value as a function of thickness. For good data set
and reasonable starting model, this curve should have a minimum at a meaningful
thickness, and all curves should indicate approximately the same thickness.
Check the checkbox next to “Thickness” and click
Jana calls Dyngo for each zone. When the calculation is finished, each selected zone has next
to the zone number in the dialog the best R-value, e.g.:
All zones’ R-values are lower than 20%, some of them even lower than 10%. This is a pretty
good start.
Click “Show thickness plots”
Now you can view the curves R-value vs. thickness for all zones. In most cases the curve has
a well-defined minimum, only in zones 30 and 70 the minimum is at the lowest thickness.
The best thickness varies in a relatively broad interval, but the minima are broad, and our
structure model is far from perfect at this stage, so this is not surprising.
Quit the viewer of the thickness curves
We want the thickess to be the same for all zones in the refinement. Now the thicknesses
are all different and we should reset them to a common most appropriate value. The best
practice seems to be a visual estimation of the “typical” thickness in the thickness plots. In
this case 500 angstroms seems to be a good guess. Do not worry about the inaccuracy of
such guess – the thickness is a very robust parameter and it will refine to the correct value
even from a quite distant start.
Click “Select zones for editing”; “Select all”; OK
Double click the yellow field EDThick and type “500” in the field
This procedure reset the thickness of all zones to 500 angstroms.
We have to recalculate the scale factors for each zone, because the scale changes with
changing thickness:
Click the button without any other checkbox checked (neither thickness nor
Orientation)
When the calculation finishes,click OK to close the electron diffraction dialog.
In the next step we will optimize the orientation of each zone. Due to the inaccuracy of
goniometer positioning and orientation matrix, the positional angles of each diffraction
pattern are not exactly correct, and a correction should be found.
Go to “Parameters → Electron diffraction”
Check the checkbox next to “Orientation”
532 Jana2006 Cookbook, version October 2017
Click the button
An optimization procedure starts that finds the best zone-axis orientation for each pattern.
Wait until the calculation finishes. Now you can now run a few more cycles of structure
refinement.
Click OK to close the form “Electron diffraction”
Double click the icon Refine
The refinement converges after 10 cycles to an R-value of 5.10%. Inspection of the structure
model shows that the occupancies of Fe1 and Fe2 refined to 0.33(3) and 0.06(2),
respectively. This is still relatively far from the values from x-ray refinement, which are
0.432(3) and 0.157(3). However, there are only 188 observed reflections in this limited
dataset for the refinement of 52 parameters. The data to parameter ratio is thus poor. If all
90 zones are included in the refinement and the whole procedure is repeated with them, the
final result is:
Robs=5.89%
occ(Fe1)=0.389(7)
occ(Fe2)=0.148(7)
4. Refinement
Run refinement
It converges with R=3.2%
Open refinement listing, GoTo Statistics
Higher R value for the strongest reflections indicates that extinction correction is needed
Close the refinement listing
Start “Parameters → Extinction”
Activate Isotropic extinction with Gaussian model
Save changes and leave to the basic window
Repeat refinement
In converges with R=2.75%. The reason extinction has not been imported from the ins file is
that SHELX and Jana use different parameters for extinction correction
4. Refinement
Run refinement
It converges with R=4.1% for 60 parameters. There were very large changes of the twin
fractions during the refinement because probably the order of twin fractions in SHELX and
Jana differs.
Open refinement listing, GoTo Statistics
There is no indication of extinction and at the bottom of the listing no warnings.
As Jana2006 allows refinement against more data sets, PRF file may consist of several data blocks:
Block1 begin
.........................
.........................
Block1 end
Block2 begin
.........................
.........................
Block2 end
In the case that there is just one data block, the labels are omitted.
Each block consists of three parts: header, Bragg positions and profile points:
2 0 0 2 3 3
1 -1 2 3. 1 10.9743 0.0237 0.0380 0.547364E+01
1 0 -2 3. 1 10.9743 0.0237 0.0380 0.547364E+01
1 -1 -4 3. 1 13.7339 0.0236 0.0414 0.698639E+01
1 0 4 3. 1 13.7339 0.0236 0.0414 0.698639E+01
0 0 -6 1. 1 14.2831 0.0236 0.0421 0.130568E+01
0 0 6 1. 1 14.2831 0.0236 0.0421 0.130568E+01
2 -1 0 6. 1 17.1683 0.0235 0.0456 0.668545E+01
2 -1 -3 6. 1 18.6097 0.0235 0.0474 0.768642E-01
2 -2 4 3. 1 22.0509 0.0234 0.0517 0.238082E+01
2 0 -4 3. 1 22.0509 0.0234 0.0517 0.238082E+01
2 -1 6 6. 1 22.4023 0.0234 0.0522 0.573109E+00
2 -1 -6 6. 1 22.4023 0.0234 0.0522 0.573109E+00
1 0 1 8. 2 25.3002 0.0232 0.0598 0.194769E+02
....................................
....................................
999
6.000000 0.209300E+04 0.249015E+04 0.457493E+02 6.023784 0 0.000000E+00 0.000000E+00 0.249015E+04
6.016766 0.209500E+04 0.247949E+04 0.457712E+02 6.040549 0 0.000000E+00 0.000000E+00 0.247949E+04
6.033532 0.217100E+04 0.246891E+04 0.465940E+02 6.057315 0 0.000000E+00 0.000000E+00 0.246891E+04
6.050298 0.206100E+04 0.245839E+04 0.453982E+02 6.074081 0 0.000000E+00 0.000000E+00 0.245839E+04
6.067064 0.204200E+04 0.244793E+04 0.451885E+02 6.090847 0 0.000000E+00 0.000000E+00 0.244793E+04
6.083830 0.217300E+04 0.243755E+04 0.466154E+02 6.107612 0 0.000000E+00 0.000000E+00 0.243755E+04
6.100596 0.210500E+04 0.242723E+04 0.458803E+02 6.124378 0 0.000000E+00 0.000000E+00 0.242723E+04
...............
...............
10.593884 0.102600E+04 0.986969E+03 0.320312E+02 10.617599 1 0.000000E+00 0.000000E+00 0.986969E+03
10.610650 0.108400E+04 0.101817E+04 0.329242E+02 10.634364 1 0.340796E+02 0.000000E+00 0.984355E+03
10.627417 0.991000E+03 0.101999E+04 0.314802E+02 10.651131 1 0.385393E+02 0.000000E+00 0.981752E+03
10.644182 0.105600E+04 0.102260E+04 0.324962E+02 10.667895 1 0.437751E+02 0.000000E+00 0.979160E+03
10.660948 0.104400E+04 0.102623E+04 0.323110E+02 10.684661 1 0.500343E+02 0.000000E+00 0.976579E+03
...............
...............
24.710857 0.361000E+03 0.316545E+03 0.190000E+02 24.734123 1 0.000000E+00 0.000000E+00 0.316545E+03
24.727623 0.347000E+03 0.316369E+03 0.186279E+02 24.750887 1 0.000000E+00 0.000000E+00 0.316369E+03
24.744389 0.341000E+03 0.353730E+03 0.184662E+02 24.767653 1 0.000000E+00 0.378276E+02 0.316194E+03
24.761154 0.353000E+03 0.356299E+03 0.187883E+02 24.784416 1 0.000000E+00 0.405929E+02 0.316019E+03
24.777920 0.385000E+03 0.359047E+03 0.196214E+02 24.801182 1 0.000000E+00 0.435381E+02 0.315845E+03
24.878517 0.375000E+03 0.383491E+03 0.193649E+02 24.901773 1 0.000000E+00 0.692188E+02 0.314806E+03
...............
...............
999.
Bragg positions
h,k,l,... n diffraction indices, where n = max(NDim(1), …, NDim(NPhases)).
For phases having the dimension smaller than n the indices are completed
by zeros.
Multiplicity multiplicity of the reflection following from the symmetry
KPhase defines the phase the reflection belongs to
The following four items are repeated twice if KADoublet=1. The first block corresponds to Kα1
TwoTheta/Time angular position/time (or for TOF d value/energy) of the peak as calculated
/d/energy from cell parameters and modulation vectors
Shift the angular shift of the peak position due to a zero shift, sycos and sysin
parameters (it need not be constant for all reflections). The actual of
reflections with respected to the profile should be calculated as:
[TwoTheta/Time/d/energy] -Shift
FWHM in the same units are previous values
Icalc(h,k,l) intensity as calculated from Rietveld (corrected for multiplicity and
preference orientation) or as extracted from the profile by le Bail fit.
The list of Bragg peaks ends with the integer 999 at the position of the h index.
Profile points
TwoTheta/Time/d/energy diffraction angle/time for TOF/d value/energy of the profile
point
Iobs observed intensity
Icalc calculated intensity
sig(Iobs) standard uncertainty of Iobs
TwoTheta-corr/Time-corr/ diffraction angle/time for TOF/d value/energy of the profile
d-corr/energy-corr point corrected for zero shift, sysin and sycos
Icalc calculated intensity
Skip flag 0 if the profile point was used in the refinement
1 if the profile point belongs to one of excluded region
Icalc(1) calculated intensity of the first phase
....
Icalc(NPhases) calculated intensity of the first phase
Ib background intensity
The list of peak points ends with the real number “999.”.
Refinement on F and on F2
The old programs used usually the refinement on F. This means that the minimized function has
the following form:
S F = ∑ wi Fi ,obs − Fi ,calc ( )2
i (1)
wi = 1 σ (Fi ,obs )
2
where . The refinement based on F2 uses an analogical formula:
i (2)
wi = 1 σ (F ). According to the error propagation theorem standard uncertainties
2 2
but now i ,obs
σ (F ) = 2 F σ ( F ) 2
are related by the equation which leads to:
i ,obs i ,obs i ,obs
SF2 = ∑
1
(F −F ) =∑
2 (F i ,obs + Fi ,calc ) (F
2
− Fi ,calc ) 2
(F ) (F )
2 2
σ σ
2 2 i ,obs i ,calc 2 2 i ,obs
4F 4F
(3)
i i ,obs i ,obs i i ,obs i ,obs
Fi ,obs ≈ Fi ,calc
In final stages of the refinement and therefore as an approximation we have:
2
SF2 ≈ ∑
4 Fi ,obs
(F − Fi ,calc ) = ∑ σ ( F1 ) ( F
2
− Fi ,calc )
2
= SF
i
2
4 Fi ,obs σ 2 Fi ,obs ( ) i ,obs
i
2
i ,obs
i ,obs
This means that the results for most cases should not give so different results.
However more precise analysis gives:
SF2 = ∑
(2 F i ,obs − Fi ,obs + Fi ,calc ) (F − F ) =
2
2
) ( F − F ) + ∑ 4F )( F )
1 1
SF − ∑
3 4
− Fi ,calc
i Fi ,obs σ 2
(F i ,obs
i ,obs i ,calc
i
2
i ,obs
2
(
σ Fi ,obs
i ,obs
S
This means analogically that difference between the minimized functions F and S F is more
2
∑ 4F
1
(F − Fi ,2calc )
2
∑ σ ( F ) (F
1
− Fi ,calc )
2
∑ σ ( F ) (F
1
− Fi ,calc )
2
(F )
2
σ
2 2 i , obs 2 i , obs 2 i , obs
i i , obs i , obs i i , obs i i , obs
RwF 2 = = =
1 1 1
∑ 4F
i
2
σ 2
(F ) Fi ,2obs ∑ 4F
i
2
σ 2
(F ) Fi ,4obs ∑ 4σ ( F ) F
i
2
2
i , obs
i , obs i , obs i , obs i , obs i , obs
∑ σ ( F ) (F − Fi ,calc )
1 2
2 i ,obs
i
RwF 2 ≈ 2 = 2 RwF
i ,obs
1
∑i σ 2 F Fi,2obs
( )
i ,obs
which means that the weighted R values based on F2 is about two times larger than that based
on F.
As obvious from the previous the weighted R values are closely related to the minimized
function and therefore they should reach their minimal values mutually.
On the other hand the non-weighted R values are defined:
∑F −F i ,obs i ,clac ∑ F −F 2
i ,obs
2
i ,calc
RF = i
RF 2 = i
∑F i
i ,obs ∑F i
2
i ,obs
and they need not to reach their minimal values at the same time as the minimized function.
RwF 2 RF 2 ≈ 2 RF
Similarly as for and RwF we can get
RwF 2
In the output of Jana2006 you can find either RwF or depending on the type of refinement.
As a non-weighted R value we use always RF .
The volume fraction vi is defined as a ratio of the individual domain Vi volumes to the sample
volume V :
V
vi = i
V
In the case of the twinned crystal with complete overlaps of diffraction spots (merohedric twins
or pseudo-merohedric twins with lattice fulfilling higher symmetry within experimental errors)
this volume fractions can directly be used to express a “combined structure factor” F
2
(h )
N
F 2 (hTi )
F 2
(h ) = ∑ vi 2
i =1 Vcell
where F (h ) is the structure factor, Vcell is the cell volume and Ti are the twinning matrices
and N is the number of diffractionally independent twin domains.
The formula can be generalized in the cace for general pseudo-merohedric twin:
(h ) = ∑ vi pi F (hTi)
N 2
F 2
2
i =1 Vcell
where the factor pi is equal 1 if the diffraction spot of i-th domain, having indices hTi , is
contributing to the spot localized at the point h and otherwise it is equal zero. Note that in this
case the diffraction vectors hTi need not have integer components.
Another generalization can be made for multiphase samples. Then the previous formula is
changed to:
N F j2 (hTi )
F 2
(h ) = ∑ vi 2
i =1 Vcell ,j
where the index j defines the j-th phase domain.
Similar formulas are used in the Rietveld refinement. Then a contribution to the profile point is:
N F j2 (h i )
ycalc (2ϑ ) = Lp (ϑ )R(ϑ )∑ vi Ai (ϑ ).∑ Pi (h i )mi (h i ) 2
Ω i (h i ,ϑ )
i =1 hi Vcell ,i
Lp (ϑ ) Lorentz-polarization factor
R (ϑ ) roughness correction
Ai (ϑ ) absorption factor
Pi (h i ) preference orientation function
mi (h i ) multiplicity factor
Ωi (h i ,ϑ ) profile function
The volume fractions can be transformed to the relative phase amounts in mass by the formula:
vi ρ i
ri= N
∑v ρ
i =1
i i
The both modes can be combined with the so call “fixed” commands which can fix different
groups of parameters. The following fixed command fixes for example coordinates of all atoms:
which means that cos modulation wave of all first harmonics for atoms having name started with
the character “C” with be set to zero and kept unrefined. For this is important that every
parameter is Jana2006 has its own name which is used in all output and input procedures.
Abstract
A full solution and refinement of modulated structures, having strong occupational modulation described
by crenel functions, is possible only if modulations of other parameters are described with respect to a set
of functions orthogonalized in the crenel interval. It is shown how such modulation can be recognized
from so called “de Wolff’s sections”. Moreover, two new sets of functions (Legendre polynomials and x-
harmonics) were introduced and tested on simulated and real structures with a new version of program
Jana2006.
1. Introduction
The superspace theory as developed by de Wollf et al., (1981) proved to be the best way how to generalize
methods for solution and refinement of regular structures to modulated and composite ones. Any structural
parameter (occupancy, atomic coordinates and ADP) can be expressed as a periodic modulation function
of the internal coordinate x4 :
p ( x4 ) = p ( x4 + 1) (1)
The actual value of the structural parameter p of the atoms located at the position r0 + n , where n defines
the position of the actual cell, is calculated from the equation (1) with x4 = q ⋅ (r0 + n ) + t . Here q stands
for the modulation vector and t defines a phase factor defining a "first" cell.
The fact that the modulation function is periodic allows its expansion into a Fourier series:
But as has been already shown (Petříček et al., 1990 and Petříček et al., 1995) there are cases where the
number of harmonics, necessary to describe a modulation function, is very large. This is true especially of
modulation functions exhibiting some discontinuity - such as positional saw-tooth function and
occupational crenel function.
It turns out that with growing number of successfully solved structures the number of cases, having some
type of discontinuity in modulation functions, increases. These functions are also often used to describe by
a unified way families of related structure ( e.g. Evain et al., 1998; Peréz-Mato et al., 1999). Therefore in
the first part of the paper it is shown how such discontinuities modulations can be correctly recognized by
analyzing of de Wolff’s sections through periodic (3+1) dimensional Fourier maps.
The basic problem in application of the crenel and saw-tooth functions is that harmonic functions (2)
defined in the interval (0,1) are not anymore orthogonal when used in the shorter interval reduced by the
crenel function. We shall call it in following text the crenel interval. This leads to strong correlations
between coefficients of positional or ADP modulations. On the other hand harmonic functions, defined in
the full interval, allow an expansion of the structure factors into a series of products of Bessel functions
(Peréz-Mato, et al.,1986 and Petříček & Coppens, 1988) which was the only method used in the old version
of the Jana program. For this reason a combined method of selection of harmonics, defined in the full
interval, and their re-orthogonalization into the crenel interval has been developed (Petříček et al., 1995).
This method has been successfully applied by different authors to various modulated structures. The main
disadvantage of this method is that depending on the crenel interval the set of used functions (linear
combinations of basic harmonics) are generally different and therefore they must be always listed together
with set of refined coefficients to allow a full reproduction of refined modulation curves.
As the Bessel expansion method was much less efficient for more complicated modulated structures than
the integral method introduced in the program REMOS by Yamamoto, (1982) the newer versions of Jana
program are offering both methods for calculation of structure factors. This opened possibility to reconsider
way how to define modulation function within a crenel interval. In the section 3 two new orthogonal sets
of functions, defined directly in the crenel interval, are introduced and tested with simulated structures.
As it has been presented in the old paper by de Wolff (1974), the generalized (3+1) dimensional maps based
on the already know model can be used to prove that the solution is correct or to find more details about
modulation functions. Similarly like for regular structure we can use Fourier, difference Fourier or Patterson
maps. Probably the most commonly used are 2d sections (x1-x4, x2-x4 or x3-x4) of (3+1) dimensional
However there are some effects which can flutter of such maps in so strong way that a wrong conclusion
can be drawn. In many cases a positional modulation affects more than one individual coordinate and the
modulation looks like for example helically or elliptically modulated distortions. Then the simple de Wolff's
sections can look like we can see in Fig. 1a. The fact that the atom is modulated in two directions leads to
artificial effects which could make a false impression of discontinuity in modulation. This can be solved
by calculating sections in other directions within some vicinity of the selected atom and projecting it (see
Fig 1b).
Another problem is connected with resolution caused by a termination effect in the Fourier summation. For
modulated structures the crucial point is if all detectable satellites were measured and used in the Fourier
calculation. This can be demonstrated by a simulated structure in which the atom has a strong harmonic
positional modulation along c axis with amplitude about 1 Å. Then the de Wolff's section calculated with
only first order satellites gives an impression that the atom is predominantly located in two distinct positions
as shown in Fig 2a. However the Fourier maps made with satellites up to fifth order (Fig 2b) shows clearly
that the modulation has purely harmonic character.
A possible smearing of a Fourier map can be estimated from so called shape function which is a Fourier
summation: S (r ) = ∑ A(H ) exp(− 2π r ⋅ H ) with A(H ) = 1 for all measured diffractions. For all others
H
diffraction points A(H ) = 0 . The calculated map is a convolution of an ideal map with the shape function.
The more indices are measured the more δ-like character the function S (r ) has and the more truthful the
map is. Fig 2c and 2d show the shape function for two simulated cases mentioned above. A strong
prolongation of the shape function along x4 axis, for the case in which only main reflections and first order
satellites are used, leads to severe smearing as visible in resulting map (Fig 2a). This is considerably
changing the character of the modulation curve in points where derivation of the ideal modulation curve
with respect to x4 is large.
On the other hand the refinement process is much less sensitive to the maximal order of used satellites. This
can be demonstrated by a simulated refinement against a data set made only from main reflections and first
order satellites. Then the refinement of a crenel-split model, as can be concluded from a Fourier maps based
on the same set reflections (see Fig 2b), gives a very bad fit. The reason for this disagreement can be
understood from differences between contributions to structure factors for a crenel-split model and the
model based on a strong positional modulation. In the simulation we used the modulation vector of form
sin πm∆
F (h, k , l , m ) ~ exp(2πiH ⋅ r0 ) sin[2π ( x3b − x3a )H 3 ] (3a)
πm
Where r0 is an average position of the split atoms, x3a and x3b actual z-coordinates of both atomic positions
and ∆ is crenel length, here equal to 1/2.
On the other hand the positional modulation leads to the following structure factor:
The fig. 3 shows that these two expressions, as a function of (x3b − x3a )H 3 and u3 H 3 , have a quite
different form even for first order satellites.
The refinement of the crenel-split model is to be combined with a positional modulation described by a set
orthogonalized modulation function as will be described later. Then one modulation wave is sufficient to
get a perfect fit and the modulation curve, as demostrated on the fig 4, which is almost identical with the
original harmonic wave.
From what was demonstrated in this chapter we can see that the de Wollf's sections can help to recognize
possible discontinuities in the modulation curves and get more realistic description of the modulated
structures. Such an analysis of (3+1) dimensional density maps can be also used to extract a correct solution
from maps supplied by Superflip (Palatinus & Chapuis, 2007).
As it was mentioned in the introduction the crenel occupancy modulation leads to severe numerical
problems when combined with modulations of other parameters (position and ADP) if the standard
description with harmonic functions (2), defined in the interval 0,1 , are used. This problem has been
solved by Petříček et al. (1995) by an orthogonalization method combined with a selection of harmonics.
Here we would like to repeat the main ideas from this paper and give two more alternative methods for
description of modulations within the crenel interval. Finally all three methods are compared on a simulated
data set.
{g i (x ), i = 1,2n + 1} = {1, sin (2πx ), cos(2πx ),sin (2πnx ), cos(2πnx )} for atomic positions and harmonic
ADP's, defined in the interval 0,1 , leads to analytical expressions for a structure factor based on sum of
products of Bessel or modified Bessel functions which is true even for cases when such modulations are
x40 − ∆ 2, x40 + ∆ 2 .
The Gram-Schmidt orthogonalizion procedure is used to construct a new set of functions {γ i (x )}, already
orthogonalized in the crenel interval, as a linear combination of the original harmonic functions to
{g i (x ), i = 1,2n + 1}:
2 n +1
γ j (x ) = ∑ T ji g i (x ) (4)
i =1
The Gram-Schmidt orthogonalizion procedure leads to the lower triangular matrix T which is generally
different for different crenel intervals. This means for each particular used crenel function the set of
coefficients forming the matrix T must be presented to allow full reconstruction of the modulated
functions.
Another point is that for numerical reasons the othogonalization method is to be combined with a selection
procedure in which each new function say g j (x ) from the sequence {g i (x ), i = 1,2n + 1} is first tested if it
makes an angle α j to already chosen subset of functions larger than a chosen limit. For practical application
As already mentioned in the original paper this leads to some incompleteness of the selected set of functions
but experience with application of this method to various compounds approved that our default value
λ = 0.95 is reasonable.
As already mentioned in the introduction, the analytical method for calculation of structure factors based
on Bessel functions was becoming less and less effective for complicated modulated structures and
integration method (Yamamoto, 1982) is now used almost exclusively. This called for reconsideration of
used set of orthogonalized sets defined directly in the crenel interval. Such function has two main
advantages to the original orthogonalized set. These sets are clearly and uniquely defined for any crenel
interval and therefore they have always the same mathematical form. Moreover, they form mathematically
complete sets of functions.
functions are not mutually orthogonalized again the Gram-Schmidt orthogonalizion procedure has to be
used. But this leads to a standard set of function depending on definition of the scalar product:
(g , g ) = ∫ w(x )g (x )g (x )dx
1
i j i j
−1
For different choice of the weight function w( x ) we can get e.g. Legendre, Chebyshev, Laguerre or Hermite
polynomials (Abramowitz & Stegun, 1972). The most natural choice for incommensurately structures is
the unit weight as all points in the interval where atom is present have the same probability. This condition
leads to the Legendre polynomials. For comparison of modulation amplitudes it important to take into
account that Legendre polynomials are not normalized but that the following equation holds:
( fn , fm ) = 2
δ nm
2+n
(
x = 2 x4 − x40 ∆ )
Another candidate for orthogonalized set of functions in the crenel interval is the set of harmonic functions
{1, sin (πx ), cos(πx ), sin (2πx ), cos(2πx ),sin (nπx ), cos(nπx )} . But as all these functions give the same
value at points -1 and +1 they cannot fully describe all possible modulations within the crenel interval as
generally the modulations need not be identical in the first and the last point of the crenel interval. For this
reason the set of functions is to be completed with a function which would correct for this. The simplest
way is to add to this sequence a linear term:
{1, x, cos(πx ), sin (πx ), cos(2πx ), sin (2πx ),, cos(nπx ), sin (nπx )}
Then in the sequence of functions all functions with odd order are even and vice versa. While the original
set of harmonic functions in the crenel interval is mutually orthonormal, the new set of functions is not as
the following formula holds:
1
(x, sin nπx ) = ∫ x sin (nπ )dx = (− 1)n+1 2
(6)
−1
nπ
On the other hand all other mixed scalar products are equal zero. But this sill means that the Gram-Schmidt
orthogonalizion procedure is to be used to get fully orthonormalized set. In the Fig.5 the orthonormalized
functions of order from 2 till 6 are visualized. The set of coefficients for first 20 functions and analytical
expressions as derived by the program Maple 12 are supplied as supplementary materials. To our
The three sets functions - harmonics in 0,1 orthogonalized to the crenel interval (in the following shortly
"Ortho-harmonics"), Legendre polynomials and x-harmonics can give generally a different fit when used
to describe modulations in the crystal. The advantage of the Legendre polynomials and x-harmonics is that
they make a complete set of functions and therefore they can describe any modulation function if a sufficient
number of terms is used. In order to demonstrate different situations two simulated structures with structural
parameters as summarised in the Table 1 have been used. In the simulation procedure noise from a random-
number generator was applied to ensure conditions similar to real experiment.
Both cases exhibit very strong positional modulation within the crenel interval and they are very far from
what is usual for real structures. However such tests can show differences between using of different set of
functions. The results of three refinements, in which the same number of positional parameters was used as
for simulations, are summarized in the Table 2. R values for the case x40 = 0.2 is almost the same for all
sets but for the second case the x-harmonics makes considerably better fits.
A comparison how refined modulation curves fit to those from the simulated model are shown in the Fig 6
and 7 for the case x40 = 0.2 and x40 = 0.6 , respectively. For refined R values and graphic curves it is
obvious that the x-harmonics gives the best fit. However after making more tests, also with some already
refined and published structures, it turned out that differences are not so pronounced as for the selected
simulated structures and that all three methods could give similar results.
4. Conclusions
The original method for refining of modulations of structural parameters combined with crenel occupancy
modulations used the “ortho-harmonics”, as was described by Petříček et al (1995), have two main
disadvantages to the newly introduced orthogonalized sets of functions. Their coefficients are to be always
updated when the crenel interval is changed. This is especially crucial during the refinement process. The
second disadvantage is that modulations functions can be reproduced only if together with refined
amplitudes a set of actual orthogonalization coefficients is supplied.
The newly introduced sets of functions, Legendre polynomials and x-harmonics, are defined in the crenel
interval in a unique mathematical form and therefore they make the refinement smoother. While the
Legendre polynomials are standard functions x-harmonics, to our knowledge, were not yet used and
tabulated. For this reason all details how to use them are given as supplementary materials. From made
tests we found that the refinement with x-harmonics gives usually better fit than compared to Legendre
polynomials and that the normalized modulated amplitudes are smaller.
Acknowledgements The work has been supported from Praemium Academiae of the Czech Academy
of Sciences and by the institutional research plan No. AVOZ10100521 of the Institute of Physics. One of
the authors (J.F.) acknowledges a financial support of the Grant Agency of the Academy of Sciences of the
Czech Republic under the grant reg. no. IAA100100918. One of the authors (V.E.) acknowledges a
financial support of Ministry of Education, Youth and Sports of the Czech Republic under grant reg. no.
MSM6046137302).
4.1. References
Formula: Ta
Basic unit cell: a = 3.018, b = 4.021, c = 3.815 Å, α = 80, ß = 75, γ = 80
Modulation vector: q = 0.3412a* + 0.4091b* + 0.2752c*
Superspace group: P 1 (αβγ )0
hklm range: all main reflections and satellites up to fifth order ( sin θ λ ≤ 1 Å-1)
n Ax ,s ,n Ay , s , n Az ,s ,n Ax ,c ,n Ay ,c ,n Az ,c ,n
0 0.2500 0.2500 0.2500
1 0.0500 -0.0500 0.0500 -0.0250 -0.0350 0.0350
2 0.0250 0.0100 -0.0150 0.0500 0.0500 -0.0250
3 0.0350 0.0150 0.0300 0.0250 0.0150 0.000
U 11 U 22 U 33 U 12 U 13 U 23
0.0150 0.0150 0.0200 0.0050 -0.0050 0.0050
Legendre
Ortho-harmonics x-harmonics
polynomial
Case RF RwF RF RwF RF RwF
x40 = 0.2 0.0703 0.0889 0.0752 0.0939 0.0675 0.0810
2
0
2
0
x
0
x
x
0
0
,
x
x
4
x
4
6
6
2
2
0
8
0
8
0
4
0
4
0
0
0
0
-
0
2
0
-
0
2
0
0
0
0
x
x
(a) (b)
(a) x1-x4 de Wolff's section through the atom having helical distortion in the (x1, x2) plane. The electron
density shows an artificial discontinuity.
(b) x1-x4 de Wolff's section through the atom having helical distortion in the (x1, x2) plane with projected
section in x2-x4. The artificial discontinuity disappears.
1
.
0
1
.
0
x
4
x
4
x
0
0
0
0
,
x
2
0
0
0
0
x
0
0
0
0
,
x
2
0
0
0
0
0
.
8
0
.
8
0
.
6
0
.
6
0
.
4
0
.
4
0
.
2
0
.
2
0
.
0
0
.
0
-
0
.
2
0
-
0
.
2
0
0
0
0
0
0
0
x
x
(a) (b)
x
4
x
4
0
.
4
0
.
4
x
0
0
0
0
,
x
2
0
0
0
0x
0
0
0
0
,
x
2
0
0
0
0
0
.
2
0
.
2
0
.
0
0
.
0
-
0
.
2
-
0
.
2
-
0
.
4
-
0
.
4
-
0
.
2
0
-
0
.
2
0
0
0
0
0
0
0
x
x
(c) (d)
(a) Fourier deWolff’s sections made from main and first order satellites
(b) Fourier deWolff’s sections made from main and satellites up to fifth satellite order
(c) Shape function made from main and first order satellites
(d) Shape function made from main and satellites up to fifth satellite order
Comparison of contribution to the structure factor for the split and harmonic model
1
.
0
x
4
x
0
0
0
0
,
x
2
0
0
0
0
0
.
8
0
.
6
0
.
4
0
.
2
0
.
0
-
0
.
2
0
0
0
0
0
0
Modulation curves refined with the split atom model described by crenel functions combined with
positional modulation.
Figure 7 Positional parameters x,y and z as refined for the case x40 = 0.6 . The red curve shows the
refined modulation curve while the black line shows the modulated wave as used in the simulation
procedure.
Dear Pavel,
I wish you all the best in 2014. I am very sorry for so long silence but
during the Christmas time I was just relaxing.
Would you be so kind to give me a piece of advice for SizeL, StrainL, SizeG,
StrainG values handling? The question, namely, is the units of these values:
is it angstroms for Size values and dimensionless epsilon for Starin values,
or any coefficients are used (as for LX,LY and GV,GU,GW,GP in simple PV
model)?
In the case that the fundamental approach is used the number of profile
parametrs is reduced to (LX, LY, GU and GP):
But instead of refining LX and GP we are using parametrs SizeL and SizeG.
They are related to LX and GP by equations:
Vaclav
Sincerely
Sebastian
E-mail: [email protected]
Dear Sebastien,
Yes, you are right the definition in Jana is different than originally
introduced by Stephens. The reason for this is that in Jana we are using
the
procedure for applying symmetry restrictions of symmetrical 4th order
tensor
as used also for anharmonic ADP's. So I am not planning to modify it.
Vaclav
Thank you for your letter and intent in Jana2006. The asymmetry
correction named in Jana as "correction by divergence" is based on the
theoretical paper by Finger, Cox and Jephcoat attached to this letter.
The two parameters labeled in Jana as "H/L" and "S/L" are based on the
symbols used in that paper - 2H is the detector opening and 2S is the
sample height. Both are divided by "L" which is the sample to detector
distance.
equation : H/L=S/L
(H-S)/L "in Jana "HmS/L and H+S/L in Jana "HpS/L". Than the condition
HmS/L=0 means H/L=S/L.
Vaclav