1
|
Xiang C, Jia W, Fang WH, Li Z. Distributed Multi-GPU Ab Initio Density Matrix Renormalization Group Algorithm with Applications to the P-Cluster of Nitrogenase. J Chem Theory Comput 2024; 20:775-786. [PMID: 38198503 DOI: 10.1021/acs.jctc.3c01228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
The presence of many degenerate d/f orbitals makes polynuclear transition-metal compounds, such as iron-sulfur clusters in nitrogenase, challenging for state-of-the-art quantum chemistry methods. To address this challenge, we present the first distributed multi-graphics processing unit (GPU) ab initio density matrix renormalization group (DMRG) algorithm suitable for modern high-performance computing (HPC) infrastructures. The central idea is to parallelize the most computationally intensive part─the multiplication of O(K2) operators with a trial wave function, where K is the number of spatial orbitals, by combining operator parallelism for distributing the workload with a batched algorithm for performing contractions on GPU. With this new implementation, we are able to reach an unprecedentedly large bond dimension D = 14,000 on 48 GPUs (NVIDIA A100 80 GB SXM) for an active space model (114 electrons in 73 active orbitals) of the P-cluster, which is nearly 3 times larger than the bond dimensions reported in previous DMRG calculations for the same system using only central processing units (CPUs).
Collapse
Affiliation(s)
- Chunyang Xiang
- State Key Lab of Processors, Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100190, China
| | - Weile Jia
- State Key Lab of Processors, Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100190, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Zhendong Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
2
|
Zhai H, Larsson HR, Lee S, Cui ZH, Zhu T, Sun C, Peng L, Peng R, Liao K, Tölle J, Yang J, Li S, Chan GKL. Block2: A comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond. J Chem Phys 2023; 159:234801. [PMID: 38108484 DOI: 10.1063/5.0180424] [Citation(s) in RCA: 13] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 11/16/2023] [Indexed: 12/19/2023] Open
Abstract
block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matrix renormalization group, such as dynamical correlation theories. The code is designed with an emphasis on flexibility, extensibility, and efficiency and to support integration with external numerical packages. Here, we explain the design principles and currently supported features and present numerical examples in a range of applications.
Collapse
Affiliation(s)
- Huanchen Zhai
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Henrik R Larsson
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Zhi-Hao Cui
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Tianyu Zhu
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Chong Sun
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Linqing Peng
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Ruojing Peng
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Ke Liao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Johannes Tölle
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Junjie Yang
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Shuoxue Li
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
3
|
Guo Y, Zhang N, Liu W. SOiCISCF: Combining SOiCI and iCISCF for Variational Treatment of Spin-Orbit Coupling. J Chem Theory Comput 2023; 19:6668-6685. [PMID: 37728243 DOI: 10.1021/acs.jctc.3c00789] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/21/2023]
Abstract
It has recently been shown that the SOiCI approach [Zhang, N.; J. Phys.: Condens. Matter 2022, 34, 224007], in conjunction with the spin-separated exact two-component relativistic Hamiltonian, can provide very accurate fine structures of systems containing heavy elements by treating electron correlation and spin-orbit coupling (SOC) on an equal footing. Nonetheless, orbital relaxations/polarizations induced by SOC are not yet fully accounted for due to the use of scalar relativistic orbitals. This issue can be resolved by further optimizing the still real-valued orbitals self-consistently in the presence of SOC, as done in the spin-orbit coupled CASSCF approach [Ganyushin, D.; et al. J. Chem. Phys. 2013, 138, 104113] but with the iCISCF algorithm [Guo, Y.; J. Chem. Theory Comput. 2021, 17, 7545-7561] for large active spaces. The resulting SOiCISCF employs both double group and time reversal symmetries for computational efficiency and the assignment of target states. The fine structures of p-block elements are taken as showcases to reveal the efficacy of SOiCISCF.
Collapse
Affiliation(s)
- Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| | - Ning Zhang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| |
Collapse
|
4
|
Nyvang A, Olsen J. A relativistic configuration interaction method with general expansions and initial applications to electronic g-factors. J Chem Phys 2023; 159:044102. [PMID: 37486047 DOI: 10.1063/5.0152655] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Accepted: 06/28/2023] [Indexed: 07/25/2023] Open
Abstract
A new implementation of the orbital-based two-component relativistic configuration interaction approach is reported and applied to calculations of the electronic g-shifts of three diatomic radicals: AlO, HgF, and PdH. The new implementation augments efficient routines for the calculation of nonrelativistic Hamiltonians with new vectorized routines for the calculation of the action of the one-electron spin-orbit operator and allows efficient calculations for the expansion of generalized active space type. The program makes full use of double group as well as time-reversal symmetry. Particle-hole reorganization of the operators is used to improve the efficiency for expansions with nearly fully occupied orbital spaces. The flexibility of the algorithm and program is used to investigate the convergence of electronic g-shifts for the three diatomic radicals as functions of the active space, states included in the orbital optimization, and excitation levels. It was possible to converge to the valence limits within a few percent using expansions containing up to quadruple excitations. However, when excitations from the core orbitals were added, it was not possible to demonstrate convergence to within a few percent with expansions containing at most 10 × 109 determinants.
Collapse
Affiliation(s)
- Andreas Nyvang
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK 8000 Aarhus C, Denmark
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK 8000 Aarhus C, Denmark
| |
Collapse
|
5
|
Zhao L, Zou W. A general method for locating stationary points on the mixed-spin surface of spin-forbidden reaction with multiple spin states. J Chem Phys 2023; 158:2895244. [PMID: 37290081 DOI: 10.1063/5.0151630] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2023] [Accepted: 05/25/2023] [Indexed: 06/10/2023] Open
Abstract
Some chemical reactions proceed on multiple potential energy surfaces and are often accompanied by a change in spin multiplicity, being called spin-forbidden reactions, where the spin-orbit coupling (SOC) effects play a crucial role. In order to efficiently investigate spin-forbidden reactions with two spin states, Yang et al. [Phys. Chem. Chem. Phys. 20, 4129-4136 (2018)] proposed a two-state spin-mixing (TSSM) model, where the SOC effects between the two spin states are simulated by a geometry-independent constant. Inspired by the TSSM model, we suggest a multiple-state spin-mixing (MSSM) model in this paper for the general case with any number of spin states, and its analytic first and second derivatives have been developed for locating stationary points on the mixed-spin potential energy surface and estimating thermochemical energies. To demonstrate the performance of the MSSM model, some spin-forbidden reactions involving 5d transition elements are calculated using the density functional theory (DFT), and the results are compared with the two-component relativistic ones. It is found that MSSM DFT and two-component DFT calculations may provide very similar stationary-point information on the lowest mixed-spin/spinor energy surface, including structures, vibrational frequencies, and zero-point energies. For the reactions containing saturated 5d elements, the reaction energies by MSSM DFT and two-component DFT agree very well within 3 kcal/mol. As for the two reactions OsO+ + CH4 → OOs(CH2)+ + H2 and W + CH4 → WCH2 + H2 involving unsaturated 5d elements, MSSM DFT may also yield good reaction energies of similar accuracy but with some counterexamples. Nevertheless, the energies may be remarkably improved by a posteriori single point energy calculations using two-component DFT at the MSSM DFT optimized geometries, and the maximum error of about 1 kcal/mol is almost independent of the SOC constant used. The MSSM method as well as the developed computer program provides an effective utility for studying spin-forbidden reactions.
Collapse
Affiliation(s)
- Long Zhao
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710127, People's Republic of China
- Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an, Shaanxi 710127, People's Republic of China
| | - Wenli Zou
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710127, People's Republic of China
- Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an, Shaanxi 710127, People's Republic of China
| |
Collapse
|
6
|
Wang X, Sharma S. Relativistic Semistochastic Heat-Bath Configuration Interaction. J Chem Theory Comput 2023; 19:848-855. [PMID: 36700783 DOI: 10.1021/acs.jctc.2c01025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
In this work we present the extension of semistochastic heat-bath configuration interaction (SHCI) to work with any two-component and four-component Hamiltonian. The vertical detachment energy (VDE) of AuH2- and zero-field splitting (ZFS) of NpO22+ is calculated by correlating more than 100 spinors in both cases. This work demonstrates the capability of SHCI to treat problems where both relativistic effect and electron correlation are important.
Collapse
Affiliation(s)
- Xubo Wang
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado80309, United States
| | - Sandeep Sharma
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado80309, United States
| |
Collapse
|
7
|
Li Z. Time-reversal symmetry adaptation in relativistic density matrix renormalization group algorithm. J Chem Phys 2023; 158:044119. [PMID: 36725514 DOI: 10.1063/5.0127621] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
In the nonrelativistic Schrödinger equation, the total spin S and spin projection M are good quantum numbers. In contrast, spin symmetry is lost in the presence of spin-dependent interactions, such as spin-orbit couplings in relativistic Hamiltonians. Therefore, the relativistic density matrix renormalization group algorithm (R-DMRG) only employing particle number symmetry is much more expensive than nonrelativistic DMRG. In addition, artificial breaking of Kramers degeneracy can happen in the treatment of systems with an odd number of electrons. To overcome these issues, we propose time-reversal symmetry adaptation for R-DMRG. Since the time-reversal operator is antiunitary, this cannot be simply achieved in the usual way. We introduce a time-reversal symmetry-adapted renormalized basis and present strategies to maintain the structure of basis functions during the sweep optimization. With time-reversal symmetry adaptation, only half of the renormalized operators are needed, and the computational costs of Hamiltonian-wavefunction multiplication and renormalization are reduced by half. The present construction of the time-reversal symmetry-adapted basis also directly applies to other tensor network states without loops.
Collapse
Affiliation(s)
- Zhendong Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
8
|
Hoyer CE, Lu L, Hu H, Shumilov KD, Sun S, Knecht S, Li X. Correlated Dirac-Coulomb-Breit multiconfigurational self-consistent-field methods. J Chem Phys 2023; 158:044101. [PMID: 36725503 DOI: 10.1063/5.0133741] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
The fully correlated frequency-independent Dirac-Coulomb-Breit Hamiltonian provides the most accurate description of electron-electron interaction before going to a genuine relativistic quantum electrodynamics theory of many-electron systems. In this work, we introduce a correlated Dirac-Coulomb-Breit multiconfigurational self-consistent-field method within the frameworks of complete active space and density matrix renormalization group. In this approach, the Dirac-Coulomb-Breit Hamiltonian is included variationally in both the mean-field and correlated electron treatment. We also analyze the importance of the Breit operator in electron correlation and the rotation between the positive- and negative-orbital space in the no-virtual-pair approximation. Atomic fine-structure splittings and lanthanide contraction in diatomic fluorides are used as benchmark studies to understand the contribution from the Breit correlation.
Collapse
Affiliation(s)
- Chad E Hoyer
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Lixin Lu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Hang Hu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Kirill D Shumilov
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Shichao Sun
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Stefan Knecht
- Algorithmiq Ltd., Kanavakatu 3C, FI-00160 Helsinki, Finland
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
9
|
Liu W. Perspective: Simultaneous treatment of relativity, correlation, and
QED. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1652] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science Shandong University Qingdao Shandong China
| |
Collapse
|
10
|
Kanemaru K, Watanabe Y, Yoshida N, Nakano H. Solvent effects in four-component relativistic electronic structure theory based on the reference interaction-site model. J Comput Chem 2022; 44:5-14. [PMID: 36190170 DOI: 10.1002/jcc.27009] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2022] [Revised: 08/29/2022] [Accepted: 09/04/2022] [Indexed: 11/07/2022]
Abstract
A combined method of the Dirac-Hartree-Fock (DHF) method and the reference interaction-site model (RISM) theory is reported; this is the initial implementation of the coupling of the four-component relativistic electronic structure theory and an integral equation theory of molecular liquids. In the method, the DHF and RISM equations are solved self-consistently, and therefore the electronic structure of the solute, including relativistic effects, and the solvation structure are determined simultaneously. The formulation is constructed based on the variational principle with respect to the Helmholtz energy, and analytic free energy gradients are also derived using the variational property. The method is applied to the iodine ion (I- ), methyl iodide (CH3 I), and hydrogen chalcogenide (H2 X, where X = O-Po) in aqueous solutions, and the electronic structures of the solutes, as well as the solvation free energies and their component analysis, solvent distributions, and solute-solvent interactions, are discussed.
Collapse
Affiliation(s)
- Kodai Kanemaru
- Department of Chemistry, Graduate School of Science, Kyushu University, Fukuoka, Japan
| | - Yoshihiro Watanabe
- Department of Chemistry, Graduate School of Science, Kyushu University, Fukuoka, Japan
| | - Norio Yoshida
- Department of Chemistry, Graduate School of Science, Kyushu University, Fukuoka, Japan.,Department of Complex Systems Science, Graduate School of Informatics, Nagoya University, Nagoya, Japan
| | - Haruyuki Nakano
- Department of Chemistry, Graduate School of Science, Kyushu University, Fukuoka, Japan
| |
Collapse
|
11
|
Hoyer CE, Hu H, Lu L, Knecht S, Li X. Relativistic Kramers-Unrestricted Exact-Two-Component Density Matrix Renormalization Group. J Phys Chem A 2022; 126:5011-5020. [PMID: 35881436 DOI: 10.1021/acs.jpca.2c02150] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this work we develop a variational relativistic density matrix renormalization group (DMRG) approach within the exact-two-component (X2C) framework (X2C-DMRG), using spinor orbitals optimized with the two-component relativistic complete active space self-consistent field. We investigate fine-structure splittings of p- (Ga, In, Tl) and d-block (Sc, Y, La) atoms and excitation energies of monohydride molecules (GeH, SnH, and TlH) with X2C-DMRG calculations using an all-electron relativistic Hamiltonian in a Kramers-unrestricted basis. We find that X2C-DMRG yields accurate 2P and 2D splittings compared to multireference configuration interaction with singles and doubles (MRCISD). We also investigated the degree of symmetry breaking in the atomic multiplets and convergence of electron correlation in the total energies. Symmetry breaking can be large in some cases (∼30 meV); however, increasing the number of renormalized block states m for the DMRG optimization recovers the symmetry breaking by several orders of magnitude. Encouragingly, we find the convergence of electron correlation to be close to MRCISDTQ5 quality. Relativistic X2C-DMRG approaches are important for cases where spin-orbit coupling is significant and the underlying reference wave function requires a large determinantal space. We are able to obtain quantitatively correct fine-structure splittings for systems up to 1019 number of determinants with traditional CI approaches, which are currently unfeasible to converge for the field.
Collapse
Affiliation(s)
- Chad E Hoyer
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hang Hu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Lixin Lu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Stefan Knecht
- Algorithmiq Ltd., Kanavakatu 3C, FI-00160 Helsinki, Finland.,Abteilung SHE Chemie, GSI Helmholtzzentrum für Schwerionenforschung, DE-64291 Darmstadt, Germany.,Department Chemie, Johannes-Gutenberg Universität Mainz, DE-55128 Mainz, Germany
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
12
|
Grofe A, Li X. Relativistic nonorthogonal configuration interaction: application to L 2,3-edge X-ray spectroscopy. Phys Chem Chem Phys 2022; 24:10745-10756. [PMID: 35451435 DOI: 10.1039/d2cp01127a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
In this article, we develop a relativistic exact-two-component nonorthogonal configuration interaction (X2C-NOCI) for computing L-edge X-ray spectra. This article to our knowledge is the first time NOCI has been used for relativistic wave functions. A set of molecular complexes, including SF6, SiCl4 and [FeCl6]3-, are used to demonstrate the accuracy and computational scaling of the X2C-NOCI method. Our results suggest that X2C-NOCI is able to satisfactorily capture the main features of the L2,3-edge X-ray absorption spectra. Excitations from the core require a large amount of orbital relaxation to yield reasonable energies and X2C-NOCI allows us to treat orbital optimization explicitly. However, the cost of computing the nonorthogonal coupling is higher than in conventional CI. Here, we propose an improved integral screening using overlap-scaled density combined with a continuous measure of the generalized Slater-Condon rules that allows us to estimate if an element is zero before attempting a two-electron integral contraction.
Collapse
Affiliation(s)
- Adam Grofe
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA.
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA.
| |
Collapse
|
13
|
Lu L, Hu H, Jenkins AJ, Li X. Exact-Two-Component Relativistic Multireference Second-Order Perturbation Theory. J Chem Theory Comput 2022; 18:2983-2992. [PMID: 35481362 DOI: 10.1021/acs.jctc.2c00171] [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/30/2022]
Abstract
As the relativistic corrections become stronger for late-row elements, the fully perturbative treatment of spin-orbit coupling and dynamic correlation may become inadequate for accurate descriptions of chemical properties. In this work, we introduce a determinant-based Kramers-unrestricted exact-two-component multireference second-order perturbation (X2C-MRPT2) method which variationally includes relativistic corrections with a perturbative dynamic correlation. The restricted active space partitioning scheme is employed to provide an adjustable correlation space for the second-order perturbation treatment. The multistate perturbation theory is also developed to improve the descriptions of ground and excited states. Benchmark studies of atomic fine-structure splittings and spectroscopic constants of molecular monohydrides using X2C-MRPT2 are compared to the other perturbative and variational approaches. The results suggest that X2C-MRPT2 is a highly accurate alternative to the fully variational multireference configuration interaction method at only a small fraction of the computational cost.
Collapse
Affiliation(s)
- Lixin Lu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hang Hu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Andrew J Jenkins
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
14
|
Sharma P, Jenkins AJ, Scalmani G, Frisch MJ, Truhlar DG, Gagliardi L, Li X. Exact-Two-Component Multiconfiguration Pair-Density Functional Theory. J Chem Theory Comput 2022; 18:2947-2954. [PMID: 35384665 DOI: 10.1021/acs.jctc.2c00062] [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/28/2022]
Abstract
Molecules containing late-row elements exhibit large relativistic effects. To account for both relativistic effects and electron correlation in a computationally inexpensive way, we derived a formulation of multiconfiguration pair-density functional theory with the relativistic exact-two-component Hamiltonian (X2C-MC-PDFT). In this new method, relativistic effects are included during variational optimization of a reference wave function by exact-two-component complete active-space self-consistent-field (X2C-CASSCF) theory, followed by an energy evaluation using pair-density functional theory. Benchmark studies of excited-state and ground-state fine-structure splitting of atomic species show that X2C-MC-PDFT can significantly improve the X2C-CASSCF results by introducing additional state-specific electron correlation.
Collapse
Affiliation(s)
- Prachi Sharma
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Andrew J Jenkins
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Giovanni Scalmani
- Gaussian Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Michael J Frisch
- Gaussian Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455, United States
| | - Laura Gagliardi
- Department of Chemistry, University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois 60637, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
15
|
Zhang N, Xiao Y, Liu W. SOiCI and iCISO: combining iterative configuration interaction with spin-orbit coupling in two ways. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:224007. [PMID: 35287124 DOI: 10.1088/1361-648x/ac5db4] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Accepted: 03/14/2022] [Indexed: 06/14/2023]
Abstract
The near-exact iCIPT2 approach for strongly correlated systems of electrons, which stems from the combination of iterative configuration interaction (iCI, an exact solver of full CI) with configuration selection for static correlation and second-order perturbation theory (PT2) for dynamic correlation, is extended to the relativistic domain. In the spirit of spin separation, relativistic effects are treated in two steps: scalar relativity is treated by the infinite-order, spin-free part of the exact two-component (X2C) relativistic Hamiltonian, whereas spin-orbit coupling (SOC) is treated by the first-order, Douglas-Kroll-Hess-like SOC operator derived from the same X2C Hamiltonian. Two possible combinations of iCIPT2 with SOC are considered, i.e., SOiCI and iCISO. The former treats SOC and electron correlation on an equal footing, whereas the latter treats SOC in the spirit of state interaction, by constructing and diagonalizing an effective spin-orbit Hamiltonian matrix in a small number of correlated scalar states. Both double group and time reversal symmetries are incorporated to simplify the computation. Pilot applications reveal that SOiCI is very accurate for the spin-orbit splitting (SOS) of heavy atoms, whereas the computationally very cheap iCISO can safely be applied to the SOS of light atoms and even of systems containing heavy atoms when SOC is largely quenched by ligand fields.
Collapse
Affiliation(s)
- Ning Zhang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, People's Republic of China
| | - Yunlong Xiao
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, People's Republic of China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, People's Republic of China
| |
Collapse
|
16
|
Jenkins AJ, Hu H, Lu L, Frisch MJ, Li X. Two-Component Multireference Restricted Active Space Configuration Interaction for the Computation of L-Edge X-ray Absorption Spectra. J Chem Theory Comput 2021; 18:141-150. [PMID: 34908414 DOI: 10.1021/acs.jctc.1c00564] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
X-ray absorption spectroscopy is a powerful probe of local electronic and nuclear structures, providing insights into chemical processes. The theoretical prediction and interpretation of metal L-edge X-ray absorption spectra are complicated by both relativistic effects, including spin-orbit coupling and the multiconfigurational nature of the states involved. This work details an exact two-component multireference restricted active space (RAS) configuration interaction scheme that uses an exact two-component state-averaged complete active space self-consistent-field method, which includes the spin-orbit coupling in a variational manner, for the accurate description of the electronic structure before using a RAS configuration interaction method to describe the core excited states of the X-ray spectrum. Benchmark calculations are presented for a series of iron-containing complexes, with results showing key features of the spectrum being reproduced, including ligand-to-metal charge transfer and shake-up excitations.
Collapse
Affiliation(s)
- Andrew J Jenkins
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hang Hu
- Molecular Engineering and Sciences Institute, University of Washington, Seattle, Washington 98195, United States
| | - Lixin Lu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Michael J Frisch
- Gaussian Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
17
|
Freitag L, Baiardi A, Knecht S, González L. Simplified State Interaction for Matrix Product State Wave Functions. J Chem Theory Comput 2021; 17:7477-7485. [PMID: 34860525 PMCID: PMC8675135 DOI: 10.1021/acs.jctc.1c00674] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
We present an approximation
to the state-interaction approach for
matrix product state (MPS) wave functions (MPSSI) in a nonorthogonal
molecular orbital basis, first presented by Knecht et al. [J. Chem. Theory Comput.,2016, 28, 5881], that allows for a significant reduction of the computational
cost without significantly compromising its accuracy. The approximation
is well-suited if the molecular orbital basis is close to orthogonality,
and its reliability may be estimated a priori with a single numerical
parameter. For an example of a platinum azide complex, our approximation
offers up to 63-fold reduction in computational time compared to the
original method for wave function overlaps and spin–orbit couplings,
while still maintaining numerical accuracy.
Collapse
Affiliation(s)
- Leon Freitag
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Street 17, 1090 Vienna, Austria
| | - Alberto Baiardi
- Laboratory for Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Stefan Knecht
- GSI Helmholtz Centre for Heavy Ion Research, Planckstr. 1, 64291 Darmstadt, Germany
| | - Leticia González
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Street 17, 1090 Vienna, Austria
| |
Collapse
|
18
|
Khedkar A, Roemelt M. Modern multireference methods and their application in transition metal chemistry. Phys Chem Chem Phys 2021; 23:17097-17112. [PMID: 34355719 DOI: 10.1039/d1cp02640b] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Transition metal chemistry is a challenging playground for quantum chemical methods owing to the simultaneous presence of static and dynamic electron correlation effects in many systems. Wavefunction based multireference (MR) methods constitute a physically sound and systematically improvable Ansatz to deal with this complexity but suffer from some conceptual difficulties and high computational costs. The latter problem partially arises from the unfavorable scaling of the Full Configuration Interaction (Full-CI) problem which in the majority of MR methods is solved for a subset of the molecular orbital space, the so-called active space. In the last years multiple methods such as modern variants of selected CI, Full-CI Quantum Monte Carlo (FCIQMC) and the density matrix renormalization group (DMRG) have been developed that solve the Full-CI problem approximately for a fraction of the computational cost required by conventional techniques thereby significantly extending the range of applicability of modern MR methods. This perspective review outlines recent advancements in the field of MR electronic structure methods together with the resulting chances and challenges for theoretical studies in the field of transition metal chemistry. In light of its emerging importance a special focus is put on the selection of adequate active spaces and the concomitant development of numerous selection aides in recent years.
Collapse
Affiliation(s)
- Abhishek Khedkar
- Lehrstuhl für theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany.
| | | |
Collapse
|
19
|
Barcza G, Ivády V, Szilvási T, Vörös M, Veis L, Gali Á, Legeza Ö. DMRG on Top of Plane-Wave Kohn-Sham Orbitals: A Case Study of Defected Boron Nitride. J Chem Theory Comput 2021; 17:1143-1154. [PMID: 33435672 DOI: 10.1021/acs.jctc.0c00809] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
In this paper, we analyze the numerical aspects of the inherent multireference density matrix renormalization group (DMRG) calculations on top of the periodic Kohn-Sham density functional theory using the complete active space approach. The potential of the framework is illustrated by studying hexagonal boron nitride nanoflakes embedding a charged single boron vacancy point defect by revealing a vertical energy spectrum with a prominent multireference character. We investigate the consistency of the DMRG energy spectrum from the perspective of sample size, basis size, and active space selection protocol. Results obtained from standard quantum chemical atom-centered basis calculations and plane-wave based counterparts show excellent agreement. Furthermore, we also discuss the spectrum of the periodic sheet which is in good agreement with extrapolated data of finite clusters. These results pave the way toward applying the DMRG method in extended correlated solid-state systems, such as point defect qubit in wide band gap semiconductors.
Collapse
Affiliation(s)
- Gergely Barcza
- Wigner Research Centre for Physics, P.O. Box 49, Budapest H-1525, Hungary.,J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Prague CZ-18223, Czechia
| | - Viktor Ivády
- Wigner Research Centre for Physics, P.O. Box 49, Budapest H-1525, Hungary.,Department of Physics, Chemistry and Biology, Linköping University, Linköping SE-581 83, Sweden
| | - Tibor Szilvási
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States.,Department of Chemical and Biological Engineering, The University of Alabama, Tuscaloosa, Alabama 35487, United States
| | - Márton Vörös
- Material Sciences Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Prague CZ-18223, Czechia
| | - Ádám Gali
- Wigner Research Centre for Physics, P.O. Box 49, Budapest H-1525, Hungary.,Department of Atomic Physics, Budapest University of Technology and Economics, Budapest H-1111, Hungary
| | - Örs Legeza
- Wigner Research Centre for Physics, P.O. Box 49, Budapest H-1525, Hungary
| |
Collapse
|
20
|
Abstract
We present a Perspective on what the future holds for full configuration interaction (FCI) theory, with an emphasis on conceptual rather than technical details. Upon revisiting the early history of FCI, a number of its key contemporary approximations are compared on as equal a footing as possible, using a recent blind challenge on the benzene molecule as a testbed [Eriksen et al., J. Phys. Chem. Lett., 2020 11, 8922]. In the process, we review the scope of applications for which FCI continues to prove indispensable, and the required traits in terms of robustness, efficacy, and reliability its modern approximations must satisfy are discussed. We close by conveying a number of general observations on the merits offered by the state-of-the-art alongside some of the challenges still faced to this day. While the field has altogether seen immense progress over the years-the past decade, in particular-it remains clear that our community as a whole has a substantial way to go in enhancing the overall applicability of near-exact electronic structure theory for systems of general composition and increasing size.
Collapse
Affiliation(s)
- Janus J Eriksen
- School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
21
|
Brabec J, Brandejs J, Kowalski K, Xantheas S, Legeza Ö, Veis L. Massively parallel quantum chemical density matrix renormalization group method. J Comput Chem 2020; 42:534-544. [DOI: 10.1002/jcc.26476] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2020] [Revised: 11/22/2020] [Accepted: 11/26/2020] [Indexed: 01/19/2023]
Affiliation(s)
- Jiri Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic Prague Czech Republic
| | - Jan Brandejs
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic Prague Czech Republic
- Faculty of Mathematics and Physics Charles University Prague Czech Republic
| | - Karol Kowalski
- Pacific Northwest National Laboratory Richland Washington USA
| | | | - Örs Legeza
- Strongly Correlated Systems “Lendület” Research group, Wigner Research Centre for Physics Budapest Hungary
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic Prague Czech Republic
| |
Collapse
|
22
|
Anderson RJ, Booth GH. Four-component full configuration interaction quantum Monte Carlo for relativistic correlated electron problems. J Chem Phys 2020; 153:184103. [DOI: 10.1063/5.0029863] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Robert J. Anderson
- Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom
| | - George H. Booth
- Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom
| |
Collapse
|
23
|
Wehrli D, Génévriez M, Knecht S, Reiher M, Merkt F. Complete characterization of the 3p Rydberg complex of a molecular ion: MgAr+. I. Observation of the Mg(3pσ)Ar+ B+ state and determination of its structure and dynamics. J Chem Phys 2020; 153:074310. [DOI: 10.1063/5.0015603] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Dominik Wehrli
- Laboratory of Physical Chemistry, ETH Zurich, CH-8093 Zurich, Switzerland
| | - Matthieu Génévriez
- Laboratory of Physical Chemistry, ETH Zurich, CH-8093 Zurich, Switzerland
| | - Stefan Knecht
- Laboratory of Physical Chemistry, ETH Zurich, CH-8093 Zurich, Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, CH-8093 Zurich, Switzerland
| | - Frédéric Merkt
- Laboratory of Physical Chemistry, ETH Zurich, CH-8093 Zurich, Switzerland
| |
Collapse
|
24
|
Saue T, Bast R, Gomes ASP, Jensen HJA, Visscher L, Aucar IA, Di Remigio R, Dyall KG, Eliav E, Fasshauer E, Fleig T, Halbert L, Hedegård ED, Helmich-Paris B, Iliaš M, Jacob CR, Knecht S, Laerdahl JK, Vidal ML, Nayak MK, Olejniczak M, Olsen JMH, Pernpointner M, Senjean B, Shee A, Sunaga A, van Stralen JNP. The DIRAC code for relativistic molecular calculations. J Chem Phys 2020; 152:204104. [PMID: 32486677 DOI: 10.1063/5.0004844] [Citation(s) in RCA: 118] [Impact Index Per Article: 29.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree-Fock, Kohn-Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.
Collapse
Affiliation(s)
- Trond Saue
- Laboratoire de Chimie et Physique Quantique, UMR 5626 CNRS-Université Toulouse III-Paul Sabatier, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Radovan Bast
- Department of Information Technology, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - André Severo Pereira Gomes
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Hans Jørgen Aa Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Lucas Visscher
- Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, NL-1081HV Amsterdam, The Netherlands
| | - Ignacio Agustín Aucar
- Instituto de Modelado e Innovación Tecnológica, CONICET, and Departamento de Física-Facultad de Ciencias Exactas y Naturales, UNNE, Avda. Libertad 5460, W3404AAS Corrientes, Argentina
| | - Roberto Di Remigio
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Kenneth G Dyall
- Dirac Solutions, 10527 NW Lost Park Drive, Portland, Oregon 97229, USA
| | - Ephraim Eliav
- School of Chemistry, Tel Aviv University, Ramat Aviv, Tel Aviv 69978, Israel
| | - Elke Fasshauer
- Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus, Denmark
| | - Timo Fleig
- Laboratoire de Chimie et Physique Quantique, UMR 5626 CNRS-Université Toulouse III-Paul Sabatier, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Loïc Halbert
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Erik Donovan Hedegård
- Division of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Miroslav Iliaš
- Department of Chemistry, Faculty of Natural Sciences, Matej Bel University, Tajovského 40, 974 01 Banská Bystrica, Slovakia
| | - Christoph R Jacob
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstr. 17, 38106 Braunschweig, Germany
| | - Stefan Knecht
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Jon K Laerdahl
- Department of Microbiology, Oslo University Hospital, Oslo, Norway
| | - Marta L Vidal
- Department of Chemistry, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Malaya K Nayak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Trombay, Mumbai 400085, India
| | - Małgorzata Olejniczak
- Centre of New Technologies, University of Warsaw, S. Banacha 2c, 02-097 Warsaw, Poland
| | - Jógvan Magnus Haugaard Olsen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | | | - Bruno Senjean
- Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, NL-1081HV Amsterdam, The Netherlands
| | - Avijit Shee
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Ayaki Sunaga
- Department of Chemistry, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji-city, Tokyo 192-0397, Japan
| | - Joost N P van Stralen
- Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, NL-1081HV Amsterdam, The Netherlands
| |
Collapse
|
25
|
Affiliation(s)
- Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| |
Collapse
|
26
|
Brandejs J, Višňák J, Veis L, Maté M, Legeza Ö, Pittner J. Toward DMRG-tailored coupled cluster method in the 4c-relativistic domain. J Chem Phys 2020; 152:174107. [DOI: 10.1063/1.5144974] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jan Brandejs
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
- Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic
| | - Jakub Višňák
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
- Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic
- Czech Academic City in Erbil, Yassin Najar Street, Kurani Ankawa, Erbil, Kurdistan, Region of Iraq
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Mihály Maté
- Strongly Correlated Systems “Lendület” Research Group, Institute for Solid State Physics and Optics, MTA Wigner Research Centre for Physics, Konkoly-Thege Miklós út 29-33, H-1121 Budapest, Hungary
- Department of Physics of Complex Systems, Eötvös Loránd University, Pf. 32, H-1518 Budapest, Hungary
| | - Örs Legeza
- Strongly Correlated Systems “Lendület” Research Group, Institute for Solid State Physics and Optics, MTA Wigner Research Centre for Physics, Konkoly-Thege Miklós út 29-33, H-1121 Budapest, Hungary
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
27
|
Hu H, Jenkins AJ, Liu H, Kasper JM, Frisch MJ, Li X. Relativistic Two-Component Multireference Configuration Interaction Method with Tunable Correlation Space. J Chem Theory Comput 2020; 16:2975-2984. [DOI: 10.1021/acs.jctc.9b01290] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hang Hu
- Molecular Engineering and Sciences Institute, University of Washington, Seattle, Washington 98195, United States
| | - Andrew J. Jenkins
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hongbin Liu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Joseph M. Kasper
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Michael J. Frisch
- Gaussian Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
28
|
Baiardi A, Reiher M. The density matrix renormalization group in chemistry and molecular physics: Recent developments and new challenges. J Chem Phys 2020; 152:040903. [DOI: 10.1063/1.5129672] [Citation(s) in RCA: 85] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
29
|
Freitag L, Ma Y, Baiardi A, Knecht S, Reiher M. Approximate Analytical Gradients and Nonadiabatic Couplings for the State-Average Density Matrix Renormalization Group Self-Consistent-Field Method. J Chem Theory Comput 2019; 15:6724-6737. [PMID: 31670947 DOI: 10.1021/acs.jctc.9b00969] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
We present an approximate scheme for analytical gradients and nonadiabatic couplings for calculating state-average density matrix renormalization group self-consistent-field wave function. Our formalism follows closely the state-average complete active space self-consistent-field (SA-CASSCF) ansatz, which employs a Lagrangian, and the corresponding Lagrange multipliers are obtained from a solution of the coupled-perturbed CASSCF (CP-CASSCF) equations. We introduce a definition of the matrix product state (MPS) Lagrange multipliers based on a single-site tensor in a mixed-canonical form of the MPS, such that a sweep procedure is avoided in the solution of the CP-CASSCF equations. We apply our implementation to the optimization of a conical intersection in 1,2-dioxetanone, where we are able to fully reproduce the SA-CASSCF result up to arbitrary accuracy.
Collapse
Affiliation(s)
- Leon Freitag
- Laboratorium für Physikalische Chemie , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| | - Yingjin Ma
- Computer Network Information Center , Chinese Academy of Sciences , Beijing 100190 , China.,Center of Scientific Computing Applications & Research, Chinese Academy of Sciences , Beijing 100190 , China
| | - Alberto Baiardi
- Laboratorium für Physikalische Chemie , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| | - Stefan Knecht
- Laboratorium für Physikalische Chemie , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| | - Markus Reiher
- Laboratorium für Physikalische Chemie , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| |
Collapse
|
30
|
Cheng L. A study of non-iterative triples contributions in relativistic equation-of-motion coupled-cluster calculations using an exact two-component Hamiltonian with atomic mean-field spin-orbit integrals: Application to uranyl and other heavy-element compounds. J Chem Phys 2019; 151:104103. [DOI: 10.1063/1.5113796] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
31
|
Fdez. Galván I, Vacher M, Alavi A, Angeli C, Aquilante F, Autschbach J, Bao JJ, Bokarev SI, Bogdanov NA, Carlson RK, Chibotaru LF, Creutzberg J, Dattani N, Delcey MG, Dong SS, Dreuw A, Freitag L, Frutos LM, Gagliardi L, Gendron F, Giussani A, González L, Grell G, Guo M, Hoyer CE, Johansson M, Keller S, Knecht S, Kovačević G, Källman E, Li Manni G, Lundberg M, Ma Y, Mai S, Malhado JP, Malmqvist PÅ, Marquetand P, Mewes SA, Norell J, Olivucci M, Oppel M, Phung QM, Pierloot K, Plasser F, Reiher M, Sand AM, Schapiro I, Sharma P, Stein CJ, Sørensen LK, Truhlar DG, Ugandi M, Ungur L, Valentini A, Vancoillie S, Veryazov V, Weser O, Wesołowski TA, Widmark PO, Wouters S, Zech A, Zobel JP, Lindh R. OpenMolcas: From Source Code to Insight. J Chem Theory Comput 2019; 15:5925-5964. [DOI: 10.1021/acs.jctc.9b00532] [Citation(s) in RCA: 399] [Impact Index Per Article: 79.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Ignacio Fdez. Galván
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
- Department of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
| | - Morgane Vacher
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Ali Alavi
- Max Planck Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Celestino Angeli
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università di Ferrara, Via Luigi Borsari 46, 44121 Ferrara, Italy
| | - Francesco Aquilante
- Département de Chimie Physique, Université de Genève, 30 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Jie J. Bao
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Sergey I. Bokarev
- Institut für Physik, Universität Rostock, Albert-Einstein-Straße 23-24, 18059 Rostock, Germany
| | - Nikolay A. Bogdanov
- Max Planck Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Rebecca K. Carlson
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Liviu F. Chibotaru
- Theory of Nanomaterials Group, University of Leuven, Celestijnenlaan 200F, Leuven 3001, Belgium
| | - Joel Creutzberg
- Department of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Nike Dattani
- Harvard Smithsonian Center for Astrophysics, Cambridge, Massachusetts 02138, United States
| | - Mickaël G. Delcey
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Sijia S. Dong
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205 A, 69120 Heidelberg, Germany
| | - Leon Freitag
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Luis Manuel Frutos
- Departamento de Química Analítica, Química Física e Ingeniería Química, and Instituto de Investigación Química “Andrés M. del Río”, Universidad de Alcalá, E-28871 Alcalá de Henares, Madrid, Spain
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Frédéric Gendron
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Angelo Giussani
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
- Instituto de Ciencia Molecular, Universitat de València, Apartado 22085, ES-46071 Valencia, Spain
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Gilbert Grell
- Institut für Physik, Universität Rostock, Albert-Einstein-Straße 23-24, 18059 Rostock, Germany
| | - Meiyuan Guo
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Chad E. Hoyer
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Marcus Johansson
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Sebastian Keller
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Stefan Knecht
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Goran Kovačević
- Division of Materials Physics, Ruđer Bošković Institute, P.O.B. 180, Bijenička 54, HR-10002 Zagreb, Croatia
| | - Erik Källman
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Giovanni Li Manni
- Max Planck Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Marcus Lundberg
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Yingjin Ma
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - João Pedro Malhado
- Department of Chemistry, Imperial College London, London SW7 2AZ, United Kingdom
| | - Per Åke Malmqvist
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Stefanie A. Mewes
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205 A, 69120 Heidelberg, Germany
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study (NZIAS), Massey University Albany, Private Bag
102904, Auckland 0632, New Zealand
| | - Jesper Norell
- Department of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Massimo Olivucci
- Department of Biotechnology, Chemistry and Pharmacy, University of Siena, via A. Moro 2, 53100 Siena, Italy
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
- USIAS and Institut de Physique et Chimie des Matériaux de Strasbourg, Université de Strasbourg-CNRS, 67034 Strasbourg, France
| | - Markus Oppel
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Quan Manh Phung
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven 3001, Belgium
| | - Kristine Pierloot
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven 3001, Belgium
| | - Felix Plasser
- Department of Chemistry, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Andrew M. Sand
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Igor Schapiro
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Prachi Sharma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Christopher J. Stein
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Lasse Kragh Sørensen
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Mihkel Ugandi
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Liviu Ungur
- Department of Chemistry, National University of Singapore, 117543 Singapore
| | - Alessio Valentini
- Theoretical Physical Chemistry, Research Unit MolSys, Allée du 6 Août, 11, 4000 Liège, Belgium
| | - Steven Vancoillie
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Valera Veryazov
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Oskar Weser
- Max Planck Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Tomasz A. Wesołowski
- Département de Chimie Physique, Université de Genève, 30 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Per-Olof Widmark
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Sebastian Wouters
- Brantsandpatents, Pauline van Pottelsberghelaan 24, 9051 Sint-Denijs-Westrem, Belgium
| | - Alexander Zech
- Département de Chimie Physique, Université de Genève, 30 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - J. Patrick Zobel
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Roland Lindh
- Department of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
- Uppsala Center for Computational Chemistry (UC3), Uppsala University, P.O. Box 596, SE-751 24 Uppsala, Sweden
| |
Collapse
|
32
|
Brandejs J, Veis L, Szalay S, Barcza G, Pittner J, Legeza Ö. Quantum information-based analysis of electron-deficient bonds. J Chem Phys 2019; 150:204117. [PMID: 31153207 DOI: 10.1063/1.5093497] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
Recently, the correlation theory of the chemical bond was developed, which applies concepts of quantum information theory for the characterization of chemical bonds, based on the multiorbital correlations within the molecule. Here, for the first time, we extend the use of this mathematical toolbox for the description of electron-deficient bonds. We start by verifying the theory on the textbook example of a molecule with three-center two-electron bonds, namely, diborane(6). We then show that the correlation theory of the chemical bond is able to properly describe the bonding situation in more exotic molecules which have been synthesized and characterized only recently, in particular, the diborane molecule with four hydrogen atoms [diborane(4)] and a neutral zerovalent s-block beryllium complex, whose surprising stability was attributed to a strong three-center two-electron π bond stretching across the C-Be-C core. Our approach is of high importance especially in the light of a constant chase after novel compounds with extraordinary properties where the bonding is expected to be unusual.
Collapse
Affiliation(s)
- Jan Brandejs
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Szilárd Szalay
- Strongly Correlated Systems "Lendület" Research Group, Institute for Solid State Physics and Optics, MTA Wigner Research Centre for Physics, Konkoly-Thege Miklós út 29-33, H-1121 Budapest, Hungary
| | - Gergely Barcza
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Örs Legeza
- Strongly Correlated Systems "Lendület" Research Group, Institute for Solid State Physics and Optics, MTA Wigner Research Centre for Physics, Konkoly-Thege Miklós út 29-33, H-1121 Budapest, Hungary
| |
Collapse
|
33
|
Jenkins AJ, Liu H, Kasper JM, Frisch MJ, Li X. Variational Relativistic Two-Component Complete-Active-Space Self-Consistent Field Method. J Chem Theory Comput 2019; 15:2974-2982. [DOI: 10.1021/acs.jctc.9b00011] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Andrew J. Jenkins
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hongbin Liu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Joseph M. Kasper
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Michael J. Frisch
- Gaussian Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
34
|
Łachmańska A, Tecmer P, Legeza Ö, Boguslawski K. Elucidating cation–cation interactions in neptunyl dications using multi-reference ab initio theory. Phys Chem Chem Phys 2019; 21:744-759. [DOI: 10.1039/c8cp04267e] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Understanding the binding mechanism in neptunyl clusters formed due to cation–cation interactions is of crucial importance in nuclear waste reprocessing and related areas of research.
Collapse
Affiliation(s)
- Aleksandra Łachmańska
- Institute of Physics
- Faculty of Physics
- Astronomy and Informatics
- Nicolaus Copernicus University in Toruń
- 87-100 Toruń
| | - Paweł Tecmer
- Institute of Physics
- Faculty of Physics
- Astronomy and Informatics
- Nicolaus Copernicus University in Toruń
- 87-100 Toruń
| | - Örs Legeza
- Strongly Correlated Systems “Lendület” Research Group
- Wigner Research Center for Physics
- H-1525 Budapest
- Hungary
| | - Katharina Boguslawski
- Institute of Physics
- Faculty of Physics
- Astronomy and Informatics
- Nicolaus Copernicus University in Toruń
- 87-100 Toruń
| |
Collapse
|
35
|
Taffet EJ, Scholes GD. Peridinin Torsional Distortion and Bond-Length Alternation Introduce Intramolecular Charge-Transfer and Correlated Triplet Pair Intermediate Excited States. J Phys Chem B 2018; 122:5835-5844. [DOI: 10.1021/acs.jpcb.8b02504] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Elliot J. Taffet
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Gregory D. Scholes
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
36
|
Veis L, Antalík A, Legeza Ö, Alavi A, Pittner J. The Intricate Case of Tetramethyleneethane: A Full Configuration Interaction Quantum Monte Carlo Benchmark and Multireference Coupled Cluster Studies. J Chem Theory Comput 2018; 14:2439-2445. [DOI: 10.1021/acs.jctc.8b00022] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Andrej Antalík
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
- Faculty of Mathematics and Physics, Charles University, 11636 Prague, Czech Republic
| | - Örs Legeza
- Strongly Correlated Systems “Lendület” Research group, Wigner Research Centre for Physics, H-1525 Budapest, Hungary
| | - Ali Alavi
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
- Max Planck Institüt für Festkörperforschung, 70569 Stuttgart, Germany
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
37
|
Battaglia S, Keller S, Knecht S. Efficient Relativistic Density-Matrix Renormalization Group Implementation in a Matrix-Product Formulation. J Chem Theory Comput 2018; 14:2353-2369. [DOI: 10.1021/acs.jctc.7b01065] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Stefano Battaglia
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sebastian Keller
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Stefan Knecht
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
38
|
|
39
|
Freidzon AY, Kurbatov IA, Vovna VI. Ab initio calculation of energy levels of trivalent lanthanide ions. Phys Chem Chem Phys 2018; 20:14564-14577. [DOI: 10.1039/c7cp08366a] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A fully ab initio computational scheme employing CASSCF/XMCQDPT2/SO-CASSCF for the absorption and emission spectra of trivalent lanthanide complexes is presented.
Collapse
Affiliation(s)
- Alexandra Ya. Freidzon
- Photochemistry Center
- Russian Academy of Sciences
- Moscow
- Russia
- National Research Nuclear University MEPhI (Moscow Engineering Physics Institute)
| | | | | |
Collapse
|
40
|
Mussard B, Sharma S. One-Step Treatment of Spin–Orbit Coupling and Electron Correlation in Large Active Spaces. J Chem Theory Comput 2017; 14:154-165. [DOI: 10.1021/acs.jctc.7b01019] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Bastien Mussard
- Department of Chemistry and
Biochemistry, University of Colorado Boulder, Boulder, Colorado 80302, United States
| | - Sandeep Sharma
- Department of Chemistry and
Biochemistry, University of Colorado Boulder, Boulder, Colorado 80302, United States
| |
Collapse
|
41
|
Ronca E, Li Z, Jimenez-Hoyos CA, Chan GKL. Time-Step Targeting Time-Dependent and Dynamical Density Matrix Renormalization Group Algorithms with ab Initio Hamiltonians. J Chem Theory Comput 2017; 13:5560-5571. [DOI: 10.1021/acs.jctc.7b00682] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Enrico Ronca
- Division of Chemistry and
Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Zhendong Li
- Division of Chemistry and
Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Carlos A. Jimenez-Hoyos
- Division of Chemistry and
Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Garnet Kin-Lic Chan
- Division of Chemistry and
Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
42
|
Yanai T, Saitow M, Xiong XG, Chalupský J, Kurashige Y, Guo S, Sharma S. Multistate Complete-Active-Space Second-Order Perturbation Theory Based on Density Matrix Renormalization Group Reference States. J Chem Theory Comput 2017; 13:4829-4840. [DOI: 10.1021/acs.jctc.7b00735] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Takeshi Yanai
- Department
of Theoretical and Computational Molecular Science, Institute for Molecular Science, Okazaki, 444-8585 Aichi Japan
- The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | - Masaaki Saitow
- The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | - Xiao-Gen Xiong
- Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China
| | - Jakub Chalupský
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nám. 2, 16110 Prague 6, Czech Republic
| | - Yuki Kurashige
- Department
of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwake-cho, Sakyou-ku, Kyoto 606-8520, Japan
- Japan Science and Technology Agency, PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Sheng Guo
- Division
of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Sandeep Sharma
- Department
of Chemistry and Biochemistry, University of Colorado, Boulder, Boulder, Colorado 80302, United States
| |
Collapse
|
43
|
Li Z, Chan GKL. Spin-Projected Matrix Product States: Versatile Tool for Strongly Correlated Systems. J Chem Theory Comput 2017; 13:2681-2695. [PMID: 28467847 DOI: 10.1021/acs.jctc.7b00270] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a new wave function ansatz that combines the strengths of spin projection with the language of matrix product states (MPS) and matrix product operators (MPO) as used in the density matrix renormalization group (DMRG). Specifically, spin-projected matrix product states (SP-MPS) are constructed as [Formula: see text], where [Formula: see text] is the spin projector for total spin S and |ΨMPS(N,M)⟩ is an MPS wave function with a given particle number N and spin projection M. This new ansatz possesses several attractive features: (1) It provides a much simpler route to achieve spin adaptation (i.e., to create eigenfunctions of Ŝ2) compared to explicitly incorporating the non-Abelian SU(2) symmetry into the MPS. In particular, since the underlying state |ΨMPS(N,M)⟩ in the SP-MPS uses only Abelian symmetries, one does not need the singlet embedding scheme for nonsinglet states, as normally employed in spin-adapted DMRG, to achieve a single consistent variationally optimized state. (2) Due to the use of |ΨMPS(N,M)⟩ as its underlying state, the SP-MPS can be closely connected to broken-symmetry mean-field states. This allows one to straightforwardly generate the large number of broken-symmetry guesses needed to explore complex electronic landscapes in magnetic systems. Further, this connection can be exploited in the future development of quantum embedding theories for open-shell systems. (3) The sum of MPOs representation for the Hamiltonian and spin projector [Formula: see text] naturally leads to an embarrassingly parallel algorithm for computing expectation values and optimizing SP-MPS. (4) Optimizing SP-MPS belongs to the variation-after-projection (VAP) class of spin-projected theories. Unlike usual spin-projected theories based on determinants, the SP-MPS ansatz can be made essentially exact simply by increasing the bond dimensions in |ΨMPS(N,M)⟩. Computing excited states is also simple by imposing orthogonality constraints, which are simple to implement with MPS. To illustrate the versatility of SP-MPS, we formulate algorithms for the optimization of ground and excited states, develop perturbation theory based on SP-MPS, and describe how to evaluate spin-independent and spin-dependent properties such as the reduced density matrices. We demonstrate the numerical performance of SP-MPS with applications to several models typical of strong correlation, including the Hubbard model, and [2Fe-2S] and [4Fe-4S] model complexes.
Collapse
Affiliation(s)
- Zhendong Li
- Division of Chemistry and Chemical Engineering, California Institute of Technology , Pasadena, California 91125, United States
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology , Pasadena, California 91125, United States
| |
Collapse
|
44
|
|
45
|
Wang X, Valiev RR, Ohulchanskyy TY, Ågren H, Yang C, Chen G. Dye-sensitized lanthanide-doped upconversion nanoparticles. Chem Soc Rev 2017. [DOI: 10.1039/c7cs00053g] [Citation(s) in RCA: 217] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Dye sensitized lanthanide upconversion entails optical upconversion with unprecedented luminescence brightness and broadband excitation wavelength.
Collapse
Affiliation(s)
- Xindong Wang
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage
- School of Chemistry and Chemical Engineering & Key Laboratory of Micro-systems and Micro-structures
- Ministry of Education
- Harbin Institute of Technology
- 150001 Harbin
| | - Rashid R. Valiev
- Department of Theoretical Chemistry and Biology
- Royal Institute of Technology
- S-10691 Stockholm
- Sweden
- Tomsk State University
| | - Tymish Y. Ohulchanskyy
- College of Optoelectronic Engineering
- Key Laboratory of Optoelectronic Devices and Systems of Ministry of Education and Guangdong Province
- Shenzhen University
- 518060 Shenzhen
- People's Republic of China
| | - Hans Ågren
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage
- School of Chemistry and Chemical Engineering & Key Laboratory of Micro-systems and Micro-structures
- Ministry of Education
- Harbin Institute of Technology
- 150001 Harbin
| | - Chunhui Yang
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage
- School of Chemistry and Chemical Engineering & Key Laboratory of Micro-systems and Micro-structures
- Ministry of Education
- Harbin Institute of Technology
- 150001 Harbin
| | - Guanying Chen
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage
- School of Chemistry and Chemical Engineering & Key Laboratory of Micro-systems and Micro-structures
- Ministry of Education
- Harbin Institute of Technology
- 150001 Harbin
| |
Collapse
|
46
|
Knecht S, Keller S, Autschbach J, Reiher M. A Nonorthogonal State-Interaction Approach for Matrix Product State Wave Functions. J Chem Theory Comput 2016; 12:5881-5894. [PMID: 27951678 DOI: 10.1021/acs.jctc.6b00889] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
We present a state-interaction approach for matrix product state (MPS) wave functions in a nonorthogonal molecular orbital basis. Our approach allows us to calculate, for example, transition and spin-orbit coupling matrix elements between arbitrary electronic states, provided that they share the same one-electron basis functions and size of the active orbital space, respectively. The key element is the transformation of the MPS wave functions of different states from a nonorthogonal to a biorthonormal molecular orbital basis representation, by exploiting a sequence of nonunitary transformations, following a proposal by Malmqvist [Int. J. Quantum Chem. 1986, 30, 479]. This is well-known for traditional wave function parametrizations but has not yet been exploited for MPS wave functions.
Collapse
Affiliation(s)
- Stefan Knecht
- Laboratorium für Physikalische Chemie, ETH Zürich , Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sebastian Keller
- Laboratorium für Physikalische Chemie, ETH Zürich , Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York , New York 14260-3000, United States
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich , Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
47
|
Chan GKL, Keselman A, Nakatani N, Li Z, White SR. Matrix product operators, matrix product states, and ab initio density matrix renormalization group algorithms. J Chem Phys 2016; 145:014102. [DOI: 10.1063/1.4955108] [Citation(s) in RCA: 101] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Garnet Kin-Lic Chan
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Anna Keselman
- Department of Condensed Matter Physics, Weizmann Institute of Science, Rehovot 76100, Israel
- Department of Physics and Astronomy, University of California, Irvine, California 92697-4575, USA
| | - Naoki Nakatani
- Institute for Catalysis, Hokkaido University, Kita 21 Nishi 10, Sapporo, Hokkaido 001-0021, Japan
| | - Zhendong Li
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Steven R. White
- Department of Physics and Astronomy, University of California, Irvine, California 92697-4575, USA
| |
Collapse
|
48
|
Sayfutyarova ER, Chan GKL. A state interaction spin-orbit coupling density matrix renormalization group method. J Chem Phys 2016; 144:234301. [DOI: 10.1063/1.4953445] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Garnet Kin-Lic Chan
- Department of Chemistry, Princeton University, Princeton, New Jersey 08540, USA
| |
Collapse
|
49
|
Evidence Supporting a Lymphatic Endothelium Origin for Angiomyolipoma, a TSC2(-) Tumor Related to Lymphangioleiomyomatosis. THE AMERICAN JOURNAL OF PATHOLOGY 2016; 186:1825-1836. [PMID: 27289491 DOI: 10.1016/j.ajpath.2016.03.009] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/05/2016] [Revised: 03/02/2016] [Accepted: 03/07/2016] [Indexed: 12/27/2022]
Abstract
Angiomyolipoma (AML) is a tumor closely related to lymphangioleiomyomatosis (LAM). Both entities are characterized by the proliferation of smooth muscle actin and melanocytic glycoprotein 100 (recognized by antibody HMB-45)-positive spindle-shaped and epithelioid cells. AML and LAM are etiologically linked to mutations in the tsc2 and tsc1 genes in the case of LAM. These genes encode the proteins tuberous sclerosis complex (TSC)-1 and TSC2, which are directly involved in suppressing the mechanistic target of rapamycin cell growth signaling pathway. Although significant progress has been made in characterizing and pharmacologically slowing the progression of AML and LAM with rapamycin, our understanding of their pathogenesis lacks an identified cell of origin. We used an AML-derived cell line to determine whether TSC2 restitution brings about the cell type from which AML arises. We found that AML cells express lymphatic endothelial cell markers consistent with lymphatic endothelial cell precursors in vivo and in vitro. Moreover, on TSC2 correction, AML cells mature into adult lymphatic endothelial cells and have functional attributes characteristic of this cell lineage, suggesting a lymphatic endothelial cell of origin for AML. These effects are dependent on TSC2-mediated mechanistic target of rapamycin inactivation. Finally, we demonstrate the in vitro effectiveness of norcantharidin, a lymphangiogenesis inhibitor, as a potential co-adjuvant therapy in the treatment of AML.
Collapse
|
50
|
Shiozaki T, Mizukami W. Relativistic Internally Contracted Multireference Electron Correlation Methods. J Chem Theory Comput 2015; 11:4733-9. [DOI: 10.1021/acs.jctc.5b00754] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Toru Shiozaki
- Department
of Chemistry, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, United States
| | - Wataru Mizukami
- Department
of Energy and Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka, 816-8580, Japan
| |
Collapse
|