<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

* Test 8 concerns the LDA+U method (for NiO)

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

* Tests 20-39 concern response functions

* Tests 40-59 concern ground state and geometry optimization

* Tests 60-... concern electronic excitations
------------------------------------------------------------------
***************

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)
     actualy using "userid" as a keyword
     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)
     actualy using "userid" as a keyword
     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)   
     
 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

 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 
     core-level shift of level 1s, thanks to the
     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.

 43. 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

 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

</pre>

