Hamilton

Figure Al.2.7. Trajectory of two coupled stretches, obtained by integrating Hamilton s equations for motion on a PES for the two modes. The system has stable anhamionic synmretric and antisyimnetric stretch modes, like those illustrated in figrne Al.2.6. In this trajectory, semiclassically there is one quantum of energy in each mode, so the trajectory corresponds to a combination state with quantum numbers nj = [1, 1]. The woven pattern shows that the trajectory is regular rather than chaotic, corresponding to motion in phase space on an invariant torus.

Hamilton C E, Bierbaum V M and Leone S R 1985 Product vibrational state distributions of thermal energy charge transfer reactions determined by laser-induced fluorescence in a flowing afterglow Ar" + CC -> CC (v= 0-6) + Ar J. Chem. Rhys. 83 2284-92... 

To generalize what we have just done to reactive and inelastic scattering, one needs to calculate numerically integrated trajectories for motions in many degrees of freedom. This is most convenient to develop in space-fixed Cartesian coordinates. In this case, the classical equations of motion (Hamilton s equations) are given... 

To develop an additional equation, we simply make the ansatz that the first temi on the left-hand side of equation (3.11.215h) equals the first temi on the right-hand side and similarly with the second temi. This innnediately gives us Hamilton s equations... 

Wright J C, Labuda M J, Zilian A, Chen P C and Hamilton J P 1997 New selective nonlinear vibrational spectroscopies J. Luminesc. 72-74 799-801... 

Chen P C, Hamilton J P, Zilian A, Labuda M J and Wright J C 1998 Experimental studies for a new family of infrared four wave mixing spectroscopies Appl. Spectrosc. 52 380-92... 

Hamilton D C, Mitchell A C and Nellis W J 1986 Electrical conductivity measurements in shock compressed liquid nitrogen Shock M/aves in Condensed Matter (Proc. 4th Am. Phys. Soc. Top. Conf.) p 473... 

Hamilton T P and Pulay P 1986 Direct Inversion In the Iterative subspace (DNS) optimization of open-shell, excited-state and small multiconfiguratlonal SCF wavefunctlons J. Chem. Phys. 84 5728... 

Sun Y-P, Wang P and Hamilton N B 1993 Fluorescence spectra and quantum yields of Buckminsterfullerene (Cgg) in room-temperature solutions. No excitation wavelength dependence J. Am. Chem. Soc. 115 6378-81... 

The method of molecular dynamics (MD), described earlier in this book, is a powerful approach for simulating the dynamics and predicting the rates of chemical reactions. In the MD approach most commonly used, the potential of interaction is specified between atoms participating in the reaction, and the time evolution of their positions is obtained by solving Hamilton s equations for the classical motions of the nuclei. Because MD simulations of etching reactions must include a significant number of atoms from the substrate as well as the gaseous etchant species, the calculations become computationally intensive, and the time scale of the simulation is limited to the... 

The total effective Hamiltonian H, in the presence of a vector potential for an A + B2 system is defined in Section II.B and the coupled first-order Hamilton equations of motion for all the coordinates are derived from the new effective Hamiltonian by the usual prescription [74], that is. 

Section VI shows the power of the modulus-phase formalism and is included in this chapter partly for methodological purposes. In this formalism, the equations of continuity and the Hamilton-Jacobi equations can be naturally derived in both the nonrelativistic and the relativistic (Dirac) theories of the electron. It is shown that in the four-component (spinor) theory of electrons, the two exha components in the spinor wave function will have only a minor effect on the topological phase, provided certain conditions are met (nearly nonrelativistic velocities and external fields that are not excessively large). 

The aim of this section is to show how the modulus-phase formulation, which is the keytone of our chapter, leads very directly to the equation of continuity and to the Hamilton-Jacobi equation. These equations have formed the basic building blocks in Bohm s formulation of non-relativistic quantum mechanics [318]. We begin with the nonrelativistic case, for which the simplicity of the derivation has... 

Again, the summation convention is used, unless we state otherwise. As will appear below, the same strategy can be used upon tbe Dirac Lagrangean density to obtain the continuity equation and Hamilton-Jacobi equation in the modulus-phase representation. 

Variationally deriving with respect to a leads to the Hamilton-Jacobi equation... 

The result of interest in the expressions shown in Eqs. (160) and (162) is that, although one has obtained expressions that include corrections to the nonrelativistic case, given in Eqs. (141) and (142), still both the continuity equations and the Hamilton-Jacobi equations involve each spinor component separately. To the present approximation, there is no mixing between the components. 

The terms before the square brackets give the nonrelativistic part of the Hamilton-Jacobi equation and the continuity equation shown in Eqs. (142) and (141), while the term with the squaie brackets contribute relativistic corrections. All terms from are of the nonmixing type between components. There are further relativistic terms, to which we now turn. 

In Eq. (168), the first, magnetic-field term admixes different components of the spinors both in the continuity equation and in the Hamilton-Jacobi equation. However, with the z axis chosen as the direction of H, the magnetic-field temi does not contain phases and does not mix component amplitudes. Therefore, there is no contribution from this term in the continuity equations and no amplitude mixing in the Hamilton-Jacobi equations. The second, electric-field term is nondiagonal between the large and small spinor components, which fact reduces its magnitude by a further small factor of 0 particle velocityjc). This term is therefore of the same small order 0(l/c ), as those terms in the second line in Eqs. (164) and (166) that refer to the upper components. 

This part of our chapter has shown that the use of the two variables, moduli and phases, leads in a direct way to the derivation of the continuity and Hamilton-Jacobi equations for both scalar and spinor wave functions. For the latter case, we show that the differential equations for each spinor component are (in the nearly nomelativistic limit) approximately decoupled. Because of this decoupling (mutual independence) it appears that the reciprocal relations between phases and moduli derived in Section III hold to a good approximation for each spinor component separately, too. For velocities and electromagnetic field strengths that ate nomrally below the relativistic scale, the Berry phase obtained from the Schrddinger equation (for scalar fields) will not be altered by consideration of the Dirac equation. 

As is well known. Molecular Dynamics is used to simulate the motions in many-body systems. In a typical MD simulation one first starts with an initial state of an N particle system F = xi,..., Xf,pi,..., pf) where / = 3N is the number of degrees of freedom in the system. After sampling the initial state one numerically solves Hamilton s equations of motion ... 

To perform MD simulation of a system with a finite number of degrees of freedom the Hamilton equations of motion... 

The explicit symplectic integrator can be derived in terms of free Lie algebra in which Hamilton equations (5) are written in the form... 

Caustics The above formulae can only be valid as long as Eq. (9) describes a unique map in position space. Indeed, the underlying Hamilton-Jacobi theory is only valid for the time interval [0,T] if at all instances t [0, T] the map (QOi4o) —> Q t, qo,qo) is one-to-one, [6, 19, 1], i.e., as long as trajectories with different initial data do not cross each other in position space (cf. Fig. 1). Consequently, the detection of any caustics in a numerical simulation is only possible if we propagate a trajectory bundle with different initial values. Thus, in pure QCMD, Eq. (11), caustics cannot be detected. 

The intriguing point about the second set of equations is that q is now kept constant. Thus the vector ip evolves according to a time-dependent Schrddinger equation with time-independent Hamilton operator H[q) and the update of the classical momentum p is obtained by integrating the Hellmann-Feynman forces [3] acting on the classical particles along the computed ip t) (plus a constant update due to the purely classical force field). 

Th c total en ergy of th e system. called the Hamilton iari, is ih c sum of th e kin elic an d poten tial energies (equation 24). 

---