[Pw_forum] esm and phonons

Angelo Peronio peronio at tasc.infm.it
Wed Jun 22 16:40:26 CEST 2011


Dear QE users,

is it possible to perform a phonon calculation after a SCF one which
employed the Effective Screening Medium by Otani and Sugino
(assume_isolated = 'esm')?

I tried (NO molecule between two electrodes) and did not succeed: ph.x
outputs
[...]
     Representation #  1 mode #   1

     Self-consistent Calculation
     kpoint   1 ibnd   8 solve_linter: root not converged       NaNQ
 
     Pert. #  1: Fermi energy shift (Ry) =           NaNQ           NaNQ

      iter #   1 total cpu time :    22.5 secs   av.it.: 200.0
      thresh= 0.100E-01 alpha_mix =  0.300 |ddv_scf|^2 =       NaNQ
[...]

and so on, up to a crash: "from broyden : error # 1    factorization".

Attached you can find pw.x and ph.x inputs and outputs, the run was done
with the QE 4.3.1 pre-built on CINECA's SP6 machine.

This was only a preliminary test: what I am trying to do is to study the
effect of an electric field on the vibrational modes of molecule adsorbed
on a metal surface. I already did it with the "tefield" implementation of
electric field and since the effect is tiny I would like to have a
consistency check.

Thanks in advance for any help, comment or suggestion.
Angelo
-- 
_____________________
Angelo Peronio
IOM CNR laboratorio TASC
AREA Science Park, edificio MM
strada statale 14, km 163.5
34149 Basovizza (TS)
Italy

laboratory: +39 040 3756452
office:     +39 040 3756458
fax:        +39 040 226767
e-mail:     peronio at tasc.infm.it
-------------- next part --------------
&CONTROL
    calculation = 'relax',
    prefix      = 'no-gas-ph-esm-02',
    pseudo_dir  = '/gpfs/scratch/userexternal/aperonio/pseudo',
    outdir      = '/gpfs/scratch/userexternal/aperonio/outdirs/no-gas-ph-esm-02',
    disk_io     = 'low',
    wf_collect  = .false.,
!    max_seconds = 120,
!    restart_mode = 'restart',
/

&SYSTEM
    ibrav   = 1,
    celldm(1) = 20,
    nat     = 2,
    ntyp    = 2,
    ecutwfc = 30,
    ecutrho = 240,
    occupations = 'smearing',
    smearing = 'marzari-vanderbilt',
    degauss = 0.02,
    nspin = 2,
    starting_magnetization(1)=0.7,
    starting_magnetization(2)=0.3,
    assume_isolated = 'esm',
    esm_bc = 'pbc',
!    esm_w =
    esm_efield = 0.0d0,
!    esm_nfit = 
/

&ELECTRONS
   conv_thr = 1.0d-10,
/

&IONS
/

ATOMIC_SPECIES
N     14.0067   N.pbe-van_bm.UPF
O     15.9994   O.pbe-van_bm.UPF

ATOMIC_POSITIONS {bohr}
N       0.0   0.0   -1.1013d0        0 0 0
O       0.0   0.0    1.1013d0

