Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aebn_v3_0.tar.gz(1519 Kbytes)
Manuscript Title: Introducing PROFESS 3.0: An advanced program for orbital-free density functional theory molecular dynamics simulations
Authors: Mohan Chen, Junchao Xia, Chen Huang, Johannes M. Dieterich, Linda Hung, Ilgyou Shin, Emily A. Carter
Program title: PROFESS
Catalogue identifier: AEBN_v3_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 190(2015)228
Programming language: Fortran 90.
Computer: UNIX/Linux with ifort or gfortran.
Operating system: Linux.
Has the code been vectorised or parallelized?: Yes
RAM: Problem dependent, but 2 GB is sufficient for up to 10,000 ions using default settings.
Keywords: Orbital-free density functional theory, Kinetic energy density functional, First-principles methods, Electronic structure, Molecular dynamics.
PACS: 71.15.Mb, 71.15.Pd.
Classification: 7.3.

External routines: FFTW 3, Libxc, Lapack

Does the new version supersede the previous version?: Yes

Nature of problem:
Given a set of coordinates describing the initial ion positions under periodic boundary conditions, the problem is to determine the ground state energy, electron density, ion positions, and cell lattice vectors predicted by orbital-free density functional theory. The computation of all terms is effectively linear scaling. Parallelization is implemented through domain decomposition, and up to around 10,000 ions may be included in the calculation on just a single processor, limited by RAM.

Solution method:
Two solutions based on computing energies as described in the text. The first method is to minimize the energy with respect to the electron density, ion positions, and cell lattice vectors. The second method is to perform molecular dynamics with different statistical mechanical ensembles.

Reasons for new version:
Conventional OFDFT performs excellently for light metals, such as aluminum, lithium, and magnesium. However, for semiconductors or transition metals, the required approximate non-interacting KEDFs and local pseudopotentials (LPSs) limit the accuracy of OFDFT. Here we present an updated version of PROFESS with a set of new KEDFs designed to model semiconductors or transition metals. Other newly developed functions, such as molecular dynamics and spin-polarized optimizers, are included.

Summary of revisions:

a) Streamlining PROFESS: most module files in PROFESS 2.0 [5] are reconstructed, including splitting large modules into several small ones but with clear functions. These reconstructions improve the readability of PROFESS, reduce the compilation time, and make implementation and contribution of new functions to PROFESS easier for external developers
b) The format of the output file is improved to be more informative and clear.
c) FFTW 2.1.5 is replaced with FFTW3-API interfaces, F95-API in serial mode, allowing for linking against either FFTW3 or a recent MKL and F2003-API in parallel mode.
d) The spin polarized Perdew-Burke-Ernzerhof (PBE) [6] functional is implemented through the Libxc library [7].
e) A subset of the latest Numerical Analysis Library (NMS) with improvements for F95 standard compliance is imported.
f) The HC KEDF is fully supported (with force and stress evaluations) [1].
g) A density decomposition method (DDM) with fixed localized electron density is supported (no force or stress yet) [2].
h) The WGC [8] decomposition (WGCD) KEDF is supported (with force but no stress yet) [3].
i) The EvW-WGC KEDF is supported (with force but no stress yet) [4].
j) A new KEDF based on a point-wise Kohn-Sham [9, 10] kinetic energy density and electron localization function is supported (with force but no stress yet) [11].
k) Spin-polarized calculations are supported [12]. Note that the reference density and therefore the kernels of the non-local KEDFs (except for the HC KEDF where the kernel for each spin channel is computed on-the-fly explicitly) are still simply derived from the total density, not from the electron density of each spin channel.
l) The Limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) method [13, 14] is added for density optimization, which yields better stability and efficiency for spin-polarized and vacuum-containing systems.
m) Molecular dynamics methods with three ensembles are implemented [15] and tested [16]: the microcanonical NVE (constant number of particles N, constant volume V, and constant energy E), the canonical NVT (constant number of particles N, constant volume V, and constant temperature T) with the Nosé-Hoover thermostat [17], and the isothermal-isobaric NPT (constant number of particles N, constant pressure P, and constant temperature T) with the Parinello-Rahman thermostat [18].
n) An option to set the initial density as a superposition of atomic densities is added.
o) A bug in PROFESS 2.0 related to cutoff functions used to treat the vacuum region is fixed.
p) The ion-ion interaction using the Ewald method (without the cubic b-spline method) is now parallelized.

Restrictions:
PROFESS cannot use non-local (such as ultrasoft) pseudopotentials. A variety of LPS files are available at the Carter group website (http://www.princeton.edu/carter/research/local-pseudopotentials).

Running time:
Problem dependent. Timing results for large scale problems are given in the PROFESS paper [19].

References:
[1] C. Huang, E. A. Carter, Phys. Rev. B 81 (2010) 045206.
[2] C. Huang, E. A. Carter, Phys. Rev. B 85 (2012) 045126.
[3] J. Xia, E. A. Carter, Phys. Rev. B 86 (2012) 235109.
[4] I. Shin, E. A. Carter, J. Chem. Phys. 140 (2014) 18A531.
[5] L. Hung, C. Huang, I. Shin, G. S. Ho, V. L. Ligneres, E. A. Carter, Computer Physics Communications 181 (2010) 2208.
[6] J. P. Perdew, K. Burke, M. Ernzehof, Phys. Rev. Lett 77 (1996) 3865.
[7] M. A. Marques, M. J. Oliveira, T. Burnus, Computer Physics Communications 183 (10) (2012) 2272 - 2281.
[8] Y. A. Wang, N. Govind, E. A. Carter, Phys. Rev. B 60 (24) (1999) 16350.
[9] P. Hohenberg, W. Kohn, Phys. Rev. 136 (1964) B864.
[10] W. Kohn, L. J. Sham, Phys. Rev. 140 (1965) A1133.
[11] J. Xia, E. A. Carter, submitted.
[12] J. Xia, E. A. Carter, J. Chem. Phys. 136 (2012) 084102.
[13] J. Nocedal, Math. Comp. 35 (1989) 773.
[14] D. C. Liu, J. Nocedal, Mathematical Programming B 45 (1989) 503.
[15] G. J. Martyna, M. E. Tuckerman, D. J. Tobias, M. L. Klein, Mol. Phys. 87 (5) (1996) 1117.
[16] M. Chen, L. Hung, C. Huang, J. Xia, E. A. Carter, Mol. Phys. 111 (2013) 3448.
[17] S. Nose, J. Chem. Phys. 81 (1984) 511.
[18] W. G. Hoover, Phys. Rev. A 31 (1985) 1695.
[19] G. S. Ho, V. L. Ligneres, E. A. Carter, Comp. Phys. Comm. 179 (2008) 839.