Tutorial: Energy evaluations
Qcore provides efficient energy calculations using a variety of standard electronic structure methods, including DFT, GFN-xTB, and MP2. Moreover, entos provides quantum embedding approaches, including embedded mean-field theory (EMFT) and ONIOM, to target computational resources at the high-accuracy description of important regions. In this section, we will demonstrate single-point energy calculations using DFT, GFN-XTB, and EMFT.
1. Density functional theory
The minimal input for running a DFT calculation looks like this:
This input corresponds to a DFT energy calculation on the built-in molecule water using the PBE exchange-correlation functional and the 6-31G* AO basis set. All other calculation details such as charge, spin multiplicity, SCF convergence assume default values.
Here is another input file with a greater level of detail:
In this example, we highlight the following options:
The molecular structure can be specified in-line using option xyz; alternatively, it can be read from an XYZ file, using the file option within the structure()subcommand.
The total charge of the molecule is specified using the charge option.
The spin multiplicity of the molecule is specified using the multiplicity option.
The maximum number of SCF iterations is specified using the max_iter option.
The print_level option controls the verbosity of output.
Additional options/subcommands for performing DFT calculations can be found in the user manual.
The minimal input for running a GFN-xTB calculation looks like this:
As for DFT, the total charge, spin multiplicity and print level can also be specified within the xtb() command.
In cases where the SCF is difficult to converge, we can choose different values for the Fock matrix damping:
EMFT is a quantum-embedding method in which the atoms are partitioned into two subsystems, with the "active" subsystem typically treated at a higher level of mean-field theory and the "environment" subsystem treated at a lower level of mean-field theory. See this paper for a detailed description of the method.
The minimal input for an EMFT energy calculation looks like this:
The EMFT energy calculation is run through the emft() command, which includes the following options/subcommands:
The active subsystem is specified using the active option, whose values are a list of atom indices; the environment subsystem includes the remaining atoms.
The high-level and low-level of theory are specified by the dft()subcommand, each specifies the exchange-correlation functional and AO basis set for each subsystem.
In the above example, the high-level subsystem corresponds to the –OH group of the methanol molecule (i.e. atoms 1 and 2 in the molecular structure specification); the high-level DFT method is B3LYP with Def2-TZVP basis set, and the low-level DFT method is PBE with Def2-SVP basis set.
Examples on specific topics: