Properties

Functions for calculating/extracting various properties.

Each function can be used to calculate the respective property as is, or to extract it from a passed qmflows.Result instance.

>>> from FOX.properties import get_bulk_modulus
>>> from qmflows.packages import Result
>>> import numpy as np

>>> # Calculate the bulk modulus from a set of arrays
>>> pressure: np.ndarray = ...
>>> volume: np.ndarray = ...
>>> get_bulk_modulus(pressure, volume)  
array([[[ 0.,  1.,  2.],
        [ 3.,  4.,  5.]],

       [[ 6.,  7.,  8.],
        [ 9., 10., 11.]]])

>>> # Calculate the bulk modulus from a qmflows.Result instance
>>> result: Result = ...
>>> get_bulk_modulus.from_result(result)  
array([[[ 0.,  1.,  2.],
        [ 3.,  4.,  5.]],

       [[ 6.,  7.,  8.],
        [ 9., 10., 11.]]])

An example for how get_bulk_modulus() can be used in conjunction with the ARMC yaml input. Note that additional CP2K print keys are required in order for it to export the necessary properties.

job:
    type: FOX.armc.PackageManager
    molecule: mol.xyz

    md:
        template: qmflows.md.specific.cp2k_mm
        settings:
            cell_parameters: [50, 50, 50]
            input:
                motion:
                    print:
                        cell on:
                            filename: ''
                        forces on:
                            filename: ''
                    md:
                        ensemble: NVE
                        thermostat:
                            print:
                                temperature on:
                                    filename: ''

pes:
    bulk_modulus:
        func: FOX.properties.get_bulk_modulus.from_result
        ref: [1.0]
        kwargs:
            reduce: mean

Index

`get_pressure` (forces, coords, volume[…])	Calculate the pressure from the passed forces.
`get_bulk_modulus` (pressure, volume[…])	Calculate the bulk modulus via differentiation of pressure w.r.t. volume.
`get_attr` (obj, name[, default, reduce, axis])	`getattr()` with support for additional keyword argument.
`call_method` (obj, name, *args[, reduce, axis])	Call the name method of obj.
`FromResult` (func[, result_func])	A class for wrapping `FunctionType` objects.

API

FOX.properties.get_pressure(forces, coords, volume, temp=298.15, *, forces_unit='ha/bohr', coords_unit='bohr', volume_unit='bohr', return_unit='ha/bohr^3')[source]

Calculate the pressure from the passed forces.

\[P = \frac{Nk_{B}T}{V} + \frac{1}{6V} \sum_i^N \sum_j^N {\boldsymbol{r}_{ij} \cdot \boldsymbol{f}_{ij}}\]

Parameters

forces (np.ndarray[np.float64], shape \((n_{\text{mol}}, n_{\text{atom}}, 3)\)) – A 3D array containing the forces of all molecules within the trajectory.
coords (np.ndarray[np.float64], shape \((n_{\text{mol}}, n_{\text{atom}}, 3)\)) – A 3D array containing the coordinates of all molecules within the trajectory.
volume (np.ndarray[np.float64], shape \((n_{\text{mol}},)\)) – A 1D array containing the cell volumes across the trajectory.
temp (np.ndarray[np.float64], shape \((n_{\text{mol}},)\)) – A 1D array of the temperatures across the trajectory.
forces_unit (str) – The unit of the forces.
coords_unit (str) – The unit of the coords.
volume_unit (str) – The unit of the volume. The passed unit will automatically cubed, e.g. Angstrom -> Angstrom**3.
return_unit (str) – The unit of the to-be returned pressure.

Returns

A 1D array with all pressures across the trajectory.

Return type

np.ndarray[np.float64], shape \((n_{\text{mol}},)\)

get_pressure.from_result(result, *, reduce=None, axis=None, return_unit='ha/bohr^3', **kwargs)

Call get_pressure() using argument extracted from result.

Parameters

result (qmflows.CP2K_Result) – The Result instance that self should operator on.
reduce (str or Callable[[Any], Any], optional) – A callback for reducing the output of self. Alternativelly, one can provide on of the string aliases from REDUCTION_NAMES.
axis (int or Sequence[int], optional) – The axis along which the reduction should take place. If None, use all axes.
return_unit (str) – The unit of the to-be returned quantity.
**kwargs (Any) – Further keyword arguments for get_pressure().

Returns

The output of get_pressure().

Return type

Any

FOX.properties.get_bulk_modulus(pressure, volume, *, pressure_unit='ha/bohr^3', volume_unit='bohr', return_unit='ha/bohr^3')[source]

Calculate the bulk modulus via differentiation of pressure w.r.t. volume.

\[B = -V * \frac{\delta P}{\delta V}\]

Parameters

pressure (np.ndarray[np.float64]) – A 1D array of pressures used for defining \(\delta P\). Must be of equal length as volume.
volume (np.ndarray[np.float64]) – A 1D array of volumes used for defining \(\delta V\). Must be of equal length as pressure.
pressure_unit (str) – The unit of the pressure.
volume_unit (str) – The unit of the volume. The passed unit will automatically cubed, e.g. Angstrom -> Angstrom**3.
return_unit (str) – The unit of the to-be returned pressure.

Returns

The bulk modulus \(B\). Returend as either a scalar or array, depending on the dimensionality volume_ref.

Return type

np.float64 or np.ndarray[np.float64]

get_bulk_modulus.from_result(result, *, reduce=None, axis=None, return_unit='ha/bohr^3', **kwargs)

Call get_bulk_modulus() using argument extracted from result.

Parameters

result (qmflows.CP2K_Result) – The Result instance that self should operator on.
reduce (str or Callable[[Any], Any], optional) – A callback for reducing the output of self. Alternativelly, one can provide on of the string aliases from REDUCTION_NAMES.
axis (int or Sequence[int], optional) – The axis along which the reduction should take place. If None, use all axes.
return_unit (str) – The unit of the to-be returned quantity.
**kwargs (Any) – Further keyword arguments for get_bulk_modulus().

Returns

The output of get_bulk_modulus().

Return type

Any

FOX.properties.get_attr(obj, name, default=<null>, reduce=None, axis=None)[source]

gettattr() with support for keyword argument.

Parameters

obj (object) – The object in question.
name (str) – The name of the to-be extracted attribute.
default (Any) – An object that is to-be returned if obj does not have the name attribute.
reduce (str or Callable[[Any], Any], optional) – A callback for reducing the extracted attribute. Alternativelly, one can provide on of the string aliases from FromResult.REDUCTION_NAMES.
axis (int or Sequence[int], optional) – The axis along which the reduction should take place. If None, use all axes.

Returns

The extracted attribute.

Return type

Any