These days, I'm having a trouble in atomic relaxing of "EuTiO3" with tetragonal symmetry using Quantum espresso 7.1 for a period of more that a week. This material have a transverse optical mode to z-axis("001" direction). So, I try to relax this material along the z-axis. My order of relaxing the EuTiO3 is as follows:
1. I download "cif.file" including the atomic position & lattice from "Materials Explorer( https://next-gen.materialsproject.org/materials ) "
2. At Quantum espresso, "vc-relax(lattice realxation)" this material.
3. "Relax" this material along the z-axis(= the direction of transverse optical mode,001).
At "2" process, total energy is decreased. But at "3" process, total energy is oddly increased!!. By other paper related EuTiO3, because this material is unstable to the direction of "z-axis", total energy should be also decreased at "3" process. I don't know what the reason is.
At similar structure "SrTiO3" with tetragonal symmetry , this problem didn't occur. At "SrTiO3", the same procedure as above was performed, but total energy was decreased at both "2" & "3" process.
Both "EuTiO3" & "SrTiO3" are perovskite material. One difference is that "EuTiO3" have antiferromagnetic property(=spin polarization) unlike "SrTiO3". Hence, I specified the "magnetization" and "Hubbard potential" to the input data of "EuTiO3"( In another paper, I know that magnetization and Hubbard potential are specified when relaxing EuTiO3 material.)
Did I typed the "magnetization" & "Hubbard potential" data wrong? Please help me. I need sincere help from experts.
P.S I add my input data of "EuTiO3" entered into Quatum Espresso.
2. process : "vc-relax"
&control
calculation = 'vc-relax'
restart_mode='from_scratch',
prefix='ETO',
pseudo_dir ='/home/dydtlr5506/pseudo',
outdir='./wavefunc'
tstress = .true.
forc_conv_thr = 1.0d-7
dt = 20
disk_io = 'none'
nstep =9910
/
&system
ibrav= 0,
nat=20, ntyp= 4,
ecutwfc = 80.0, ecutrho = 640.0,
occupations = 'smearing',
smearing = 'gaussian'
degauss = 0.0010
nspin = 2
starting_magnetization(1)=7.0
starting_magnetization(2)=-7.0
/ &electrons
mixing_mode = 'plain'
electron_maxstep = 100
conv_thr = 1.0d-13
/
&ions
ion_dynamics = 'bfgs'
/ &cell
cell_dofree='xyz'
cell_dynamics = 'bfgs'
/
ATOMIC_SPECIES
Eu1 151.964 Eu.pbe-spdfn-kjpaw_psl.1.0.0.UPF
Eu2 151.964 Eu.pbe-spdfn-kjpaw_psl.1.0.0.UPF
Ti 47.86 Ti.pbe-spn-kjpaw_psl.1.0.0.UPF
O 15.999 O.pbe-n-kjpaw_psl.1.0.0.UPF
K_POINTS {automatic}
6 6 4 0 0 0
CELL_PARAMETERS (angstrom)
5.510044859 0.000000000 0.000000000
0.000000000 5.510044859 0.000000000
0.000000000 0.000000000 7.871680636
ATOMIC_POSITIONS (crystal)
Eu1 0.0000000000 0.5000000000 0.2500000000 1 1 0
Eu2 0.5000000000 0.0000000000 0.2500000000 1 1 0
Eu1 0.5000000000 0.0000000000 0.7500000000 1 1 0
Eu2 0.0000000000 0.5000000000 0.7500000000 1 1 0
Ti 0.5000000000 0.5000000000 0.5000000000 1 1 0
Ti 0.0000000000 0.0000000000 0.5000000000 1 1 0
Ti 0.0000000000 0.0000000000 0.0000000000 1 1 0
Ti 0.5000000000 0.5000000000 0.0000000000 1 1 0
O 0.5000000000 0.5000000000 0.2500000000 1 1 0
O 0.0000000000 0.0000000000 0.2500000000 1 1 0
O 0.7183286073 0.2183286073 0.0000000000 1 1 0
O 0.7816713921 0.2816713922 0.5000000000 1 1 0
O 0.2816713922 0.2183286079 0.5000000000 1 1 0
O 0.2183286073 0.2816713927 0.0000000000 1 1 0
O 0.0000000000 0.0000000000 0.7500000000 1 1 0
O 0.5000000000 0.5000000000 0.7500000000 1 1 0
O 0.2183286079 0.7183286078 0.5000000000 1 1 0
O 0.2816713927 0.7816713927 0.0000000000 1 1 0
O 0.7816713927 0.7183286073 0.0000000000 1 1 0
O 0.7183286078 0.7816713921 0.5000000000 1 1 0
HUBBARD (ortho-atomic)
U Eu1-4f 5.9
U Eu2-4f 5.9
2. process: output total energy = -5712.21091180 Ry
3. process : "relax"
&control
calculation = 'relax'
restart_mode='from_scratch',
prefix='ETO',
pseudo_dir ='/home/dydtlr5506/pseudo',
outdir='./wavefunc'
tstress = .true.
forc_conv_thr = 1.0d-7
dt = 20
disk_io = 'none'
nstep =9910
/
&system
ibrav= 0,
nat=20, ntyp= 4,
ecutwfc = 80.0, ecutrho = 640.0,
occupations = 'smearing',
smearing = 'gaussian'
degauss = 0.0010
nspin = 2
starting_magnetization(1)=0.5
starting_magnetization(2)=-0.5
/ &electrons
mixing_mode = 'plain'
electron_maxstep = 100
conv_thr = 1.0d-9
/
&ions
ion_dynamics = 'bfgs'
/ &cell
cell_dofree='xyz'
cell_dynamics = 'bfgs'
/
ATOMIC_SPECIES
Eu1 151.964 Eu.pbe-spdfn-kjpaw_psl.1.0.0.UPF
Eu2 151.964 Eu.pbe-spdfn-kjpaw_psl.1.0.0.UPF
Ti 47.86 Ti.pbe-spn-kjpaw_psl.1.0.0.UPF
O 15.999 O.pbe-n-kjpaw_psl.1.0.0.UPF
K_POINTS {automatic}
6 6 4 0 0 0
CELL_PARAMETERS (angstrom)
5.505334661 -0.000000000 0.000000000
0.000000000 5.505334661 0.000000000
-0.000000000 -0.000000000 7.873593873
ATOMIC_POSITIONS (crystal)
Eu1 0.0000000000 0.5000000000 0.1739600153 0 0 1
Eu2 0.5000000000 0.0000000000 0.1738399582 0 0 1
Eu1 0.5000000000 0.0000000000 0.8220612019 0 0 1
Eu2 0.0000000000 0.5000000000 0.8225353983 0 0 1
Ti 0.5000000000 0.5000000000 0.4904621300 0 0 1
Ti 0.0000000000 0.0000000000 0.4904621300 0 0 1
Ti 0.0000000000 0.0000000000 -0.0080196857 0 0 1
Ti 0.5000000000 0.5000000000 -0.0080196857 0 0 1
O 0.5000000000 0.5000000000 0.2458119409 0 0 1
O 0.0000000000 0.0000000000 0.2458119409 0 0 1
O 0.7179836696 0.2179836696 -0.0002906743 0 0 1
O 0.7820163277 0.2820163278 0.5028007560 0 0 1
O 0.2820163278 0.2179836723 0.4996782263 0 0 1
O 0.2179836696 0.2820163304 0.0012584417 0 0 1
O 0.0000000000 0.0000000000 0.7710861746 0 0 1
O 0.5000000000 0.5000000000 0.7710861746 0 0 1
O 0.2179836723 0.7179836722 0.4872825726 0 0 1
O 0.2820163304 0.7820163304 0.5132563164 0 0 1
O 0.7820163304 0.7179836696 0.0012584417 0 0 1
O 0.7179836722 0.7820163277 0.4996782263 0 0 1
HUBBARD (ortho-atomic)
U Eu1-4f 5.9
U Eu2-4f 5.9
3. process: output total energy = -5711.56708785 Ry
Thank you.
- Badges
- Science topic
- Similar topics
- Physical Sciences
- Materials Science
- Materials Chemistry
- Polymer Chemistry