<html>
<head>
<title>v5 : README</title>
</head>
<body>

<h5>Copyright (C) 2002-2007 ABINIT group (FJollet, MTorrent, GZerah, XGonze, DHamann, MVeithen)
<br> This file is distributed under the terms of the GNU General Public License, see
~abinit/COPYING or <a href="http://www.gnu.org/copyleft/gpl.txt">
http://www.gnu.org/copyleft/gpl.txt </a>.
<br> For the initials of contributors, see ~abinit/doc/developers/contributors.txt .
</h5>

<pre>

This directory, v5, contains tests which 
are related to the development of the various development in 
version 5.x.x of Abinit: New build system, new files format, 
real space wavelet calculation, real space preconditioner ...
This file gives first a very brief description
of the tests, then a much longer one.

A classification of the tests is provided now. Later, each test is described
in more detail, with accompanying notes and references (some of the notes
might have to be updated, due to the small possible drift related to improvement
of the accuracy of the code).

------------------------------------------------------------------
(the following list is to be clarified later)

* Tests 1-5 concern the SCF preconditioners

* Test 7 concerns the Poisson solver in real space, using wavelets

* Tests 8 and 11 concern the LDA+U method (for NiO)

* Tests 9-10 concern the calculation of the conductivity with PAW (for Al)

* Tests 20-29 concern response functions

* Tests 30-59 concern ground state and geometry optimization

* Tests 60-79 concern electronic excitations

* Tests 80-99 concern molecular dynamics
------------------------------------------------------------------
***************

To run these tests, do the following:

0. Be sure that the perl script "run-standard-tests", in the main directory,
   has been generated from the primitive file "run-standard-tests.pl",
   thanks to the command "make perl"
   issued in ~abinit (see the installation notes on the Web).

1. Submit the "run-standard-tests" script, specifying a machine, and
   the keyword "v5" (for Tests_v5) and either
   the index of a test case, i.e.
   (run-standard-tests name_of_machine v5 22) >& log_file
   or a whole range of test cases (two indices), i.e.
   (run-standard-tests name_of_machine v5 02 08) >& log_file
   or, if you want to run all the test cases of this directory
   (run-standard-tests name_of_machine v5) >& log_file
   This will send stdout and stderr to log_file.
   The script 'run-standard-tests', will create a subdirectory of ~abinit/tests/v5,
   with the name_of_machine and the
   date, where all the results will be placed.

2. In that directory, you will find for each test case that you have
   run, a log file (with the name of the test case), an output
   file, but also a 'diff.xxx' file, automatically created by making
   a 'diff' with respect to the "Refs" subdirectory output files.
   It contains output files from a recent version of the ABINIT code.
   There may be large differences in timing but there should only
   be minor differences in the output of physical quantities.

3. There is also a global report file, generated by the use of the
   fldiff script. Its name is fldiff.report . See the last
   version of the  ~abinit/doc/install_notes/install** file
   in the Infos directory for information about the use of this file.
   This file is the most convenient for a quick look at the correctness
   of results. When the results are not correct, one has often to
   rely on the 'diff.xxx' file to understand what was going wrong.


**********

