ISSN:
0020-7608
Keywords:
Computational Chemistry and Molecular Modeling
;
Atomic, Molecular and Optical Physics
Source:
Wiley InterScience Backfile Collection 1832-2000
Topics:
Chemistry and Pharmacology
Notes:
A more general application of the self-consistent field iteration is coupled with a finite-difference Newton-Raphson algorithm to solve the set of coupled second-order integro-differential equations with split boundary conditions which constitutes the Hartree-Fock problem for diatomic molecules. The N orbitals are assumed to be of the form ψα = Lα(λ) Mα (μ)eimαφ (2π)-1/2, (α = 1, …, N), where λ, μ, and φ are the usual confocal elliptical coordinates. Requiring the expectation value of the electronic Hamiltonian to be stationary with respect to independent variations of the functions Lα and Mα, subject to constraints of orthonormality, leads to a set of coupled one-dimensional differential equations for the functions Lα and Mα. In the new method a corresponding set of finite-difference equations including the split boundary conditions for each function, as well as the Lagrange multipliers and associated constraints on normalization and orthogonality, are incorporated into a large system of nonlinear algebraic equations which is solved by means of a coupled self-consistent field-generalized Newton-Raphson iteration. As examples, calculations of the (1sσg)2 1Σg+ and (1sσg) (2pσu) 3Σu+ states of H2 are presented. The calculated energy for the 1Σg+ state of H2 is 99.985% of the three-dimensional Hartree-Fock limit. The discrepancy is due to the assumed factored form of the orbitals ψα, and a generalization of the finite-difference method is suggested to improve the results.
Additional Material:
2 Tab.
Type of Medium:
Electronic Resource
URL:
http://dx.doi.org/10.1002/qua.560090210
Permalink