##### Copyright (C) 2008 ABINIT group (XG)
This file is distributed under the terms of the GNU General Public License, see ~abinit/COPYING or http://www.gnu.org/copyleft/gpl.txt .
For the initials of contributors, see ~abinit/doc/developers/contributors.txt .

    Many thanks to the contributors to the ABINIT project
    between December 2007 and February 2009. These release notes
    are relative to improvements of ABINITv5.6 with respect to v5.5.

    The list of contributors includes :
    B. Amadon, E. Arras, J.-M. Beuken, F. Bottin, F. Bruneval, D. Caliste, Chunping Hu, M. Cote,
    T. Deutsch, M. Giantomassi, D. Hamann, A. Jacques, J. Janssen,
    F. Jollet, T. Jourdan, S. Leroux, C. Motta, S. Pesant,
    Y. Pouillon, T. Rangel, G.-M. Rignanese, R. Shaltaf, L. Sindic, M. Torrent,
    M. Verstraete, G. Zerah, J. Zwanziger

    It is worth to read carefully
    all the modifications that are mentioned in the present file,
    and examine the links to help files or test cases ...
    This might take some time ...
    Please note the WARNINGS !

    Xavier

    __________________________________________________________________

    Version 5.6, released on September 3, 2008.

    Changes with respect to version 5.5 :

    A. WARNINGS.
    B. Most noticeable achievements (for users)
    C. Changes in the package, for developers
    D. Other changes (or on-going developments, not finalized)

    __________________________________________________________________

    A.  WARNINGS

    A.1 The input variable ratsph is now an array ratsph(ntypat).

    A.2 The input variable exchn2n3 has been replaced by exchn2n3d,
        because the former was terminated with a digit, which is incompatible
        with the indexing of the input variables by the dataset number.

    A.3 The interface to the configure system for the production of the parallel version of abinit (abinip),
        has changed. In brief, if you were using the "manual" generation previously,
        you likely have to switch to the "auto-detection" (enable_mpi=yes).
        Also, with enable_mpi=yes, the user should provide the path to
        the correct MPI directory prefix, e.g.
        with_mpi_prefix="/usr/local/mpich-1.2.7p1-intel-9.1-64".
        The "manual" generation of the MPI support (enable_mpi=manual)
        needs many configure variables to be define explicitely, unlike
        previously.

        More details :
        Interface changes to the configure system :
        configure: reporting user interface changes:
        configure:  * new option --enable-error-handlers is available
        configure:  * new option --enable-gw-dpc is available
        configure:  * --enable-mpi has a new meaning
        configure:    please check that enable-mpi="yes" is OK for you
        configure:  * new option --enable-scalapack is available
        configure:  * --with-fc-name renamed --with-fc-vendor
        configure:  * new option --with-mpi-level is available
        configure:  * new option --with-plugins-prefix is available
        configure:  * new option --with-plugins-srcdir is available
        configure:  * new option --with-scalapack-includes is available
        configure:  * new option --with-scalapack-libs is available
        configure:  * --enable-tricks has a new meaning
        configure:    please check that enable-tricks="yes" is OK for you

    __________________________________________________________________

    B.  Most noticeable achievements

    B.1 The PAW can be used to perform GW calculations.
        Compared to the norm-conserving case, this will allow
        to get closer to all-electron calculations,
        and also will allow to gain CPU (when the planewave energy cut-off to be used is smaller)
        This new capability of ABINIT is exemplified in v5#66 .
        Work done by Matteo Giantomassi.

    B.2 The convergence of the GW calculations with respect to the number of bands
        has been improved dramatically, thanks to the extrapolar method,
        described in Phys.Rev.B78, 085125 (2008).
        This should results in a 3-5 speed up of GW calculations for realistic systems,
        thanks to a large lowering of the number of unoccupied bands needed to obtain the same accuracy.
        See the new input variables gwcomp and gwencomp, and the new
        automatic tests v5#65, 69 and paral#V .
        Work by Fabien Bruneval. Test v5#69 provided by Matteo Giantomassi.

    B.3 Non-collinear magnetism (with spinors) is available for PAW calculations.
        Also, spin-orbit interaction is available for PAW calculations.
        For the spin-orbit calculation, the new input variable pawnsporb has
        been defined.
        Automatic test v5#17
        Work done by Marc Torrent, with participation of T. Jourdan and E. Arras.

    B.4 LDA+U is now available in production (always need PAW).
        See the (already existing) input variables usepawu, upawu, jpawu
        as well as some new ones : dmatudiag, lpawu, dmatpuopt, usedmatpu
        (and the internal variable natpawu).
        Automatic tests v5#08, 11, 19, 53
        Work done by Bernard Amadon.

        LDA+U is available with non-collinear magnetism and with spin-orbit interaction
        Automatic test v5#16 for non-collinear magnetism
        Work done by M. Torrent

        A tutorial is available to explain the use of LDA+U,
        see the menu of the tutorials, work by Simon Pesant, helped
        by M. Coté and XG. In this tutorial, several additional examples
        are given.

        A lm-dependent energy shift can be applied atom by atom,
        see the input variables natvshift, atvshift. This would e.g. allow
        to compute the value of the U according to the method by
        Coccochioni and de Gironcoli.
        Automatic tests v5#14 and 15 are available. Work done by XG.

    B.5 In the framework of band/FFT/k-point parallelisation,
        spin-polarized (nsppol=2 or nspden=2) as well as PAW method calculations
        are now allowed.
        PAW case is shown in the test paral#R .
        Precondition of density and/or correction of forces (see iprcch) and time-reversal
        symmetry (istwfk=2) are now available for beta tests.
        Work done by F. Bottin, with M. Torrent, F. Jollet and G. Zerah.

        Also allowed within band/FFT/k-point parallelization : optics_paw.
        Work by M. Torrent.

        The computation of the static susceptibility is now compatible with
        band/FFT/k-point parallelisation (done in sequential on the nfftdiel grid).
        Work done by G. Zerah.

    B.6 The static susceptibility can now be computed within PAW, opening access
        to the extrapolar algorithm to speed up and/or stabilize the SCF cycles.
        The iprcel input variable can thus be used with the following values : 4x, 6x, 14x .
        Automatic test paral#R.
        Work done by M. Torrent.

    B.7 The Wannier90 library, interfaced with ABINIT, can now
        be used in production,
        work by Tonatiuh Rangel, Bernard Amadon, and Don Hamann, helped by G.-M. Rignanese, as well as
        Y. Pouillon for the build system.
        New test has been provided in the tests/wannier90 directory (B. Amadon).
        Related input variables w90iniprj, w90lplot, w90nplot, w90prtunk.
        However, the testing procedure is still rather incomplete. Only the .out is tested.

    B.8 Anti ferromagnetic symmetries can be used for GW calculations,
        providing significant speed up for antiferromagnetic materials.
        Test v5#68. Work done by Matteo Giantomassi.
        See also section D.15 .

    B.9 A cut-off for the Coulomb interaction, to be used for GW calculations
        of nanosystems (clusters and wires) has been introduced.
        See the new input variable vcutgeo, as well as icutcoul, and the related rcut.
        One test is provided, paral#V , for a finite system.
        More tests should be provided for the case of the cylinder geometry.
        Work done by M. Giantomassi. Test paral#V provided by F. Bruneval.

    B.10 Local exchange (=atom-centered-only part of the exchange)
         has been coded within the PAW.
         Input variables : useexexch, lexexch.
         Automatic test v5#18 is available, testing local part of PBE0 hybrid functional.
         Work done by Francois Jollet.

    B.11 The ability to compute the Fermi contact interaction at each nuclear
        site, that is, the electronic density at each nuclear site, has been coded (PAW).
        This is needed for the computation of isomer shifts to model Mossbauer spectra.
        Input variable prtfc.  Automatic test v5#34.
        Also, the input variable prtefg, coupled with ptcharge,
        allows to compute the Electric Field Gradient using an ionic point charge model .
        Automatic test number v5#35.
        Work done by Joe Zwanziger, with some help from M. Torrent.

    B.12 For very high electronic temperatures, the direct computation
        of the density based on the recursion method (instead of the
        iterative determination of the eigenfunctions) has been coded.
        Input variables recnpath, recnrec, recptrott, rectesteg,
        recrcut, recefermi, rectolden.
        Still to be documented. Automatic test paral T (to be activated !).
        Reference ?
        Work by Stephane Leroux and Gilles Zerah.

    B.13 Make the ABINIT parser independent for third party codes
         through a libabinis.a library.
         Introduction of language bindings for C and Python.
         See src/bindings/README .
         Based on the libabinis.a library, one can parse an ABINIT input
         file and access all scalar and array data of all datasets on this file through
         get routines using the name of the attribute as a keyword.
         A simple pyGtk interface has been created to read ABINIT input
         file to demonstrate the capabilities of Python bindings.
         Further developments will be done later to access the geometry tools
         (symmetry finder ...) with the binding.
         The binding can be used by V_Sim (http://www-drfmc.cea.fr/sp2m/L_Sim/V_Sim/index.en.html)
         Work by Damien Caliste.

    B.14 The m-projected DOS is now available via the input variable prtdosm,
         and has been parallelised (including band/FFT), and made compatible with mkmem=0
         (among other improvements). Tests to be provided.
         Work by M. Torrent

    B.15 The algorithm of D. Alfe (Comp. Phys. Comm. 118, 31 (1999)) for the
         prediction of the density, especially useful for the parallel case,
         has been implemented. It is available with iprcch 5 or 6.
         Tests to be provided.
         Work by M. Torrent

    __________________________________________________________________

    C. Changes for the developers

    C.1 The PAW part of ABINIT has been largely rewritten :
        creation of methods to manipulate the objects, routines to modularise the
        operations, etc. This was needed to gather the new developements
        related to GW as well as response-function.
        Some speed-up of PAW calculations thanks to the storage in memory
        of the coefficients of the scalar product between wavefunctions and
        projectors (cprj). Input variable pawusecp. Default value is such
        that the cprj coefficients are now always stored
        Also, new DEV variables for PAW : e.g. pawnhatxc .
        Work by M. Torrent and M. Giantomassi.

    C.2 The robustness of the build system (bug fixing) has been
        improved a lot by Y. Pouillon . Also, several problems have been
        observed, isolated, fixed by Y. Pouillon, J.-M. Beuken, D. Caliste and A. Jacques.

    C.3 Preparation of a testing farm (platforms, compilers) by JM Beuken.
        At present, the machine chum has been customized with 6 compilers.
        New config/build-examples/*cf files have been prepared.

    C.4 The GW part of the code can now use double precision complex
        or single precision complex. This is to allow better portability and testing.
        At present, only accessible through modification by hand of the
        meaning of gwp and gwpc in src/defs_basis.F90 .
        Work by Matteo Giantomassi.

        abiruls.pl has been updated to accomodate this modification (L. Sindic)

    C.5 Cleaning : lib has been replaced by plugins ; lib/light
        has been suppressed ; lib/numeric transferred to src
    __________________________________________________________________

    D.  Other changes
    (or on-going developments, not yet finalized).

    D.1 PAW + Response function  (phonons) is under development, by M. Torrent.

    D.2 For GW calculations, more choice is offered concerning
        memory management. Input variable gwmem.
        Test v5#65. Work done by Matteo Giantomassi.

    D.3 In order to restart calculations of Sigma, one will be able to rely on the susceptibility,
        thanks to the input variables irdsuscep and getsuscep .
        Doc and automatic tests missing.
        Work in progress, by Matteo Giantomassi.

    D.4 The input variable prtgkk can be used to ease the preparation of the
        computation of electron-phonon matrix elements.
        See automatic test v5#85,86,87,88.
        Associated input variable prepgkk.
        Work done by Matthieu Verstraete.

    D.5 The response to a magnetic field (with PAW) is under development.
        Work by Joe Zwanziger. Input variable rfmgfd.

    D.6 The ability to define the number of bands per processor for the lobpcg
        algorithm has been coded. Input variable bandpp.
        No automatic test available yet.
        Work by Francois Bottin.

    D.7 New test v5#67, for the input variable soenergy, namely the
        scissor operator possibly used in GW calculations.
        Contributed by Riad Shaltaf.

    D.8 New test paral#U, to test the parallelism of GW calculations without
        plasmon pole models. Contributed by F. Bruneval.

    D.9 Cut3D output of wavefunction : option 14, to output wfs using the Gaussian/cube format.
        Work by F. Bruneval.

    D.10 New input variable prtspcur, only for development, by MVerstraete.
        Associated tests v5#95,96,97,98

    D.11 New test v5#91, for the development of the lobpcgII algorithm.
         Work by S. Leroux and G. Zerah .

    D.12 Maintenance work on abilint by T. Deutsch

    D.13 Decomposition of the frequency-dependent response of a slab
         into a bulk contribution and a surface contribution.
         Work by C. Motta. Tests and documentation to be provided.

    D.14 For spin polarized systems, density and magnetization are mixed,
         instead of up and down density.
         Work by M. Torrent

    D.15 Work on GW for spin-polarized and spinorial cases.
         In particular, for spin-polarized case, COH and contour integration method are treated.
         Work by M. Giantomassi

    D.16 Work on parallelism of GW, by F. Bruneval.

    D.17 Work on LDA + DMFT, by B. Amadon.

    D.18 Allowing GW for charged system, by F. Bruneval. Test to be provided ?!

    D.19 In addition, miscellaneous bug fixes by many developers (see the list of contributors),
         many bug fixes in documentation, improvements in the tutorials ...

    D.20 Allowed to use '*' in the input file, without specifying the number of times the value has to be repeated.