# [Pw_forum] Dielectric constant calculations for InAs

Zhiping Xu xuzp at tsinghua.edu.cn
Tue Sep 8 11:48:30 CEST 2015

```Dear all:

I am using pw.x and ph.x to calculate the dielectric constant of InAs, and here are some problems with the calculations.

(1) After structural relaxation using vc-relax, the stress is now zero (confirmed using scf calculations for the final structure).

convergence has been achieved in  10 iterations

Forces acting on atoms (Ry/au):

atom    1 type  1   force =     0.00000000    0.00000000    0.00000000
atom    2 type  2   force =     0.00000000    0.00000000    0.00000000

Total force =     0.000000     Total SCF correction =     0.000000

entering subroutine stress ...

total   stress  (Ry/bohr**3)                   (kbar)     P=    0.00
0.00000001  -0.00000000  -0.00000000          0.00     -0.00     -0.00
-0.00000000   0.00000001   0.00000000         -0.00      0.00      0.00
-0.00000000   0.00000000   0.00000001         -0.00      0.00      0.00

However, the following phonon calculation predicts not only negative phonon frequency at q = 0, but also very strange dielectric constant.

Dielectric constant in cartesian axis

(    -230.720680077      -0.000000000       0.000000000 )
(      -0.000000000    -230.720680077       0.000000000 )
(       0.000000000       0.000000000    -230.720680077 )

Effective charges (d Force / dE) in cartesian axis

atom      1   In
Ex  (        9.55076        0.00000        0.00000 )
Ey  (        0.00000        9.55076        0.00000 )
Ez  (        0.00000        0.00000        9.55076 )
atom      2   As
Ex  (       -9.43886        0.00000        0.00000 )
Ey  (       -0.00000       -9.43886       -0.00000 )
Ez  (        0.00000        0.00000       -9.43886 )

Diagonalizing the dynamical matrix

q = (    0.000000000   0.000000000   0.000000000 )

**************************************************************************
freq (    1) =      -6.064621 [THz] =    -202.293974 [cm-1]
freq (    2) =      -6.064621 [THz] =    -202.293974 [cm-1]
freq (    3) =      -6.064621 [THz] =    -202.293974 [cm-1]
freq (    4) =       0.164141 [THz] =       5.475159 [cm-1]
freq (    5) =       0.164141 [THz] =       5.475159 [cm-1]
freq (    6) =       0.164141 [THz] =       5.475159 [cm-1]
**************************************************************************

The input files are attached below for your reference.

(2) One more question, the vc-relax calculation converges at zero stress, which was calculated using the cut-off set for initial cell parameters, so one has to redo the calculations using the new cell parameters following the suggestion at http://www.quantum-espresso.org/faq/self-consistency/#6.11 <http://www.quantum-espresso.org/faq/self-consistency/#6.11>. Is there any better (more automatic) way to obtain a final relaxed cell/atom structures?

Could any one help me on these two questions?

Thank you,

Zhiping

input files:

—scf.in—
&control
calculation     = 'scf',
restart_mode    = 'from_scratch',
pseudo_dir      = '/home/xuzp/bin/lib/upf_files',
outdir          = '.',
prefix          = 'InAs',
tstress         = .true.
tprnfor         = .true.
wf_collect      = .true.
/
&system
ibrav           = 2,
celldm(1)       = 11.69675,
nat             = 2,
ntyp            = 2,
ecutwfc         = 90,
/
&electrons
diagonalization = 'cg'
conv_thr        =  1.0d-10
/

ATOMIC_SPECIES
In  111.818  In.pbe-dn-rrkjus_psl.0.2.2.UPF
As   74.920  As.pbe-n-rrkjus_psl.0.2.UPF

ATOMIC_POSITIONS crystal
In 0.00 0.00 0.00
As 0.25 0.25 0.25

K_POINTS automatic
9 9 9 0 0 0

—ph.in—
phonons of InAs
&inputph
tr2_ph   = 1.0d-16,
niter_ph = 100,
prefix   = 'InAs',
epsil    = .true.
amass(1) = 111.818,
amass(2) = 74.920,
outdir   = '.',
fildyn   = 'InAs.matdyn',
/
0.0 0.0 0.0
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20150908/5209f624/attachment.html>
```
{"ft_lang_label":"__label__en","ft_lang_prob":0.6185084,"math_prob":0.99357396,"size":3444,"snap":"2019-26-2019-30","text_gpt3_token_len":1245,"char_repetition_ratio":0.21453488,"word_repetition_ratio":0.057613168,"special_character_ratio":0.49941927,"punctuation_ratio":0.22537112,"nsfw_num_words":0,"has_unicode_error":false,"math_prob_llama3":0.97711337,"pos_list":[0],"im_url_duplicate_count":[null],"WARC_HEADER":"{\"WARC-Type\":\"response\",\"WARC-Date\":\"2019-07-18T08:52:41Z\",\"WARC-Record-ID\":\"<urn:uuid:ea8fc9a8-5428-49ea-b588-bd0f462af530>\",\"Content-Length\":\"7222\",\"Content-Type\":\"application/http; msgtype=response\",\"WARC-Warcinfo-ID\":\"<urn:uuid:2c35fc73-0b70-4573-829e-2c1b20f5a8da>\",\"WARC-Concurrent-To\":\"<urn:uuid:c8ef5118-9369-47ec-b2e6-1ea2bca50b85>\",\"WARC-IP-Address\":\"147.122.3.40\",\"WARC-Target-URI\":\"https://lists.quantum-espresso.org/pipermail/users/2015-September/033220.html\",\"WARC-Payload-Digest\":\"sha1:4H6ZR7GSO4JQC4ZMN6GNLTBYEKYR2ZVD\",\"WARC-Block-Digest\":\"sha1:34OLXXIGNMFOGP3MEGGGYBS2R47FMJT7\",\"WARC-Identified-Payload-Type\":\"text/html\",\"warc_filename\":\"/cc_download/warc_2019/CC-MAIN-2019-30/CC-MAIN-2019-30_segments_1563195525587.2_warc_CC-MAIN-20190718083839-20190718105839-00431.warc.gz\"}"}