[Pw_forum] NEB calculation

Layla Martin-Samos lmartinsamos at gmail.com
Mon Apr 18 13:13:46 CEST 2011


Dear Savita, atention the NEB input has changed in version 4.3. There is an
independent executable that parses a file, see example17. For identifying
the last image configuration, maybe you should make a simple relaxation of
CO at a fcc site and O very close to CO. Then you can start a NEB
calculation with  first image as CO and O at two near fcc sites and the last
image as the relaxed configuration obtained from the previous relaxation.

cheers

Layla

2011/4/18 <psavita at crlindia.com>

> Dear Layla and Prasenjeet,
>
> I think I should take Layla's suggestion and first consider CO2 formation.
> But then can I assume that O moves towards CO and CO2 forms at atop?
> Also as per Prasenjeet's advice, I should reduce vacuum layer to 25 bohr.
> I am giving here my input file for you to comment more.
> Thank you so much.
>
>
> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
> &CONTROL
>   calculation = 'neb',
>   restart_mode = 'from_scratch'
>   prefix='pt_co',
>   pseudo_dir  =
> "/home/nano/sivaram/simulation_zone/dft_simulation/year_2011_work/NEB",
>   outdir      =
> "/home/nano/sivaram/simulation_zone/dft_simulation/year_2011_work/NEB/relax",
> !  nstep       = 50,
> /
> &SYSTEM
>   ibrav       = 0,
>   nat         = 15,
>   ntyp        = 3,
>   ecutwfc     = 60.00,
>   ecutrho     = 360.00,
>   occupations = "smearing",
>   smearing    = "mp",
>   degauss     = 0.014699,
>   nspin       = 2,
>   starting_magnetization = 0.5
> /
> &ELECTRONS
>   electron_maxstep = 100,
>   conv_thr    = 1.D-7,
>   mixing_beta = 0.2,
>   diagonalization = 'cg'
> /
> &IONS
>   ds                = 2.0d0,
>   opt_scheme        = "broyden",
>   num_of_images     = 5,
>   k_max             = 0.4,
>   k_min             = 0.3,
>   CI_scheme         = "no-CI"
>   pot_extrapolation = "second_order",
>   wfc_extrapolation = "second_order",
>   path_thr          = 0.1d0
> /
> CELL_PARAMETERS
> 10.476 0 0
> 0  9.0727 0
> 0 0 40.215
> ATOMIC_SPECIES
> Pt  195.08  Pt.pw91-n-van.UPF
> O   15.999  O.pw91-van_ak.UPF
> C   12.011  C.pw91-van_ak.UPF
> ATOMIC_POSITIONS {bohr}
> first_image
> Pt       0.000000000   0.000000000   0.000000000
> Pt       5.238100000   0.000000000   0.000000000
> Pt       2.619050000   4.536328000   0.000000000
> Pt       7.857150000   4.536328000   0.000000000
> Pt       2.619088571   1.512240622   4.513968243
> Pt       7.857077843   1.512255365   4.514105935
> Pt       5.238032011   6.048542598   4.513881620
> Pt      10.476075613   6.048536230   4.513989805
> Pt       0.000134895   3.029013846   9.028523678
> Pt       5.238095443   3.029034864   9.028457281
> Pt       2.619088295   7.565397431   9.028126889
> Pt       7.857080929   7.565358308   9.028145500
> C       7.857080929   7.565358308    11.5281455
> O       7.857080929   7.565358308    13.7013455
> O       2.619050000   4.536328000    11.5281455
> last_image
> Pt       0.000000000   0.000000000   0.000000000
> Pt       5.238100000   0.000000000   0.000000000
> Pt       2.619050000   4.536328000   0.000000000
> Pt       7.857150000   4.536328000   0.000000000
> Pt       2.619088571   1.512240622   4.513968243
> Pt       7.857077843   1.512255365   4.514105935
> Pt       5.238032011   6.048542598   4.513881620
> Pt      10.476075613   6.048536230   4.513989805
> Pt       0.000134895   3.029013846   9.028523678
> Pt       5.238095443   3.029034864   9.028457281
> Pt       2.619088295   7.565397431   9.028126889
> Pt       7.857080929   7.565358308   9.028145500
> C        3.9286        5.2972       20.000
> O        3.9286        3.1240       20.0
> O        3.9286        7.4704       20.0
> K_POINTS {automatic}
>  6 6 3 0 0 0
>
> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
>
> Savita Pundlik
> Computational Materials Applied Research Group
> Computational Research Laboratories Ltd.,
> Taco House, Damle Path, Off Law College Road
> Pune - 411004, India.
>
>
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110418/bf7adfd6/attachment.htm 


More information about the Pw_forum mailing list