Fukui Functions and the Dual Descriptor¶

Step 1: Setting up the calculation¶

Start ADFinput

Build the carbon monoxide molecule using the search function

Set the XC functional to GGA BP

Open the Fukui tab in Properties → Fukui Functions

Check both Fukui options

Here you can specify which Fukui functions you want to calculate. Both can be calculated at the same time, and is also required for the Dual Descriptor. The default way to calculate the Fukui functions is to add or subtract a whole electron. However, it is possible to use fractional charges by modifying the Charge change parameter.

Open the Bader tab in Properties → Other: Etot, Bader, Charge Transport, ...

Check the Bader option for Critical points, bond paths and atomic properties

The condensed Fukui Functions and local softness are calculated using atomic charges. By default the Hirshfeld, Voronoi and Mulliken charges are used. The Bader charges can be used as well if they are calculated.

Use File → Save As to save the file as ‘CO_Fukui’

Run the calculation

Step 2: The output¶

The condensed Fukui Functions and the local softness can be found in the output file of the main calculation.

Open the output file with ADFoutput

Go to the Condensed Fukui section by Properties → Condensed Fukui

The condensed Fukui functions are given per atom, with the various atomic charge methods. The same applies for the local softness.

Step 3: Visualizing the Fukui functions and Dual Descriptor¶

The Fukui functions can be visualized in ADFview, using the electron density.

Set the grid size to medium Fields → Grid → Medium

Add a double isosurface Add → Isosurface: Double (+/-)

At the bottom of the screen choose the Fukui Plus Function

The instruction on the screen asks for the result file of the Fukui Plus calculation

Choose CO_Fukui.plus.t21

Enable the second checkbox for the field to use the wireframe model

ADFview will create some calculated fields to combine the spin densities for the unrestricted calculations, and calculates the Fukui plus function. The wireframe model allows us to see any hidden isosurface. Another way to do this is by lowering the opacity in the field details.

Disable the field

Create a new double isosurface

At the bottom of the screen choose the Dual Descriptor in the Fukui menu

The fukui plus function is already loaded, the instruction now asks for the missing Fukui Minus calculation

Choose CO_Fukui.minus.t21

Enable the second checkbox for the field to use the wireframe model