No abstract

The recursion method of Haydock, Heine and Kelly is a powerful tool for calculating diagonal matrix elements of the resolvent of quantum-mechanical Hamiltonian operators by elegantly expressing them in terms of continued fractions. In this paper we extend the recursion method to off-diagonal matrix elements of general (possibly non-Hermitian) operators and apply it to the simulation of molecular optical absorption and photoemission spectra within time-dependent density-functional and many-body perturbation theories, respectively. This method is demonstrated with a couple of applications to the optical absorption and photoemission spectra of the caffeine molecule.

We investigate lattice thermal conductivity κ of MgSiO3 perovskite (pv) by ab initio lattice dynamics calculations combined with exact solution of linearized phonon Boltzmann equation. At room temperature, κ of pristine MgSiO3 pv is found to be 10.7 W/(m · K) at 0 GPa. It increases linearly with pressure and reaches 59.2 W/(m · K) at 100 GPa. These values are close to multi-anvil press measurements whereas about twice as large as those from diamond anvil cell experiments. The increase of k with pressure is attributed to the squeeze of weighted phase-spaces phonons get emitted or absorbed. Moreover, we find κ exhibits noticeable anisotropy, with κ zz being the largest component and being about 25%. Such extent of anisotropy is comparable to those of upper mantle minerals such as olivine and enstatite. By analyzing phonon mean free paths and lifetimes, we further show that the weak temperature dependence of κ observed in experiments should not be caused by phonons reaching ‘minimum’ mean free paths. These results clarify the microscopic mechanism of thermal transport in MgSiO3 pv, and provide reference data for understanding heat conduction in the Earth’s deep interior.

The two-phase thermodynamic (2PT) model [S.-T. Lin et al., J. Chem. Phys. 119, 11792-11805 (2003)] provides a promising paradigm to efficiently determine the ionic entropies of liquids from molecular dynamics. In this model, the vibrational density of states (VDoS) of a liquid is decomposed into a diffusive gas-like component and a vibrational solid-like component. By treating the diffusive component as hard sphere (HS) gas and the vibrational component as harmonic oscillators, the ionic entropy of the liquid is determined. Here we examine three issues crucial for practical implementations of the 2PT model: (i) the mismatch between the VDoS of the liquid system and that of the HS gas; (ii) the excess entropy of the HS gas; (iii) the partition of the gas-like and solid-like components. Some of these issues have not been addressed before, yet they profoundly change the entropy predicted from the model. Based on these findings, a revised 2PT formalism is proposed and successfully tested in systems with Lennard-Jones potentials as well as many-atom potentials of liquid metals. Aside from being capable of performing quick entropy estimations for a wide range of systems, the formalism also supports fine-tuning to accurately determine entropies at specific thermal states.