The VMFCI Method explained
VMFCI is the acronym for Vibrational Mean Field Configuration Interactions. The method consists in performing vibrational configuration interactions of some degrees of freedom (dof) in the field of the other. The aim of the method is to keep the size of the finite basis sets used in successive VCI manageable by contracting groups of dof together, as in the traditional contraction method [1]. However, the power of the method comes from the mean field term added to the group Hamiltonians not present in previous approaches.
A VMFCI step starts with a partition, , of the nvib vibrational degrees of freedom into np subsets :
Using partition the vibrational Hamiltonian can be written as :
where
denotes a vibrational operator which only depends upon operators acting on the degrees in subsets I.
Then, one defines a possibly coarser partition, :
For such a step,
we call contractions the subsets J and components of contraction J
, the subsets I
such that :
Notes
- The components of one step are the contractions of the previous step.
is allowed: Iterating the same partition until self-consistency lead to VSCFCI methods which generalizes
the well-known VSCF method [2,3]. Such a generalization has been considered by Bowman and Gazdi [4] but not in the frame of an iterated VCI approach.
------
Let us consider one given contraction J . We assume that contraction J
has
components :
Having a basis set spanning an Hilbert subspace of dimension, say d, for each component I
:
we build for contraction J, a so-called product basis set, spanning an Hilbert subspace of dimension D
,
by constructing product functions of the form:
or more explicitly using variables dependencies:
Notes
- In CONVIV, the process is initialized with an initial basis set of modals, that is to say, functions of a single vibrational degree of freedom, but this constraint can be walked around if groups of dof have to be contracted from the start. The modals can be eigenfunctions of one dimensional Schrödinger equations with harmonic potential (with arbitrary center and frequency), Morse potential, Trigonometric Pösch-Teller potential, Kratzer potential, or can be Chebychev polynomials.
- By convention, the eigenstates are numbered in increasing order with natural numbers starting with 0, so the vector index M
=(0, ... ,0) always corresponds to product of ground state functions.
- The dimension of the basis set for the contraction J
can be different from the product of dimensions of its component's basis sets because of possible basis function truncations, usually performed according to some energy criteria (see infra).
------
For the contraction J, we define a partial Hamiltonian, H
, by grouping all the terms in Hvib involving the degrees in components Ijl ofJ
:
In contrast with the original contraction method, a term accounting for the mean field of all the modes not in contraction J, called spectator modes, is added to this partial Hamiltonian:
The VMFCI step consists in performing a vibrational configuration interaction [5] (Galerkin method) for each mean field Hamiltonians, Eq.(3), in the product basis sets, Eq.(2).
Thereby, we obtain new basis sets of dimension D made of eigenvectors of the mean field Hamiltonians:
,
to construct the product basis sets of the next VMFCI step. Their associated eigenvalues can be used to truncate the product basis sets according to energy criteria: either the individual component basis functions are selected only if their associated eigenvalue is less than a given threshold, or a product function is selected only if the sum of its component basis function eigenvalues is less than a given threshold. Both criteria can be applied together.
Notes
- In CONVIV, it is possible to control the number of eigenpairs calculated by the diagonalizer when
solving the mean field Hamiltonian eigenvalue problem. So, in fact, the dimension of the new basis set for contraction J, can be less than D
.