Single-point calculation#

In this type of project, it is calculated the wave function and charge density, and hence the energy, of the particular arrangement of nuclei defined by the atomistic model.
We refer the user to section Structure Modeling in ASAP for further information on ASAP structure builder/viewer interface.
ASAP sets single-point calculation as the default project type.
Workflow single select
Workflow single select widget

Click on the Parameters icon to open the Single-point Parameter widget,

Workflow single parameters widget

Please refet to the Variables section to know how to edit Atomic positions, Magnetic moments, Atomic masses and Velocities in ASAP.

When Single-point calculation is selected the user can decide to compute additional properties, please refer to the section Properties for a detailed explanation.

Click on the Calculator icon to select the computational engine to be used during the single point calculation.

Workflow single calc select
We refer the user to chapter Calculators for further information on ASAP available calculators.
Click on the Run icon to open Run widget. Then click on the Run button to start the the single-point calculation.
Workflow single run select
Workflow single run widget
We refer the user to chapter Runners for further information on computational resources configuration in ASAP.
After submitting a single-point calculation, the complete calculation output in real-time is shown in the Run widget.
Workflow single run output calc

Analysis#

The information provided by the analysis widget also depends on the properties that were selected on the Parameters widget, see section Properties.

When the calculation is completed, press the Exit and analyse button to open the single-point analysis widget.

Workflow single run exit analyse

The following information extracted from the single-point calculation output is available:

Input#

It shows the complete input file of the single-point calculation. This option is only available for SIESTA calculator.

Workflow single analyse input

Output#

It shows the complete output file of the single-point calculation.

Workflow single analyse output

Energies#

It shows the information on: Total energy, Fermi energy minimum, maximum and frontier orbital energies HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital). As well as energy gap and Eigenvalue energies.

Workflow single analyse energy

We refer the user to chapter Analysis for further information on ASAP analysis widget.

Density of States#

Density Of States describes the number of states per interval of energy, at each energy level, that are available to be occupied.
Check the SCF DOS tick box to visualise the density of states of the system originated from the single-particle (Kohn-Sham) eigenvalues computed during self-consistent field loop.
Workflow single analyse dos
You can adjust the energy range (x-axis) by tuning the minimum and maximum energies (in eV) shown at the top of the figure.
You can also tune the following parameters:
  • Smearing. Gaussian or Lorentzian representations of the Dirac delta function.

  • Full width at half maximum (FWHM) The broadening is defined by the full-width-at-half-maximum in both representations (Gaussian and Lorentzian) of the Dirac delta function.

  • Zero @ Fermi. Sets the origin of the DOS function at Fermi level.

    Workflow single analyse zero fermi

Right-click on the density of states figure to make available the following options:

Workflow single analyse dos options
  • Save figure as…. We support the following formats: pdf, png, jpg, jpeg, ps, eps, and svg.

  • Open in Matplotlib… Open with Matplotlib the python script used to plot the figure.

  • Save a Matplotlib script as… Save the python script used to plot the figure.

  • Export data… Export the data from the density of states figure to a text file format.

Check the Post-SCF-DOS tick box to visualize the density of states computed after a self-consistent field loop with a custom BZ sampling.

Workflow single analyse post scf dos

You can adjust the energy range (x-axis) by tuning the minimum and maximum energies shown at the top of the widget, in units of eV, meV, kcal/mol, kJ/mol, Ha, mHa, Ry and mRy.

Workflow single analyse post scf dos energy

Check the Zero @ Fermi level tick box to set the origin of the post-SCF DOS function at the Fermi energy level.

Workflow single analyse post scf dos fermi

The electronic band structure of a solid describes the range of energies that an electron within the solid may have (energy bands) and the ranges of energies it may not have (band gaps). Check the Band structure tick box to visualise the corresponding graph. Please note that this option is only available for infinite-sized systems.

Workflow analyse band structure

Notice that you can adjust the band structure energy range (y-axis) by tuning the minimum and maximum energies shown at the top of the figure, in units of eV, meV, kcal/mol, kJ/mol, Ha, mHa, Ry and mRy.

Workflow analyse band structure energy

Check the Zero @ Fermi level tick box to center the DOS at the Fermi energy level.

Workflow analyse band structure fermi

Right-click on the band structure figure to make available the following options:

Workflow analyse band structure options
  • Save figure as…. We support the following formats: pdf, png, jpg, jpeg, ps, eps and svg.

  • Open in Matplotlib…. Matplolib library allows the user an interactive visualisation of the figure (scale axis, zoom on speci1c X and Y values…).

  • Save a Matplotlib script as…. Save the Matplotlib script used to plot the figure.

  • Export band structure to JSON…

  • Export band structure to .bands…

  • Export DOS…. Export the data from the band structure figure to a text file format.

Partial DOS#

Check the PDOS option to visualise the partial density of states of the system.

Analysis pdos widget

To visualise the projected DOS on specific atoms, select them by their chemical species or choose a range of atom indices.

Analysis pdos widget atom selection

You can also select the angular momenta of the atomic orbitals that will be included in the PDOS visualisation.

Analysis pdos widget angular selection

You can change the label at your convenience. Press the Add curve button to incorporate each of your selections.

