Copyright (C) 2013-2015 Sebastian Wouters sebastianwouters@gmail.com

This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.

CheMPS2 is a scientific library which contains a spin-adapted implementation of the density matrix renormalization group (DMRG) for ab initio quantum chemistry. This method allows to obtain numerical accuracy in active spaces beyond the capabilities of full configuration interaction (FCI):

• up to 40 electrons in 40 orbitals for general active spaces
• up to 100 electrons in 100 orbitals for one-dimensional active spaces, such as the pi-system of all-trans polyenes

In addition, DMRG allows to obtain the 2-RDM of the active space efficiently. The method is therefore ideal to replace the FCI solver in the complete active space configuration interaction (CASCI) and complete active space self consistent field (CASSCF) methods, when the active space sizes become prohibitively expensive for FCI. The corresponding methods are called DMRG-CI and DMRG-SCF, respectively. Because DMRG can handle the abovementioned active space sizes, it allows to obtain FCI energies for small systems such as dimers, while for larger systems it is ideal to treat the static/strong correlation in a large active space.

The design philosophy for CheMPS2 is to be a lightweight, efficient, and stable library. For an input Hamiltonian and targeted symmetry sector, the library performs successive DMRG sweeps according to a user-defined convergence scheme. As output, the library returns the minimal encountered energy as well as the 2-RDM of the active space. With the latter, various molecular properties can be calculated, as well as the gradient and Hessian for orbital rotations or nuclear displacements. In addition, several correlation functions can be obtained to investigate the electronic structure in the active space.

To gain a better understanding of how to perform DMRG calculations for large active spaces, you are encouraged to read the user manual and the three papers listed below.

Incorporation of the library into other codes is very simple due a minimal API, as well as a python interface. Direct usage of the library is illustrated in the c++ and python tests, which can be easily adapted to your needs. The interfaces to psi4 and pyscf are described in the user manual.

In the future, the parallelization of CheMPS2 for shared memory architectures will be extended to a hybrid scheme with both shared (OpenMP) and distributed (MPI) memory parallelization. In addition, the calculation of the 2-RDM of the active space will be extended to the 3-RDM.

To acknowledge CheMPS2, please cite

1. S. Wouters, W. Poelmans, P.W. Ayers and D. Van Neck, CheMPS2: a free open-source spin-adapted implementation of the density matrix renormalization group for ab initio quantum chemistry, Computer Physics Communications 185 (6), 1501-1514 (2014), doi:10.1016/j.cpc.2014.01.019

    @article{CheMPS2cite1,
        author = {Sebastian Wouters and Ward Poelmans and Paul W. Ayers and Dimitri {Van Neck}},
        title = {CheMPS2: a free open-source spin-adapted implementation of the density matrix renormalization group for ab initio quantum chemistry},
        journal = {Computer Physics Communications},
        year = {2014},
        volume = {185},
        number = {6},
        pages = {1501-1514},
        doi = {10.1016/j.cpc.2014.01.019}
    }
   

2. S. Wouters and D. Van Neck, The density matrix renormalization group for ab initio quantum chemistry, European Physical Journal D 68 (9), 272 (2014), doi:10.1140/epjd/e2014-50500-1

    @article{CheMPS2cite2,
        author = {Sebastian Wouters and Dimitri {Van Neck}},
        title = {The density matrix renormalization group for ab initio quantum chemistry},
        journal = {European Physical Journal D},
        year = {2014},
        volume = {68},
        number = {9},
        pages = {272},
        doi = {10.1140/epjd/e2014-50500-1}
    }
   

