Random-phase approximation and its applications in computational chemistry and materials science

The random-phase approximation (RPA) as an approach for computing the electronic correlation energy is reviewed. After a brief account of its basic concept and historical development, the paper is devoted to the theoretical formulations of RPA, and its applications to realistic systems. With several illustrating applications, we discuss the implications of RPA for computational chemistry and materials science. The computational cost of RPA is also addressed which is critical for its widespread use in future applications. In addition, current correction schemes going beyond RPA and directions of further development will be discussed.Comment: 25 pages, 11 figures, published online in J. Mater. Sci. (2012

The Molpro quantum chemistry package

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework

Accurate Treatment of Large Supramolecular Complexes by Double-Hybrid Density Functionals Coupled with Nonlocal van der Waals Corrections

In this work, we present a thorough assessment of the performance of some representative double-hybrid density functionals (revPBE0-DH-NL and B2PLYP-NL) as well as their parent hybrid and GGA counterparts, in combination with the most modern version of the nonlocal (NL) van der Waals correction to describe very large weakly interacting molecular systems dominated by noncovalent interactions. Prior to the assessment, an accurate and homogeneous set of reference interaction energies was computed for the supramolecular complexes constituting the L7 and S12L data sets by using the novel, precise, and efficient DLPNO-CCSD(T) method at the complete basis set limit (CBS). The correction of the basis set superposition error and the inclusion of the deformation energies (for the S12L set) have been crucial for obtaining precise DLPNO-CCSD(T)/CBS interaction energies. Among the density functionals evaluated, the double-hybrid revPBE0-DH-NL and B2PLYP-NL with the three-body dispersion correction provide remarkably accurate association energies very close to the chemical accuracy. Overall, the NL van der Waals approach combined with proper density functionals can be seen as an accurate and affordable computational tool for the modeling of large weakly bonded supramolecular systems.Financial support by the “Ministerio de Economía y Competitividad” (MINECO) of Spain and European FEDER funds through projects CTQ2011-27253 and CTQ2012-31914 is acknowledged. The support of the Generalitat Valenciana (Prometeo/2012/053) is also acknowledged. J.A. thanks the EU for the FP7-PEOPLE-2012-IEF-329513 grant. J.C. acknowledges the “Ministerio de Educación, Cultura y Deporte” (MECD) of Spain for a predoctoral FPU grant

van der Waals Interactions in Material Modelling

Van der Waals (vdW) interactions stem from electronic zero-point fluctuations and are often critical for the correct description of structure, stability, and response properties of molecules and materials, including biomolecules, nanomaterials, and material interfaces. Here, we give a conceptual as well as mathematical overview of the current state of modeling vdW interactions,focusing in particular on the consequences of different approximations for practical applications. We present a systematic classification of approximate first-principles models based on the adiabatic-connection fluctuation-dissipation theorem, namely the nonlocal density functionals, interatomic methods, and methods based on the random-phase approximation. The applicability of these methods to different types of materials and material properties is discussed in connection with availability of theoretical and experimental benchmarks. We conclude with a roadmap of the open problems that remain to be solved to construct a universal, efficient, and accurate vdW model for realistic material modeling

Diffusion Monte Carlo Study of Para -Diiodobenzene Polymorphism Revisited

We revisit our investigation of the diffusion Monte Carlo (DMC) simulation of p-DIB molecular crystal polymorphism. [J. Phys. Chem. Lett. 2010, 1, 1789-1794] We perform, for the first time, a rigorous study of finite-size effects and choice of nodal surface on the prediction of polymorph stability in molecular crystals using fixed-node DMC. Our calculations are the largest which are currently feasible using the resources of the K computer and provide insights into the formidable challenge of predicting such properties from first principles. In particular, we show that finite-size effects can influence the trial nodal surface of a small (1×1×1) simulation cell considerably. We therefore repeated our DMC simulations with a 1×3×3 simulation cell, which is the largest such calculation to date. We used a DFT nodal surface generated with the PBE functional and we accumulated statistical samples with ∼6.4×105 core-hours for each polymorph. Our final results predict a polymorph stability consistent with experiment, but indicate that results in our previous paper were somewhat fortuitous. We analyze the finite-size errors using model periodic Coulomb (MPC) interactions and kinetic energy corrections, according to the CCMH scheme of Chiesa, Ceperley, Martin, and Holzmann. We investigate the dependence of the finite-size errors on different aspect ratios of the simulation cell (k-mesh convergence) in order to understand how to choose an appropriate ratio for the DMC calculations. Even in the most expensive simulations currently possible, we show that the finite size errors in the DMC total energies are far larger than the energy difference between the two polymorphs, although error cancellation means that the polymorph prediction is accurate. Finally, we found that the T-move scheme is essential for these massive DMC simulations in order to circumvent population explosions and large time-step biases.Chemistry and Chemical Biolog

Comment on "The role of electron-electron interactions in two-dimensional Dirac fermions"

Tang et al. [Science 361, 570 (2018)] report on the properties of Dirac fermions with both on-site and Coulomb interactions. The substantial decrease up to ~40% of the Fermi velocity of Dirac fermions with on-site interaction is inconsistent with the numerical data near the Gross-Neveu quantum critical point. This results from an inappropriate finite-size extrapolation.Comment: Submitted: October 24, 2018, Accepted: November 4, 201