1
|
Iyer GR, Whelpley N, Tiihonen J, Kent PRC, Krogel JT, Rubenstein BM. Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories. J Chem Theory Comput 2024; 20:7416-7429. [PMID: 39172163 DOI: 10.1021/acs.jctc.4c00214] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2024]
Abstract
The accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH3) molecule and the nucleophilic substitution (SN2) reaction F- + CH3F → FCH3 + F-. We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.
Collapse
Affiliation(s)
- Gopal R Iyer
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Noah Whelpley
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Juha Tiihonen
- Department of Physics, Nanoscience Center, University of Jyväskylä, Jyväskylä 40014, Finland
| | - Paul R C Kent
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Brenda M Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
2
|
Shin H, Gasperich K, Rojas T, Ngo AT, Krogel JT, Benali A. Systematic Improvement of Quantum Monte Carlo Calculations in Transition Metal Oxides: sCI-Driven Wavefunction Optimization for Reliable Band Gap Prediction. J Chem Theory Comput 2024. [PMID: 39213586 DOI: 10.1021/acs.jctc.4c00335] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
Accurate determination of the electronic properties of correlated oxides remains a significant challenge for computational theory. Traditional Hubbard-corrected density functional theory (DFT+U) frequently encounters limitations in precisely capturing electron correlation, particularly in predicting band gaps. We introduce a systematic methodology to enhance the accuracy of diffusion Monte Carlo (DMC) simulations for both ground and excited states, focusing on LiCoO2 as a case study. By employing a selected configuration interaction (sCI) approach, we demonstrate the capability to optimize wavefunctions beyond the constraints of single-reference DFT+U trial wavefunctions. We show that the sCI framework enables accurate prediction of band gaps in LiCoO2, closely aligning with experimental values and substantially improving traditional computational methods. The study uncovers a nuanced mixed state of t2g and eg orbitals at the band edges that is not captured by conventional single-reference methods, further elucidating the limitations of PBE+U in describing d-d excitations. Our findings advocate for the adoption of beyond-DFT methodologies, such as sCI, to capture the essential physics of excited-state wavefunctions in strongly correlated materials. The improved accuracy in band gap predictions and the ability to generate more reliable trial wavefunctions for DMC calculations underscore the potential of this approach for broader applications in the study of correlated oxides. This work not only provides a pathway for more accurate simulations of electronic structures in complex materials but also suggests a framework for future investigations of the excited states of other challenging systems.
Collapse
Affiliation(s)
- Hyeondeok Shin
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Kevin Gasperich
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Tomas Rojas
- Department of Chemical Engineering, University of Illinois at Chicago, Chicago, Illinois 60608, United States
- Material Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Anh T Ngo
- Department of Chemical Engineering, University of Illinois at Chicago, Chicago, Illinois 60608, United States
- Material Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Anouar Benali
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| |
Collapse
|
3
|
Szabó PB, Schätzle Z, Entwistle MT, Noé F. An Improved Penalty-Based Excited-State Variational Monte Carlo Approach with Deep-Learning Ansatzes. J Chem Theory Comput 2024. [PMID: 39213603 PMCID: PMC11428158 DOI: 10.1021/acs.jctc.4c00678] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
We introduce several improvements to the penalty-based variational quantum Monte Carlo (VMC) algorithm for computing electronic excited states of Entwistle et al. [Nat. Commun. 14, 274 (2023)] and demonstrate that the accuracy of the updated method is competitive with other available excited-state VMC approaches. A theoretical comparison of the computational aspects of these algorithms is presented, where several benefits of the penalty-based method are identified. Our main contributions include an automatic mechanism for tuning the scale of the penalty terms, an updated form of the overlap penalty with proven convergence properties, and a new term that penalizes the spin of the wave function, enabling the selective computation of states with a given spin. With these improvements, along with the use of the latest self-attention-based ansatz, the penalty-based method achieves a mean absolute error below 1 kcal/mol for the vertical excitation energies of a set of 26 atoms and molecules, without relying on variance matching schemes. Considering excited states along the dissociation of the carbon dimer, the accuracy of the penalty-based method is on par with that of natural-excited-state (NES) VMC, while also providing results for additional sections of the potential energy surface, which were inaccessible with the NES method. Additionally, the accuracy of the penalty-based method is improved for a conical intersection of ethylene, with the predicted angle of the intersection agreeing well with both NES-VMC and multireference configuration interaction.
Collapse
Affiliation(s)
- P Bernát Szabó
- Department of Mathematics and Computer Science, FU Berlin, Arnimallee 6, Berlin 14195, Germany
| | - Zeno Schätzle
- Department of Mathematics and Computer Science, FU Berlin, Arnimallee 6, Berlin 14195, Germany
| | - Michael T Entwistle
- Department of Mathematics and Computer Science, FU Berlin, Arnimallee 6, Berlin 14195, Germany
| | - Frank Noé
- Department of Mathematics and Computer Science, FU Berlin, Arnimallee 6, Berlin 14195, Germany
- Microsoft Research AI4Science, Karl-Liebknecht Str. 32, Berlin 10178, Germany
- Department of Physics, FU Berlin, Arnimallee 14, Berlin 14195, Germany
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
| |
Collapse
|
4
|
Huber TB, Wheeler RA. Fixed-node diffusion Monte Carlo shows promise for modeling reaction thermochemistry of hydrocarbon-based radicals. J Chem Phys 2024; 161:034303. [PMID: 39007382 DOI: 10.1063/5.0211903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2024] [Accepted: 06/12/2024] [Indexed: 07/16/2024] Open
Abstract
Reliable thermodynamic and kinetic properties of free radical polymerization reactions are essential for synthesizing both primary polymeric materials and specialty polymers. The computational generation of these data from quantum chemistry requires a time-efficient method capable of capturing the essential physics. One such method, fixed-node diffusion Monte Carlo (FN-DMC) (using single Slater-Jastrow trial wavefunctions), has demonstrated the capability to recover 90%-95% of missing dynamic correlation energy for typical systems. In this study, methyl radical addition to ethylene serves as a simple model to test FN-DMC's ability to calculate enthalpies of reaction and activation energies with different time steps, antisymmetric trial wavefunctions, basis set sizes, and effective core potentials. The FN-DMC computational protocol thus defined for methyl radical addition to ethylene is subsequently benchmarked against Weizmann-1 and experimental reaction enthalpies from Lin et al.'s test set of 21 radical addition and 28 hydrogen abstraction enthalpies. Our findings reveal that FN-DMC consistently generates reaction enthalpies with chemical accuracy, exhibiting mean absolute deviation of 3.5(7) and 1.4(8) kJ/mol from the Weizmann-1 reference for radical addition and hydrogen abstraction reactions, respectively. Given its favorable computational scaling and high degree of parallelizability, we, therefore, recommend more comprehensive testing of FN-DMC with effective core potentials to address more extensive and intricate polymerization reactions and reactions with other radicals.
Collapse
Affiliation(s)
- Timothy B Huber
- Department of Chemistry and Biochemistry, Northern Illinois University, 1425 W Lincoln Hwy, Dekalb, Illinois 60115, USA
| | - Ralph A Wheeler
- Department of Chemistry and Biochemistry, Northern Illinois University, 1425 W Lincoln Hwy, Dekalb, Illinois 60115, USA
| |
Collapse
|
5
|
Wang M, Zhou Y, Wang H. Performance assessment of the effective core potentials under the fermionic neural network: First and second row elements. J Chem Phys 2024; 160:204109. [PMID: 38785290 DOI: 10.1063/5.0207853] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2024] [Accepted: 05/05/2024] [Indexed: 05/25/2024] Open
Abstract
The rapid development of deep learning techniques has driven the emergence of a neural network-based variational Monte Carlo (VMC) method (referred to as FermiNet), which has manifested high accuracy and strong predictive power in the electronic structure calculations of atoms, molecules, and some periodic systems. Recently, the implementation of the effective core potential (ECP) scheme has further facilitated more efficient calculations in practice. However, there is still a lack of comprehensive assessments of the ECP's performance under the FermiNet. In this work, we set sail to fill this gap by conducting extensive tests on the first two row elements regarding their atomic, spectral, and molecular properties. Our major finding is that, in general, the qualities of ECPs have been correctly reflected under FermiNet. Two recently built ECP tables, namely, correlation consistent ECP (ccECP) and energy consistent correlated electron pseudopotential (eCEPP), seem to prevail in terms of overall performance. In particular, ccECP performs slightly better on spectral precision and covers more elements, while eCEPP is more systematically built from both shape and energy consistency and better treats the core polarization. On the other hand, the high accuracy of the all-electron calculations is hindered by the absence of relativistic effects as well as the numerical instabilities in some heavier elements. Finally, with further in-depth discussions, we generate possible directions for developing and improving FermiNet in the near future.
Collapse
Affiliation(s)
- Mengsa Wang
- Graduate School of China Academy of Engineering Physics, Beijing 100088, China
- National Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
| | - Yuzhi Zhou
- CAEP Software Center for High Performance Numerical Simulation, Beijing 100088, China
- Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
| | - Han Wang
- National Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
- HEDPS, CAPT, College of Engineering, Peking University, Beijing, China
| |
Collapse
|
6
|
Yeh CN, Morales MA. Low-Scaling Algorithms for GW and Constrained Random Phase Approximation Using Symmetry-Adapted Interpolative Separable Density Fitting. J Chem Theory Comput 2024; 20:3184-3198. [PMID: 38597496 DOI: 10.1021/acs.jctc.4c00085] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
We present low-scaling algorithms for GW and constrained random phase approximation based on a symmetry-adapted interpolative separable density fitting (ISDF) procedure that incorporates the space-group symmetries of crystalline systems. The resulting formulations scale cubically, with respect to system size, and linearly with the number of k-points, regardless of the choice of single-particle basis and whether a quasiparticle approximation is employed. We validate these methods through comparisons with published literature and demonstrate their efficiency in treating large-scale systems through the construction of downfolded many-body Hamiltonians for carbon dimer defects embedded in hexagonal boron nitride supercells. Our work highlights the efficiency and general applicability of ISDF in the context of large-scale many-body calculations with k-point sampling beyond density functional theory.
Collapse
Affiliation(s)
- Chia-Nan Yeh
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Miguel A Morales
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|
7
|
Zhou H, Kincaid B, Wang G, Annaberdiyev A, Ganesh P, Mitas L. A new generation of effective core potentials: Selected lanthanides and heavy elements. J Chem Phys 2024; 160:084302. [PMID: 38391016 DOI: 10.1063/5.0180057] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Accepted: 11/05/2023] [Indexed: 02/24/2024] Open
Abstract
We construct correlation-consistent effective core potentials (ccECPs) for a selected set of heavy atoms and f elements that are currently of significant interest in materials and chemical applications, including Y, Zr, Nb, Rh, Ta, Re, Pt, Gd, and Tb. As is customary, ccECPs consist of spin-orbit (SO) averaged relativistic effective potential (AREP) and effective SO terms. For the AREP part, our constructions are carried out within a relativistic coupled-cluster framework while also taking into account objective function one-particle characteristics for improved convergence in optimizations. The transferability is adjusted using binding curves of hydride and oxide molecules. We address the difficulties encountered with f elements, such as the presence of large cores and multiple near-degeneracies of excited levels. For these elements, we construct ccECPs with core-valence partitioning that includes 4f subshell in the valence space. The developed ccECPs achieve an excellent balance between accuracy, size of the valence space, and transferability and are also suitable to be used in plane wave codes with reasonable energy cutoffs.
Collapse
Affiliation(s)
- Haihan Zhou
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Benjamin Kincaid
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Guangming Wang
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Abdulgani Annaberdiyev
- Center for Nanophase Materials Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Panchapakesan Ganesh
- Center for Nanophase Materials Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Lubos Mitas
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| |
Collapse
|
8
|
Shin H, Ganesh P, Kent PRC, Benali A, Bhattacharya A, Lee HN, Heinonen O, Krogel JT. DFT+ U and quantum Monte Carlo study of electronic and optical properties of AgNiO 2 and AgNi 1-xCo xO 2 delafossite. Phys Chem Chem Phys 2024; 26:6967-6976. [PMID: 38334756 DOI: 10.1039/d3cp03477a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
As the only semimetallic d10-based delafossite, AgNiO2 has received a great deal of attention due to both its unique semimetallicity and its antiferromagnetism in the NiO2 layer that is coupled with a lattice distortion. In contrast, other delafossites such as AgCoO2 are insulating. Here we study how the electronic structure of AgNi1-xCoxO2 alloys vary with Ni/Co concentration, in order to investigate the electronic properties and phase stability of the intermetallics. While the electronic and magnetic structure of delafossites have been studied using density functional theory (DFT), earlier studies have not included corrections for strong on-site Coulomb interactions. In order to treat these interactions accurately, in this study we use Quantum Monte Carlo (QMC) simulations to obtain accurate estimates for the electronic and magnetic properties of AgNiO2. By comparison to DFT results we show that these electron correlations are critical to account for. We show that Co doping on the magnetic Ni sites results in a metal-insulator transition near x ∼0.33, and reentrant behavior near x ∼ 0.66.
Collapse
Affiliation(s)
- Hyeondeok Shin
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA.
| | - Panchapakesan Ganesh
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Paul R C Kent
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Anouar Benali
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA.
| | - Anand Bhattacharya
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Ho Nyung Lee
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Olle Heinonen
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA.
| |
Collapse
|
9
|
Ditte M, Barborini M, Medrano Sandonas L, Tkatchenko A. Molecules in Environments: Toward Systematic Quantum Embedding of Electrons and Drude Oscillators. PHYSICAL REVIEW LETTERS 2023; 131:228001. [PMID: 38101380 DOI: 10.1103/physrevlett.131.228001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/29/2023] [Revised: 07/26/2023] [Accepted: 10/20/2023] [Indexed: 12/17/2023]
Abstract
We develop a quantum embedding method that enables accurate and efficient treatment of interactions between molecules and an environment, while explicitly including many-body correlations. The molecule is composed of classical nuclei and quantum electrons, whereas the environment is modeled via charged quantum harmonic oscillators. We construct a general Hamiltonian and introduce a variational Ansatz for the correlated ground state of the fully interacting molecule-environment system. This wave function is optimized via the variational Monte Carlo method and the ground state energy is subsequently estimated through the diffusion Monte Carlo method. The proposed scheme allows an explicit many-body treatment of electrostatic, polarization, and dispersion interactions between the molecule and the environment. We study solvation energies and excitation energies of benzene derivatives, obtaining excellent agreement with explicit ab initio calculations and experiments.
Collapse
Affiliation(s)
- Matej Ditte
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Matteo Barborini
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Leonardo Medrano Sandonas
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
10
|
Al-Hamdani YS, Zen A, Alfè D. Unraveling H2 chemisorption and physisorption on metal decorated graphene using quantum Monte Carlo. J Chem Phys 2023; 159:204708. [PMID: 38018756 DOI: 10.1063/5.0174232] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2023] [Accepted: 10/24/2023] [Indexed: 11/30/2023] Open
Abstract
Molecular hydrogen has the potential to significantly reduce the use of carbon dioxide emitting energy processes. However, hydrogen gas storage is a major bottleneck for its large-scale use as current storage methods are energy intensive. Among different storage methods, physisorbing molecular hydrogen at ambient pressure and temperatures is a promising alternative-particularly in light of the advancements in tunable lightweight nanomaterials and high throughput screening methods. Nonetheless, understanding hydrogen adsorption in well-defined nanomaterials remains experimentally challenging and reference information is scarce despite the proliferation of works predicting hydrogen adsorption. We focus on Li, Na, Ca, and K, decorated graphene sheets as substrates for molecular hydrogen adsorption, and compute the most accurate adsorption energies available to date using quantum diffusion Monte Carlo (DMC). Building on our previous insights at the density functional theory (DFT) level, we find that a weak covalent chemisorption of molecular hydrogen, known as Kubas interaction, is feasible on Ca decorated graphene according to DMC, in agreement with DFT. This finding is in contrast to previous DMC predictions of the 4H2/Ca+ gas cluster (without graphene) where chemisorption is not favored. However, we find that the adsorption energy of hydrogen on metal decorated graphene according to a widely used DFT method is not fully consistent with DMC. The reference adsorption energies reported herein can be used to find better work-horse methods for application in large-scale modeling of hydrogen adsorption. Furthermore, the implications of this work affect strategies for finding suitable hydrogen storage materials and high-throughput methods.
Collapse
Affiliation(s)
- Yasmine S Al-Hamdani
- Department of Earth Sciences, University College London, London WC1E 6BT, United Kingdom
- Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
- London Centre for Nanotechnology, University College London, London WC1E 6BT, United Kingdom
| | - Andrea Zen
- Department of Earth Sciences, University College London, London WC1E 6BT, United Kingdom
- Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
| | - Dario Alfè
- Department of Earth Sciences, University College London, London WC1E 6BT, United Kingdom
- Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
- London Centre for Nanotechnology, University College London, London WC1E 6BT, United Kingdom
- Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
| |
Collapse
|
11
|
Ichibha T, Nikaido Y, Bennett MC, Krogel JT, Hongo K, Maezono R, Reboredo FA. Locality error free effective core potentials for 3d transition metal elements developed for the diffusion Monte Carlo method. J Chem Phys 2023; 159:164114. [PMID: 37888761 DOI: 10.1063/5.0175381] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Accepted: 10/05/2023] [Indexed: 10/28/2023] Open
Abstract
Pseudopotential locality errors have hampered the applications of the diffusion Monte Carlo (DMC) method in materials containing transition metals, in particular oxides. We have developed locality error free effective core potentials, pseudo-Hamiltonians, for transition metals ranging from Cr to Zn. We have modified a procedure published by some of us in Bennett et al. [J. Chem. Theory Comput. 18, 828 (2022)]. We carefully optimized our pseudo-Hamiltonians and achieved transferability errors comparable to the best semilocal pseudopotentials used with DMC but without incurring in locality errors. Our pseudo-Hamiltonian set (named OPH23) bears the potential to significantly improve the accuracy of many-body-first-principles calculations in fundamental science research of complex materials involving transition metals.
Collapse
Affiliation(s)
- Tom Ichibha
- School of Information Science, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Yutaka Nikaido
- School of Information Science, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
| | - M Chandler Bennett
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Kenta Hongo
- Research Center for Advanced Computing Infrastructure, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
| | - Ryo Maezono
- School of Information Science, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
| | - Fernando A Reboredo
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| |
Collapse
|
12
|
Tyagi R, Zen A, Voora VK. Quantifying the Impact of Halogenation on Intermolecular Interactions and Binding Modes of Aromatic Molecules. J Phys Chem A 2023. [PMID: 37406194 DOI: 10.1021/acs.jpca.3c02291] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/07/2023]
Abstract
Halogenation of aromatic molecules is frequently used to modulate intermolecular interactions with ramifications for optoelectronic and mechanical properties. In this work, we accurately quantify and understand the nature of intermolecular interactions in perhalogenated benzene (PHB) clusters. Using benchmark binding energies from the fixed-node diffusion Monte Carlo (FN-DMC) method, we show that generalized Kohn-Sham semicanonical projected random phase approximation (GKS-spRPA) plus approximate exchange kernel (AKX) provides reliable interaction energies with mean absolute error (MAE) of 0.23 kcal/mol. Using the GKS-spRPA+AXK method, we quantify the interaction energies of several binding modes of PHB clusters ((C6X6)n; X = F, Cl, Br, I; n = 2, 3). For a given binding mode, the interaction energies increase 3-4 times from X = F to X = I; the X-X binding modes have energies in the range of 2-4 kcal/mol, while the π-π binding mode has interaction energies in the range of 4-12 kcal/mol. SAPT-DFT-based energy decomposition analysis is then used to show that the equilibrium geometries are dictated primarily by the dispersion and exchange interactions. Finally, we test the accuracy of several dispersion-corrected density functional approximations and show that only the r2SCAN-D4 method has a low MAE and correct long-range behavior, which makes it suitable for large-scale simulations and for developing structure-function relationships of halogenated aromatic systems.
Collapse
Affiliation(s)
- Ritaj Tyagi
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Andrea Zen
- Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
| | - Vamsee K Voora
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| |
Collapse
|
13
|
Kolesár V, Dubecký M. Accuracy of Noncovalent Interactions Involving d-Elements by the 1-Determinant Fixed-Node Diffusion Monte Carlo Method with Effective Core Potentials. J Chem Theory Comput 2023; 19:1170-1176. [PMID: 36751996 DOI: 10.1021/acs.jctc.2c00872] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Abstract
A critical assessment of effective core potential (ECP)-based single-determinant (SD) fixed-node diffusion quantum Monte Carlo (FNDMC) accuracy in prototypical noncovalent closed-shell systems involving d-elements is presented. Careful analysis of biases and elimination of possible bias sources leads to two findings of practical importance for SD FNDMC in these systems. First, in some systems (HCu:HCu, HCu:CuH), SD FNDMC reveals large biases of interaction energy differences (significantly exceeding the target 2% relative error) vs a reliable coupled-cluster CCSD(T)/CBS (complete basis set) reference. Second, the leading error of SD FNDMC with ECPs was attributed to a higher nuclear charge Z of d-group (pseudo) atoms, when compared to sp elements, in line with a previously reported finding that aggregate SD FNDMC bias tends to increase in systems with higher electronic densities. Therefore, SD FNDMC should only be used with caution in systems with a large Z.
Collapse
Affiliation(s)
- Vladimír Kolesár
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
| | - Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic.,ATRI, Faculty of Materials Science and Technology in Trnava, Slovak University of Technology in Bratislava, J. Bottu 25, 917 24 Trnava, Slovakia
| |
Collapse
|
14
|
Dubecký M, Minárik S, Karlický F. Benchmarking fundamental gap of Sc 2C(OH) 2 MXene by many-body methods. J Chem Phys 2023; 158:054703. [PMID: 36754808 DOI: 10.1063/5.0140315] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
Sc2C(OH)2 is a prototypical non-magnetic member of MXenes, a promising transition-metal-based 2D material family, with a direct bandgap. We provide here a benchmark of its fundamental gap Δ obtained from many-body GW and fixed-node diffusion Monte Carlo methods. Both approaches independently arrive at a similar value of Δ ∼ 1.3 eV, suggesting the validity of both methods. Such a bandgap makes Sc2C(OH)2 a 2D semiconductor suitable for optoelectronic applications. The absorbance spectra and the first exciton binding energy (0.63 eV), based on the Bethe-Salpeter equation, are presented as well. The reported results may serve to delineate experimental uncertainties and enable selection of reasonable approximations such as density functional theory functionals, for use in modeling of related MXenes.
Collapse
Affiliation(s)
- Matúš Dubecký
- Department of Physics, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
| | - Stanislav Minárik
- ATRI, Slovak University of Technology in Bratislava, J. Bottu 25, 917 24 Trnava, Slovakia
| | - František Karlický
- Department of Physics, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
| |
Collapse
|
15
|
Wines D, Choudhary K, Tavazza F. Systematic DFT+U and Quantum Monte Carlo Benchmark of Magnetic Two-Dimensional (2D) CrX 3 (X = I, Br, Cl, F). THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:10.1021/acs.jpcc.2c06733. [PMID: 36727030 PMCID: PMC9888057 DOI: 10.1021/acs.jpcc.2c06733] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
The search for two-dimensional (2D) magnetic materials has attracted a great deal of attention because of the experimental synthesis of 2D CrI3, which has a measured Curie temperature of 45 K. Often times, these monolayers have a higher degree of electron correlation and require more sophisticated methods beyond density functional theory (DFT). Diffusion Monte Carlo (DMC) is a correlated electronic structure method that has been demonstrated to be successful for calculating the electronic and magnetic properties of a wide variety of 2D and bulk systems, since it has a weaker dependence on the Hubbard parameter (U) and density functional. In this study, we designed a workflow that combines DFT +U and DMC in order to treat 2D correlated magnetic systems. We chose monolayer CrX3 (X = I, Br, Cl, F), with a stronger focus on CrI3 and CrBr3, as a case study due to the fact that they have been experimentally realized and have a finite critical temperature. With this DFT+U and DMC workflow and the analytical method of Torelli and Olsen, we estimated a maximum value of 43.56 K for the Tc of CrI3 and 20.78 K for the Tc of CrBr3, in addition to analyzing the spin densities and magnetic properties with DMC and DFT+U. We expect that running this workflow for a well-known material class will aid in the future discovery and characterization of lesser known and more complex correlated 2D magnetic materials.
Collapse
Affiliation(s)
- Daniel Wines
- Materials Science and Engineering Division, National Institute of Standards and Technology (NIST), Gaithersburg, Maryland 20899, United States
| | - Kamal Choudhary
- Materials Science and Engineering Division, National Institute of Standards and Technology (NIST), Gaithersburg, Maryland 20899, United States; Theiss Research, La Jolla, California 92037, United States
| | - Francesca Tavazza
- Materials Science and Engineering Division, National Institute of Standards and Technology (NIST), Gaithersburg, Maryland 20899, United States
| |
Collapse
|
16
|
Wang G, Kincaid B, Zhou H, Annaberdiyev A, Bennett MC, Krogel JT, Mitas L. A new generation of effective core potentials from correlated and spin-orbit calculations: selected heavy elements. J Chem Phys 2022; 157:054101. [DOI: 10.1063/5.0087300] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We introduce new correlation consistent effective core potentials (ccECPs) for the elements I, Te, Bi, Ag, Au, Pd, Ir, Mo, and W with $4d$, $5d$, $6s$ and $6p$ valence spaces. These ccECPs are given as a sum of spin-orbit averaged relativistic effective potential (AREP) and effective spin-orbit (SO) terms. The construction involves several steps with increasing refinements from more simple to fully correlated methods. The optimizations are carried out with objective functions that include weighted many-body atomic spectra, norm-conservation criteria, and spin-orbit splittings. Transferability tests involve molecular binding curves of corresponding hydride and oxide dimers. The constructed ccECPs are systematically better and in a few cases on par with previous effective core potential (ECP) tables on all tested criteria and provide a significant increase in accuracy for valence-only calculations with these elements. Our study confirms the importance of the AREP part in determining the overall quality of the ECP even in the presence of sizable spin-orbit effects. The subsequent quantum Monte Carlo (QMC) calculations point out the importance of accurate trial wave functions which in some cases (mid series transition elements) require treatment well beyond single-reference.
Collapse
Affiliation(s)
| | | | - Haihan Zhou
- NC State University, United States of America
| | | | | | - Jaron T. Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, United States of America
| | - Lubos Mitas
- North Carolina State University, United States of America
| |
Collapse
|
17
|
Zhou H, Scemama A, Wang G, Annaberdiyev A, Kincaid B, Caffarel M, Mitas L. A quantum Monte Carlo study of systems with effective core potentials and node nonlinearities. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2021.111402] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
18
|
Bennett MC, Reboredo FA, Mitas L, Krogel JT. High Accuracy Transition Metal Effective Cores for the Many-Body Diffusion Monte Carlo Method. J Chem Theory Comput 2022; 18:828-839. [PMID: 35001633 DOI: 10.1021/acs.jctc.1c00992] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Practical applications of the real-space diffusion Monte Carlo (DMC) method require the removal of core electrons, where currently localization approximations of semilocal potentials are generally used in the projector. Accurate calculations of complex solids and large molecules demand minimizing the impact of approximated atomic cores. Prior works have shown that the errors from such approximations can be sizable in both finite and periodic systems. In this work, we show that a class of differential pseudopotentials, known as pseudo-Hamiltonians, can be constructed for the 3d transition metal atoms, entirely removing the need for any localization scheme in the DMC projector. As a proof of principle, we demonstrate the approach for the case of Co. In order to minimize errors in the pseudo-Hamiltonian at the many-body level, we generalize the recently proposed correlation-consistent pseudopotential generation scheme to successively close semilocal representations of the differential potentials. Our generation scheme successfully produces potentials tailored specifically for real space projector quantum Monte Carlo methods with low error at the many-body level, i.e., with many-body scattering properties very close to relativistic all-electron results. In particular, we show that the agreement with respect to atomic and molecular quantities reach chemical accuracy in many cases─on par with the most accurate semilocal pseudopotentials available. Further, our pseudo-Hamiltonian generation scheme utilizes standard quantum chemistry codes designed only to work with semilocal pseudopotentials, enabling straightforward generation of pseudo-Hamiltonians for additional elements in future works.
Collapse
Affiliation(s)
- M Chandler Bennett
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Fernando A Reboredo
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Lubos Mitas
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| |
Collapse
|
19
|
Wang X, Berkelbach TC. Absorption Spectra of Solids from Periodic Equation-of-Motion Coupled-Cluster Theory. J Chem Theory Comput 2021; 17:6387-6394. [PMID: 34559525 DOI: 10.1021/acs.jctc.1c00692] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present ab initio absorption spectra of six three-dimensional semiconductors and insulators calculated using Gaussian-based periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). The spectra are calculated efficiently by solving a system of linear equations at each frequency, giving access to an energy range of tens of electronvolts without explicit enumeration of excited states. We assess the impact of cost-saving approximations associated with Brillouin zone sampling, frozen orbitals, and the partitioned EOM-CCSD approximation. Although our most converged spectra exhibit line shapes that are in good agreement with experimental spectra, they are uniformly shifted to higher energies by about 1 eV, which is not explained by the remaining finite-size errors. We tentatively attribute this discrepancy to a combination of vibrational effects and the remaining electron correlation, i.e., triple excitations and above.
Collapse
Affiliation(s)
- Xiao Wang
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Timothy C Berkelbach
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States.,Department of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
20
|
Fanta R, Dubecký M. Noncovalent Interactions by the Quantum Monte Carlo Method: Strong Influence of Isotropic Jastrow Cutoff Radii. J Chem Theory Comput 2021; 17:4242-4249. [PMID: 34169721 DOI: 10.1021/acs.jctc.1c00467] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a paradigmatic example of a strong effect of Jastrow cutoff radii setup on the accuracy of noncovalent interaction energy differences within one-determinant Slater-Jastrow fixed-node diffusion Monte Carlo (1FNDMC) simulations using isotropic Jastrow terms and effective-core potentials. Analysis of total energies, absolute and relative errors, and local energy variance of energy differences vs the reference results suggests a simple procedure to marginalize the related biases. The presented data showcase improvements in dispersion-bounded systems within such a 1FNDMC method.
Collapse
Affiliation(s)
- Roman Fanta
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
| | - Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic.,ATRI, Faculty of Materials Science and Technology in Trnava, Slovak University of Technology in Bratislava, J. Bottu 25, 917 24 Trnava, Slovakia
| |
Collapse
|
21
|
Tiihonen J, Clay RC, Krogel JT. Toward quantum Monte Carlo forces on heavier ions: Scaling properties. J Chem Phys 2021; 154:204111. [PMID: 34241166 DOI: 10.1063/5.0052266] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Quantum Monte Carlo (QMC) forces have been studied extensively in recent decades because of their importance with spectroscopic observables and geometry optimization. Here, we benchmark the accuracy and computational cost of QMC forces. The zero-variance zero-bias (ZVZB) force estimator is used in standard variational and diffusion Monte Carlo simulations with mean-field based trial wavefunctions and atomic pseudopotentials. Statistical force uncertainties are obtained with a recently developed regression technique for heavy tailed QMC data [P. Lopez Rios and G. J. Conduit, Phys. Rev. E 99, 063312 (2019)]. By considering selected atoms and dimers with elements ranging from H to Zn (1 ≤ Zeff ≤ 20), we assess the accuracy and the computational cost of ZVZB forces as the effective pseudopotential valence charge, Zeff, increases. We find that the costs of QMC energies and forces approximately follow simple power laws in Zeff. The force uncertainty grows more rapidly, leading to a best case cost scaling relationship of approximately Zeff 6.5(3) for diffusion Monte Carlo. We find that the accessible system size at fixed computational cost scales as Zeff -2, insensitive to model assumptions or the use of the "space warp" variance-reduction technique. Our results predict the practical cost of obtaining forces for a range of materials, such as transition metal oxides where QMC forces have yet to be applied, and underscore the importance of further developing force variance-reduction techniques, particularly for atoms with high Zeff.
Collapse
Affiliation(s)
- Juha Tiihonen
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Raymond C Clay
- Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| |
Collapse
|
22
|
Zhou Q, Liu ZF, Marks TJ, Darancet P. Electronic Structure of Metallophthalocyanines, MPc (M = Fe, Co, Ni, Cu, Zn, Mg) and Fluorinated MPc. J Phys Chem A 2021; 125:4055-4061. [PMID: 33961423 DOI: 10.1021/acs.jpca.0c10766] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We compute the electronic structure and optical excitation energies of metal-free and transition-metal phthalocyanines (H2Pc and MPc for M = Fe, Co, Ni, Cu, Zn, Mg) using density functional theory with optimally tuned range-separated hybrid functionals (OT-RSH). We show that the OT-RSH approach provides photoemission spectra in quantitative agreement with experiments as well as optical band gaps within 10% of their experimental values, capturing the interplay of localized d-states and delocalized π-π* states for these organometallic compounds. We examine the tunability of MPcs and H2Pc through fluorination, resulting in quasi-rigid shifts of the molecular orbital energies by up to 0.7 eV. Our comprehensive data set provides a new computational benchmark for gas-phase phthalocyanines, significantly improving upon other density-functional-theory-based approaches.
Collapse
Affiliation(s)
- Qunfei Zhou
- Materials Research Science and Engineering Center, Northwestern University, Evanston, Illinois 60208, United States.,Center for Nanoscale Materials, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Zhen-Fei Liu
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Tobin J Marks
- Materials Research Science and Engineering Center, Northwestern University, Evanston, Illinois 60208, United States.,Department of Chemistry and Department of Materials Science and Engineering, Northwestern University, Evanston, Illinois 60208, United States
| | - Pierre Darancet
- Center for Nanoscale Materials, Argonne National Laboratory, Argonne, Illinois 60439, United States.,Northwestern Argonne Institute for Science and Engineering, Evanston, Illinois 60208, United States
| |
Collapse
|
23
|
Lau BTG, Knizia G, Berkelbach TC. Regional Embedding Enables High-Level Quantum Chemistry for Surface Science. J Phys Chem Lett 2021; 12:1104-1109. [PMID: 33475362 DOI: 10.1021/acs.jpclett.0c03274] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Compared to common density functionals, ab initio wave function methods can provide greater reliability and accuracy, which could prove useful when modeling adsorbates or defects of otherwise periodic systems. However, the breaking of translational symmetry necessitates large supercells that are often prohibitive for correlated wave function methods. As an alternative, this paper introduces the regional embedding approach, which enables correlated wave function treatments of only a target fragment of interest through small, fragment-localized orbital spaces constructed using a simple overlap criterion. Applications to the adsorption of water on lithium hydride, hexagonal boron nitride, and graphene substrates show that regional embedding combined with focal-point corrections can provide converged CCSD(T) (coupled-cluster) adsorption energies with very small fragment sizes.
Collapse
Affiliation(s)
- Bryan T G Lau
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Gerald Knizia
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Timothy C Berkelbach
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
24
|
Benali A, Shin H, Heinonen O. Quantum Monte Carlo benchmarking of large noncovalent complexes in the L7 benchmark set. J Chem Phys 2020; 153:194113. [DOI: 10.1063/5.0026275] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Anouar Benali
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Hyeondeok Shin
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Olle Heinonen
- Material Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| |
Collapse
|
25
|
Dubecký M, Karlický F, Minárik S, Mitas L. Fundamental gap of fluorographene by many-body GW and fixed-node diffusion Monte Carlo methods. J Chem Phys 2020; 153:184706. [DOI: 10.1063/5.0030952] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
- ATRI, Faculty of Materials Science and Technology in Trnava, Slovak University of Technology in Bratislava, J. Bottu 25, 917 24 Trnava, Slovakia
| | - František Karlický
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
| | - Stanislav Minárik
- ATRI, Faculty of Materials Science and Technology in Trnava, Slovak University of Technology in Bratislava, J. Bottu 25, 917 24 Trnava, Slovakia
| | - Lubos Mitas
- Department of Physics and CHiPS, North Carolina State University, Raleigh, North Carolina 27695, USA
| |
Collapse
|
26
|
Wines D, Saritas K, Ataca C. A first-principles Quantum Monte Carlo study of two-dimensional (2D) GaSe. J Chem Phys 2020; 153:154704. [DOI: 10.1063/5.0023223] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Affiliation(s)
- Daniel Wines
- Department of Physics, University of Maryland Baltimore County, Baltimore, Maryland 21250, USA
| | - Kayahan Saritas
- Department of Applied Physics, Yale University, New Haven, Connecticut 06520, USA
| | - Can Ataca
- Department of Physics, University of Maryland Baltimore County, Baltimore, Maryland 21250, USA
| |
Collapse
|
27
|
Wang G, Annaberdiyev A, Mitas L. Binding and excitations in SixHy molecular systems using quantum Monte Carlo. J Chem Phys 2020; 153:144303. [DOI: 10.1063/5.0022814] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Affiliation(s)
- Guangming Wang
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Abdulgani Annaberdiyev
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Lubos Mitas
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| |
Collapse
|
28
|
Krogel JT, Reboredo FA. Hybridizing pseudo-Hamiltonians and non-local pseudopotentials in diffusion Monte Carlo. J Chem Phys 2020; 153:104111. [DOI: 10.1063/5.0016778] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jaron T. Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Fernando A. Reboredo
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| |
Collapse
|
29
|
Malone FD, Zhang S, Morales MA. Accelerating Auxiliary-Field Quantum Monte Carlo Simulations of Solids with Graphical Processing Units. J Chem Theory Comput 2020; 16:4286-4297. [DOI: 10.1021/acs.jctc.0c00262] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Shuai Zhang
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Miguel A. Morales
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| |
Collapse
|
30
|
Kent PRC, Annaberdiyev A, Benali A, Bennett MC, Landinez Borda EJ, Doak P, Hao H, Jordan KD, Krogel JT, Kylänpää I, Lee J, Luo Y, Malone FD, Melton CA, Mitas L, Morales MA, Neuscamman E, Reboredo FA, Rubenstein B, Saritas K, Upadhyay S, Wang G, Zhang S, Zhao L. QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo. J Chem Phys 2020; 152:174105. [DOI: 10.1063/5.0004860] [Citation(s) in RCA: 43] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- P. R. C. Kent
- Center for Nanophase Materials Sciences Division and Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Abdulgani Annaberdiyev
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Anouar Benali
- Computational Science Division, Argonne National Laboratory, 9700 S. Cass Avenue, Lemont, Illinois 60439, USA
| | - M. Chandler Bennett
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Edgar Josué Landinez Borda
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| | - Peter Doak
- Center for Nanophase Materials Sciences Division and Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Hongxia Hao
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Kenneth D. Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Jaron T. Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Ilkka Kylänpää
- Computational Physics Laboratory, Tampere University, P.O. Box 692, 33014 Tampere, Finland
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Ye Luo
- Computational Science Division, Argonne National Laboratory, 9700 S. Cass Avenue, Lemont, Illinois 60439, USA
| | - Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| | - Cody A. Melton
- Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| | - Lubos Mitas
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Miguel A. Morales
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Fernando A. Reboredo
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Brenda Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Kayahan Saritas
- Department of Applied Physics, Yale University, New Haven, Connecticut 06520, USA
| | - Shiv Upadhyay
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Guangming Wang
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Shuai Zhang
- Laboratory for Laser Energetics, University of Rochester, 250 E River Rd., Rochester, New York 14623, USA
| | - Luning Zhao
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
31
|
Annaberdiyev A, Melton CA, Bennett MC, Wang G, Mitas L. Accurate Atomic Correlation and Total Energies for Correlation Consistent Effective Core Potentials. J Chem Theory Comput 2020; 16:1482-1502. [DOI: 10.1021/acs.jctc.9b00962] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Abdulgani Annaberdiyev
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, United States
| | - Cody A. Melton
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, United States
- Sandia National Laboratories, Albuquerque, New Mexico 87123, United States
| | - M. Chandler Bennett
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, United States
| | - Guangming Wang
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, United States
| | - Lubos Mitas
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, United States
| |
Collapse
|
32
|
Phung QM, Hagai M, Xiong XG, Yanai T. Polarization consistent basis sets using the projector augmented wave method: a renovation brought by PAW into Gaussian basis sets. Phys Chem Chem Phys 2020; 22:27037-27052. [DOI: 10.1039/d0cp05229a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A new family of polarization consistent basis set, combined with the projector augmented wave method, was introduced. The basis sets are compact and have good performance as compared to conventional all-electron basis sets in DFT calculations.
Collapse
Affiliation(s)
- Quan Manh Phung
- Department of Chemistry
- Graduate School of Science
- Nagoya University
- Nagoya
- Japan
| | - Masaya Hagai
- Department of Chemistry
- Graduate School of Science
- Nagoya University
- Nagoya
- Japan
| | - Xiao-Gen Xiong
- Sino-French Institute of Nuclear Engineering and Technology
- Sun Yat-Sen University
- Zhuhai
- China
| | - Takeshi Yanai
- Department of Chemistry
- Graduate School of Science
- Nagoya University
- Nagoya
- Japan
| |
Collapse
|