<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&nbsp; <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 &gt; Ti.scf.in &lt;&lt; EOF<br>&nbsp;&amp;control<br>&nbsp;&nbsp;&nbsp; calculation = 'scf'<br>&nbsp;&nbsp;&nbsp; restart_mode='from_scratch',<br>&nbsp;&nbsp;&nbsp; prefix='Ti',<br>&nbsp;&nbsp;&nbsp; pseudo_dir = '$PSEUDO_DIR',<br>&nbsp;&nbsp;&nbsp; outdir='$TMP_DIR'<br>&nbsp;/<br>&nbsp;&amp;system<br>&nbsp;&nbsp;&nbsp; ibrav=&nbsp; 0, celldm(1)=20,
 nat=&nbsp; 1, ntyp= 1,<br>&nbsp;&nbsp;&nbsp; ecutwfc = 45 , ecutrho=400 , <br>&nbsp;&nbsp;&nbsp; lda_plus_u = .true.<br>&nbsp;&nbsp;&nbsp; hubbard_u(1) = 1.D-20, <br>&nbsp;&nbsp;&nbsp; starting_magnetization(1)= 1.0,<br>&nbsp;&nbsp;&nbsp; occupations='smearing', smearing='gauss', degauss=0.01,<br>&nbsp;&nbsp;&nbsp; nspin=2,<br>&nbsp;/<br>&nbsp;&amp;electrons<br>&nbsp;&nbsp;&nbsp; mixing_beta = 0.7<br>&nbsp;&nbsp;&nbsp; conv_thr =&nbsp; 1.0d-6<br>&nbsp;/<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&nbsp;&nbsp; 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 &lt; Ti.scf.in &gt; Ti.scf.out <br><br><br>for a in 0.0 -0.1 0.1&nbsp; -0.2 0.2&nbsp;&nbsp; ; do<br>rm -rf $TMP_DIR/*<br><br># perturbed calculations (finite alpha)<br>rm -f Ti.scf.in<br>cat &gt; Ti.scf.in &lt;&lt;
 EOF<br>&nbsp;&amp;control<br>&nbsp;&nbsp;&nbsp; calculation = 'scf'<br>&nbsp;&nbsp;&nbsp; restart_mode='from_scratch',<br>&nbsp;&nbsp;&nbsp; prefix='Ti',<br>&nbsp;&nbsp;&nbsp; pseudo_dir = '$PSEUDO_DIR',<br>&nbsp;&nbsp;&nbsp; outdir='$TMP_DIR'<br>&nbsp;/<br>&nbsp;&amp;system<br>&nbsp;&nbsp;&nbsp; ibrav=&nbsp; 0, celldm(1)=20, nat=&nbsp; 1,&nbsp; ntyp= 1,<br>&nbsp;&nbsp;&nbsp; ecutwfc = 45.0 , ecutrho=400 , <br>&nbsp;&nbsp;&nbsp; lda_plus_u = .true.<br>&nbsp;&nbsp;&nbsp; hubbard_u(1) = 1.D-20, <br>&nbsp;&nbsp;&nbsp; starting_magnetization(1)= 1.0,<br>&nbsp;&nbsp;&nbsp; occupations='smearing', smearing='gauss', degauss=0.01,<br>&nbsp;&nbsp;&nbsp; nspin=2,<br>&nbsp;&nbsp;&nbsp; hubbard_alpha(1)=$a<br>&nbsp;/<br>&nbsp;&amp;electrons<br>&nbsp;&nbsp;&nbsp; diagonalization = 'cg'<br>&nbsp;&nbsp;&nbsp; mixing_beta = 0.7<br>&nbsp;&nbsp;&nbsp; conv_thr =&nbsp; 1.0d-6<br>&nbsp;&nbsp;&nbsp; startingpot = 'file'<br>&nbsp;&nbsp;&nbsp; startingwfc =
 'file'<br>&nbsp;&nbsp;&nbsp; diago_thr_init = 1.d-12<br>&nbsp;/<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&nbsp;&nbsp; 47&nbsp; 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 &lt; Ti.scf.in &gt; Ti.pert_$a.out<br><br>done<br><br><br><br>thanks&nbsp; 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>