Coulomb Schrodinger Hamiltonian

Abstract When considering the work of Carl Ballhausen on vibrational spectra, it is suggested that his use of the Born-Oppenheimer approximation is capable of some refinement and extension in the light of later developments. A consideration of the potential energy surface in the context of a full Coulomb Schrodinger Hamiltonian in which translational and rotational motions are explicitly considered would seem to require a reformulation of the Born-Oppenheimer approach. The resulting potential surface for vibrational motion should be treated, allowing for the rotational motion and the nuclear permutational symmetry of the molecule. 

The relativistic many-electron Hamiltonian cannot be written in closed form it may be derived perturbatively from quantum electrodynamics [1]. The simplest form is the Dirac-Coulomb (DC) Hamiltonian, where the nonrelativistic one-electron terms in the Schrodinger equation are replaced by the one-electron Dirac operator hj). 

The corrections of order (Za) are just the first order matrix elements of the Breit interaction between the Coulomb-Schrodinger eigenfunctions of the Coulomb Hamiltonian Hq in (3.1). The mass dependence of the Breit interaction is known exactly, and the same is true for its matrix elements. These matrix elements and, hence, the exact mass dependence of the contributions to the energy levels of order (Za), beyond the reduced mass, were first obtained a long time ago [2]... 

The intrinsic failure of the Foldy-Wouthuysen protocol is therefore without doubt related to the ill-defined 1/c expansion of the kinetic term Ep, which does not bear any reference to the external potential V. However, in the literature the ill-defined behavior of the Foldy-Wouthuysen transformation has sometimes erroneously been attributed to the singular behavior of the Coulomb potential near the nucleus, and even the existence of the correct nonrelativistic limit of the Foldy-Wouthuysen Hamiltonian is sometimes the subject of dispute. Because of Eqs. (11.82) and (11.83) and the analysis given above, the nonrelativistic limit c —> oo, i.e., X —> 0 is obviously well defined, and for positive-energy solutions given by the Schrodinger Hamiltonian /nr = / 2me + V. 

The solute-solvent system, from the physical point of view, is nothing but a system that can be decomposed in a determined collection of electrons and nuclei. In the many-body representation, in principle, solving the global time-dependent Schrodinger equation with appropriate boundary conditions would yield a complete description for all measurable properties [47], This equation requires a definition of the total Hamiltonian in coordinate representation H(r,X), where r is the position vector operator for all electrons in the sample, and X is the position vector operator of the nuclei. In molecular quantum mechanics, as it is used in this section, H(r,X) is the Coulomb Hamiltonian[46]. The global wave function A(r,X,t) is obtained as a solution of the equation ... 

Technically, the time-independent Schrodinger equation (2) is solved for clamped nuclei. The Hamiltonian is broken into its electronic part, He, including the nuclear Coulomb repulsion energy, and the nuclear Hamiltonian HN. At this level, mass polarization effects are usually neglected. The wave function is therefore factorized as usual (r,X)= vP(r X)g(X). Formally, the electronic wave function d lnX) and total electronic energy, E(X), are obtained after solving the equation for each value of X ... 

Consequently, from the density the Hamiltonian can be readily obtained, and then every property of the system can be determined by solving the Schrodinger equation to obtain the wave function. One has to emphasize, however, that this argument holds only for Coulomb systems. By contrast, the density functional theory formulated by Hohenberg and Kohn is valid for any external potential. Kato s theorem is valid not only for the ground state but also for the excited states. Consequently, if the density n, of the f-th excited state is known, the Hamiltonian H is also known in principle and its eigenvalue problem ... 

The simple calculation of the matrix element of this Hamiltonian between the nonrelativistic Schrodinger-Coulomb wave functions gives the Fermi result [2] for the splitting between the l S i and states ... 

In principle, it should be possible to obtain the electronic energy levels of the molecules as a solution of the Schrodinger equation, if inter-electronic and internuclear cross-coulombic terms are included in the potential energy for the Hamiltonian. But the equation can be solved only if it can be broken up into equations which are functions of one variable at a time. A simplifying feature is that because of the much larger mass of the nucleus the motion of the electrons can be treated as independent of that of the nucleus. This is known as the Bom-Oppen-heimer approximation. Even with this simplification, the exact solution has been possible for the simplest of molecules, that is, the hydrogen molecule ion, H + only, and with some approximations for the H2 molecule. 

On the other hand, we sometimes discuss the Schrodinger equation without an absorption potential, for example, the Coulomb few-body Schrodinger equation with the Hamiltonian (3), assuming the energy to take complex values. This is analytic continuation of the quantal problem into fake, complex energies. This way, we depart from what actually occurs in nature for real energies. By regarding the nature as a special... 

It leaves intact the fermion operators related to the /1-th group itself. By virtue of this the two-electron operators WBA result in a renormalization of one-electron terms in the Hamiltonians for each group. <4 = 1,..., M. The expectation values ((b+b ))B are the one-electron densities. The Schrodinger equation eq. (1.193) can be driven close to the standard HFR form. This can be done if one defines generalized Coulomb and exchange operators for group A by their matrix elements in the carrier space of group A ... 

For the computational investigation of molecular systems containing heavy atoms, such as transition metals, lanthanides, and actinides, we could neglect neither relativity nor electron correlation. Relativistic effects, both spin-free and spin-orbit, increase with the nuclear charge of atoms. Therefore, instead of the nonrelativistic Schrodinger equation, we must start with the Dirac equation, which has four-component solutions. For many-electron systems, the four-component Hamiltonian is constructed from the one-electron Dirac operator with an approximated relativistic two-electron operator, such as the Coulomb, Breit, or Gaunt operator, within the nopair approximation. The four-component method is relativistically rigorous, which includes both spin-free and spin-orbit effects in a balanced way. However it requires much computational time since it contains more variational parameters than the approximated, one or two-component method. 

Solving the electronic Schrodinger equation for anything but the simplest systems is an extremely complex numerical problem, which relies heavily on the use of approximations. One of the most important of these is already present in the choice of the Hamiltonian, which in almost all computational studies is taken to include only the kinetic energy and Coulombic charge interaction terms, because these are usually quantitatively the most important (3) ... 

---