Dynamical method

More informative and precision instruments are devices based on the dynamic methods of indentation. The impedance of a vibrating probe perturbing the medium is related to the... 

The surface tension of a pure liquid should and does come out to be the same irrespective of the method used, although difficulties in the mathematical treatment of complex phenomena can lead to apparent discrepancies. In the case of solutions, however, dynamic methods, including detachment ones, often tend... 

It was made clear in Chapter II that the surface tension is a definite and accurately measurable property of the interface between two liquid phases. Moreover, its value is very rapidly established in pure substances of ordinary viscosity dynamic methods indicate that a normal surface tension is established within a millisecond and probably sooner [1], In this chapter it is thus appropriate to discuss the thermodynamic basis for surface tension and to develop equations for the surface tension of single- and multiple-component systems. We begin with thermodynamics and structure of single-component interfaces and expand our discussion to solutions in Sections III-4 and III-5. 

Both the Monte Carlo and the molecular dynamics methods (see Section III-2B) have been used to obtain theoretical density-versus-depth profiles for a hypothetical liquid-vapor interface. Rice and co-workers (see Refs. 72 and 121) have found that density along the normal to the surface tends to be a... 

Two parameters must be measured to apply the BET equation, the pressure at the sample and the amount adsorbed at this pressure. There are tlnee conmron methods for measuring the amount of gas adsorbed, called the volumetric method, the gravimetric method and the dynamic method, of which the volumetric method is the connnonest [21],... 

Nos e S 1984 A molecular dynamics method for simulations In the canonical ensemble Mol. Phys. 52 255-68... 

Nose S 1984 A unified formulation of the constant-temperature molecular dynamics methods J. Chem. Phys. 81 511-19... 

Parrinello M and Rahman A 1981 Polymorphic transitions In single crystals a new molecular dynamics method J. Appl. Phys. 52 7182-90... 

A comprehensive introduction to the field, covering statistical mechanics, basic Monte Carlo, and molecular dynamics methods, plus some advanced techniques, including computer code. 

Many experimental techniques now provide details of dynamical events on short timescales. Time-dependent theory, such as END, offer the capabilities to obtain information about the details of the transition from initial-to-final states in reactive processes. The assumptions of time-dependent perturbation theory coupled with Fermi s Golden Rule, namely, that there are well-defined (unperturbed) initial and final states and that these are occupied for times, which are long compared to the transition time, no longer necessarily apply. Therefore, truly dynamical methods become very appealing and the results from such theoretical methods can be shown as movies or time lapse photography. 

Final state analysis is where dynamical methods of evolving states meet the concepts of stationary states. By their definition, final states are relatively long lived. Therefore experiment often selects a single stationary state or a statistical mixture of stationary states. Since END evolution includes the possibility of electronic excitations, we analyze reaction products in terms of rovibronic states. 

In this chapter, we look at the techniques known as direct, or on-the-fly, molecular dynamics and their application to non-adiabatic processes in photochemistry. In contrast to standard techniques that require a predefined potential energy surface (PES) over which the nuclei move, the PES is provided here by explicit evaluation of the electronic wave function for the states of interest. This makes the method very general and powerful, particularly for the study of polyatomic systems where the calculation of a multidimensional potential function is an impossible task. For a recent review of standard non-adiabatic dynamics methods using analytical PES functions see [1]. 

A further model Hamiltonian that is tailored for the treatment of non-adiabatic systems is the vibronic coupling (VC) model of Koppel et al. [65]. This provides an analytic expression for PES coupled by non-adiabatic effects, which can be fitted to ab initio calculations using only a few data points. As a result, it is a useful tool in the description of photochemical systems. It is also very useful in the development of dynamics methods, as it provides realistic global surfaces that can be used both for exact quantum wavepacket dynamics and more approximate methods. 

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. 

