&CONTROL                                                                      
                     title = 'LaMnO3-Pero',                                    
               calculation = 'vc-relax' ,                                      
               max_seconds = 84000 ,                                           
              restart_mode = 'restart' ,                                  
                    outdir = './OUT',
                pseudo_dir = '../',               
                    prefix = 'LaMnO3' ,
             etot_conv_thr = 1.0d-5 ,                                          
             forc_conv_thr = 1.0d-4 ,                                          
                   tstress = .true. ,                                          
                   tprnfor = .true. ,
		        dt = 80                              
 /                                                                             
 &SYSTEM                                                                       
                     ibrav =  8 ,
                 celldm(1) =   10.858559 ,
                 celldm(2) =   1.33372100, 
                 celldm(3) =   0.9629650 ,
                       nat = 20 ,                                              
                      ntyp = 4 ,                                               
                   ecutwfc = 25 ,                                              
                   ecutrho = 250 ,                                             
                     nspin = 2,
    starting_magnetization(1)=0.0,    starting_magnetization(2)=1.0,
    starting_magnetization(3)=-1.0,   starting_magnetization(4)=0.0,
               occupations = 'smearing' ,
                  smearing = 'methfessel-paxton',                              
                 degauss = 0.01                                              
 /                                                                             
 &ELECTRONS                                                                    
               mixing_beta = 0.5 ,                                             
           diagonalization = 'david' ,                                         
            diago_thr_init = 1.D-2 ,                                           
          diago_david_ndim = 10                                                
 /                                                                             
 &IONS                                                                         
      ion_dynamics = 'damp',
 / 
 &CELL
             cell_dynamics = 'damp-w',
                     press = 0.0D0 ,
             press_conv_thr=0.5D0,
	             wmass = 0.01
/                                                                             
ATOMIC_SPECIES                                                                 
  La   138.9055    La.pbe-nsp-van.UPF
  Mn1   54.938049  Mn.pbe-sp-van.UPF
  Mn2   54.938049  Mn.pbe-sp-van.UPF
  O     15.9994    O.pbe-van_ak.UPF
ATOMIC_POSITIONS (crystal)                                                     
La       0.0513  0.2500  0.9814
La       0.4487  0.7500  1.4814 
La      -0.0513  0.7500 -0.9814
La       0.5513  0.2500 -0.4814
Mn1      0.0000  0.0000  0.5000 
Mn1      0.5000  0.0000  0.0000
Mn2      0.0000  0.5000  0.5000
Mn2      0.5000  0.5000  0.0000
O        0.4849  0.2500  0.0778
O        0.0151  0.7500  0.5778
O       -0.4849  0.7500 -0.0778 
O        0.9849  0.2500  0.4222
O        0.3085  0.0408  0.7227
O        0.1915 -0.0408  1.2227
O       -0.3085  0.5408 -0.7227
O        0.8085  0.4592 -0.2227 
O       -0.3085 -0.0408 -0.7227
O        0.8085  0.0408 -0.2227 
O        0.3085  0.4592  0.7227
O        0.1915  0.5408  1.2227
K_POINTS automatic                                                             
 4  4  4  1 1 1