Commit Graph

338 Commits

Author SHA1 Message Date
giannozz b95d77d6c7 First attempt to include new functionals contributed by Ikutaro Hamada.
If I understand correctly: vdw-DF4 and optB86b-vdW are the same; there 
are two implementations of Becke 86, one in b86b and one in becke86b;
two implementation of optimized Becke 86, one in pbex and one in b86b.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@10838 c92efa57-630b-4861-b058-cf58834340f0
2014-04-02 16:21:34 +00:00
giannozz 13a55144da Added Engel-Vosko XC, courtesy of Guido Roma
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@10800 c92efa57-630b-4861-b058-cf58834340f0
2014-03-13 17:37:20 +00:00
giannozz 1d385beca9 SOGGA wasn't really working
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@10656 c92efa57-630b-4861-b058-cf58834340f0
2014-01-02 18:54:17 +00:00
giannozz 3b294825c2 Added new vdw-DF functionals, labelled vdw-DF3 and vdW-DF4, by Klimes et al.
Contributed by Guido Roma. Still untested.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@10583 c92efa57-630b-4861-b058-cf58834340f0
2013-11-06 17:46:39 +00:00
ccavazzoni eabe5f67fd - commit of the latest Axel patch
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@10533 c92efa57-630b-4861-b058-cf58834340f0
2013-10-12 22:28:12 +00:00
dalcorso 49690bc820 Some change to the logic of generate_k_along_lines.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@10464 c92efa57-630b-4861-b058-cf58834340f0
2013-09-21 09:24:27 +00:00
aoterodelaroza 2a6b38c368 Added PW86 (unrevised, PRB 33, 8800 (1986)) and B86B (JCP 85 (1986)
7184) exchange functionals. 'pw86pbe' and 'b86bpbe' are shorthand for
pw86 exchange + pbe correlation and b86b exchange + pbe
correlation. 



git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@10454 c92efa57-630b-4861-b058-cf58834340f0
2013-09-17 01:34:16 +00:00
sjoshua 9d7e462e0e new hybrid functional gau-pbe is added
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@10228 c92efa57-630b-4861-b058-cf58834340f0
2013-05-07 11:55:59 +00:00
giannozz 1e87110327 Poor-man parallelization on atoms of some terms in Hubbard contribution to
force - surprising useful when there are hundreds of Hubbard atoms


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@10200 c92efa57-630b-4861-b058-cf58834340f0
2013-04-29 11:33:04 +00:00
giannozz 2f1d6fe097 Numeric constant given a name (ry_to_ghz) and moved to constants.f90;
make.depend updated; misc. minor things


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@10015 c92efa57-630b-4861-b058-cf58834340f0
2013-02-26 10:03:47 +00:00
dalcorso 1c5ee82018 Added the possibility to give a list of q points in phonon input.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9990 c92efa57-630b-4861-b058-cf58834340f0
2013-02-19 10:35:12 +00:00
giannozz 8d6bdd1d33 Remove blanks at the beginning of directory names. The presence of a white
space at the beginning of "pseudo_dir" was causing a funny MD5 hash,
because C function fopen didin't like it and couldn't open the file 


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9979 c92efa57-630b-4861-b058-cf58834340f0
2013-02-16 18:02:54 +00:00
giannozz 9be6f7c79c Added note about case ibrav=-5 as suggested by Giovanni Pizzi and Andrea
Cepellotti


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9952 c92efa57-630b-4861-b058-cf58834340f0
2013-02-07 16:40:14 +00:00
giannozz 2d58494592 1) Routine "read_input_file" modified so that it accepts the input file name
as an argument instead of reading it from command line. Input file name is
   read instead from command line at startup by mp_startup. 
2) Removed #ifdef __MPI before mp_startup (command line must be read in serial 
   execution as well)
3) Replaced calls to lower-level routines read_namelists, read_cards etc by a
   call to read_input_file in NEB
