Atomic Calculations

In the case of atomic calculations, this choice of normalization does not need to be imposed. It results spontaneously from the form of the generahzed Sturmian basis set. 

To obtain the generalized Sturmian secular equations, we substitute the superposition (6) into the many-particle Schrodinger equation (5)  

Finally, multiplying from the left by a conjugate function in the basis set, integrating over space and spin coordinates, and making use of the potential-weighted orthonormality relations (13), we obtain the set of secular equations  

Notice that the kinetic energy term has vanished This remarkable feature of equation (16) results from the fact that we have chosen the energy of our isoenergetic basis set to be the same as the energy of the state that we are trying to represent. 

In the case of atoms, equation (3) can be solved exactly if we set Vq equal to the Coulomb attraction potential of the bare nucleus  

As our first example we will present several atomic calculations. These simple systems will allow us to gain a fist impression of the capabilities and limitations of DFT. To solve the Kohn-Sham equations we used the code of J. L. Martins [77]. The results are then compared to Hartree-Fock calculations performed with GAMESS [78]. As an approximation to the xc potential, we 

Our first many-electron example is argon. Argon is a noble gas with the closed shell configuration ls 2s 2p 3s 3p , so its ground-state is spherical. In Fig. 6.5 we plot the electron density for this atom as a function of the distance to the nucleus. The function n r) decays monotonically, with very little structure, and is therefore not a very elucidative quantity to behold. However, if we choose to represent r n(r), we can clearly identify the shell structure of the atom Three maxima, corresponding to the center of the three shells, 

Fernando Nogueira, Alberto Castro, and Miguel A.L. Marques 

The problem of the exponential decay can yet be seen from a different perspective. For a many-electron atom the Hartree energy can be written, in 

Note that in the sum the term with i = y is not excluded. This diagonal represents the interaction of one electron with itself, and is therefore called the self-interaction term. It is clearly a spurious term, and is exactly canceled by the diagonal part of the exchange energy. It is easy to see that neither the LDA nor the GGA exchange energy cancel exactly the self-interaction. This is, however, not the case in more sophisticated functionals like the exact exchange or the self-interaction-corrected LDA. 

In a free atom, however, the situation is much simpler owing to the spherical symmetry only a scalar component is relevant and it is gauge invariant since the position of the nucleus refers to the natural gauge origin. The magnetic susceptibility reduces to 

Some atomic calculations of magnetic susceptibilities based on the analytical Hartree-Fock functions were collected in [28] and are presented in Table 5.4. 

For heavy atoms, however, it is generally accepted that relativistic corrections are essential. The most important relativistic Hamiltonian terms appear to be the mass-velocity correction, the spin-orbit coupling and the Darwin 

Calculated diamagnetic susceptibilities —xdia (10-12 m3 mol-1) for neutral and charged atomic systems in ground states8 

It can be seen that the contributions of valence orbitals to the total second momentum (r2) are essential the inner shells contribute insignificantly. For example, within the series of Cu, Ag and Au atoms the contributions of valence nsi/2 orbitals are 40, 44 and 18% and those of (n — l)dy2 plus (n — l)d5/2 are 43, 37 and 50% of the total, respectively. The reduced contribution of the 6sl/2 orbital for the Au atom has roots in the relativistic contraction of this orbital, as is seen from the comparison with the non-relativistic limit. 

---

Basis sets can be constructed using an optimisation procedure in which the coefficients and the exponents are varied to give the lowest atomic energies. Some complications can arise when this approach is applied to larger basis sets. For example, in an atomic calculation the diffuse functions can move towards the nucleus, especially if the core region is described... 

Assuming a 2sf 2pf electron distribution for the carbon atoms, calculate the energy of fomiation of ethylene from the gaseous atoms. 

There are several types of basis functions listed below. Over the past several decades, most basis sets have been optimized to describe individual atoms at the EIF level of theory. These basis sets work very well, although not optimally, for other types of calculations. The atomic natural orbital, ANO, basis sets use primitive exponents from older EIF basis sets with coefficients obtained from the natural orbitals of correlated atom calculations to give a basis that is a bit better for correlated calculations. The correlation-consistent basis sets have been completely optimized for use with correlated calculations. Compared to ANO basis sets, correlation consistent sets give a comparable accuracy with significantly fewer primitives and thus require less CPU time. 

The contracted basis set created from the procedure above is listed in Figure 28.3. Note that the contraction coefficients are not normalized. This is not usually a problem since nearly all software packages will renormalize the coefficients automatically. The atom calculation rerun with contracted orbitals is expected to run much faster and have a slightly higher energy. 

The Schrodinger equation is a nonreiativistic description of atoms and molecules. Strictly speaking, relativistic effects must be included in order to obtain completely accurate results for any ah initio calculation. In practice, relativistic effects are negligible for many systems, particularly those with light elements. It is necessary to include relativistic effects to correctly describe the behavior of very heavy elements. With increases in computer capability and algorithm efficiency, it will become easier to perform heavy atom calculations and thus an understanding of relativistic corrections is necessary. 

Formal charge (Section 1 6) The charge either positive or negative on an atom calculated by subtracting from the number of valence electrons in the neutral atom a number equal to the sum of its unshared electrons plus half the elec trons in its covalent bonds... 