Test cases:

  1. Simple test for TFvW preconditionner (not already effective)
     simple scf calculation
     following part of the code are tested:
     04rsprc/prctfw.F90
     defs/functofrho_tfw.F90
     01cg/*
     (from PMAnglade)

  2. Simple test for TFvW2 preconditionner (not already effective)
     simple scf calculation
     following part of the code are tested:
     04rsprc/prctfw2.F90
     defs/func_tfvw2.F90     
     01cg/*
     THIS TEST DOES NOT WORK AT PRESENT
     (from PMAnglade)

  3. Simple test for rsprc preconditionner (not already effective)
     simple scf calculation
     following part of the code are tested:
     04rsprc/prcrscgres.F90
     defs/func_rscgres.F90     
     01cg/*     
     (from PMAnglade)

  4. Simple test for rsprc preconditionner (not already effective)
     simple scf calculation
     following part of the code are tested:
     04rsprc/prcrscgres2.F90
     defs/func_rscgres2.F90     
     01cg/*          
     (from PMAnglade)
     
  7. Simple test for real space computation (using a wavelet based poisson solver), 
     no forces, just an H atom
     in an empty box. The cut-off is unrealistic. This test tests the
     following part of the code :
      * psp spline generation for real space ;
      * local part of potential from pseudo ;
      * ion-ion interaction computation( instead of Ewald) ;
      * Hartree potential with Poisson's solver.
     (from D. Caliste)
     
  8. Simple test for LDA+U implementation.
     no forces, 2 Ni and 2 O atoms in an fcc supercell.
     The cut-off is unrealistic. This test tests the
     following part of the code :
      * paw LDA+U (pawpupot, pawpuenergy, pawpuinit)
     (from B. Amadon)
     
  9. Test to prepare a conductivity calculation (test 10)
     Ground state calculation for 4 Al atoms
     The number of bands is unrealistic. This test produces files that are inputs for test 10. It tests the
     following part of the code :
     *optic_paw
     (from S. Mazevet) 
     
 10. Test a PAW conductivity calculation
     The code conducti calls conducti_paw
     (from S. Mazevet)   
  
 11. Simple test for LDA+U features.
     2 Ni and 2 O atoms in an fcc supercell.
     The cut-offs are unrealistic.
     The aim of this test is to verify:
      * Antiferromagnetism and PAW+U
      * Use of an initial and constant density matrix during the first SCF iterations
      * Computation of "local magnetic moment" with prtdensph keyword
     (from M. Torrent)
     
 21. PbFCl in the matlockite structure P4/nmm. 6 atoms per cell.
     Compute a set of matrix elements of the dynamical matrix,
     for q wavevector 0.5 0.5 0.5 .
     In this test case, the symmetry operations have an important
     non-symmorphic character. There was a bug in the treatment
     of this material, prior to v5.2.4 .

 22. SiO2 alpha-quartz with 9 atoms per cell. 
     Compute the atomic temperature factors. 
     Contributed by P. Boulanger

 23. AlAs in Zincblende structure.
     Test calculation of the Born effective charge and dielectric constant
     tensor at finite electric fields.
     The first dataset is for normal ground-
     state calculation at zero electric field. The second dataset is for
     ground-state calculation at finite electric field. The final dataset
     is the linear-response calculation at finite electric field.
     (contributed by Xinjie Wang)

 28. H2 molecule, using a local pseudopotential, and Wu-Cohen GGA.
     Similar to test v3#08
     Compute the interatomic force constant for the
     displacement of the second atom along x.
     RF calculation of the 2DTE gives 189.49410411 Ha
     while finite differences gives :
     with delta(xred)=0.0002, 189.4945677  Ha
     with delta(xred)=0.0001, 189.49422003  Ha .
     Combining the results gives 189.4941041 Ha. Perfect !

 31. Bismuth atom, isolated, in a supercell.
     Compares different spin treatments :
     - spin-polarized collinear (nspinor=1, nspden=2, nsppol=2, no spin-orbit)
     - spinor (nspinor=2, nspden=4, nsppol=1, no spin-orbit)
     Both with fixed occupation numbers and metallic occupations.
     Also, consider mkmem=0
     This test has been contributed by LEDiaz

 32. Bismuth atom, isolated, in a supercell.
     Analyse the collinear spin-density obtained in test 31 , using cut3d

 33. Bismuth atom, isolated, in a supercell.
     Analyse the non-collinear spin-density obtained in test 31 , using cut3d

 41. Al2 molecule in a big box. Treat 8 bands, occupied using occopt=4.
     Compare tolrff and toldff stopping criteria
     Contributed M. Verstraete

 42. NH3 and NH2-CH3 molecules , supercell geometry
     Show the computations needed to obtain the differential 
     pseudopotentials 7n.pspnc (just usual atomic configuration)
     and 7n.1s_psp.mod (with a hole in state 1s).
     The differential core-level shift is obtained by comparing
     difference in the total energies of the NH3 molecule with both pseudopotentials
     with the one of the NH2CH3 molecule : 
     0.41522 Ha for the NH3 molecule 
     0.41768 Ha for the NH2-CH3 molecule, making for a shift of 2.46 mHa=0.0669 eV
     due to the replacement of one H atom of the NH3 molecule by a CH3 group.
     Note that these are completely non converged values. Parameters
     like ecut and acell should be the subject of a convergence study.
     For applications of this technique, see
     Rignanese GM et al, Phys. Rev. Lett. 79, 5174 (1997),
     or Travaly Y et al, Phys. Rev. B 61, 7716 (2000). 
     This test has been contributed by GMRignanese.

 44. Bismuth atom, isolated, in a supercell.
     Try alchemical pseudopotentials with spin-orbit.
     Spinor calculation (nspinor=2, nspden=1, nsppol=1)
     Metallic occupations.
     Also compare the halved spin-orbit interaction alchemical mixing
     with a pseudopotential directly prepared with halved spin-orbit interaction.
     The agreement is perfect.

 48. H2 molecule, using a local pseudopotential,
     and the Wu-Cohen GGA.
     Similar to Test v3#48
     Compute the total energy and force as a
     displacement of the second atom along x.
     Direct calculation of the gradient with respect
     to reduced coordinates gives -3.469435369425Ha
     while finite differences gives :
     with delta(xred)=0.0002, -3.4694479488 Ha
     with delta(xred)=0.0001, -3.4694385143 Ha .
     Combining the results gives -3.4694353695 Ha.
     Excellent.

 49. Si2 diatomic molecule. Spin-polarized (ferromagnetic).
     Test the Wu-Cohen GGA functional.
     Similar to Test v3#49
     Direct calculation of the gradient with respect
     to reduced coordinates gives  0.344935232634 Ha
     while finite differences gives :
     with delta(xred)=0.002,  0.34493521   Ha

 51. Stishovite, high-pressure polymorph of SiO2. Computation of the electric field
     gradient at each atomic site. This tests the implementation of the PAW-based
     efg algorithm as outlined in Profeta, Mauri, and Pickard, J. Am. Chem. Soc.
     125, 541-548 (2003). Coded by J. Zwanziger and M. Torrent, test contributed
     by J. Zwanziger.

 54. Crystalline Silicon to test the CP2K pseudopotential implementation with
     the new pseudo type: pspcod 10 (same formalism as HGH pseudos, pspcod 3).
     Same input file as t35.in of tests/tutorial, the only difference is that
     in the test the used pseudo is taken from the CP2K library.
     (http://cvs.berlios.de/cgi-bin/viewcvs.cgi/cp2k/potentials/Goedecker/abinit/)
     The bandstructure is in agreement with that calculated in the tutorial
     example with the pspcod=1 pseudopotential.
     Tests provided by S. Caravati .

 55. Crystalline rhombohedral Antimonium. CP2K PBE pseudopotential.
     Optimization of the atomic positions with ionmov 2. In the input file
     several datasets are defined but actually only one of them is chosen for
     the calculation. The provided xred (one for each dataset) are already the
     equilibrium position, so the calculation should use only one optimization
     step per dataset.
     Tests provided by S. Caravati .

 56. Bulk jellium (the slab fills up the entire cell): computation of the total
     energy. Tested 3 Wigner-Seitz radii: 1.5, 2 (~aluminum), 4. For each
     Wigner-Seitz radius the energy per particle agrees with that analytically
     calculated using the Teter-Pade functional form.
     The number of k points has been much reduced, for the purpose of the automatic testing.
     Test provided by S. Caravati .

 57. Jellium slab (50 jellium/ 50 vacuum) with Wigner-Seitz radius=4.
     Computation of the total energy
     The number of k points has been much reduced
     Test provided by S. Caravati .

 58. Jellium thin slab with Wigner-Seitz radius=2.56 (~copper) + 1 Na adatom.
     Computation of the total energy and force on the Na atom. Its position
     with respect to the slab already is the equilibrium position so the
     calculation should use only one optimization step.
     Test provided by S. Caravati .

 59. Jellium thin slab with Wigner-Seitz radius=2.56 (~copper) + 1 Na adatom.
     Computation of the total energy at the Na equilibrium distance from the
     jellium edge and two other ones, slightly farther and slightly closer.
     The total energies should indicate as more stable the equilibrium distance.
     The number of k points has been much reduced
     Test provided by S. Caravati .

 60. Crystalline Silicon 
     Computation of the non-self-consistent kinetic energy at fixed density.

 61. N2 molecule non-spin-polarized, but treated with nsppol=2, in a big box.
     Similar to test 55 of tests/v3. There should be only
     small numerical differences with respect to datasets 1 and 2 of this test,
     due to the fact that the tetragonal symmetry was here slightly broken
     for sake of portability of the automatic test.
     Note however that the singlet and triplet states cannot be disentagled,
     so one gets a big list of excitations, the merge of singlet and triplet
     excitations.

 62. BeH spin-polarized, in a supercell
     Test the spin-polarized TDDFT implementation. 
     nband 20 , ecut 8 , acell 14 14 14 .
     For the first two excitations, one finds :
     2.40967E+00 eV , doubly degenerate  (Pi excitation)
     3.94165E+00 eV  (Sigma+ excitation)
     They come mainly from the 2nd occupied Kohn-Sham state, spin up.
     The third excitation, at 5.11592E+00 eV (Sigma+ excitation) comes mainly from 
     the 1st occupied Kohn-Sham state, spin down.
     Then, come excitations that are in the continuum, or close to it 
     (depending on the parameters of the calculation), so, harder to converge,
     In order to have a calculation converged
     at the level of 0.1 eV, one needs parameters like
     nband 40 , ecut 10 , acell 30.
     In that case, the three above-mentioned excitations become
     2.388 eV (Pi), 4.432 eV (Sigma+), 5.357 eV (Sigma+). 
     In this case, the latter is no more the third excitation, as other ones,
     coupled to the vacuum have a lower energy. 
     These numbers can be compared to the values given by 
     J. Guan, ME Casida, DR Salahub, J. of Mol. Structure - Theochem 527,
     229, Sp. Iss. SI (2000) :
     2.391 eV (1Pi), 4.593 eV (2Sigma+), 5.418 eV (5Sigma+),
     and to the experimental values (see the above paper)
     2.56 eV/2.484 eV , 5.51 eV , 6.71 eV.
     Test provided by D. Sangalli

 63. Bcc Hydrogen, r_s=2.6 a.u AF
     To check a GW+spin calculation with and without symmetries
     Test provided by Matteo Giantomassi

 64. Crystalline alpha-quartz
     Calculation of the GW corrections with and without non-symmorphic operations

 81. This is a test, in the norm conserving case, of the Molecular Dynamics flag
     ionmov==13. With the values of optcell=0, 1, 2 (isothermal (0) , isothermal/isenthalpic
     with homogeneous cell deformations (1) or full cell deformations (2).)
     This is the reversible integrator of Martyna, Tuckerman, Tobias and Klein.
     (Mol. Phys., 1996, p 1117).

</pre>

