The random phase approximation (RPA) systematically overestimates the magnitude of the correlation energy and generally underestimates cohesive energies. This originates in part from the complete lack of exchange terms that would otherwise cancel Pauli exclusion principle violating (EPV) contributions. The uncanceled EPV contributions also manifest themselves in form of an unphysical negative pair density of spin-parallel electrons close to electron–electron coalescence. We follow considerations of many-body perturbation theory to propose an exchange correction that corrects the largest set of EPV contributions, while having the lowest possible computational complexity. The proposed method exchanges adjacent particle/hole pairs in the RPA diagrams, considerably improving the pair density of spin-parallel electrons close to coalescence in the uniform electron gas (UEG). The accuracy of the correlation energy is comparable to other variants of second-order screened exchange (SOSEX) corrections although it is slightly more accurate for the spin-polarized UEG. Its computational complexity scales as O(N5) or O(N4) in orbital space or real space, respectively. Its memory requirement scales as O(N2).
No abstract is provided for this article.
No abstract is provided for this article.
The structure of Brønsted acid sites and the effects of the Al → Si and Si → P substitution in chabazite and AlPO4-34 frameworks, respectively, have been studied by ab initio plane-wave calculations on periodic models. Geometry distortions remain well localized to the nearest neighborhood of the Brønsted acid site. The Si−O(H) and Al−O(H) distances of the common Si−O(H)−Al moiety are significantly different in the two kinds of materials. Calculated OH stretching frequencies of acidic protons, which could be assigned to the low- and high-frequency OH bands, in agreement with experimental infrared spectra of H-SSZ-13 (high silica chabazite) and HSAPO-34, correlate with OH bond lengths and Si−O(H)−Al angles.
I will discuss recent developments in ab initio molecular dynamics (MD) based on the density functional theory (DFT). Presently, most such calculations are performed using the pseudopotential approximation. With the introduction of the projector augmented-wave method, the exact valence wave functions can be applied during MD without worsening the efficiency of the calculations significantly. Hence, calculations are now exact at least in the framework established by DFT. This allows attention to be turned to other limitations, such as system size or the present semilocal approximations used in DFT. I will show results of a comprehensive study for the liquid metallic (and semiconducting) p elements Al, Ga, In, Tl, Si, Ge, Sn, Pb, As, Sb, Bi, Se and Te. Compared to experiment, discrepancies in the structure factor are observed in particular for heavy elements and for As, Se and Te. Possible reasons for these discrepancies are discussed.
No abstract is provided for this article.
The high efficiency of lead organo-metal-halide perovskite solar cells has raised many questions about the role of the methylammonium (MA) molecules in the Pb-I framework. Experiments indicate that the MA molecules are able to ``freely'' spin around at room temperature even though they carry an intrinsic dipole moment. We have performed large supercell (2592 atoms) finite-temperature ab initio molecular dynamics calculations to study the correlation between the molecules in the framework. An underlying long-range antiferroelectric ordering of the molecular dipoles is observed. The dynamical correlation between neighboring molecules shows a maximum around room temperature in the mid-temperature phase. In this phase, the rotations are slow enough to (partially) couple to neighbors via the Pb-I cage. This results in a collective motion of neighboring molecules in which the cage acts as the mediator. At lower and higher temperatures, the motions are less correlated.
Determining the finite-temperature structure of the hybrid perovskite MAPbI3 is a challenge for both experimental and theoretical methods. A very powerful computational method that can resolve the atomic structure is molecular dynamics (MD). The resulting structure depends on the density functional approximation (DFA) in the case of ab initio MD and the force field in classical MD. We compare the structure between 250 and 400 K obtained with different DFAs and force fields in one consistent manner. The symmetry of the PbI3 framework is analyzed as well as the relative ordering of the neighboring organic molecules inside the framework. The distribution function of the molecules is used to map out an effective energy surface for the rotation of a single molecule. This surface is accurately modeled by a pair of cubic harmonics. Available experimental data in literature are discussed and compared to the structure obtained with the different methods. The spread in these data is still too large to uniquely determine the method that “best” describes the perovskite, however, promising candidates and outliers have been identified.