BondAngleMap
- class atomman.cluster.BondAngleMap(model: str | IOBase | DataModelDict | None = None, rmin: float | None = None, rmax: float | None = None, rnum: int | None = None, thetamin: float | None = None, thetamax: float | None = None, thetanum: int | None = None, r_ij: _SupportsArray[dtype[Any]] | _NestedSequence[_SupportsArray[dtype[Any]]] | bool | int | float | complex | str | bytes | _NestedSequence[bool | int | float | complex | str | bytes] | None = None, r_ik: _SupportsArray[dtype[Any]] | _NestedSequence[_SupportsArray[dtype[Any]]] | bool | int | float | complex | str | bytes | _NestedSequence[bool | int | float | complex | str | bytes] | None = None, theta: _SupportsArray[dtype[Any]] | _NestedSequence[_SupportsArray[dtype[Any]]] | bool | int | float | complex | str | bytes | _NestedSequence[bool | int | float | complex | str | bytes] | None = None, energy: _SupportsArray[dtype[Any]] | _NestedSequence[_SupportsArray[dtype[Any]]] | bool | int | float | complex | str | bytes | _NestedSequence[bool | int | float | complex | str | bytes] | None = None, symbols: str | list | None = None)
Bases:
objectClass for generating and analyzing energies of three atom clusters that explore a range of interatomic distances and bond angles. Can be used to characterize the 3 atom bond nature of interatomic potentials.
- cumulative_pdf(nbins: int = 301, energymin: float = -15.0, energymax: float = 15.0) Tuple[ndarray, ndarray]
Returns the cumulative probability density function for the energy.
- Parameters:
- Returns:
cum_pdf (numpy.NDArray) – The cumulative probability density function associated with each bin.
centers (numpy.NDArray) – The center values for each bin.
- itercoords() Generator[Tuple[float, float, float, float], None, None]
Iterates through the three-body coordinates, which can be used as inputs for computing energies.
- Yields:
r_ij (float) – The radial distance between atoms i and j.
r_ik (float) – The radial distance between atoms i and k.
r_jk (float) – The radial distance between atoms j and k.
theta (float) – The angle between i-j and i-k in degrees.
- itersystem(symbols: str | list | None = None, copy: bool = False) Generator[System, None, None]
Iterates through the three-body coordinates and returns a System for each. Useful for generating configuration files for simulators. The atom coordinates will be set such that atom 0 is at [0,0,0], atom 1 at [r_ij,0,0] and atom 2 is in the xy plane based on r_ik and theta.
- Parameters:
symbols (str or list, optional) – The element model symbols to assign to the atoms. Can either be one value for all atoms, or three values for each atom individually. If not given here, will use the values set during class initialization.
copy (bool, optional) – If False (default), then the yielded system is the same object with the coordinates shifted. If True, each yielded system is a new object.
- Yields:
atomman.System – The atomic system containing the three-body cluster.
- load_table(filename: str, length_unit: str = 'angstrom', energy_unit: str = 'eV')
Loads a tabulated representation of the coordinates and energy values from a file.
- model(model: str | IOBase | DataModelDict | None = None, length_unit: str = 'angstrom', energy_unit: str = 'eV') DataModelDict | None
Loads or generates a bond angle map data model.
Note: Generating data models is currently limited to regular values, i.e. ones which the coordinates correspond to embedded loops with r_ij iterating in the outside loop, r_ik in the middle loop and theta in the inside loop.
- Parameters:
model (str, file-like object or DataModelDict, optional) – The data model content or file containing the bond angle map data. If given, the content will be read in and set to the current object. If not given, then a data model will be returned for the object.
length_unit (str, optional) – The unit of length to save the rmin and max values in when generating a data model. Default is ‘angstrom’. value
energy_unit (str, optional) – The unit of energy to save the energy values in when generating a data model. Default value is ‘eV’.
- Returns:
The data model containing the bond angle map coordinate information and measured energies. Only returned if model is not given as a parameter.
- Return type:
DataModelDict.DataModelDict
- Raises:
ValueError – If the data is irregular, i.e. coordinates do not conform to embedded loops with r_ij in the outer loop, r_ik in the middle loop and theta in the inside loop.
- pdf(nbins: int = 301, energymin: float = -15.0, energymax: float = 15.0) Tuple[ndarray, ndarray]
Returns the probability density function for the energy
- Parameters:
- Returns:
pdf (numpy.NDArray) – The probability density function associated with each bin.
centers (numpy.NDArray) – The center values for each bin.
- plot_cumulative_pdf(nbins: int = 301, energymin: float = -15.0, energymax: float = 15.0, matplotlib_axes: axes | None = None, **kwargs) figure | None
Generates a plot of the cumulative probability density function of the energy.
- Parameters:
nbins (int, optional) – The number of histogram bins to use. Default value is 301.
energymin (float, optional) – The minimum energy bound to consider. Default value is -15.0.
energymax (float, optional) – The maximum energy bound to consider. Default value is 15.0.
matplotlib_axes (matplotlib.Axes.axes, optional, optional) – An existing plotting axis to add the pdf plot to. If not given, a new figure object will be generated.
**kwargs (any, optional) – Any additional key word arguments will be passed to matplotlib.pyplot.figure for generating a new figure object (if axis is not given).
- Returns:
The generated figure. Not returned if matplotlib_axes is given.
- Return type:
matplotlib.Figure
- plot_pdf(nbins: int = 301, energymin: float = -15.0, energymax: float = 15.0, matplotlib_axes: axes | None = None, **kwargs) figure | None
Generates a plot of the probability density function of the energy.
- Parameters:
nbins (int, optional) – The number of histogram bins to use. Default value is 301.
energymin (float, optional) – The minimum energy bound to consider. Default value is -15.0.
energymax (float, optional) – The maximum energy bound to consider. Default value is 15.0.
matplotlib_axes (matplotlib.Axes.axes, optional, optional) – An existing plotting axis to add the pdf plot to. If not given, a new figure object will be generated.
**kwargs (any, optional) – Any additional key word arguments will be passed to matplotlib.pyplot.figure for generating a new figure object (if axis is not given).
- Returns:
The generated figure. Not returned if matplotlib_axes is given.
- Return type:
matplotlib.Figure
- property rmax: float | None
The maximum value used for the r_ij and r_ik spacings.
- Type:
float or None
- property rmin: float | None
The minimum value used for the r_ij and r_ik spacings.
- Type:
float or None
- property rnum: int | None
The number of values used for the r_ij and r_ik spacings.
- Type:
int or None
- save_table(filename: str, include_header: bool = True)
Saves a tabulated representation of the coordinates and energy values to a file.
- set(rmin: float | None = None, rmax: float | None = None, rnum: int | None = None, thetamin: float | None = None, thetamax: float | None = None, thetanum: int | None = None, r_ij: _SupportsArray[dtype[Any]] | _NestedSequence[_SupportsArray[dtype[Any]]] | bool | int | float | complex | str | bytes | _NestedSequence[bool | int | float | complex | str | bytes] | None = None, r_ik: _SupportsArray[dtype[Any]] | _NestedSequence[_SupportsArray[dtype[Any]]] | bool | int | float | complex | str | bytes | _NestedSequence[bool | int | float | complex | str | bytes] | None = None, theta: _SupportsArray[dtype[Any]] | _NestedSequence[_SupportsArray[dtype[Any]]] | bool | int | float | complex | str | bytes | _NestedSequence[bool | int | float | complex | str | bytes] | None = None, energy: _SupportsArray[dtype[Any]] | _NestedSequence[_SupportsArray[dtype[Any]]] | bool | int | float | complex | str | bytes | _NestedSequence[bool | int | float | complex | str | bytes] | None = None)
Sets the bond angle coordinates and the associated energies, if given.
- Parameters:
energy (array-like, optional) – All measured energies. If r_ij, r_ik and theta are given then all should be the same length. If the range parameters are given then the length should be thetanum*rnum*rnum. If not given, then a new array of nan values will be constructed.
rmin (float, optional) – The minimum value used for the r_ij and r_ik spacings.
rmax (float, optional) – The maximum value used for the r_ij and r_ik spacings.
rnum (float, optional) – The number of values used for the r_ij and r_ik spacings.
thetamin (float, optional) – The minimum value used for the theta angles.
thetamax (float, optional) – The maximum value used for the theta angles.
thetanum (float, optional) – The number of values used for the theta angles.
r_ij (array-like, optional) – All r_ij values used. If given, the lengths of r_ij, r_ik and theta need to be the same.
r_ik (array-like, optional) – All r_ik values used. If given, the lengths of r_ij, r_ik and theta need to be the same.
theta (array-like, optional) – All theta values used. If given, the lengths of r_ij, r_ik and theta need to be the same.
energy – All measured energies. If r_ij, r_ik and theta are given then all should be the same length. If the coordinate range parameters are given, then the energies should be of length rnum*rnum*thetanum and ordered to correspond to three embedded loops with r_ij iterating in the outside loop, r_ik in the middle and theta in the inside. If energy is not given, then all values will initially be set to np.nan.