Defranceschi presentation

Numerical methods vs. basis sets in quantum chemistry: 

Numerical methods vs. basis sets in quantum chemistry M. Defranceschi CEA-Saclay

Various molecular ab initio models: 

Various molecular ab initio models Minimization problem where

For sake of simplicity: 

For sake of simplicity non relativistic equations, time-independent model, nuclei are points at fixed known positions, real-valued functions, spin not explicitly considered, simplifications to make it more convenient

Vector space considered: 

Vector space considered Physical functions  are square integrable functions (three dimensional measure in the Lebesgue sense) Hilbert space Reduced to a subspace

Overwhelming numerical difficulties: 

Overwhelming numerical difficulties Problem too difficult to be solved numerically Vector space too large Non linear terms

Two classes of simplification: 

Two classes of simplification Rigourous energy/approximated wavefunction Hartree-Fock approx. Restriction to a set of functions Rigourous density/approximated energy Density functional approx.

Hartree-Fock settings: 

Hartree-Fock settings

Mathematical fundation: 

Mathematical fundation Define the energy functional E() on a set X of functions  (the set of all the possible states of the molecule). Then find a function (the ground-state) satisfying some given constraint (i.e. constant number of electrons) and minimize the energy E on the convenient set of states :

Notion of physical space: 

Notion of physical space What are the variables of ? Physical notion : coordinates can be either position or momentum (or both) or any other quantity.

First ideas in position space: 

First ideas in position space Analytical solutions Numerical solutions Radial function of in a one-center approximation Spheroidal cooordinates for diatomic molecules Complete numerical integration for diatomic molecules In the case of atoms numerical integration are reliable

The quantum chemist procedure: 

The quantum chemist procedure Molecules are not considered as a whole but as constructed from atoms. Use of atomic basis sets Slater type Gaussian type Any functions which contain the correct physical information. The procedure most widely used consists in writing the molecular orbitals as LCAO which belong to a given complete set of the Sobolev space

Manageable approximate solutions: 

Manageable approximate solutions Infinite basis sets are impracticable Truncated basis sets Large expansions but often to small Tendency towards linear dependence Inherent deficiencies for GTO cusp problem wrong asymptotic decay

Some attempts of numerical solutions: 

Some attempts of numerical solutions Integration over a numerical grid Finite element method Momentum space direct numerical integration Numerical solution using a wavelet basis

Finite element method: 

Finite element method Very accurate results for even time-independent problems for simple systems. Large storage requirement for the FE matrices for extended three-dimensional systems Removal of the singularities inherent in the nuclear potentials.

Momentum space approach(1): 

Momentum space approach(1) In position space HF equations are integro-differential FT of operators and not of functions

Momentum space approach(2): 

Momentum space approach(2) In momentum space HF equations are first order integral equations

Momentum space approach(3): 

Momentum space approach(3) The solutions for bound states ( ) can be obtained by an iterative procedure starting with a LCAO in momentum space ( a modified Lanczos procedure). Enables to recover basis functions, and then basis sets not limited in size Enables to recover the asymptotic behavior. Removal of the singularities inherent in the nuclear and interelectronic potentials.

Momentum space approach(4): 

Momentum space approach(4)

Momentum space approach(5): 

Momentum space approach(5) Disadvantage of a FT of a function is that all information about its support or its singularities is lost. A function with high variations of momenta is hardly interpretable A compactly supported function requires a lot of sinusoidal functions

A wavelet approach(1): 

A wavelet approach(1) The idea is to realize a decomposition with vanishing functions which leads to a momentum representation involving a position parameter Functions depending on two variables linked to momentum and position are used

A wavelet approach(2): 

A wavelet approach(2) A representation is obtained by means of a decomposition of the Schrödinger operator onto an orthonormal wavelet basis. scaling function wavelet mother

A wavelet approach(3): 

A wavelet approach(3) The approach is related to multiresolution analysis, which is a decomposition of the Hilbert space into a chain of closed subspaces. The family defined by the scaling function constitutes an o.n. basis set for Vj. Let Wj be the space containing nthe difference in information between Vj-1 and Vj. It allows to decompose

A wavelet approach (4): 

A wavelet approach (4) The two part of a Fock operator has to be treated in two differents ways: The NS form of the Laplacian operator is solved iteratively The NS form of the potential term is obtained by a quadrature formula

A wavelet approach(5): 

A wavelet approach(5) The matrix representation of an operator applied to a vector may be depicted

A wavelet approach (6): 

A wavelet approach (6)

Conclusion: 

Conclusion The numerical development is far from the state of the art of the current quantum chemistry practice based on the use of atomic basis sets.