[Pw_forum] graphen with adatoms input
yelena
yelena at ipb.ac.rs
Thu May 30 16:39:58 CEST 2013
I tried to make input for graphene with adatoms in two ways, but
everytime when i open it in crysden to check if it looks ok and I got
this image (in attachment).
Here are inputs:
&CONTROL
prefix='gradat',
calculation='scf',
restart_mode='from_scratch',
outdir = '/home/ipb167/Desktop/espresso/gradat' ,
pseudo_dir = '/home/ipb167/Desktop/espresso/pseudo' ,
/
&SYSTEM
ibrav= 0, celldm(1) = 4.66, nat=3, ntyp= 2,
ecutwfc = 65.D0,
occupations='smearing',
smearing='methfessel-paxton',
degauss=0.01,
/
&ELECTRONS
conv_thr=1.D-12,
diagonalization='david',
mixing_mode='plain'
/
ATOMIC_SPECIES {alat}
C 12.0107 C.pz-vbc.UPF
Li 6.941 Li.pz-n-vbc.UPF
CELL_PARAMETERS (alat=4.66000000){hexagonal}
0.999984012 0.000000000 0.000000000
-0.499992006 0.866011558 0.000000000
0.000000000 0.000000000 12.188495864
ATOMIC_POSITIONS crystal
C 0.000000 0.000000 0.000000
C 0.333333 -0.33333 0.000000
Li 0.666666 0.333333 0.000000
K_POINTS automatic
32 32 1 0 0 0
And other:
&CONTROL
prefix='gradat',
calculation='scf',
restart_mode='from_scratch',
outdir = '/home/ipb167/Desktop/espresso/gradat' ,
pseudo_dir = '/home/ipb167/Desktop/espresso/pseudo' ,
/
&SYSTEM
ibrav= 4, celldm(1) = 4.66, celldm(3)= 6.0, nat=3, ntyp= 2,
ecutwfc = 65.D0,
occupations='smearing',
smearing='methfessel-paxton',
degauss=0.01,
/
&ELECTRONS
conv_thr=1.D-12,
diagonalization='david',
mixing_mode='plain'
/
ATOMIC_SPECIES
C 12.0107 C.pz-vbc.UPF
Li 6.941 Li.pz-n-vbc.UPF
ATOMIC_POSITIONS crystal
C 0.000000 0.000000 0.000000
C 0.333333 -0.33333 0.000000
Li 0.666666 0.333333 0.000000
K_POINTS automatic
32 32 1 0 0 0
Both give same image and I'm not sure if Li and C should be connected
it that way.
Any suggestions?
Jelena
On 08 May 2013 15:03, Stefano de Gironcoli wrote:
> you can use the cell you have created for the pure graphene, or a
> 2x2,
> 3x3 version of it and then simply add the adatom in the desired place
> .
> the code will recognize the reduced symmetry and will adapt the list
> of
> kpoints accordingly.
>
> stefano
>
> On 05/08/2013 02:23 PM, yelena wrote:
>> Hello everyone!
>> I am asking for suggestion. I'm trying to make input for scf
>> calculation (and then for electron-phonon... ) for graphene with
>> adatoms.
>> I've been doing calculation with monolayer graphene and I used
>> hexagonal lattice to describe it and large c/a with vacuum to
>> simulate
>> monolayer. I've been using GPU version of code and it works
>> perfectly
>> (I'll write something more about that on GPU message board)
>> Now I'd like to add adatoms in the middle of lattice and I'm not
>> sure
>> how to do that. Is there any way to make input for hexagonal lattice
>> with adatom or I have to make free lattice ibrav=0 and then use CELL
>> PARAMETER card?
>>
>> Looking forward to your answers,
>>
>> Jelena Pesic
>> PhD Student,
>> Institute of Physics Belgrade, Serbia
>> _______________________________________________
>> Pw_forum mailing list
>> Pw_forum at pwscf.org
>> http://pwscf.org/mailman/listinfo/pw_forum
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://pwscf.org/mailman/listinfo/pw_forum
-------------- next part --------------
A non-text attachment was scrubbed...
Name: adatoms.png
Type: image/png
Size: 106675 bytes
Desc: not available
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20130530/5a648dae/attachment.png>
More information about the users
mailing list