4) Unused routines removed
Note that the above changes affect only pw.x, cp.x, neb.x; in all other cases,
mp_startup is not executed in serial case (will be fixed sooner or later) and
the command line is still read (by a different routine) at a later stage.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9919 c92efa57-630b-4861-b058-cf58834340f0
2013-02-02 21:12:16 +00:00
giannozz b6fa893c20 Unused routines removed
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9881 c92efa57-630b-4861-b058-cf58834340f0
2013-01-28 17:45:47 +00:00
giannozz fac96b9402 Added As
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9880 c92efa57-630b-4861-b058-cf58834340f0
2013-01-28 17:44:26 +00:00
giannozz 2e7ac5706b Misc I/O problems: some versions of gfortran still do not like iargc
declared as external; a file is accepted as xml if it ends with .xml
or .XML. What is the usage of "attr" by the way?


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9865 c92efa57-630b-4861-b058-cf58834340f0
2013-01-25 13:25:30 +00:00
giannozz c6c82dc0f5 No need any longer for g95 special case
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9833 c92efa57-630b-4861-b058-cf58834340f0
2013-01-21 21:45:52 +00:00
spigafi 2239d9e10b PGI needs iargc properly decalred.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9824 c92efa57-630b-4861-b058-cf58834340f0
2013-01-21 08:41:12 +00:00
giannozz 632047ae09 Parsing of command-line options in neb partly merged with the rest of the code,
simplified


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9823 c92efa57-630b-4861-b058-cf58834340f0
2013-01-20 21:02:07 +00:00
giannozz 652f460dfa There is no reason to have command-line options so looooong. I have introduced
a few shorthands (-i, -ni, -nk, -nb, -nt, -nd) with obvious meaning and updated
the documentation. The old syntax will still work. Cleanup of routines catching
command-line options. More serious cleanup of the entire stuff is anyway needed


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9819 c92efa57-630b-4861-b058-cf58834340f0
2013-01-19 21:52:32 +00:00
giannozz b5b42e892c Safer test: accepts anything starting with <?xml ... > or <xml> in
lowercase or capital letters as valid xml file


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9774 c92efa57-630b-4861-b058-cf58834340f0
2013-01-15 20:07:39 +00:00
giannozz 1dfc944014 "smarter" xml check introduced yesterday wasn't working
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9766 c92efa57-630b-4861-b058-cf58834340f0
2013-01-09 21:35:35 +00:00
giannozz 6dc558d2d6 Subroutine test_input_file converted to function, makes a more careful check
whether a file is xml or not; cleanup pf neb.f90


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9761 c92efa57-630b-4861-b058-cf58834340f0
2013-01-08 21:14:37 +00:00
giannozz fef27547bb Mistake in the new functional pbeq2d
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9424 c92efa57-630b-4861-b058-cf58834340f0
2012-09-12 17:01:06 +00:00
giannozz 30cd3a76b0 Minor cleanup; integration routine prepared for dealing with even
number of grid point (still commented out). I think we should figure 
out which integration routine is the best and stick to it: there are
two simpson-style routines that yield slightly different results


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9416 c92efa57-630b-4861-b058-cf58834340f0
2012-09-10 15:36:39 +00:00
giannozz 866702e9bf Minor changes: source "normalization", documentation
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9411 c92efa57-630b-4861-b058-cf58834340f0
2012-09-05 20:19:51 +00:00
marsamos 31f09febfd added relax cell with constrained area for 2D systems
contribution from Richards Charles Andrew,
Department of Physics, University of Pretoria
South Africa



git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9410 c92efa57-630b-4861-b058-cf58834340f0
2012-09-05 14:45:57 +00:00
riccardo 11e6ab0b2a New definition of rPW86 functional, where the LDA part has been removed
from the GGA definition. The inclusion caused some integration errors



git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9380 c92efa57-630b-4861-b058-cf58834340f0
2012-08-27 11:14:29 +00:00
dalcorso 33b724aac3 Atomic weights updated to element 109.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9239 c92efa57-630b-4861-b058-cf58834340f0
2012-07-28 13:23:05 +00:00
giannozz 54da68f234 List of elements extended to element 109 (now you can study Meitnerium-Hassium
alloys)


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9169 c92efa57-630b-4861-b058-cf58834340f0
2012-07-10 08:49:53 +00:00
giannozz 359f01fc56 First attempt to write the input data for generation into the
pseudopotential file - sems to work 


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9138 c92efa57-630b-4861-b058-cf58834340f0
2012-06-28 22:10:45 +00:00
giannozz dbe8fdcf3c New functional PBEQ2D (PRL 108, 126402 (2012)) courtesy of Letizia
Chiodo - untested for now


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9108 c92efa57-630b-4861-b058-cf58834340f0
2012-06-19 16:54:48 +00:00
ccavazzoni 3072a444e0 - yet another communicator, between images and pools, to parallelize
over perturbations in linear response codes like recomb
  Contributed by Marco Govoni



