<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Dear all<br>By using script in wiki I obtain U for isolated Ti atom. but the value of  <br>#!/bin/bash<br>PSEUDO_DIR=/home/kazempou/WORK_FHI/TiO2/u-linear-response-atom<br>PW_ROOT=/home/akazem/espresso-4.0.2/<br>TMP_DIR=/home/kazempou/tmp/Ucalc<br>TMP_DIR1=/home/kazempou/tmp/Ucalc1<br>mkdir -p $TMP_DIR<br>mkdir -p $TMP_DIR1<br><br>WORK=$( pwd )<br><br>rm -rf $TMP_DIR/*<br>rm -rf $TMP_DIR1/*<br><br>mkdir $WORK/results_Mn_ucalc<br>cd $WORK/results_Mn_ucalc<br># unperturbed self-consistent calculation<br>cat > Ti.scf.in << EOF<br> &control<br>    calculation = 'scf'<br>    restart_mode='from_scratch',<br>    prefix='Ti',<br>    pseudo_dir = '$PSEUDO_DIR',<br>    outdir='$TMP_DIR'<br> /<br> &system<br>    ibrav=  0, celldm(1)=20,
 nat=  1, ntyp= 1,<br>    ecutwfc = 45 , ecutrho=400 , <br>    lda_plus_u = .true.<br>    hubbard_u(1) = 1.D-20, <br>    starting_magnetization(1)= 1.0,<br>    occupations='smearing', smearing='gauss', degauss=0.01,<br>    nspin=2,<br> /<br> &electrons<br>    mixing_beta = 0.7<br>    conv_thr =  1.0d-6<br> /<br>CELL_PARAMETERS<br>1.0 0.0 0.0<br>0.0 1.0 0.0<br>0.0 0.0 1.0<br>ATOMIC_SPECIES<br>Ti   47 ti.optgga2.fhi.UPF<br>ATOMIC_POSITIONS {crystal}<br>Ti 0.5 0.5 0.5<br>K_POINTS {gamma}<br>EOF<br><br>/home/kazempou/Desktop/espresso-4.0.5/bin/pw.x < Ti.scf.in > Ti.scf.out <br><br><br>for a in 0.0 -0.1 0.1  -0.2 0.2   ; do<br>rm -rf $TMP_DIR/*<br><br># perturbed calculations (finite alpha)<br>rm -f Ti.scf.in<br>cat > Ti.scf.in <<
 EOF<br> &control<br>    calculation = 'scf'<br>    restart_mode='from_scratch',<br>    prefix='Ti',<br>    pseudo_dir = '$PSEUDO_DIR',<br>    outdir='$TMP_DIR'<br> /<br> &system<br>    ibrav=  0, celldm(1)=20, nat=  1,  ntyp= 1,<br>    ecutwfc = 45.0 , ecutrho=400 , <br>    lda_plus_u = .true.<br>    hubbard_u(1) = 1.D-20, <br>    starting_magnetization(1)= 1.0,<br>    occupations='smearing', smearing='gauss', degauss=0.01,<br>    nspin=2,<br>    hubbard_alpha(1)=$a<br> /<br> &electrons<br>    diagonalization = 'cg'<br>    mixing_beta = 0.7<br>    conv_thr =  1.0d-6<br>    startingpot = 'file'<br>    startingwfc =
 'file'<br>    diago_thr_init = 1.d-12<br> /<br>CELL_PARAMETERS<br>1.0 0.0 0.0<br>0.0 1.0 0.0<br>0.0 0.0 1.0<br>ATOMIC_SPECIES<br>Ti   47  ti.optgga2.fhi.UPF<br>ATOMIC_POSITIONS {crystal}<br>Ti 0.5 0.5 0.5<br>K_POINTS {gamma}<br>EOF<br><br>/home/kazempou/Desktop/espresso-4.0.5/bin/pw.x < Ti.scf.in > Ti.pert_$a.out<br><br>done<br><br><br><br>thanks  a lot<br><br><br><br>Ali Kazempour<br>
Physics department, Isfahan University of Technology<br>
84156 Isfahan, Iran.            Tel-1:  +98 311 391 3733<br>
Fax: +98 311 391 2376      Tel-2:  +98 311 391 2375</td></tr></table><br>