3. S. Wouters, T. Bogaerts, P. Van Der Voort, V. Van Speybroeck and D. Van Neck, Communication: DMRG-SCF study of the singlet, triplet, and quintet states of oxo-Mn(Salen), Journal of Chemical Physics 140 (24), 241103 (2014), doi:10.1063/1.4885815

    @article{CheMPS2cite3,
        author = {Sebastian Wouters and Thomas Bogaerts and Pascal {Van Der Voort} and Veronique {Van Speybroeck} and Dimitri {Van Neck}},
        title = {Communication: DMRG-SCF study of the singlet, triplet, and quintet states of oxo-Mn(Salen)},
        journal = {Journal of Chemical Physics},
        year = {2014},
        volume = {140},
        number = {24},
        pages = {241103},
        doi = {10.1063/1.4885815}
    }
   

CheMPS2 can be built with CMake and depends on

• BLAS
• LAPACK
• GSL (GNU Scientific library)
• HDF5 (Hierarchical Data Format Release 5)

It is parallelized for shared memory architectures with the Open Multi-Processing (OpenMP) API. In the future, the parallelization of CheMPS2 for shared memory architectures will be extended to a hybrid scheme with both shared (OpenMP) and distributed (MPI) memory parallelization.

It is advised to clone the CheMPS2 git repository from github. In your terminal, run:

> cd /sourcefolder
> git clone 'https://github.com/sebwouters/chemps2'
> cd chemps2

That way, future updates and bug fixes can be easily pulled in:

> cd /sourcefolder/chemps2
> git pull

The files

/sourcefolder/chemps2/CMakeLists.txt
/sourcefolder/chemps2/CheMPS2/CMakeLists.txt
/sourcefolder/chemps2/tests/CMakeLists.txt
/sourcefolder/chemps2/sphinx/CMakeLists.txt

provide a minimal compilation. In your terminal, run:

> cd /sourcefolder/chemps2
> mkdir build
> cd build

CMake generates makefiles based on the user’s specifications:

> CXX=option1 cmake .. -DMKL=option2 -DCMAKE_INSTALL_PREFIX=option3 -DBUILD_DOXYGEN=option4 -DBUILD_SPHINX=option5

1. Option1 is the c++ compiler; typically g++, icpc, or clang++ on Linux.
2. Option2 can be ON or OFF and is used to switch on the intel math kernel library.
3. Option3 is the prefix of the installation directory; typically /usr or /usr/local on Linux. On my computer, libchemps2 is then installed in /prefix/lib/x86_64-linux-gnu and the headers in /prefix/include/chemps2.
4. Option4 can be ON or OFF and is used to switch on the possibility to compile the doxygen documentation.
5. Option5 can be ON or OFF and is used to switch on the possibility to compile the user manual with sphinx.

If one or more of the required libraries are not found, use the command

> CMAKE_INCLUDE_PATH=option6 CMAKE_LIBRARY_PATH=option7 CXX=option1 cmake .. -DMKL=option2 -DCMAKE_INSTALL_PREFIX=option3 -DBUILD_DOXYGEN=option4 -DBUILD_SPHINX=option5

instead, where option6 and option7 are respectively the missing colon-separated include and library paths:

CMAKE_INCLUDE_PATH=/my_libs/lib1/include:/my_libs/lib2/include
CMAKE_LIBRARY_PATH=/my_libs/lib1/lib:/my_libs/lib2/lib

To compile, run:

> make

To install, run:

> make install

For non-standard installation directories, please remember to append the library path to LD_LIBRARY_PATH in your .bashrc.

To test libchemps2, run:

> cd /sourcefolder/chemps2/build
> make test

The tests only require a very limited amount of memory (order 10-120 MB).

PyCheMPS2, a python interface to libchemps2, can be built with Cython. The installation is independent of CMake and assumes that you have installed the CheMPS2 library with make install. For non-standard installation directories of CheMPS2, please remember to append the library path to LD_LIBRARY_PATH in your .bashrc. In addition, the include path should be appended to CPATH:

> export CPATH=${CPATH}:option3/include

where option3 is the option provided to CMake with -DCMAKE_INSTALL_PREFIX=option3 above. Then the python wrapper can be installed with:

