Equilibrium atomic positions

Our intention is to give a brief survey of advanced theoretical methods used to detennine the electronic and geometric stmcture of solids and surfaces. The electronic stmcture encompasses the energies and wavefunctions (and other properties derived from them) of the electronic states in solids, while the geometric stmcture refers to the equilibrium atomic positions. Quantities that can be derived from the electronic stmcture calculations include the electronic (electron energies, charge densities), vibrational (phonon spectra), stmctiiral (lattice constants, equilibrium stmctiires), mechanical (bulk moduli, elastic constants) and optical (absorption, transmission) properties of crystals. We will also report on teclmiques used to study solid surfaces, with particular examples drawn from chemisorption on transition metal surfaces. 

In the expansion (A2.32), the first term is merely a constant, while the second one renormalizes equilibrium atomic positions but gives no contribution to the interaction of the atom C with a thermostat (provided a symmetric disposition of atoms, the term linear in r vanishes). The third term contains small corrections to... 

Owing to the complexity of zeolitic systems, most computational studies are still performed with the help of classical models. These methods use a set of potential functions to describe the potential energy surface (PES) in a manydimensional space of geometrical parameters of the system. Although the PES can be tested in terms of observables such as equilibrium atom positions, vibrational frequencies, heats of formation, and other experimental information, the PES itself is not an observable quantity. Because of that, there is no unique representation of the PES, and several coordinate systems and parameteriza-... 

Notably, the anisotropic thermal displacement factors form the elements of a 3x3 symmetric matrix. The physically meaningful form of this matrix when it is positive-definite is that of an ellipsoidal probability surface centered at the equilibrium atom position. An alternative form for Equation (22) frequently used in crystallography ... 

There is a problem with the rotational term. A rigid body (the equilibrium atomic positions Ooi are used), such as the benzene molecule, rotates, but due to symmetry, it may have some special axes characterizing the moments of inertia. The moment of inertia represents a tensor of rank 3 with the following components ... 

All of the techniques discussed so far indicate that the solid surface is ordered on an atomic scale. Most of the surface atoms occupy equilibrium atomic positions that are located in well-defined rows separated by equal interatomic positions. This atomic order is predominant despite the fact that there are large numbers of atomic positions on the surface where atoms have different numbers of neighbors. A pictorial representation of the topology of a monatomic crystal on an atomic scale is shown in... 

Broadening of spots can result from themial diffuse scattering and island fomiation, among other causes. The themial effects arise from the disorder in atomic positions as they vibrate around their equilibrium sites the sites themselves may be perfectly crystalline. 

HyperChem models the vibrations of a molecule as a set of N point masses (the nuclei of the atoms) with each vibrating about its equilibrium (optimized) position. The equilibrium positions are determined by solving the electronic Schrodinger equation. 

Fig. 6. Radiation damage in graphite showing the induced crystal dimensional strains. Impinging fast neutrons displace carbon atoms from their equilibrium lattice positions, producing an interstitial and vacancy. The coalescence of vacancies causes contraction in the a-direction, whereas interstitials may coalesce to form dislocation loops (essentially new graphite planes) causing c-direction expansion.

Lateral density fluctuations are mostly confined to the adsorbed water layer. The lateral density distributions are conveniently characterized by scatter plots of oxygen coordinates in the surface plane. Fig. 6 shows such scatter plots of water molecules in the first (left) and second layer (right) near the Hg(l 11) surface. Here, a dot is plotted at the oxygen atom position at intervals of 0.1 ps. In the first layer, the oxygen distribution clearly shows the structure of the substrate lattice. In the second layer, the distribution is almost isotropic. In the first layer, the oxygen motion is predominantly oscillatory rather than diffusive. The self-diffusion coefficient in the adsorbate layer is strongly reduced compared to the second or third layer [127]. The data in Fig. 6 are qualitatively similar to those obtained in the group of Berkowitz and coworkers [62,128-130]. These authors compared the structure near Pt(lOO) and Pt(lll) in detail and also noted that the motion of water in the first layer is oscillatory about equilibrium positions and thus characteristic of a solid phase, while the motion in the second layer has more... 

