QuEmb

QuEmb is a robust framework designed to implement the Bootstrap Embedding (BE) method, efficiently treating electron correlation in molecules, surfaces, and solids. This repository contains the Python implementation of the BE methods, including periodic bootstrap embedding. The code leverages PySCF library for quantum chemistry calculations and utlizes Python's multiprocessing module to enable parallel computations in high-performance computing environments.

QuEmb includes two libraries: molbe and kbe. The molbe library implements BE for molecules and supramolecular complexes, while the kbe library is designed to handle periodic systems such as surfaces and solids using periodic BE.

Features

• Fragment-based quantum embedding: Utilizes flexible system partioning with overlapping regions to improve quantum embedding techniques.
• Periodic Bootstrap Embedding: Extends BE method to treat periodic systems (1D & 2D systems) using reciprocal space sums.
• High accuracy and efficiency: Capable of recovering ~99.9% of electron correlation energy.
• Parallel computing: Employ's Python multiprocessing module to perform parallel computations across multiple processors.

Installation

Prerequisites

• Python >=3.10
• PySCF library
• Numpy
• Scipy
• libDMET (required for periodic BE)
• Wannier90^## (to use Wannier functions)

^## Wannier90 code is optional and only necessary to use Wannier functions in periodic code.

The required dependencies, with the exception of the optional Wannier90, are automatically installed by pip.

Installation

One can just pip install directly from the Github repository:

pip install git+https://https://github.com/troyvvgroup/quemb

Alternatively one can manually clone and install as in:

git clone https://https://github.com/troyvvgroup/quemb
cd quemb
pip install .

Basic Usage

# Molecular
from quemb.molbe import fragpart
from quemb.molbe import BE

# Periodic
#from quemb.kbe import fragpart
#from quemb.kbe import BE

# Perform pyscf HF/KHF calculations
# get mol: pyscf.gto.M or pyscf.pbc.gto.Cell
# get mf: pyscf.scf.RHF or pyscf.pbc.KRHF

# Define fragments
myFrag = fragpart(be_type='be2', mol=mol)

# Initialize BE
mybe = BE(mf, myFrag)

# Perform density matching in BE
mybe.optimize(solver='CCSD')

See documentation and quemb/example for more details.

Documentation

Comprehensive documentation for QuEmb is available at quemb/docs. The documentation provides detailed infomation on installation, usage, API reference, and examples. To build the documentation locally, simply navigate to docs and build using make html or make latexpdf.

Alternatively, you can view the latest documentation online here.

References

The methods implemented in this code are described in details in the following papers:

• OR Meitei, T Van Voorhis, Periodic bootstrap embedding, JCTC 19 3123 2023
• OR Meitei, T Van Voorhis, Electron correlation in 2D periodic systems, arXiv:2308.06185
• HZ Ye, HK Tran, T Van Voorhis, Bootstrap embedding for large molecular systems, JCTC 16 5035 2020