> cd /sourcefolder/chemps2/PyCheMPS2
> python setup.py build_ext -L ${LD_LIBRARY_PATH}
> python setup.py install --prefix=option3

On my machine, the python wrapper is installed to the folder option3/lib/python2.7/site-packages, but the folder lib and the distribution of python can vary.

Compilation of PyCheMPS2 occurs by linking to the c++ library in the installation directory. The installation of PyCheMPS2 will fail if that library is not properly installed. If you have pulled a newer version of CheMPS2, please remember to reinstall the c++ library first, before reinstalling PyCheMPS2!

To test PyCheMPS2 (remember that the python site-packages folder can vary), run:

> cd /sourcefolder/chemps2/PyCheMPS2/tests
> export PYTHONPATH=${PYTHONPATH}:option3/lib/python2.7/site-packages
> python test1.py
> python test2.py
> python test3.py
> python test4.py
> python test5.py
> python test6.py
> python test7.py
> python test8.py
> python test9.py

If you compiled the c++ library with -DMKL=ON, you might get the error

Intel MKL FATAL ERROR: Cannot load libmkl_avx.so or libmkl_def.so.

This issue of using Intel’s MKL inside python is known and reported. To get the python tests to run, you can set the variable LD_PRELOAD in order to preload libmkl_rt.so. On my system, this is done with

> export LD_PRELOAD=/opt/intel/mkl/lib/intel64/libmkl_rt.so

The python tests do exactly the same thing as the c++ tests above, and illustrate the usage of the python interface to libchemps2. The tests should end with a line stating whether or not they succeeded. They only require a very limited amount of memory (order 10-120 MB).

To build the doxygen manual, the BUILD_DOXYGEN flag should have been on: -DBUILD_DOXYGEN=ON. In your terminal, run:

> cd /sourcefolder/chemps2/build
> make doc
> cd LaTeX-documents
> make
> evince refman.pdf &
> cd ../html
> firefox index.html &

To build the sphinx user manual, the BUILD_SPHINX flag should have been on: -DBUILD_SPHINX=ON. In your terminal, run:

> cd /sourcefolder/chemps2/build
> make sphinx
> cd sphinx/html
> firefox index.html &

For information on how to perform DMRG and DMRG-SCF calculations with CheMPS2, please consult the

1. publications
2. user manual
3. doxygen html output

The user manual contains elaborate information on

• the DMRG and DMRG-SCF algorithms
• the symmetries which are exploited in CheMPS2
• how to generate matrix elements with plugins to psi4
• how to perform DMRG and DMRG-SCF calculations
• the interfaces of CheMPS2 to psi4 and pyscf

CheMPS2/CASSCF.cpp contains the constructor and destructor of the CASSCF class, as well as the functions which allow to update, get or set its variables.

CheMPS2/CASSCFdebug.cpp contains two functions: one for calculating the ROHF energy; and one for fetching FCI coefficients to determine the point group symmetry of certain electronic states of the carbon dimer.

CheMPS2/CASSCFnewtonraphson.cpp contains all DMRG-SCF functions which are specific to the augmented Hessian Newton-Raphson update scheme, including the functions to calculate the gradient and Hessian.

CheMPS2/ConvergenceScheme.cpp contains all functions of the ConvergenceScheme class, which contains the instructions for the subsequent DMRG sweeps.

CheMPS2/Correlations.cpp contains all the functionality to calculate the spin, density, and spin-flip correlation functions as well as the two-orbital mutual information.

CheMPS2/Davidson.cpp is an implementation of Davidson's algorithm.

CheMPS2/DIIS.cpp contains a DIIS convergence speed-up for DMRG-SCF.

CheMPS2/DMRG.cpp contains the constructor and destructor of the DMRG class, as well as the top-level sweep functions.

CheMPS2/DMRGmpsio.cpp contains the store and load functions for the DMRG checkpoint file (the MPS).

