First principles dynamics

Figure 3.31. (a) Experimental dissociation probability S0(= S) for D2 on Pt(l 11) as a function of Et and 0,. From Ref. [292]. (b) Points connected by lines are some of the experimental results of (a) re-plotted at fixed parallel incident energy Epar(= ). The pure solid, dashed and dotted lines are the equivalent results from 6D first principles dynamics. From Ref. [300]. 

Figure 3.36. Nitrogen dissociation on W(100). (a) Experimental measurements of the dissociation probability S as a function of En and Ts. (b) Experimental measurements of only the direct component of dissociation probability S as a function of Et and 6f. (a) and (b) from Ref. [339]. (c) Dissociation probability S from first principles classical dynamics, separated into a dynamic trapping fraction and a direct dissociation fraction, (d) Approximate reaction path for dynamic trapping mediated dissociation from the first principles dynamics. The numbers indicate the temporal sequence, (c) and (d) from Ref. [343].

In this author s opinion, the experimental study of gas-surface dynamics has unfortunately dropped off precipitously in the last several years. There are probably many reasons for this funding vagaries, the dominance of the STM in surface science today, the success of first principles dynamics, etc. However, it is absolutely essential to have the most detailed experiments to benchmark first principles theory. Comparisons for even the simple systems considered here demonstrate that theory is far from fully understood. Hopefully, a new generation of experimentalists/theorists will accept this challenge so that the next status report will make this one obsolete. 

An approach briefly presented here is based on a combination of the eikonal (or short wavelength) approximation for nuclei, and time-dependent Hartree-Fock states for the many-electron system, in what we have called the Eikonal/ TDHF approach.[13] A similar description can be obtained with narrow wavepackets for the nuclear motions. Several other approaches have recently been proposed for doing first principles dynamics, a very active area of current research. [39, 11, 15]... 

To consider next the partitioning approach for many-atom systems, we start again with the full density operator f(t), describing both electrons and nuclei. We describe a procedure useful in a first principles dynamics where the nuclear motions are described in an eikonal approximation with Q = Q(t), so that the electronic part of the problem involves the calculation of an electronic density operator r[Q(f), Q t), f] = pei t). This satisfies an electronic L-vN equation... 

This overview has also considered the advantages and disadvantages of a description using first principles dynamics. Its application to many-atom systems undergoing electronic transitions is a very active and challenging subject of molecular quantum dynamics. 

The demand for the research will cover the development of the novel algorithms utilizing parallel computation methods. The development of a hierarchical multi-scale paradigm will consolidate theoretical analysis and will lead to large-scale decision-making criteria of the process level design based on the first-principle dynamics. 

We have constructed a rigorous nonlinear first-principles dynamic model of this process with TMODS. We have used the model to test the control strategy and show that it does provides effective control of the vinyl acetate monomer process. 

In this paper we gave a dynamic extension of the DFT, by deriving a L-D equation (11) with the fluctuation-dissipation theorem (9). We showed that the stochastic equation correctly samples the density field according to the probability exp —jflf [n], (17), based on the second H-theorem (16). At this point we note however that our TO-DFT is phenomenological md it is desirable to have a first-principle dynamics generalization of DFT. 

Molecular dynamics is a true first principles dynamic molecular model. It simply solves the equations of motion. Given an intermolecular potential, MD provides the exact spatial and temporal evolution of the system. The stiffness caused by fast vibrations compared with slow molecular relaxations demands relatively small time steps and challenges current simulations. As an example, the time scale associated with vibrations is a fraction of a picosecond, whereas those associated with diffusion or reaction may easily be in the seconds to hours range depending on the activation energy. Consequently, MD on a single processor is usually limited to short time and length scales (e.g., pico- to nanoseconds and 1-2 nm). 

Fig 11.39 Structural changes observed during the reaction of a chlorine molecule with a silicon surface, (Figure redrawn from Stich I, A De Vita, M C Payne, M ] Gillan and L ] Clarke 1994 Surface Dissociation from First Principles Dynamics and Chemistry Physical Review B49 8076-8085 )... 

Ong MT, Leiding J, Tao H, Virshup AM, Martinez TJ (2009) First principles dynamics and minimum energy pathways for mechanochemical ring opening of cyclobutene. J Am Chem... 

Chapter 3 treats nuclear motions on the adiabatic potential energy surfaces (PES). One of the most powerful and simplest means to study chemical dynamics is the so-called ab initio molecular dynamics (or the first principle dynamics), in which nuclear motion is described in terms of the Newtonian d3mamics on an ab initio PES. Next, we review some of the representative time-dependent quantum theory for nuclear wavepackets such as the multiconfigurational time-dependent Hartree approach. Then, we show how such nuclear wavepacket d3mamics of femtosecond time scale can be directly observed with pump>-probe photoelectron spectroscopy. 

With respect to the dynamical properties of the hydrated electron in cluster systems, the first principle dynamics using ab initio molecular dynamics and so on have been extensively applied. [135, 180, 371, 408, 446] They revealed information about the structure and relative stabilities of the isomer clusters. Nonadiabatic dynamics of a solvated electron in various photochemical processes has also been studied experimentally. [62, 123, 294, 329] Rossky and co-workers [327, 468] also studied the relaxation dynamics of excess electrons using quantum molecular dynamics simulation techniques. Here the nonadiabatic interactions were taken into account basically within the scheme of surface hopping technique. [444]... 

One of the first approaches introduced and used in condensed phases, also in the sense of the historical development, was a biasing technique known as the Blue Moon [55]. This approach was better refined over the years [56] to become more user-friendly and easily implementable in a computer code. For all the details, we refer the reader to the cited original publications. Just to summarize the essential points, let us recall that in the case of first principles dynamical simulation, the method relies on the identification of a reaction coordinate, or order parameter, = (R/) of a given subset of atomic coordinates (Lagrangean variables) R/ able to track the activated process or chemical reaction on which one wants to focus. The simplest example is represented by the distance = R/—Ry between two atoms that are expected to form or break a chemical bond. This analytical function is added to, e.g., a Car-Parrinello Lagrangean as a holonomic constraint. 

---