K_POINTS {automatic}
1 1 1   0 0 0
-------------- next part --------------

     Program PWSCF v.4.3.1      starts on 22Jun2011 at 15:42:26 

     This program is part of the open-source Quantum ESPRESSO suite
     for quantum simulation of materials; please cite
         "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);
          URL http://www.quantum-espresso.org", 
     in publications or presentations arising from this work. More details at
     http://www.quantum-espresso.org/wiki/index.php/Citing_Quantum-ESPRESSO

     Parallel version (MPI), running on    16 processors
     K-points division:     npool     =    2
     R & G space division:  proc/pool =    8

     Current dimensions of program PWSCF are:
     Max number of different atomic species (ntypx) = 10
     Max number of k-points (npk) =  40000
     Max angular momentum in pseudopotentials (lmaxx) =  3
     Waiting for input...
          Reading PseudoPotential from 
               /gpfs/scratch/userexternal/aperonio/pseudo/N.pbe-van_bm.UPF
          Reading PseudoPotential from 
               /gpfs/scratch/userexternal/aperonio/pseudo/O.pbe-van_bm.UPF

     Subspace diagonalization in iterative solution of the eigenvalue problem:
     parallel, distributed-memory algorithm (size of sub-group:  2*  2 procs)

 
   Stick Mesh
   ----------
   nst =  7645,  nstw =   949, nsts =  3809
               n.st   n.stw   n.sts    n.g    n.gw   n.gs
   min         955     118     476   62807    2762   22180
   max         956     119     477   62812    2767   22202
       7645     949    3809  502479   22119  177509
 


     bravais-lattice index     =            1
     lattice parameter (alat)  =      20.0000  a.u.
     unit-cell volume          =    8000.0000 (a.u.)^3
     number of atoms/cell      =            2
     number of atomic types    =            2
     number of electrons       =        11.00
     number of Kohn-Sham states=           10
     kinetic-energy cutoff     =      30.0000  Ry
     charge density cutoff     =     240.0000  Ry
     convergence threshold     =      1.0E-10
     mixing beta               =       0.7000
     number of iterations used =            8  plain     mixing
     Exchange-correlation      =  SLA  PW   PBE  PBE ( 1 4 3 4 0)
     EXX-fraction              =        0.00
     nstep                     =           50


     Effective Screening Medium Method
     =================================
     field strength (Ry/a.u.)         =       0.00 
     ESM offset from cell edge (a.u.) =       0.00 
     grid points for fit at edges     =          4 
     Ordinary Periodic Boundary Conditions
 

     celldm(1)=  20.000000  celldm(2)=   0.000000  celldm(3)=   0.000000
     celldm(4)=   0.000000  celldm(5)=   0.000000  celldm(6)=   0.000000

     crystal axes: (cart. coord. in units of alat)
               a(1) = (   1.000000   0.000000   0.000000 )  
               a(2) = (   0.000000   1.000000   0.000000 )  
               a(3) = (   0.000000   0.000000   1.000000 )  

     reciprocal axes: (cart. coord. in units 2 pi/alat)
               b(1) = (  1.000000  0.000000  0.000000 )  
               b(2) = (  0.000000  1.000000  0.000000 )  
               b(3) = (  0.000000  0.000000  1.000000 )  


     PseudoPot. # 1 for N  read from file N.pbe-van_bm.UPF
     MD5 check sum: 90d10febe142ad51ac9d5be1deb811c4
     Pseudo is Ultrasoft, Zval =  5.0
     Generated by new atomic code, or converted to UPF format
     Using radial grid of  729 points,  4 beta functions with: 
                l(1) =   0
                l(2) =   0
                l(3) =   1
                l(4) =   1
     Q(r) pseudized with  8 coefficients,  rinner =    0.800   0.800   0.800


     PseudoPot. # 2 for O  read from file O.pbe-van_bm.UPF
     MD5 check sum: ca2828d9fb1f294809e918903b7648d1
     Pseudo is Ultrasoft, Zval =  6.0
     Generated by new atomic code, or converted to UPF format
     Using radial grid of  737 points,  4 beta functions with: 
                l(1) =   0
                l(2) =   0
                l(3) =   1
                l(4) =   1
     Q(r) pseudized with  8 coefficients,  rinner =    0.800   0.800   0.800


     atomic species   valence    mass     pseudopotential
        N              5.00    14.00670     N ( 1.00)
        O              6.00    15.99940     O ( 1.00)

     Starting magnetic structure 
     atomic species   magnetization
        N            0.700
        O            0.300

      8 Sym.Ops. (no inversion)


   Cartesian axes

     site n.     atom                  positions (alat units)
         1           N   tau(  1) = (   0.0000000   0.0000000  -0.0550650  )
         2           O   tau(  2) = (   0.0000000   0.0000000   0.0550650  )

     number of k points=     2  Marzari-Vanderbilt smearing, width (Ry)=  0.0200
                       cart. coord. in units 2pi/alat
        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   1.0000000
        k(    2) = (   0.0000000   0.0000000   0.0000000), wk =   1.0000000

     G cutoff = 2431.7084  ( 502479 G-vectors)     FFT grid: (110,110,110)
     G cutoff = 1215.8542  ( 177509 G-vectors)  smooth grid: ( 70, 70, 70)

     Largest allocated arrays     est. size (Mb)     dimensions
        Kohn-Sham Wavefunctions         0.42 Mb     (   2767,  10)
        NL pseudopotentials             0.68 Mb     (   2767,  16)
        Each V/rho on FFT grid          5.17 Mb     ( 169400,   2)
        Each G-vector array             0.48 Mb     (  62812)
        G-vector shells                 0.02 Mb     (   2016)
     Largest temporary arrays     est. size (Mb)     dimensions
        Auxiliary wavefunctions         1.69 Mb     (   2767,  40)
        Each subspace H/S matrix        0.02 Mb     (     40,  40)
        Each <psi_i|beta_j> matrix      0.00 Mb     (     16,  10)
        Arrays for rho mixing          20.68 Mb     ( 169400,   8)

     Initial potential from superposition of free atoms
     Check: negative starting charge=(component1):   -0.002084
     Check: negative starting charge=(component2):   -0.001185

     starting charge   10.99999, renormalised to   11.00000

     negative rho (up, down):  0.208E-02 0.118E-02
     Starting wfc are    8 atomic +    2 random wfc

     total cpu time spent up to now is      2.39 secs

     per-process dynamical memory:   197.4 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations =  6.0

     negative rho (up, down):  0.252E-02 0.199E-02

     total cpu time spent up to now is      4.03 secs

     total energy              =     -51.65926069 Ry
     Harris-Foulkes estimate   =     -51.44041938 Ry
     estimated scf accuracy    <       0.21844163 Ry

     total magnetization       =     1.00 Bohr mag/cell
     absolute magnetization    =     1.14 Bohr mag/cell

     iteration #  2     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.99E-03,  avg # of iterations =  2.0

     negative rho (up, down):  0.305E-02 0.251E-02

     total cpu time spent up to now is      5.53 secs

     total energy              =     -51.65146146 Ry
     Harris-Foulkes estimate   =     -51.73300180 Ry
     estimated scf accuracy    <       0.21580913 Ry

     total magnetization       =     1.00 Bohr mag/cell
     absolute magnetization    =     1.05 Bohr mag/cell

     iteration #  3     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.96E-03,  avg # of iterations =  2.0

     negative rho (up, down):  0.276E-02 0.253E-02

     total cpu time spent up to now is      7.02 secs

     total energy              =     -51.71705264 Ry
     Harris-Foulkes estimate   =     -51.72770453 Ry
     estimated scf accuracy    <       0.07169371 Ry

     total magnetization       =     1.00 Bohr mag/cell
     absolute magnetization    =     1.05 Bohr mag/cell

     iteration #  4     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  6.52E-04,  avg # of iterations =  1.5

     negative rho (up, down):  0.274E-02 0.265E-02

     total cpu time spent up to now is      8.51 secs

     total energy              =     -51.72855187 Ry
     Harris-Foulkes estimate   =     -51.72799460 Ry
     estimated scf accuracy    <       0.00069238 Ry

     total magnetization       =     1.00 Bohr mag/cell
     absolute magnetization    =     1.03 Bohr mag/cell

     iteration #  5     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  6.29E-06,  avg # of iterations = 16.5

     negative rho (up, down):  0.270E-02 0.269E-02

     total cpu time spent up to now is     10.26 secs

     total energy              =     -51.72905395 Ry
     Harris-Foulkes estimate   =     -51.72895663 Ry
     estimated scf accuracy    <       0.00010165 Ry

     total magnetization       =     1.00 Bohr mag/cell
     absolute magnetization    =     1.03 Bohr mag/cell

     iteration #  6     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  9.24E-07,  avg # of iterations = 13.5

     negative rho (up, down):  0.268E-02 0.271E-02

     total cpu time spent up to now is     11.98 secs

     total energy              =     -51.72910574 Ry
     Harris-Foulkes estimate   =     -51.72909107 Ry
     estimated scf accuracy    <       0.00001217 Ry

     total magnetization       =     1.00 Bohr mag/cell
     absolute magnetization    =     1.03 Bohr mag/cell

     iteration #  7     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.11E-07,  avg # of iterations =  2.5

     negative rho (up, down):  0.268E-02 0.271E-02

     total cpu time spent up to now is     13.56 secs

     total energy              =     -51.72911799 Ry
     Harris-Foulkes estimate   =     -51.72911738 Ry
     estimated scf accuracy    <       0.00000073 Ry

     total magnetization       =     1.00 Bohr mag/cell
     absolute magnetization    =     1.03 Bohr mag/cell

     iteration #  8     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  6.60E-09,  avg # of iterations =  2.5

     negative rho (up, down):  0.268E-02 0.271E-02

     total cpu time spent up to now is     15.11 secs

     total energy              =     -51.72911884 Ry
     Harris-Foulkes estimate   =     -51.72911878 Ry
     estimated scf accuracy    <       0.00000003 Ry

     total magnetization       =     1.00 Bohr mag/cell
     absolute magnetization    =     1.03 Bohr mag/cell

     iteration #  9     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.72E-10,  avg # of iterations =  2.5

     negative rho (up, down):  0.268E-02 0.271E-02

     total cpu time spent up to now is     16.68 secs

     total energy              =     -51.72911887 Ry
     Harris-Foulkes estimate   =     -51.72911887 Ry
     estimated scf accuracy    <          6.6E-10 Ry

     total magnetization       =     1.00 Bohr mag/cell
     absolute magnetization    =     1.03 Bohr mag/cell

     iteration # 10     ecut=    30.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  5.97E-12,  avg # of iterations =  2.0

     negative rho (up, down):  0.268E-02 0.271E-02

     total cpu time spent up to now is     18.03 secs

     End of self-consistent calculation

 ------ SPIN UP ------------


          k = 0.0000 0.0000 0.0000 ( 22119 PWs)   bands (ev):

   -31.2947 -16.8855 -12.6879 -12.6879 -11.5761  -4.4472  -4.4472  -0.4962
     0.7319   0.7723

 ------ SPIN DOWN ----------


          k = 0.0000 0.0000 0.0000 ( 22119 PWs)   bands (ev):

   -30.8017 -16.0562 -11.7974 -11.7974 -10.9813  -3.0901  -3.0901  -0.4338
     0.7944   0.8554

     the Fermi energy is    -4.4093 ev