Analysis pdos widget add

You can change units on the top right side of the widget. The options are eV, meV, kcal/mol, kJ/mol, Ry, mRy, Ha and mHa.

Analysis pdos widget updated
You can adjust the energy range (x-axis) by tuning the minimum and maximum energies shown at the top of the figure.
Click on the Zero at Fermi level check box to center DOS at the fermi level.
Analysis pdos widget zero fermi

You can change the colors of the curves by tuning the value on the Color column in the table on the bottom left of the widget. Please note that you must use Hexadecimal code.

Analysis pdos widget colors

Additionally, you can change the style of the lines by tuning the value on the column Style on the bottom left of the widget. Supported values are ‘-’, ‘–’, ‘-.’, ‘:’, ‘None’, ‘ ‘ , , ‘solid’, ‘dashed’, ‘dashdot’, ‘dotted’.

Analysis pdos widget style

Charge#

Option only available if at least one method has been selected at Charge Analysis. See subsection Charge.
Click on the Charge Analysis… button to visualise the atomic charge.
Analysis properties charge

You can visualise the charge of each atom of the project calculated by the selected methods. These charges are given in units of electrons.

Bader analysis#

When the project is successfully completed, click the Exit & analyse button.

Analysis properties bader exit analyse

Then select Bader analysis in the Analysis widget.

Analysis properties bader bader analysis

Click the View in 3D button to view the structure in the advanced 3D viewer. You can see the equipotential surface induced by the charge distribution for a determined isovalue. See section Advanced 3D viewer for further information on the 3D advanced viewer.

Analysis properties bader 3d
Analysis properties bader 3d view

In the previous figure you can see the structure of the project in 3D and the Bader charge isosurface.

Press the Export to cube button to export the volumetric data of the Bader analysis to a cube file.

Analysis properties bader export

Electrostatic potential#

When the calculation is completed, select Exit and analyse… to open the analysis widget.
Check the tick-box Potential to visualise the electrostatic potential as a function of the distance. By default, the widget shows the average electrostatic potential (in eV units) in cross-section (z direction in Å).
The planar averaged electrostatic potential (\(V(z)\)) is calculated across the system as:
\[V(z)=1/S_{xy}\int \int V(x, y, z) dx dy,\]
where \(S_{xy}\) is the cross-section of the system.
Analysis properties electrostatic potential
In this example, we have computed the electrostatic potential for the frozen bulk geometry of freestanding PbZrTiO\(_3\)(001). It presents a tetragonal phase characterised by alternating columnar arrangement of Ti and Zr cations along the polar axis, \(z\).
The planar averaged of the electrostatic potential (red curve in the figure above) exhibits fast oscillations correlated with the atomic arrangement in the material. (Massimiliano Stengel and Nicola A. Spaldin, Phys. Rev. B 75, 20512, 2007). On the other hand, you can observe that the potential is flat at the vacuum region (0-18 Å  and 34 to 50 Å).
The potential drop observed at 0 Å  is due to the fact that “Slab.DipoleCorrection True” was requested to compensate for the electric field across the whole system.
Check the tick-box Macroaverage to visualise the macroscopically averaged potential. It is a convenient approach to evaluate information related to the variation of the electrostatic potential but discarding its fast oscillation related to atomic arrangement. (Massimiliano Stengel and Nicola A. Spaldin, Phys. Rev. B 75, 20512, 2007)
Analysis properties electrostatic potential macroaverage
When Macroaverage is activated, you can modify two related parameters: the Length for averaging (in Units of Å, Bohr or nm) and the damping (dimensionless). The damping factor controls the local oscillations. The selection of a damping factor that is close to the lattice parameter smoothens local oscillations and makes the slope of the electrostatic potential easier to analyse.
You can also tune the following settings in the analysis widget:
  • Axis: The axis (X, Y or Z) in which the electrostatic potential is evaluated.

  • Potential: The units of the electrostatic potential: eV, meV, Rydberg or Hartree.

  • Range: The range used as the x-axis. The units can be Å, Bohr and nm.

In this second example, we have computed the electrostatic potential of the tetracyanoquinodimethane molecule (TCNQ). TCNQ has a planar structure that consists of a central quinoid ring (C6H4) and four cyano groups (CN) attached to the corners.
Analysis properties electrostatic potential tcnq macro
The planar averaged of the electrostatic potential (red curve in the figure above) exhibits the variation of the electrostatic potential along the X axis correlated with the atomic arrangement in the material.
The analysis widget also offers the possibility of visualising the electrostatic potential by selecting View in 3D.
Analysis properties electrostatic potential tcnq 3d appear
The cyano groups draw electron density away from the quinoid ring. Consequently it becomes electron-defficient, leading to a negative electrostatic potential around the carbon atoms of the ring.
On the other hand, nitrogen atoms in the cyano groups are electron-rich since they have they have lone pairs of electrons that are not shared with other atoms in the molecule. As a consecuence, the regions around the nitrogen atoms present positive electrostatic potential, approximately in [0,5.4] Ang and [14.5,20] Ang intervals.
Press the Export to cube button to export the volumetric data of the computed electrostatic potential to a cube file.
Analysis properties electrostatic potential tcnq export