CheMPS2/DMRGoperators.cpp contains all functions related to the DMRG renormalized operators: saving to disk, loading from disk, and updating.

CheMPS2/DMRGSCFindices.cpp contains the index conversions for the DMRG-SCF algorithm.

CheMPS2/DMRGSCFintegrals.cpp is a container class for two-body matrix elements with at most two virtual indices.

CheMPS2/DMRGSCFmatrix.cpp is a container class for orbital matrices which are blockdiagonal in the irreps.

CheMPS2/DMRGSCFoptions.cpp is a container class to pass the DMRGSCF options to the augmented Hessian Newton-Raphson routine.

CheMPS2/DMRGSCFunitary.cpp contains the storage and handling of the unitary matrix and its nonredundant skew-symmetric parametrization required for the DMRG-SCF algorithm.

CheMPS2/DMRGSCFVmatRotations.cpp performs the two-body matrix element rotations for the DMRGSCF and Edmiston-Ruedenberg classes.

CheMPS2/DMRGSCFwtilde.cpp is a container class to store a tensor which is required to compute the DMRG-SCF Hessian.

CheMPS2/DMRGtechnics.cpp contains the functions related to the 2-RDM and the excited-state calculations.

CheMPS2/EdmistonRuedenberg.cpp contains an orbital localization function based on the Edmiston-Ruedenberg cost function and an augmented Hessian Newton-Raphson optimizer.

CheMPS2/FCI.cpp contains a full configuration interaction solver based on Davidson's algorithm. It also contains the functionality to calculate Green's functions.

CheMPS2/FourIndex.cpp contains all functions of the FourIndex container class for the two-body matrix elements.

CheMPS2/Hamiltonian.cpp contains all functions of the Hamiltonian class, including functions to get or set specific variables, as well as the save and load functions to store the Hamiltonian on disk.

CheMPS2/Heff.cpp contains the constructor and destructor of the Heff class, as well as the top-level functions to perform the effective Hamiltonian times guess-vector multiplication and Davidson's algorithm.

CheMPS2/HeffDiagonal.cpp contains the functions to calculate the diagonal elements of the effective Hamiltonian. These are required as preconditioner in Davidson's algorithm.

CheMPS2/HeffDiagrams1.cpp contains a subset of functions to perform the effective Hamiltonian times guess-vector multiplication.

CheMPS2/HeffDiagrams2.cpp contains a subset of functions to perform the effective Hamiltonian times guess-vector multiplication.

CheMPS2/HeffDiagrams3.cpp contains a subset of functions to perform the effective Hamiltonian times guess-vector multiplication.

CheMPS2/HeffDiagrams4.cpp contains a subset of functions to perform the effective Hamiltonian times guess-vector multiplication.

CheMPS2/HeffDiagrams5.cpp contains a subset of functions to perform the effective Hamiltonian times guess-vector multiplication.

CheMPS2/Initialize.cpp allows to set the seed of the random number generator and cout.precision (added for PyCheMPS2).

CheMPS2/Irreps.cpp contains the psi4 symmetry conventions.

CheMPS2/PrintLicense.cpp contains a function which prints the license disclaimer.

CheMPS2/Problem.cpp contains all Problem class functions. This wrapper class allows to set the desired symmetry sector for the DMRG algorithm.

CheMPS2/Sobject.cpp contains all Sobject class functions. This class constructs, stores, and decomposes the reduced two-site object.

CheMPS2/SyBookkeeper.cpp contains all SyBookkeeper functions. This class keeps track of the FCI and DMRG virtual dimensions of all symmetry sectors at all boundaries.

CheMPS2/TensorA.cpp contains the TensorA functions. This class stores and handles the complementary reduced spin-0 S0-tensors.

CheMPS2/TensorB.cpp contains the TensorB functions. This class stores and handles the complementary reduced spin-1 S1-tensors.

