3

ABCAP : advanced — qsub csh —Ferromagnetic Fe

Prepare a directory for the calculation:

August 30, 2017----------------------------cd ~/abc
mkdir Fe
cd Fe../Setnew.csh-----------------------------

Purpose:(1) Study the symmetry of crystal (symmetry operations).
(2) Study the symmetry of wavefunction (irreduciblerepresentations).

Enter a command, Hqc, which shows a procedure of thecalculation.
The procedure consists of(3) Study the optical transition between valence andconduction bands (dipole allowed).

• Edit a file, ???.data.
(4) Calculate the optical absorption spectrum.

• Execute a shell-script file, ???.?sh.

(5) Treat magnetic materials.

[1] Edit the file, ab_prp.data by using the following
information:(6) Use the +U method.

1>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
Fea=2.87A (229, Oh9, Im-3m)
body-centered (il=3) cubic latticegenerators
5 (0, 0, 0)
19 (0, 0, 0)
25 (0, 0, 0)
atomic position (0.0, 0.0, 0.0)
magnetic (jmag=2)>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

Environment variables

Let’s use tcsh.
In order to set some environment variables, executethe following commands:---------------------------------------------cat /home/CMD/teac02/.envs >> ~/.tcshrc
source ~/.tcshrc--------------------------------------------------------------------------(ab_prp.data)
ab_prp.qc-----------------------------

The content of this .envs file is######### .envs ###########
setenv ABCAP /home/CMD/teac02/abcap1707

Execute ab_prp.qc: This is used for executing
ab_prp.csh.
It makes a file, ab_input.data, which is a real input file for the ABCAP calculation.

setenv FC_TYPE gen
setenv FC ifort
setenv FO "-save"setenv MPIFC "mpif90"

The file, ab_input.data, is made from the informations of ab_prp.data and atom.data. The
atom.data file is a data base for atoms from H toU.

set path=(. $path)######### .envs ###########[2] Calculate the initial charge density:

2Preparation----------------------------(ab_input.data) ab_in.qc-----------------------------

Make a directory, e.g., abc, for the ABCAP calculation.
Copy a shell-script file, Setnew.csh, as follows.

In ab_in.qc, the following programs run.

----------------------------cd ~mkdir abc
cp $ABCAP/samples/Setnew.csh abc/-----------------------------• ab_in.exe checks the crystal structure, andmakes totally-symmetric basis functions.
• ab_inch.exe makes the initial charge density
as the superposition of atomic charge densities.

The shell-script file, Setnew.csh, is used on starting a
new calculation: Files needed for ABCAP calculation
are copied from the directory,$ABCAP/samples/LaMnO3c_f_6/.

• ab_kpgn.exe generates the k points used in
the calculation.
• ab_size.exe estimates the size of calculation.
1• ab_ospw.exe calculates the overlap integral
among the totally symmetric basis functions.

[7] Drawing of the density of states (DOS)• Edit bn_pdos.data. Give the energy region
to be drawn.

• fl_dm00.exe sets the initial value of the local
density matrix. (diagonal elements are 1/2)

• Run bn_pdos.sh.

[3] Iteration for the selfconsistent calculation. Edit

• Edit p2_dos.data. Select the total DOS, thes, p, d, and f partial DOS.

fl06.csh and set the number of iteration at theparameter, ITER_MAIN.

• Run p2_dos.sh----------------------------(ab_input.data) fl06.qc-----------------------------•

In fl06.qc, the following programs run:•fl_pot.exe makes the potential from thecharge densty.

4• fl_potuj.exe makes the +U potential fromthe local density matrix.

The k group[1] Irreducible representations for the k group• fl_bnd.exe calculates eigen functions
through a standard diagonalization.

• Edit ab_irep.data: Give k points to be calculated.

• fl_chg.exe determines the Fermi energy andcalculate the charge desity from the eigen functions.

