qat.fermion.chemistry.wrapper.MolecularHamiltonian

class qat.fermion.chemistry.wrapper.MolecularHamiltonian(one_body_integrals: ndarray, two_body_integrals: ndarray, constant_coeff: ndarray)

MolecularHamiltonian helper class. It represents the electronic-structure Hamiltonian defined using one- and two-body integrals.

This electronic-structure Hamiltonian is defined by:

$$H=\sum _{uv\sigma }{I}_{uv}{c}_{u\sigma }^{†}{c}_{v\sigma }+\frac{1}{2}\sum _{uvwx}\sum _{\sigma {\sigma }^{\prime }}{I}_{uvwx}{c}_{u\sigma }^{†}{c}_{v{\sigma }^{\prime }}^{†}{c}_{w{\sigma }^{\prime }}{c}_{x\sigma }+r\mathbb{I}$$

with $r$ the core repulsion constant, and with $I_{uv}$ and $I_{uvwx}$ the one- and two-body integrals defined by:

$$\begin{array}{r}\begin{array}{c}{I}_{uv}=\int dr{\varphi }_u^{\ast }(r)h_1[{\varphi }_v(r)]\\ I_{uvwx}=\int drd{r}^{\prime }{\varphi }_u^{\ast }(r){\varphi }_v^{\ast }(r^{\prime })v[{\varphi }_w(r){\varphi }_x(r^{\prime })]\end{array}\end{array}$$

Here, $\{{\varphi }_i(r){\}}_{i=0...N-1}$ is the single-particle basis, with $N$ the size, which depends on the basis chosen. $h_1=h_{kin}+h_{pot}$ is the one-body Hamiltonian, and $v$ the Coulomb operator.

Note

This electronic-structure Hamiltonian definition is different than the one used in ElectronicStructureHamiltonian.

Parameters:

• one_body_integrals (np.ndarray) – One-body integral $I_{uv}$.

• two_body_integrals (np.ndarray) – Two-body integral $I_{uvwx}$.

• constant_coeff (np.ndarray) – Constant coefficient $r$ (core repulsion).

nqbits

The total number of qubits.

Type:

int

one_body_integrals

One-body integral $I_{uv}$.

Type:

np.ndarray

two_body_integrals

Two-body integral $I_{uvwx}$.

Type:

np.ndarray

constant_coeff

Constant coefficient $r$ (core repulsion).

Type:

np.ndarray

Example

import numpy as np
from qat.fermion.chemistry import MolecularHamiltonian

# Initialize random one- and two-body integrals, and a constant
one_body_integral = np.random.randn(2, 2)
two_body_integral = np.random.randn(2, 2, 2, 2)
constant = np.random.rand()

# Define the MolecularHamiltonian
mol_h = MolecularHamiltonian(one_body_integral, two_body_integral, constant)

print(mol_h)

 MolecularHamiltonian(
 - constant_coeff : 0.9796315264586001
 - integrals shape
    * one_body_integrals : (2, 2)
    * two_body_integrals : (2, 2, 2, 2)
)

get_electronic_hamiltonian() → ElectronicStructureHamiltonian

Converts the MolecularHamiltonian to an ElectronicStructureHamiltonian. To do so, it converts from $I_{uv},I_{uvwx}$ to $h_{pq},h_{pqrs}$, with

$$\begin{array}{r}\begin{array}{c}{h}_{u\sigma ,v{\sigma }^{\prime }}=I_{u,v}{\delta }_{\sigma ,{\sigma }^{\prime }}\\ h_{u{\sigma }_1,v{\sigma }_2,w{\sigma }_2^{\prime },x{\sigma }_1^{\prime }}=I_{uvwx}((1-{\delta }_{\sigma ,{\sigma }^{\prime }})+{\delta }_{\sigma ,{\sigma }^{\prime }}(1-{\delta }_{u,v})(1-{\delta }_{w,x}))\end{array}\end{array}$$

and where the one- and two-body integrals are defined as:

$$I_{uv}\equiv (u|h|v)=\int \mathrm{d}r{\varphi }_u^{\ast }(r)T{\varphi }_v(r)$$

$$I_{uvwx}\equiv (ux|vw)=\iint \mathrm{d}{r}_1\mathrm{d}{r}_2{\varphi }_u^{\ast }(r_1){\varphi }_x(r_1)v(r_{12}){\varphi }_v^{\ast }(r_2){\varphi }_w(r_2)$$

with $T$ (resp. $v$) the one- (resp. two-) body potentials, and ${\varphi }_u(r)$ is the molecular orbital wavefunction.

The $h$ integrals are used to construct hamiltonians of the ElectronicStructureHamiltonian type.

Returns:

ElectronicStructureHamiltonian Electronic structure hamiltonian.

select_active_space(noons: List[float], n_electrons: int, threshold_1: float | None = 0.02, threshold_2: float | None = 0.001) → Tuple[MolecularHamiltonian, List[int], List[int]]

Selects the right active space and freezes core electrons according to their NOONs $n_i$.

This function is an implementation of the Complete Active Space (CAS) approach. It divides orbital space into sets of active and inactive orbitals, the occupation number of the latter remaining unchanged during the computation.

The active space indices are defined as:

$$\mathcal{A}=\{i,n_i\in [{\epsilon }_2,2-{\epsilon }_1[\}\cup \{i,n_i\ge 2-{\epsilon }_1,2(i+1)\ge N_e\}$$

The inactive occupied orbitals are defined as:

$$\mathcal{O}=\{i,n_i\ge 2-{\epsilon }_1,2(i+1)<N_e\}$$

The restriction of the one- and two-body integrals (and update of the core energy) is then carried out according to:

$$\mathrm{\forall }u,v\in \mathcal{A},I_{uv}^{(a)}=I_{uv}+\sum _{i\in \mathcal{O}}2I_{i,u,v,i}-I_{i,u,i,v}$$

$$\mathrm{\forall }u,v,w,x\in \mathcal{A},I_{uvwx}^{(a)}=I_{uvwx}$$

$$E_{\mathrm{core}}^{(a)}=E_{\mathrm{core}}+\sum _{i\in \mathcal{O}}{I}_{ii}+\sum _{ij\in \mathcal{O}}2I_{ijji}-I_{ijij}$$

Parameters:

• noons (List[float]) – the natural-orbital occupation numbers $n_i$, sorted in descending order (from high occupations to low occupations)

• n_electrons (int) – The number of electrons $N_e$.

• threshold_1 (Optional[float]) – The upper threshold ${\epsilon }_1$ on the NOON of an active orbital.

• threshold_2 (Optional[float]) – The lower threshold ${\epsilon }_2$ on the NOON of an active orbital.

Returns:

• the molecular Hamiltonian in active space $H^{(a)}$

• the list of indices corresponding to the active orbitals, $\mathcal{A}$

• the list of indices corresponding to the occupied orbitals, $\mathcal{O}$

Return type:

Tuple[MolecularHamiltonian, List[int], List[int]]

transform_basis(transformation_matrix: ndarray) → MolecularHamiltonian

Change one and two body integrals (indices p, q…) to new basis (indices i, j…) using transformation U such that

$${\hat{c}}_i=\sum _q{U}_{qi}{c}_q$$

i.e

$$\begin{array}{r}\begin{array}{c}{\hat{I}}_{ij}=\sum _{pq}{U}_{pi}{I}_{pq}{U}_{jq}^{†}\\ {\hat{I}}_{ijkl}=\sum _{pqrs}{U}_{pi}{U}_{qj}{I}_{pqrs}{U}_{kr}^{†}{U}_{ls}^{†}\end{array}\end{array}$$

Parameters:

transformation_matrix (np.array) – transformation matrix $U$

Returns:

MolecularHamiltonian updated to the new basis.

Return type:

molecular_hamiltonian (MolecularHamiltonian)