Hello,

I am trying to study the Li adsorption on 2d g-GeC using Quantum Espresso. According to this paper "Article Thermodynamics and kinetics of 2D g-GeC monolayer as an anod...

" when the Li starts from the hexagon centre it will adsorb at the same location with the minimum adsorption energy.

But when I am trying this it moves to adsorb above Ge atom. I tried different: psedopotentials, Li heights, K-points, cutoff.

The input and output pictures are attached.

The QE input file:

&CONTROL

calculation = "vc-relax"

title='1-site',

prefix='1-site',

nstep = 1000,

verbosity='high',

restart_mode='from_scratch'

pseudo_dir = "........"

outdir='......',

tprnfor = .TRUE.

tstress = .TRUE.

forc_conv_thr=0.002

etot_conv_thr=0.002

/

&SYSTEM

celldm(1)= 25.212280,

celldm(3)= 1.499052,

degauss = 3.70000e-04,

ecutrho = 5.00000e+02,

ecutwfc = 5.00000e+01,

ibrav = 0,

nat = 33,

ntyp = 3,

vdw_corr = 'dft-d3',

occupations = "smearing",

smearing = "gaussian",

/

&ELECTRONS

electron_maxstep = 500,

conv_thr = 1.0d-6,

mixing_mode = 'plain',

mixing_beta = 0.7d0,

/

&IONS

ion_dynamics = "bfgs",

upscale=20.0

/

&CELL

press_conv_thr = 0.5D0

cell_dynamics = 'bfgs',

cell_dofree = '2Dxy'

cell_factor = 5D0

/

K_POINTS {automatic}

7 7 1 0 0 0

ATOMIC_SPECIES

Ge 72.61000 Ge.pbe-n-kjpaw_psl.1.0.0.UPF

C 12.01070 C.pbe-n-kjpaw_psl.1.0.0.UPF

Li 6.94100 Li.pbe-sl-kjpaw_psl.1.0.0.UPF

CELL_PARAMETERS (alat= 25.21228009)

0.978856540 -0.000000286 0.000000000

-0.489428518 0.847715729 0.000000000

0.000000000 0.000000000 1.499052149

ATOMIC_POSITIONS (angstrom)

Ge -0.0685088180 1.7663353790 10.0000000000

C 1.5639746730 0.8238744820 10.0000000000

Ge -1.7009846540 4.5938866400 10.0000000000

C -0.0685020540 3.6513305800 10.0000000000

Ge 3.1964211260 1.7663085260 10.0000000000

C 4.8288204320 0.8237634790 10.0000000000

Ge 1.5639541240 4.5938440260 10.0000000000

C 3.1963711390 3.6512803200 10.0000000000

Ge -3.3333384800 7.4214564210 10.0000000000

C -1.7009145670 6.4789401420 10.0000000000

Ge -4.9657855320 10.2488869080 10.0000000000

C -3.3333244220 9.3063797920 10.0000000000

Ge -0.0684177790 7.4213840550 10.0000000000

C 1.5640458680 6.4788762370 10.0000000000

Ge -1.7009381880 10.2488567360 10.0000000000

C -0.0684876070 9.3063873220 10.0000000000

Ge 6.4612804110 1.7663669610 10.0000000000

C 8.0938562570 0.8239638020 10.0000000000

Ge 4.8287702250 4.5939908570 10.0000000000

C 6.4612927710 3.6514479780 10.0000000000

Ge 9.7263013060 1.7663954530 10.0000000000

C 11.3587095330 0.8239202840 10.0000000000

Ge 8.0938150710 4.5940050230 10.0000000000

C 9.7262456720 3.6513471870 10.0000000000

Ge 3.1964821670 7.4214834880 10.0000000000

C 4.8289711400 6.4790431520 10.0000000000

Ge 1.5639788520 10.2489855620 10.0000000000

C 3.1964935670 9.3064864880 10.0000000000

Ge 6.4614298680 7.4214916730 10.0000000000

C 8.0938940990 6.4790567610 10.0000000000

Ge 4.8289447090 10.2489942550 10.0000000000

C 6.4613771570 9.3064460800 10.0000000000

Li 3.196371139 5.5363819000 12.500000000

Can anyone help me in this?

More Sarah Gamal's questions See All
Similar questions and discussions