git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9035 c92efa57-630b-4861-b058-cf58834340f0
2012-05-25 06:46:17 +00:00
giannozz de7b82c60c makedeps.sh wasn't working as expected, make.depend updated
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@9016 c92efa57-630b-4861-b058-cf58834340f0
2012-05-20 17:34:00 +00:00
giannozz f231e7f158 Added SO-GGA and M06L meta-gga, contributed Yves Ferro. Missing 3/5
factor added to potinit. For some obscure reason this produces a worse
starting point in the current metaGGA tests. M06L is numerically very
unstable: there is still something to be understood about the correct
implementation of metaGGA


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@8993 c92efa57-630b-4861-b058-cf58834340f0
2012-05-14 09:26:40 +00:00
spigafi 3878febb32 OpenMP improvements for stresses and forces calculations.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@8685 c92efa57-630b-4861-b058-cf58834340f0
2012-02-18 14:49:55 +00:00
giannozz af20c5cc65 CP can now read PW data in the DFT+U case without crashing. It's a hack
because the way ns occupancies are written to file is very questionable
and should be reconsidered. Cleanup: set_hubbard_l merged and moved to flib/

There are still incompatibilities between PW and CP data (e.g. wfcs 
in LSDA case cannot be read), and there will ever be until all I/O 
routines are completely merged (and not partially as they are now) 


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@8656 c92efa57-630b-4861-b058-cf58834340f0
2012-02-15 15:32:00 +00:00
paulatto e8df497cbc Safe lower bound for string check, or it may crash in phonon (phq_readin, line 488) if fildyn is too short
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@8595 c92efa57-630b-4861-b058-cf58834340f0
2012-02-09 11:30:40 +00:00
paulatto 40b51af447 For a high enough value of n (~about 120 on my system) this function overflows and gives NaN. This should never happen in principle, but can happen in projwfc or similar if the compile
r misreads the value of ngauss (e.g. if it is input a a float ngauss=0.d0).



git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@8427 c92efa57-630b-4861-b058-cf58834340f0
2012-01-10 17:07:41 +00:00
dalcorso 13a671641f Problem fix: hse + lsda was giving NaN in some cases.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@8264 c92efa57-630b-4861-b058-cf58834340f0
2011-11-21 16:41:56 +00:00
giannozz a0203f0e7c Added yet another ibrav=-5 : trigonal and 3fold axis along <111>
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@8030 c92efa57-630b-4861-b058-cf58834340f0
2011-08-01 07:34:00 +00:00
riccardo 556e9b846f Added support for the C09 GGA Exchange functional, contribution by
Ikutaro Hamada - ikutaro@wpi-aimr.tohoku.ac.jp



git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@8013 c92efa57-630b-4861-b058-cf58834340f0
2011-07-27 12:38:37 +00:00
giannozz 4b4b144691 make.depend updated
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@7930 c92efa57-630b-4861-b058-cf58834340f0
2011-07-04 17:13:58 +00:00
marsamos abef4fab91 moved test_input_file for testing xml from NEB to flib
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@7928 c92efa57-630b-4861-b058-cf58834340f0
2011-07-04 15:35:18 +00:00
marsamos 9566f7a717 trimcheck and find_free_unit now in flib
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@7925 c92efa57-630b-4861-b058-cf58834340f0
2011-07-04 15:24:48 +00:00
giannozz 16c44b5d85 Added trigonal R with threfold axis along 111 (ibrav=-5)
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@7920 c92efa57-630b-4861-b058-cf58834340f0
2011-06-30 15:38:44 +00:00
giannozz 59c8ac8b8a Fixed-volume variable cell optimization did not work with hexagonal lattices.
The "corrected' version is likely sub-optimal, but it shouldn't make such a
big difference in practical calculations. Subtraction of isostatic pressure 
from stress moved out of the stress subroutine. Note that the call to the 
stress routine has changed (it has an argument now; since it is called by
pwscf.f90 only, not a big deal)


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@7888 c92efa57-630b-4861-b058-cf58834340f0
2011-06-21 16:53:36 +00:00
giannozz 0a4d5a9225 ibrav=-12 (as 12 with a and c axis orthogonal) seems to work
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@7882 c92efa57-630b-4861-b058-cf58834340f0
2011-06-16 15:43:28 +00:00
giannozz 01a1070dbf Added ibrav=-12 for monoclinic with b as unique axis (i.e. orthogonal to a).
ibrav=12 produces lattice vectors with b tilted axis and c unique axis.
This is the most common choice in crystallography. Documentation updated.
Untested.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@7876 c92efa57-630b-4861-b058-cf58834340f0
2011-06-16 13:08:14 +00:00