!    total energy              =     -51.72911887 Ry
     Harris-Foulkes estimate   =     -51.72911887 Ry
     estimated scf accuracy    <          3.1E-11 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =    -103.00652293 Ry
     hartree contribution      =      52.78098676 Ry
     xc contribution           =     -11.64948453 Ry
     ewald contribution        =      10.15246548 Ry
     smearing contrib. (-TS)   =      -0.00656366 Ry

     total magnetization       =     1.00 Bohr mag/cell
     absolute magnetization    =     1.03 Bohr mag/cell

     ESM Charge and Potential
     ========================

    z (A)      rho (e)      Avg v_hartree        Avg v_local  Avg v_hart+v_loc
                                   (eV/A)             (eV/A)            (eV/A)
    ==========================================================================
   -5.196     -0.00002        -13.9264976         14.2712723         0.3447748
   -5.099      0.00002        -13.8974847         14.2433913         0.3459066
   -5.003     -0.00002        -13.8391001         14.1861335         0.3470333
   -4.907      0.00002        -13.7513427         14.0995026         0.3481599
   -4.811      0.00000        -13.6342133         13.9835046         0.3492912
   -4.714      0.00001        -13.4877116         13.8381248         0.3504132
   -4.618      0.00001        -13.3118376         13.6633852         0.3515476
   -4.522      0.00000        -13.1065917         13.4592583         0.3526666
   -4.426      0.00003        -12.8719735         13.2257750         0.3538015
   -4.330     -0.00001        -12.6079840         12.9629033         0.3549193
   -4.233      0.00004        -12.3146223         12.6706738         0.3560516
   -4.137      0.00001        -11.9918898         12.3490595         0.3571696
   -4.041      0.00006        -11.6397859         11.9980820         0.3582961
   -3.945      0.00004        -11.2583121         11.6177264         0.3594143
   -3.849      0.00008        -10.8474683         11.2080000         0.3605316
   -3.752      0.00010        -10.4072561         10.7689034         0.3616473
   -3.656      0.00013         -9.9376768         10.3004282         0.3627514
   -3.560      0.00020         -9.4387319          9.8025901         0.3638582
   -3.464      0.00023         -8.9104248          9.2753669         0.3649421
   -3.367      0.00038         -8.3527585          8.7187860         0.3660276
   -3.271      0.00046         -7.7657393          8.1328164         0.3670771
   -3.175      0.00071         -7.1493735          7.5174909         0.3681175
   -3.079      0.00093         -6.5036727          6.8727765         0.3691038
   -2.983      0.00135         -5.8286503          6.1987049         0.3700546
   -2.886      0.00188         -5.1243279          5.4952473         0.3709193
   -2.790      0.00266         -4.3907346          4.7624278         0.3716932
   -2.694      0.00377         -3.6279124          4.0002286         0.3723161
   -2.598      0.00529         -2.8359209          3.2086597         0.3727388
   -2.502      0.00756         -2.0148430          2.3877204         0.3728774
   -2.405      0.01067         -1.1647998          1.5374006         0.3726009
   -2.309      0.01530         -0.2859609          0.6577228         0.3717620
   -2.213      0.02172          0.6214258         -0.2513497         0.3700761
   -2.117      0.03105          1.5570118         -1.1897644         0.3672474
   -2.020      0.04421          2.5202955         -2.1575914         0.3627041
   -1.924      0.06331          3.5105641         -3.1547418         0.3558223
   -1.828      0.09084          4.5267870         -4.1813248         0.3454622
   -1.732      0.13041          5.5674778         -5.2372107         0.3302672
   -1.636      0.18694          6.6305011         -6.3225495         0.3079516
   -1.539      0.26664          7.7128117         -7.4371723         0.2756394
   -1.443      0.37971          8.8101061         -8.5812648         0.2288413
   -1.347      0.53804          9.9162954         -9.7546247         0.1616708
   -1.251      0.75679         11.0228510        -10.9574670         0.0653841
   -1.155      1.03962         12.1180768        -12.1894531        -0.0713763
   -1.058      1.38402         13.1868805        -13.4524122        -0.2655317
   -0.962      1.77275         14.2109278        -14.7561650        -0.5452372
   -0.866      2.12333         15.1698382        -16.0731229        -0.9032847
   -0.770      2.33264         16.0453749        -17.3176295        -1.2722546
   -0.673      2.37679         16.8264704        -18.3889554        -1.5624850
   -0.577      2.32638         17.5104278        -19.2037165        -1.6932887
   -0.481      2.25285         18.0995999        -19.7309585        -1.6313585
   -0.385      2.15719         18.5979300        -20.0092095        -1.4112795
   -0.289      2.03191         19.0105744        -20.1254095        -1.1148351
   -0.192      1.91772         19.3440342        -20.1806255        -0.8365912
   -0.096      1.87358         19.6041279        -20.2555020        -0.6513741
    0.000   -287.51570         19.7931416        -20.3708587        -0.5777171
    0.096   -286.04631         19.9088070        -20.5172867        -0.6084797
    0.192   -282.40604         19.9441223        -20.7063017        -0.7621794
    0.289   -276.81219         19.8861134        -20.9127531        -1.0266398
    0.385   -269.61565         19.7166175        -21.0431534        -1.3265359
    0.481   -261.15867         19.4166405        -20.9730835        -1.5564430
    0.577   -251.64729         18.9712676        -20.5956396        -1.6243720
    0.673   -241.16305         18.3727057        -19.8673076        -1.4946019
    0.770   -229.75048         17.6226661        -18.8227858        -1.2001197
    0.866   -217.37361         16.7353553        -17.5637779        -0.8284226
    0.962   -203.86939         15.7371031        -16.2150613        -0.4779582
    1.058   -189.03428         14.6601590        -14.8763509        -0.2161920
    1.155   -175.48462         13.5344885        -13.5804525        -0.0459640
    1.251   -635.36735         12.3835347        -12.3161422         0.0673925
    1.347      0.45944         11.2246977        -11.0808868         0.1438109
    1.443      0.30658         10.0703485         -9.8752192         0.1951293
    1.539      0.20477          8.9287758         -8.6987326         0.2300432
    1.636      0.13743          7.8054861         -7.5518013         0.2536848
    1.732      0.09166          6.7041357         -6.4340751         0.2700607
    1.828      0.06126          5.6271936         -5.3458742         0.2813194
    1.924      0.04093          4.5763084         -4.2869143         0.2893941
    2.020      0.02758          3.5525788         -3.2574376         0.2951411
    2.117      0.01872          2.5567291         -2.2572474         0.2994816
    2.213      0.01265          1.5892399         -1.2864924         0.3027474
    2.309      0.00865          0.6504379         -0.3450711         0.3053668
    2.405      0.00577         -0.2594581          0.5669603         0.3075022
    2.502      0.00402         -1.1402955          1.4496171         0.3093216
    2.598      0.00267         -1.9919762          2.3029196         0.3109434
    2.694      0.00192         -2.8144302          3.1268188         0.3123886
    2.790      0.00127         -3.6076143          3.9213848         0.3137705
    2.886      0.00091         -4.3714954          4.6865343         0.3150390
    2.983      0.00062         -5.1060528          5.4223560         0.3163032
    3.079      0.00042         -5.8112711          6.1287631         0.3174919
    3.175      0.00032         -6.4871400          6.8058341         0.3186941
    3.271      0.00019         -7.1336534          7.4535039         0.3198504
    3.367      0.00018         -7.7508052          8.0718206         0.3210154
    3.464      0.00007         -8.3385935          8.6607551         0.3221616
    3.560      0.00011         -8.8970144          9.2203171         0.3233027
    3.656      0.00003         -9.4260677          9.7505154         0.3244478
    3.752      0.00006         -9.9257508         10.2513250         0.3255742
    3.849      0.00002        -10.3960640         10.7227839         0.3267198
    3.945      0.00002        -10.8370061         11.1648450         0.3278389
    4.041      0.00002        -11.2485770         11.5775601         0.3289831
    4.137      0.00000        -11.6307764         11.9608773         0.3301009
    4.233      0.00003        -11.9836035         12.3148445         0.3312410
    4.330     -0.00002        -12.3070592         12.6394213         0.3323621
    4.426      0.00003        -12.6011421         12.9346378         0.3334957
    4.522     -0.00002        -12.8658534         13.2004760         0.3346226
    4.618      0.00002        -13.1011919         13.4369412         0.3357493
    4.714     -0.00001        -13.3071586         13.6440407         0.3368821
    4.811      0.00001        -13.4837526         13.8217561         0.3380035
    4.907      0.00000        -13.6309744         13.9701151         0.3391407
    5.003     -0.00001        -13.7488238         14.0890832         0.3402594
    5.099      0.00001        -13.8373007         14.1786982         0.3413975
    5.196     -0.00002        -13.8964055         14.2389220         0.3425165
    5.292      0.00002        -13.9261374         14.2697901         0.3436527
 

     convergence has been achieved in  10 iterations

     Forces acting on atoms (Ry/au):

     atom   1 type  1   force =     0.00000000    0.00000000   -0.00086155
     atom   2 type  2   force =     0.00000000    0.00000000    0.00086155

     Total force =     0.000862     Total SCF correction =     0.000012

     BFGS Geometry Optimization

     bfgs converged in   1 scf cycles and   0 bfgs steps
     (criteria: energy < 0.10E-03, force < 0.10E-02)

     End of BFGS Geometry Optimization

     Final energy   =     -51.7291188729 Ry
