Atomic boundary

Choose the atoms of interest for the sem i-empirical calculation, then use the Bctend to sp option on the. Select menu to establish the appropriate atomic boundaries for the c uantnm mechanics calculation. TTyperChem substitutes pararmeteri/ed pseudo-fluorine atom s for th e portion s of the molecule n ot included directly in the calculation (see the second part of this book, Theory and Methods). 

The density is a maximum in all directions perpendicular to the bond path at the position of a bond CP, and it thus serves as the terminus for an infinite set of trajectories, as illustrated by arrows for the pair of such trajectories that lie in the symmetry plane shown in Fig. 7.2. The set of trajectories that terminate at a bond-critical point define the interatomic surface that separates the basins of the neighboring atoms. Because the surface is defined by trajectories of Vp that terminate at a point, and because trajectories never cross, an interatomic surface is endowed with the property of zero-flux - a surface that is not crossed by any trajectories of Vp, a property made clear in Fig. 7.2. The final set of diagrams in Fig. 7.1 depict contour maps of the electron density overlaid with trajectories that define the interatomic surfaces and the bond paths to obtain a display of the atomic boundaries and the molecular structure. 

These definitions apply to any atomic system, molecule or crystal. Fig. 7.3 a illustrates their application to the charge distribution of the guanine-cytosine base-pair. Fig. 7.3 b shows the molecular structure defined by the bond paths and the associated CPs that clearly and uniquely define the three hydrogen bonds that link the two bases. Fig. 7.3 c shows the atomic boundaries and bond paths overlaid on the electron density in the plane of the nuclei. All properties of the atoms can be determined, enabling one, for example, to determine separately the energy of formation of each of the three hydrogen bonds. 

The thus obtained HF equations are solved with the additional (compared to the free atom) boundary conditions imposed on radial wavefunctions Pni(r) ... 

Total Energy for Heavy Neutral Atoms.—The fact that the neutral atom solution has the form I in Figure 1 implies that (x) - 0 as x - oo, in fact as 144/x which is readily verified to be an exact solution of the dimensionless TF equation (10), not however satisfying the atomic boundary condition (11). Since V(r) - 0 at infinity, it follows from equation (7) that, for this neutral case with N=Z, we must have p=0. The condition that, in the simplest density description of neutral atoms, the chemical potential is zero is important for the arguments which follow. We shall see below that one of the objectives of more sophisticated density descriptions must be to find p. Equations (25) and (26) can be rewritten in the form, using E= — T and p=0,... 

Such analysis can reveal bond paths (lines of maximum ED finking two atoms), exact location of the bond critical points, bond elhpticity (deviation from cylindrical symmetry), and so on. The real atomic boundaries can be found and the effective atomic charges integrated. 

In summary, electrons can be localized within an atomic boundary in a closed-shell interaction. They can also be separately localized within the inner or valence shells, but not to bonded regions which cross interatomic surfaces, nor, in general, to non-bonded regions which occupy only a portion of a valence shell, to yield the Lewis model of the localized electron pair. 

HBoundary is the boundary Hamiltonian and it will consist of two parts — an operator for the interaction with the QM atoms, Boundary(QMp and an energy for the interaction with the MM atoms, Boundary(MM)-... 

Replacing the covalent radii for Tq (see Sect. 1.4.1) and multiplying the right part of Eq. 1.24 by fc = 0.9 + 0.05n [37] in order to convert the radius of the maximum electron density Tq to the radius of its minimum, r, which can be regarded as the atomic boundary. The results are listed in Table SI.6. The difference between these approaches for calculating radii of free atoms to day allows to present the average values to within 0.1 A (Table 1.10). 

On the Atomic Scale. The most fundamental of structural composites are metal alloys. These materials do not appear to have the material interfaces that normally characterize structural composites. One could argue that these interfaces are provided by the atomic boundaries themselves, but the point is academic. The important feature is that the principle of additive strength begins at the atomic and molecular level and carries over to the bulk properties of the composite material. In alloys, the combination of two or more dis-tincdy different materials at the atomic level produces a new material with distinctive properties of its own. 

There is evidently a grave problem here. The wavefiinction proposed above for the lithium atom contains all of the particle coordinates, adheres to the boundary conditions (it decays to zero when the particles are removed to infinity) and obeys the restrictions = P23 that govern the behaviour of the... 

For many-electron systems such as atoms and molecules, it is obviously important that approximate wavefiinctions obey the same boundary conditions and symmetry properties as the exact solutions. Therefore, they should be antisynnnetric with respect to interchange of each pair of electrons. Such states can always be constmcted as linear combinations of products such as... 

A related advantage of studying crystalline matter is that one can have synnnetry-related operations that greatly expedite the discussion of a chemical bond. For example, in an elemental crystal of diamond, all the chemical bonds are equivalent. There are no tenninating bonds and the characterization of one bond is sufficient to understand die entire system. If one were to know the binding energy or polarizability associated with one bond, then properties of the diamond crystal associated with all the bonds could be extracted. In contrast, molecular systems often contain different bonds and always have atoms at the boundary between the molecule and the vacuum. 

Besides the expressions for a surface derived from the van der Waals surface (see also the CPK model in Section 2.11.2.4), another model has been established to generate molecular surfaces. It is based on the molecular distribution of electronic density. The definition of a Limiting value of the electronic density, the so-called isovalue, results in a boundary layer (isoplane) [187]. Each point on this surface has an identical electronic density value. A typical standard value is about 0.002 au (atomic unit) to represent electronic density surfaces. 

Here, the component of the autocorrelation vector a for the distance interval between the boundaries dj (lower) and (upper) is the sum of the products of property p for atoms i and j, respectively, having a Euclidian distance d within this interval. 

---