[Pw_forum] Problem with ibrav=0 in supercell calculations
Mark Mazar
mazar83 at gmail.com
Fri Jun 4 22:49:11 CEST 2010
Your CELL PARAMETERS flag is missing its underscore.
On Fri, Jun 4, 2010 at 2:57 PM, DHRUV SINGH <dhruv.singh at gmail.com> wrote:
> Hi everyone,
> I am trying to do a supercell calculation of InAs (zinc blende structure)
> using the following input file(ibrav=0 and specifying cell parameters with
> cubic symmetry)
> &control
> calculation='scf'
> restart_mode='from_scratch',
> tstress = .true.
> tprnfor = .true.
> prefix='inas',
> pseudo_dir = '/home/ba01/u110/singh36/espresso-4.1.3/pseudo/',
> outdir='./tmp'
> /
> &system
> ibrav=0,
> celldm(1)=11.05,
> nat= 8,
> ntyp= 2,
> ecutwfc =25.0
> ecutrho=250.0
> /
> &electrons
> conv_thr = 1.0d-8
> mixing_beta = 0.7
> /
> ATOMIC_SPECIES
> In 114.818 In.pz-bhs.UPF
> As 74.921 As.pz-bhs.UPF
> CELL PARAMETERS{cubic}
> 1.0 0.0 0.0
> 0.0 1.0 0.0
> 0.0 0.0 1.0
> ATOMIC_POSITIONS {crystal}
> In 0.00 0.00 0.00
> As 0.25 0.25 0.25
> In 0.50 0.50 0.00
> As 0.75 0.75 0.25
> In 0.00 0.50 0.50
> As 0.25 0.75 0.75
> In 0.50 0.00 0.50
> As 0.75 0.25 0.75
> K_POINTS {automatic}
> 4 4 4 0 0 0
> Strangely enough an scf run does not read the cell parameters card and I get
> an error
> Waiting for input...
> Warning: card CELL PARAMETERS{CUBIC} ignored
> Warning: card 1.0 0.0 0.0 ignored
> Warning: card 0.0 1.0 0.0 ignored
> Warning: card 0.0 0.0 1.0 ignored
> %%%%%%%%%%%%%%%%%%%%%%%%%
> from read_cards : error # 1
> ibrav=0: must read cell parameters
> %%%%%%%%%%%%%%%%%%%%%%%%%
> stopping ...
> --------
> I have run a few calculations in the past with the cubic supercell
> specifying different values of cell parameters and they work fine. Also the
> calculation runs fine with ibrav=1 (and not specifying the cell parameters).
> However, I am not able to pin down the mistake I am making here (Are we not
> allowed to specify a free lattice if it has equal cell parameters?).
> Thanks,
> Dhruv Singh
> School of Mechanical Engineering, Purdue University
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
>
More information about the users
mailing list