CheMPS2/TensorC.cpp contains the TensorC functions. This class stores and handles the complementary reduced spin-0 F0-tensors.

CheMPS2/TensorD.cpp contains the TensorD functions. This class stores and handles the complementary reduced spin-1 F1-tensors.

CheMPS2/TensorDiag.cpp contains the TensorDiag functions. This storage class handles reduced tensors which are diagonal in the symmetry sectors.

CheMPS2/TensorF0Cbase.cpp contains the common storage and handling functions of the TensorF0 and TensorC classes.

CheMPS2/TensorF0.cpp contains all TensorF0 functions. This class stores and handles the reduced spin-0 part of two sandwiched second quantized operators, of which the particle symmetry sectors are equal.

CheMPS2/TensorF1.cpp contains all TensorF1 functions. This class stores and handles the reduced spin-1 part of two sandwiched second quantized operators, of which the particle symmetry sectors are equal.

CheMPS2/TensorF1Dbase.cpp contains the common storage and handling functions of the TensorF1 and TensorD classes.

CheMPS2/TensorGYZ.cpp contains the contruct and update functions for the G-, Y-, and Z-tensors. They are required for the two-orbital mutual information.

CheMPS2/TensorK.cpp contains the contruct and update functions for the K-tensor. It is required for the two-orbital mutual information.

CheMPS2/TensorL.cpp contains all TensorL functions. This class stores and handles the reduced spin-1/2 part of a single sandwiched second quantized operator.

CheMPS2/TensorM.cpp contains the contruct and update functions for the M-tensor. It is required for the two-orbital mutual information.

CheMPS2/TensorO.cpp implements the storage and handling of the partial terms which are required to calculate the overlap between two MPSs.

CheMPS2/TensorQ.cpp contains all TensorQ functions. This class stores and handles the complementary reduced spin-1/2 part of three sandwiched second quantized operators.

CheMPS2/TensorS0Abase.cpp contains the common storage and handling functions of the TensorS0 and TensorA classes.

CheMPS2/TensorS0.cpp contains all TensorS0 functions. This class stores and handles the reduced spin-0 part of two sandwiched second quantized operators, of which the particle symmetry sectors differ by 2.

CheMPS2/TensorS1Bbase.cpp contains the common storage and handling functions of the TensorS1 and TensorB classes.

CheMPS2/TensorS1.cpp contains all TensorS1 functions. This class stores and handles the reduced spin-1 part of two sandwiched second quantized operators, of which the particle symmetry sectors differ by 2.

CheMPS2/TensorSwap.cpp contains the common storage and handling functions of the TensorL and TensorQ classes.

CheMPS2/TensorT.cpp contains all TensorT functions. This class stores and handles the reduced part of an MPS site-tensor.

CheMPS2/TensorX.cpp contains all TensorX functions. This class stores and handles the complementary reduced spin-0 part of four sandwiched second quantized operators, which is of course diagonal in the symmetry sectors.

CheMPS2/TwoDM.cpp contains all functions to calculate the 2-RDM from the DMRG-optimized MPS.

CheMPS2/TwoDMstorage.cpp contains all functions to store the 2-RDM from the DMRG-optimized MPS.

CheMPS2/TwoIndex.cpp contains all functions of the TwoIndex container class for the one-body matrix elements.

CheMPS2/include/CASSCF.h contains the definitions of the CASSCF class.

CheMPS2/include/ConvergenceScheme.h contains the definitions of the ConvergenceScheme class.

CheMPS2/include/Correlations.h contains the definitions of the Correlations class.

CheMPS2/include/Davidson.h contains the definitions of the Davidson class.

CheMPS2/include/DIIS.h contains the definitions of the DIIS class.

CheMPS2/include/DMRG.h contains the definitions of the DMRG class.

CheMPS2/include/DMRGSCFindices.h contains the definitions of the DMRGSCFindices class.

