"""Bond topology features from ``Multiwfn``."""
from __future__ import annotations
import os
from typing import TYPE_CHECKING, Any, Dict, List, Optional, Tuple, Union
import numpy as np
from bonafide.utils.base_featurizer import BaseFeaturizer
from bonafide.utils.constants import PROGRAM_ENVIRONMENT_VARIABLES
from bonafide.utils.driver import multiwfn_driver
from bonafide.utils.helper_functions_chemistry import align_coordinates
from bonafide.utils.io_ import read_xyz_file
from bonafide.utils.multiwfn_properties import read_prop_file
if TYPE_CHECKING:
from numpy.typing import NDArray
from rdkit import Chem
[docs]
class _Multiwfn3DBondTopology(BaseFeaturizer):
"""Parent feature factory for the 3D bond topology Multiwfn features.
For details, please refer to the Multiwfn manual (http://sobereva.com/multiwfn/, last accessed
on 12.09.2025).
"""
neighbor_distance_cutoff: float
step_size: float
def __init__(self) -> None:
self.extraction_mode = "multi"
super().__init__()
[docs]
def calculate(self) -> None:
"""Calculate the Multiwfn bond topology features."""
self._run_multiwfn()
self._read_cp_prop_file()
[docs]
def _run_multiwfn(self) -> None:
"""Run Multiwfn.
Returns
-------
None
"""
# Select topology analysis
multiwfn_commands: List[Union[str, int, float]]
multiwfn_commands = [2]
# Add path generation options
multiwfn_commands.extend([-2, 2, self.step_size, 3, self.neighbor_distance_cutoff, 0])
# Search critical points
multiwfn_commands.extend([2, 3])
# Generate paths
multiwfn_commands.append(8)
# Print results to separate text file
multiwfn_commands.extend([7, 0])
# Exit program
multiwfn_commands.extend([-10, "q"])
# Set up environment variables
environment_variables = {
var: getattr(self, var, None) for var in PROGRAM_ENVIRONMENT_VARIABLES["multiwfn"]
}
# Run Multiwfn
multiwfn_driver(
cmds=multiwfn_commands,
input_file_path=str(self.electronic_struc_n),
output_file_name=f"Multiwfn3DBondTopology__{self.conformer_name}",
environment_variables=environment_variables,
namespace=self.conformer_name[::-1].split("__", 1)[-1][::-1],
)
# The atom coordinates used by Multiwfn are also needed as they can differ from the ones in
# the mol object. Therefore, they are written to a file for later retrieval.
_output_file_name = f"Multiwfn3DAtomXyzCoordinates__{self.conformer_name}"
multiwfn_driver(
cmds=[300, 7, -1, f"_xyz_from_multiwfn__{self.conformer_name}.xyz", -10, 0, "q"],
input_file_path=str(self.electronic_struc_n),
output_file_name=_output_file_name,
environment_variables=environment_variables,
namespace=self.conformer_name[::-1].split("__", 1)[-1][::-1],
)
# Directly remove the Multiwfn output as it does not contain any relevant data.
# The coordinates are written to a separate file.
if os.path.isfile(_output_file_name) is True:
os.remove(_output_file_name)
[docs]
def _read_cp_prop_file(self) -> None:
"""Read the output file from Multiwfn and write the results to the ``results`` dictionary.
Returns
-------
None
"""
# Check if the output file exists
_opath = f"CPprop__{self.conformer_name}.txt"
if os.path.isfile("CPprop.txt") is False:
self._err = (
f"Multiwfn output file '{_opath}' not found; probably the calculation "
"did not run. Check your input"
)
return
# Check if the additional atom coordinates files exists
_opath2 = f"_xyz_from_multiwfn__{self.conformer_name}.xyz"
if os.path.isfile(_opath2) is False:
self._err = (
f"Multiwfn output file '{_opath2}' not found; probably the calculation "
"did not run. Check your input"
)
return
# Process the additional xyz file
multiwfn_xyz_coords, error_message = self._process_xyz_file(file_path=_opath2)
if error_message is not None:
self._err = (
f"Error processing Multiwfn output file '{_opath2}': {error_message}. The "
"coordinates need to be separately processed successfully for this feature to "
"be calculated"
)
return
os.remove(_opath2)
assert self.coordinates is not None # for type checker
assert multiwfn_xyz_coords is not None # for type checker
rotation_matrix, translation_vector, error_message = align_coordinates(
reference_coords=self.coordinates, to_be_aligned_coords=multiwfn_xyz_coords
)
if error_message is not None:
self._err = (
f"Error aligning coordinates from Multiwfn output to the coordinates of the mol "
f"object: {error_message}. The coordinates need to be successfully aligned for "
"this feature to be calculated"
)
return
# Rename CPprop.txt file
os.rename(src="CPprop.txt", dst=f"CPprop__{self.conformer_name}.txt")
# Open output file and read data
with open(_opath, "r") as f:
cp_prop_output = f.readlines()
all_data = read_prop_file(
file_content=cp_prop_output,
prefix="bcp_",
rotation_matrix=rotation_matrix,
translation_vector=translation_vector,
)
for data in all_data:
assert isinstance(data["_atoms"], tuple) # for type checker
atom_idx_1 = data["_atoms"][0] - 1
atom_idx_2 = data["_atoms"][1] - 1
# Find the bond in the mol object and save the data
for bond in self.mol.GetBonds():
bond_idx = bond.GetIdx()
if any(
[
(
bond.GetBeginAtomIdx() == atom_idx_1
and bond.GetEndAtomIdx() == atom_idx_2
),
(
bond.GetBeginAtomIdx() == atom_idx_2
and bond.GetEndAtomIdx() == atom_idx_1
),
]
):
# Add data that needs to be separately calculated
distance_begin_atom, begin_atom_coordinates = self._get_distance(
data=data, to="start", bond_obj=bond, mol_obj=self.mol
)
distance_end_atom, end_atom_coordinates = self._get_distance(
data=data, to="end", bond_obj=bond, mol_obj=self.mol
)
relative_bond_position = self._get_relative_bond_position(
begin_atom_coordinates=begin_atom_coordinates,
end_atom_coordinates=end_atom_coordinates,
distance_begin_atom=distance_begin_atom,
)
data["bcp_distance_begin_atom"] = distance_begin_atom
data["bcp_distance_end_atom"] = distance_end_atom
data["bcp_relative_bond_position"] = relative_bond_position
# Correct data types are ensured -> typing is ignored here
self.results[bond_idx] = {
f"multiwfn3D-bond-topology_{feature_id_str}": value # type: ignore[misc]
for feature_id_str, value in data.items()
if feature_id_str not in ["_atoms", "_not_found"]
}
break
[docs]
def _process_xyz_file(
self, file_path: str
) -> Tuple[Optional[NDArray[np.float64]], Optional[str]]:
"""Read and process a xyz file and ensure that it is compatible with the data of the initial
mol object.
Parameters
----------
file_path : str
The path to the xyz file to be processed.
Returns
-------
Tuple[Optional[NDArray[np.float64]], Optional[str]]
A tuple containing:
* The coordinates read from the xyz file as a numpy array, or ``None`` if an error
occurred.
* An error message if an error occurred, or ``None`` if the processing was successful.
"""
# Read file
coords_list, error_message = read_xyz_file(file_path=file_path)
if error_message is not None:
return None, error_message
assert coords_list is not None # for type checker
coords = [c for c in coords_list[0].split("\n")[2:] if c.strip() != ""]
_splitted = np.array([line.split() for line in coords])
atom_symbols = _splitted[:, 0]
# Check if the data was read correctly and matches the mol object
if len(atom_symbols) != len(self.elements):
_errmsg = (
"number of atoms in Multiwfn output does not match number of atoms in mol vault"
)
return None, _errmsg
if not all(atom_symbols == self.elements):
_errmsg = "atom symbols in Multiwfn output do not match atom symbols in mol vault"
return None, _errmsg
# Get the coordinates
atom_coordinates = _splitted[:, 1:4].astype(float)
return atom_coordinates, None
[docs]
@staticmethod
def _get_distance(
data: Dict[str, Any], to: str, bond_obj: Chem.rdchem.Bond, mol_obj: Chem.rdchem.Mol
) -> Tuple[float, NDArray[np.float64]]:
"""Calculate the distance of a bond critical point to one of the atoms forming the bond.
Parameters
----------
data : Dict[str, Any]
The dictionary containing the coordinates of the bond critical point.
to : str
Whether to calculate the distance to the "start" or "end" atom of the bond.
bond_obj : Chem.rdchem.Bond
The RDKit bond object of the bond under consideration.
mol_obj : Chem.rdchem.Mol
The RDKit molecule object of the entire molecule.
Returns
-------
Tuple[float, NDArray[np.float64]]
The distance of the bond critical point to the respective atom and that atom's
cartesian coordinates.
"""
bcp_coordinates = np.array([float(x) for x in data["bcp_coordinates"].split(",")])
if to == "start":
atom_idx = bond_obj.GetBeginAtomIdx()
if to == "end":
atom_idx = bond_obj.GetEndAtomIdx()
pos = mol_obj.GetConformer().GetAtomPosition(atom_idx)
atom_coordinates = np.array([pos.x, pos.y, pos.z])
distance = np.linalg.norm(bcp_coordinates - atom_coordinates)
return float(distance), atom_coordinates
[docs]
@staticmethod
def _get_relative_bond_position(
begin_atom_coordinates: NDArray[np.float64],
end_atom_coordinates: NDArray[np.float64],
distance_begin_atom: float,
) -> float:
"""Calculate the relative position of a bond critical point along its respective bond.
It is always ensured that the ratio is between 0 and 0.5, meaning that it is always
calculated with respect to the closest atom of the bond.
Parameters
----------
begin_atom_coordinates : NDArray[np.float64]
The cartesian coordinates of the start atom of the bond.
end_atom_coordinates : NDArray[np.float64]
The cartesian coordinates of the end atom of the bond.
distance_begin_atom : float
The distance of the bond critical point to the start atom of the bond.
Returns
-------
float
The ratio between the bond length and the distance of the bond critical point to
the bond atom that is closest to it.
"""
bond_length = np.linalg.norm(end_atom_coordinates - begin_atom_coordinates)
ratio = float(distance_begin_atom / bond_length)
if ratio > 0.5:
ratio = 1 - ratio
return ratio
[docs]
class Multiwfn3DBondTopologyBcpAverageLocalIonizationEnergy(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_average_local_ionization_energy",
calculated with multiwfn.
The index of this feature is 446 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpCoordinates(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_coordinates", calculated with
multiwfn.
The index of this feature is 447 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpDeltaGHirshfeld(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_delta_g_hirshfeld", calculated with
multiwfn.
The index of this feature is 448 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpDistanceEndAtom(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_distance_end_atom", calculated with
multiwfn.
The index of this feature is 450 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpDistanceBeginAtom(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_distance_begin_atom", calculated
with multiwfn.
The index of this feature is 451 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpElectronDensity(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_electron_density", calculated with
multiwfn.
The index of this feature is 452 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpElectronDensityAlpha(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_electron_density_alpha", calculated
with multiwfn.
The index of this feature is 453 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpElectronDensityBeta(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_electron_density_beta", calculated
with multiwfn.
The index of this feature is 454 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpElectronDensityEllipticity(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_electron_density_ellipticity",
calculated with multiwfn.
The index of this feature is 455 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpElectronLocalizationFunction(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_electron_localization_function",
calculated with multiwfn.
The index of this feature is 456 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpElectrostaticPotential(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_electrostatic_potential",
calculated with multiwfn.
The index of this feature is 457 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpElectrostaticPotentialElectrons(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_electrostatic_potential_electrons",
calculated with multiwfn.
The index of this feature is 458 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpElectrostaticPotentialNuclearCharges(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature
"topology_bcp_electrostatic_potential_nuclear_charges", calculated with multiwfn.
The index of this feature is 459 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpEnergyDensity(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_energy_density", calculated with
multiwfn.
The index of this feature is 460 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpEtaIndex(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_eta_index", calculated with
multiwfn.
The index of this feature is 461 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpGradientComponentsXYZ(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_gradient_components_x_y_z",
calculated with multiwfn.
The index of this feature is 462 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpGradientNorm(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_gradient_norm", calculated with
multiwfn.
The index of this feature is 463 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpHamiltonianKineticEnergy(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_hamiltonian_kinetic_energy",
calculated with multiwfn.
The index of this feature is 464 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpHessianDeterminant(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_hessian_determinant", calculated
with multiwfn.
The index of this feature is 465 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
[docs]
class Multiwfn3DBondTopologyBcpHessianEigenvalues(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_hessian_eigenvalues", calculated
with multiwfn.
The index of this feature is 466 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpInteractionRegionIndicator(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_interaction_region_indicator",
calculated with multiwfn.
The index of this feature is 467 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpLagrangianKineticEnergy(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_lagrangian_kinetic_energy",
calculated with multiwfn.
The index of this feature is 468 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpLaplacianComponentsXYZ(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_laplacian_components_x_y_z",
calculated with multiwfn.
The index of this feature is 469 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpLaplacianOfElectronDensity(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_laplacian_of_electron_density",
calculated with multiwfn.
The index of this feature is 470 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpLaplacianTotal(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_laplacian_total", calculated with
multiwfn.
The index of this feature is 471 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpLocalizedOrbitalLocator(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_localized_orbital_locator",
calculated with multiwfn.
The index of this feature is 473 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpPotentialEnergyDensity(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_potential_energy_density",
calculated with multiwfn.
The index of this feature is 474 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpReducedDensityGradient(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_reduced_density_gradient",
calculated with multiwfn.
The index of this feature is 475 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpRelativeBondPosition(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_relative_bond_position", calculated
with multiwfn.
The index of this feature is 477 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpSignSecondLargestEigenvalueElectronDensityHessian(
_Multiwfn3DBondTopology
):
"""Feature factory for the 3D bond feature
"topology_bcp_sign_second_largest_eigenvalue_electron_density_hessian", calculated with
multiwfn.
The index of this feature is 478 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpSpinDensity(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_spin_density", calculated with
multiwfn.
The index of this feature is 480 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyBcpVdwPotential(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_bcp_vdw_potential", calculated with
multiwfn.
The index of this feature is 481 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyCpIndex(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_cp_index", calculated with multiwfn.
The index of this feature is 482 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology
[docs]
class Multiwfn3DBondTopologyCpType(_Multiwfn3DBondTopology):
"""Feature factory for the 3D bond feature "topology_cp_type", calculated with multiwfn.
The index of this feature is 483 (see the ``list_atom_features()`` and
``list_bond_features()`` method). The corresponding configuration settings can be found
under "multiwfn.topology" in the _feature_config.toml file.
"""
def __init__(self) -> None:
super().__init__()
# This feature is automatically calculated in _Multiwfn3DBondTopology