Periodic local MP2 method employing orbital specific virtuals

Abstract

Articles you may be interested inA new ab initio potential energy surface for the collisional excitation of HCN by para-and ortho-H2 J. Chem. Phys. 139, 224301 (2013) We introduce orbital specific virtuals (OSVs) to represent the truncated pair-specific virtual space in periodic local Møller-Plesset perturbation theory of second order (LMP2). The OSVs are constructed by diagonalization of the LMP2 amplitude matrices which correspond to diagonal Wannier-function (WF) pairs. Only a subset of these OSVs is adopted for the subsequent OSV-LMP2 calculation, namely, those with largest contribution to the diagonal pair correlation energy and with the accumulated value of these contributions reaching a certain accuracy. The virtual space for a general (non diagonal) pair is spanned by the union of the two OSV sets related to the individual WFs of the pair. In the periodic LMP2 method, the diagonal LMP2 amplitude matrices needed for the construction of the OSVs are calculated in the basis of projected atomic orbitals (PAOs), employing very large PAO domains. It turns out that the OSVs are excellent to describe short range correlation, yet less appropriate for long range van der Waals correlation. In order to compensate for this bias towards short range correlation, we augment the virtual space spanned by the OSVs by the most diffuse PAOs of the corresponding minimal PAO domain. The Fock and overlap matrices in OSV basis are constructed in the reciprocal space. The 4-index electron repulsion integrals are calculated by local density fitting and, for distant pairs, via multipole approximation. New procedures for determining the fit-domains and the distant-pair lists, leading to higher efficiency in the 4-index integral evaluation, have been implemented. Generally, and in contrast to our previous PAO based periodic LMP2 method, the OSV-LMP2 method does not require anymore great care in the specification of the individual domains (to get a balanced description when calculating energy differences) and is in that sense a black box procedure. Discontinuities in potential energy surfaces, which may occur for PAO-based calculations if one is not careful, virtually disappear for OSV-LMP2. Moreover, due to much increased compactness of the pair-specific virtual spaces, the OSV-LMP2 calculations are faster and require much less memory than PAO-LMP2 calculations, despite the noticeable overhead of the initial OSV construction procedure. C 2015 AIP Publishing LLC.[http://dx

Understanding and Calibrating Density-Functional-Theory Calculations Describing the Energy and Spectroscopy of Defect Sites in Hexagonal Boron Nitride

“…The recent development [35][36][37][38][39][40] of efficiently scaling techniques of this type gives potential for expanded applications. However, these methods are still not as accurate as is desired, with variations of on average 0.3 eV being found between different approaches.…”

“…However, CCSD and MRCI based methods can only be applied to model compounds, 33 not materials such as 3-D h-BN or 2-D h-BN nanoflakes. A range of improved high-level methods is available [34][35][36][37][38][39][40][41][42][43][44][45][46][47][48][49][50][51][52][53] but unfortunately rarely deployed.…”

Applying the Coupled-Cluster Ansatz to Solids and Surfaces in the Thermodynamic Limit

“…Significant progress has been achieved for manyelectron wave-function-based theories by exploiting the locality of electronic correlation in large molecules and solids. The development of so-called local correlation methods and embedding theories has improved their computational efficiency considerably [15][16][17][18][19][20][21][22][23][24]. However, theories that approximate long-range correlation effects such as van der Waals interactions must carefully be checked for convergence with respect to the employed cutoff parameters to allow for accurate and predictive ab initio studies of real materials.…”

Range-separated double-hybrid density-functional theory applied to periodic systems

“…In particular, the most relevant one is the selection of excitation domains assigned to each occupied orbital, since in the present calculations the recently implemented orbital-specific-virtual method 110 was not adopted. Domains for rare-gas and ionic solids were chosen in order to consider PAOs belonging to two coordinated shells of atoms around the occupied orbital, while for molecular crystals the domains coincide with the molecule to which the orbital belongs.…”