# 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 *n _{vib}* vibrational degrees of freedom into

*n*subsets :

_{p}

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*always corresponds to product of ground state functions._{}=(0, ... ,0)

- 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

_{}*H*involving the degrees in components

_{vib}*I*of

_{jl}*J*:

_{}

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

_{}