# Open SourceMolecular Dynamics

### Sidebar

#### For Developers

exercises:2016_uzh_cmest:faq

## Which basis sets and/or pseudopotentials should I use?

When in doubt, always use the DZVP-MOLOPT-GTH basis set and the GTH-PBE pseudopotential.

In some cases (for example for tin) there is only a short-range basis set available, called DZVP-MOLOPT-SR-GTH.

## How do I obtain a charge density difference?

In Exercise 2 a series of *.cube files was already produced for MOs. To obtain the electronic charge density in a *.cube, file the following snippet of input file can be used (see also the reference manual):

&FORCE_EVAL
&DFT
&PRINT
&E_DENSITY_CUBE
&END E_DENSITY_CUBE
...
&END PRINT
... ...
&END DFT
&END FORCE_EVAL

The cubecruncher tool can be used to obtain charge density differences. A precompiled executable is available on the server tcopt3, its path being /users/scaravat/bin/cubecruncher.x.

The basic usage to obtain a charge density difference is:

## How do I obtain the stress tensor?

You have to add the following snippet to the input file:

&FORCE_EVAL
...
STRESS_TENSOR ANALYTICAL
&PRINT
&STRESS_TENSOR
&END STRESS_TENSOR
&END PRINT
...
&END FORCE_EVAL

Your calculation should be set up in such a way that forces are being calculated: e.g. RUN_TYPE set to GEO_OPT or ENERGY_FORCE; instead RUN_TYPE ENERGY won't work!

## How do I restart a calculation?

There are two restart files usually produced by a run:

• a binary file, typically <PROJECT>-RESTART.wfn, which contains the WFs coefficients.
• a ASCII file, typically <PROJECT>-1.restart, which has the same structure as the input file with all the keyword specified (in case with the default values).

To employ the WFs restart file it's enough to change in the input file SCF_GUESS ATOMIC or SCF_GUESS RANDOM into SCF_GUESS RESTART (see the manual).

Instead, to restart e.g. the positions of the atoms along a GEO_OPT run, you should add a new section to the input file:

&EXT_RESTART
RESTART_FILE_NAME <PROJECT>-1.restart
&END EXT_RESTART

where you should replace <PROJECT> with the content of your PROJECT_NAME (eventually aliased to PROJECT) keyword in section &GLOBAL.