Stationary interatomic force constant matrix method with ultrasoft pseudopotentials

2020 ◽  
Vol 102 (17) ◽  
Author(s):  
Kazutoshi Miwa ◽  
Ryoji Asahi
2003 ◽  
Vol 81 (1-2) ◽  
pp. 115-122 ◽  
Author(s):  
A S Cote ◽  
I Morrison ◽  
X Cui ◽  
S Jenkins ◽  
D K Ross

We present the results of first-principles computational studies of the dynamical properties of hexagonal ice using both the ab-initio pseudopotential method and the full-potential augmented plane-wave method. Properties obtained using both the generalized gradient approximation (GGA) and the meta-GGA in density-functional theory are compared. The lattice-dynamical properties of the structures are obtained using a finite-difference evaluation of the dynamical matrix and force-constant matrix from atomic forces. Phonon dispersion is evaluated by the direct determination of the force-constant matrix in supercells derived from the primitive molecule unit cells with the assumption that force constants are zero beyond the second molecular nearest neighbors. The k-dependent phonon frequencies are then obtained from the force-constant matrix and dispersion relations, and the Brillouin-zone integrated density of states is evaluated. The importance of phonon dispersion in the various regions of the phonon spectra is then assessed and compared to existing neutron-scattering data. Frozen-phonon calculations are used to compare phonon frequencies evaluated in both the GGA and meta-GGA. PACS Nos.: 61.12Ex, 63.20Dj


It is shown that the force constants of a solid cannot be deduced from the lattice frequencies alone. The fact that a set of force constants agrees with the frequencies exactly is no guarantee that the force constants are even approximately correct. Detailed eigenvector measurements would resolve the uncertainty. A solid has many more force constants than frequencies. The usual procedure gets round this by retaining only the largest force constants and then using at least as many frequencies as the number of retained force constants to calculate the latter. This procedure is incorrect, because neglecting the small force constants limits the accuracy to which the frequency data can be used to calculate the retained constants. In fact, the number of data that can be extracted from the frequency measurements is always less than the number of force constants one wants to calculate. Therefore, even if all the lattice frequencies were known exactly they could still be satisfied with a very wide range of very different sets of force constants. A large proportion of these sets cannot be rejected on the basis of physical criteria alone. Using the methods of continuous transformation theory, for diamond all the ways are constructed in which the force constant matrix can be changed continuously without altering the agreement with the frequencies. A numerical example is included.


1991 ◽  
Vol 46 (5) ◽  
pp. 426-432
Author(s):  
Zdenek Slanina

AbstractVarious refined potentials describing the intra- and inter-molecular force fields of water molecules arc used to calculate the properties of the gas-phase water dimer. The intra-molecular parts have been taken from spectroscopic or quantum-chemical sources. The minimum energy structure was found iteratively using the first derivatives of the potential; the force-constant matrix was constructed by numerical difierentation. A quite close agreement between the Bopp-Jancso-Heinzinger and the Matsuoka-Clementi-Yoshimine-Lie potentials is found. The treatment is applied to seven observed water-dimer isotopomeric isomerizations


Sign in / Sign up

Export Citation Format

Share Document