Ehrenfest

The initial classification of phase transitions made by Ehrenfest (1933) was extended and clarified by Pippard [1], who illustrated the distmctions with schematic heat capacity curves. Pippard distinguished different kinds of second- and third-order transitions and examples of some of his second-order transitions will appear in subsequent sections some of his types are unknown experimentally. Theoretical models exist for third-order transitions, but whether tiiese have ever been found is unclear. 

To add non-adiabatic effects to semiclassical methods, it is necessary to allow the trajectories to sample the different surfaces in a way that simulates the population transfer between electronic states. This sampling is most commonly done by using surface hopping techniques or Ehrenfest dynamics. Recent reviews of these methods are found in [30-32]. Gaussian wavepacket methods have also been extended to include non-adiabatic effects [33,34]. Of particular interest here is the spawning method of Martinez, Ben-Nun, and Levine [35,36], which has been used already in a number of direct dynamics studies. 

Direct dynamics attempts to break this bottleneck in the study of MD, retaining the accuracy of the full electronic PES without the need for an analytic fit of data. The first studies in this field used semiclassical methods with semiempirical [66,67] or simple Hartree-Fock [68] wave functions to heat the electrons. These first studies used what is called BO dynamics, evaluating the PES at each step from the elech onic wave function obtained by solution of the electronic structure problem. An alternative, the Ehrenfest dynamics method, is to propagate the electronic wave function at the same time as the nuclei. Although early direct dynamics studies using this method [69-71] restricted themselves to adiabatic problems, the method can incorporate non-adiabatic effects directly in the electionic wave function. 

The center of the wavepacket thus evolves along the trajectory defined by classical mechanics. This is in fact a general result for wavepackets in a hannonic potential, and follows from the Ehrenfest theorem [147] [see Eqs. (154,155) in Appendix C]. The equations of motion are straightforward to integrate, with the exception of the width matrix, Eq. (44). This equation is numerically unstable, and has been found to cause problems in practical applications using Morse potentials [148]. As a result, Heller inboduced the P-Z method as an alternative propagation method [24]. In this, the matrix A, is rewritten as a product of matrices... 

The standard semiclassical methods are surface hopping and Ehrenfest dynamics (also known as the classical path (CP) method [197]), and they will be outlined below. More details and comparisons can be found in [30-32]. The multiple spawning method, based on Gaussian wavepacket propagation, is also outlined below. See [1] for further infomiation on both quantum and semiclassical non-adiabatic dynamics methods. 

Both the BO dynamics and Gaussian wavepacket methods described above in Section n separate the nuclear and electronic motion at the outset, and use the concept of potential energy surfaces. In what is generally known as the Ehrenfest dynamics method, the picture is still of semiclassical nuclei and quantum mechanical electrons, but in a fundamentally different approach the electronic wave function is propagated at the same time as the pseudoparticles. These are driven by standard classical equations of motion, with the force provided by an instantaneous potential energy function... 

The MMVB force field has also been used with Ehrenfest dynamics to propagate trajectories using mixed-state forces [84]. The motivation for this is... 

Quantum chemical methods, exemplified by CASSCF and other MCSCF methods, have now evolved to an extent where it is possible to routinely treat accurately the excited electronic states of molecules containing a number of atoms. Mixed nuclear dynamics, such as swarm of trajectory based surface hopping or Ehrenfest dynamics, or the Gaussian wavepacket based multiple spawning method, use an approximate representation of the nuclear wavepacket based on classical trajectories. They are thus able to use the infoiination from quantum chemistry calculations required for the propagation of the nuclei in the form of forces. These methods seem able to reproduce, at least qualitatively, the dynamics of non-adiabatic systems. Test calculations have now been run using duect dynamics, and these show that even a small number of trajectories is able to produce useful mechanistic infomiation about the photochemistry of a system. In some cases it is even possible to extract some quantitative information. 

In a final step, we follow the ideas of Ehrenfest [252], who first looked for classical structures in the equations of quantum mechanics, and look at the time... 

Solving the Eqs. (C.6-C.8,C.12,C.13) comprise what is known as the Ehrenfest dynamics method. This method has appealed under a number of names and derivations in the literatnre such as the classical path method, eilconal approximation, and hemiquantal dynamics. It has also been put to a number of different applications, often using an analytic PES for the electronic degrees of freedom, but splitting the nuclear degrees of freedom into quantum and classical parts. 

Ehrenfest, R Bemerkung fiber die angenaherte Giiltigkeit der klassischen Mechanik innerhalb der Quantenmechanik. Z. fiir Physik 45 (1927) 455-457. 

For a second-order transition Eq. (4.53) is the analog of Eq. (4.5), which is useful for first-order transitions. Equation (4.53) and the first- and second-order terminology are due to Ehrenfest. 

The transition from a ferromagnetic to a paramagnetic state is normally considered to be a classic second-order phase transition that is, there are no discontinuous changes in volume V or entropy S, but there are discontinuous changes in the volumetric thermal expansion compressibility k, and specific heat Cp. The relation among the variables changing at the transition is given by the Ehrenfest relations. 

Compressibility and pressure dependence of Curie temperature are directly measured changes in specific heat and thermal expansion are calculated from the Ehrenfest relation. 

At the basis of the building up procedure used by Bohr there lies the adiabatic principle introduced by Ehrenfest [1917]. 

In some cases the hypothesis fixes completely which special motions are to be considered as allowed. This occurs if the new class of motions are derived by means of an adiabatic transformation from some class of motions already known (Ehrenfest [1917]). 

There are however some stringent restrictions on the applicability of the adiabatic principle (Mehra and Rechenberg [1982]). Ehrenfest himself [1917] showed that it was applicable to simply periodic systems. These are systems having two or more frequencies which are rational fractions of each other. In such systems the motion will necessarily repeat itself after a fixed interval of time. Burgers [1917], a student of Ehrenfest, showed that it was also applicable... 

Ehrenfest, P. [1917] Adiabatic Invariant and the Theory of Quanta , Philosophical Magazine, 33, 500. 

Classical dynamics is studied as a special case by analyzing the Ehrenfest theorem, coherent states (16) and systems with quasi classical dynamics like the rigid rotor for molecules (17) and the oscillator (18) for various particle systems and for EM field in a laser. 

Ehrenfest trajectory for three-dimensional D + H2 generated by the RWP method, that is, the modified Hamiltonian operator f H). Dotted curves in (b) correspond to the Ehrenfest trajectory determined by the usual Schrddinger equation. See text for further details. 

---