1
|
Pokhilko P, Zgid D. Natural orbitals and two-particle correlators as tools for the analysis of effective exchange couplings in solids. Phys Chem Chem Phys 2023; 25:21267-21279. [PMID: 37548912 DOI: 10.1039/d3cp01975f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/08/2023]
Abstract
Using generalizations of spin-averaged natural orbitals and two-particle charge correlators for solids, we investigate the electronic structure of antiferromagnetic transition-metal oxides with a fully self-consistent, imaginary-time GW method. Our findings disagree with the Goodenough-Kanamori (GK) rules that are commonly used for the qualitative interpretation of such solids. First, we found a strong dependence of the natural orbital occupancies on momenta, contradicting GK assumptions. Second, along the momentum path, the character of natural orbitals changes. In particular, the contributions of oxygen 2s orbitals are important, which has not been considered in the GK rules. To analyze the influence of the electronic correlation on the values of effective exchange coupling constants, we use both natural orbitals and two-particle correlators and show that electronic screening modulates the degree of superexchange by stabilizing the charge-transfer contributions, which greatly affects these coupling constants. Finally, we give a set of predictions and recommendations regarding the use of density functional, Green's function, and wave-function methods for evaluating effective magnetic couplings in molecules and solids.
Collapse
Affiliation(s)
- Pavel Pokhilko
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA.
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
2
|
Shee A, Yeh CN, Peng B, Kowalski K, Zgid D. Triple Excitations in Green's Function Coupled Cluster Solver for Studies of Strongly Correlated Systems in the Framework of Self-Energy Embedding Theory. J Phys Chem Lett 2023; 14:2416-2424. [PMID: 36856741 DOI: 10.1021/acs.jpclett.2c03616] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Embedding theories became important approaches used for accurate calculations of both molecules and solids. In these theories, a small chosen subset of orbitals is treated with an accurate method, called an impurity solver, capable of describing higher correlation effects. Ideally, such a chosen fragment should contain multiple orbitals responsible for the chemical and physical behavior of the compound. Handling a large number of chosen orbitals presents a very significant challenge for the current generation of solvers used in the physics and chemistry community. Here, we develop a Green's function coupled cluster singles doubles and triples (GFCCSDT) solver that can be used for a quantitative description in both molecules and solids. This solver allows us to treat orbital spaces that are inaccessible to other accurate solvers. At the same time, GFCCSDT maintains high accuracy of the resulting self-energy. Moreover, in conjunction with the GFCCSD solver, it allows us to test the systematic convergence of computational studies. Developing the CC family of solvers paves the road to fully systematic Green's function embedding calculations in solids. In this paper, we focus on the investigation of GFCCSDT self-energies for a strongly correlated problem of SrMnO3 solid. Subsequently, we apply this solver to solid MnO showing that an approximate variant of GFCCSDT is capable of yielding a high accuracy orbital resolved spectral function.
Collapse
Affiliation(s)
- Avijit Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Chia-Nan Yeh
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Bo Peng
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Karol Kowalski
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
3
|
Rudshteyn B, Weber JL, Coskun D, Devlaminck PA, Zhang S, Reichman DR, Shee J, Friesner RA. Calculation of Metallocene Ionization Potentials via Auxiliary Field Quantum Monte Carlo: Toward Benchmark Quantum Chemistry for Transition Metals. J Chem Theory Comput 2022; 18:2845-2862. [PMID: 35377642 PMCID: PMC9123894 DOI: 10.1021/acs.jctc.1c01071] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
The accurate ab initio prediction of ionization energies is essential to understanding the electrochemistry of transition metal complexes in both materials science and biological applications. However, such predictions have been complicated by the scarcity of gas phase experimental data, the relatively large size of the relevant molecules, and the presence of strong electron correlation effects. In this work, we apply all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing multideterminant trial wave functions to six metallocene complexes to compare the computed adiabatic and vertical ionization energies with experimental results. We find that ph-AFQMC yields mean absolute errors (MAEs) of 1.69 ± 1.02 kcal/mol for the adiabatic energies and 2.85 ± 1.13 kcal/mol for the vertical energies. We also carry out density functional theory (DFT) calculations using a variety of functionals, which yields MAEs of 3.62-6.98 kcal/mol and 3.31-9.88 kcal/mol, as well as one variant of localized coupled cluster calculations (DLPNO-CCSD(T0) with moderate PNO cutoffs), which has MAEs of 4.96 and 6.08 kcal/mol, respectively. We also test the reliability of DLPNO-CCSD(T0) and DFT on acetylacetonate (acac) complexes for adiabatic energies measured in the same manner experimentally, and we find higher MAEs, ranging from 4.56 to 10.99 kcal/mol (with a different ordering) for DFT and 6.97 kcal/mol for DLPNO-CCSD(T0). Finally, by utilizing experimental solvation energies, we show that accurate reduction potentials in solution for the metallocene series can be obtained from the AFQMC gas phase results.
Collapse
Affiliation(s)
- Benjamin Rudshteyn
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Dilek Coskun
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Pierre A Devlaminck
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 Fifth Avenue, New York, New York 10010, United States
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
4
|
Van Benschoten WZ, Shepherd JJ. Piecewise Interaction Picture Density Matrix Quantum Monte Carlo. J Chem Phys 2022; 156:184107. [DOI: 10.1063/5.0094290] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact $N$-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC method (IP-DMQMC) which overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.
Collapse
|
5
|
Abstract
Quantum embedding schemes are a promising way to extend multireference computations to large molecules with strong correlation effects localized on a small number of atoms. This work introduces a second-order active-space embedding theory [ASET(2)] which improves upon mean-field frozen embedding by treating fragment-environment interactions via an approximate canonical transformation. The canonical transformation employed in ASET(2) is formulated using the driven similarity renormalization group. The ASET(2) scheme is benchmarked on the N═N bond dissociation in pentyldiazene, the S0 to S1 excitation in 1-octene, and the interaction energy of the O2-benzene complex. The ASET(2) explicit treatment of fragment-environment interactions beyond the mean-field level generally improves the accuracy of embedded computations, and it becomes necessary to achieve an accurate description of excitation energies of 1-octene and the singlet-triplet gap of the O2-benzene complex.
Collapse
Affiliation(s)
- Nan He
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chenyang Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
6
|
Cheng HP, Deumens E, Freericks JK, Li C, Sanders BA. Application of Quantum Computing to Biochemical Systems: A Look to the Future. Front Chem 2020; 8:587143. [PMID: 33330375 PMCID: PMC7732423 DOI: 10.3389/fchem.2020.587143] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2020] [Accepted: 10/15/2020] [Indexed: 12/24/2022] Open
Abstract
Chemistry is considered as one of the more promising applications to science of near-term quantum computing. Recent work in transitioning classical algorithms to a quantum computer has led to great strides in improving quantum algorithms and illustrating their quantum advantage. Because of the limitations of near-term quantum computers, the most effective strategies split the work over classical and quantum computers. There is a proven set of methods in computational chemistry and materials physics that has used this same idea of splitting a complex physical system into parts that are treated at different levels of theory to obtain solutions for the complete physical system for which a brute force solution with a single method is not feasible. These methods are variously known as embedding, multi-scale, and fragment techniques and methods. We review these methods and then propose the embedding approach as a method for describing complex biochemical systems, with the parts not only treated with different levels of theory, but computed with hybrid classical and quantum algorithms. Such strategies are critical if one wants to expand the focus to biochemical molecules that contain active regions that cannot be properly explained with traditional algorithms on classical computers. While we do not solve this problem here, we provide an overview of where the field is going to enable such problems to be tackled in the future.
Collapse
Affiliation(s)
- Hai-Ping Cheng
- Quantum Theory Project, Department of Physics, University of Florida, Gainesville, FL, United States
| | - Erik Deumens
- Quantum Theory Project, Department of Physics, University of Florida, Gainesville, FL, United States
| | - James K. Freericks
- Department of Physics, Georgetown University, Washington, DC, United States
| | - Chenglong Li
- Department of Medicinal Chemistry, University of Florida, Gainesville, FL, United States
| | - Beverly A. Sanders
- Department of Computer and Information Science and Engineering, University of Florida, Gainesville, FL, United States
| |
Collapse
|
7
|
Rudshteyn B, Coskun D, Weber JL, Arthur EJ, Zhang S, Reichman DR, Friesner RA, Shee J. Predicting Ligand-Dissociation Energies of 3d Coordination Complexes with Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2020; 16:3041-3054. [DOI: 10.1021/acs.jctc.0c00070] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Benjamin Rudshteyn
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Dilek Coskun
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L. Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Evan J. Arthur
- Schrodinger Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, United States
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, United States
| | - David R. Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
8
|
Dong X, Zgid D, Gull E, Strand HUR. Legendre-spectral Dyson equation solver with super-exponential convergence. J Chem Phys 2020; 152:134107. [DOI: 10.1063/5.0003145] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Xinyang Dong
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Dominika Zgid
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Emanuel Gull
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Hugo U. R. Strand
- Department of Physics, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden
- Center for Computational Quantum Physics, The Flatiron Institute, New York, New York 10010, USA
| |
Collapse
|
9
|
Shee A, Zgid D. Coupled Cluster as an Impurity Solver for Green’s Function Embedding Methods. J Chem Theory Comput 2019; 15:6010-6024. [DOI: 10.1021/acs.jctc.9b00603] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Avijit Shee
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
10
|
Tran LN, Shea JAR, Neuscamman E. Tracking Excited States in Wave Function Optimization Using Density Matrices and Variational Principles. J Chem Theory Comput 2019; 15:4790-4803. [DOI: 10.1021/acs.jctc.9b00351] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Lan Nguyen Tran
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Ho Chi Minh City Institute of Physics, VAST, Ho Chi Minh City 700000, Vietnam
| | - Jacqueline A. R. Shea
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
11
|
On Achieving High Accuracy in Quantum Chemical Calculations of 3d Transition Metal-Containing Systems: A Comparison of Auxiliary-Field Quantum Monte Carlo with Coupled Cluster, Density Functional Theory, and Experiment for Diatomic Molecules. J Chem Theory Comput 2019; 15:2346-2358. [DOI: 10.1021/acs.jctc.9b00083] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
|
12
|
Rusakov AA, Iskakov S, Tran LN, Zgid D. Self-Energy Embedding Theory (SEET) for Periodic Systems. J Chem Theory Comput 2018; 15:229-240. [DOI: 10.1021/acs.jctc.8b00927] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Alexander A. Rusakov
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Sergei Iskakov
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Lan Nguyen Tran
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, United States
- Ho Chi Minh City Institute of Physics, Vietnam Academy of Science and Technology (VAST), Ho Chi Minh City 70000, Vietnam
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, United States
- Center for Computational Quantum Physics, The Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|