CheMPS2/include/DMRGSCFintegrals.h contains the definitions of the DMRGSCFintegrals class.

CheMPS2/include/DMRGSCFmatrix.h contains the definitions of the DMRGSCFmatrix class.

CheMPS2/include/DMRGSCFoptions.h contains the definitions of the DMRGSCFoptions container class.

CheMPS2/include/DMRGSCFunitary.h contains the definitions of the DMRGSCFunitary class.

CheMPS2/include/DMRGSCFVmatRotations.h contains the definitions of the DMRGSCFVmatRotations class.

CheMPS2/include/DMRGSCFwtilde.h contains the definitions of the DMRGSCFwtilde class.

CheMPS2/include/EdmistonRuedenberg.h contains the definitions of the EdmistonRuedenberg class.

CheMPS2/include/FCI.h contains the definitions of the FCI class.

CheMPS2/include/FourIndex.h contains the definitions of the FourIndex class.

CheMPS2/include/Gsl.h contains the definitions of the external GSL routines.

CheMPS2/include/Hamiltonian.h contains the definitions of the Hamiltonian class.

CheMPS2/include/Heff.h contains the definitions of the Heff class.

CheMPS2/include/Initialize.h contains the definitions of the Initialize class.

CheMPS2/include/Irreps.h contains the definitions of the Irrep class.

CheMPS2/include/Lapack.h contains the definitions of the external BLAS and LAPACK routines.

CheMPS2/include/MyHDF5.h forces the use of the HDF5 1.8 API, e.g. H5Gcreate2 instead of H5Gcreate1, a known issue in Ubuntu 12.04.

CheMPS2/include/Options.h contains all the options of the CheMPS2 namespace. Here the checkpoint storage names and folders can be set, as well as parameters related to memory usage and convergence.

CheMPS2/include/Problem.h contains the definitions of the Problem class.

CheMPS2/include/Sobject.h contains the definitions of the Sobject class.

CheMPS2/include/SyBookkeeper.h contains the definitions of the SyBookkeeper class.

CheMPS2/include/TensorA.h contains the definitions of the TensorA class.

CheMPS2/include/TensorB.h contains the definitions of the TensorB class.

CheMPS2/include/TensorC.h contains the definitions of the TensorC class.

CheMPS2/include/TensorD.h contains the definitions of the TensorD class.

CheMPS2/include/TensorDiag.h contains the definitions of the TensorDiag class.

CheMPS2/include/TensorF0Cbase.h contains the definitions of the TensorF0Cbase class.

CheMPS2/include/TensorF0.h contains the definitions of the TensorF0 class.

CheMPS2/include/TensorF1Dbase.h contains the definitions of the TensorF1Dbase class.

CheMPS2/include/TensorF1.h contains the definitions of the TensorF1 class.

CheMPS2/include/TensorGYZ.h contains the definitions of the TensorGYZ class.

CheMPS2/include/Tensor.h contains the definitions of the virtual Tensor class.

CheMPS2/include/TensorK.h contains the definitions of the TensorK class.

CheMPS2/include/TensorL.h contains the definitions of the TensorL class.

CheMPS2/include/TensorM.h contains the definitions of the TensorM class.

CheMPS2/include/TensorO.h contains the definitions of the TensorO class.

CheMPS2/include/TensorQ.h contains the definitions of the TensorQ class.

CheMPS2/include/TensorS0Abase.h contains the definitions of the TensorS0Abase class.

CheMPS2/include/TensorS0.h contains the definitions of the TensorS0 class.

CheMPS2/include/TensorS1Bbase.h contains the definitions of the TensorS1Bbase class.

CheMPS2/include/TensorS1.h contains the definitions of the TensorS1 class.

CheMPS2/include/TensorSwap.h contains the definitions of the TensorSwap class.

CheMPS2/include/TensorT.h contains the definitions of the TensorT class.

CheMPS2/include/TensorX.h contains the definitions of the TensorX class.