Begin final coordinates

ATOMIC_POSITIONS (bohr)
N        0.000000000   0.000000000  -1.101300000    0   0   0
O        0.000000000   0.000000000   1.101300000
End final coordinates



     Writing output data file no-gas-ph-esm-02.save
 
     init_run     :      0.98s CPU      2.00s WALL (       1 calls)
     electrons    :      7.53s CPU     15.84s WALL (       1 calls)
     forces       :      0.32s CPU      0.64s WALL (       1 calls)

     Called by init_run:
     wfcinit      :      0.03s CPU      0.09s WALL (       1 calls)
     potinit      :      0.41s CPU      0.82s WALL (       1 calls)

     Called by electrons:
     c_bands      :      0.96s CPU      1.99s WALL (      10 calls)
     sum_band     :      1.90s CPU      3.84s WALL (      10 calls)
     v_of_rho     :      3.07s CPU      6.23s WALL (      11 calls)
     newd         :      1.13s CPU      2.29s WALL (      11 calls)
     mix_rho      :      0.34s CPU      0.69s WALL (      10 calls)

     Called by c_bands:
     init_us_2    :      0.02s CPU      0.05s WALL (      21 calls)
     cegterg      :      0.90s CPU      1.87s WALL (      10 calls)

     Called by *egterg:
     h_psi        :      0.73s CPU      1.46s WALL (      65 calls)
     s_psi        :      0.02s CPU      0.03s WALL (      65 calls)
     g_psi        :      0.00s CPU      0.02s WALL (      54 calls)
     cdiaghg      :      0.09s CPU      0.18s WALL (      64 calls)

     Called by h_psi:
     add_vuspsi   :      0.01s CPU      0.03s WALL (      65 calls)

     General routines
     calbec       :      0.00s CPU      0.04s WALL (      76 calls)
     fft          :      1.99s CPU      3.99s WALL (     335 calls)
     ffts         :      0.03s CPU      0.12s WALL (      42 calls)
     fftw         :      0.67s CPU      1.30s WALL (     746 calls)
     interpolate  :      0.34s CPU      0.71s WALL (      42 calls)
     davcio       :      0.00s CPU      0.00s WALL (       1 calls)
 
     Parallel routines
     fft_scatter  :      0.96s CPU      1.95s WALL (    1123 calls)
 
     PWSCF        :     8.97s CPU        19.08s WALL

 
   This run was terminated on:  15:42:45  22Jun2011            

