Theorem virial

The density scaling relations to be presented in Sect. 1.4, which constitute important constraints on the density functionals, are rooted in the same wave function scaling that will be used here to derive the virial theorem [26]. 

The density corresponding to the scaled wavefunction is the scaled density 

7 1 leads to densities n.y(r) that are higher (on average) and more contracted than n(r), while 7 1 produces densities that are lower and more expanded. 

Now consider what happens to H) = T- -F) under scaling. By definition of F, 

For example, n = —1 for the Hamiltonian of (1.36) in the presence of a single nucleus, or more generally when the Hellmann-Feynman forces of (1.37) vanish for the state F. 

---

The pressure follows from the virial theorem, or from the characteristic thennodynamic equation and the PF. It is given by... 

The pressure is usually calculated in a computer simulation via the virial theorem ol Clausius. The virial is defined as the expectation value of the sum of the products of the coordinates of the particles and the forces acting on them. This is usually written iV = X] Pxi where x, is a coordinate (e.g. the x ox y coordinate of an atom) and p. is the first derivative of the momentum along that coordinate pi is the force, by Newton s second law). The virial theorem states that the virial is equal to —3Nk T. 

The energy obtained from a calculation using ECP basis sets is termed valence energy. Also, the virial theorem no longer applies to the calculation. Some molecular properties may no longer be computed accurately if they are dependent on the electron density near the nucleus. 

In spite of its simplicity and the visual similarity of this equation to Eq. (7), we would like to note that Eq. (11) leads to a nontrivial thermodynamics of a partially quenched system in terms of correlation functions, see, e.g.. Ref. 25 for detailed discussion. Evidently, the principal route for and to the virial theorem is to exploit the thermodynamics of the replicated system. However, special care must be taken then, because the V and s derivatives do not commute. Moreover, the presence of two different temperatures, Pq and P, requires attention in taking temperature derivatives, setting those temperatures equal, if appropriate, only at the end of the calculations. 

In a simulation [19] the pressure tensor is obtained from the virial theorem [78]... 

The energy of a Is-electron in a hydrogen-like system (one nucleus and one electron) is —Z /2, and classically this is equal to minus the kinetic energy, 1/2 mv, due to the virial theorem E — —T = 1/2 V). In atomic units the classical velocity of a Is-electron is thus Z m= 1). The speed of light in these units is 137.036, and it is clear that relativistic effects cannot be neglected for the core electrons in heavy nuclei. For nuclei with large Z, the Is-electrons are relativistic and thus heavier, which has the effect that the 1 s-orbital shrinks in size, by the same factor by which the mass increases (eq. (8.2)). 

Any changes in the potential energy because of the Coulomb correlation must therefore also influence the kinetic energy. The virial theorem will be further discussed below. 

If a trial function 9 leads to a kinetic energy 1 and a potential energy Vx which do not fulfill the virial theorem (Eq. 11.15), the total energy (7 +Ei) is usually far from the correct result. Fortunately, there exists a very simple scaling procedure by means of which one can construct a new trial function which not only satisfies the virial theorem but also leads to a considerably better total energy. The scaling idea goes back to a classical paper by Hylleraas (1929), but the connection with the virial theorem was first pointed out by Fock.5 It is remarkable how many times this idea has been rediscovered and published in the modern literature. 

Let us first assume that (px is the exact solution. This implies that the variation principle must be fulfilled ior r] = 1, and substitution of this value into Eq. 11.24 leads then to the virial theorem in the form of Eq. 11.15. 

In the derivation above, we have included the kinetic energy of the nuclei in the Hamiltonian and considered a stationary state. In Eq. II.3, this term has been neglected, and we have instead assumed that the nuclei have given fixed positions. It has been pointed out by Slater34 that, if the nuclei are not situated in the proper equilibrium positions, the virial theorem will appear in a slightly different form. (A variational derivation has been given by Hirschfelder and Kincaid.11)... 

By multiplying Eq. 11.32 by 77 and by using Eq. 11.30 and Eq. 11.31, it is then easily checked that the virial theorem (Eq. 11.33) is satisfied for the scaled function internuclear distance R — rj xp. The distance R is here a simple function of p, and, after establishing the relationship in the form of a graph or a table, we can also solve the reverse problem of finding the properly scaled func-... 