Since the basis set is obtained from atomic calculations, it is still desirable to scale exponents for the molecular environment. This is accomplished by defining an inner valence scale factor and an outer valence scale factor ( double zeta ) and multiplying the corresponding inner and outer a s by the square of these factors. Only the valence shells are scaled. 

Fig. 11. Abundance mass spectra of differently charged hot CgoLL clusters evaporating atoms calculated with a Monte-Carlo simulation (the Li and Cgg isotope distributions are included). Energies required to remove Li atoms were calculated using the MNDO method. The peaks at x = 12 and at x = 6 + n (where n is the cluster charge) observed in experiment (Fig. 9) are well reproduced.

The concentration dependence of the average electron transfer from the Zn to the Cu atoms calculated with the LSMS, SCF-KKR-CPA, and CPA-LSMS are shown in Fig. 2. The values obtained with the CPA-LSMS are almost the same as those from the LSMS, but... 

Use Figure 17-11 to estimate the resistivities of two metal samples, one made of pure copper and the other of a copper-manganese alloy containing one atom of manganese for every one hundred copper atoms. Calculate the ratio of the cost due to power loss from wire of the impure material to the cost due to the power loss from wire of the pure material. 

The volume per mole of atoms of some fourth-row elements (in the solid state) are as follows K, 45.3 Ca, 25.9 Sc, 18.0 Br, 23.5 and Kr, 32.2 ml/mole of atoms. Calculate the atomic volumes (volume per mole of atoms) for each of the fourth-row transition metals. Plot these atomic volumes and those of the elements given above against atomic numbers. 

I Meanwhile others object to the suggestion that the optimization of basis sets are carried out by reference to experimental data. While accepting that the exponents and contraction coefficients are generally optimized in atomic calculations, they insist that these optimizations are in themselves ab initio. 

For atoms with more than two electrons, it is very difficult to obtain such a small absolute error in the energy as in the helium case, but, within an isoelectronic sequence, the relative error will, of course, go down rapidly with increasing atomic number Z. The method of superposition of configurations has been used successfully in a number of applications, particularly by Boys (1950-) and Jucys (1947-), and, for a more detailed survey of the work on atoms, we will refer to the special table on atomic calculations in the bibliography. This is a field of rapid development, where one can expect important new results within the next few years. 

The starting point to obtain a PP and basis set for sulphur was an accurate double-zeta STO atomic calculation4. A 24 GTO and 16 GTO expansion for core s and p orbitals, respectively, was used. For the valence functions, the STO combination resulting from the atomic calculation was contracted and re-expanded to 3G. The radial PP representation was then calculated and fitted to six gaussians, serving both for s and p valence electrons, although in principle the two expansions should be different. Table 3 gives the numerical details of all these functions. 

The result of the kinetic energy operator on the basis function < )j is known from the atomic calculations. The remaining integrals, the overlap integrals... 

Following the most commonly used approach, that of truncating the cluster with hydrogen atoms, calculations were performed on an Si50i6Hi2 molecule, using a 3-2IG basis set. Four variants of this system were tested, and the... 

Let us now consider the velocity autocorrelation function (VACF) obtained from the MCYL potential, (namely, with the inclusion of vibrations). Figure 3 shows the velocity autocorrelation function for the oxygen and hydrogen atoms calculated for a temperature of about 300 K. The global shape of the VACF for the oxygen is very similar to what was previously determined for the MCY model. Very notable are the fast oscillations for the hydrogens relative to the oxygen. 

C07-0090. It requires 496 kJ/mol to break O2 molecules into atoms and 945 kJ/mol to break N2 molecules into atoms. Calculate the maximum wavelengths of light that can break these molecules apart. What part of the electromagnetic spectrum contains these photons ... 

Yakobi, H., Eliav, E. and Kaldor, U. (2007) Nuclear quadrupole moment of 197Au from high-accuracy atomic calculations. Journal of Chemical Physics, 126, 184305-1-184305-4. 

What is obviously needed is a generally accepted recipe for how atomic states should be dealt with in approximate density functional theory and, indeed, a few empirical rules have been established in the past. Most importantly, due to the many ways atomic energies can be obtained, one should always explicitly specify how the calculations were performed to ensure reproducibility. From a technical point of view (after considerable discussions in the past among physicists) there is now a general consensus that open-shell atomic calculations should employ spin polarized densities, i. e. densities where not necessarily... 

A further simplication often used in density-functional calculations is the use of pseudopotentials. Most properties of molecules and solids are indeed determined by the valence electrons, i.e., those electrons in outer shells that take part in the bonding between atoms. The core electrons can be removed from the problem by representing the ionic core (i.e., nucleus plus inner shells of electrons) by a pseudopotential. State-of-the-art calculations employ nonlocal, norm-conserving pseudopotentials that are generated from atomic calculations and do not contain any fitting to experiment (Hamann et al., 1979). Such calculations can therefore be called ab initio, or first-principles. ... 

The oxidation number, or oxidation state, is the formal charge on an atom calculated on the basis that it is in a wholly ionic compound. Oxidation numbers are assigned according to several rules. 

---