=------------------------------------------------------------------------------=
   JOB DONE.
=------------------------------------------------------------------------------=
-------------- next part --------------
NO gas phase phonons at Gamma.
&inputph
  prefix='no-gas-ph-esm-02',
  outdir='/gpfs/scratch/userexternal/aperonio/outdirs/no-gas-ph-esm-02',
  fildyn='no-gas-ph-esm-02.dynmat',
  alpha_mix(1) = 0.3,
  tr2_ph = 1d-14,
/
0.0 0.0 0.0
-------------- next part --------------

     Program PHONON v.4.3.1     starts on 22Jun2011 at 15:43:12 

     This program is part of the open-source Quantum ESPRESSO suite
     for quantum simulation of materials; please cite
         "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);
          URL http://www.quantum-espresso.org", 
     in publications or presentations arising from this work. More details at
     http://www.quantum-espresso.org/wiki/index.php/Citing_Quantum-ESPRESSO

     Parallel version (MPI), running on    16 processors
     K-points division:     npool     =    2
     R & G space division:  proc/pool =    8

     Ultrasoft (Vanderbilt) Pseudopotentials

   Info: using nr1, nr2, nr3 values from input

   Info: using nr1s, nr2s, nr3s values from input
          Reading PseudoPotential from 
               /gpfs/scratch/userexternal/aperonio/outdirs/no-gas-ph-esm-02/no-gas-ph-esm-02.save/N.pbe-van_bm.UPF
          Reading PseudoPotential from 
               /gpfs/scratch/userexternal/aperonio/outdirs/no-gas-ph-esm-02/no-gas-ph-esm-02.save/O.pbe-van_bm.UPF
 
   Stick Mesh
   ----------
   nst =  7645,  nstw =   949, nsts =  3809
               n.st   n.stw   n.sts    n.g    n.gw   n.gs
   min         955     118     476   62807    2762   22180
   max         956     119     477   62812    2767   22202
       7645     949    3809  502479   22119  177509
 

     negative rho (up, down):  0.268E-02 0.271E-02

     Calculation of q =    0.0000000   0.0000000   0.0000000

                                                                                

     bravais-lattice index     =            1
     lattice parameter (alat)  =      20.0000  a.u.
     unit-cell volume          =    8000.0000 (a.u.)^3
     number of atoms/cell      =            2
     number of atomic types    =            2
     kinetic-energy cut-off    =      30.0000  Ry
     charge density cut-off    =     240.0000  Ry
     convergence threshold     =      1.0E-14
     beta                      =       0.3000
     number of iterations used =            4
     Exchange-correlation      =  SLA  PW   PBE  PBE ( 1 4 3 4 0)
     EXX-fraction              =        0.00


     celldm(1)=   20.00000  celldm(2)=    0.00000  celldm(3)=    0.00000
     celldm(4)=    0.00000  celldm(5)=    0.00000  celldm(6)=    0.00000

     crystal axes: (cart. coord. in units of alat)
               a(1) = (  1.0000  0.0000  0.0000 )  
               a(2) = (  0.0000  1.0000  0.0000 )  
               a(3) = (  0.0000  0.0000  1.0000 )  

     reciprocal axes: (cart. coord. in units 2 pi/alat)
               b(1) = (  1.0000  0.0000  0.0000 )  
               b(2) = (  0.0000  1.0000  0.0000 )  
               b(3) = (  0.0000  0.0000  1.0000 )  


     Atoms inside the unit cell: 

   Cartesian axes

     site n.  atom      mass           positions (alat units)
        1        N   14.0067   tau( 1) = (    0.00000    0.00000   -0.05506  )
        2        O   15.9994   tau( 2) = (    0.00000    0.00000    0.05506  )

     Computing dynamical matrix for 
                    q = (   0.0000000   0.0000000   0.0000000 )
 
      9 Sym.Ops. (with q -> -q+G )


     G cutoff = 2431.7084  (  62812 G-vectors)     FFT grid: (110,110,110)
     G cutoff = 1215.8542  (  22190 G-vectors)  smooth grid: ( 70, 70, 70)

     number of k points=     2  Marzari-Vanderbilt smearing, width (Ry)=  0.0200
                       cart. coord. in units 2pi/alat
        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   1.0000000
        k(    2) = (   0.0000000   0.0000000   0.0000000), wk =   1.0000000

     PseudoPot. # 1 for N  read from file N.pbe-van_bm.UPF
     MD5 check sum: 643e1b9210c5596d3e8220b7f5116234
     Pseudo is Ultrasoft, Zval =  5.0
     Generated by new atomic code, or converted to UPF format
     Using radial grid of  729 points,  4 beta functions with: 
                l(1) =   0
                l(2) =   0
                l(3) =   1
                l(4) =   1
     Q(r) pseudized with  8 coefficients,  rinner =    0.800   0.800   0.800


     PseudoPot. # 2 for O  read from file O.pbe-van_bm.UPF
     MD5 check sum: 7ee51cec6b4d96081c02596cd828ca05
     Pseudo is Ultrasoft, Zval =  6.0
     Generated by new atomic code, or converted to UPF format
     Using radial grid of  737 points,  4 beta functions with: 
                l(1) =   0
                l(2) =   0
                l(3) =   1
                l(4) =   1
     Q(r) pseudized with  8 coefficients,  rinner =    0.800   0.800   0.800


     k=gamma and q=gamma tricks are used


     Atomic displacements:
     There are   6 irreducible representations

     Representation     1      1 modes - To be done

     Representation     2      1 modes - To be done

     Representation     3      1 modes - To be done

     Representation     4      1 modes - To be done

     Representation     5      1 modes - To be done

     Representation     6      1 modes - To be done



     Alpha used in Ewald sum =   2.5000
     PHONON       :     6.66s CPU        13.95s WALL



     Representation #  1 mode #   1

     Self-consistent Calculation
     kpoint   1 ibnd   8 solve_linter: root not converged       NaNQ
 
     Pert. #  1: Fermi energy shift (Ry) =           NaNQ           NaNQ

      iter #   1 total cpu time :    22.5 secs   av.it.: 200.0
      thresh= 0.100E-01 alpha_mix =  0.300 |ddv_scf|^2 =       NaNQ
     kpoint   1 ibnd   8 solve_linter: root not converged       NaNQ
 
     Pert. #  1: Fermi energy shift (Ry) =           NaNQ           NaNQ

      iter #   2 total cpu time :    31.1 secs   av.it.: 200.0
      thresh=      NaNQ alpha_mix =  0.300 |ddv_scf|^2 =       NaNQ
     kpoint   1 ibnd   8 solve_linter: root not converged       NaNQ
 
     Pert. #  1: Fermi energy shift (Ry) =           NaNQ           NaNQ

      iter #   3 total cpu time :    39.7 secs   av.it.: 200.0
      thresh=      NaNQ alpha_mix =  0.300 |ddv_scf|^2 =       NaNQ
     kpoint   1 ibnd   8 solve_linter: root not converged       NaNQ
 
     Pert. #  1: Fermi energy shift (Ry) =           NaNQ           NaNQ
-------------- next part --------------

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     task #         9
     from broyden : error #         1
     factorization
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


More information about the Pw_forum mailing list