By its nature, the application of direct dynamics requires a detailed knowledge of both molecular dynamics and quantum chemistry. This chapter is aimed more at the quantum chemist who would like to use dynamical methods to expand the tools at theh disposal for the study of photochemistry, rather than at the dynamicist who would like to learn some quantum chemishy. It hies therefore to introduce the concepts and problems of dynamics simulations, shessing that one cannot strictly think of a molecule moving along a trajectory even though this is what is being calculated. 

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 simplest way to add a non-adiabatic correction to the classical BO dynamics method outlined above in Section n.B is to use what is known as surface hopping. First introduced on an intuitive basis by Bjerre and Nikitin [200] and Tully and Preston [201], a number of variations have been developed [202-205], and are reviewed in [28,206]. Reference [204] also includes technical details of practical algorithms. These methods all use standard classical trajectories that use the hopping procedure to sample the different states, and so add non-adiabatic effects. A different scheme was introduced by Miller and George [207] which, although based on the same ideas, uses complex coordinates and momenta. 

The multiple spawning method described in Section IV.C has been applied to a number of photochemical systems using analytic potential energy surfaces. As well as small scattering systems [36,218], the large retinal molecule has been treated [243,244]. It has also been applied as a direct dynamics method. 

Full quantum wavepacket studies on large molecules are impossible. This is not only due to the scaling of the method (exponential with the number of degrees of freedom), but also due to the difficulties of obtaining accurate functions of the coupled PES, which are required as analytic functions. Direct dynamics studies of photochemical systems bypass this latter problem by calculating the PES on-the-fly as it is required, and only where it is required. This is an exciting new field, which requires a synthesis of two existing branches of theoretical chemistry—electronic structure theory (quantum chemistiy) and mixed nuclear dynamics methods (quantum-semiclassical). 

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. 

Nose, S. A molecular dynamics method for simulations in the canonical ensemble. Mol. Phys. 52 (1984) 255-268 ibid. A unified formulation of the constant temperature molecular dynamics method. J. Chem. Phys. 81 (1984) 511-519. 

Procacci, P., Darden, T., Marchi, M., A very fast molecular dynamics method to simulate biomolecular systems with realistic electrostatic interactions. J. Phys. Chem. 100 (1996) 10464-10468. 

Martyna, G.J. Adiabatic path integral molecular dynamics methods. I. Theory. J. Chem. Phys. 104 (1996) 2018-2027. 

R. Zhou and B. J. Berne. A new molecular dynamics method combining the reference system propagator algorithm with a fast multipole method for simulating proteins and other complex systems. J. Phys. Chem., 103 9444-9459, 1995. 

The niolcciilar dynamics method is useful for calculating the tint e-dependent properties of an isolated inoleciile. However, more often, one Is interested in th e properties of a molecule that is in ler-aclin. with other molecules. With IlyperC hem, yon can add solvent molecules to the simulation explicitly, but the addition of many solven t molecu les will make the sun u lation much slower. A faster so In Lion is to sim n late them otion of th e m olecu le of in Lercst n sin g Lan gevin dyn am ics. 

One of the main advantages of the stochastic dynamics methods is that dramatic tirn savings can he achieved, which enables much longer stimulations to he performed. Fc example, Widmalm and Pastor performed 1 ns molecular dynamics and stochastic dynamic simulations of an ethylene glycol molecule in aqueous solution of the solute and 259 vvatc jnolecules [Widmalm and Pastor 1992]. The molecular dynamics simulation require 300 hours whereas the stochastic dynamics simulation of the solute alone required ju 24 minutes. The dramatic reduction in time for the stochastic dynamics calculation is du not only to the very much smaller number of molecules present hut also to the fact the longer time steps can often he used in stochastic dynamics simulations. 

Lolecular dynamics methods that we have discussed in this chapter, and the examples ave been used to illustrate them, fall into the category of atomistic simulations, in... 

Zhou R and B J Berne 1995. A New Molecular Dynamics Method Combining the Reference Sys Propagator Algorithm with a Fast Multipole Method for Simulating Proteins and Ol Complex Systems. Journal of Chemical Physics 103 9444-9459. 

---