Effective Screening Medium (ESM) Calculation¶
In this tutorial, we demonstrate how to create a Job in order to extract the potential/charge profiles via the Effective Screening Medium (ESM) approach for simulating surfaces and interfaces, based on Density Functional Theory.
We consider a water (H2O) molecule in the present example, and use Quantum ESPRESSO as our main simulation engine.
Quantum ESPRESSO version considered in this tutorial
The present tutorial is written for Quantum ESPRESSO at versions 5.2.1, 5.4.0, 6.0.0 or 6.3.
Workflow (Quantum ESPRESSO)¶
Expand to view ...
Examples on how ESM is enabled and supported in Quantum ESPRESSO are offered in Ref. 1. Here, we will offer a brief review of the most important input keywords, that are required to be included in Quantum ESPRESSO input scripts in the context of ESM calculations.
Quantum ESPRESSO ESM Input Parameters¶
assume_isolated = 'esm'¶
This parameter is used to perform the calculation assuming the system to be isolated, such as in the cases of a molecule or a cluster, as opposed to regular periodic boundary conditions.
For polarized or charged slab calculation, the 'esm' option embeds the simulation cell within an effective semi- infinite medium in the perpendicular direction (along z). Embedding regions can be vacuum or semi-infinite metal electrodes.If between two electrodes, an optional electric field may be applied via the 'esm_efield' keyword described in what follows.
This requires a simulation cell with the c lattice vector along z, normal to the xy plane, with the slab centered around z=0.
This option determines the boundary conditions used for either side of the slab. The available possibilities are listed in this page.
This keyword determines the position offset of the start of the effective screening region, measured relative to the edge of the simulation cell (of total vertical thickness L_z). The ESM region begins at (assuming the slab to be centered around z=0):
This other option gives the magnitude of the electric field to be applied between semi-infinite ESM electrodes (metals). It is applicable only in the case of the metal-slab-metal (bc2) boundary condition.
lfcpopt option is set to ".TRUE.", it performs a constant bias potential (constant-mu) calculation 2 for a static system with ESM method. This option is subject to the following two conditions:
- calculation must be of type 'relax'.
- Boundary conditions can be of type "bc2" or "bc3" only.
Using the constant-mu method, one can control the Fermi energy, that is the applied bias, during a simulation.
fcp_mu tag in the Quantum ESPRESSO input script sets the target Fermi energy for the simulation, if the aforementioned
lfcpopt input parameter has been set to ".TRUE.".
SCF vs Relax ESM Calculations¶
Two different flavors of ESM workflow calculations are offered on our platform, the first one performing a basic ground state energy self-consistent field (SCF) calculation, whereas the second affording also for the relaxation of the inter-atomic forces, within the structure under consideration, during the course of the ESM computation. The latter option is enabled via the
calculation = 'relax' Quantum ESPRESSO input tag.
Prepare Water Molecule¶
In the present example, we shall opt for the "Vacuum-Slab-Vacuum" (bc1) boundary condition option. The vacuum boundaries should be shifted by half of the lattice c constant, by leaving the "Offset" option of the "Set Boundary Conditions" dialog to its default zero value.
After finishing setting up the boundary conditions for our water molecule structure, the user should Save the changes to the structure into the account-owned materials collection, and then exit Materials Designer.
Import Water Molecule in Job Designer¶
Copy ESM Workflow from Bank¶
Workflows for performing Effective Screening Medium (ESM) computations with Quantum ESPRESSO can readily be imported from the Workflows Bank into the account-owned collection. The user should search for the "ESM" keyword whilst performing a search within the Bank.
Change Important Settings¶
- Type of boundary conditions
- Electric Field
- Target Fermi Energy
In the present example, we shall keep the previously-defined 'bc1' boundary conditions, and leave the remaining three options to their default zero values.
In addition, the user should set the size of the grid of k-points to 1x1x1 in this case, also under "Important Settings", since we are dealing with a water molecule as opposed to a periodic crystal.
Before submitting the job, the user should click on the "Compute" tab of Job Designer and examine the compute parameters included therein. Water is a small structure, so 4 CPUs and a few minutes of calculation runtime should be sufficient.
Examine Final Results¶
Potential Energy Profile¶
When the ESM computation is complete at the end of Job execution, switching to the Results tab of Job Viewer will show the Potential Energy profile of our water-vacuum system, plotted as an energy curve (in eV) as a function of the distance along the vertical perpendicular direction (the "z" coordinate), away from the central water slab. The "local" and "Hartree" contributions to the Potential energy are also given separately.
Charge Density Profile¶
Similarly, the Charge Density profile is also displayed under the Results tab, showing the evolution of the charge density (in electron charge units/Angstrom) as a function of the same vertical "z" coordinate mentioned previously.
We demonstrate the above-mentioned steps involved in the creation and execution of an ESM computation on a water molecule, using the Quantum ESPRESSO simulation engine, in the following animation. Here, we shall make use of the "Relax" variant of the Quantum ESPRESSO ESM workflow.