Schrodinger equation nuclear

In Section II, molecular dynamics within the BO approximation was introduced. As shown in Appendix A, the full nuclear Schrodinger equation is, however. 

The familiar BO approximation is obtained by ignoring the operators A completely. This results in the picture of the nuclei moving over the PES provided by the electrons, which are moving so as to instantaneously follow the nuclear motion. Another common level of approximation is to exclude the off-diagonal elements of this operator matrix. This is known as the Bom-Huang, or simply the adiabatic, approximation (see [250] for further details of the possible approximations and nomenclature associated with the nuclear Schrodinger equation). 

By using this expression for G, it is possible to write the nuclear Schrodinger equation (A.8) in matrix form [54,179] as... 

This last equation is the nuclear Schrodinger equation describing the motion of nuclei. The electronic energy computed from solving the electronic Schrodinger equation (3) on page 163 plus the nuclear-nuclear interactions Vjjjj(R,R) provide a potential for nuclear motion, a Potential Energy Surface (PES). 

Solutions to a Schrodinger equation for this last Hamiltonian (7) describe the vibrational, rotational, and translational states of a molecular system. This release of HyperChem does not specifically explore solutions to the nuclear Schrodinger equation, although future releases may. Instead, as is often the case, a classical approximation is made replacing the Hamiltonian by the classical energy ... 

This Hamiltonian is used in the Schrodinger equation for nuclear motion, describing the vibrational, rotational, and translational states of the nuclei. Solving the nuclear Schrodinger equation (at least approximately) is necessary for predicting the vibrational spectra of molecules. 

If we want to calculate the potential energy curve, then we have to change the intemuclear separation and rework the electronic problem at the new A-B distance, as in the H2 calculation. Once again, should we be so interested, the nuclear problem can be studied by solving the appropriate nuclear Schrodinger equation. This is a full quantum-mechanical equation, not to be confused with the MM treatment. 

PES), which is different for each electronic state of the system (i.e. each eigenfunction of the BO Schrodinger equation). Based on these PESs, the nuclear Schrodinger equation is solved to define, for example, the possible nuclear vibrational levels. This approach will be used below in the description of the nuclear inelastic scattering (NIS) method. 

Most of the AIMD simulations described in the literature have assumed that Newtonian dynamics was sufficient for the nuclei. While this is often justified, there are important cases where the quantum mechanical nature of the nuclei is crucial for even a qualitative understanding. For example, tunneling is intrinsically quantum mechanical and can be important in chemistry involving proton transfer. A second area where nuclei must be described quantum mechanically is when the BOA breaks down, as is always the case when multiple coupled electronic states participate in chemistry. In particular, photochemical processes are often dominated by conical intersections [14,15], where two electronic states are exactly degenerate and the BOA fails. In this chapter, we discuss our recent development of the ab initio multiple spawning (AIMS) method which solves the elecronic and nuclear Schrodinger equations simultaneously this makes AIMD approaches applicable for problems where quantum mechanical effects of both electrons and nuclei are important. We present an overview of what has been achieved, and make a special effort to point out areas where further improvements can be made. Theoretical aspects of the AIMS method are... 

A set of coupled equations for the evolution of the basis function coefficients is obtained by substituting the wavefunction ansatz of Eqs. (2.5)-(2.7) into the nuclear Schrodinger equation... 

Gaussian wavepacket propagation, 377-381 initial condition selection, 373-377 nuclear Schrodinger equation, 363-373 Adiabatic-to-diabatic transformation (ADT). 

In this section, the basic theory of molecular dynamics is presented. Starting from the BO approximation to the nuclear Schrodinger equation, the picture of nuclear dynamics is that of an evolving wavepacket. As this picture may be unusual to readers used to thinking about nuclei as classical particles, a few prototypical examples are shown. 

In a diabatic representation, the electronic wave functions are no longer eigenfunctions of the electronic Hamiltonian. The aim is instead that the functions are so chosen that the (nonlocal) non-adiabatic coupling operator matrix, A in Eq. (52), vanishes, and the couplings are represented by (local) potential operators. The nuclear Schrodinger equation is then written... 

The potential energy function U(R) that appears in the nuclear Schrodinger equation is the sum of the electronic energy and the nuclear repulsion. The simplest case is that of a diatomic molecule, which has one internal nuclear coordinate, the separation R of the two nuclei. A typical shape for U(R) is shown in Fig. 19.1. For small separations the nuclear repulsion, which goes like 1 /R, dominates, and liniR >o U(R) = oo. For large separations the molecule dissociates, and U(R) tends towards the sum of the energies of the two separated atoms. For a stable molecule in its electronic ground state U(R) has a minimum at a position Re, the equilibrium separation. 

Even for a diatomic molecule the nuclear Schrodinger equation is generally so complicated that it can only be solved numerically. However, often one is not interested in all the solutions but only in the ground state and a few of the lower excited states. In this case the harmonic approximation can be employed. For this purpose the potential energy function is expanded into a Taylor series about the equilibrium separation, and terms up to second order are kept. For a diatomic molecule this results in ... 

The first derivative vanishes since U(R) has a minimum at Re. Within this approximation the nuclear Schrodinger equation reduces to that of a harmonic oscillator, whose frequency to is given by ... 

Once this equation is solved for all relevant regions of the nuclear configuration space, in the BO framework, the nuclear motion can be treated either via a classical mechanical analysis with the help of computer simulations [6], or it can be treated quantum mechanically for simple models [54], In the latter scheme, the nuclear Schrodinger equation must be solved ... 

In order to compute the polarizability for the lowest rovibrational level of the EF state we have to solve the radial nuclear Schrodinger equation. 

---