VESTA Manual
VESTA Manual
VESTA Manual
1
Koichi MOMMA
National Museum of Nature and Science,
4-1-1 Amakubo, Tsukuba, Ibaraki 305-0005, Japan
2
Fujio IZUMI
Graduate School of Engineering, Kyoto University
1
E-mail: vesta.dev@gmail.com
2
E-mail: fizumi3776@gmail.com
Contents
LICENSE AGREEMENT xi
1 INTRODUCTION AND BACKGROUND 1
3.1
Minimum requirements of hardware . . . . . . . . . . . . . . . . . . . . . . . . . 15
3.2
Windows . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
3.3
macOS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
3.4
Linux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
3.4.1 32-bit version . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
3.4.2 64-bit version . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
3.5 Trouble Shooting . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3.6 Loading Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
4 MAIN WINDOW 19
ii
4.2.6 Help menu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
4.3 Tools in the Horizontal Toolbar . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
4.3.1 Alignment . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
4.3.2 Rotation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
4.3.3 Translation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
4.3.4 Scaling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
4.4 Tools in the Vertical Toolbar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
4.5 Text Area . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
4.5.1 Output tab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
4.5.2 Summary tab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
4.5.3 Comment tab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
5 DISPLAY STYLES 27
iii
8 CREATING BONDS AND POLYHEDRA 57
8.1 Specifications of Searching for Bonds . . . . . . . . . . . . . . . . . . . . . . . . . 58
8.1.1 Search mode . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
8.1.2 Boundary mode . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
8.1.3 Options . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
8.1.4 A pair of atoms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
8.2 Operating Instructions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
8.2.1 Creating coordination polyhedra . . . . . . . . . . . . . . . . . . . . . . . 59
8.2.2 Searching for molecules and clusters . . . . . . . . . . . . . . . . . . . . . 60
8.2.3 Applications of the “Search molecules” mode to inorganic crystals . . . . . 60
8.2.4 Searching for hydrogen bonds . . . . . . . . . . . . . . . . . . . . . . . . . 61
8.2.5 Visualizing cage-like structures . . . . . . . . . . . . . . . . . . . . . . . . 61
9 ADDITIONAL OBJECTS 63
9.1 Vectors on Atoms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
9.1.1 Creation and editing of a vector . . . . . . . . . . . . . . . . . . . . . . . 64
9.1.2 Attachment of vectors to crystallographic sites . . . . . . . . . . . . . . . 66
9.1.3 Attachment of vectors to individual atoms . . . . . . . . . . . . . . . . . . 66
9.2 Lattice Planes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
9.2.1 Appearance of lattice planes . . . . . . . . . . . . . . . . . . . . . . . . . . 67
9.2.2 Calculate the best plane for selected atoms . . . . . . . . . . . . . . . . . 68
11 INTERACTIVE MANIPULATIONS 73
11.1 Rotate . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
11.1.1 Drag mode . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
11.1.2 Animation mode . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
11.2 Magnify . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
11.3 Translate . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
11.4 Select . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
11.4.1 Atom . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
11.4.2 Bond . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
11.4.3 Coordination polyhedron . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
11.5 Distance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
11.6 Bond angle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
11.7 Dihedral angle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
11.8 Interfacial angle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
iv
12 PROPERTIES OF OBJECTS 82
12.1 Properties Dialog Box . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
12.1.1 General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
12.1.2 Atoms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
12.1.3 Bonds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
12.1.4 Polyhedra . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
12.1.5 Isosurfaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
12.1.6 Sections . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
12.2 Objects Tab in the Side Panel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
12.2.1 List of phases and objects . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
12.2.2 Atoms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
12.2.3 Bonds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
12.2.4 Polyhedra . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
12.2.5 Slices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
12.2.6 Shapes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
13 OVERALL APPEARANCE 97
13.1 Background . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
13.2 Lighting . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
13.3 Projection Mode . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
13.4 Depth-Cueing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
14 UTILITIES 101
14.1 Equivalent Positions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
14.2 Geometrical Parameters . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
14.3 Standardization of Crystal Data . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
14.4 Niggli-Reduced Cell . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
14.5 Site Potentials and Madelung Energy . . . . . . . . . . . . . . . . . . . . . . . . . 106
14.6 Powder Diffraction Pattern . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
14.7 Structure Factors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
14.8 Fourier Synthesis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
14.9 Model Electron Densities . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
14.10Model Nuclear Densities . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
14.11Patterson Densities . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
14.122D Data Display . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
14.13Line Profile . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
14.14Peak Search . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
14.15Conversion of Electron Densities . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
v
15.5.1 (hkl) plane in the bounding box . . . . . . . . . . . . . . . . . . . . . . . 120
15.5.2 (hkl) plane defined by two vectors . . . . . . . . . . . . . . . . . . . . . . 120
15.5.3 Project along [hkl] axis . . . . . . . . . . . . . . . . . . . . . . . . . . . . 120
15.6 Controlling Properties of a 2D Image . . . . . . . . . . . . . . . . . . . . . . . . . 122
15.7 Exporting 2D data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124
16 PREFERENCES 125
vi
List of Figures
5.1 The Structural model frame box in the Style tab of the Side Panel . . . . . . . . . 27
5.2 Crystal structure of quartz represented as stick model with dot surfaces . . . . . 27
5.3 The Volumetric data frame box in the Style tab of the Side Panel . . . . . . . . . . 29
5.4 The Crystal shapes frame box in the Style tab of the Side Panel . . . . . . . . . . 30
6.1 New Data dialog box . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
vii
6.13 Volumetric data page in the Edit Data dialog box . . . . . . . . . . . . . . . . . . 48
6.14 A dialog box to choose operations for volumetric data . . . . . . . . . . . . . . . 49
6.15 Distributions of electron densities and effective spin densities in an O2 molecule . 50
6.16 Crystal shape page in the Edit Data dialog box . . . . . . . . . . . . . . . . . . . . 51
6.17 Crystal morphologies of anatase-type TiO2 crystals . . . . . . . . . . . . . . . . . 51
6.18 The Objects tab of the Side Panel showing a list of crystallographic faces . . . . . 52
6.19 Crystal morphologies of anatase having faces that do not follow symmetry operations 53
7.1 Phase page in the Edit Data dialog box . . . . . . . . . . . . . . . . . . . . . . . . 54
7.2 Molecules of L-aspartic acid and 2,4,6-trichlorobenzoic acid on the (1014) surface
of calcite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
7.3 A moire pattern of graphene on the (111) surface of Ir . . . . . . . . . . . . . . . 56
7.4 Crystal structures of Ca-olivine (Ca2SiO4) and isostrucural monticellite (CaMgSiO4)
overlapped with each other. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
8.1 Bonds dialog . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
11.2 Calculation of a dihedral angle for four carbon atoms in an aromatic ring in
3-[4-(dimethylamino)phenyl]-1-(2-hydroxyphenyl)-prop-2-en-1-one . . . . . . . . . 80
11.3 Calculation of an interfacial angle between (100) and (511) faces in quartz. . . . 81
12.1 General page in the Properties dialog box . . . . . . . . . . . . . . . . . . . . . . 83
light enabled . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
viii
12.6 Bonds page in the Properties dialog . . . . . . . . . . . . . . . . . . . . . . . . . 87
12.7 Polyhedra page in the Properties dialog . . . . . . . . . . . . . . . . . . . . . . . 88
12.8 Isosurfaces page in the Properties dialog box . . . . . . . . . . . . . . . . . . . . . 89
12.9 Comparison between the two modes of rendering isosurfaces . . . . . . . . . . . . 90
4 91
12.10The 64a1g orbital for the {Cd(S4Mo3(Hnta)3)2} ion with a ball-and-stick model
12.11Composite images of isosurfaces and a ball-and-stick model for albatrossene . . . 91
12.12Relations between orientation of isosurface and the two parameters for opacity
setting . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
12.13Sections page in the Properties dialog box . . . . . . . . . . . . . . . . . . . . . . 93
12.14Distribution of nuclear densities obtained for KOD at 580 K by MEM from single-
crystal neutron diffraction data . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
12.15Relations among saturation levels, data values, and colors of sections . . . . . . . 94
12.16The Objects tab of the Side Panel . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
12.17A list of bond specifications in the Objects tab of the Side Panel . . . . . . . . . . 95
12.18A list of polyhedra in the Objects tab of the Side Panel . . . . . . . . . . . . . . . 95
12.19A list of slices in the Objects tab of the Side Panel . . . . . . . . . . . . . . . . . . 96
12.20A list of forms and faces of crystal morphologies in the Objects tab of the Side Panel 96
13.1 Overall Appearance dialog box . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
13.2 Crystal structure of milarite drawn with light effects (A) enabled and (B) disabled
dialog box . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
13.3 Crystal structure of post-perovskite, MgSiO3, rendered by parallel and perspec-
tive projections . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
13.4 Crystal structure of mordenite rendered with and without depth-cueing . . . . . 100
14.1 Equivalent Positions dialog box . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
14.2 The Equivalent Positions dialog box showing magnetic data . . . . . . . . . . . . . 102
14.3 Geometrical Parameters dialog box . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
14.4 Simulation of an X-ray powder diffraction pattern of YBa2Cu4O8 . . . . . . . . . 109
14.5 Structure Factors dialog box . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
14.6 Fourier Synthesis dialog box . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
15.1 2D Data Display window . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 117
15.2 A dialog box for three different modes to draw a 2D image . . . . . . . . . . . . . 120
15.3 Three modes of creating a 2D image . . . . . . . . . . . . . . . . . . . . . . . . . 121
15.4 Schematic image of the relation between the center of the projection and the
center of the slice . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 122
15.5 Bird’s eye view of a (001) slice of rutile-type TiO2 . . . . . . . . . . . . . . . . . 123
15.6 Electron-density distribution on the (001) plane in rutile-type TiO2 with and
without grid edges . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 123
15.7 Contours page in the Side Panel . . . . . . . . . . . . . . . . . . . . . . . . . . . 124
16.1 Preferences dialog box . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125
ix
List of Tables
x
LICENSE AGREEMENT
VESTA LICENSE
Version 3
Copyright © 2006–2019, Koichi Momma and Fujio Izumi
VESTA is currently copyrighted with its source code not open to the public; we wish to
control the development and future of VESTA by ourselves. This software is distributed free
of charge for academic, scientific, educational, and non-commercial users. Users belonging
to commercial enterprises may also use this software at no cost until a license for business
users is established.
Permission to use this software is hereby granted under the following conditions:
1. Drawings produced by VESTA may be used in any publications (original and review
articles) provided that its use is explicitly acknowledged. A suitable reference for
VESTA is:
K. Momma and F. Izumi, “VESTA 3 for three-dimensional visualization of crystal, volu-
metric and morphology data,” J. Appl. Crystallogr., 44, 1272–1276 (2011).
2. You should not redistribute any copy of the distributed files including MADEL and
STRUCTURE TIDY unless you have a written permission from us.
Part of these terms may be changed without any prior announcement. This software is
provided “as is” without any expressed or implied warranty.
xi
Chapter 1
INTRODUCTION AND
BACKGROUND
1. VICS (VIsualization of Crystal Structures) for displaying and manipulating crystal struc-
tures [15, 16],
1
3. PRIMA (PRactice Iterative MEM Analyses) for MEM analysis from X-ray and neutron
diffraction data [5],
4. Dysnomia [7, 8, 9] which is a successor to PRIMA,
5. ALBA (After Le Bail Analysis) for the maximum-entropy Patterson method [17],
6. Alchemy: A file converter to make it possible to analyze observed structure factors and
their estimated standard uncertainty, which result from Rietveld analysis using GSAS
[18] and FullProf [19], by MEM with Dysnomia [7, 8, 9].
VICS and VEND are programs for 3D visualization of crystal and electronic structures,
2
respectively. They were written in the ANSI C language with the full use of the OpenGL
technology. The graphical user interfaces (GUIs) of VICS and VEND were constructed by
using GLUT and GLUI libraries.
VICS and VEND saw the light of day at the end of 2002, and since then, continued their
growth to be used widely in a variety of studies. However, we never get full satisfaction from
their usability and performance. First, the combined use of VICS and VEND to visualize both
crystal and electronic structures via text files is rather troublesome; on-the-fly visualization of
these two kinds of images is highly desired. Second, their GUIs are not very user-friendly
because they are based on the old-fashioned toolkits, GLUT and GLUI, which have been no
longer upgraded. Above all things, they lack scalability and require large system resources
owing to unrefined programming.
To overcome the faults described above, we at first upgraded VICS to VICS-II employing
3
a modern C++ GUI framework wxWidgets [20]. Later, we further integrated VICS-II and
VEND into VESTA.
• Superimposition of multiple structural models, volumetric data, and crystal faces on the
same Graphic Area
• Visualization of isosurfaces with multiple levels
• Output information about principal axes and mean square displacements for anisotropic
atomic displacements
2
http://www.opengl.org/
3 http://www.wxwidgets.org/
2
• Determination of the best plane for selected atoms
The text box supports four types of basic arithmetic operations: +, , , and /, which means
that you can input, for example, “1/3” instead of “0.333333.” Pressing the <Tab> or <Enter>
key will focus a next control, i.e., text box, button, ratio button, or check box, in a dialog box.
Press <Shift> + <Tab> to focus the preceding control.
3
Chapter 2
4
A B
Figure 2.1: Screenshots of VESTA running on the three different operating systems: (A) Win-
dows, (B) Linux, and (C) macOS.
5
2+
Figure 2.2: Crystal structure of a variant of mica group minerals, masutomilite K(Li,Al,Mn )3-
[(Si,Al)4O10](F,OH)2 [24].
6
Figure 2.4: Crystal structure of beryl, Be3Al2Si6O18 [26].
7
Figure 2.6: Crystal structure of the tetragonal form of melanophlogite,
14 12 14 12
46SiO2 6M 2M (M = N2, CO2; M = CH4, N2) [28]. Bright-blue and
14 12
pink spheres in cages of the SiO4 framework represent M and M
sites for guest molecules, respectively.
8
2.2.1 A variety of structural information derived on selection of objects
Selection of objects (atoms, bonds, and coordination polyhedra) by clicking with a mouse
pro-vides us with a variety of structural information:
• fractional coordinates,
• site symmetries,
1 http://fujioizumi.verse.jp/rietan/book/book.html#BVP
9
B
Mg
b a
Figure 2.8: Electron-density distribution in MgB2. Four hexagonal unit cells are
3
shown with an isosurface level of 0.11a0 (a0: Bohr radius).
0.5
0.45
0.4
0.35
0.3
0.25
0.2
0.15
0.1
0.05
0.0a0-3
10
b
O
H
0 O a
c
Figure 2.10: Nuclear-density distribution in the
paraelectric phase of KH2PO4.
Figure 2.10 shows isosurfaces of scattering-length densities determined from neutron pow-
der diffraction data of KH PO (paraelectric phase, space group: ) at room temperature by
2 4 I42d
MEM-based pattern fitting [43]. In this way, two different colors are assigned to positive and
3
negative isosurfaces. Blue isosurfaces (density: 2:5 fm/Å ) for H atoms are elongated toward
3
yellow ones (density: 2.5 fm/Å ) for O atoms because of double minimum potential with an
H–H distance of about 0.36 Å.
Crystal morphologies can be drawn by inputting Miller indices of faces (see section 6.5). Crystal
faces can be overlapped with structural models as well as isosurfaces of volumetric data, as
exemplified by Fig. 2.11. Thus, VESTA facilitates understanding of the relationship between
crystal morphologies and chemical bonds, for example, by application of the periodic bond chain
(PBC) theory [44]. For example, Fig. 2.12 illustrates an external morphology of anatase-type TiO 2
superimposed on its ball-and-stick model. According to Hartman [44], the f101g face of anatase is
an F face, where two PBCs are parallel to the plane. All the other fh0lg faces are S faces, where
only one PBC is parallel to the plane. The PBC along [010] is intuitively recognized by interactive
rotation of the model in three dimensions.
Morphologies of twinned crystals or epitaxial intergrowths of two or more crystals are vi-
sualized by inputting multiple crystal data (Fig. 2.13). After morphological data have been
given, information about Miller indices, distances from the center of the crystal to faces, and
the surface area of each face are output to the text area.
11
Figure 2.11: An external morphology of an Al2O3 crystal
composed of {001}, {110}, and {113} faces.
12
A B C
Figure 2.13: Crystal morphologies of adularia, a variant of orthoclase (KAlSi 3O8). (A) A single
crystal. (B) A Hypothetical morphology composed of four individual crystals twinned after the
Babeno law. (C) An idealized morphology of repeated twins, which were found from Kobushi
mine, Japan, composed of eight individual crystals twinned after the Babeno and Manebach
laws.
13
2.7.1 Graphical user interface
The GUI of VESTA is built on top of a cross-platform application framework (toolkit), wxWid-
gets, written in the C++ language. The wxWidgets libraries are one of the best toolkits for
cross-platform GUI programming. It provides us with a consistent look-and-feel inherent in
2
each operating system. The license agreement of wxWidgets, an LGPL-like license with
some exceptions allowing binary distribution without source code and copyright, is flexible
enough to permit us to develop any types of applications incorporating wxWidgets.
2 http://www.wxwidgets.org/about/newlicen.htm
14
Chapter 3
GETTING STARTED
3.2 Windows
Both 32- and 64-bit versions of VESTA are available for Windows. They were tested on
Windows 10, 8, 7, Vista, and XP. VESTA does not run on Windows ME or older.
The Windows versions of VESTA are archived in the zip format. To get the 32-bit version
of VESTA for Windows, download
http://jp-minerals.org/vesta/archives/VESTA.zip
To get the 64-bit version of VESTA for Windows, download
http://jp-minerals.org/vesta/archives/VESTA-win64.zip
Note that the 32-bit version of VESTA can also be used on 64-bit Windows. Extract the whole
contents of the archive file in the same directory. VESTA can be launched in the following
four fashions:
1. Double click the icon of VESTA.exe.
This command will start VESTA and open specified files, file1, file2, etc. Input files may
be omitted. When the current directory is different from the directory of VESTA, you
1 http://software.intel.com/en-us/articles/intel-sdp-products/
15
must input the absolute path of the executable binary file of VESTA. To execute VESTA
by simply typing VESTA regardless of the current directory, set environment variable
PATH to include absolute path of the directory where VESTA is placed.
3. Drag & drop file icons of supported file types on the icon of VESTA.exe.
To associate an extension with VESTA, right click a file with the extension and select
“Proper-ties.” Then change the “Program” section in the Properties dialog box.
Beware that VESTA cannot be executed unless the zip file is decompressed. By default,
double-clicking a zip file on Windows XP or later will open it with Explore, i.e., the standard
file manager of Windows. However, the contents of the archive file seen in the Explore
window are actually not extracted but just ‘previewed’. In such a case, copy the whole
contents of the ‘previewed’ archive into an appropriate directory, typically in C:\Program
Files. Then, VESTA.exe in the VESTA folder can be executed as described above.
3.3 macOS
The macOS version of VESTA runs on OS X 10.9 or later using Macs equipped with Intel
CPUs that support 64 bit instruction sets. For those using older machines, VESTA v3.4.4 or
earlier runs on OS X 10.6 or later and Intel 64 bit CPU. VESTA 2.1.6 or earlier runs on Macs
(Mac OS X 10.4) equipped with PowerPC or Intel 32 bit CPUs.
The macOS version is contained in VESTA.dmg with a disc image format. Download
http://jp-minerals.org/vesta/archives/VESTA.dmg
and mount its disc image by double-clicking the icon of VESTA.dmg. Copy the whole
contents of the disc image into an appropriate directory, typically /Applications/VESTA.
VESTA can be launched in the following four fashions:
1. Double click the icon of VESTA (or VESTA.app).
3. Drag and drop file icons of supported file types on the icon of VESTA.
3.4 Linux
Both 32- and 64-bit versions are available for Linux platforms.
16
2. Distributions where VESTA is known to work
Linux versions of VESTA are archived in the tar.bz2 format. For 32-bit OSs running on
PCs equipped with Intel x86 CPUs, download
http://jp-minerals.org/vesta/archives/VESTA-i686.tar.bz2
For 64-bit OSs running on PCs equipped with AMD 64 CPUs, download
http://jp-minerals.org/vesta/archives/VESTA-x86_64.tar.bz2
Extract the whole contents of the archive file to a directory. Then, execute VESTA in that
directory by double-clicking the “VESTA” file on a file manager or by typing the following
command in a command line:
VESTA file1 file2 ...
This command will start VESTA to open the specified files, file1, file2, etc. Of course, the
input files may be omitted. When the current directory is different from that of VESTA, you
must input the absolute path of the executable binary file of VESTA. To launch VESTA by
simply typing VESTA regardless of the current directory, create a symbolic link of VESTA in a
directory where environment variable path is set. You can optionally add VESTA to a panel,
dock, or the application menu.
VESTA always uses “.” (period) as decimal points regardless of locale settings of OSs,
even though some languages use “,” (comma) as decimal points. Such behavior is achieved
by internally changing locale settings on execution of the VESTA program. This is necessary
because most of file formats, such as Crystallographic Information Framework (CIF), use
period as decimal points.
VESTA crashes on Ubuntu 6.06 with some locale settings. In this case, try the following
command:
GTK_IM_MODULE=xim VESTA
17
3.5 Trouble Shooting
VESTA may not function properly with part of graphics cards; most of such troubles seem to
arise from bugs of their drivers. Updating drivers to latest ones may solve the problems. If
some troubles are encountered on the use of VESTA, please try to run it on another machine
to check whether or not the trouble is caused by a bug in it.
If you fail in executing VESTA on a Linux platform, try to execute it in a command line.
Then you will probably get some information about the cause of the problem. Most likely, the
failure is caused by missing libraries. In such a case, the names of missing libraries are
output to the terminal.
18
Chapter 4
MAIN WINDOW
1. Menu Bar
2. Horizontal Toolbar
7.Te xtAre a
19
1. Menu bar: “File”, “Edit”, “View”, “Objects”, “Utilities”, and “Help” menus are placed on
the menu bar. The menu bar is placed at the top of the Main Window in Windows and
Linux while the menus are displayed at the top of the screen in macOS.
4. Side Panel: Composed of three tab pages. The Tools tab contains options for GUI tools
in the Vertical Toolbar (see chapter 11). The Style tab contains GUI tools for display
styles (see chapter 5). The Objects tab shows a list of phases and objects in each
phase (see chapter 12.2).
6. Graphics Area: An area where various objects are displayed in three dimensions.
7. Text Area: Two types of texts are displayed here (see section 4.5). A variety of
information related to users’ operations is displayed in the Output tab. In the Comment
tab, comments on data displayed currently can be input.
4.2 Menus
4.2.1 File menu
• New Structure...: Open a Edit Data dialog box to input crystal data of a new structure.
• Save...: Save current data to a file, *.vesta, with the VESTA format. If the current data
has been input from *.vesta or once saved as *.vesta, *.vesta is overwritten. Otherwise,
a file selection dialog box will appear to prompt you to input a new file name.
• Save As...: Save current data to a new file, *.vesta, with the VESTA format. A file
selection dialog box will appear to prompt you to input a new file name.
• Export Data...: Export current data to a file with formats other than the VESTA format.
• Export Raster Image...: Export a graphic image to a file with a raster (pixel-based)
format.
• Export Vector Image...: Export a graphic image to a file with a vector format.
• Save Output Text...: Save text data in the Text Area as a text file.
• Close: Close the current page. When no data are displayed, this menu is practically the
same as the Exit menu.
20
4.2.2 Edit menu
• Undo: Undo the last action.
• Edit Data:
Phase…: Open an Edit Data dialog box to select the Phase page.
Unit Cell…: Open an Edit Data dialog box to select the Unit Cell page.
Structure Parameters…: Open an Edit Data dialog box to select the Structure
parameters page.
Volumetric Data…: Open an Edit Data dialog box to select the Volumetric data
page.
Crystal Shape…: Open an Edit Data dialog box to select the Crystal shape page.
• Bonds...: Open a Bonds dialog box to create or edit data related to bonds.
• Fit to Screen: Scale and center justify objects to fit to the Graphics Area.
• Clear Text Are: Clear all the text in the Text Area.
21
• Volumetric Data: items of this menu can be also selected in the Side Panel.
• Properties
• Site Potentials and Madelung Energy...: Calculate site potentials and the Madelung
energy of the currently displayed crystal by an external program, MADEL.
• Powder Diffraction Pattern...: Simulate a powder diffraction pattern with RIETAN-FP
[12] and display the result with an external graph-plotting program.
• Structure Factors...: Calculate structure factors.
• Patterson Densities:
22
From Volumetric Data…: Calculate and display density distributions of Patterson
function from the currently displayed volumetric data.
• Line Profile...: Calculate a line profile of volumetric data between two positions and
output them in a text file.
2
• Manual...: Open the user’s manual of VESTA in a PDF viewer.
3
• Check for Updates...: Open the web page of VESTA with a browser to check whether
or not a new version of VESTA has been released.
4
• About VESTA... : Show information about VESTA.
4.3.1 Alignment
The above six buttons are used to align objects parallel to the a, b, or c axis, or parallel to the
a , b , or c axis, respectively.
2
The present PDF file, VESTA_Manual.pdf, need to share the same folder with a binary executable file of
VESTA.
3
http://jp-minerals.org/vesta/
4 In the macOS version, this item is placed under the VESTA menu.
23
4.3.2 Rotation
These six buttons are used to rotate objects around the x, y, or z axis. The step width of
rotation (in degrees) is specified in the text box next to the sixth button:
4.3.3 Translation
Translate upward
Translate downward
Translate leftward
Translate rightward
These four buttons are used to translate objects upward, downward, leftward, and rightward,
respectively. The step width of translation (in pixels) is specified in the text box next to the
fourth button:
4.3.4 Scaling
Zoom in
Zoom out
These three buttons are used to change object sizes. The step width of zooming (in %) is
specified in a text box next to the third button:
24
4.4 Tools in the Vertical Toolbar
Rotate
Select
Translate
Magnify
Measure distance
Measure angle
is displayed in the Text Area. SSE2 means that Streaming SIMD Extensions 2 is supported
in this graphics card. In the case of a Power Mac G5 (Dual 2.5 GHz) equipped with ATI
Radeon 9600 XT, a message
Video configuration: ATI Radeon 9600 XT OpenGL Engine
is issued in the Text Area. If the GPU of your PC does not support any hardware acceleration of
OpenGL, Video configuration would be GDI Generic on Windows, and Software Rasterizer
or Mesa GLX Indirect on Linux.
When a new file is loaded, a summary of data appears in the Output area, including the
absolute path of the file and the title of the data. For crystal-structure data, lattice parameters, a
unit-cell volume, and structure parameters (element names, site names, fractional coordinates,
occupancies, isotropic atomic displacement parameters, multiplicities plus Wyckoff letters, and
site symmetries) are output. In the case of volumetric data, lattice parameters, number of grids
along each axis, and number of polygons and unique vertices on slices and isosurfaces are
displayed. For instance, a CIF of PbSO4 affords the following output:
Title O4 Pb1 S1
Lattice type P
25
Space group number 62
Setting number 1
Lattice parameters
Structure parameters
Such a type of a list is also output when pressing [OK] button in the Edit Data dialog box.
26
Chapter 5
DISPLAY STYLES
Display styles of structural models, volumetric data, and crystal shapes are controlled either
from the Style tab of the Side Panel or from the “Objects” menu. Both of the methods allow
us to use the same options.
This option controls the visibility of a structural model. When Figure 5.1: The Structural
this option is checked (default), a structural model is visible; model frame box in the
otherwise, no structural model is shown. Uncheck this option Style tab of the Side Panel.
when you want to see only isosurfaces and sections for data
containing both structural and volumetric ones.
27
5.1.2 Styles
As described in 2.2, VESTA represents crystal structures by the five different styles: ball-and-
stick, space-filling, polyhedral, wireframe, and stick models. When atoms are drawn as
spheres, they are rendered with radii corresponding to 40 % of actual atomic radii in all but
the space-filling model, where atoms are rendered with the actual atomic radii. Default radii
of atoms are selected from three types: atomic, ionic, and van der Waals radii. The radius of
each element and a type of radii are specified at the Atoms page in the Properties dialog box
(see section 12.1.2). Features in each structural model are described below with screenshots
of the structure for quartz [49] on the right side.
“Ball-and-stick”
In the “Ball-and-stick” model, all the atoms are expressed as solid
spheres or displacement ellipsoids. Bonds are expressed as either
cylinders or lines.
“Space-filling”
In the “Space-filling” model, atoms are drawn as interpenetrating
solid spheres, with radii specified at the Atoms tab in the
Properties dialog box (see section 12.1.2). This model is useful for
understanding how atoms are packed together in the structure.
“Polyhedral”
In the polyhedral model, crystal structures are represented by co-
ordination polyhedra where central atoms, bonds, and apex atoms
may also be included. Bonds between central and apex atoms
have to be searched with the Bonds dialog box to display
coordination polyhedra comprising them. Atoms are expressed as
solid spheres or displacement ellipsoids. Bonds are expressed as
either cylinders or lines. Needless to say, the transparency of the
coordination polyhe-dra must be high enough to make it possible
to see the central atoms and bonds. One of six different styles for
representing polyhedra is specified at the Polyhedra tab in the
Properties dialog box (see section 12.1.4).
28
“Wireframe”
In the “Wireframe” model, atoms having no bonds are drawn as
wire-frame spheres whereas those bonded to other atoms are
never drawn. All the bonds are presented as lines with gradient
colors. This model is useful for seeing and manipulating complex
and/or large structures because this is usually the fastest model for
rendering structures on the screen.
“Stick”
In the “Stick” model, atoms with no bonds are drawn as solid spheres
while atoms bonded to other atoms are never drawn. All the bonds
are expressed as cylinders, whose properties can be changed at the
Bonds tab in the Properties dialog box (see section 12.1.3). This
model serves to see frameworks or molecular geometry.
Displacement ellipsoids
In the ball-and-stick and polyhedral models, atoms can be
rendered as displacement ellipsoids. The probability for atoms to
be included in the ellipsoids is also specified in the Properties
dialog box (see section 12.1.2).
• Show sections
This option controls whether or not sections of isosurfaces
are visible. When this option is checked (default), sections Figure 5.3: The Volumetric
are visible, otherwise they are not shown. data frame box in the Style
tab of the Side Panel.
• Show isosurfaces
This option controls whether or not isosurfaces of volu-
metric data are visible. When this option is checked (default), isosurfaces are visible,
otherwise they are not shown. To see only a structural model for data containing both
29
structural and volumetric ones, uncheck this option in addition to the “Show sections”
option.
• Surface coloring
“Surface coloring” means that colors of isosurfaces drawn from one data set are
determined by a secondary data set. A typical example is coloring of electron-density
isosurfaces on the basis of electrostatic potentials (see Fig. 12.9). This option is
enabled only when the secondary volumetric data for surface coloring has been loaded
with the Volumetric data tab in the Edit Data dialog box (see 6.4.2).
• Smooth shading
• Wireframe
• Dot surface
Isosurfaces are drawn as solid surfaces with variable opacity in the Smooth shading mode
whereas isosurfaces are represented by lines and points, respectively, in the Wireframe and
Dot surface modes.
5.3.1 Styles
“Unicolor”
All the faces of a crystal are filled with a single color, as the right
figure (quartz) illustrates.
30
“Custom color”
Each face is filled with a color assigned to symmetrically
equivalent fhklg faces specified in the Edit Data dialog box (see
section 6.5) or in the Side Panel.
“Wireframe”
Only edges of a crystal shape are drawn with lines.
31
Chapter 6
To create a new structure, choose the “File” menu ) “New Structure…”. To edit current data,
choose “Edit” menu ) “Edit Data” ) “Phase…”. In both of the cases, the same dialog box
named New Data or Edit Data appears (Fig. 6.1). This dialog box consists of the following
five tab pages:
• Phase
• Unit cell
• Structure parameters
• Volumetric data
• Crystal shape
At the top of the dialog box, a serial number and a title are displayed for the selected phase.
32
6.1 Defining Phases
In the Phase tab (Fig. 6.1) of the Edit Data dialog box, you can add, delete, copy, or import
phase data that are overlaid on the same Graphics Area. To edit the title of a phase, select a
row in the list and click on the Title column. When dealing with a single phase, you have
nothing to do with this tab except for editing the title. The Positioning and Orientation frame
boxes are used when superimposing two or more phases. See chapter 7 for how to visualize
multiple phase data in the same Graphics Area.
Figure 6.2: Unit cell page in the Edit Data dialog box.
• Custom
33
• Triclinic
• Monoclinic
• Orthorhombic
• Tetragonal
• Trigonal
• Hexagonal
• Cubic
• Rhombohedral
34
Table 6.1: Non-standard settings in two triclinic space groups.
Space groups
Axis choice No. 48, 50, 59, 68, and 70 Others
Origin choice 1 Origin choice 2
abc 1 2 1
ba c 3 4 2
cab 5 6 3
c ba 7 8 4
bca 9 10 5
ac b 11 12 6
35
6.2.2 Behavior when changing a space-group setting
The list box below the Setting controls the behavior when changing a space-group setting.
There are two options:
• Update structure parameters to keep 3D geometry
On the use of the first option (default), lattice parameters and fractional coordinates are con-
verted to keep the structural geometry of a crystal when the setting number of a space group
is changed or a transformation matrix is specified. When the second option is used, lattice
parameters and fractional coordinates remain unchanged instead of the geometry.
If a correct structure is displayed in the Graphic Area but another space-group setting is
preferred to the current one, choose Update structure parameters to keep 3D geometry and
change the space-group setting. The second option is typically used to change a setting
number when the correct setting number is not recognized by VESTA after reading in a file
such as CIF. In such a case, lattice parameters and fractional coordinates are correct
whereas symmetry operations are incorrect. Then, choose Keep structure parameters
unchanged and select the correct setting number. When a setting number cannot be uniquely
determined from data recorded in a structural-data file, VESTA assumes setting number 2 (a
setting with the inversion center at the origin) in centrosymmetric space groups with two
origin choices, and setting number 1 in all the other space groups.
Note that even if structure parameters remain unchanged, some lattice parameters may be
changed in conformity with constraints imposed on them in the current space-group setting.
36
Figure 6.3: Equivalent Positions dialog box with an editing mode.
be removed by clicking the [Remove duplicate atoms] button in the Structure parameters tab
(see page 48).
37
Non-conventional space-group settings
Note that the orientation of the c axis is inverted in this case, or otherwise the transformed
structure becomes enantiomorphic with respect to the original structure in non-
centrosymmetric space groups. In general, when the determinant of P , det(P ), is negative,
the crystal-coordinate system is transformed from right-handed to left-handed, and vice
versa. When such a matrix is specified as P , a warning dialog box will appear to confirm
whether or not you proceed to the next step.
′
The unit-cell volume, V , is transformed to V by
′
V = det(P )V: (6.2)
P P P p2C
B 32 33
C : (6.3) p3A
@ 0 0 0
1
38
Then each of the following crystallographic variables, i.e., primitive translation vectors a, b, c,
reciprocal lattice vectors a , b , c , metric tensor G, metric tensor of reciprocal lattice G , atomic
coordinate X, symmetry operation W, and anisotropic atomic displacement parameters
, are transformed as described below.
Primitive translation vectors a, b and c are transformed by P as
′ ′ ′
(a ; b ; c ) = (a; b; c)P
0 1
P P P
B 11 12 13C
(6.4)
= (a; b; c) @P21 P22 P23A
P P P
31 32 33
The shift of the origin, i.e., (0, 0, 0), is defined with a shift vector t as
t = (a; b; c)p
0 1
p
B 1C (6.5)
= (a; b; c) @p2A
p
3
= p1a + p2b + p3c:
Let A be the representation matrix of primitive translation vectors, a, b, and c, on the basis
of a set of orthonormal vectors, x, y, and z:
A = (a; b; c) = 0a b c 1 : (6.6)
x x x
@ A
ay by cy
B b C
az z cz
Then the transformation by Eq. (6.4) can be rewritten as
′
A = AP : (6.7)
Alternatively, the use of the metric tensor, G, is often preferable to the direct use of A
owing to the unique definition of G. On the other hand, the number of representations for A is
infinite because the orthonormal bases, x, y, and z, may be arbitrarily chosen. The metric
tensor (matrix) of the direct lattice is defined as
t 0a a a b a c 1 0 ab cos ac cos1 (6.8)
G=AA= a2
2
b a b b b c = ba cos b bc cos ;
cb c
@ A @ A
B C B cb cos C
ca c ca cos c2
′
where t denotes transposition. Then, G is transformed into G by
′ ′ t ′ t t t t (6.9)
G = (A ) A = (AP ) (AP ) = P (A A)P = P GP :
Similarly, the representation matrix of the reciprocal basis vectors a , b , and c is defined as
A = (a ; b ; c ) = Ba bx cx C (6.10)
x
0 ay by cy 1 :
@a bz cz A
z
39
The direct and reciprocal basis vectors are related by
1 t t 1 (6.11)
A = (A ) = (A )
and
t t t t (6.12)
(A ) A = A(A ) = A A = A A = E;
where E is the identity matrix:
B 1 0 0C (6.13)
0 0 1
0
E= 0 1 01 :
@ A
b c cos A:
B C B c b cos
c a c b c c c a cos
(6.14)
′
The transformation of G into G is represented by
′ ′ t ′ 1 t 1 1 t 1 t t (6.15)
G = (A ) A = (AP ) ((AP ) ) = P ((A ) A )(P ) = QG Q ;
1
where Q = P .
The coordinates x, y, and z in the direct space are transformed by the 4 4 transformation
matrix Q as
0 y′′ 1 0 y1
Bx C B xC
′ = Q z ;
z (6.16)
B 1 C B 1 C
BC BC
@A @A
where
( )
Q= 1 1
P P p
o 1
( )
= Qqo
1 1 (6.17)
0 q
1C
q
B 11 2 CC
Q Q Q
= BB 21 22 23 @Q31 q3A
Q12 Q13 Q
Q32Q33 1
0 0 0
1
=P :
′
X = QX: (6.18)
40
a symmetrically equivalent point X with the coordinates, x~, y~, and z~, by
B x~ f W W W w CBxC
11 12 13
1
1 0 W W W
0 1 0 y1
1 C B 0 0 0 1 1
y~ 21 22 23 w2
BC B W W W CB C
B z~C =
B 31 32 33 w3 C B z C
@A @ A@ A (6.19)
0 1
x+W y+W z+w
B 11 12 13 1C
= BBW 21
x + W y + W z + w
22 23 2CC ;
@W31x + W32y + W33z + w3A
W
1
W = QWP:
′ (6.23)
t
Let Q be the transposed matrix of Q defined in Eq. (6.17). Then, anisotropic atomic
displacement parameters, , are transformed by
0 1
′ ′ ′ ′
B 11 12 13 C = @
′ ′ ′
12 22 23 A
′ ′ ′
13 23 33
@Q Q Q A@ 11 12 13A @Q Q Q (6.24)
11 12 13 11 21 31
Q Q Q Q Q Q
= 0
21 22 23 10
12 22 23 10
12 22 32A 1
BQ 13 23 33C
Q Q CB BQ Q Q C
31 32 33 13 23 33
t
= QQ:
41
Creation of superstructures and substructures
The transformation matrix is also used to convert (primitive lattice)⇄(complex lattice) (Fig.
6.6), and to create superstructures and averaged substructures. To convert a structure, you
have to select mode Update structure parameters to keep 3D geometry in the Unit cell tab of
the Edit Data dialog box before opening the Unit Cell Transformation dialog box.
When the determinant of P , det(P ) is larger than 1, the transformed unit-cell becomes
larger than the original one, and when det(P ) is smaller than 1, the transformed unit-dell
becomes smaller than the original one. In such cases, a dialog box appears to confirm
whether you intend to change the unit cell volume (Fig. 6.7).
When a rotation matrix with det(P ) > 1 is specified, or the transformed unit cell contains
additional lattice points even when det(P ) = 1, a dialog box appears to ask you how to
convert the structure (Fig. 6.8). There are the following three options:
• “Add new equivalent positions to a list of symmetry operations”
• “Search atoms in the new unit cell and add them as new sites”
• “Do nothing”
In the first and second options, the structural geometry of the original data will remain
unchanged. On the other hand, on selection of the third option, some atoms in the new unit
cell will be missing. On selection of the first option, a superstructure will be created by
′
examining the following equation to find additional symmetry operations W having a
translational component in between (0, 0, 0) and (1, 1, 1):
0 W ′ W ′ W ′ w′1 0 W W W 1
11 12 13 1 11 12 13 w +n
1 xC
′ ′ ′ ′ ′ W W W w +n
W = B
W21
W31′
W22
W32′
W23
W33′
w2
w 3′ C
=Q B W
31
21
W
32
22
W
33
23
2 yCC P;
B C B
@ A @ w3 + nz A (6.25) 1
B 0 0 0 1 C B 0 0 0 1; 2; : : : ):
(nx; ny; nz = 0;
42
On selection of the second option, a superstructure will be created by examining the
following equation to find additional sites lying in between (0, 0, 0) and (1, 1, 1):
01 0 1
′
x x + nx
B y ′ C = Q B y + ny C : (nx; ny; nz = 0; 1; 2; : : : ) (6.26)
1 1
B z′ C B z + nz C
B C B C
@A @ A
Note that the transformed new unit-cell should be symmetrically consistent with the space
group currently in use. Otherwise, reduce symmetry to P 1 by clicking [Remove symmetry]
before creating a superstructure, and input a new space group after the superstructure has
been created. If det(P ) < 1, the same position may result from two or more sites. In such a
case, the redundant sites can be removed by clicking [Remove duplicate atoms] button in the
Structure parameters tab (see page 48).
Figure 6.9: The Unit cell page in the Edit Data dialog box, with
the “Magnetic structure” option checked.
43
6.3 Structure Parameters
The Structure parameters tab in the Edit Data dialog box (Fig. 6.10) is used to input and edit
parameters related to atomistic structures. The following structural information is input or
edited by using text boxes at the upper half of the page:
• [Symbol...]: Symbol of an element (up to two characters).
Figure 6.10: Structure parameter page in the Edit Data dialog box.
A list of sites in the asymmetric unit is shown at the lower half of the page. Text boxes are
disabled when no site is selected in the list. When an item in the list is selected, the text
boxes are updated to show data relevant to the selected site. In-place editing of a list item is
also supported, as Fig. 6.10 shows.
44
To add a new site, click the [New] button at first, and edit texts in the text boxes. A
modification in a text box is applied to the selected site immediately after the focus (caret)
has moved to another text box or another GUI control. In the case of text boxes for numerical
values, modifications will be discarded if an input text is not a valid numerical value, or if the
text box is empty. Therefore, if you have accidentally changed a value in a text box, the
original value can be restored by erasing the entire text in the text box before moving the
focus to another GUI control.
To delete a site, select an item in the list, and click the [Delete] button or press the
<Delete> key. Clicking the [Clear] button will delete all the sites in the list.
45
Figure 6.12: Crystal structure of Cs 6C60 [53], with C60 rep-
resented by translucent polyhedra. A virtual site X with an
occupancy of g = 0 was added at the center of C 60. Polyhedra of
C60 were created by specifying X–C bonds in the Bonds dia-log
box. C–C bonds are also searched by a “Search molecules”
mode to represent polyhedral edges by solid cylinders. See
chapter 8 for details in creating bonds and polyhedra.
6.3.4 Occupancy
The occupancy is unity for full occupation and zero for virtual sites. If the occupancy of a site
is less than unity, atoms occupying there are displayed as circle graphs for occupancies. If
unicolor balls are preferred to bicolor ones, change the occupancy to unity for convenience.
If the occupancy of a site is zero, the site is treated as a virtual site that is not occupied by
any atoms in the real structure. Neither the virtual sites nor bonds connecting them are
displayed on the screen. Virtual sites are used to visualize, for example, cage structures of
porous crystals as solid polyhedra without creating a large number of unnecessary bonds
(Fig. 6.12).
The Debye–Waller factor, Tj (h), which is often referred to as the temperature factor, is
included in formulae for structure factors, F (h), to represent the effect of static and dynamic
displacement of atom j (see Eqs. 14.8 and 14.11). The displacement of atom is formulated in
two different ways: anisotropic and isotropic atomic displacement.
Anisotropic models
The anisotropic Debye–Waller factor, T (h), is defined as )]
[
= exp h (
2 2 2 2 2 2 2
T (h) = exp 2 h a U11 + k b U22 + l c U33 + 2hka b U12 + 2hla c U13 + 2klb c U23
[ ( 2+ k2 +l
2 + 2hk + 2hl + 2kl )] ;
11 22 33 12 13 23
(6.27)
46
where a ; b ; c ; ; , and are reciprocal lattice parameters.
The type of anisotropic atomic displacement parameters must be specified in list box
{Anisotropic}: ‘None’ to omit anisotropic atomic displacement parameters, ‘U’ to input Uij (U11,
U22, U33, U12, U13, and U23), and ‘beta’ to input ij ( 11, 22, 33, 12, 13, and 23). Using this list box, we
can convert Uij into ij and vice versa.
Anisotropic atomic displacement parameters defined in different ways must be converted
into Uij or ij defined above.
On refinement of ij ’s by the Rietveld method with RIETAN-FP [12], ij’s for each site have
to satisfy linear constraints imposed on them [54, 55].
Isotropic models
[ ( sin 2 ]
T (h) = exp B
)
( )
= exp B (6.28)
2
4d
= exp
( 2
2 U
) 2
;d
where B and U are the isotropic atomic displacement parameters (U = B/8 2), is the Bragg
angle, is the X-ray or neutron wavelength, and d is the lattice-plane spacing. B and U are
related to the mean square displacement, u2 , along the direction perpendicular to the reflection
plane with (6.29)
2 2
B=8 U=8 u2 :
For example, a B value of 3 Å 2 corresponds to a displacement of about 0.2 Å.
The type of isotropic atomic displacement parameters must be specified in list box {Isotropic}.
Using this list box, we can convert U into B and vice versa.
Isotropic atomic displacement parameters defined in different ways must be converted
into B or U defined above.
In crystallographic sites for which anisotropic atomic displacement parameters, Uij or ij,
have been input, equivalent isotropic atomic displacement parameters, Beq or Ueq, are
calculated from Uij, a , b , c , and the metric tensor G [56, 57]:
B eq = 8 2U eq
(6.30)
with
Ueq = 1 ∑ ∑ Uija a ai a
j
2 2
i j
1 2
3[ i j
(6.31)
= 3 U11 (aa ) + U22 (bb ) + U33 (cc )
+ 2U12a b ab cos ]
+ 2U13a c ac cos + 2U23b c bc cos :
Beq and Ueq are regarded as B and U in the Structure parameters tab of the Edit Data dialog
box after reopening it.
47
Option “Link” specifies the manner of outputting current structure data in a file, *.vesta, with
the VESTA format. Structure data are usually recorded in *.vesta. On the other hand, for
volumetric data, relative paths (including a file name) of data files are recorded in *.vesta instead
of recording the volumetric data directly in it. When option “Link” is checked, structural data are
also recorded as the relative path of the data file. Then, even if the structural-data file is changed
after saving *.vesta, the changes are reflected in VESTA when *.vesta is reopened. This option is
useful in preparing a file, *.vesta, with the VESTA format as a template file for setting objects and
overall appearances. This option disables text boxes and list boxes in the Structure parameters
tab to prohibit users from editing structure data because any modifications to structure data will
not be saved to a file when this option is enableld.
Figure 6.13: Volumetric data page in the Edit Data dialog box.
48
6.4.1 Volumetric data to draw isosurfaces
VESTA enables us to deal with more than two volumetric data sets. Clock [Import...] to select
a file in the file selection dialog box. Only files with extensions of supported formats are
visible in the file selection dialog box. After selecting a volumetric data file, a Choose
operations dialog box appears (Fig. 6.14).
Operation
In the Operation radio box, one of the following five data operations can be selected in
addition to conversion of data units.
“Multiply to current data” is convenient when squaring wave functions to obtain electron
densities (existing probabilities for electrons). The last two options are not displayed when no
volumetric data are contained in the current data. In such a case, the first and third options
have also the same effect.
49
A B
Figure 6.15: Distributions of electron densities and effective spin densities in an O2 molecule.
(A) up-spin electron densities, ", (B) down-spin electron densities, #, and (C) effective spin
densities, ∆ = " ", calculated with VESTA. Both " and # were calculated with DVSCAT [58].
3 3
Isosurface levels were set at 0.01a0 in (A) and (B), and at 0.001a0 in (C). where a0 is the
Bohr radius.
Interpolation
The spacial resolution of volumetric data can be increased using the algorithm of cubic spline
interpolation. The interpolation level is specified in the {Interpolation} text box.
50
Figure 6.16: Crystal shape page in the Edit Data dialog box.
A B
Figure 6.17: Crystal morphologies of anatase-type TiO 2 crystals. In the inner crystal
in (B), distances from the origin are set at 3.5 Å for both {101} and {103}. In the outer
crystal, distances from the origin are set at 3.8 Å for {101} and 5 Å for {103}.
51
6.5.1 Examples
52
A B C
Figure 6.19: Crystal morphologies of anatase having faces that do not follow symmetry opera-
tions. (A) Option “Apply symmetry operations” is unchecked for a (103) face. (B) The distance
from the origin was set at 3.8 Å for {101} and 5 Å for {103}, and a (103) lattice plane was
inserted by setting a distance from the origin at 3.8 Å. (C) A (103) face, which overlaps with
the (103) lattice plane in (B), was further added.
53
Chapter 7
To display multiple-phase data in a single Graphics Area, set positioning and orientation of
each phase after clicking the Phase tab in the Edit Data dialog box (Fig. 7.1). The positioning
and orientation of a phase are specified relative to another phase or the Cartesian coordinate
system that is commonly used as internal representation of all phases. By default, the origin,
(0; 0; 0), of a phase is placed at that of the internal coordinate system. The orientation of a
phase is set such that the [1 0 0] axis of a phase is parallel to the x axis of the internal
coordinate system with the [0 1 0] axis parallel to the x–y plane.
To avoid the circular reference of phases, positioning and orientation must be set in
descend-ing order of the phase list.
d
e f g
54
the internal Cartesian coordinate system. Then, input (x; y; z) coordinates of the selected and
reference layers; they are exactly overlapped with each other (Fig. 7.1, c).
7.3 Examples
Figure 7.2 shows two kinds of organic molecules approaching a surface of calcite-type
CaCO3. Figure 7.3 illustrates graphene on a (111) surface of Ir [59]. An example of
overlaying isostu-ractural crystals is displayed in Fig. 7.4.
Figure 7.2: Molecules of L-aspartic acid and 2,4,6-trichlorobenzoic acid near a (1014) surface of
step.
calcite (CaCO3) having a [441]
55
Figure 7.3: A moire pattern of graphene on the (111) surface of Ir.
56
Chapter 8
To search for bonds and atoms connected by them, and to create coordination polyhedra,
choose the “Edit” menu ) “Bonds…”. At the top of the Bonds dialog box (Fig. 8.1), select a
phase to edit. A list of bond specifications is shown at the lower half of the dialog box. GUI
controls in the Search bonds and atoms frame box are used to edit a bond specification.
They are disabled when no bond specification is selected in the list.
To add a new bond specification, click the [New] button at first, select atoms relevant to the
bond, and edit minimum and maximum lengths. Modifications in the GUI control are applied to the
list either immediately or after the focus (caret) has moved from a text box to another GUI control.
To delete a bond specification, select it in the list, and then click the [Delete] button or press the
<Delete> key. Clicking the [Clear] button deletes all the bond specifications in the list. Press the
[OK] or [Apply] button to reflect editing results in the Graphics Area.
57
8.1 Specifications of Searching for Bonds
8.1.1 Search mode
Bonds and, in turn, atoms connected by them are searched in one of the following three
search modes:
• “Search molecules”
All the pairs of atoms are searched on the basis of user-specified {Min. length} and
{Max. length}. Neither A1 nor A2 is specified in this option. This mode best meets
searching for atoms and bonds in molecular crystals.
8.1.3 Options
• “Search by label”
On selection of this option, a pair of atoms is specified by labels of sites. When this
option is unchecked, a pair of atoms is specified by elemental symbols.
• “Show polyhedra”
This option specifies that A1 atoms are central atoms of coordination polyhedra. This
option has no effect on bond specifications in the “Search molecules” mode because no
coordination polyhedra are searched in that mode.
58
8.1.4 A pair of atoms
Atoms A1 and A2 may be specified by entering either elemental symbols or site names. The {A1}
and {A2} list boxes list elemental symbols of atoms if option “Search by label” is unchecked
(default). On the other hand, {A1} and {A2} list site labels if option “Search by label” is checked.
The minimum and maximum bond lengths, {Min. length} and {Max. length}, are input in
the unit of Å. Though {Min. length} is usually set at zero, it may be positive when dealing with
a disordered structure whose split-atom model gives seemingly very short bonds.
To build up coordination polyhedra, either the “Search A2 bonded to A1” or “Search atoms
bonded to A1” mode should be selected. The central atoms of the coordination polyhedra
must be specified as A1. In the “Search molecules” mode, no coordination polyhedra are
created because VESTA does not have any information about central atoms. In the first two
bond search modes, the Boundary mode 2 (“Search additional atoms if A1 is included in the
boundary”) enables us to search for all the A2 atoms bonded to A1 atoms so that no
coordination polyhedra are truncated even if A2 atoms lie outside the boundary. A1 atoms
lying outside the boundary are never searched.
To search for additional atoms bonded to A2 atoms, use the Boundary mode 3 (“Search
additional atoms recursively if either A1 or A2 is visible”). For example, if hydrogen atoms are
coordinated to some A2 atoms, those hydrogen atoms can be searched even if A2 atoms lie
outside the boundary (Fig. 8.2).
A B
59
A B
Figure 8.3: Crystal structure of anthraquinone [61]. Bonds are searched in the “Search mole-
cules” mode with the boundary modes (A) “Do not search atoms beyond the boundary” and
(B) “Search additional atoms recursively if either A1 or A2 is visible”. In both cases, the
bounding box has a range of coordinates from (0, 0, 0) to (1, 1, 1).
60
Figure 8.4: Crystal structure of tetrakis(di-4-pyridylsulfane)dinitratocopper(II) [63]. Bonds are
searched in the “Search molecules” mode with a maximum distance of 1.6 Å in combination with
the “Search A2 bonded to A1” mode for Cu–N and S–C bonds. The Boundary mode 3 was used.
because searching within a certain area ensures that all the atoms contained in molecules in
the boundary are searched in cases of periodic structures.
There are exceptions where the combination of such modes is well-suited to search for bonds
in inorganic crystals, e.g., searching for O H bonds. Usually, no bonds are shorter than O H ones;
thus, the use of this mode is safe with a small value of {Max. length}. The maximum bond length
of ca. 1.1 Å leads to a search for all the H atoms bonded to O atoms in the Graphics Area.
61
Figure 8.5: Crystal structure of the tetragonal variant of chibaite repre-
sented by a polyhedral model [64].
62
Chapter 9
ADDITIONAL OBJECTS
Figure 9.1: Vectors dialog box showing a list of atoms in SrFeO 2 [65], with the
atom list filtered by elements. Atoms selected in the dialog box are highlighted in
the underlying Graphics Area displaying the magnetic structure of SrFeO2.
63
Vectors can be assigned to either crystallographic sites or individual atoms in the Graphic
Area. When a vector is attached to a crystallographic site, all of its symmetrically equivalent
atoms will bear the same vector, which is properly rotated in accordance with symmetry opera-
tions applied to the atoms. When a vector is assigned to each atom, only that atom will have the
vector whereas no other symmetrically equivalent atoms will bear the vector. Multiple numbers of
vectors can be attached to atoms and sites, and a vector can be attached to more than two atoms
and sites. The length of arrows, which are vectors rendered in the Graphics Area, is scaled by
moduli of the vectors multiplied by the input value of the {Scale factor for modulus}. For more
details about the scaling of vectors, see 9.1.1 and Fig. 9.4.
Two types of vectors can be specified: {Axial vector} and {Polar vector}. Because these
two differ from each other in the mode of applying symmetry operations, the vector type
affects only vectors attached to crystallographic sites. When a vector is attached to individual
atoms, no symmetry operation is applied to the vector, with no differences between the two
types of the vectors.
Typical examples of axial and polar vectors are magnetic moment, , and electrostatic po-
larization vector, p, respectively. The magnetic moment is induced by tiny electrical currents
resulting from electrons moving around the nucleus of an atom. The direction of the moment vec-
tor is determined by that of the current. For example, counter-clockwise and clockwise currents
generate upward and downward moment vectors, respectively. A mirror plane perpendicular to
the moment vector does not change the current direction, while a mirror plane parallel to the
moment vector inverses the current direction. As a result, magnetic moments are transformed by
a mirror plane as Fig. 9.3A illustrates. In the case of the electrostatic polarization vector,
64
A B
μ p
Figure 9.3: Results of applying a mirror operation to (A) magnetic moments and (B) polar-
ization vectors p
on the other hand, a mirror plane perpendicular to the moment vector inverses the
polarization vector, while a mirror plane parallel to the polarization vector does not change
′
the polarization vector (Fig. 9.3B). To be more specific, p is transformed to p by a 3 3
operator matrix, M, as
′
p = Mp (9.1)
′
while is transformed to as
′
=PM; (9.2)
where P is the determinant of M and referred to as the “parity” of the operator M. Note that
the term P M is a scalar-matrix multiplication.
Aside from conventional symmetry operations M, there is another important symmetry
element for vectors, i.e. time reversal. In view of time reversal, the current direction around
an atom is changed, with a result that the direction of is reversed. Let the time reversal be T
(= 1 or 1), then the transformation of is represented by
′
=TPM: (9.3)
65
A s×μ B s×μ
C s×μ D s×μ
Figure 9.4: Effects of two options for vector styles. Option “Penetrate
atoms”: on (A, B) and off (C, D), and option “Add atom radius to the
vector modulus”: off (A, C) and on (B, D). When the latter option is on,
the length of an arrow appearing above an atom sphere is linearly
scaled by multiplied by the {Scale factor for modulus}, s.
66
Figure 9.5: Lattice Plane dialog box.
When volumetric data are included in the current data, lattice planes are colored according to
volumetric data on the lattice planes. Saturation levels of colors are specified at Sections
page in the Properties dialog box (Fig. 9.6; see 12.1.6). When only structural information is
included in the current data, lattice planes are drawn with colors specified in this dialog box.
To draw lattice planes with specified colors for data having both structural and volumetric
information, volumetric data should be deleted at the Edit Data dialog box (see 6.4).
Material settings of lattice planes, i.e. specular color and shininess are input in the Material
67
frame box. Drawing of edges for lattice planes can be controlled in the Edge frame box.
These settings are common to all the lattice planes.
Figure 9.7: The best plane calculated for an aromatic ring in anthraquinone [61].
68
Chapter 10
DEFINING DRAWING
BOUNDARIES AND VIEW
DIRECTIONS
To change the size of a drawing boundary (box), select the [Boundary] button in the Side Panel,
or choose “Objects” menu ) “Boundary…”. Then, the Boundary dialog box appears as Fig. 10.1
illustrates. This dialog box can also be opened with a keyboard shortcut of <Ctrl> + <Shift>
+ <B>. At the top of the dialog box, select a phase to edit. Press the [OK] or [Apply] button to
reflect editing results in the Graphics Area.
Changing the boundary regenerates all the atoms in the Graphics Area and reset all the
properties of objects to default values. Selected or hidden states of atoms, bonds, and
polyhedra are reset to the default states.
69
Figure 10.2: Electron density distribution determined for D-sorbitol
[66] by the maximum-entropy method from synchrotron X-ray
powder diffraction data. Drawn in ranges from (0, 0, 0) to (1, 1, 1)
with a cutoff plane (1 1 0).
70
Figure 10.3: Crystal structure of diamond drawn in ranges from ( 8; 8; 8) to (8; 8; 8) with and
without a cutoff plane f111g. The drawing boundaries are displayed as transparent faces.
71
Figure 10.4: Orientation dialog box.
hu + kv + lw = 0: (10.1)
This condition must, therefore, be satisfied in order to specify the upward direction on the
screen, otherwise the upward direction on the screen is automatically determined by VESTA.
To view the contents of Graphics Area along basis vectors of a unit cell or a reciprocal cell,
simply press one of the above buttons in the Horizontal Toolbar. When two or more phases
are visualized in the same Graphics Area, the above buttons set the viewing direction relative
to the first phase. To set a viewing direction relative to a phase other than the first one, use
the Orientation dialog box.
72
Chapter 11
INTERACTIVE MANIPULATIONS
11.1 Rotate
In the Rotation mode, mouse behavior can be changed after clicking the Tools tab in the Side
Panel (Fig. 11.1).
“Click” mode
In the “Click” mode, click the left mouse button to rotate the objects. In the “Click” plus [Free]
rotation modes, the rotation axis is perpendicular to a line connecting the clicked position and
the central point. To restrict the rotation axis, select [Around X axis], [Around Y axis], or
[Around Z axis] in the pull down menu of the Rotation modes frame box. To stop rotating the
objects, select either the “Drag” or “Push” mode. The animation speed (intervals between
frames in ms) is specified in the Preferences dialog box.
“Push” mode
In the “Push” mode, press the left mouse button at point 1 and drag the mouse to point 2. In the
[Free rotation] mode, objects are rotated around an axis perpendicular to a line connecting points
1 and 2. The rotation speed is proportional to the speed of moving the mouse. The objects stop
rotating immediately after releasing the left mouse button. To restrict the rotation
73
axis, select [Around X axis], [Around Y axis], or [Around Z axis] in the pull down menu of the
Rotation modes frame box.
“Random” mode
In the “Random” mode, the rotation axis is automatically set and it changes dynamically. To
stop the rotation of the objects, select either the “Drag” or “Push” mode.
11.2 Magnify
Objects are magnified in proportion to the distance of dragging the mouse upward. On the other
hand, they are shrunk in inverse proportion to the distance of dragging the mouse downward.
11.3 Translate
Drag the mouse in the Graphics Area to translate objects. When a lattice plane with a color
specified in the Lattice Planes dialog box is selected and dragged in this mode, it is interactively
moved. Otherwise the entire objects are translated with the mouse along the same direction.
11.4 Select
Several ways of selecting objects can be used. On the use of the select mode in the Vertical
Toolbar, left click on an object to select it. To select all the objects in a certain area, press the
left mouse button and drag the mouse to specify an area. On selection of new objects,
objects that have previously been selected are reset to the normal state. To select additional
objects while keeping the present objects alive, press <Shift> while clicking or dragging the
Graphics Area. Regardless of the current manipulation mode, two or more objects can be
selected by clicking or dragging on them while pressing the <Shift> key. A single object is
selectable by double-clicking on it. Objects other than atoms, bonds, and polyhedra cannot
be selected. For example, atoms behind isosurfaces can be selected.
After a single object has been selected, a variety of information about it is output in the
Text Area. The estimated standard uncertainty is also displayed for interatomic distance,
bond angle, and dihedral angle if those of lattice parameters and fractional coordinates have
been supplied. However, note that VESTA gives only rough estimates of standard
uncertainties by neglecting off-diagonal elements of the variance-covariance matrix because
no off-diagonal elements are included in crystal-data files.
Objects can also be selected by using the Objects tab in the Side Panel (see section
12.2), the Vectors dialog box (see section 11.4), and the Geometrical Parameters dialog box
(see section 14.2).
Press the <Delete> key to hide selected objects. The hidden objects are not actually deleted
but just made invisible. To restore all the hidden objects, press the <Esc> key. By hiding part of
coordination polyhedra, you can easily mix polyhedra with a ball-and-stick model.
11.4.1 Atom
On selection of an atom, site number, site name, symbol of the element, fractional coordi-
nates (x; y; z), translation vector, symmetry operations (coordinate triplet), occupancy,
isotropic atomic displacement parameter, site multiplicity plus Wyckoff symbol, and site
symmetry are displayed for the atom. For example, selection of Al and O atoms connected
with each other in a ball-and-stick model of -Al2O3 gives the following lines:
74
Atom: 1 Al Al 0.00000 1.00000 -0.15000 ( 0, 1,-1)+ -y, -x, z+1/2
Occ. = 1.000 Ueq = 0.17370 12c 3.
Atom: 2 O O 0.33333 0.97667 -0.08333 ( 0, 0,-1)+ -y+1/3, x-y+2/3, z+2/3
Occ. = 1.000 Ueq = 0.19915 18e .2
When atoms are rendered as displacement ellipsoids, information about principal axes
and mean square displacements are also output:
Atom: 2 O O 0.97667 0.64333 0.08333 ( 0, 0, 0)+ x-y+2/3, x+1/3, -z+1/3
Occ. = 1.000 Ueq = 0.19915 18e .2
Principal axes of the anisotropic atomic displacement parameters
MSD (Å^2) x y z u v w
1 0.003223 -0.043986 0.027435 -0.023150 -0.005924 0.006668 -0.001785
2 0.002471 0.027516 0.041255 -0.003390 0.010805 0.010027 -0.000261
3 0.001873 -0.013219 0.012054 0.039403 -0.001318 0.002930 0.003038
11.4.2 Bond
When a bond is selected, its length is displayed with its estimated standard uncertainty, if
any, in the Text Area as well as the Status Bar. Site number, site name, symbol of the
element, fractional coordinates (x; y; z), translation vector, and symmetry operations are also
displayed for two atoms connected with the bond. For example, on selection of an Al O bond
in the ball-and-stick model of -Al2O3, the following three lines are output:
Polyhedral volume
The volume of the selected coordination polyhedron is calculated as a fundamentally
important physical quantity [31]. For example, distortion in perovskite-type compounds, ABO3
can be quantified with polyhedral volume ratios, VA/VB, rather than with tilting angles [68].
Distortion index
D= 1 n
jli lavj ; (11.1)
n ∑i l
=1 av
where li is the distance from the central atom to the ith coordinating atom, and lav is the
average bond length.
75
Quadratic elongation
The quadratic elongation, < > [33], is defined only for tetrahedra, octahedra, cubes, dodecahe-dra, and
icosahedra:
n
1 ( li ) 2 (11.2)
∑
<>=
n i=1 l0 ;
where l0 is the center-to-vertex distance of a regular polyhedron of the same volume. < > is
dimensionless, giving a quantitative measure of polyhedral disortion which is independent of
the effective size of the polyhedron.
where m is (number of faces in the polyhedron) 3/2 (i.e., number of bond angles), ϕi is the ith
◦
bond angle, and ϕ0 is the ideal bond angle for a regular polyhedron (for example, 90 for an
◦ ′
octahedron or 109 28 for a tetrahedron).
Effective coordination number
The coordination number denotes the number of atoms coordinated to a central atom in a
coordination polyhedron. However, stating the coordination of the central atom as a single
number is somewhat difficult in relatively distorted coordination polyhedra. Several proposals
are made for calculation of a mean or ‘effective’ coordination number (ECoN) by adding all
surrounding atoms with a weighting scheme, where the atoms are not counted as full atoms
but as fractional atoms with numbers between 0 and 1. This number gets closer to zero with
an increase in distance from the central atom to a surrounding atom.
VESTA adopts ECoN [33, 34, 35, 36] defined as
∑ w i; (11.4)
ECoN = i
i exp [ 1 ]
(li/lmin)
where lmin is the smallest bond length in the coordination polyhedron.
Charge distribution
Let qX be the formal charge (oxidation number) of the central atom, X. Then, the fraction of the charge
received by an ion at a corner of a coordination polyhedron is calculated as
wq :
∆qi = i X (11.7)
ECoNX
76
The total charge, QA, received by an ion A is obtained by summing the relevant charge fractions,
∆qi’s, over its n bonds. Similarly, the charge, QX, received by the ion at the center of a
coordination polyhedron is calculated by
wi(qA/QA)i ] q :
X
Q = [∑ (11.8)
X i ECoNX
Calculation of the charge distribution [35, 36, 37] in the crystal structure depends on the current
bond specifications. In other words, the fraction of the charge is not given to or received by ions
that are not bonded to each other. Formal charges of ions are read in from data files when ICSD
files or CIFs (see 17.4.1) including oxidation states (oxidation numbers) are opened. For other file
formats, formal charges of ions have to be input in the Structure parameters tab of the Edit Data
dialog box to calculate the charge distribution. Note that charge distribution in nonstoichiometric
compounds can be calculated because occupancies are used in the calculation.
In addition to the above physical quantities, the bond valence sum, V , defined as
n
∑ ( )
l0 b li (11.9)
V = i=1 exp
is also obtainable from the bond valence parameter, l0, of the central atom [38, 39, 40]. In
VESTA, the empirical constant, b, in Eq. (11.9) is fixed at a typical value of 0.37 Å.
The bond valence model, which is a development of Pauling’s rules, has been
theoretically described in terms of classical electrostatic theory without resorting to quantum
mechanics. Nevertheless, V serves, in practice, to estimate the oxidation state of the central
atom only from bond lengths determined by X-ray or neutron diffraction.
A CIF, bvparm2011.cif, storing bond valence parameters of most chemical species is
available for download from
http://www.iucr.org/resources/data/datasets/bond-valence-parameters
This CIF is also included in the folder where the executable binary file of VESTA is located so
as to refer to it when l0 values are required. In bvparm2009.cif, an l0 value for a pair of a
3+ 2 3+
cation and an anion is obtained, e.g., l0 = 1:620 Å for Al and O and l0 = 2:172 Å for La
2
and O (b = 0:37 Å). Never select an l0 value for which b ≠ 0:37 Å.
If a coordination polyhedron is clicked in the Graphic Area while pressing the <Ctrl> key
in the manipulation mode of Select, you are asked to input a bond valence parameter, l0, for
the central metal of the selected coordination polyhedron. After entering l0, the bond valence
sum [38, 39, 40] for the central atom is calculated and displayed in the Text Area from all the
bond lengths, li, for the current polyhedron and l0.
Unless oxidation states of sites have been input in the Structure parameters tab of the Edit
Data dialog box or from some kinds of structural files, VESTA prompts you to enter an
oxidation number corresponding to V in Eq. (11.9); pressing the <Enter> key skips the
subsequent cal-culation. After entering the value, an bond length expected with Eq. (11.9) is
calculated and given in the Text Area.
77
An example of getting information on a coordination polyhedron
Suppose that a TiO6 octahedron is selected in a structural model of perovskite (CaTiO 3) [69]
4+
while pressing the <Ctrl> key. You are asked to enter the bond valence parameter of Ti (l0
= 1:815 Å) and the oxidation state of Ti (= +4) unless it has already been input in the Structure
parameters tab of the Edit Data dialog box. Then, the following data including the physical
quantities described above are output to the Text Area:
POLYHEDRON:
1 Ti1 Ti 1.00000 0.50000 0.00000 ( 1, 0, 0)+ x, y, z
----------------------------------------------------------------------------
4 O2 O 0.71030 0.71120 -0.03730 ( 1, 1, 0)+ -x, -y, -z
4 O2 O 0.78970 0.21120 -0.03730 ( 0, 0, 0)+ x+1/2, -y+1/2, -z
3 O1 O 0.92890 0.48390 0.25000 ( 1, 0, 0)+ -x+1/2, y+1/2, z
3 O1 O 1.07110 0.51610 -0.25000 ( 0, 0, 0)+ x+1/2, -y+1/2, -z
4 O2 O 1.28970 0.28880 0.03730 ( 1, 0, 0)+ x, y, z
4 O2 O 1.21030 0.78880 0.03730 ( 1, 0, 0)+ -x+1/2, y+1/2, z
----------------------------------------------------------------------------
l(Ti1-O2) = 1.9571(11) Å
l(Ti1-O2) = 1.9569(11) Å
l(Ti1-O1) = 1.9497(4) Å
l(Ti1-O1) = 1.9497(4) Å
l(Ti1-O2) = 1.9571(11) Å
l(Ti1-O2) = 1.9569(11) Å
---------------------------------------
Average bond length = 1.9546 Å
Polyhedral volume = 9.9546 Å^3
Distortion index (bond length) = 0.00168
Quadratic elongation = 1.0001
Bond angle variance = 0.3773 deg.^2
Effective coordination number = 5.9994
Charge distribution
----------------------------------------------------------------------------
delta_q: Fraction of the charge received by the ion
Q: Total charge received by the ion
q: Formal charge (oxidation number)
----------------------------------------------------------------------------
x y z delta_q Q q
4 O2 O 0.71030 0.71120 -0.03730 0.661 -1.994 -2.000
4 O2 O 0.78970 0.21120 -0.03730 0.662 -1.994 -2.000
3 O1 O 0.92890 0.48390 0.25000 0.677 -2.013 -2.000
3 O1 O 1.07110 0.51610 -0.25000 0.677 -2.013 -2.000
4 O2 O 1.28970 0.28880 0.03730 0.661 -1.994 -2.000
4 O2 O 1.21030 0.78880 0.03730 0.662 -1.994 -2.000
----------------------------------------------------------------------------
1 Ti1 Ti 1.00000 0.50000 0.00000 4.000 4.000
All the fractional coordinates, translations, equivalent positions, bond lengths with their esti-mated
standard uncertainties for the TiO6 octahedron are output before the polyhedral volume.
78
11.5 Distance
To calculate an interatomic distance, select the fifth button in the Vertical Toolbar. In this
mode, only atoms can be selected. Select a pair of atoms, A and B, to obtain the interatomic
distance between A and B. Then, the two atoms are highlighted, and a dashed line
connecting them is plotted on the screen. The A–B distance is displayed on the Status Bar
with its estimated standard uncertainty [70], if any, enclosed by a pair of parentheses. The
estimated standard uncertainty is calculated only when those of lattice parameters and
fractional coordinates have been input.
To obtain a bond length in a ball-and-stick model, selecting the relevant bond (see 11.4.2) is
faster than clicking two atoms. More information on the interatomic distance is displayed in the
Text Area, where site number, site name, symbol of the element, fractional coordinates (x; y; z),
symmetry operations, and translation vector are displayed for each of atoms A and B.
′ ′′
For example, an interatomic distance for A = C7 and B = O2 in 2 -Hydroxyl-4 -
dimethylami-nochalcone [71] is output as follows:
l(C7-O2) = 1.254(4) Angstrom
15 C7 C 0.75640 0.11620 0.34510 ( 1,-1, 0)+ -x, y+1/2, -z+1/2
4 O2 O 0.80530 0.04440 0.43470 ( 1,-1, 0)+ -x, y+1/2, -z+1/2
79
More information about the dihedral angle is output in the Text Area, where site number,
site name, symbol of the element, fractional coordinates (x; y; z), symmetry operations, and
translation vector are displayed for each of atoms A, B, C, and D. For examples, in the case
of C1 (= A), C2 (= B), C3 (= C), and C4 (= D) atoms contained in an aromatic ring of 3-[4-
(dimethylamino)phenyl]-1-(2-hydroxyphenyl)prop-2-en-1-one [71] (Fig. 11.2), the following
five lines are output in the text area:
Figure 11.2: Calculation of a dihedral angle for four carbon atoms in an aromatic ring in
3-[4-(dimethylamino)phenyl]-1-(2-hydroxyphenyl)-prop-2-en-1-one
80
omega(C1-C2-C3-C4) = 2.36(8) deg.
5 C1 C 0.58190 0.79990 0.21900 ( 0, 0, 0)+ x, y, z
7 C2 C 0.52280 0.72240 0.11990 ( 0, 0, 0)+ x, y, z
9 C3 C 0.41240 0.66560 0.09830 ( 0, 0, 0)+ x, y, z
10 C4 C 0.36010 0.68130 0.17850 ( 0, 0, 0)+ x, y, z
Data in lines No. 2 5 (selected lines in the Text Area in Fig. 11.2) can be used when imposing
nonlinear restraints on the dihedral angle in Rietveld refinement with RIETAN-FP [12].
Figure 11.3: Calculation of an interfacial angle between (100) and (511) faces in quartz.
81
Chapter 12
PROPERTIES OF OBJECTS
Properties of various objects are edited in the Properties dialog box and the Objects tab of
the Side Panel.
12.1.1 General
The first page in the Properties dialog box is General (Fig. 12.1).
Unit cell
The display of Unit cell edges in the Graphics Area is controlled with the following three radio
buttons:
• “All unit cells”: Show all the edges of unit cells within the drawing boundary.
Line styles of the unit cell edges are selected from three styles: “Solid lines”, “Dotted lines”,
and “Dashed lines”. The width of lines is input in text box {Line width}. The color of lines is
specified either (A) by entering R, G, and B values ranging from 0 to 255 or (B) by picking a
color from a color selection dialog box opened after clicking the [Select] button.
Axes
“Show Compass” is used to turn on or off a display of three arrows indicating a, b, and c axes
(or x, y, and z axes in the case of Cartesian coordinates). “Show Axis Labels” is used to turn
on or off a display of axis labels, ‘a’, ‘b’, and ‘c’ (or ‘x’, ‘y’, and ‘z’ in the case of Cartesian
coordinates).
82
Figure 12.1: General page in the Properties dialog box.
Shapes
The visibility, line width, and color of edges of crystal morphologies are input in the Shapes
frame box.
12.1.2 Atoms
The second page in the Properties dialog box is Atoms (Fig. 12.2).
Material
• {Specular} is a color of reflected light that bounces sharply in a particular direction in the
manner of a mirror. A highly specular light tends to cause a bright spot on the surface it
shines upon, which is called the specular highlight.
• {Shininess} is a property, which specifies how small and focused the specular highlight.
A value of 0 specifies an unfocused specular highlight.
Resolution
{Stacks} and {Slices} are parameters common to all the atoms, allowing you to change the
reso-lution (quality) of atoms displayed on the screen. {Stacks} denotes the numbers of
subdivisions along the Z axis (similar to lines of latitude) while {Slices} denotes the number of
subdivisions around the Z axis (similar to lines of longitude). {Stacks} and {Slices} should be
equal to each other in ball-and-stick and space-filling models where atoms are represented
by perfect spheres. In general, decreasing {Stacks} and {Slices} accelerates the rendering of
objects in the Graphics Area.
83
Figure 12.2: Atoms page in the Properties dialog box.
Atom style
Select either of the following two radio buttons for the mode of displaying atoms on the
Graphics Area: (1) “Show as balls” or (2) “Show as displacement ellipsoids”.
• “Show as balls”: Atoms are rendered as spheres. List box {Radii type} is used to select
a type of default atomic radii from the following three:
– Atomic: Metallic or covalent radii, whose values were mostly taken from Refs. [74,
75, 76].
– Ionic: Effective ionic radii compiled by Shannon [77] for representative oxidation
states and coordination numbers.
– van der Waals: van der Waals radii [78].
The user may modify default values of atomic, ionic, and van der Waals radii by editing
a text file, elements.ini, in the program folder of VESTA.
84
A B
If one or more of principal axes have negative mean square displacements, atoms are
represented by cuboids so that the unusual atomic displacement parameters can easily
be recognized (Fig. 12.4). Each cuboid is oriented in accordance with principal axes
with the dimension of the cuboid scaled according to the absolute value of the mean-
square displacement.
• “Hide non-bonding atoms”: This option hides atoms that are connected by no bonds but
the crystallographic site for those atoms have a coordination number larger than 0. In
other words, atoms are made invisible if all the coordinated atoms are lying outside of
the drawing boundary.
Select a symbol of an element from list box {Symbol} and then specify its {Radius}. The color
of an atom (element) is specified either (a) by entering R, G, and B values ranging from 0 to
255 or (b) by picking a color from a color selection dialog box opened when clicking the
[Select] button.
85
A B
Labels
Select either of the following two types for atom labels: (1) “Names of elements” or (2) “Names of
sites”. Labels are displayed near atoms with an offset along the z axis, which is specified in the
unit of Å. When option “Mark enantiomorphic sites” is checked, atoms generated by symmetry
operations of inversion, rotoinversion, or mirror are marked with a ‘,’ symbol. These positions are
enantiomorphs with respect to non-marked positions. This option can be used to draw general
position diagrams (Fig. 12.5) similar to those illustrated in “International Tables for
Crystallography,” Vol. A [30]. When an atom position coincides with an a enantiomorphic sym-
metry element, the original and enantiomorphic positions coincide. Such an atom is represented
by a ‘,’ mark on the left half of the atom separated with ‘j’ at its center.
12.1.3 Bonds
The third page in the Properties dialog box is Bonds (Fig. 12.6).
Material
Resolution
Bond style
One of the following six types of bonds are specified in this frame box:
• “Unicolor cylinder”: Each bond is drawn as a cylinder, whose color can be changed in
the Color panel.
86
Figure 12.6: Bonds page in the Properties dialog box.
• “Bicolor cylinder”: Each bond is drawn as a cylinder with colors of two atoms connected
with each other. The two colors of the bond are just the same as those of atoms
connected by bonds.
• “Color line”: Each bond is drawn as a straight line, whose color is changed in the Color
panel.
• “Gradient line”: Each bond is drawn as a line connecting two atoms with gradient distri-
bution of colors. Colors of each bond at the ends of the line are the same as those of
the atoms.
• “Dotted line”
• “Dashed line”
The radius of each cylindrical bond in the stick model is changed in text box {Radius
(cylinder)}. The actual radius of cylinders is 40 % of the input value. The line width of line-
style bonds is changed in text box, {Width (line)}. The color of bonds in the single color style
is changed in the Color tool below these two text boxes.
12.1.4 Polyhedra
The fourth page in the Properties dialog box is Polydedra (Fig. 12.7).
87
Figure 12.7: Polyhedra page in the Properties dialog box.
Material
{Specular} and {Shininess} are just the same as those in Atoms page (see 12.1.2). The
{Opacity} of coordination polyhedra is input with 255 corresponding to opaque planes (the
inside of each coordination polyhedron is invisible) and 0 corresponding to fully transparent
planes (coordina-tion polyhedra are invisible).
Polyhedral style
Styles of coordination polyhedra are selected from the following six types of expressions that
are schematically illustrated in the Polyhedral style frame box:
Planes
The surface color of each coordination polyhedron is specified here. Select an element from
list box {Central atom} to specify its color. Default colors of coordination polyhedra are the
same as those of the central atoms.
88
Edges
The visibility, line width, and color of edges of coordination polyhedra are input in the Edges
frame box.
12.1.5 Isosurfaces
The fifth page in the Properties dialog box is Isosurface (Fig. 12.8).
Material
It is the same as that in Atoms (see 12.1.2).
Isosurfaces
At the top of the Isosurfaces frame box, the minimum and maximum data values are
displayed in the unit of the raw data. Data values such as electron and nuclear densities, and
wave functions on isosurfaces are equal to {Isosurface level}, d(iso). All the points with
densities larger than d(iso) lie inside the isosurfaces whereas those with densities smaller
than d(iso) are situated outside the isosurfaces.
At the lower half of the frame box, a list of isosurface levels is shown. To add a new
isosurface level, add the [New] button, and edit the isosurface level and color. To delete an
isosurface level, select an item in the list, and press the [Delete] button. Pressing the [Clear]
button deletes all the isosurfaces.
Order of rendering polygons The sequence of drawing polygons can be changed with option
“Render from front to back.” By default, polygons are rendered from behind so that isosurfaces
89
behind the nearest surface are visible through translucent isosurfaces. In some cases with com-
plex isosurfaces, it may become difficult to understand the shape of isosurfaces because they
heavily overlap each other. When option “Render from front to back” is checked, only the nearest
surfaces are rendered. The use of this option may sometimes improve the visibility of complex
isosurfaces because neither back surfaces nor internal ones are drawn.
Figure 12.9 illustrates a difference between the two rendering modes in 3D visualization
4
of results of an electronic-state calculation for the [Cd{S4Mo3(Hnta)3}2] ion (H3nta: nitrilotri-
acetic acid) [79] by a discrete variational X method with DVSCAT [58].
A B
Figure 12.9: Comparison between the two modes of rendering isosurfaces for a molybdenum-
4
cadmium cluster [Cd{S4Mo3(Hnta)3}2] . (A) “Render from behind (default)” and (B) “Render from
front to back.” Composite images of electron-density isosurfaces colored according to
3
electrostatic potentials and a ball-and-stick model are shown with an isosurface level of 0:03a0
and electrostatic potentials on the isosurface ranging from 0:814 Ry (blue) to 0.174 Ry (red).
Kinds of isosurfaces For volumetric data having both positive and negative values, you can
select which surfaces are visible:
• {Positive},
• {Negative}.
Figure 12.10 exemplifies isosurfaces of wave functions calculated with DVSCAT [58] for a
complex ion with a ball-and-stick model superimposed on the isosurfaces. Yellow and blue
surfaces show positive and negative values, respectively.
90
4
Figure 12.10: The 64a1g orbital for the {Cd[S4Mo3(Hnta)3]2} ion
[79] with a ball-and-stick model. The isosurface levels of the
3/2 3/2
wave function were set at 0.01a0 (yellow) and 0:01a0
(blue), where a0 is the Bohr radius.
A B
Surface coloring
After volumetric data for surface coloring have been loaded, isosurfaces can be colored according to
those data. The saturation level of colors is specified as (a) a value normalized between 0 and 100,
and (b) a value corresponding to raw data. This is mostly equal to color settings for lattice planes and
sections of isosurfaces (see 12.1.6), except that the normalized values of 0 and 100 corresponds to
the minimum and maximum data values on the current isosurfaces, not the
91
data set itself.
O2
Isosurface
O1
Z
Screen
The color index, T , to determine the color of a data point with a value of d is calculated
from the unit of the raw data: the minimum saturation level, Smin, and the maximum one, Smax,
in the unit of the raw data:
S
T= d min : (12.1)
S S
max min
Data points with values larger than Smax and smaller than Smin are given the same colors as
those assigned to Smin and Smax, respectively.
12.1.6 Sections
The seventh (final) page in the Properties dialog box is Sections (Fig. 12.13).
For volumetric data, both lattice planes and sections of isosurfaces are colored according
to numerical values on them (Fig. 12.14). The saturation level of colors is specified as (a) a
value normalized between 0 and 100 corresponding the minimum and maximum data values,
respectively, and (b) a value corresponding to raw data. The color index, T , for a data point
with a value of d is calculated from the minimum saturation level, Smin, and the maximum one,
Smax, in the unit of the raw data:
T= d Smin : (12.3)
S S
max min
Data points with numerical values larger than Smax and smaller than Smin are given the same
colors as those assigned to Smin and Smax, respectively. The color of each point is
determined from T , depending on color modes. One of six color modes, that is, B-G-R, R-G-
B, C-M-Y, Y-M-C, gray scale (from black to white), and inverse gray scale, is selected from a
list box. For example, blue is assigned to 0, green to 0.5, and red to 1 in the B-G-R mode, as
shown in Fig. 12.15.
When option “Absolute values” is checked, colors are assigned on the basis of absolute
values of data.
92
Figure 12.13: Sections page in the Properties dialog box.
Option “Assign colors recursively” assigns rainbow colors recursively to data points with
values smaller or larger than the saturation levels, which affords an effect similar to contour
lines with gradient colors in between them. The opacity of sections are specified in the
{Opacity of drawn sections} box, and that of lattice planes in the Lattice Planes dialog box
(see 9.2). When {Cutoff level of lattice plane} is larger than 0, some parts of lattice planes
with data values smaller than the cutoff level are omitted.
93
D
O
D D
Data range
0% 100 %
Saturation level
Min. Max.
Figure 12.15: Relations among saturation levels, data values, and colors of sections.
94
12.2 Objects Tab in the Side Panel
12.2.2 Atoms
12.2.3 Bonds
12.2.4 Polyhedra
95
12.2.5 Slices
When a list of slices is displayed in the lower half of the
Objects page (Fig. 12.19), the first column (Slice) gives
Miller indices of slices. The second column (d (Å)) is dis-
tances from the origin to the slices, and the third column
(C) shows colors and opacities of slices.
Figure 12.19: A list of slices
in the Objects tab of the Side
Panel.
12.2.6 Shapes
When a list of forms and faces of crystal morphologies is
displayed in the lower half of the Objects page (Fig. 12.20),
the first row is assigned to a special item to control a color of
shape in the “Unicolor” style and the visibility of labels all at
once. The first column (Face) gives Miller indices of faces
grouped by crystallographically equivalent forms. The second
column (d (Å)) is distances from the origin to faces, and the
third one (C) shows colors and opacities of faces. The fourth
column (L) controls the visibility of labels for faces. If
properties of a form are edited, the modifications are applied
to all the faces that are crystallographically equivalent. Figure 12.20: A list of forms
and faces of crystal morpholo-
gies in the Objects tab of the
Side Panel.
96
Chapter 13
OVERALL APPEARANCE
The Overall Appearance dialog box (Fig. 13.1) appears on selection of the “Overall
appearance...” item under the “View” menu.
13.1 Background
The background color is specified with three values in between 0 to 255, or selected from a
color selection dialog box after clicking the button at the right of the text boxes.
97
A B
Figure 13.2: Crystal structure of milarite drawn with light effects (A) enabled and (B) disabled.
13.2 Lighting
When option “Enable lighting” is checked, the following light properties are set. When lighting
is disabled, the resulting figure looks like a traditional 2D plot with no gradients (Fig. 13.2).
• Ambient: A light that comes from all directions equally and is scattered in all directions
equally by objects. Specified by a value in between 0 to 100.
• Diffuse: A light that comes from a particular direction and hits objects with an intensity
that depends on the orientation of their surfaces. However, once the light hits a surface,
it reflects evenly off a surface and radiates in all directions. Specified by a value in
between 0 to 100.
• Lighting direction: This can be changed by dragging the “track ball” placed below the
“Lighting” frame box.
Even though the ambient and diffuse lights are evenly reflected, some part of lights are reflected
more in the manner of a mirror where most of the light bounces off in a particular direction. A light
having such a reflection component is called “specular” light. In VESTA, the intensity of the
specular light is fixed at 100 %, and a color of a reflected specular light, i.e., “specular color” is
controlled by two properties of objects, {Specular} and {Shininess}. The final color of a surface is
the sum of all three components of lights. A object surface appears to be brighter than the color of
object if the specular light is reflected in the direction of viewpoint.
Figure 13.3 shows the structure of post-perovskite, MgSiO3 [82], visualized in the parallel
and perspective modes.
f= end z : (13.1)
end start
98
A B
13.4 Depth-Cueing
Depth-cueing blends a “fog” color with the original color of each object using the blending
factor f. The factor f at depth z is computed by Both of the starting depth, start, and the
ending depth, end, are input by the user. VESTA automatically assigns the background color
of the Graphics Area to the fog color Cf . Then, the color of a rendering object, Cr, is replaced by
′ (13.2)
C =f C + (1 f) C :
r r f
When VESTA renders objects in the Graphics Area, internal coordinates of the OpenGL
scene are normalized in such a way that a radius of the bounding sphere for the scene
becomes 0.9 and that the center of the scene is placed at 0. Objects at z < start are clearly
rendered without any fog whereas objects at z > end are completely invisible.
Depth-cueing can be enabled or disabled in check box Enable depth-cueing. The effect of
depth-cueing is schematically displayed below the two text boxes Starting depth and Ending
depth.
Figure 13.4 illustrates the effect of depth-cueing on images of mordenite [83] viewed
along the c axis in the perspective-projection mode.
99
A B
Figure 13.4: Crystal structure of mordenite rendered with and without depth-cueing. (A)
Depth-cueing enabled by setting start and end parameters at 1:2 and 1:0, respectively. (B)
Depth-cueing disabled.
100
Chapter 14
UTILITIES
The Equivalent Positions dialog box (Fig. 14.1) appears on selection of the “Equivalent Posi-
tions...” item under the “Utilities” menu.
A list of general equivalent positions is displayed in this dialog box. When one of the
equivalent positions in the list is selected, the corresponding symmetry operation is displayed
in a matrix form at the upper left of the dialog box. In the right side of this dialog box,
symmetry operation W and transformation of fractional coordinates (x, y, z) with it are
explicitly described (see 6.2.6).
In the case of magnetic space groups, the time reversal term t and the resultant magnetic
moments mx, my, and mz are also displayed in the list as well as the x, y, z notation of general
positions.
101
Figure 14.2: The Equivalent Positions dialog box showing magnetic data.
This dialog box lists interatomic distances and bond angles recorded in a file *.ffe output by
ORFFE [45]. ORFFE calculates geometrical parameters from crystal data in file *.xyz created
by RIETAN-FP [12], outputting them in file *.ffe. When reading in input and/or output files of
RIETAN-FP (*.ins, *.lst), VESTA also inputs *.ffe automatically provided that *.ffe shares the
same folder with *.lst and/or *.ins. Otherwise, *.ffe can be input by clicking the [Read *.ffe]
button in the Geometrical Parameters dialog box.
VESTA allows us to locate the bonds and bond angles displayed in the Geometrical
Parameters dialog box in the Graphic Area. On selection of a bond (2 atoms) or a bond angle
(3 atoms) in this dialog box, the corresponding objects in a ball-and-stick model is selected
(highlighted), and vice versa. Thus, atom pairs and triplets associated with geometrical
parameters on which restraints are imposed in Rietveld analysis with RIETAN-FP are easily
recognized in the ball-and-stick model.
Figure 14.3 exemplifies visualization of a bond in a ball-and-stick model of fluorapatite; a
P O3 bond (grey line) selected in the dialog box is highlighted in the structural model in the
graphic window. The upper part of the dialog box displays detailed information on the P O3
bond.
Because ORFFE calculates standard uncertainties of geometrical parameters from both
di-agonal and off-diagonal terms in the variance-covariance matrix output by RIETAN-FP, the
resulting standard uncertaities are more accurate than those evaluated by VESTA from only
the diagonal terms. Accordingly, the standard uncertaities output by ORFFE should be de-
scribed in papers rather than those calculated by VESTA.
102
Figure 14.3: The Geometrical Parameters dialog box showing a list of bonds recorded for
fluora-patite [84] in FapatiteJ.ffe.
103
Standardized crystal data are generally suitable for structure refinement because atoms
in the asymmetric unit are confined in a relatively narrow region. VESTA is capable of
launching STRUCTURE TIDY [46] for this purpose.
This feature of standardization of crystal-structure data is particularly useful when
searching for compounds with similar structures using the Inorganic Crystal Structure
Database (ICSD) [87].
In both STRUCTURE TIDY and LAZY PULVERIX, the following three standard settings
are preferred to other settings:
1. monoclinic system: b-axis unique setting ( ̸= 90°),
Let n be the number of atoms in the asymmetric unit, and (xj ; yj ; zj; j = 1; 2; ; n) their
fractional coordinates. Then, the standardization parameter, , is defined as
n 1
∑
= (x 2 + y 2 + z 2 ) (14.1)
j=1 j j j
:
Note that this equation does not contain occupancies, gj. STRUCTURE TIDY selects a set of
xj , yj and zj (j = 1; 2; ; n) minimizing the value. For better distinction between different
structure-type branches, STRUCTURE TIDY further outputs another standardization
parameter, CG, depending also on lattice parameters:
2 ( ) ( ) ( )
n 2 n 2 n 2
1
a x + b j=1 y + c z
CG = nV j j j
j=1 =1
4 ∑ ∑ ∑
j
1 (14.2)
+2ab cos ( n ) + 2ac cos ( n ) 32 ;
x jy j xjzj + 2bc cos( n yj zj )
∑ ∑j ∑ 5
j=1 =1 j=1
104
conduction sheet and space group Ammm (No. 65) [88]. However, the standard setting
described in International Tables for Crystallography, volume A [30] is Cmmm; Rietveld
analysis with RIETAN-FP has to be carried out on the basis of Cmmm.
Running STRUCTURE TIDY by selecting the “Standardization of Crystal Data” item
under the “Utilities” menus in VESTA, we obtain optimum crystal data based on space group
Ammm; they are listed at the tail of data.sto:
Axes changed to : b,c,a
-------------------------------------------------------------------------
DATA YBa2Cu4O8 2.9785 0.7284
CELL 3.8708 27.2309 3.8402 *
ATOM O1 4(j) 0 0.05253 1/2 O2
ATOM Ba 4(j) 0 0.36498 1/2 Ba
ATOM Cu1 4(i) 0 0.06138 0 Cu2
ATOM O2 4(i) 0 0.14562 0 O1
ATOM Cu2 4(i) 0 0.21296 0 Cu1
ATOM O3 4(i) 0 0.28178 0 O4
ATOM O4 4(i) 0 0.44786 0 O3
ATOM Y 2(c) 1/2 0 1/2 Y
TRANS b,c,a
REMARK Transformed from setting A m m m.
STRUCTURE TIDY converted a, b, and c axes in Ammm into b, c, and a axes in Cmmm,
respectively. The value of 2.9785 is the standardization parameter ‘gamma’, the last data,
0.7284, in the ‘DATA’ line is CG, and the last data of each site is the site name input by the user.
• In the monoclinic system, the lattice parameters should be a c with the b-axis unique
setting.
• In the triclinic system, the reduction becomes more complicated owing to possible
multiple choices of basis vectors in the lattice.
Niggli’s approach to unit-cell reduction [89, 90] defines the reduced cell in terms of the shortest
vectors:
+ 2 + 3 = minimum:
1 (14.3)
STRUCTURE TIDY [46] transforms the given unit cell to a Niggli-reduced cell to describe the
structure in space group or when ‘*’ precedes a space-group symbol in an input file. If P 1 P 1
the Niggli-reduction is successful, the solution is output in the Text Area.
From the reduced form, we can obtain the metric symmetry of the lattice, which is the
highest symmetry possible for the lattice on the basis of only geometric considerations. A
105
second application of the Niggli-reduced cell is that different deformation variants of a basis
type (arstotype), probably conventionally described with different space groups and probably
non-comparable unit cells, can be compared if they have been standardized in their
corresponding Niggli-reduced cells.
The output in the Text Area on transformation to a Niggli-reduced cell for a spinel-type oxide
CoAl2O4 is given below:
Lattice type F
Space group name F d -3 m
Space group number 227
Setting number 1
Lattice parameters
Structure parameters
x y z g B Site Sym.
1 Co Co1 0.00000 0.00000 0.00000 1.000 1.000 8a -43m
====================================================================================
( 2) P -1 - i aP2
-------------------------------------------------------------------------
DATA 0.2165 0.3437
CELL 5.7240 5.7240 5.7240 60.000 60.000 60.000 *
ATOM Co 2(i) 0.12500 0.12500 0.12500 Co1
TRANS 0.5a+0.5b,0.5b+0.5c,0.5a+0.5c
REMARK Niggli-reduced cell
VESTA utilizes an external program, MADEL [91], to calculate electrostatic site potentials, ϕi,
and the Madelung energy, EM, of a crystal. Three methods are used to calculate Madelung
en-ergies: Ewald, Evjen, and Fourier methods; the Fourier method is adopted in MADEL.
MADEL was originally written by Katsuo Kato (old National Institute for Research in Inorganic
Mate-rials) and slightly modified later by one of the authors (F.I.). An advantage of using
MADEL in VESTA is that troublesome inputting of formatted data, particularly symmetry
operations, is avoidable.
The electrostatic potential, ϕi, for site i is computed by
Z
j
ϕi = ∑ ; (14.4)
j 4 ϵ0lij
where Zj is the valence (oxidation state) of the jth ion in the unit of the elementary charge, e
19 12 1
(= 1:602177 10 C), ϵ0 is the vacuum permittivity (= 8:854188 10 F m ), and lij is the
106
distance between ions i and j; the summation is carried out over all the ions j (i ≠ j) in the
crystal. In case site j is partially occupied, Zj should be multiplied by its occupancy, gj . EM
per asymmetric unit is calculated by using the formula
1∑ ϕZW; (14.5)
i i i i
EM = 2
with (occupncy) (number of equivalent positions)
Wi = : (14.6)
(number of general equivalent positions)
The summation in Eq. (14.5) is carried out over all the sites in the asymmetric unit. To obtain
the Madelung energy for the unit cell, EM must be multiplied by the number of general
equivalent positions .
Prior to the execution of MADEL, the oxidation numbers of atoms in the asymmetric unit
must be input in the Structure parameters tab of the Edit Data dialog box. Just after
launching MADEL, you are prompted to input two parameters, RADIUS and REGION:
RADIUS: Radius of an ionic sphere, s, in Å. The charge-density distribution, r, is given
by [ 2 3 4] (14.7)
1 6(r/s) + 8(r/s) 3(r/s) ;
(r) = 0
where r is the distance from the center of the ionic sphere (r < s and (r) = 0
for r s). When lines for interstitial sites are not given in the input file, set
RADIUS at a value that is large enough but less than the smallest
interatomic distance (not half of it!).
1
REGION: Reciprocal-space range (in Å ) within which Fourier coefficients are
summed up. MADEL sums up the Fourier coefficients with respect to all
hkl’s within a sphere having a radius equal to RADIUS. Choose an
1 1
appropriate value ranging from 2 Å to 4 Å according to the desired
precision of calculation. Also, check whether or not a curve for Madelung
energy versus REGION is nearly flat around the selected value of REGION.
The standard output of MADEL is displayed in the Text Area. The unit of ϕi is e/Å (1 e/Å=
14:39965 V). The accuracy of ϕi and EM obtained using MADEL is limited to 3 or 4 digits.
The following lines give part of an output file when this feature is applied to investigating
distribution of hole carriers in a high-Tc superconductor, YBa2Cu4O8 [88, 92]:
Potentials of sites in the asymmetric unit
Charge W x y z phi
Ba 2.000000 0.250000 0.500000 0.135000 0.500000 -1.330304E+00
Y 3.000000 0.125000 0.500000 0.000000 0.500000 -1.639974E+00
Cu1 2.000000 0.250000 0.000000 0.213000 0.000000 -2.164424E+00
Cu2 2.500000 0.250000 0.000000 0.061400 0.000000 -1.798502E+00
O1 -2.000000 0.250000 0.000000 0.145600 0.000000 1.310975E+00
O2 -2.000000 0.250000 0.000000 0.052500 0.500000 1.931017E+00
O3 -2.000000 0.250000 0.500000 0.052100 0.000000 1.925720E+00
O4 -2.000000 0.250000 0.500000 0.218200 0.000000 1.133873E+00
Electrostatic energy per asymmetric unit
The electrostatic energy per mole is calculated by multiplying a factor of the Avogadro constant
23 1
(6.02214 10 mol ) and EM together. In the above calculation, all the holes are assumed
107
to be doped into the Cu2 atoms on the CuO2 conduction sheet with the Cu1 atoms having an
oxidation state of +2.
Standard input files, *.pme, (see p. 141) for MADEL can be output by selecting the Export
Data... item under the File menu. This file is helpful in appending fractional coordinates of
interstitial sites to obtain their site potentials using the RIETAN-FP–VENUS integrated
assistance environment.
• Wavelength (neutron and synchrotron X-ray diffraction) or target (diffraction using char-
acteristic X rays).
• Geometry (Bragg–Brentano, transmission, or Debye–Scherrer geometry).
• Profile cutoff.
Profile parameters refined from typical diffraction data measured on a powder diffractometer,
which is often used on a routine basis, are preferable to those in one of template files
included in the distribution files of the RIETAN-FP–VENUS system. Further, a variable, NPAT,
must be given in the template file to specify a graph plotting program:
1. NPAT = 1: a pair of files, *.plt including commands to plot a graph and *.gpd storing data to
1
be graphed, for gnuplot (default),
2 3
2. NPAT = 2: an Igor text file, *.itx, for Igor Pro or WinPLOTR,
4
3. NPAT = 3: a text file, *.itx, for RietPlot.
The use of gnuplot is highly recommended because it is free software supporting many
features; commands may be modified and added in *.plt to change the appearance of a
graph. For details in the content of *.plt, refer to RIETAN-FP_manual.pdf. If NPAT = 1, no
graph plotting program need to be specified in the Preferences dialog box.
Stable versions of gnuplot are included in the distribution files of the RIETAN-FP–VENUS
5
system for Windows and macOS.
Gnuplot can output a PDF file, *.pdf, whose margins are automatically cut, which is very
convenient when exporting the resulting PDF file to an application such as T EX, Microsoft
PowerPoint and Word, and Adobe Illustrator. (Windows) or ps2pdf (macOS).
1
http://www.gnuplot.info/
2
http://www.wavemetrics.com/
3
http://www-llb.cea.fr/winplotr/winplotr.htm
4
http://www.nims.go.jp/denzai/publ_rietplot.html
5 http://fujioizumi.verse.jp/download/download_Eng.html
108
Figure 14.4: Simulation of an X-ray powder diffraction pattern of a high-Tc superconductor,
YBa2Cu4O8 [88], which was plotted with Igor Pro.
On selection of the “Powder Diffraction Pattern...” item under the “Utilities” menu, a series
of procedures, i.e., (a) generation of an input file, *.ins, for RIETAN-FP, (b) execution of
RIETAN-FP in the simulation mode (NMODE = 1), and (c) representation of the resulting text
file(s) by a graphing program are executed by VESTA as if the two programs were
implemented in VESTA (Fig. 14.4). Text files, *.ins, *.lst (standard output), and a pair of files
(*.plt and *.gpd for gnuplot) or *.itx, output by RIETAN-FP are saved in folder tmp under the
settings directory (see 17.2).
For those who have not any program to plot data in *.itx with the Igor text format, a small
6
Java program named PowderPlot, which was developed by one of the authors (K.M.), is
included in the archive file of VESTA.
VESTA automatically standardizes lattice settings when outputting *.ins to satisfy con-
ditions required by LAZY PULVERIX [86], which is embedded in RIETAN-FP to generate
6 http://jp-minerals.org/powderplot/
109
diffraction indices, hkl, and their multiplicities, m. Standard lattice settings adopted in LAZY
PULVERIX are just the same as those in STRUCTURE TIDY (see 14.3) [85, 46]: second
settings in centrosymmetric space groups with more than one origin choice and first settings
in the other space groups. Then, beware that diffraction indices displayed in PowderPlot are
incompatible with the coordinate system adopted in VESTA if the lattice setting in VESTA is
other than a standard one.
A list of reflections output at the tail of *.lst if NPRINT > 0 provides us with a variety of
information about hkl reflections included in a 2 range specified by the user. For instance,
simulation of an X-ray diffraction pattern for fluorapatite gives the following output where line
heads were deleted for convenience:
h k l Code 2-theta d Ical |F(crys)| |F(magn)| POF FWHM m Dd/d
1 0 0 +1 10.895 8.11382 12287 16.7800 - 1.001 0.07451 6 0.013637
1 0 0 +2 10.923 8.11382 6113 16.7800 - 1.001 0.07451 6 0.013602
1 0 1 +1 16.877 5.24919 4465 11.1545 - 0.999 0.07443 12 0.008757
1 0 1 +2 16.919 5.24919 2221 11.1545 - 0.999 0.07443 12 0.008734
1 1 0 +1 18.929 4.68452 2079 12.0931 - 1.001 0.07445 6 0.007795
1 1 0 +2 18.976 4.68452 1034 12.0931 - 1.001 0.07445 6 0.007775
2 0 0 +1 21.891 4.05691 9679 30.2953 - 1.001 0.07453 6 0.006726
2 0 0 +2 21.946 4.05691 4814 30.2953 - 1.001 0.07453 6 0.006709
1 1 1 +1 22.945 3.87283 8580 21.1855 - 1.000 0.07457 12 0.006413
1 1 1 +2 23.003 3.87283 4268 21.1855 - 1.000 0.07458 12 0.006397
.....
Selection of “Utilities” menu ) “Structure Factors...” pops up the Equivalent Positions dialog
box (Fig. 14.5), whose use makes it possible to calculate structure factors from structural models
(structure parameters) or volumetric data (electron or coherent-scattering length densities, ).
110
If multiple phase data are included in the data set, select a phase, for which structure
factors are calculated, at the top of the dialog box. If the phase data contain both structural
and volumetric data, select {Structural model} or {Volumetric data} in the list box to specify
the type of data from which structure factors are calculated. In the Source list box, {X-ray}
and {Neutron} beam sources can be switched, and some typical characteristic X-rays can be
selected at a list box next to the Source list box. For the neutron source and the X-ray source
with its target specified as {Custom}, the wavelength can be edited in the text box. The
minimum lattice-plane spacing, d, of reflections whose structure factors are calculated is also
specified in the text box.
In X-ray diffraction, structure factor F of reflection h with diffraction indices hkl is formu-lated
as
∑j
n
[ ( )]
F (h) = (14.8)
gj fj(h)Tj (h) exp 2 i hxj + kyj + lzj ;
=1
where j is the atom number, n is the total number of atoms in the unit cell, gj is the
occupancy, fj(h) is the (X-ray) atomic form factor, Tj(h) is the Debye–Waller factor, and xj, yj ,
and zj are the fractional coordinates. The atomic form factor f(h) is a complex number that
can be divided into four separate components as
′
f(h) = f0(h) + f + if + fNT;
′′ (14.9)
′ ′′
where f0(h) is the coherent scattering factor, f and f are the real and imaginary components
of the dispersion correction, and fNT is the nuclear Thomson scattering [93],. The coherent
atomic scattering factor f0, which is a function of sin / , is usually analytically approximated by
a sum of exponentials plus a constant:
N 2
∑ sin
f0(sin / ) = a ( −bi 2 ) (14.10)
i=1 i exp + c:
VESTA utilizes a set of parameters (ai, bi, and c) with N = 5 tabulated in [94]. These
′ ′′
parameters can be used in the range of 0 < sin / < 6. The real and imaginary parts, f and f ,
of the dispersion correction are calculated by linear interpolation of theoretically calculated
form factors in energy ranges of E = 1–10 eV and E = 0:4–1:0 MeV [95, 96]. The nuclear
Thompson scattering term fNT is small and negative in the phase relative to the electronic
form factor, with small (negligible) dependences on and sin / .
In the case of neutron diffraction, fj(h) has to be replaced by the coherent-scattering
length [10], bc:
∑
n
[ ( )]
(14.11)
j
111
14.182259 5.936858 0.112726 34.958481 0.390240 5.803000
X-ray dispersion coefficients for λ = 0.154059 nm
If all the parameters of some elements are 0, the element names are not recognized by the
program; consequently, such elements do not contribute to calculated structure factors at all.
The element names can be edited in the Structure parameters page of the Edit Data dialog
′ ′′
box. When an option of “Ignore the effect of X-ray (anomalous) dispersion” is enabled, f , f ,
and fNT are set at 0, and none of these values are output in the Text Area.
Structure factors, F (h), can be calculated by the Fourier transform of (x; y; z):
∫∫∫
F (h) = (x; y; z) exp[2 i(hx + ky + lz)]dxdydz: (14.12)
Only unique reflections should be given in the file because all the equivalent reflections
are automatically generated according to the current space group. If parts of reflections in the
list do not satisfy conditions limiting possible reflections in the space group, a warning
message will appear, and their structure factors will be reset to 0. F (000) is separately input
in the text box placed above the list of structure factors. In addition, the spacial resolution in
real space must be specified in the unit of Å. After all the data have been input, press the
[Calculate] button to carry out Fourier synthesis from them.
112
Figure 14.6: Fourier Synthesis dialog box.
electron densities is specified in the unit of Å. The number of grids along each
crystallographic axis, Nx, is automatically set such that Nx gives a resolution close to the
specified value and that Nx satisfies symmetrical constraints. Then, the electron density, (x; y;
z), at the coordinates of (x; y; z) is calculated by
N
Nx y Nz
∑2 k ∑2 ∑2
2 2 2
1 F (h) exp[ 2 i(hx + ky + lz)]: (14.14)
(x; y; z) =
V
N N N
h= x = y l= z
If the resolution is not sufficiently high, spurious noises appear in the resulting densities because
of the so-called termination effect in Fourier synthesis. On the other hand, atomic scattering
1
factors compiled in Ref. [94] are reliable only in a range of sin / < 6 Å , which corresponds to a
spacial resolution of approximately 0.042 Å. VESTA simply extrapolates the atomic scattering
factors given in Ref. [94] even when the specified resolution is smaller than 0.042 Å.
1. “From Model Electron Densities” submenu: from structure parameters and atomic
scattering factors of free atoms [94].
113
2. “From Model Nuclear Densities” submenu: from structure parameters and coherent-
scattering lengths of atoms.
3. “From Volumetric Data” submenu: by convolution of volumetric data that are currently
displayed.
On selection of the first two submenu, structure factors are calculated with Eqs. 14.8 and
14.11, respectively, with the spacial resolution of electron or coherent-scattering-length
densities spec-ified in the unit of Å. Then, the Patterson density, P (x; y; z), at the coordinates
of (x; y; z) is calculated by
Nx Ny Nz
∑ ∑
k ∑2
2
2
2 2
2
1 2
jF (h)j exp[ 2 i(hx + ky + lz)]; (14.15)
P (x; y; z) =
V
N N N
h= x = y l= z
where Nx, Ny, and Nz are the numbers of grids along the x, y, and z axes, respectively (see
section 14.9).
When Patterson functions are calculated from the current volumetric data, structure factors are calculated
by
∑∑∑ (14.16)
F (h; k; l) = (x; y; z) exp[2 i(hx + ky + lz)];
x y z
and then P (x; y; z) is computed by Eq. (14.15). Therefore, the spacial resolution of the
Patterson function is the same as that of the volumetric data.
114
8.095806E-002 2.272991E+002
1.214371E-001 1.836505E+002
1.619161E-001 1.328268E+002
.....
The first line gives a title, followed by two lines for the fractional coordinates of both ends,
those of Ti and O atoms in this case. The fourth line is the total number of data points (= 50),
followed by 50 lines. In each of the 50 lines, the first and second columns give a distance
from “Point 1” in the unit of Å and an interpolated value, respectively.
The integrated values represent numbers of electrons per each site, which is calculated by
sum-mation of numbers of electrons over the Voronoi volume of the site.
m [ ] ] 2 (14.17)
(r) (r)
115
2
where ℏ is Dirac’s constant, and m is the electron mass; the gradient ℏ-expansion of the
Green function around the classical Thomas-Fermi approximation [98] yields = 1/72 and k =
1/6. The electronic-energy density equals the sum of g(r) and (r):
116
Chapter 15
TWO-DIMENSIONAL DATA
DISPLAY
• Menu bar: On Windows and Linux, the menu bar is placed at the top of the Main
Window. On macOS, menus are displayed at the top of the screen. “File” and “View”
menus are available in this window.
117
• Side Panel: Colors, contours, and other properties of a 2D image are controlled in the
Side Panel. The minimum and maximum values in the 2D slice/projection data are
displayed at the top of the panel.
15.2 Menus
• File menu
– Export Raster Image...: Export a graphic image as a file with a raster (pixel-
based) format.
– Close: Close the 2D Data Display window.
• View menu
15.3.1 Rotation
These six buttons are used to rotate objects around one of the x, y, and z axes. The step
width of rotation (in degrees) is specified in a text box next to the sixth button:
118
15.3.2 Translation
Translate upward
Translate downward
Translate leftward
Translate rightward
These four buttons are used to translate objects toward up, down, left, and right, respectively.
The step width of translation (in pixels) is specified in a text box next to the fourth button:
15.3.3 Scaling
Zoom in
Zoom out
These three buttons is used to change object sizes. The step width of zooming (in %) is
specified in a text box next to the third button:
Rotate
Translate
Magnify
119
Figure 15.2: A dialog box for three different modes to draw a 2D image
Vectors parallel to the (hkl) plane: A pair of vectors, V1 and V2, is specified in the expres-
sion of the lattice vector, ua + vb + wc, on the (hkl) plane. The two lattice vectors should
be parallel to the (hkl) plane but should not be parallel to each other.
Center of projection: Projection of this point along the [hkl] direction intersects the slice at
the center of the slice.
In this mode, two-dimensional data on the slice are recalculated by linear interpolation.
When the numbers of data points along x, y, and z directions in the volumetric data are N[0],
N[1], and N[2], respectively, the number of data points, Nx, along the [u1v1w1] direction is
computed by
Nx = (int)(N[0] u1 + N[1] v1 + N[2] w1): (15.1)
The same procedure is applied to the number of data points along the [u2v2w2] direction.
120
A
B
[u2 v2 w2]
[u1 v1 w1]
Figure 15.3: Electron-density distribution drawn with the three different modes for the (111)
plane in rutile-type TiO2. (A) “(hkl) plane in the bounding box” mode, (B) “(hkl) plane defined
by two vectors” mode, and (C) “Project along [hkl] axis” mode. The electron densities were
determined by MEM from X-ray powder diffraction data.
121
Distance from the center of projection
Center of projection
Projection vector [ ]
Figure 15.4: Schematic image of the relation between the center of the projection
and the center of the slice.
h, k, and l: Three values in an equation to represent the lattice vector, R = ha + kb + lc, where
a, b, and c denote fundamental lattice vectors. The [hkl] direction is defined with R.
Range of projection: A pair of values specified by distances from the Center of projection.
They are input in the unit of either its lattice-plane spacing, d, or Å.
Vectors parallel to the (hkl) plane: A pair of vectors, V1 and V2, is specified in the expres-
sion of the lattice vector, ua + vb + wc, on the (hkl) plane. The two lattice vectors should
be parallel to the (hkl) plane but should not be parallel to each other.
Center of projection: Projection of this point along the [hkl] direction intersects the slice at
the center of the slice.
122
Figure 15.5: Bird’s eye view of a (001) slice of rutile-type TiO2.
A B
Figure 15.6: Electron-density distribution on the (001) plane in rutile-type TiO 2. (A) without
grid edges and (B) with grid edges drawn with solid lines.
In the Contours page (Fig. 15.7), various properties of contour lines are specified. Check
“Draw contour lines” to draw contour lines, which are plotted in two different modes: linear or
logarithmic.
Contours are plotted in linear and logarithmic modes as solid lines (style L1) and dashed
ones (style L2). The width and colors of these lines are specified in the text boxes and
buttons in the Style frame box.
In the linear mode, lines are drawn at every Interval in data ranging from Min. to Max. The
numerical value at the Nth line, F (N), is given by
In this mode, line style L1 (solid line) is applied to lines of positive values, and line style L2
(dashed line) to lines of negative values.
123
A B
Figure 15.7: Contours page in the Side Panel. (A) Linear mode and (B) logarithmic mode.
F (N) = A B
N/Step
; (15.3)
where A and B are constants specified by user. Line style L1 is applied to every “step” lines
starting from the first one, and for other lines the L2 style is applied. For example, when Nmin
is set to an integer value, the L1 style is used for lines with integer values of N/Step, and the
L2 style is used for lines with non-integer values of N/Step.
124
Chapter 16
PREFERENCES
Default settings for the behavior of VESTA are changed in the Preferences dialog box (Fig.
16.1).
125
parameters are copied from the template file. Beware that the format of the template file must
be compatible with the version of RIETAN-FP used by VESTA. The use of *.ins for an older
version of RIETAN-FP may cause an error on execution of its latest version. The default
template file distributed with VESTA is that for the latest version of RIETAN-FP.
16.4 Animation
Two settings for animation speed, i.e., the rotation angle per frame and the interval between
frames are specified in frame box Animation. The rotation angle per frame is specified in text
box {Step (deg./frame)}. The interval between frames is specified in milliseconds in text box
{Interval (ms)}. The number of frames per second is
1/ (t +t ); (16.1)
draw interval
where tdraw is the time required to draw a single frame of an image on the Graphics Area,
and tinterval is the interval between frames.
Rg = Rd s f; (16.2)
where Rd is the number of stacks or slices for display on the screen, s is the scale of the image,
and f is the factor to increase Rg specified in text box {Increasing factor for stacks/slices}.
126
16.7 Default isosurface level
The default isosurface level, d(iso), after opening a file storing volumetric data, , is calculated
as
d(iso) = j j + n (j j); (16.3)
where j j is the average of j j’s, (j j) is the standard deviation of j j, and n is a parameter to
adjust d(iso).
127
Chapter 17
A B
128
17.2 Directories for User Settings
VESTA loads and saves two files, VESTA.ini and style/default.ini in a directory for user
settings to store user settings for the program and for graphics, respectively. A directory, tmp,
to store temporary files is also created in the same directory. The location of the settings
directory depends on operating systems, and is determined in the following order of priority.
17.2.1 Windows
• When environmental variable VESTA_PREF is defined, the specified directory is used.
• When the user has permission to write in the program directory of VESTA, that
directory is used so as not to put the home directory in disorder.
• A directory %HOMEPATH%\AppData\Roaming\VESTA\ is created on the first execution of
VESTA to use it.
17.2.2 macOS
• When environmental variable VESTA_PREF is defined, the specified directory is used. To
define environmental variables for GUI applications, they must be described in a text
file named environment.plist under a hidden directory, ~/.MacOSX. To make it easier to
define VESTA_PREF, a utility called set_VESTA_PREF.app is included in the RIETAN-
1
VENUS package ; refer to Readme_mac.pdf for details in set_VESTA_PREF.app.
• A directory ~/Library/Application Support/VESTA/ is created on the first execution
of VESTA, and it is used. Note that ~/Library is a hidden directory under the home
directory (~/).
17.2.3 Linux
• When environmental variable VESTA_PREF is defined, the specified directory is used.
• A hidden directory ~/.VESTA/ is created on the first execution of VESTA to use it.
1 http://fujioizumi.verse.jp/download/download_Eng.html
129
17.3 Files Used by VESTA
130
3. American Mineralogist Crystal Structure Database [99] (*.amc)
http://rruff.geo.arizona.edu/AMS/amcsd.php
In part of AMCSD text files, extra characters are attached to space group names, and
non-standard space-group symbols are used. Some of such non-standard space-group
settings may not be read in correctly. In such a case, modify the space-group symbol
appropriately and then change the setting number in the Unit cell tab of the Edit Data
dialog box if necessary. Sometimes, you have to convert fractional coordinates by
yourself if a non-standard setting that is not described in International Tables for
Crystallography, volume A [30] is adopted.
4. asse (*.asse)
http://www.nims.go.jp/cmsc/staff/arai/asse/
5. Chem3D (*.cc1)
http://openbabel.org/docs/2.3.1/FileFormats/3D_viewer_Formats.html
6. Crystallographic Information File (CIF; *.cif) [100]
http://www.iucr.org/resources/cif/
CIF has a variety of formats for crystal data. You can get detailed information about CIF
from the above Web site. For example, CIF files may contain Cartesian coordinates, but
VESTA cannot input them. Note that VESTA does not support all the formats allowed in
the CIF format. For example, Cartesian coordinates included in CIFs cannot be input
with VESTA; only fractional coordinates should be given in CIFs. For readable formats,
refer to example files of *.cif in VENUS/examples/VICS/CIF.
In the case of *.cif containing multiphase data, all the data are input in the same tab and
overlapped with each other. To visualize only one phase in such a case, select “Edit”
menu ) “Edit Data” ) “Phase…”. Select an unnecessary phase in the list, and press
[Delete] button.
The CIF specification presents multiple ways of inputting entries for space-group
symme-try. VESTA searches the entries in the following order:
_symmetry_equiv_pos_as_xyz
_symmetry_space_group_number
_symmetry_space_group_name_H-M
If one of the above entries are given in *.cif, VESTA recognizes space-group symmetry.
For example, CIFs created by SHELX-97 [97] through WinGX [101] contain neither
space-group number nor space-group symbol but have a list of symmetry operations.
If none of the above entries are correctly given, the space group is regarded as P 1. In this
case, set the space group in the Unit cell tab of the Edit Data dialog box after opening such
a file, or modify a line relevant to the space group in the CIF in the following manner:
_symmetry_space_group_number 12
or
_symmetry_space_group_name_H-M 'P 21/n'
131
8. Crystal Structure Search and Retrieval (CSSR; *.cssr)
http://www.maciejharanczyk.info/Zeopp/input.html
In a file with the CSSR format, a setting numbers may be given after ‘OPT =’. Unfortu-
nately, we have no information about setting numbers in this format. Then, please
change it in the Unit cell tab of the Edit Data dialog box if necessary.
9. Cambridge Structural Database [102] (CSD/FDAT; *.csd, *.fdt))
http://www.ccdc.cam.ac.uk/Solutions/CSDSystem/Pages/CSD.aspx
10. DL_POLY [103] format (CONFIG, REVCON, *.config)
https://www.scd.stfc.ac.uk/Pages/DL_POLY.aspx
11. FEFF input file (feff.inp)
http://www.feffproject.org/
FEFF [104, 105] is an automated program for ab initio multiple-scattering calculations of
X-ray Absorption Fine Structure (XAFS) and X-ray Absorption Near-Edge Structure
(XANES) spectra for clusters of atoms. The names of input files for FEFF must be
either feff.inp or FEFF.INP.
12. FHI-aims input file (*.in)
https://aimsclub.fhi-berlin.mpg.de/
FHI-aims is an accurate all-electron, full-potential electronic structure code package for
computational materials science.
13. Output file of the Elk FP-LAPW Code (GEOMETRY.OUT)
http://elk.sourceforge.net/
Elk is an all-electron full-potential linearised augmented-planewave (FP-LAPW) code
for determining the properties of crystalline solids.
14. GSAS [18] format (*.EXP)
http://www.ncnr.nist.gov/xtal/software/gsas.html
15. Inorganic Crystal Structure Database [87] (ICSD; *.ics)
http://www2.fiz-karlsruhe.de/icsd_home.html
Two retrieval programs, RETRIEVE for MS-DOS and FindIt for Windows, of ICSD
output text files with quite different formats. VESTA is capable of reading in both types
of the crystal data files.
In these *.ics files, extra characters are sometimes attached to space group names,
e.g., ’P 42/n m c S’, which should be ’P 42/n m c’ (P42/nmc). In addition, full Hermann-
Mauguin space-group symbols are sometimes given in ICSD text files. In such a case,
an error message appears in both the Text Area and a message box. Read it carefully
to proceed to a next operation. When you encounter this type of an error, it is strongly
recommended to output a CIF instead of *.ics.
16. ICSD-CRYSTIN: (*.cry)
132
In part of MINCRYST text files, extra characters are attached to space group names,
and non-standard space-group symbols are used. In such a case, an error message
appears in the Text Area. Such space group names need to be appropriately changed.
Change the setting number in the Unit cell tab of the Edit Data dialog box if necessary.
19. MOLDA [106] (*.mld) http://www3.u-
toyama.ac.jp/kihara/cc/mld/readme.html
The Web site of MOLDA has been closed because the author, Hiroshi Yoshida passed
away in 2005. Then, the MODRAST/MOLDA format for *.mld is briefly explaned here.
This format consists of the following lines:
(a) Line No. 1: A comment about the compound, e.g., its name
(b) Line No. 2: Number of atoms, na, in the compound
(c) Lines No. 3 (3 + na): Cartesian coordinates (x, y, and z) and an atomic number
(d) Line No. (4 + na): Number of bonds, nb, in the compound
(e) Lines No. (5 + na) (5 + na + nb): A pair of atom numbers
For example, in the case of ethylene, where na = 6 and nb = 5, the following lines are
required:
"Ethylene (CH2=CH2)"
6
.66958, 0, 0, 6
1.23873, -.94397, 0, 1
1.23873, .94397, 0, 1
-.66958, 0, 0, 6
-1.23873, -.94397, 0, 1
-1.23873, .94397, 0, 1
5
1, 4
1, 2
1, 3
4, 5
4, 6
133
of space group Fd m. The input of ’I-230-2’ causes an error. Lattice parameters must be 3
given within one line in the following way:
XMol developed at the Minnesota Supercomputer Center is a utility for creating and
viewing graphic images of molecules.
28. F01 for SCAT [58, 111] and C04D for
contrd http://www.dvxa.org/
VESTA need to read in c04d, which is an input file of contrd, in addition to f01 if a
structural model is to be overlap with volumetric data. For this purpose, the dimensions
of a boundary box (an area where volumetric data are output to text files) in c04d are
required. For details in contrd, refer to Readme_contrd.txt in the VENUS package [15,
16]. Of course, c04d and f01 should be placed in the same folder. If c04d has not been
input by VESTA, atomic coordinates are treated in the Cartesian coordinate system, as
is the case for *.xyz files. In this case, volumetric data, *.scat and *.sca, cannot be
overlapped with the structural model.
All the atoms recorded in f01 must be included within the above boundary box.
Otherwise, atomic coordinates are normalized within the boundary box by assuming
periodicity, which leads to the appearance of incorrect structures in the Graphics Area.
2
The assistance environment for the DV-X method using Hidemaru Editor is very con-
venient when carrying out a series of electronic-state calculations.
The Web site of Genta Sakane of Okayama University of
Science: http://www.chem.ous.ac.jp/~gsakane/
is very useful for those who like to visualize physical quantities calculated with contrd.
3
“Introduction to the assistance environment for the DV-X method ” distributed there is a
detailed Japanese document suitable for beginners in the DV-X method as well as
VESTA.
2 http://fujioizumi.verse.jp/visualization/VENUS.html#assistance_environment
3
http://www.chem.ous.ac.jp/~gsakane/dvxa_assistance_environment_1.html
134
29. MXDORTO/MXDTRICL [112, 113] (FILE06.DAT, FILE07.DAT)
http://kats-labo.jimdo.com/mxdorto-mxdtricl/
MXDORTO and MXDTRICL are Fortran programs for molecular dynamics simulation.
135
followed by elements of a three-dimensional array, D
with any number of data in each line and at least one space between two real data (free
format). Note that voxels at N1+1, N2+1, and N3+1 lie on x = 1, y = 1, and z = 1,
respectively. A top part of an example of *.ted is given below:
TiO2 (rutile)
4.59393 4.59393 2.95886 90.00000 90.00000 90.00000
65 65 65
-3.700291E+004 -3.355551E+004 -2.515189E+004 -1.581991E+004 -8.551052E+003 -4.096268E+003
-1.802126E+003 -7.564262E+002 -3.147750E+002 -1.347187E+002 -6.129182E+001 -3.045389E+001
-1.683229E+001 -1.043066E+001 -7.226008E+000 -5.529954E+000 -4.588889E+000 -4.036150E+000
-3.669965E+000 -3.365605E+000 -3.043458E+000 -2.664371E+000 -2.229370E+000 -1.771544E+000
-1.337361E+000 -9.658093E-001 -6.762105E-001 -4.680259E-001 -3.280885E-001 -2.393913E-001
.....
VESTA allows you to save a file, *.?ed, with the general volumetric-data format.
Accord-ingly, VESTA serves as a converter from a binary file to a text one.
35. Periodic volumetric-data (text format) (*.grd)
http://www.ncnr.nist.gov/xtal/software/gsas.html
Text files for 3D Fourier maps output by GSAS [18]. A program, 3DBVSMAPPER [118],
for automatically generating bond-valence sum images is also capable of outputting
files with this format. To create *.grd files that can be input by VESTA, select option “C
– Select section (X, Y, or Z) selection” in Fourier calculation setup and then enter “X” for
prompt “Enter section desired (X,Y,Z - choose Z for DSN6 maps).”
The file format is essentially the same as the general volumetric-data format, except
that the three-dimensional data array, D, are output in the following range:
136
as a part of package [15, 16]. Three-dimensional numerical data recorded in *.SCA or
*.SCAT are drawn without any conversion.
Units for electron densities, electrostatic potentials, and wave functions are,
3 3/2
respectively, bohr , Ry (rydberg), and bohr , where bohr is the atomic unit for length,
11
i.e., 1 bohr = a0 = 5.29177211 10 m = 0.529177211 Å (a0: Bohr radius), and 1 Ry =
18
Eh/2 = 2.179 871 9 10 J (Eh: hartree).
Refer to No. 23 in 17.4.1 for the assistance environment for the DV-X method and its
detailed document written in Japanese.
40. WIEN2k (*.rho) http://www.wien2k.at/ (WIEN2k)
http://www.nims.go.jp/cmsc/staff/arai/wien/venus.html (wien2venus.py)
5
A script, wien2venus.py, coded in Python by Masao Arai of NIMS makes it possible to
export electron densities calculated with WIEN2k [41] to a text file, *.rho, which is in turn
3
visualized by VESTA. The unit of electron densities stored in this file is bohr .
41. WinGX 3D Fourier map (*.fou)
http://www.chem.gla.ac.uk/~louis/software/wingx/
Text files of 3D Fourier Maps output by WinGX [101]. To create *.fou that can be input
with VESTA, open “FOURIER MAP Control Panel” of WinGX from Maps menu !
FOURIER MAP ! Slant plane. Select “3D Fouier (Beevers-Lipson)” and “Write
MarchingCubes File” options and set “Projection” at Z axis. The minimum and
maximum of “Summation limits” should be set at 0 and 1 for all of X, Y, and Z axes. A
resolution for each axis should be carefully set in view of the following matter.
In this format, data points along each axis are not uniformly distributed but just placed
at intervals of given resolutions. When length of an axis is L and resolution is set at d,
the number of data points, NPIX, is set at the integer part of L/d + 1. The output files
have a general grid format with nearly correct periodicity only if L/d is close to an
integer. It is recommended that special positions lie exactly on data grid. For example, if
mirror planes lie at x = 1/4 and x = 3/4, L/d should be multiples of 4.
42. X-PLOR/CNX (*.xplor)
http://en.wikipedia.org/wiki/X-PLOR (CNX)
http://superflip.fzu.cz/ (Superflip)
Superflip [119] is a computer program for ab initio structure solution of crystal structures
by the charge-flipping method [120, 121, 122]. Superflip outputs electron densities in
the unit cell in file *.x-plor with the X-PLOR format [123], which can be directly
visualized by VESTA.
137
44. GAMESS [126] input and volumetric data files output by MacMolPlt [127]
http://www.msg.ameslab.gov/GAMESS/GAMESS.html (GAMESS)
https://brettbode.github.io/wxmacmolplt/ (MacMolPlt)
A GAMESS input file, *.inp, can readily be obtained from a GAMESS log file, *.log, with
MacMolPlt. At first, examine the unit of final Cartesian coordinates of atoms after
convergence in this file using a text editor. Then, run MacMolPlt to open *.log. In the
Windows menu, select ‘Coordinates’ through the name of *.log and check whether the
unit of Cartesian coordinates is Å or bohr (au). The unit should be Å, as is usual with
VESTA. If the unit is Bohr, select ‘Convert to Angstroms’ under the Molecule menu.
Then, select ‘Input Builder’ through the name of *.log in the Windows menu and click
[Write File] to create *.inp storing atomic symbols and Cartesian coordinates.
Next, a volumetric data file, *.mmp, in which the origin of the 3D grid is recorded, has to
be output. In the Windows menu, select ‘Surfaces’ through the name of *.log. Specify
an item from ‘3D Orbital’, ‘3D Total Electron Density’, and ‘3D Molecular Electrostatic
Potential’. In the dialog box that appears subsequently, change the number of grid
points and grid size appropriately, select an orbital (in the case of ‘3D Orbital’), and click
[Update]. An isosurface plus a ball-and-stick model appear in the *.log window. The
number of grid points, origin, and grid increment are displayed by clicking
[Parameters...]. Then, click [Export...]. Specify the name and location of *.mmp with the
same basename as *.inp. Note that *.inp and *.mmp must share the same folder.
45. Gaussian Cube format (*.cube, *.cub)
http://www.gaussian.com/
Text files storing electron densities, spin densities, electrostatic potentials, wave
6
functions, and so forth calculated with Gaussian [128, 129] with a keyword of ‘Cube’. .
7
Cube files can also be created by Firefly (previously known as PC GAMESS).
46. VASP (*.vasp, CHG, CHGCAR, PARCHG, LOCPOT, ELFCAR, POSCAR, CONTCAR)
http://www.vasp.at/
http://www.materialsdesign.com/medea (commercial software MedeA including VASP
as its component)
All the above files are text files storing crystal-structure and volumetric data output by
VASP [130, 131].
CHG stores lattice vectors, atomic coordinates, and total charge densities multiplied by
the unit-cell volume, V . PAW one-center occupancies are added to them in CHGCAR.
Though both CHG and CHGCAR provide us with the same information on the valence
charges, the file size of CHG is smaller than that of CHGCAR owing to the lower
accuracy of the numerical data. PARCHG has the same format as CHG has, storing
partial charge densities of a particular k-point and/or band. When these files are read in
3
to visualize isosurfaces and sections, data values are divided by V in the unit of bohr .
3
The unit of charge densities input by VESTA is, therefore, bohr .
LOCPOT contains lattice vectors, atomic coordinates, and Coulomb potentials (unit: eV),
i.e., total potentials without exchange-correlation contributions (unless the line of LEXCHG=
1 is commented out in main.F). ELFCAR, which has the same format as CHG has, stores
dimensionless electron localization functions (ELF). POSCAR and CON-TCAR include
lattice vectors, atomic coordinates, and optionally starting velocities and
6 Units used in Gaussian 03 are described in http://www.gaussian.com/g_tech/g_ur/k_constants.htm
7
http://classic.chem.msu.su/gran/gamess/
138
predictor-corrector coordinates for molecular-dynamics calculation. POSCAR and CON-
TCAR, respectively, correspond to the initial structure and the final one output by VASP
at the end of a job; CONTCAR can be used to continue the job.
Owing to the absence of symbols of elements or atomic numbers in these files, they
must accompany another file, OUTCAR, to display structural models. OUTCAR may be
re-named *.out having the same basename as *.vasp has. Only lines following
‘POTCAR:’ at the top part of OUTCAR are required to obtain symbols of elements.
When files are read in for surface coloring, no data values are scaled unless the file
name is CHGCAR or PARCHG.
The ABINIT package [132] for the ab initio computation of material properties has a
feature to output binary files storing electron densities, electrostatic potentials, and
wave functions. They can be converted into text files, *.xsf, having the XSF format of
XCryS-Den [133, 134] using a converter named Cut3D. The unit of electron densities is
3
bohr . Cut3D supports data Type 13 (XCrySDen/VENUS wavefunction real data),
whereby wave functions can be directly output in *.xsf.
Quantum ESPRESSO [135] for quantum simulations of materials has also a feature to
output files with the XSF format.
Detailed information about the XSF format can be obtained from
http://www.xcrysden.org/doc/XSF.html
In general, *.xsf consists of some sections that begins with keywords. VESTA reads in
volumetric data from a BEGIN_BLOCK_DATAGRID section. To make it possible to overlap
isosurfaces with a structural model, *.xsf should also contain (1) PRIMEVEC and
PRIMCOORD sections or (2) CONVVEC and CONVCOORD sections. Furthermore, lattice vectors
in the PRIMVEC or CONVVEC section must agree with spanning vectors in the
BEGIN_DATAGRID section. The unit of primitive lattice vectors (PRIMVEC) and Cartesian
coordinates is Å in the XSF format.
The entire information about data displayed currently are saved in a text file, *.vesta,
with the VESTA format. Files with the VESTA format contains all the structural data and
graphic settings whereas volumetric data are not directly recorded in *.vesta but
imported from external files. Directories and names of files storing volumetric data are
recorded in *.vesta as relative paths to the volumetric-data files from the directory of
*.vesta, which makes it possible to save the entire information about the current data in
a small file without duplicating huge volumetric data.
139
Structural data can be optionally imported from an external file in the same manner. To
prevent the structural data from being recorded directly in *.vesta, import structural data
using the Import Data dialog box with option “Link” (see 6.3.6). A name of file *.ffe
output by ORFFE [45] is also recorded in *.vesta so that geometrical parameters
recorded in *.ffe are automatically listed in the Geometrical Parameters dialog box after
reopening *.vesta (see 14.2). After Rietveld analysis with RIETAN-FP [12], lattice and
structure parameters in *.vesta are automatically updated provided that *.vesta and a
standard input file, *.ins, share the same folder.
2. Chem3D (*.cc1)
http://openbabel.org/docs/2.3.1/FileFormats/3D_viewer_Formats.html
140
8. VRML (*.wrl)
http://www.web3d.org/x3d/specifications/vrml/
9. XMol XYZ (*.xyz) http://en.wikipedia.org/wiki/XYZ_file_format
http://openbabel.org/docs/2.3.1/FileFormats/XYZ_cartesian_coordinates_format.html
Volumetric data
14. PRIMA binary format (*.pri)
141
19. Compressed volumetric-data format (*.m3d)
Two-dimensional data
Suppose that the 2D Data Display window has been opened to choose option “(hkl) plane
defined by two vectors” or “project along [hkl] axis” in the Create New Slice dialog box (see
15.5). Then, two-dimensional data shown in 2D Data Display window can be output by
selecting the 2D Data Display item from the File menu. Note that this function is disabled on
selection of the first option: “(hkl) plane in the bounding box.”
For example, electron densities on the (100) plane with a distance of 0.5d from the origin
in rutile-type TiO2 are output as follows:
TiO2
(hkl) : 1 0 0
Vector 1: 0 1 0
Vector 2: 0 0 1
Distance from the origin: 2.29696 2.29696
65 65
0.000000 0.000000 1.771768E-02
0.000000 0.046232 1.777942E-02
0.000000 0.092464 1.796342E-02
0.000000 0.138697 1.826639E-02
0.000000 0.184929 1.868208E-02
.....
The two integers (65 and 65) in the sixth line are numbers of divisions along the x and y axes.
Three data in each data line are X, Y , and density, where X and Y are, respectively, x and y
coordinates (unit: Å) from the origin in the two-dimensional plot.
2. EPS
9
3. JPEG
10
4. JPEG 2000
11
5. PNG
6. PPM
7. RAW
8. RGB (SGI)
9. TGA
10. TIFF
9
http://www.jpeg.org/jpeg/index.html
10 http://www.jpeg.org/jpeg2000/index.html
11 http://www.libpng.org/pub/png/
142
Images of structural and volumetric data displayed in the Graphics Area and the 2D Data
Display window can be recorded in files with a variety of graphics formats. When saving files
storing bit-mapped images, select “Export Raster Image...” in the “File” menu. The image
size of the exported graphics files is scaled up as a multiple of the size of the Graphics Area.
The scale factor is input after the file name has been entered in a file selection dialog box.
VESTA can export huge images that exceed the maximum width and height of the viewport,
which are displayed in the Text Area just after launching VESTA, e.g.,
Maximum supported width and height of the viewport: 4096 x 4096
The resolution of objects, i.e., the number of stacks and slices for atoms and bonds, are
also scaled by two factors {Scale} and {Increasing factor for stacks/slices} specified in the
Preferences dialog box (see 12.1.2, 12.1.3, and 16).
3. PostScript (PS)
Images of structural and volumetric data displayed in the Graphics Area of the Main Window
can be saved as vector images. When saving files storing vector images, select “Export
Vector Image...” under the “File” menu. There are some limitations on exporting vector
images. Colors of isosurfaces and sections colored according to volumetric data cannot be
represented with vector image formats. Translucent polygons of polyhedra, lattice planes,
and isosurfaces become opaque in all the formats except for PDF.
143
CLOSING REMARKS
Acknowledgments
We wish to thank Ruben A. Dilanian of Melbourne University for his dedication to a great deal
of programming for VICS and VEND, the two predecessors of VESTA. Akiji Yamamoto of
NIMS coded a highly efficient compressed volumetric-data format (*.m3d), and Takashi Ida of
Nagoya Institute of Technology wrote Voronoi integration code for volumetric data. Masao
Arai of NIMS and Masataka Mizuno of Osaka University kindly coded file converters,
wien2venus.py for WIEN2k and contrd for SCAT, respectively. Hiroshi Ohki of Shinshu
University helped us improve VICS, VEND and VICS-II significantly. Genta Sakane of
Okayama University of Science, Takuji Ikeda of AIST, and Kazuki Komatsu of the University
of Tokyo gave us many suggestions and helped us debug the program.
Finally, we thank for all who have helped us improve the program by reporting bugs and
offering valuable ideas.
144
Appendix A
Keyboard Shortcuts
Global (Menu)
Ctrl + N File menu - New Structure…
Ctrl + O File menu - Open…
Ctrl + S File menu - Save…
Ctrl + Shift + S File menu - Save as
Ctrl + W File menu - Close
Ctrl + Q File menu - Exit
Ctrl + E Edit menu - Edit Data - Phase…
Ctrl + B Edit menu - Bonds…
Ctrl + + View menu - Zoom In
Ctrl + – View menu - Zoom Out
Ctrl + Shift + A View menu - Overall Appearance…
Ctrl + Shift + P Objects menu - Properties…
Ctrl + Shift + B Objects menu - Boundary…
Ctrl + Shift + O Objects menu - Orientation…
When the Graphics Area has focus
r Change Manipulation to the Rotate mode
m Change Manipulation to the Magnify mode
t Change Manipulation to the Translate mode
s Change Manipulation to the Select mode
d Change Manipulation to the Distance mode
a Change Manipulation to the Angle mode
Left Rotate along y
Right Rotate along +y
Up Rotate along x
Down Rotate along +x
Ctrl + Up Rotate along z
Ctrl + Down Rotate along +z
Shift + Left Rotate along y (1/10 step)
Shift + Right Rotate along +y (1/10 step)
Shift + Up Rotate along x (1/10 step)
145
Table A.1: Continued from the previous page.
146
Bibliography
[2] M. Sakata and M. Sato, Acta Crystallogr., Sect. A: Found. Crystallogr., 46, 263 (1990).
[3] G. Bricogne and C. J. Gilmore, Acta Crystallogr., Sect. A: Found. Crystallogr., 46, 284
(1991).
[4] R. J. Papoular, Acta Crystallogr., Sect. A: Found. Crystallogr., 47, 293 (1991).
[5] F. Izumi and R. A. Dilanian, “Recent Research Developments in Physics,” Vol. 3, Part
II, Transworld Research Network, Trivandrum (2002), pp. 699–726.
[6] F. Izumi, Solid State Ionics, 172, 1 (2004).
[7] F. Izumi and K. Momma, IOP Conf. Ser.: Mater. Sci. Eng., 18, 022001 (2011).
[9] K. Momma, T. Ikeda, A. A. Belik, and F. Izumi, Powder Diffr., 28, 184 (2013).
[10] V. F. Sears, “International Tables for Crystallography,” Vol. C, 3rd ed., ed. by E. Prince,
Kluwer, Dordrecht (2004), pp. 444–452.
[11] K. Momma and F. Izumi, Commission on Crystallographic Computing, IUCr Newsletter,
No. 7, 106 (2006).
[12] F. Izumi and K. Momma, Solid State Phenom., 130, 15 (2007).
[16] F. Izumi and R. A. Dilanian, Commission on Powder Diffraction, IUCr Newsletter, No.
32, 59 (2005).
[17] F. Izumi, “Spectroscopy and Diffraction III,” Vol. 11: the 5th Series of Chemistry,
Maruzen, Tokyo (2006), p. 220.
[18] A. C. Larson and R. B. V. Dreele, “General Structure Analysis System (GSAS),” Report
LAUR 86-748, Los Alamos National Laboratory, Los Alamos (2004).
[19] J. Rodriguez-Carvajal, Phys. B (Amsterdam, Neth.), 192, 55 (1993).
[20] J. Smart, K. Hock, and S. Csomor, “Cross-Platform GUI Programming with wxWidgets,”
Prentice Hall (2005).
147
[21] W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, “Numerical Recipes:
the Art of Scientific Computing,” 3rd ed., Vol. 11: the 5th Series of Chemistry, Maruzen,
Tokyo (2007), pp. 1142–1146.
[22] B. Quentrec and C. Brot, J. Comput. Phys., 13, 430 (1973).
[27] I. Hassan, S. M. Antao, and J. B. Parise, Am. Mineral., 89, 359 (2004).
[28] T. Nakagawa, K. Kihara, and K. Harada, Am. Mineral., 86, 1506 (2001).
[32] W. H. Baur, Acta Crystallogr., Sect. B: Struct. Sci., 30, 1195 (1974).
[33] K. Robinson, G. V. Gibbs, and P. H. Ribbe, Science (Washington, DC, U. S.), 172, 567
(1971).
[34] R. Hoppe, Z. Kristallogr., 150, 23 (1979).
[39] N. E. Brese and M. O’Keeffe, Acta Crystallogr., Sect. B: Struct. Sci., 47, 192 (1991).
[40] I. D. Brown, “The Chemical Bond in Inorganic Chemistry,” Oxford University Press,
Oxford (2001).
[41] P. Blaha, K. Schwarz, G. Madsen, D. Kvasnicka, and J. Luitz, “WIEN2k, an Augmented
Plane Wave + Local Orbitals Program for Calculating Crystal Properties,” Technische
Universität Wien, Wien (2001).
[42] J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, and J. Akimitsu, Nature
(London, U. K.), 410, 63 (2001).
148
[43] S. Kumazawa, S. Yamamura, E. Nishibori, M. Takata, M. Sakata, F. Izumi, and Y. Ishii,
J. Phys. Chem. Solids, 60, 1407 (1999).
[44] P. Hartman, “Crystal Growth: An Introdution,” North-Holland Publishing Co., Amster-
dam (1973), pp. 367–402.
[45] W. R. Busing, K. O. Martin, and H. A. Levy, “A FORTRAN Crystallographic Function
and Error Program,” Report ORNL-TM-306, Oak Ridge National Laboratory, Oak Ridge
(1964).
[46] L. M. Gelato and E. Parthé, J. Appl. Crystallogr., 20, 139 (1987).
[48] V. G. Tsirelson, Acta Crystallogr., Sect. B: Struct. Sci., 58, 632 (2002).
[50] H. Arnold, “International Tables for Crystallography,” Vol. A, 5th ed., ed. by Th. Hahn,
Kluwer, Dordrecht (2006), pp. 78–85.
[51] N. V. Belov, N. N. Neronova, and T. S. Smirnova, Sov. Phys. Crystallogr., 2, 311 (1957).
[52] W. Opechowski and R. Guccione, “Magnetism,” ed. by G. T. Rado and H. Shut, Vol. II,
Part A, Academic Press, New York (1965), pp. 105–165.
[53] O. Zhou and D. E. Cox, J. Phys. Chem. Solids, 53, 1373 (1992).
[57] F. Izumi, “Spectroscopy and Diffraction III,” Vol. 11: the 5th Series of Chemistry,
Maruzen, Tokyo (2006), p. 161.
[58] H. Adachi, M. Tsukada, and C. Satoko, J. Phys. Soc. Jpn., 45, 875 (1978).
[59] A. T. N’Diaye, S. Bleikamp, P. J. Feibelman, and T. Michely, Phys. Rev. Lett., 97,
215501 (2006).
[60] K. Komatsu, T. Kuribayashi, A. Sano, E. Ohtani, and Y. Kudoh, Acta Crystallogr., Sect.
E: Struct. Rep. Online, 62, i216 (2006).
[61] A. Prakash, Acta Crystallogr., 22, 439 (1967).
[62] M. N. Burnett and C. K. Johnson, “ORTEP-III: Oak Ridge Thermal Ellipsoid Plot Pro-
gram for Crystal Structure Illustrations,” Report ORNL-6895, Oak Ridge National Labo-
ratory, Oak Ridge (1996).
[63] R. H. Gao, Acta Crystallogr., Sect. E: Struct. Rep. Online, 66, m487 (2010).
149
[66] M. Rakiah, J. Refebvre, O. Hernandez, W. van Beek, and M. Serpelloni, J. Appl.
Crystal-logr., 37, 766 (2004).
[67] R. M. Hazen, R. T. Downs, and C. T. Prewitt, “High-Temperature and High-Pressure
Crystal Chemistry,” Vol. 41: Reviews in Mineralogy, Mineral. Soc. Am., Washington,
D.C. (2000), Chap. 1.
[68] M. Avdeev, E. N. Caspi, and S. Yakovlev, Acta Crystallogr., Sect. B: Struct. Sci., 63,
363 (2007).
[69] A. Beran, E. Libowitzky, and T. Armbruster, Can. Mineral., 34, 803 (1996).
[85] E. Parthé and L. M. Gelato, Acta Crystallogr., Sect. A: Found. Crystallogr., 40, 169
(1984).
[86] K. Yvon, W. Jeitschko, and E. Parthé, J. Appl. Crystallogr., 10, 73 (1977).
150
[87] A. Belsky, M. Hellenbrandt, V. L. Karen, and P. Luksch, Acta Crystallogr., Sect. B:
Struct. Sci., B58, 364 (2002).
[88] P. Lightfoot, S. Pei, J. D. Jorgensen, Y. Yamada, T. Matsumoto, F. Izumi, and Y.
Kodama, Acta Crystallogr., Sect. C: Cryst. Struct. Commun., 47, 1143 (1991).
[89] V. K. Pecharsky and P. Y. Zavalij, “Fundamentals of Powder Diffraction and Structural
Characterization of Materials,” 2nd ed., Springer, New York (2009), pp. 447-450.
[90] C. Giacovazzo, “Fundamentals of Crystallography,” 3rd ed., ed. by C. Giacovazzo,
Oxford University Press, Oxford (2011), pp. 84–86.
[91] K. Kato and F. Izumi, unpublished work (2008).
[93] Z. Berant, R. Moreh, and S. Kahane, Phys. Lett. B, 69, 281 (1977).
[94] D. Waasmaier and A. Kirfel, Acta Crystallogr., Sect. A: Found. Crystallogr., 51, 416
(1995).
[95] C. T. Chantler, J. Phys. Chem. Ref. Data, 24, 71 (1995).
[97] G. M. Sheldrick, Acta Crystallogr., Sect. A: Found. Crystallogr., 64, 112 (2008).
[100] S. R. Hall, F. H. Allen, and I. D. Brown, Acta Crystallogr., Sect. A: Found. Crystallogr.,
47, 655 (1991).
[101] L. J. Farrugia, J. Appl. Crystallogr., 45, 849 (2012).
[102] F. H. Allen, Acta Crystallogr., Sect. B: Struct. Sci., 58, 380 (2002).
[104] J. Rehr and R. Albers, Phys. Rev. B: Condens. Matter Mater. Phys., 41, 8139 (1990).
[105] J. Rehr and R. Albers, Rev. Mod. Phys., 72, 621 (2000).
[108] F. Izumi and T. Ikeda, Mater. Sci. Forum, 321, 198 (2000).
[110] C. Glass, A. Oganov, and N. Hansen, Comp. Phys. Comm., 175, 713 (2006).
151
[113] K. Kawamura, Jpn. Chem. Prog. Exchange Newsletter, 6, No. 4, 91 (1995).
[117] H. Tanaka, M. Takata, E. Nishibori, K. Kato, T. Iishi, and M. Sakata, J. Appl. Crystal-
logr., 35, 282 (2002).
[118] M. Sale and M. Avdeev, J. Appl. Crystallogr., 45, 1054 (2012).
[120] G. Oszlányi and A. Sütő, Acta Crystallogr., Sect. A: Found. Crystallogr., 60, 134 (2004).
[121] G. Oszlányi and A. Sütő, Acta Crystallogr., Sect. A: Found. Crystallogr., 61, 147 (2005).
[122] G. Oszlányi and A. Sütő, Acta Crystallogr., Sect. A: Found. Crystallogr., 63, 156 (2006).
[123] A. T. Brünger, “X-Plor Version 3.1: A System for X-Ray Crystallography and NMR,”
Yale University Press, New Haven (1993).
[124] M. D. Segall, P. J. D. Lindan, M. J. Probert, C. J. Pickard, P. J. Hasnip, S. J. Clark, and
M. C. Payne, J. Phys.: Condens. Matter, 14, 2717 (2002).
152
[131] G. Kresse and J. Furthmüller, Comput. Mater. Sci., 6, 15 (1996).
[132] X. Gonze, J.-M. Beuken, R. Caracas, F. Detraux, M. Fuchs, G.-M. Rignanese, L. Sindic,
M. Verstraete, G. Zerah, F. Jollet, M. Torrent, A. Roy, M. Mikami, P. Ghosez, J.-Y.
Raty, and D. C. Allan, Comput. Mater. Sci., 25, 478 (2002).
[133] A. Kokalj, J. Mol. Graph. Mod., 17, 176 (1999).
153
Index
154
GAMESS, 138 MINCRYST, 132, 133
gamma, 105 MOLDA, 133
Gaussian, 138 MS-DOS, 132
general equivalent position, 22, 37, MXDORTO, 135
101 general grid, 128 MXDTRICL, 135
general position diagram, 86
general volumetric data, 141 neutron diffraction, 111
general volumetric-data format, 141 Niggli-reduced cell, 105
GLUI, 2 nuclear Thomson scattering, 111
GLUT, 2
occupancy, 46, 74, 111
gnuplot, 108, 114
opacity, 4, 9, 30, 88, 90, 91, 93
Graphic Area, 102
OpenGL, 2, 15–17, 99
GSAS, 132, 135, 136
ORFFE, 102, 140
GULP, 135
ORTEP-III, 60
hexagonal lattice, 104 OUTCAR, 139
oxidation state, 13, 77
ICSD, 13, 77, 104, 132
ICSD-CRYSTIN, 132 Igor PARCHG, 138, 139
Pro, 108, 114 interatomic Patterson function, 9, 135
distance, 9 interfacial Patterson-function density, 2
angle, 81 inversion PDB, 13, 133, 140
center, 103, 104 PDF, 143
isotropic atomic displacement parameter, 47 peak search, 115
periodic grid, 128
JPEG, 142 periodic volumetric-data format, 141
JPEG 2000, 142 PNG, 142
polar vector, 64
Laplacian, 115 polyhedra, 74
Laplacian of electron density, 115, polyhedral model, 4, 21, 28, 29,
135 lattice parameter, 36, 104 lattice 61 polyhedral volume ratio, 75
plane, 66 portable document format, 143
lattice-plane spacing, 47 POSCAR, 138, 140 POTCAR,
LAZY PULVERIX, 103, 104, 109 139
Line Profile, 114 powder diffraction pattern, 22, 108, 125, 130
Linux, 4, 16–18, 20, 117 PowderPlot, 109
LOCPOT, 138 PPN, 142
PRIMA, 2, 135
MacMolPlt, 138 macOS, 4,
primitive lattice, 42
16, 21, 23 MADEL, 22, 106,
primitive translation vector, 39
141 Madelung energy, 45,
principal axis, 2 protein data
106 magnetic moment, 4,
bank, 133
63 maximum-entropy
PS, 143
method, 1
maximum-entropy Patterson method, quadratic elongation, 9, 76
2 MDL Molfile, 132
mean square displacement, 2, RAW, 142
47 MEED, 135 reciprocal basis vector, 39
MEM-based pattern fitting, 1 reciprocal-lattice, 47, 71, 72
MEND, 135 redo, 3
metric tensor, 39, 40, 47 RETRIEVE, 132
155
RGB (SGI), 142 wxWidgets, 2, 14
rhombohedral lattice, 103
RIETAN-2000, 133, 140 X-PLOR, 137
RIETAN-FP, 13, 22, 81, 102, 108, 109, 125, XAFS, 132
130, 133, 134, 140 XANES, 132
RietPlot, 108 XCrySDen, 139
rotation matrix, 37, 38 XMol XYZ, 134, 141
XSF format, 139
SCAT, 134, 136, 144 XTL format, 135
SHELX-97, 112, 131
shift of the origin, 39
shininess, 67, 83, 88, 98 space-filling
model, 4, 6, 21, 27, 28, 83 specular,
67, 83
standard setting, 104, 105
standardization of crystal data,
103 standardization parameter,
104 static displacement, 63
stick model, 4, 21, 27–29, 87
STRUCTURE TIDY, 104, 105, 110
substructure, 42
Superflip, 137
superstructure, 42
surface coloring, 9, 22, 30, 49, 91,
139 SVG, 143
symmetry operation, 9, 37, 38, 40, 41, 74,
75, 79, 80, 101
temperature factor, 46
TGA, 142
TIFF, 142
transformation matrix, 38, 40
translation vector, 74
undo, 3
unit-cell volume, 38
156
USER’S MANUAL FOR MADEL
1. WHAT IS MADEL?
This manual was translated from Japanese into English with some modifications and additions
by Dr. Fujio Izumi of NIRIM during his short stay at the Argonne National Laboratory in
March, 1991.
1
Line 3. Three integers related with the crystal structure (3I3)
NS: number of symmetry operations. Inverted positions (–x, –y, –z) in a centrosymmetric
cell are not required, whereas translated positions in a complex lattice must be input.
NA: number of atoms contained in the asymmetric unit.
IZ: = 0, non-centrosymmetric cell.
= 1, centrosymmetric cell.
Column
1 - 11: translational part of xj or blank.
12 - 13: 1, 2, 3, –1, –2, –3, or blank for x, y, z, –x, –y, –z, or blank, respectively, as used in
the expression for the transformed xj.
14 - 15: 1, 2, 3, –1, –2, –3, or blank for x, y, z, –x, –y, –z, or blank, respectively, as used in
the expression for the transformed xj. Columns 12 and 13 are exactly equivalent to
Columns 14 and 15. Also, note that an expression such as xj = 2x must be treated as
xj = x + x.
16 - 26: translational part of yj or blank.
27 - 30: integers representing plus or minus x, y, or z in the expression for transformed yj as
described above.
31 - 41: translational part of zj or blank.
42 - 45: integers representing plus or minus x, y, or z in the expression for transformed zj as
described above.
These lines have the same format as those in ORFLS (least-squares structure-refinement
program) and ORFFE (program to calculate interatomic distances, bond angles, etc.)
developed at the Oak Ridge National Laboratory. Therefore, they can be copied easily from
input files for these two programs.
2
Lines 6. Sites in the asymmetric unit (A8, 1X, 5F9.4)
Input NA lines (I = 1, NA).
Lines 7. Fractional coordinates of vacant sites whose potentials are to be calculated (3F9.6)
Input lines as many as you like (I = NA+1, NA+2, .....). These lines are optional.
4. OUTPUT
The output of MADEL is described in German in the present version. The unit of the
electrostatic site potential, φi , for the ith site is e/Å (1 e/Å = 14.399652 V), where e is the
elementary charge (= 1.6021773×10 –19 C). The precision of results obtained using MADEL is
limited to 3 or 4 digits.
When the site input in Lines 7 is located within an ionic sphere, MADEL prints out a potential
(POBBA) excluding contribution of the sphere and, in addition, a potential (PMBBA)
calculated by substituting an original point charge for the sphere. If the site is not contained in
any ionic sphere, a potential in which contribution of the nearest neighbor (either one if two or
more nearest neighbors are present) is subtracted is output as POBBA.
The Madelung energy for the asymmetric unit, EM, is calculated by using the formula
3
1
E = ZW
M 2∑i φi i i
where Zi is the valence (oxidation state) of the ith site in the unit of e, and Wi is W(I). EM must
be multiplied by the number of general equivalent positions to obtain the Madelung energy for
the unit cell.
4
µ.5 2. 31% #.-20$ "-$ ,"*%#
$
¤&
#.-20
,"*%#
$ #.-20$ #
$
#.-20$ ¤
$4$)0
!¤6!
#.-20$ &
&
#
$
¤
&
¤
&
#
$ ¤ &
#
$ ¤ &
#
$
&
&
&
#.-20$ ¤ &
&
&
$.-20% $
%
$
%
#.)0
¤¤
!
3/ %.4-
§ '
-
-
¤
¤
3,#&0 .' /*5&+1 .32 .' 0"-(& $.-20%
'
µ !
µ !
# $-
#
$-
µ
-"! -
µ§
%§¤
"
1*#&. ,' "0,*/ 0, #& "+")(4&%
")) "0,*/
$
%
$,+0.%
§µ /$"0
§µ !/$"0 -
§µ !/$"0 ,3+
/$"0
/$"0
/$"0 1-
/$"0 %,3+
/$"0 /$"
/,/1.'"$&/
$,+0.%
' $
%
$,+0.%
µ§¤ %-"
!-%$# $