User's Guide For Quantum ESPRESSO (v.6.6)
User's Guide For Quantum ESPRESSO (v.6.6)
Contents
1 Introduction 1
1.1 People . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2 Contacts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.3 Guidelines for posting to the mailing list . . . . . . . . . . . . . . . . . . . . . . 4
1.4 Terms of use . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2 Installation 6
2.1 Download . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.2 Prerequisites . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.3 configure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.3.1 Manual configuration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.4 Libraries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.5 Libxc library . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.5.1 Linking in Quantum ESPRESSO . . . . . . . . . . . . . . . . . . . . . 13
2.5.2 Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.5.3 Differences between Libxc and internal functionals . . . . . . . . . . . . . 15
2.5.4 Special cases . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.5.5 XC test . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.6 Compilation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.7 Running tests and examples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.7.1 Test-suite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.7.2 Examples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.8 Installation tricks and problems . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.8.1 All architectures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.8.2 Linux PC . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.8.3 Linux PC clusters with MPI . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.8.4 Microsoft Windows . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.8.5 Mac OS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.9 Cray machines . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.9.1 Obsolescent architectures . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
1
3 Parallelism 24
3.1 Understanding Parallelism . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
3.2 Running on parallel machines . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
3.3 Parallelization levels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
3.3.1 Understanding parallel I/O . . . . . . . . . . . . . . . . . . . . . . . . . 26
3.4 Tricks and problems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
1 Introduction
This guide gives a general overview of the contents and of the installation of Quantum
ESPRESSO (opEn-Source Package for Research in Electronic Structure, Simulation, and Op-
timization), version 6.6.
Important notice: due to the lack of time and of manpower, this manual does not cover
many important aspects and may contain outdated information.
The Quantum ESPRESSO distribution contains the core packages PWscf (Plane-Wave
Self-Consistent Field) and CP (Car-Parrinello) for the calculation of electronic-structure prop-
erties within Density-Functional Theory (DFT), using a Plane-Wave (PW) basis set and pseu-
dopotentials. It also includes other packages for more specialized calculations:
• PWneb: energy barriers and reaction pathways through the Nudged Elastic Band (NEB)
method.
• GWL: electronic excitations within the GW approximation and with the Bethe-Salpeter
Equation
• PWgui: a Graphical User Interface, producing input data files for PWscf and some PostProc
codes.
A copy of required external libraries is also included. Finally, several additional packages that
exploit data produced by Quantum ESPRESSO or patch some Quantum ESPRESSO
routines can be automatically installed using make:
2
• WanT: quantum transport properties with Wannier functions.
• YAMBO: electronic excitations within Many-Body Perturbation Theory, GW and Bethe-
Salpeter equation.
• PLUMED (v.1.3 only): calculation of free-energy surface through metadynamics.
• GIPAW (Gauge-Independent Projector Augmented Waves): NMR chemical shifts and EPR
g-tensor.
For Quantum ESPRESSO with the self-consistent continuum solvation (SCCS) model, aka
“Environ”, see http://www.quantum-environment.org/.
Documentation on single packages can be found in the Doc/ directory of each package.
A detailed description of input data is available for most packages in files INPUT *.txt and
INPUT *.html.
The Quantum ESPRESSO codes work on many different types of Unix machines, includ-
ing parallel machines using both OpenMP and MPI (Message Passing Interface). Quantum
ESPRESSO also runs on Mac OS X and MS-Windows machines (see section 2.2). The GPU-
enabled version is available on https://gitlab.com/QEF/q-e-gpu/-/release.
Further documentation, beyond what is provided in this guide, can be found in:
• the Doc/ and examples/ directories of the Quantum ESPRESSO distribution;
• the web site www.quantum-espresso.org;
• the archives of the mailing list: See section 1.2, “Contacts”, for more info.
People who want to contribute to Quantum ESPRESSO should read the Developer Manual:
Doc/developer man.pdf.
This guide does not explain the basic Unix concepts (shell, execution path, directories etc.)
and utilities needed to run Quantum ESPRESSO; it does not explain either solid state
physics and its computational methods. If you want to learn the latter, you should first read a
good textbook, such as e.g. the book by Richard Martin: Electronic Structure: Basic Theory
and Practical Methods, Cambridge University Press (2004); or: Density functional theory: a
practical introduction, D. S. Sholl, J. A. Steckel (Wiley, 2009); or Electronic Structure Calcula-
tions for Solids and Molecules: Theory and Computational Methods, J. Kohanoff (Cambridge
University Press, 2006). Then you should consult the documentation of the package you want
to use for more specific references.
All trademarks mentioned in this guide belong to their respective owners.
1.1 People
The maintenance and further development of the Quantum ESPRESSO distribution is pro-
moted by the Quantum ESPRESSO Foundation under the coordination of Paolo Giannozzi
(Univ. Udine and IOM-CNR, Italy) and Pietro Delugas (SISSA Trieste) with a strong support
from the MaX - Materials design at the Exascale EU Centre of Excellence.
Contributors to Quantum ESPRESSO, beyond the authors of the papers mentioned in
Sec.1.4, include:
• Ye Luo (Argonne) for improved FFT threading and miscellaneous contributions and op-
timizations;
3
• Fabio Affinito (CINECA) for ELPA support, for contributions to the FFT library, and
for various parallelization improvements;
• Sebastiano Caravati for direct support of GTH pseudopotentials in analytical form, San-
tana Saha and Stefan Goedecker (Basel U.) for improved UPF converter of newer GTH
pseudopotentials;
• Axel Kohlmeyer for libraries and utilities to call Quantum ESPRESSO from external
codes (see the COUPLE sub-directory), made the parallelization more modular and usable
by external codes;
• Daniel Forrer (Padua Univ.) and Michele Pavone (Naples Univ. Federico II) for disper-
sions interaction in the framework of DFT-D;
• Filippo Spiga (University of Cambridge, UK) for mixed MPI-OpenMP parallelization and
for the GPU-enabled version;
• Costas Bekas and Alessandro Curioni (IBM Zurich) for the initial BlueGene porting.
Åke Sandgren, Audrius Alkauskas, Alain Allouche, Francesco Antoniella, Uli As-
chauer, Francesca Baletto, Gerardo Ballabio, Mauro Boero, Scott Brozell, Claudia
Bungaro, Paolo Cazzato, Gabriele Cipriani, Jiayu Dai, Stefano Dal Forno, Cesar Da
Silva, Alberto Debernardi, Gernot Deinzer, Alin Marin Elena, Francesco Filipponi,
Prasenjit Ghosh, Marco Govoni, Thomas Gruber, Martin Hilgeman, Yosuke Kanai,
Konstantin Kudin, Nicolas Lacorne, Hyungjun Lee, Stephane Lefranc, Sergey Lisenkov,
Kurt Maeder, Andrea Marini, Giuseppe Mattioli, Nicolas Mounet, William Parker,
Pasquale Pavone, Mickael Profeta, Chung-Yuan Ren, Kurt Stokbro, David Strubbe,
Sylvie Stucki, Paul Tangney, Pascal Thibaudeau, Davide Tiana, Antonio Tilocca,
Jaro Tobik, Malgorzata Wierzbowska, Vittorio Zecca, Silviu Zilberman, Federico
Zipoli,
4
1.2 Contacts
The web site for Quantum ESPRESSO is http://www.quantum-espresso.org/. Releases
and patches can be downloaded from the links contained in it. The main entry point for
developers is the GitLab web site: https://gitlab.com/QEF/q-e.
The recommended place where to ask questions about installation and usage of Quantum
ESPRESSO, and to report problems, is the mailing list [email protected].
Here you can obtain help from the developers and from knowledgeable users. You have to be
subscribed (see the “Contacts” section of the web site) in order to post to the users’ list. Please
check your spam folder if you do not get a confirmation message when subscribing.
Please read the guidelines for posting, section 1.3! PLEASE NOTE: only messages that
appear to come from the registered user’s e-mail address, in its exact form, will be accepted.
In case of trouble, carefully check that your return e-mail is the correct one (i.e. the one you
used to subscribe). Please also note that since some time the correct functioning of mailing list
is frequently disrupted by policy changes at large commercial providers. There is not much we
can do about that.
If you need to contact the developers for specific questions about coding, proposals, offers
of help, etc., you may either post an “Issue” to GitLab, or send a message to the developers’
mailing list [email protected]. Please do not post general questions
there: they will be ignored.
• Before posting, please: browse or search the archives – links are available in the “Contacts”
section of the web site. Most questions are asked over and over again. Also: make an
attempt to search the available documentation, notably the FAQs and the User Guide(s).
The answer to most questions is already there.
• Reply to both the mailing list and the author or the post, using “Reply to all”.
• Be short: no need to send 128 copies of the same error message just because this is what
came out of your 128-processor run. No need to send the entire compilation log for a
single error appearing at the end.
• Do not post large attachments: post a link to a place where the attachment(s) can be
downloaded from, such as e.g. DropBox, GoogleDocs, or one of the various web temporary
storage spaces.
• Avoid excessive or irrelevant quoting of previous messages. Your message must be imme-
diately visible and easily readable, not hidden into a sea of quoted text.
5
• Remember that even experts cannot guess where a problem lies in the absence of sufficient
information. One piece of information that must always be provided is the version number
of Quantum ESPRESSO.
• Remember that the mailing list is a voluntary endeavor: nobody is entitled to an answer,
even less to an immediate answer.
• Finally, please note that the mailing list is not a replacement for your own work, nor is
it a replacement for your thesis director’s work.
and
Users of the GPU-enabled version should also cite the following paper:
Note the form Quantum ESPRESSO for textual citations of the code. Please also see
package-specific documentation for further recommended citations. Pseudopotentials should
be cited as (for instance)
6
2 Installation
2.1 Download
Quantum ESPRESSO is distributed in source form, but selected binary packages and virtual
machines are also available. Stable and development releases of the Quantum ESPRESSO
source package (current version is 6.6), as well as available binary packages, can be downloaded
from the links listed in the “Download” section of www.quantum-espresso.org.
The Quantum Mobile virtual machine for Windows/Mac/Linux/Solaris provides a com-
plete Ubuntu Linux environment, containing Quantum ESPRESSO and more. Link and
description in https://www.materialscloud.org/work/quantum-mobile.
For source compilation, uncompress and unpack compressed archives in the typical .tar.gz
format using the command:
(a hyphen before ”zxvf” is optional) where X.Y.Z stands for the version number. If your version
of tar doesn’t recognize the ”z” flag:
7
PW/ PWscf package
EPW/ EPW package
NEB/ PWneb package
PP/ PostProc package
PHonon/ PHonon package
PWCOND/ PWcond package
CPV/ CP package
atomic/ atomic package
GUI/ PWGui package
HP/ HP package
Finally, directory COUPLE/ contains code and documentation that is useful to call Quantum
ESPRESSO programs from external codes; directory LR Modules/ contains source files for
modules that are common to all linear-response codes.
2.2 Prerequisites
To install Quantum ESPRESSO from source, you need first of all a minimal Unix environ-
ment, that is: a command shell (e.g., bash, sh) and utilities make, awk, sed. For MS-Windows,
see Sec.2.8.4.
Note that the scripts contained in the distribution assume that the local language is set to
the standard, i.e. ”C”; other settings may break them. Use export LC ALL=C (sh/bash) or
setenv LC ALL C (csh/tcsh) to prevent any problem when running scripts (including installa-
tion scripts).
Second, you need C and Fortran compilers, compliant with C89 and F2003 standards1 . For
parallel execution, you will also need MPI libraries and a parallel (i.e. MPI-aware) compiler.
For massively parallel machines, or for simple multicore parallelization, an OpenMP-aware
compiler and libraries are also required. For GPU’s, you will need the CUDA Fortran version
and to compile with the PGI compiler and run with CUDA run-time environment.
As a rule, Quantum ESPRESSO tries to keep compatibility with older compilers, avoiding
nonstandard extensions and newer features that are not widespread or stabilized. If however
your compiler is very old (say older than ∼ 5 years or so) it is quite likely that something will
not work. The same applies to mathematical and MPI libraries.
Big machines with specialized hardware (e.g. IBM SP, CRAY, etc) typically have a For-
tran compiler with MPI and OpenMP libraries bundled with the software. Workstations or
“commodity” machines, using PC hardware, may or may not have the needed software. If not,
you need either to buy a commercial product (e.g Intel, NAG, PGI) or to use an open-source
compiler like gfortran from the gcc distribution. Some commercial compilers (e.g. PGI) may
be available free of charge under some conditions (e.g. academic or personal usage).
2.3 configure
To install the Quantum ESPRESSO source package, run the configure script. This is ac-
tually a wrapper to the true configure, located in the install/ subdirectory (configure -h
for help). configure will (try to) detect compilers and libraries available on your machine,
and set up things accordingly. Presently it is expected to work on all ”common” computers,
1
since v.6.4 a standard 2008 feature is used: if unallocated pointers are passed as optional arguments, they
are interpreted as not present
8
that is: based on Intel, AMD, ARM CPUs, running Linux, Mac OS X, MS-Windows. Quan-
tum ESPRESSO is known go work on many more kinds of machines but may requires some
tweaking, especially for the hardware of large HPC centers. Detailed but sometimes outdated
installation instructions for specific HPC machines may be found in files install/README.sys,
where sys is the machine name.
Note: support for cMake is being added.
Instructions for the impatient:
cd qe-X.Y.Z/
./configure
make all
This will (try to) produce parallel (MPI) executable if a proper parallel environment is detected,
serial executables otherwise. For OpenMP executables, specify ./configure --enable-openmp.
Symlinks to executable programs will be placed in the bin/ subdirectory. Note that both C
and Fortran compilers must be in your execution path, as specified in the PATH environment
variable.
configure generates the following files:
make.inc compilation rules and flags (used by Makefile)
install/configure.msg a report of the configuration run (not needed for compilation)
install/config.log detailed log of the configuration run (useful for debugging)
include/qe cdefs.h (previously: include/c defs.h) a few definitions used by C files
include/configure.h optionally: info on compilation flags (uncomment #define __HAVE_CO
NOTA BENE: if you interrupt make, it may fail when you start it again later (this will happen
for instance if make is interrupted while unpacking and compiling the FoX library). If so, run
make clean before running make again.
NOTA BENE 1: it is convenient to use ”parallel make” to speed up compilation: make -jN
compiles in parallel on N processors.
NOTA BENE 2: configure no longer updates files make.depend, containing dependencies
upon modules, in the various subdirectories. If you modify the sources, run ./install/makedeps.sh
or type make depend to update files make.depend.
NOTA BENE 3: make.inc used to be called make.sys until v.6.0. The change of name is due
to frequent problems with mailers assuming that whatever ends in .sys is a suspect virus.
You should always be able to compile the Quantum ESPRESSO suite of programs without
having to edit any of the generated files. However you may have to tune configure by specifying
appropriate environment variables and/or command-line options. Usually the tricky part is to
get external libraries recognized and used: see Sec.2.4 for details and hints.
Environment variables may be set in any of these ways:
As a rule: do not define environment variables for configure unless you need it. Always try
configure with no options as a first step. Some environment variables that are relevant to
configure are:
9
ARCH label identifying the machine type (see below)
F90, F77, CC names of Fortran, Fortran-77, and C compilers
MPIF90 name of parallel Fortran 90 compiler (using MPI)
CPP source file preprocessor (defaults to $CC -E)
LD linker (defaults to $MPIF90)
(C,F,F90,CPP,LD)FLAGS compilation/preprocessor/loader flags
LIBDIRS extra directories where to search for libraries
(note that F90 is an “historical” name – we actually use Fortran 2003 – and that it should
be used only together with option --disable-parallel. In fact, the value of F90 must be
consistent with the parallel Fortran compiler which is determined by configure and stored in
the MPIF90 variable).
For example, the following command line:
instructs configure to use mpif90 as Fortran compiler with flags -O2 -assume byterecl, gcc
as C compiler with flags -O3, and to link with flag -static. Note that the value of FFLAGS must
be quoted, because it contains spaces. NOTA BENE: passing the complete path to compilers
(e.g., F90=/path/to/f90xyz) may lead to obscure errors during compilation.
If your machine type is unknown to configure, you may use the ARCH variable to suggest
an architecture among supported ones. Some large parallel machines using a front-end (e.g.
Cray XT) will actually need it, or else configure will correctly recognize the front-end but not
the specialized compilation environment of those machines. In some cases, cross-compilation
requires to specify the target machine with the --host option. This feature has not been
extensively tested, but we had at least one successful report (compilation for NEC SX6 on a
PC). Currently supported architectures are:
x86 64 Intel and AMD 64-bit running Linux
arm ARM machines (with gfortran or armflang)
crayxt4 Cray XT4/XT5/XE machines
mac686 Apple Intel machines running Mac OS X
mingw32 Cross-compilation for MS-Windows, using mingw, 32 bits
mingw64 As above, 64 bits
cygwin MS-Windows PCs with Cygwin
ppc64 Linux PowerPC machines, 64 bits
ppc64-mn as above, with IBM xlf compiler
ppc64-bg* IBM BlueGene
necsx* NEC SX-6 and SX-8 machines
ia32* Intel 32-bit machines (x86) running Linux
ia64* Intel 64-bit (Itanium) running Linux
Note: x86 64 replaces amd64 since v.4.1. Cray Unicos machines, SGI machines with MIPS
architecture, HP-Compaq Alphas are no longer supported since v.4.2; PowerPC Macs are no
longer supported since v.5.0; IBM machines with AIX are no longer supported since v.6.0.
Architectures marked with a * are to be considered obsolescent or obsolete.
Finally, configure recognizes the following command-line options:
10
--enable-parallel compile for parallel (MPI) execution if possible (default: yes)
--enable-openmp compile for OpenMP execution if possible (default: no)
--enable-shared use shared libraries if available (default: yes;
”no” is implemented, untested, in only a few cases)
--enable-debug compile with debug flags (only for some compilers; default: no)
--enable-pedantic compile with pedantic flags (only for gfortran; default: no)
--enable-signals enable signal trapping (default: disabled)
and the following optional packages:
--with-scalapack (yes|no|intel) Use scalapack if available.
Set to intel to use Intel MPI and blacs (default: use OpenMPI)
--with-elpa-include Specify full path of ELPA include and modules headers (default: no)
--with-elpa-lib Specify full path of the ELPA library (default: no)
--with-elpa-version Specify ELPA version, only year (2015 or 2016, default: 2016)
--with-hdf5 (no | | yes | <path>)
Use HDF5, if yes configure assumes
that a valid installation with version ¿= 1.8.16 is
available, and h5cc and h5fc are in the default
executable search path; ¡path¿ must be the root
folder of a standalone hdf5 installation. (default: no).
--with-hdf5-libs= Specify the linker options needed by HDF5 when
configure fails to detect them by itself. As value
to specify here is usually composed by many
substrings it should be enclosed by quotes so to
prevent configure failures. (default: no)
--with-hdf5-include Specify full path the HDF5 include folder containing
module and headers files. Use it if configure fails
to detect the path by itself. (default: no)
--with-libxc Link the libxc library (default:no)
--with-libxc-prefix directory where libxc is installed
--with-libxc-include directory where libxc Fortran headers reside
The following options are available for the CUDA Fortran GPU version:
--with-cuda=value enables compilation of the CUDA Fortran
accelerated subroutines.
value should point the path where the CUDA toolkit
is installed, e.g. /opt/cuda (default: no)
--with-cuda-cc=value sets the compute capabilities for the compilation
of the accelerated subroutines.
value must be consistent with the hardware and the
NVidia driver installed on the workstation or on the
compute nodes of the HPC facility (default: 35)
--with-cuda-runtime=value sets the version of the CUDA toolkit used
for the compilation of the accelerated code.
value must be consistent with the
CUDA Toolkit installed on the workstation
or available on the compute nodes of the HPC facility.
PGI compilers currently accept 7.5, 8.0 or 9.0 (default: 8.0)
11
Please remember that in order to compile the CUDA code for GPU’s you need ... the
CUDA code (it is aligned to but separated from the main distribution)! you also need a recent
version of the PGI Fortran compiler. OpenMP must be enabled, and you may want to use a
CUDA-aware MPI distribution to optimize the data transfer between the processes.
If you want to modify configure (advanced users only!), see the Developer Manual.
2.4 Libraries
Quantum ESPRESSO contains a copy of some needed external libraries:
Optimized vendor-specific libraries often yield huge performance gains with respect to com-
piled libraries and should be used whenever possible. configure always try to locate the best
mathematical libraries.
BLAS and LAPACK Quantum ESPRESSO can use any architecture-optimized BLAS
and LAPACK replacements, like those contained e.g. in the MKL for Intel and AMD CPUs,
or ESSL for IBM Power machines.
If no optimized libraries are available, one may try the ATLAS library: http://math-atlas.sourcef
Note that ATLAS is not a complete replacement for LAPACK: it contains all of the BLAS,
plus the LU code, plus the full storage Cholesky code. Follow the instructions in the ATLAS
distributions to produce a full LAPACK replacement.
Sergei Lisenkov reported success and good performances with optimized BLAS by Kazushige
Goto. The library is now available under an open-source license: see the GotoBLAS2 page at
http://www.tacc.utexas.edu/tacc-software/gotoblas2/.
12
FFT The FFTXlib of Quantum ESPRESSO contains a copy of an old FFTW library. It
also supports the newer FFTW3 library and some vendor-specific FFT libraries. configure
will first search for vendor-specific FFT libraries; if none is found, it will search for an external
FFTW v.3 library; if none is found, it will fall back to the internal copy of FFTW. configure
will add the appropriate preprocessing options:
to DFLAGS in the make.inc file. If you edit make.inc manually, please note that one and only
one among the mentioned preprocessing option must be set.
If you have MKL libraries, you may either use the provided FFTW3 interface (v.10 and
later), or directly link FFTW3 from MKL (v.12 and later) or use DFTI (recommended).
MPI libraries MPI libraries are usually needed for parallel execution, unless you are happy
with OpenMP-only multicore parallelization. In well-configured machines, configure should
find the appropriate parallel compiler for you, and this should find the appropriate libraries.
Since often this doesn’t happen, especially on PC clusters, see Sec.2.8.3.
Note: since v.6.1, MPI libraries implementing v.3 of the standard (notably, non-blocking
broadcast and gather operations) are required.
Libraries for accelerators The accelerated version of the code uses standard CUDA libraries
such as cublas, cufft, cusolver and the eigensolver library explicitly developed for Quan-
tum ESPRESSO by NVidia and distributed at https://github.com/NVIDIA/Eigensolver gpu.
13
The configure script is able to determine whether one is linking to a serial or parallel HDF5
library, and will set the flag -D HDF5 SERIAL in the make.inc file accordingly.
Other libraries Quantum ESPRESSO can use the MASS vector math library from IBM,
if available (only on machines with XLF compiler: likely obsolete).
If optimized libraries are not found The configure script attempts to find optimized
libraries, but may fail if they have been installed in non-standard places. You should exam-
ine the final value of BLAS LIBS, LAPACK LIBS, FFT LIBS, MPI LIBS (if needed), MASS LIBS
(IBM only), either in the output of configure or in the generated make.inc, to check whether
it found all the libraries that you intend to use.
If some library was not found, you can specify a list of directories to search in the envi-
ronment variable LIBDIRS, and rerun configure; directories in the list must be separated by
spaces. For example:
./configure LIBDIRS="/opt/intel/mkl70/lib/32 /usr/lib/math"
If this still fails, you may set some or all of the * LIBS variables manually and retry. For
example:
./configure BLAS_LIBS="-L/usr/lib/math -lf77blas -latlas_sse"
Beware that in this case, configure will blindly accept the specified value, and won’t do any
extra search.
14
• add -D LIBXC to DFLAGS
Then Quantum ESPRESSO can be compiled as usual. If the version of libxc is older than
5.0.0, the last point must be replaced by:
Note for version 5.0.0: the f03 interfaces are no longer available in libxc 5.0.0. They
have been reintroduced in the current develop version. Version 5.0.0 is still usable, but, before
compiling Quantum ESPRESSO, a string replacement is necessary, namely ‘xc f03’ must be
repalced with ‘xc f90’ everywhere in the following files: funct.f90, xc lda lsda drivers.f90,
xc gga drivers.f90, xc mgga drivers.f90, dmxc drivers.f90 and dgcxc drivers.f90 in
Modules folder and xctest qe libxc.f90 in PP/src folder.
2.5.2 Usage
In order to use libxc functionals, you must enforce them from input by including the input dft
string in the system namelist. For example, to use the libxc version of the PBE functionals
(both exchange and correlation):
Capital or lowercase is not relevant, neither the separation between one name and the other.
The complete list of libxc functionals is available at:
https://www.tddft.org/programs/libxc/functionals/
Combinations of Quantum ESPRESSO and libxc functionals are allowed in PW, but some
attention has to be paid to their reciprocal compatibility (see section below). In these cases,
again, the functionals must be enforced from input. For example, the exchange term of PBE
in Quantum ESPRESSO together with the correlation term of PBE in libxc is obtained by
putting:
Note that when using GGA functionals of Quantum ESPRESSO you must always specify
the LDA term too, while it is not the case for the libxc ones.
Outside PW (PHonon, PP, etc.), libxc functionals are still usable, but in these cases it is highly
recommended to avoid combinations of QE and libxc functionals together, since, at present,
they may lead to some issues in the dft detection.
15
2.5.3 Differences between Libxc and internal functionals
There are some differences between Quantum ESPRESSO functionals and libxc ones. In
Quantum ESPRESSO the LDA and the GGA terms are separated and must be specified
independently. In libxc the GGA functionals already include the LDA part (Slater exchange
and Perdew&Wang correlation in most of the cases with the exception, for example, of Lee
Yang Parr functionals).
The libxc metaGGA functionals may or may not need the LDA and GGA terms, depending
on the cases, therefore a good check of the chosen functionals is recommended before doing
expensive runs.
Some functionals in libxc incorporate the exchange part and the correlation one into one
term only (e.g. the ones that include the ‘ xc’ label in their name). In these cases the whole
functional is formally treated as ‘exchange only’ by Quantum ESPRESSO. This does not
imply any loss of information.
External parameters. Several functionals in libxc depend on one or more external pa-
rameters. Some of these can be recovered inside Quantum ESPRESSO, some others are
not directly available. In all these cases a direct intervention on the Quantum ESPRESSO
source code is necessary in order to be able to properly use such functionals.
The libxc routine xc f03 func set ext params allows to set the value of the external pa-
rameters, which are contained in a single array. A trace for the usage of this routine is
present, among the comments, in the XC-drivers of Quantum ESPRESSO (in Modules
folder: xc lda lsda drivers.f90, xc gga drivers.f90, etc.).
Without setting the external parameters inside the code, their default value will be assumed.
This could lead to results different from the expectations. To get the default values, the xc f03-
func info get ext params default value routine is available in libxc. For more details we
refer to the official documentation of libxc.
In any case, when external parameters are needed by the chosen functionals, a warning message
will appear in the output of Quantum ESPRESSO.
Hybrid functionals. The hybrid functionals of libxc are allowed, but, as in the previ-
ous case, their correct operation is not guaranteed by default. It is important to ensure that
the exx fraction in Quantum ESPRESSO assumes the desired value when the XC rou-
tines are called. To do this you need to act on the external parameters of the functionals, as
explained previously. Moreover it is important to ensure that the starting value of the vari-
able exx fraction in Quantum ESPRESSO , which is set in the matching libxc routine
in Modules/funct.f90, is the desired one (exx fraction may not always coincide with the
exx coef of libxc).
16
Functionals with partial output. A few libxc functional routines provides the energy and
some others the potential only, but not both of them. At present, Quantum ESPRESSO,
assumes that both the energy and the potential are provided when a XC-functional routine
is called, therefore, in order to use these functionals, some modifications on the code will be
necessary.
Other functionals. Besides exchange ( x), correlation ( c) and exchange plus correlation
( xc), a fourth kind of functionals is available in libxc, the kinetic functionals ( k). At present,
they are not usable in Quantum ESPRESSO.
2.5.5 XC test
A simple testing program, xctest qe libxc.f90, for the XC part of Quantum ESPRESSO
is available in the PP/src folder.
The program tests the output (energies and potentials) of the XC-driver routines by comparing
the results obtained with the Quantum ESPRESSO internal library of XC functionals with
those obtained with the libxc library.
It also tests the XC-derivative part (mainly used in PHonon) by comparing the QE derivative
of the XC potential (for LDA and GGA) with the libxc one.
2.6 Compilation
The compiled codes can run with any input: almost all variables are dinamically allocated at
run time. Only a few variables have fixed dimensions, set in file Modules/parameters.f90:
These values should work for the vast majority of cases. In case you need more atomic types
or more k-points, edit this file and recompile.
At your choice, you may compile the complete Quantum ESPRESSO suite of programs
(with make all), or only some specific programs. All executables are linked in main bin
directory.
make with no arguments yields a list of valid compilation targets:
17
• make pwall all of the above
• make hp HP package.
• make w90 downloads wannier90, unpacks it, copies an appropriate make.inc file, pro-
duces all executables in W90/wannier90.x and in bin/
• make want downloads WanT, unpacks it, runs its configure, produces all executables for
WanT in WANT/bin.
• make yambo downloads yambo, unpacks it, runs its configure, produces all yambo exe-
cutables in YAMBO/bin
• make gipaw downloads GIPAW, unpacks it, runs its configure, produces all GIPAW exe-
cutables in GIPAW/bin and in main bin directory.
• make plumed unpacks PLUMED, patches several routines in PW/, CPV/ and clib/, recom-
piles PWscf and CP with PLUMED support
For the setup of the GUI, refer to the PWgui-X.Y.Z /INSTALL file, where X.Y.Z stands for the
version number of the GUI (should be the same as the general version number). If you are
using sources from the git repository, see the GUI/README file instead.
If make refuses for some reason to download additional packages, manually download them
into subdirectory archive/, not unpacking or uncompressing them, and try make again. Also
see Sec.(2.1).
2.7.1 Test-suite
Automated tests give a ”pass/fail” answer. All tests run quickly (less than a minute at most),
but they are not meant to be realistic, just to test a specific case. Many features are tested but
only for the following codes: pw.x, cp.x, ph.x, epw.x. Instructions for the impatient:
cd test-suite
make run-tests
18
Instructions for all others: go to the test-suite/ directory, read the README file, or at least,
type make. You may need to edit the run-XX.sh shells, defining variables PARA PREFIX and
PARA POSTFIX (see below for their meaning).
2.7.2 Examples
There are many examples and reference data for almost every piece of Quantum ESPRESSO,
but you have to manually inspect the results.
In order to use examples, you should edit file environment variables, setting the following
variables as needed.
BIN DIR: directory where executables reside
PSEUDO DIR: directory where pseudopotential files reside
TMP DIR: directory to be used as temporary storage area
The default values of BIN DIR and PSEUDO DIR should be fine, unless you have installed
things in nonstandard places. TMP DIR must be a directory where you have read and write
access to, with enough available space to host the temporary files produced by the example
runs, and possibly offering high I/O performance (i.e., don’t use an NFS-mounted directory).
NOTA BENE: do not use a directory containing other data: the examples will clean it!
If you have compiled the parallel version of Quantum ESPRESSO (this is the default
if parallel libraries are detected), you will usually have to specify a launcher program (such
as mpirun or mpiexec) and the number of processors: see Sec.3 for details. In order to do
that, edit again the environment variables file and set the PARA PREFIX and PARA POSTFIX
variables as needed. Parallel executables will be run by a command like this:
$PARA_PREFIX pw.x $PARA_POSTFIX -i file.in > file.out
For example, if the command line is like this (as for an IBM SP):
poe pw.x -procs 4 -i file.in > file.out
you should set PARA PREFIX="poe", PARA POSTFIX="-procs 4". Furthermore, if your machine
does not support interactive use, you must run the commands specified above through the batch
queuing system installed on that machine. Ask your system administrator for instructions. For
execution using OpenMP on N threads, use PARA PREFIX="env OMP NUM THREADS=N ... ".
To run an example, go to the corresponding directory (e.g. PW/examples/example01) and
execute:
./run_example
This will create a subdirectory results/, containing the input and output files generated by
the calculation. Some examples take only a few seconds to run, while others may require up to
several minutes.
In each example’s directory, the reference/ subdirectory contains verified output files,
that you can check your results against. They were generated on a Linux PC using the Intel
compiler. On different architectures the precise numbers could be slightly different, in particular
if different FFT dimensions are automatically selected. For this reason, a plain diff of your
results against the reference data doesn’t work, or at least, it requires human inspection of the
results.
The example scripts stop if an error is detected. You should look inside the last written
output file to understand why.
19
2.8 Installation tricks and problems
2.8.1 All architectures
• Working Fortran and C compilers, compliant with F2003 and C89 standards (see Sec.2.2)
respectively, are needed in order to compile Quantum ESPRESSO. Most recent Fortran
compilers will do the job.
C and Fortran compilers must be in your PATH. If configure says that you have no
working compiler, well, you have no working compiler, at least not in your PATH, and
not among those recognized by configure.
• If you get Compiler Internal Error or similar messages: your compiler version is buggy.
Try to lower the optimization level, or to remove optimization just for the routine that
has problems. If it doesn’t work, or if you experience weird problems at run time, try
to install patches for your version of the compiler (most vendors release at least a few
patches for free), or to upgrade to a more recent compiler version.
• If you get error messages at the loading phase that look like file XYZ.o: unknown / not
recognized/ invalid / wrong file type / file format / module version, one of the following
things have happened:
1. you have leftover object files from a compilation with another compiler: run make
clean and recompile.
2. make did not stop at the first compilation error (it may happen in some software
configurations). Remove the file *.o that triggers the error message, recompile, look
for a compilation error.
If many symbols are missing in the loading phase: you did not specify the location of all
needed libraries (LAPACK, BLAS, FFTW, machine-specific optimized libraries), in the
needed order. Note that Quantum ESPRESSO is self-contained (with the exception
of MPI libraries for parallel compilation): if system libraries are missing, the problem is
in your compiler/library combination or in their usage, not in Quantum ESPRESSO.
• If you get Segmentation fault or similar errors in the provided tests and examples: your
compiler, or your mathematical libraries, or MPI libraries, or a combination thereof, is
buggy, or there is some software incompatibility. Although one can never rule out the
presence of subtle bugs in Quantum ESPRESSO that are not revealed during the
testing phase, it is very unlikely that this happens on the provided tests and examples.
• If all test fails, look into the output and error files: there is some dumb reason for failure.
• If most test pass but some fail, again: look into the output and error files. A frequent
source of trouble is complex function zdotc. See the ”Linux PCs with gfortran compilers”
paragraph, or replace zdotc with fortran intrinsic dot product.
2.8.2 Linux PC
Both AMD and Intel CPUs, 32-bit and 64-bit, are supported and work, either in 32-bit emu-
lation and in 64-bit mode. 64-bit executables can address a much larger memory space than
32-bit executable, but there is no gain in speed. Beware: the default integer type for 64-bit
20
machine is typically 32-bit long. You should be able to use 64-bit integers as well, but it is not
guaranteed to work and will not give any advantage anyway.
Currently, configure supports Intel (ifort), NAG (nagfor), PGI (pgf90) and gfortran com-
pilers. Pathscale, Sun Studio, AMD Open64, are no longer supported after v.6.2: g95, since
v.6.1.
Both Intel MKL and AMD acml mathematical libraries are supported, the former much
better than the latter.
It is usually convenient to create semi-statically linked executables (with only libc, libm,
libpthread dynamically linked). If you want to produce a binary that runs on different machines,
compile it on the oldest machine you have (i.e. the one with the oldest version of the operating
system).
Linux PCs with gfortran You need at least gfortran v.4.4 or later to properly compile
Quantum ESPRESSO.
”There is a known incompatibility problem between different calling convention for Fortran
functions that return complex values [...] If your code crashes during a call to zdotc, recompile
Quantum ESPRESSO using the internal BLAS and LAPACK routines (using configure op-
tions --with-internal-blas and --with-internal-lapack) to see if the problem disappears;
or, add the -ff2c flag” (info by Giovanni Pizzi, Jan. 2013).
If you want to use MKL libraries together with gfortran, link -lmkl gf lp64, not -lmkl intel lp64
(and the like for other architectures).
Linux PCs with Intel compiler (ifort) IMPORTANT NOTE: ifort versions earlier than
v.15 miscompile the new XML code in QE v.6.4 and later. Please install this patch:
https://gitlab.com/QEF/q-e/wikis/Support/Patch-for-old-Intel-compilers.
The Intel compiler ifort http://software.intel.com/ produces fast executables, at least
on Intel CPUs, but not all versions work as expected. In case of trouble, update your version
with the most recent patches. Since each major release of ifort differs a lot from the previous
one, compiled objects from different releases may be incompatible and should not be mixed.
The Intel compiler is no longer free for personal usage, but it is still for students and open-
source contributors (https://software.intel.com/en-us/qualify-for-free-software).
If configure doesn’t find the compiler, or if you get Error loading shared libraries at run
time, you may have forgotten to execute the script that sets up the correct PATH and library
path. Unless your system manager has done this for you, you should execute the appropriate
script – located in the directory containing the compiler executable – in your initialization files.
Consult the documentation provided by Intel.
The warning: feupdateenv is not implemented and will always fail, can be safely ignored.
Warnings on “bad preprocessing option” when compiling iotk and complains about “recom-
mended formats” may also be ignored.
Linux PCs with MKL libraries On Intel CPUs it is very convenient to use Intel MKL
libraries (freely available at https://software.intel.com/en-us/performance-libraries).
MKL libraries can be used also with non-Intel compilers. They work also for AMD CPU,
selecting the appropriate machine-optimized libraries: see “Linux PCs with AMD processors”.
configure properly detects only recent (v.12 or later) MKL libraries, as long as the $MKL-
ROOT environment variable is set in the current shell. Normally this environment variable
21
is set by sourcing the Intel MKL or Intel Parallel Studio environment script. By default the
non-threaded version of MKL is linked, unless option configure --with-openmp is specified.
In case of trouble, refer to the following web page to find the correct way to link MKL:
http://software.intel.com/en-us/articles/intel-mkl-link-line-advisor/.
For parallel (MPI) execution on multiprocessor (SMP) machines, set the environment vari-
able OMP NUM THREADS to 1 unless you know what you are doing. See Sec.3 for more info
on this and on the difference between MPI and OpenMP parallelization.
If you get a mysterious ”too many communicators” error and a subsequent crash: there is
a bug in Intel MPI and MKL 2016 update 3. See this thread and the links quoted therein:
http://www.mail-archive.com/[email protected]/msg29684.html.
Linux PCs with AMD processors For AMD CPUs you may find convenient to link AMD
acml libraries (can be freely downloaded from AMD web site). configure should recognize
properly installed acml libraries. UPDATE (February 2020): ACML have been discontinued.
There are new libraries called AMD Optimizing CPU Libraries (AOCL), tuned for the AMD
EPYC processor family. “ Recently I played around with some AMD EPYC cpus and the bad
thing is that I also saw some strange numbers when using libflame/aocl 2.1. (...) Since version
2020 the MKL performs rather well when using AMD cpus, however, if you want to get the
best performance you have to additionally set:
export MKL_DEBUG_CPU_TYPE=5
which gives an additional 10-20% speedup with MKL 2020, while for earlier versions the speedup
is greater than 200%. [...] Another note, there seem to be problems using FFTW interface of
MKL with AMD cpus. To get around this problem, one has to additionally set
export MKL_CBWR=AUTO
• configure tries to locate a parallel compiler in a logical place with a logical name, but
if it has a strange names or it is located in a strange location, you will have to instruct
configure to find it. Note that in many PC clusters (Beowulf), there is no parallel
Fortran compiler in default installations: you have to configure an appropriate script,
such as mpif90.
• configure tries to locate libraries (both mathematical and parallel libraries) in the usual
places with usual names, but if they have strange names or strange locations, you will
have to rename/move them, or to instruct configure to find them. If MPI libraries are
not found, parallel compilation is disabled.
22
• configure tests that the compiler and the libraries are compatible (i.e. the compiler may
link the libraries without conflicts and without missing symbols). If they aren’t and the
compilation fails, configure will revert to serial compilation.
Apart from such problems, Quantum ESPRESSO compiles and works on all non-buggy,
properly configured hardware and software combinations. In some cases you may have to
recompile MPI libraries: not all MPI installations contain support for the Fortran compiler of
your choice (or for any Fortran compiler at all!).
If Quantum ESPRESSO does not work for some reason on a PC cluster, try first if
it works in serial execution. A frequent problem with parallel execution is that Quantum
ESPRESSO does not read from standard input, due to the configuration of MPI libraries: see
Sec.3.4. If you are dissatisfied with the performances in parallel execution, see Sec.3 and in
particular Sec.3.4.
2.8.5 Mac OS
Mac OS-X machines with gfortran or with the Intel compiler ifort and MKL libraries should
work, but ”your mileage may vary”, depending upon the specific software stack you are using.
Parallel compilation with OpenMPI should also work.
Gfortran information and binaries for Mac OS-X here: http://hpc.sourceforge.net/ and
https://wiki.helsinki.fi/display/HUGG/GNU+compiler+install+on+Mac+OS+X.
Mysterious crashes in zdotc are due to a known incompatibility of complex functions with
some optimized BLAS. See the ”Linux PCs with gfortran” paragraph.
23
$ ./configure --enable-openmp --enable-parallel --with-scalapack
$ vim make.inc
• offload mode, executed on main CPU and offloaded onto coprocessor ”automagically”;
”You can take advantage of the offload mode using the libxphi library. This library offloads the
BLAS/MKL functions on the Xeon Phi platform hiding the latency times due to the communi-
cation. You just need to compile this library and then to link it dynamically. The library works
with any version of QE. Libxphi is available from https://github.com/cdahnken/libxphi.
Some documentation is available therein.
Instead, if you want to compile a native version of QE, you just need to add the -mmic flag
and cross compile. If you want to use the symmetric mode, you need to compile twice: with
and without the -mmic flag”. ”[...] everything, i.e. code+libraries, must be cross-compiled
with the -mmic flag. In my opinion, it’s pretty unlikely that native mode can outperform the
execution on the standard Xeon cpu. I strongly suggest to use the Xeon Phi in offload mode,
for now” (info by Fabio Affinito, March 2015).
IBM BlueGene The current configure was working on the machines at CINECA and at
Jülich. For other machines, you may need something like
where the various * LIBS and * DIR ”suggest” where the various libraries are located.
24
3 Parallelism
3.1 Understanding Parallelism
Two different parallelization paradigms are currently implemented in Quantum ESPRESSO:
1. Message-Passing (MPI). A copy of the executable runs on each CPU; each copy lives in a
different world, with its own private set of data, and communicates with other executables
only via calls to MPI libraries. MPI parallelization requires compilation for parallel
execution, linking with MPI libraries, execution using a launcher program (depending
upon the specific machine). The number of CPUs used is specified at run-time either as
an option to the launcher or by the batch queue system.
2. OpenMP. A single executable spawn subprocesses (threads) that perform in parallel spe-
cific tasks. OpenMP can be implemented via compiler directives (explicit OpenMP) or
via multithreading libraries (library OpenMP). Explicit OpenMP require compilation for
OpenMP execution; library OpenMP requires only linking to a multithreading version of
mathematical libraries, e.g.: ESSLSMP, ACML MP, MKL (the latter is natively multi-
threading). The number of threads is specified at run-time in the environment variable
OMP NUM THREADS.
1. a launcher program such as mpirun or mpiexec, with the appropriate options (if any);
Items 1) and 2) are machine- and installation-dependent, and may be different for interactive
and batch execution. Note that large parallel machines are often configured so as to disallow
interactive execution: if in doubt, ask your system administrator. Item 3) also depend on your
specific configuration (shell, execution path, etc). Item 4) is optional but it is very important
for good performances. We refer to the next section for a description of the various possibilities.
25
3.3 Parallelization levels
In Quantum ESPRESSO several MPI parallelization levels are implemented, in which both
calculations and data structures are distributed across processors. Processors are organized in
a hierarchy of groups, which are identified by different MPI communicators level. The groups
hierarchy is as follow:
• images: Processors can then be divided into different ”images”, each corresponding to a
different self-consistent or linear-response calculation, loosely coupled to others.
• pools: each image can be subpartitioned into ”pools”, each taking care of a group of
k-points.
• bands: each pool is subpartitioned into ”band groups”, each taking care of a group of
Kohn-Sham orbitals (also called bands, or wavefunctions). Especially useful for calcula-
tions with hybrid functionals.
• PW: orbitals in the PW basis set, as well as charges and density in either reciprocal or real
space, are distributed across processors. This is usually referred to as ”PW paralleliza-
tion”. All linear-algebra operations on array of PW / real-space grids are automatically
and effectively parallelized. 3D FFT is used to transform electronic wave functions from
reciprocal to real space and vice versa. The 3D FFT is parallelized by distributing planes
of the 3D grid in real space to processors (in reciprocal space, it is columns of G-vectors
that are distributed to processors).
• tasks: In order to allow good parallelization of the 3D FFT when the number of processors
exceeds the number of FFT planes, FFTs on Kohn-Sham states are redistributed to
“task” groups so that each group can process several wavefunctions at the same time.
Alternatively, when this is not possible, a further subdivision of FFT planes is performed.
Note however that not all parallelization levels are implemented in all codes.
About communications Images and pools are loosely coupled: inter-processors communi-
cation between different images and pools is modest. Processors within each pool are instead
tightly coupled and communications are significant. This means that fast communication hard-
ware is needed if your pool extends over more than a few processors on different nodes.
26
Choosing parameters : To control the number of processors in each group, command line
switches: -nimage, -npools, -nband, -ntg, -ndiag or -northo (shorthands, respectively: -ni,
-nk, -nb, -nt, -nd) are used. As an example consider the following command line:
mpirun -np 4096 ./neb.x -ni 8 -nk 2 -nt 4 -nd 144 -i my.input
This executes a NEB calculation on 4096 processors, 8 images (points in the configuration space
in this case) at the same time, each of which is distributed across 512 processors. k-points are
distributed across 2 pools of 256 processors each, 3D FFT is performed using 4 task groups (64
processors each, so the 3D real-space grid is cut into 64 slices), and the diagonalization of the
subspace Hamiltonian is distributed to a square grid of 144 processors (12x12).
Default values are: -ni 1 -nk 1 -nt 1 ; nd is set to 1 if ScaLAPACK is not compiled, it
is set to the square integer smaller than or equal to half the number of processors of each pool.
Massively parallel calculations For very large jobs (i.e. O(1000) atoms or more) or for
very long jobs, to be run on massively parallel machines (e.g. IBM BlueGene) it is crucial to use
in an effective way all available parallelization levels: on linear algebra (requires compilation
with ELPA and/or ScaLAPACK), on ”task groups” (requires run-time option ”-nt N”), and
mixed MPI-OpenMP (requires OpenMP compilation: configure–enable-openmp). Without a
judicious choice of parameters, large jobs will find a stumbling block in either memory or CPU
requirements. Note that I/O may also become a limiting factor.
• “portable”: as above, but data can be copied to and read from a different machines
(this is not guaranteed with Fortran binary files). Requires compilation with -D__HDF5
preprocessing option and HDF
There is a third format, no longer used for final data but used for scratch and restart files:
• “distributed”: each processor writes its own slice to disk in its internal format to a
different file. The “distributed” format is fast and simple, but the data so produced is
readable only by a job running on the same number of processors, with the same type of
parallelization, as the job who wrote the data, and if all files are on a file system that is
visible to all processors (i.e., you cannot use local scratch directories: there is presently
no way to ensure that the distribution of processes across processors will follow the same
pattern for different jobs).
The directory for data is specified in input variables outdir and prefix (the former can
be specified as well in environment variable ESPRESSO TMPDIR): outdir/prefix.save. A
copy of pseudopotential files is also written there. If some processor cannot access the data
directory, the pseudopotential files are read instead from the pseudopotential directory specified
27
in input data. Unpredictable results may follow if those files are not the same as those in the
data directory!
IMPORTANT: Avoid I/O to network-mounted disks (via NFS) as much as you can! Ideally
the scratch directory outdir should be a modern Parallel File System. If you do not have any,
you can use local scratch disks (i.e. each node is physically connected to a disk and writes to
it) but you may run into trouble anyway if you need to access your files that are scattered in
an unpredictable way across disks residing on different nodes.
You can use input variable disk io to vary the amount of I/O done by pw.x. Since v.5.1,
the dafault value is disk io=’low’, so the code will store wavefunctions into RAM and not
on disk during the calculation. Specify disk io=’medium’ only if you have too many k-points
and you run into trouble with memory; choose disk io=’none’ if you do not need to keep final
data files.
Trouble with input files Some implementations of the MPI library have problems with
input redirection in parallel. This typically shows up under the form of mysterious errors when
reading data. If this happens, use the option -i (or -in, -inp, -input), followed by the input
file name. Example:
Of course the input file must be accessible by the processor that must read it (only one processor
reads the input file and subsequently broadcasts its contents to all other processors).
Apparently the LSF implementation of MPI libraries manages to ignore or to confuse even
the -i/in/inp/input mechanism that is present in all Quantum ESPRESSO codes. In this
case, use the -i option of mpirun.lsf to provide an input file.
Trouble with MKL and MPI parallelization If you notice very bad parallel performances
with MPI and MKL libraries, it is very likely that the OpenMP parallelization performed by
the latter is colliding with MPI. Recent versions of MKL enable autoparallelization by default
on multicore machines. You must set the environment variable OMP NUM THREADS to 1 to
28
disable it. Note that if for some reason the correct setting of variable OMP NUM THREADS
does not propagate to all processors, you may equally run into trouble. Lorenzo Paulatto (Nov.
2008) suggests to use the -x option to mpirun to propagate OMP NUM THREADS to all
processors. Axel Kohlmeyer suggests the following (April 2008): ”(I’ve) found that Intel is now
turning on multithreading without any warning and that is for example why their FFT seems
faster than FFTW. For serial and OpenMP based runs this makes no difference (in fact the
multi-threaded FFT helps), but if you run MPI locally, you actually lose performance. Also
if you use the ’numactl’ tool on linux to bind a job to a specific cpu core, MKL will still try
to use all available cores (and slow down badly). The cleanest way of avoiding this mess is to
either link with
-lmkl intel lp64 -lmkl sequential -lmkl core (on 64-bit: x86 64, ia64)
-lmkl intel -lmkl sequential -lmkl core (on 32-bit, i.e. ia32 )
or edit the libmkl ’platform’.a file. I’m using now a file libmkl10.a with:
It works like a charm”. UPDATE: Since v.4.2, configure links by default MKL without
multithreaded support.
Trouble with compilers and MPI libraries Many users of Quantum ESPRESSO, in
particular those working on PC clusters, have to rely on themselves (or on less-than-adequate
system managers) for the correct configuration of software for parallel execution. Mysteri-
ous and irreproducible crashes in parallel execution are sometimes due to bugs in Quantum
ESPRESSO, but more often than not are a consequence of buggy compilers or of buggy or
miscompiled MPI libraries.
29