>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
abcap-ab_irep.data1! the number of k
0 0 0 1
! k point>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>• fl_dmmx.exe calculates the local density matrix.
• fl_pot.exe calculates the total energy.
•plot.ps is obtained for the nonmagnetic
case (jmag=0) or the antiferromagnetic case(jmag=1); plot1.ps and plot2.ps are obtained for the magnetic case (jmag=2).

fl_mx5.exe calculates the next input ofcharge density.

• Run ab_irep.sh.
• See ab_irep.log.

Here, fl_ptuj.exe works only with lda+u> 0.

[2] Compatibility in k space

A summary on each cycle of the iterative calculation is written in the file, iter.log. The summaryis shown on the display by the commnad,•----------------------------check.sh-----------------------------

Edit ab_compat.data: Give a series of kpoints for which you want to know the compatibility between the successive k points.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
abcap-ab_compat.data2! the number of k points
0 0 0 1
! k point
1 0 0 8! k point>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>[4] Drawing of the crystal structure• Edit bn_atps.data. Give the region of drawing in the file.
• Run bn_atps.sh and p3_atps.sh.

• Run ab_compat.sh.

• A postscript file, plot.ps, is obtained.

• See ab_compat.log.

[5] Drawing of the band structure (e-k curve)[3] Optical transition (dipole allowed)• Edit a_bnpl.data. Give paths in the k space
to be drawn.

• Edit optsel.data: Give k points to be calculated.

• Run bnpl.sh.
•plot.ps is obtained for the nonmagnetic
case (jmag=0) or the antiferromagnetic case(jmag=1); plot1.ps and plot2.ps are obtained for the magnetic case (jmag=2).

>>>>>>>>>> optsel.data >>>>>>>>>>>>>
0! jpr (std output control)
0! idoub (0: noSOi or 1: SOi)
0 0 0 1 ! k-point
6! irrep. of vector (r or p)>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>[6] Drawing of the k-space path
• Run brzone.sh.

• Run optsel.sh.

• plot.ps is obtained, which shows he k-spacepath in the e-k curve.

• See optsel.txt.
25Si>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
NiO [a(cubic)=4.17A]
a=2.948635A, c=4.17A
simple (il=1) tetragonal latticegenerators
2 (0, 0, 0)
21 (0, 0, 0)
25 (0, 0, 0)
atomic positions (0.0, 0.0, 0.0) Ni-up(0.5, 0.5, 0.5) Ni-down(0.5, 0.5, 0.0) O(0.0, 0.0, 0.5) O
antiferromagnetic (jmag=1)
1 (1/2, 1/2, 1/2)>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ab_prp.data >>>>>>>>>>>
Si (227, Oh7, Fd-3m)
a=5.4296 A
face-centered (il=2)
generators
5 (0/1, 0/1, 0/1)
19 (1/4, 1/4, 1/4)
25 (1/4, 1/4, 1/4)
atomic position (0.0, 0.0, 0.0)
nonmagnetic (jmag=0)>>>>>>>>> ab_prp.data >>>>>>>>>>>6Antiferromagnetic Cr• Type II AFM NiO>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>Cr andiferromagnetic (221, Oh1, Pm-3m)
a=2.88A (bcc)
simple (il=1) cubic lattice for AFM
generators
5 (0, 0, 0)
19 (0, 0, 0)
25 (0, 0, 0)
atomic positions (0.0, 0.0, 0.0)
(0.5, 0.5, 0.5)
antiferromagnetic (jmag=1)
operation 1 (1/2, 1/2, 1/2)>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

7>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
NiO [a(cubic)=4.17A]
a=2.948635A, c=14.445304A
Rhombohedral (il=-1) lattie
generators
3 (0, 0, 0)
10 (0, 0, 0)
13 (0, 0, 0)
atomic positions (0.0, 0.0, 0.0) Ni-up(0.0, 0.0, 0.5) Ni-down(2/3, 1/3, 1/12) O
antiferromagnetic (jmag=1)
1 (0, 0, 1/2)>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>LDA+U method
10The eﬀect of +U should be studied: The value of theband gap and the feature of valence band top.

