Molecular Schrodinger equation

Balint-Kurti G G, Dixon R N and Marston C C 1992 Grid methods for solving the Schrodinger equation and time-dependent quantum dynamics of molecular photofragmentation and reactive scattering processes/of. Rev. Phys. Chem. 11 317—44... 

Reactive atomic and molecular encounters at collision energies ranging from thermal to several kiloelectron volts (keV) are, at the fundamental level, described by the dynamics of the participating electrons and nuclei moving under the influence of their mutual interactions. Solutions of the time-dependent Schrodinger equation describe the details of such dynamics. The representation of such solutions provide the pictures that aid our understanding of atomic and molecular processes. 

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

Substitution of Eq. (3) into the molecular Schrodinger equation leads to a system of coupled equations in a coupled multistate electronic manifold... 

Various kinds of mixed quantum-classical models have been introduced in the literature. We will concentrate on the so-called quantum-classical molecular dynamics (QCMD) model, which consists of a Schrodinger equation coupled to classical Newtonian equations (cf. Sec. 2). 

It was stated above that the Schrodinger equation cannot be solved exactly for any molecular systems. However, it is possible to solve the equation exactly for the simplest molecular species, Hj (and isotopically equivalent species such as ITD" ), when the motion of the electrons is decoupled from the motion of the nuclei in accordance with the Bom-Oppenheimer approximation. The masses of the nuclei are much greater than the masses of the electrons (the resting mass of the lightest nucleus, the proton, is 1836 times heavier than the resting mass of the electron). This means that the electrons can adjust almost instantaneously to any changes in the positions of the nuclei. The electronic wavefunction thus depends only on the positions of the nuclei and not on their momenta. Under the Bom-Oppenheimer approximation the total wavefunction for the molecule can be written in the following form ... 

Much of quantum chemistry attempts to make more quantitative these aspects of chemists view of the periodic table and of atomic valence and structure. By starting from first principles and treating atomic and molecular states as solutions of a so-called Schrodinger equation, quantum chemistry seeks to determine what underlies the empirical quantum numbers, orbitals, the aufbau principle and the concept of valence used by spectroscopists and chemists, in some cases, even prior to the advent of quantum mechanics. 

Statistical mechanics is the mathematical means to calculate the thermodynamic properties of bulk materials from a molecular description of the materials. Much of statistical mechanics is still at the paper-and-pencil stage of theory. Since quantum mechanicians cannot exactly solve the Schrodinger equation yet, statistical mechanicians do not really have even a starting point for a truly rigorous treatment. In spite of this limitation, some very useful results for bulk materials can be obtained. 

Both molecular and quantum mechanics methods rely on the Born-Oppenheimer approximation. In quantum mechanics, the Schrodinger equation (1) gives the wave functions and energies of a molecule. 

The quantum mechanics methods in HyperChem differ in how they approximate the Schrodinger equation and how they compute potential energy. The ab initio method expands molecular orbitals into a linear combination of atomic orbitals (LCAO) and does not introduce any further approximation. 

These atomic orbitals, called Slater Type Orbitals (STOs), are a simplification of exact solutions of the Schrodinger equation for the hydrogen atom (or any one-electron atom, such as Li" ). Hyper-Chem uses Slater atomic orbitals to construct semi-empirical molecular orbitals. The complete set of Slater atomic orbitals is called the basis set. Core orbitals are assumed to be chemically inactive and are not treated explicitly. Core orbitals and the atomic nucleus form the atomic core. 

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 ... 

Molecular quantum mechanics finds the solution to a Schrodinger equation for an electronic Hamiltonian, Hgjg., that gives a total energy, Egjg(-(R) + V (R,R). Repeated solutions at different nuclear configurations, R, lead to some approximate potential energy sur-... 

Rather than solve a Schrodinger equation with the Nuclear Hamiltonian (above), a common approximation is to assume that atoms are heavy enough so that classical mechanics is a good enough approximation. Motion of the particles on the potential surface, according to the laws of classical mechanics, is then the subject of classical trajectory analysis or molecular dynamics. These come about by replacing Equation (7) on page 164 with its classical equivalent ... 

For small molecules, the accuracy of solutions to the Schrodinger equation competes with the accuracy of experimental results. However, these accurate ab initio calculations require enormous computation and are only suitable for the molecular systems with small or medium size. Ab initio calculations for very large molecules are beyond the realm of current computers, so HyperChem also supports semi-empirical quantum mechanics methods. Semi-empirical approximate solutions are appropriate and allow extensive chemical exploration. The inaccuracy of the approximations made in semi-empirical methods is offset to a degree by recourse to experimental data in defining the parameters of the method. Indeed, semi-empirical methods can sometimes be more accurate than some poorer ab initio methods, which require much longer computation times. 

The simplest approximation to the Schrodinger equation is an independent-electron approximation, such as the Hiickel method for Jt-electron systems, developed by E. Hiickel. Later, others, principally Roald Hoffmann of Cornell University, extended the Hiickel approximations to arbitrary systems having both n and a electrons—the Extended Hiickel Theory (EHT) approximation. This chapter describes some of the basics of molecular orbital theory with a view to later explaining the specifics of HyperChem EHT calculations. 

Now that you know the mathematical form, you can solve the independent-electron Schrodinger equation for the molecular orbitals. First substitute the LCAO form above into equation (47) on page 193, multiply on the left by and integrate to represent... 

It is not the intention that this book should be a primary reference on quantum mechanics such references are given in the bibliography at the end of this chapter. Nevertheless, it is necessary at this stage to take a brief tour through the development of the Schrodinger equation and some of its solutions that are vital to the interpretation of atomic and molecular spectra. 

A theoretical model should be uniquely defined for any given configuration of nuclei and electrons. This means that specifying a molecular structure is all that is required to produce an approximate solution to the Schrodinger equation no other parameters are needed to specify the problem or its solution. 

The orbital model would be exact were the electron repulsion terms negligible or equal to a constant. Even if they were negligible, we would have to solve an electronic Schrodinger equation appropriate to CioHs " " in order to make progress with the solution of the electronic Schrodinger equation for naphthalene. Every molecular problem would be different. 

Exact solutions to the electronic Schrodinger equation are not possible for many-electron atoms, but atomic HF calculations have been done both numerically and within the LCAO model. In approximate work, and for molecular applications, it is desirable to use basis functions that are simple in form. A polyelectron atom is quite different from a one-electron atom because of the phenomenon of shielding", for a particular electron, the other electrons partially screen the effect of the positively charged nucleus. Both Zener (1930) and Slater (1930) used very simple hydrogen-like orbitals of the form... 

In the Bom-Oppenheimer picture the nuclei move on a potential energy surface (PES) which is a solution to the electronic Schrodinger equation. The PES is independent of the nuclear masses (i.e. it is the same for isotopic molecules), this is not the case when working in the adiabatic approximation since the diagonal correction (and mass polarization) depends on the nuclear masses. Solution of (3.16) for the nuclear wave function leads to energy levels for molecular vibrations (Section 13.1) and rotations, which in turn are the fundamentals for many forms of spectroscopy, such as IR, Raman, microwave etc. 

Electronic structure methods are aimed at solving the Schrodinger equation for a single or a few molecules, infinitely removed from all other molecules. Physically this corresponds to the situation occurring in the gas phase under low pressure (vacuum). Experimentally, however, the majority of chemical reactions are carried out in solution. Biologically relevant processes also occur in solution, aqueous systems with rather specific pH and ionic conditions. Most reactions are both qualitatively and quantitatively different under gas and solution phase conditions, especially those involving ions or polar species. Molecular properties are also sensitive to the environment. 

The description of electronic distribution and molecular structure requires quantum mechanics, for which there is no substitute. Solution of the time-independent Schrodinger equation, Hip = Eip, is a prerequisite for the description of the electronic distribution within a molecule or ion. In modern computational chemistry, there are numerous approaches that lend themselves to a reasonable description of ionic liquids. An outline of these approaches is given in Scheme 4.2-1 [1] ... 

---