Calculators¶
For ASE, a calculator is a black box that can take atomic numbers and
atomic positions from an Atoms
object and calculate the
energy and forces and sometimes also stresses.
In order to calculate forces and energies, you need to attach a calculator object to your atoms object:
>>> atoms = read('molecule.xyz')
>>> e = atoms.get_potential_energy()
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/home/jjmo/ase/atoms/ase.py", line 399, in get_potential_energy
raise RuntimeError('Atoms object has no calculator.')
RuntimeError: Atoms object has no calculator.
>>> from ase.calculators.abinit import Abinit
>>> calc = Abinit(...)
>>> atoms.calc = calc
>>> e = atoms.get_potential_energy()
>>> print(e)
-42.0
Here we attached
an instance of the ase.calculators.abinit
class and then
we asked for the energy.
Supported calculators¶
The calculators can be divided in four groups:
Abacus, AMS, Asap, BigDFT, DeePMD-kit, DFTD3, DFTD4, DFTK, FLEUR, GPAW, Hotbit, TBLite, and XTB have their own native or external ASE interfaces.
ABINIT, AMBER, CP2K, CASTEP, deMon2k, DFTB+, ELK, EXCITING, FHI-aims, GAUSSIAN, Gromacs, LAMMPS, MOPAC, NWChem, Octopus, ONETEP, PLUMED, psi4, Q-Chem, Quantum ESPRESSO, SIESTA, TURBOMOLE and VASP, have Python wrappers in the ASE package, but the actual FORTRAN/C/C++ codes are not part of ASE.
Pure python implementations included in the ASE package: EMT, EAM, Lennard-Jones, Morse and HarmonicCalculator.
Calculators that wrap others, included in the ASE package:
ase.calculators.checkpoint.CheckpointCalculator
, thease.calculators.loggingcalc.LoggingCalculator
, thease.calculators.mixing.LinearCombinationCalculator
, thease.calculators.mixing.MixedCalculator
, thease.calculators.mixing.SumCalculator
, thease.calculators.mixing.AverageCalculator
, thease.calculators.socketio.SocketIOCalculator
, the Grimme-D3 potential, and the qmmm calculatorsEIQMMM
, andSimpleQMMM
.
name |
description |
---|---|
DFT supporting both pw and lcao basis |
|
Amsterdam Modeling Suite |
|
Highly efficient EMT code |
|
Wavelet based code for DFT |
|
A deep learning package for many-body potential energy representation |
|
London-dispersion correction |
|
Charge-dependent London-dispersion correction |
|
Plane-wave code for DFT and related models |
|
Full Potential LAPW code |
|
Real-space/plane-wave/LCAO PAW code |
|
DFT based tight binding |
|
Light-weight tight-binding framework |
|
Semiemprical extended tight-binding program package |
|
Plane-wave pseudopotential code |
|
Classical molecular dynamics code |
|
Plane-wave pseudopotential code |
|
DFT and classical potentials |
|
Gaussian based DFT code |
|
DFT based tight binding code |
|
DFT based tight binding |
|
Atomic orbital DFT code |
|
Embedded Atom Method |
|
elk |
Full Potential LAPW code |
Plane-wave pseudopotential code |
|
Full Potential LAPW code |
|
Numeric atomic orbital, full potential code |
|
Gaussian based electronic structure code |
|
Gaussian based electronic structure code |
|
Classical molecular dynamics code |
|
Interatomic potential code |
|
Hessian based harmonic force-field code |
|
Classical MD with standardized models |
|
Classical molecular dynamics code |
|
Combination of multiple calculators |
|
Semiempirical molecular orbital code |
|
Gaussian based electronic structure code |
|
Real-space pseudopotential code |
|
Linear-scaling pseudopotential code |
|
LCAO pseudopotential code |
|
Gaussian based electronic structure code |
|
Enhanced sampling method library |
|
Gaussian based electronic structure code |
|
Gaussian based electronic structure code |
|
LCAO pseudopotential code |
|
Fast atom orbital code |
|
Plane-wave PAW code |
|
Effective Medium Theory calculator |
|
lj |
Lennard-Jones potential |
morse |
Morse potential |
Checkpoint calculator |
|
Socket-based interface to calculators |
|
Logging calculator |
|
DFT-D3 dispersion correction calculator |
|
Explicit Interaction QM/MM |
|
Subtractive (ONIOM style) QM/MM |
Note
A Fortran implemetation of the Grimme-D3 potential, that can be used as an add-on to any ASE calculator, can be found here: https://gitlab.com/ehermes/ased3/tree/master.
The calculators included in ASE are used like this:
>>> from ase.calculators.abc import ABC
>>> calc = ABC(...)
where abc
is the module name and ABC
is the class name.
Calculator keywords¶
Example for a hypothetical ABC calculator:
- ABC(restart=None, ignore_bad_restart_file=False, label=None,
- atoms=None, parameters=None, command='abc > PREFIX.abc',
- xc=None, kpts=[1, 1, 1], smearing=None,
- charge=0.0, nbands=None, **kwargs)
Create ABC calculator
- restart: str
Prefix for restart file. May contain a directory. Default is None: don’t restart.
- ignore_bad_restart_file: bool
Ignore broken or missing restart file. By default, it is an error if the restart file is missing or broken.
- label: str
Name used for all files. May contain a directory.
- atoms: Atoms object
Optional Atoms object to which the calculator will be attached. When restarting, atoms will get its positions and unit-cell updated from file.
- command: str
Command used to start calculation. This will override any value in an
ASE_ABC_COMMAND
environment variable.- parameters: str
Read parameters from file.
- xc: str
XC-functional (
'LDA'
,'PBE'
, …).- kpts:
Brillouin zone sampling:
(1,1,1)
: Gamma-point(n1,n2,n3)
: Monkhorst-Pack grid(n1,n2,n3,'gamma')
: Shifted Monkhorst-Pack grid that includes \(\Gamma\)[(k11,k12,k13),(k21,k22,k23),...]
: Explicit list in units of the reciprocal lattice vectorskpts=3.5
: \(\vec k\)-point density as in 3.5 \(\vec k\)-points per Å\(^{-1}\).
- smearing: tuple
The smearing of occupation numbers. Must be a tuple:
('Fermi-Dirac', width)
('Gaussian', width)
('Methfessel-Paxton', width, n)
, where \(n\) is the order (\(n=0\) is the same as'Gaussian'
)
Lower-case names are also allowed. The
width
parameter is given in eV units.- charge: float
Charge of the system in units of \(|e|\) (
charge=1
means one electron has been removed). Default ischarge=0
.- nbands: int
Number of bands. Each band can be occupied by two electrons.
Not all of the above arguments make sense for all of ASE’s
calculators. As an example, Gromacs will not accept DFT related
keywords such as xc
and smearing
. In addition to the keywords
mentioned above, each calculator may have native keywords that are
specific to only that calculator.
Keyword arguments can also be set or changed at a later stage using
the set()
method:
- ase.calculators.set(key1=value1, key2=value2, ...)¶
- EAM
- Pure Python EMT calculator
- ABINIT
- Amber
- CASTEP
- CP2K
- CRYSTAL14
- Demon
- deMon-Nano
- DFTB+
- DMol3
- Espresso
- exciting
- FHI-aims
- FLEUR
- GAMESS-US
- Gaussian
- Gromacs
- GULP
- Harmonic calculator
- Communication with calculators over sockets
- Jacapo - ASE python interface for Dacapo
- KIM
- LAMMPS Calculators
- LAMMPSrun
- Mopac
- NWChem
- Octopus
- ONETEP
- OpenMX
- ORCA
- PLUMED
- psi4
- Q-Chem
- SIESTA
- Introduction
- Environment variables
- SIESTA Calculator
- Extra FDF parameters
- Example
- Defining Custom Species
- Pseudopotentials
- Restarting from an old Calculation
- Choosing the coordinate format
- Siesta Calculator Class
- Excited states calculations
- Raman Calculations with SIESTA and PyNAO
- Further Examples
- Siesta lrtddft Class
- Siesta RamanCalculatorInterface Calculator Class
- TURBOMOLE
- VASP
- QMMM
- Checkpointing
- Mixing Calculators
- Logging Calculator
- DFT-D3
- Other built-in calculators
- Stuff for testing things
- ACE-Molecule