Optical property• Nonmagnetic NiO• Calculate the dielctric function within the dipole
interband transition.

• Ferromagnetic NiO• Calculate the optical constants.
• Calculate the optical property of thin film.

• Type II AFM NiO[1] Preparation

8

Use a command, Setopt.csh, to prepare a directory
for the optical-spectrum calculation. We assume
that the band structure calculation has been doneat the directory, NiOr_a2_u6.

Nonmagnetic andferromagnetic NiO>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
NiO
a=4.17A,NaCl structureface-centered (il=2) cubic lattie
generators
5 (0, 0, 0)
19 (0, 0, 0)
25 (0, 0, 0)
atomic positions (0.0, 0.0, 0.0) Ni(0.5, 0.5, 0.5) O
nonmagnetic(jmag=0)
ferromagnetic (jmag=2)>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

9----------------------------cd ~/abc
cp $ABCAP/samples/Setopt.csh .
Setopt.csh NiOr_a2_u6----------------------------This makes a directory, NiOr_a2_u6_s, and copies
some files from $ABCAP/samples/ZnO_5_s0 and
NiOr_a2_u6 to this directory.

[2] Initial check----------------------------cd ~/abc/NiOr_a2_u6_sHbnpl.sh ; gs plot.ps

Antiferromagnetic NiO• Type I AFM
3

12bngap.sh ; cat bngap.log; cat bngap.txt-----------------------------Al>>>>>>>>> ab_prp.data >>>>>>>>>>>Al (225, Oh5, Fm-3m) fcc
a=4.05 A
face-centered (il=2)
generators
5 (0, 0, 0)
19 (0, 0, 0)
25 (0, 0, 0)atomic position (0.0, 0.0, 0.0)
nonmagnetic (jmag=0)>>>>>>>>> ab_prp.data >>>>>>>>>>>[3] Adjustment of the band gap----------------------------(bn_scis.data) bn_scis.shbngap.sh ; cat bngap.log----------------------------bn_scis.sh makes a scissors operation.

[4] Carrier doping

13----------------------------(bn_dope.data) bn_dope.sh(bnpl.data)
bnpl.sh-----------------------------Space group information• generators:
File: generator.data• Wycoﬀ positions:File: wycoff.data[5] Calculation of the dielectric function----------------------------(df01.data) df01.qc-----------------------------• To make ab_prp.data,
there is a tool, ab_crystal.sh.

[6] Contribution of plasma oscillation

14----------------------------(wpw.data) wpw.qc
bn_efp2.sh(df_total.data) df_total.sh-----------------------------

Crystal structures

NIMS database:
http://crystdb.nims.go.jp/crystdb/
search-materials• Search ’nims’ and ’atomwork’ by using internet.

[7] Drawing of optical constants• Click ”Enter AtomWork”----------------------------(p2_dielec.data) p2_dielec.sh(p2_optnk.data)
p2_optnk.sh(p2_absorb.data) p2_absorb.sh-----------------------------• login.
• Input elements of the compound(example: Al P)
Click ”Search materials”.

[8] Drawing of optical property of the thin film

• Choose a structure.
• Choose a paper from the refernce list.

----------------------------(optra.data)
optra.sh(p2_optra.data) p2_optra.sh-----------------------------

11• See crystal structure (standardized).

AlP>>>>>>>>> ab_prp.data >>>>>>>>>>>AlP (216, Td2, F-43m) zinc blendea=5.45 A
face-centered (il=2)
generators
5 (0, 0, 0)
43 (0, 0, 0)
atomic positions(0.0, 0.0, 0.0 ) Al
(0.25, 0.25, 0.25) P
nonmagnetic (jmag=0)>>>>>>>>> ab_prp.data >>>>>>>>>>>
4