CheMPS2/include/TwoDM.h contains the definitions of the TwoDM class.

CheMPS2/include/TwoDMstorage.h contains the definitions of the TwoDMstorage class.

CheMPS2/include/TwoIndex.h contains the definitions of the TwoIndex class.

Please note that these files are documented with doxygen comments. The installation section discusses how a manual can be generated from these comments.

tests/test1.cpp.in contains several DMRG ground state calculations in different symmetry sectors for the N2 molecule (d2h symmetry) in the minimal STO-3G basis set.

tests/test2.cpp.in contains a ground state DMRG calculation of the ^1A1 state of H2O (c2v symmetry) in the 6-31G basis set.

tests/test3.cpp.in contains a ground state DMRG calculation of the ^1A1 state of CH4 (c2v symmetry) in the STO-3G basis set.

tests/test4.cpp.in contains a ground state DMRG calculation of the ^6A state of a linear Hubbard chain (forced c1 symmetry) with 10 sites and open boundary conditions, containing 9 fermions (just below half-filling).

tests/test5.cpp.in contains an excited state DMRG calculation in the ^1Ag symmetry sector of N2 (d2h symmetry) in the minimal STO-3G basis set. The ground and two lowest excited states are determined.

tests/test6.cpp.in contains a state-averaged DMRG-SCF calculation of the first excited state of the ^1Ag sector of O2 (d2h symmetry) in the CC-pVDZ basis set. The two 1s core orbitals are kept frozen, and two Ag, B2g, B3g, B1u, B2u, and B3u orbitals are chosen as active space. A significant speedup is obtained with DIIS.

tests/test7.cpp.in reads in tests/matrixelements/O2.CCPVDZ.FCIDUMP, stores these matrix elements to disk, reads them back in from disk, and compares the two versions.

tests/test8.cpp.in contains a DMRG-SCF ground state calculation of the ^1Ag state of N2 (d2h symmetry) in the CC-pVDZ basis set. The two 1s core orbitals are kept frozen. The next two Ag and B1u orbitals (sigma bonding and antibonding), as well as one B2g, B3g, B2u, and B3u orbital (pi bonding and antibonding) are chosen as active space. A significant speedup is obtained with DIIS. This test is smaller than test6, and is included for debugging with valgrind.

tests/test9.cpp.in is a copy of tests/test8.cpp.in, with a slightly larger active space and which works with ordered localized orbitals instead of natural orbitals. The localization occurs by means of Edmiston-Ruedenberg, and the ordering based on the Fiedler vector with the exchange matrix as cost function.

tests/matrixelements/CH4.STO3G.FCIDUMP contains the matrix elements for test3.

tests/matrixelements/H2O.631G.FCIDUMP contains the matrix elements for test2.

tests/matrixelements/N2.STO3G.FCIDUMP contains the matrix elements for test1 and test5.

tests/matrixelements/O2.CCPVDZ.FCIDUMP contains the matrix elements for test6 and test7.

tests/matrixelements/N2.CCPVDZ.FCIDUMP contains the matrix elements for test8 and test9.

PyCheMPS2/test1.py.in is the python version of tests/test1.cpp.in

PyCheMPS2/test2.py.in is the python version of tests/test2.cpp.in

PyCheMPS2/test3.py.in is the python version of tests/test3.cpp.in

PyCheMPS2/test4.py.in is the python version of tests/test4.cpp.in

PyCheMPS2/test5.py.in is the python version of tests/test5.cpp.in

PyCheMPS2/test6.py.in is the python version of tests/test6.cpp.in

PyCheMPS2/test7.py.in is the python version of tests/test7.cpp.in

PyCheMPS2/test8.py.in is the python version of tests/test8.cpp.in

PyCheMPS2/test9.py.in is the python version of tests/test9.cpp.in

These test files illustrate how to use libchemps2. The tests only require a very limited amount of memory (order 10-120 MB).