qat.fermion.chemistry.wrapper.MoleculeInfo

class qat.fermion.chemistry.wrapper.MoleculeInfo(hamiltonian: MolecularHamiltonian, n_electrons: int, noons: ndarray | List[float], orbital_energies: ndarray)

MoleculeInfo helper class. This class is a even higher level version of the MolecularHamiltonian.

Parameters:

• hamiltonian (MolecularHamiltonian) – The MolecularHamiltonian of the studied molecule.

• n_electrons (int) – Number of electrons.

• noons (Union[np.ndarray, List[float]]) – Natural orbital occupation number.

• orbital_energies (np.ndarray) – Orbital energies.

nqbits

The total number of qubits.

Type:

int

one_body_integrals

One-body integrals \(I_{uv}\).

Type:

np.ndarray

two_body_integrals

Two-body integrals \(I_{uvwx}\).

Type:

np.ndarray

constant_coeff

Constant coefficient \(r\) (core repulsion).

Type:

np.ndarray

hamiltonian

The MolecularHamiltonian of the studied molecule.

Type:

MolecularHamiltonian

n_electrons

Number of electrons.

Type:

int

noons

Natural orbital occupation number.

Type:

Union[np.ndarray, List[float]]

orbital_energies

Orbital energies.

Type:

np.ndarray

Example

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

# For illustration purpose, 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()
noons = list(np.random.randn(10))
orbital_energies = list(np.random.randn(10))

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

# Define MoleculeInfo
molecule = MoleculeInfo(
    mol_h,
    n_electrons=4,
    noons=noons,
    orbital_energies=orbital_energies
)

print(molecule)

MoleculeInfo(
 - MolecularHamiltonian(
    * constant_coeff : 0.7860417265255741
    * integrals shape
       ** one_body_integrals : (2, 2)
       ** two_body_integrals : (2, 2, 2, 2)
   )
 - n_electrons = 4
 - noons = [np.float64(2.0096740839670644), np.float64(0.09969545342290652), np.float64(1.0347201718418875), np.float64(-0.3866663821225546), np.float64(0.41887104570565087), np.float64(-0.5124012698410086), np.float64(-0.04354229203479397), np.float64(-2.1708797602858594), np.float64(1.438845092297319), np.float64(-0.40960453230934074)]
 - orbital energies = [np.float64(-0.187207465983109), np.float64(0.2450355409971492), np.float64(-0.844017434525592), np.float64(-0.8786620538867721), np.float64(0.3967040036706545), np.float64(1.7704905125702282), np.float64(-0.022987744188976873), np.float64(-1.7332509360637418), np.float64(-1.5669216101695511), np.float64(-3.1247733521961245)]
)

restrict_active_space(threshold_1: float | None = 0.02, threshold_2: float | None = 0.001)

Same method as the MolecularHamiltonian method select_active_space(), except it also modifies all the molecule parameters accordingly (NOONs, orbital energies, and number of electrons).

For more information, see select_active_space() documentation.

Parameters:

• threshold_1 (Optional[float]) – The upper threshold \(\varepsilon_1\) on the NOON of an active orbital.

• threshold_2 (Optional[float]) – The lower threshold \(\varepsilon_2\) on the NOON of an active orbital.