Schrodinger equation basis sets

This makes it desirable to define other representations in addition to the electronically adiabatic one [Eqs. (9)-(12)], in which the adiabatic electronic wave function basis set used in the Bom-Huang expansion (12) is replaced by another basis set of functions of the electronic coordinates. Such a different electronic basis set can be chosen so as to minimize the above mentioned gradient term. This term can initially be neglected in the solution of the / -electionic-state nuclear motion Schrodinger equation and reintroduced later using perturbative or other methods, if desired. This new basis set of electronic wave functions can also be made to depend parametrically, like their adiabatic counterparts, on the internal nuclear coordinates q that were defined after Eq. (8). This new electronic basis set is henceforth refened to as diabatic and, as is obvious, leads to an electronically diabatic representation that is not unique unlike the adiabatic one, which is unique by definition. 

U(qJ is referred to as an adiabatic-to-diabatic transformation (ADT) matrix. Its mathematical sbucture is discussed in detail in Section in.C. If the electronic wave functions in the adiabatic and diabatic representations are chosen to be real, as is normally the case, U(q ) is orthogonal and therefore has n n — l)/2 independent elements (or degrees of freedom). This transformation mabix U(qO can be chosen so as to yield a diabatic electronic basis set with desired properties, which can then be used to derive the diabatic nuclear motion Schrodinger equation. By using Eqs. (27) and (28) and the orthonormality of the diabatic and adiabatic electronic basis sets, we can relate the adiabatic and diabatic nuclear wave functions through the same n-dimensional unitary transformation matrix U(qx) according to... 

When the wave function is completely general and pennitted to vary in the entire Hilbert space the TDVP yields the time-dependent Schrodinger equation. However, when the possible wave function variations are in some way constrained, such as is the case for a wave function restricted to a particular functional form and represented in a finite basis, then the corresponding action generates a set of equations that approximate the time-dependent Schrodinger equation. 

For both types of orbitals, the coordinates r, 0, and (j) refer to the position of the electron relative to a set of axes attached to the center on which the basis orbital is located. Although Slater-type orbitals (STOs) are preferred on fundamental grounds (e.g., as demonstrated in Appendices A and B, the hydrogen atom orbitals are of this form and the exact solution of the many-electron Schrodinger equation can be shown to be of this form (in each of its coordinates) near the nuclear centers), STOs are used primarily for atomic and linear-molecule calculations because the multi-center integrals < XaXbl g I XcXd > (each... 

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. 

HyperChem s ab initio calculations solve the Roothaan equations (59) on page 225 without any further approximation apart from the use of a specific finite basis set. Therefore, ab initio calculations are generally more accurate than semi-empirical calculations. They certainly involve a more fundamental approach to solving the Schrodinger equation than do semi-empirical methods. 

The cell in the lower right corner of the chart represents the exact solution of the Schrodinger equation, the limit toward which all approximate methods strive. Full Cl using an infinitely flexible basis set is the exact solution. ... 

As the basis set becomes infinitely flexible, full Cl approaches the exact solution of the time-independent, non-relativistic Schrodinger equation. 

The K-matrix method is essentially a configuration interaction (Cl) performed at a fixed energy lying in the continuum upon a basis of "unperturbed funetions that (at the formal level) includes both diserete and eontinuous subsets. It turns the Schrodinger equation into a system of integral equations for the K-matrix elements, which is then transformed into a linear system by a quadrature upon afinite L basis set. 

We should also mention that basis sets which do not actually comply with the LCAO scheme are employed under certain circumstances in density functional calculations, i. e., plane waves. These are the solutions of the Schrodinger equation of a free particle and are simple exponential functions of the general form... 

The TD wavefunction satisfying the Schrodinger equation ih d/dt) F(f) = // (/,) can be expanded in a basis set whose elements are the product of the translational basis of R, vibrational wavefunctions for r, r2, and the body-fixed (BF) total angular momentum eigenfunctions as41... 

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

We have investigated another procedure for reducing the computational expense of the AIMS method, which capitalizes on the temporal nonlocality of the Schrodinger equation and the deterministic aspect of the AIMS method. Recall that apart from the Monte Carlo procedure that we employ for selecting initial conditions, the prescription for basis set propagation and expansion is deterministic. We emphasize the deterministic aspect because the time-displaced procedure relies on this property. 

If the PES are known, the time-dependent Schrodinger equation, Eq. (1), can in principle be solved directly using what are termed wavepacket dynamics [15-18]. Here, a time-independent basis set expansion is used to represent the wavepacket and the Hamiltonian. The evolution is then carried by the expansion coefficients. While providing a complete description of the system dynamics, these methods are restricted to the study of typically 3-6 degrees of freedom. Even the highly efficient multiconfiguration time-dependent Hartree (MCTDH) method [19,20], which uses a time-dependent basis set expansion, can handle no more than 30 degrees of freedom. 

Ab initio quantum mechanics is based on a rigorous treatment of the Schrodinger equation (or equivalent matrix methods)4-7 which is intellectually satisfying. While there are a number of approximations made, it relies on a set of equations and a few physical constants.8 The use of ab initio methods on large systems is limited if not impossible, even with the fastest computers available. Since the size of an ab initio calculation is defined by the number of basis functions in the system, ab initio calculations are extremely costly for anything past the second row in the periodic table, and for all systems with more than 20 or 30 total atoms. 

A numerical solution of the Schrodinger equation in Eq. [1] often starts with the discretization of the wave function. Discretization is necessary because it converts the differential equation to a matrix form, which can then be readily handled by a digital computer. This process is typically done using a set of basis functions in a chosen coordinate system. As discussed extensively in the literature,5,9-11 the proper choice of the coordinate system and the basis functions is vital in minimizing the size of the problem and in providing a physically relevant interpretation of the solution. However, this important topic is out of the scope of this review and we will only discuss some related issues in the context of recursive diagonalization. Interested readers are referred to other excellent reviews on this topic.5,9,10... 

In the variational calculations, the expansions of the kinetic energy factors and the pseudo-potential 1/ [1] are taken to fourth order, and, as mentioned above, the potential energy V is expanded through sixth order. In the numerical integration of the inversion Schrodinger equation a grid of 1000 points is used. The basis set [1] is truncated so that... 

---