There are three different approaches to a thermodynamic theory of continuum that can be distinguished. These approaches differ from each other by the fundamental postulates on which the theory is based. All of them are characterized by the same fundamental requirement that the results should be obtained without having recourse to statistical or kinetic theories. None of these approaches is concerned with the atomic structure of the material. Therefore, they represent a pure phenomenological approach. The principal postulates of the first approach, usually called the classical thermodynamics of irreversible processes, are documented. The principle of local state is assumed to be valid. The equation of entropy balance is assumed to involve a term expressing the entropy production which can be represented as a sum of products of fluxes and forces. This term is zero for a state of equilibrium and positive for an irreversible process. The fluxes are function of forces, not necessarily linear. However, the reciprocity relations concern only coefficients of the linear terms of the series expansions. Using methods of this approach, a thermodynamic description of elastic, rheologic and plastic materials was obtained. 

A (Figure 4.9). The diameter of such a neck, 2.3 A, is sufficiently large for a linear C-C chain to pass, but too small to also be an equilibrium adsorption position. The largest compound allowed inside the pores is a linear molecule limited in length to four carbon atoms due to the distance between two subsequent necks [103]. Another example of shape-selective behavior is found in a Zn-based MOF able to encapsulate linear hexane while branched hexanes are blocked [104]. 

The atomic temperature factor, or B factor, measures the dynamic disorder caused by the temperature-dependent vibration of the atom, as well as the static disorder resulting from subtle structural differences in different unit cells throughout the crystal. For a B factor of 15 A2, displacement of an atom from its equilibrium position is approximately 0.44 A, and it is as much as 0.87 A for a B factor of 60 A2. It is very important to inspect the B factors during any structural analysis a B factor of less than 30 A2 for a particular atom usually indicates confidence in its atomic position, but a B factor of higher than 60 A2 likely indicates that the atom is disordered. 

With such calculations one can approach Hartree-Fock accuracy for a particular cluster of atoms. These calculations yield total energies, and so atomic positions can be varied and equilibrium positions determined for both ground and excited states. There are, however, drawbacks. First, Hartree-Fock accuracy may be insufficient, as correlation effects beyond Hartree-Fock may be of physical importance. Second, the cluster of atoms used in the calculation may be too small to yield an accurate representation of the defect. And third, the exact evaluation of exchange integrals is so demanding on computer resources that it is not practical to carry out such calculations for very large clusters or to extensively vary the atomic positions from calculation to calculation. Typically the clusters are too small for a supercell approach to be used. 

The statistical treatment of the vibrational degrees of freedom of crystals is far more difficult compared to gases. Let us initially consider a monoatomic crystal. An atom in a crystal vibrates about its equilibrium lattice position. In the simplest approach, three non-interacting superimposed linear harmonic oscillators represent the vibrations of each atom. The total energy, given by the sum of the kinetic and potential energies for the harmonic oscillators, is... 

A first impression of collective lattice vibrations in a crystal is obtained by considering one-dimensional chains of atoms. Let us first consider a chain with only one type of atom. The interaction between the atoms is represented by a harmonic force with force constant K. A schematic representation is displayed in Figure 8.4. The average interatomic distance at equilibrium is a, and the equilibrium rest position of atom n is thus un =na. The motion of the chain of atoms is described by the time-dependent displacement of the atoms, un(t), relative to their rest positions. We assume that each atom only feels the force from its two neighbours. The resultant restoring force (F) acting on the nth atom of the one dimensional chain is now in the harmonic approximation... 

Let us describe the displacement of the nth atom from its equilibrium rest position by a cosine-wave with amplitude w0, angular frequency and wave vector q = 2n/X 1... 

Take an N-atomic molecule with the nuclei each at their equilibrium internuclear position. Establish a Cartesian x, y, z coordinate system for each of the nuclei such that, for Xj with i = 1,.., 3N, xi is the Cartesian x displacement of nucleus 1, x2 is the Cartesian y displacement coordinate for nucleus 1, X3 is the Cartesian z displacement coordinate for nucleus 1,..., x3N is the Cartesian z displacement for nucleus N. Use of one or another quantum chemistry program yields a set of force constants I ij in Cartesian displacement coordinates... 

---