We note that the virial theorem is automatically fulfilled in the Hartree-Fock approximation. This result follows from the fact that the single Slater determinant (Eq. 11.38) built up from the Hartree-Fock functions pk x) satisfying Eq. 11.46 is the optimum wave function of this particular form, and, since this wave function cannot be further improved by scaling, the virial theorem must be fulfilled from the very beginning. If we consider a stationary state with the nuclei in their equilibrium positions, we have particularly Thf = — Fhf, and for the correlation terms follows consequently that... 

Wigner s formula is open to criticism also on another point, since he assumes the existence of a stationary electron state where the density is so low that the kinetic energy may be neglected. This is in contradiction to the virial theorem (Eq. 11.15), which tells us that the kinetic energy can never be neglected in comparison to the potential energy and that the latter quantity is compensated by the former to fifty per cent. A reexamination of the low density case would hence definitely be a problem of essential interest. 

The reason for this complication of the theory is evident the truncated set may contain certain variable parameters, and, if these are carefully adjusted to render the best possible description of a specific state, they may become rather unsuitable for the description of another state. According to Section II.C(3), a truncated set should, e.g., always contain a scale factor as a variable parameter and, if this quantity is fitted to the ground state, it may give a basic set which is rather "out of scale for even the first excited state. Since the virial theorem is not satisfied for this state, the corresponding total energy may be comparatively poorly reproduced. This implies that in treating excited states, it is desirable to have reliable criteria for the accuracy of both energies and wave functions. 

If the basic set xpk is chosen complete, the virial theorem will be automatically fulfilled and no scaling is necessary. In such a case, the wave function under consideration may certainly be expressed in the form of Eq. III. 18, but, if the basis is chosen without particular reference to the physical conditions of the problem, the series of determinants may be extremely slowly convergent with a corresponding difficulty in interpreting the results. It therefore seems tempting to ask whether there exists any basic set of spin orbitals. which leads to a most "rapid convergency in the expansion, Eq. III. 18, of the wave function for a specific state (Slater 1951). 

In the case the calculations are based on a truncated set Wlf 2,. . . containing adjustable parameters, the A splitting is of particular importance, since it permits the investigator to use different values of these parameters for different eigenvalues Xk— the relation III.95 will anyway be valid. The scale factor rj is such a parameter, and the results in Section II.C(3) and III.D(lb) show that, by means of the A splitting, it is now possible to get the virial theorem exactly fulfilled for at least one of the eigenfunctions associated with each Xk. 

The results show that it is possible to improve the Hartree-Fock energy —2.86167 at.u. considerably by means of a simple correlation factor, but also that it is essential to scale the total function W properly to fulfil the virial theorem. The parameters in the best function u of the form of Eq. III. 121 are further given below ... 

Carr, W. J., Phys. Rev. 106, 414, Use of a general virial theorem with perturbation theory. ... 

Lowdin, P.-O., Scaling problem, virial theorem and connected relations in quantum mechanics."... 

The relationship E = —T = V /2) an example of the quantum-mechanical virial theorem. 

Show explicitly for a hydrogen atom in the Is state that the total energy is equal to one-half the expectation value of the potential energy of interaction between the electron and the nucleus. This result is an example of the quantum-mechanical virial theorem. 

We now consider more interesting properties that can be extracted in our approach which cannot be extracted in a standard X-ray charge analysis. For a system at equilibrium, the virial theorem gives the total energy as... 

Lembarki, A., F. Regemont, and H. Chermette. 1995. Gradient-corrected exchange potential with the correct asymptotic behavior and the corresponding exchange-energy functional obtained from virial theorem. Phys. Rev. A 52, 3704. 

In previous work, we have mainly used the DPM model to investigate the effects of the coefficient of normal restitution and the drag force on the formation of bubbles in fluidized beds (Hoomans et al., 1996 Li and Kuipers, 2003, 2005 Bokkers et al., 2004 Van der Floef et al., 2004), and not so much to obtain information on the constitutive relations that are used in the TFMs. In this section, however, we want to present some recent results from the DPM model on the excess compressibility of the solids phase, which is a key quantity in the constitutive equations as derived from the KTGF (see Section IV.D.). The excess compressibility y can be obtained from the simulation by use of the virial theorem (Allen and Tildesley, 1990). 

---