1
|
Drwal D, Matousek M, Golub P, Tucholska A, Hapka M, Brabec J, Veis L, Pernal K. Role of Spin Polarization and Dynamic Correlation in Singlet-Triplet Gap Inversion of Heptazine Derivatives. J Chem Theory Comput 2023; 19:7606-7616. [PMID: 37864544 PMCID: PMC10653106 DOI: 10.1021/acs.jctc.3c00781] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Indexed: 10/23/2023]
Abstract
The new generation of proposed light-emitting molecules for organic light-emitting diodes (OLEDs) has raised considerable research interest due to its exceptional feature─a negative singlet-triplet (ST) gap violating Hund's multiplicity rule in the excited S1 and T1 states. We investigate the role of spin polarization in the mechanism of ST gap inversion. Spin polarization is associated with doubly excited determinants of certain types, whose presence in the wave function expansion favors the energy of the singlet state more than that of the triplet. Using a perturbation theory-based model for spin polarization, we propose a simple descriptor for prescreening of candidate molecules with negative ST gaps and prove its usefulness for heptazine-type molecules. Numerical results show that the quantitative effect of spin polarization decreases linearly with the increasing highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) exchange integral. Comparison of single- and multireference coupled-cluster predictions of ST gaps shows that the former methods provide good accuracy by correctly balancing the effects of doubly excited determinants and dynamic correlation. We also show that accurate ST gaps may be obtained using a complete active space model supplemented with dynamic correlation from multireference adiabatic connection theory.
Collapse
Affiliation(s)
- Daria Drwal
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Mikulas Matousek
- 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, 12116 Prague, Czech Republic
| | - Pavlo Golub
- J.
Heyrovský Institute of Physical Chemistry, Academy of Sciences
of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Aleksandra Tucholska
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Michał Hapka
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Jiri Brabec
- 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
| | - Katarzyna Pernal
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| |
Collapse
|
2
|
Hapka M, Krzemińska A, Modrzejewski M, Przybytek M, Pernal K. Efficient Calculation of the Dispersion Energy for Multireference Systems with Cholesky Decomposition: Application to Excited-State Interactions. J Phys Chem Lett 2023; 14:6895-6903. [PMID: 37494637 PMCID: PMC10405273 DOI: 10.1021/acs.jpclett.3c01568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2023] [Accepted: 07/17/2023] [Indexed: 07/28/2023]
Abstract
Accurate and efficient prediction of dispersion interactions in excited-state complexes poses a challenge due to the complex nature of electron correlation effects that need to be simultaneously considered. We propose an algorithm for computing the dispersion energy in nondegenerate ground- or excited-state complexes with arbitrary spin. The algorithm scales with the fifth power of the system size due to employing Cholesky decomposition of Coulomb integrals and a recently developed recursive formula for density response functions of the monomers. As a numerical illustration, we apply the new algorithm in the framework of multiconfigurational symmetry adapted perturbation theory, SAPT(MC), to study interactions in dimers with localized excitons. The SAPT(MC) analysis reveals that the dispersion energy may be the main force stabilizing excited-state dimers.
Collapse
Affiliation(s)
- Michał Hapka
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Agnieszka Krzemińska
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005 Lodz, Poland
| | - Marcin Modrzejewski
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Michał Przybytek
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Katarzyna Pernal
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005 Lodz, Poland
| |
Collapse
|
3
|
Posenitskiy E, Chilkuri VG, Ammar A, Hapka M, Pernal K, Shinde R, Landinez Borda EJ, Filippi C, Nakano K, Kohulák O, Sorella S, de Oliveira Castro P, Jalby W, Ríos PL, Alavi A, Scemama A. TREXIO: A file format and library for quantum chemistry. J Chem Phys 2023; 158:2888846. [PMID: 37144717 DOI: 10.1063/5.0148161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Accepted: 04/04/2023] [Indexed: 05/06/2023] Open
Abstract
TREXIO is an open-source file format and library developed for the storage and manipulation of data produced by quantum chemistry calculations. It is designed with the goal of providing a reliable and efficient method of storing and exchanging wave function parameters and matrix elements, making it an important tool for researchers in the field of quantum chemistry. In this work, we present an overview of the TREXIO file format and library. The library consists of a front-end implemented in the C programming language and two different back-ends: a text back-end and a binary back-end utilizing the hierarchical data format version 5 library, which enables fast read and write operations. It is compatible with a variety of platforms and has interfaces for Fortran, Python, and OCaml programming languages. In addition, a suite of tools have been developed to facilitate the use of the TREXIO format and library, including converters for popular quantum chemistry codes and utilities for validating and manipulating data stored in TREXIO files. The simplicity, versatility, and ease of use of TREXIO make it a valuable resource for researchers working with quantum chemistry data.
Collapse
Affiliation(s)
- Evgeny Posenitskiy
- Laboratoire de Chimie et Physique Quantiques - UMR5626, CNRS/Université Paul Sabatier, Bat. 3R1b4, 118 route de Narbonne, 31062 Toulouse Cedex 09, France
- Qubit Pharmaceuticals, Incubateur Paris Biotech Santé, 24 Rue du Faubourg Saint Jacques, 75014 Paris, France
| | - Vijay Gopal Chilkuri
- Laboratoire de Chimie et Physique Quantiques - UMR5626, CNRS/Université Paul Sabatier, Bat. 3R1b4, 118 route de Narbonne, 31062 Toulouse Cedex 09, France
- Institut des Sciences Moléculaires de Marseille, Service 561, Campus Scientifique de St. Jérôme, Aix Marseille Université, Centrale Marseille 13, 397 Marseille Cedex 20, France
| | - Abdallah Ammar
- Laboratoire de Chimie et Physique Quantiques - UMR5626, CNRS/Université Paul Sabatier, Bat. 3R1b4, 118 route de Narbonne, 31062 Toulouse Cedex 09, France
| | - Michał Hapka
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005 Lodz, Poland
| | - Ravindra Shinde
- MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Edgar Josué Landinez Borda
- MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Claudia Filippi
- MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Kosuke Nakano
- Research and Services Division of Materials Data and Integrated System, National Institute for Materials Science (NIMS), Tsukuba, Ibaraki 305-0047, Japan
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| | - Otto Kohulák
- Laboratoire de Chimie et Physique Quantiques - UMR5626, CNRS/Université Paul Sabatier, Bat. 3R1b4, 118 route de Narbonne, 31062 Toulouse Cedex 09, France
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| | - Sandro Sorella
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| | | | - William Jalby
- Université Paris-Saclay, UVSQ, LI-PaRAD, 9 Boulevard d'Alembert, 78280 Guyancourt, France
| | - Pablo López Ríos
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques - UMR5626, CNRS/Université Paul Sabatier, Bat. 3R1b4, 118 route de Narbonne, 31062 Toulouse Cedex 09, France
| |
Collapse
|
4
|
Abstract
A multiconfigurational adiabatic connection (AC) formalism is an attractive approach to compute the dynamic correlation within the complete active space self-consistent field and density matrix renormalization group (DMRG) models. Practical realizations of AC have been based on two approximations: (i) fixing one- and two-electron reduced density matrices (1- and 2-RDMs) at the zero-coupling constant limit and (ii) extended random phase approximation (ERPA). This work investigates the effect of removing the "fixed-RDM" approximation in AC. The analysis is carried out for two electronic Hamiltonian partitionings: the group product function- and the Dyall Hamiltonians. Exact reference AC integrands are generated from the DMRG full configuration interaction solver. Two AC models are investigated, employing either exact 1- and 2-RDMs or their second-order expansions in the coupling constant in the ERPA equations. Calculations for model molecules indicate that lifting the fixed-RDM approximation is a viable way toward improving the accuracy of existing AC approximations.
Collapse
Affiliation(s)
- Mikuláš Matoušek
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Michał Hapka
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - 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
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005 Lodz, Poland
| |
Collapse
|
5
|
Drwal D, Beran P, Hapka M, Modrzejewski M, Sokół A, Veis L, Pernal K. Efficient Adiabatic Connection Approach for Strongly Correlated Systems: Application to Singlet-Triplet Gaps of Biradicals. J Phys Chem Lett 2022; 13:4570-4578. [PMID: 35580342 PMCID: PMC9150121 DOI: 10.1021/acs.jpclett.2c00993] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2022] [Accepted: 05/09/2022] [Indexed: 06/15/2023]
Abstract
Strong electron correlation can be captured with multireference wave function methods, but an accurate description of the electronic structure requires accounting for the dynamic correlation, which they miss. In this work, a new approach for the correlation energy based on the adiabatic connection (AC) is proposed. The ACn method accounts for terms up to order n in the coupling constant, and it is size-consistent and free from instabilities. It employs the multireference random phase approximation and the Cholesky decomposition technique, leading to a computational cost growing with the fifth power of the system size. Because of the dependence on only one- and two-electron reduced density matrices, ACn is more efficient than existing ab initio multireference dynamic correlation methods. ACn affords excellent results for singlet-triplet gaps of challenging organic biradicals. The development presented in this work opens new perspectives for accurate calculations of systems with dozens of strongly correlated electrons.
Collapse
Affiliation(s)
- Daria Drwal
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Pavel Beran
- 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, Ke Karlovu 2027/3, 12116 Prague 2, Czech Republic
| | - Michał Hapka
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Marcin Modrzejewski
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Adam Sokół
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - 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
| | - Katarzyna Pernal
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| |
Collapse
|
6
|
Jangrouei MR, Krzemińska A, Hapka M, Pastorczak E, Pernal K. Dispersion Interactions in Exciton-Localized States. Theory and Applications to π-π* and n-π* Excited States. J Chem Theory Comput 2022; 18:3497-3511. [PMID: 35587598 PMCID: PMC9202351 DOI: 10.1021/acs.jctc.2c00221] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We address the problem
of intermolecular interaction energy calculations
in molecular complexes with localized excitons. Our focus is on the
correct representation of the dispersion energy. We derive an extended
Casimir-Polder formula for direct computation of this contribution
through second order in the intermolecular interaction operator V̂. An alternative formula, accurate to infinite order
in V̂, is derived within the framework of the
adiabatic connection (AC) theory. We also propose a new parametrization
of the VV10 nonlocal correlation density functional, so that it corrects
the CASSCF energy for the dispersion contribution and can be applied
to excited-state complexes. A numerical investigation is carried out
for benzene, pyridine, and peptide complexes with the local exciton
corresponding to the lowest π–π* or n– π*
states. The extended Casimir-Polder formula is implemented in the
framework of multiconfigurational symmetry-adapted perturbation theory,
SAPT(MC). A SAPT(MC) analysis shows that the creation of a localized
exciton affects mostly the electrostatic component of the interaction
energy of investigated complexes. Nevertheless, the changes in Pauli
repulsion and dispersion energies cannot be neglected. We verify the
performance of several perturbation- and AC-based methods. Best results
are obtained with a range-separated variant of an approximate AC approach
employing extended random phase approximation and CASSCF wave functions.
Collapse
Affiliation(s)
- Mohammad Reza Jangrouei
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005, Lodz, Poland
| | - Agnieszka Krzemińska
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005, Lodz, Poland
| | - Michał Hapka
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093, Warsaw, Poland
| | - Ewa Pastorczak
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005, Lodz, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005, Lodz, Poland
| |
Collapse
|
7
|
Beran P, Matoušek M, Hapka M, Pernal K, Veis L. Density Matrix Renormalization Group with Dynamical Correlation via Adiabatic Connection. J Chem Theory Comput 2021; 17:7575-7585. [PMID: 34762423 DOI: 10.1021/acs.jctc.1c00896] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The quantum chemical version of the density matrix renormalization group (DMRG) method has established itself as one of the methods of choice for calculations of strongly correlated molecular systems. Despite its great ability to capture strong electronic correlation in large active spaces, it is not suitable for computations of dynamical electron correlation. In this work, we present a new approach to the electronic structure problem of strongly correlated molecules, in which DMRG is responsible for a proper description of the strong correlation, whereas dynamical correlation is computed via the recently developed adiabatic connection (AC) technique which requires only up to two-body active space reduced density matrices. We report the encouraging results of this approach on typical candidates for DMRG computations, namely, n-acenes (n = 2 → 7), Fe(II)-porphyrin, and the Fe3S4 cluster.
Collapse
Affiliation(s)
- Pavel Beran
- 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, 121 16 Prague, Czech Republic
| | - Mikuláš Matoušek
- 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, 121 16 Prague, Czech Republic
| | - Michał Hapka
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland.,Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - 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
| |
Collapse
|
8
|
Abstract
![]()
We present a formulation
of the multiconfigurational (MC) wave
function symmetry-adapted perturbation theory (SAPT). The method is
applicable to noncovalent interactions between monomers which require
a multiconfigurational description, in particular when the interacting
system is strongly correlated or in an electronically excited state.
SAPT(MC) is based on one- and two-particle reduced density matrices
of the monomers and assumes the single-exchange approximation for
the exchange energy contributions. Second-order terms are expressed
through response properties from extended random phase approximation
(ERPA). The dispersion components of SAPT(MC) have been introduced
in our previous works [HapkaM.2019, 15, 1016−102730525591; HapkaM.2019, 15, 6712–672331670950]. SAPT(MC) is applied either with generalized valence
bond perfect pairing (GVB) or with complete active space self-consistent
field (CASSCF) treatment of the monomers. We discuss two model multireference
systems: the H2 ··· H2 dimer
in out-of-equilibrium geometries and interaction between the argon
atom and excited state of ethylene. Using the C2H4* ··· Ar complex as an example, we examine second-order
terms arising from negative transitions in the linear response function
of an excited monomer. We demonstrate that the negative-transition
terms must be accounted for to ensure qualitative prediction of induction
and dispersion energies and develop a procedure allowing for their
computation. Factors limiting the accuracy of SAPT(MC) are discussed
in comparison with other second-order SAPT schemes on a data set of
small single-reference dimers.
Collapse
Affiliation(s)
- Michał Hapka
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland.,Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Michał Przybytek
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| |
Collapse
|
9
|
Affiliation(s)
| | - Michał Hapka
- Lodz University of Technology Institute of Physics Lodz Poland
- Faculty of Chemistry University of Warsaw Warsaw Poland
| |
Collapse
|
10
|
Madajczyk K, Żuchowski PS, Brzȩk F, Rajchel Ł, Kȩdziera D, Modrzejewski M, Hapka M. Dataset of noncovalent intermolecular interaction energy curves for 24 small high-spin open-shell dimers. J Chem Phys 2021; 154:134106. [PMID: 33832261 DOI: 10.1063/5.0043793] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We introduce a dataset of 24 interaction energy curves of open-shell noncovalent dimers, referred to as the O24 × 5 dataset. The dataset consists of high-spin dimers up to 11 atoms selected to assure diversity with respect to interaction types: dispersion, electrostatics, and induction. The benchmark interaction energies are obtained at the restricted open-shell CCSD(T) level of theory with complete basis set extrapolation (from aug-cc-pVQZ to aug-cc-pV5Z). We have analyzed the performance of selected wave function methods MP2, CCSD, and CCSD(T) as well as the F12a and F12b variants of coupled-cluster theory. In addition, we have tested dispersion-corrected density functional theory methods based on the PBE exchange-correlation model. The O24 × 5 dataset is a challenge to approximate methods due to the wide range of interaction energy strengths it spans. For the dispersion-dominated and mixed-type subsets, any tested method that does not include the triples contribution yields errors on the order of tens of percent. The electrostatic subset is less demanding with errors that are typically an order of magnitude smaller than the mixed and dispersion-dominated subsets.
Collapse
Affiliation(s)
- Katarzyna Madajczyk
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, ul. Grudzia̧dzka 5, 87-100 Toruń, Poland
| | - Piotr S Żuchowski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, ul. Grudzia̧dzka 5, 87-100 Toruń, Poland
| | - Filip Brzȩk
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, ul. Grudzia̧dzka 5, 87-100 Toruń, Poland
| | - Łukasz Rajchel
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, ul. Grudzia̧dzka 5, 87-100 Toruń, Poland
| | - Dariusz Kȩdziera
- Faculty of Chemistry, Nicolaus Copernicus University, ul. Gagarina 7, Toruń, Poland
| | - Marcin Modrzejewski
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Michał Hapka
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
11
|
Hapka M, Krzemińska A, Pernal K. How Much Dispersion Energy Is Included in the Multiconfigurational Interaction Energy? J Chem Theory Comput 2020; 16:6280-6293. [PMID: 32877179 PMCID: PMC7586340 DOI: 10.1021/acs.jctc.0c00681] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2020] [Indexed: 11/30/2022]
Abstract
We demonstrate how to quantify the amount of dispersion interaction recovered by supermolecular calculations with the multiconfigurational self-consistent field (MCSCF) wave functions. For this purpose, we present a rigorous derivation which connects the portion of dispersion interaction captured by the assumed wave function model-the residual dispersion interaction-with the size of the active space. Based on the obtained expression for the residual dispersion contribution, we propose a dispersion correction for the MCSCF that avoids correlation double counting. Numerical demonstration for model four-electron dimers in both ground and excited states described with the complete active space self-consistent field (CASSCF) reference serves as a proof-of-concept for the method. Accurate results, largely independent of the size of the active space, are obtained. For many-electron systems, routine CASSCF interaction energy calculations recover a tiny fraction of the full second-order dispersion energy. We found that the residual dispersion is non-negligible only for purely dispersion-bound complexes.
Collapse
Affiliation(s)
- Michał Hapka
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Agnieszka Krzemińska
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Katarzyna Pernal
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| |
Collapse
|
12
|
Hapka M, Pernal K, Gritsenko OV. Local Enhancement of Dynamic Correlation in Excited States: Fresh Perspective on Ionicity and Development of Correlation Density Functional Approximation Based on the On-Top Pair Density. J Phys Chem Lett 2020; 11:5883-5889. [PMID: 32589027 PMCID: PMC7467739 DOI: 10.1021/acs.jpclett.0c01616] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2020] [Accepted: 06/26/2020] [Indexed: 06/11/2023]
Abstract
We discuss the interplay between the nondynamic and dynamic electron correlation in excited states from the perspective of the suppression of dynamic correlation (SDC) and enhancement of dynamic correlation (EDC) effects. We reveal that there exists a connection between the ionic character of a wave function and EDC. Following this finding we introduce a quantitative measure of ionicity based solely on local functions without referring to valence bond models. The ability to recognize both the SDC and EDC regions underlies the presented method, named CASΠDFT, combining complete active space (CAS) wave function and density functional theory (DFT) via the on-top pair density (Π) function. We extend this approach to excited states by devising an improved representation of the EDC effect in the correlation functional. The generalized CASΠDFT uses different DFT functionals for ground and excited states. Numerical demonstration for singlet π → π* excitations shows that CASΠDFT offers satisfactory accuracy at a fraction of the cost of the ab initio approaches.
Collapse
Affiliation(s)
- Michał Hapka
- Institute
of Physics, Lodz University of Technology, PL-90-924 Lodz, Poland
- Faculty
of Chemistry, University of Warsaw, 02-093 Warsaw, Poland
| | - Katarzyna Pernal
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Oleg V. Gritsenko
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
- Section
Theoretical Chemistry, VU University, NL-1081 HV Amsterdam, The Netherlands
| |
Collapse
|
13
|
Maradzike E, Hapka M, Pernal K, DePrince AE. Reduced Density Matrix-Driven Complete Active Apace Self-Consistent Field Corrected for Dynamic Correlation from the Adiabatic Connection. J Chem Theory Comput 2020; 16:4351-4360. [PMID: 32538086 DOI: 10.1021/acs.jctc.0c00324] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The recently proposed multireference adiabatic connection (AC) formalism [Pernal, Phys. Rev. Lett. 120, 013001 (2018)] is applied to recover dynamic electron correlation effects lacking in variational two-electron reduced density matrix (v2RDM)-driven complete active space self-consistent field theory (CASSCF). The AC approach is validated by computing potential energy curves for the dissociation of molecular nitrogen and the symmetric double dissociation of H2O while enforcing two sets of approximate N-representability conditions in the underlying v2RDM-driven CASSCF calculations (either two-particle or two-particle plus partial three-particle conditions). The AC yields smaller absolute errors than second-order N-electron perturbation theory (NEVPT2) at all molecular geometries for both sets of the N-representability conditions considered. The efficacy of the approach for thermochemistry is also assessed for a set of 31 small-molecule reactions. When imposing partial three-particle N-representability conditions, mean and maximum unsigned errors in reaction energies from the AC are superior to those from NEVPT2.
Collapse
Affiliation(s)
- Elvis Maradzike
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Michał Hapka
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland.,Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| |
Collapse
|
14
|
Hapka M, Pernal K, Gritsenko OV. Molecular multibond dissociation with small complete active space augmented by correlation density functionals. J Chem Phys 2020; 152:204118. [DOI: 10.1063/5.0009253] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michał Hapka
- Institute of Physics, Lodz University of Technology, PL-90-924 Lodz, Poland
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, PL-90-924 Lodz, Poland
| | - Oleg V. Gritsenko
- Institute of Physics, Lodz University of Technology, PL-90-924 Lodz, Poland
- Section Theoretical Chemistry, VU University, NL-1081 HV Amsterdam, The Netherlands
| |
Collapse
|
15
|
Hapka M, Pastorczak E, Krzemińska A, Pernal K. Long-range-corrected multiconfiguration density functional with the on-top pair density. J Chem Phys 2020; 152:094102. [DOI: 10.1063/1.5138980] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Michał Hapka
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Ewa Pastorczak
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Agnieszka Krzemińska
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| |
Collapse
|
16
|
Hapka M, Jaszuński M. The effect of weak intermolecular interactions on the nuclear magnetic resonance shielding constant in N 2. Magn Reson Chem 2020; 58:245-248. [PMID: 31680305 DOI: 10.1002/mrc.4943] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/07/2019] [Revised: 09/14/2019] [Accepted: 09/17/2019] [Indexed: 06/10/2023]
Abstract
Ab initio calculations are applied to examine the influence of the intermolecular interactions on the shielding constant in gaseous nitrogen. An accurate literature potential energy surface and the nuclear magnetic resonance shielding surface of the N2 -N2 complex calculated in this work provide results in satisfactory agreement with the available experimental estimates of the effect.
Collapse
Affiliation(s)
- Michał Hapka
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093, Warszawa, Poland
| | - Michał Jaszuński
- Institute of Organic Chemistry, Polish Academy of Sciences, Kasprzaka 44, 01-224, Warszawa, Poland
| |
Collapse
|
17
|
Abstract
We present a method for calculation of the second-order exchange-dispersion energy in the framework of the symmetry-adapted perturbation theory (SAPT) for weakly interacting monomers described with multiconfigurational wave functions. The proposed formalism is based on response properties obtained from extended random phase approximation (ERPA) equations and assumes the single-exchange (S2) approximation. The approach is applicable to closed shell systems where static correlation cannot be neglected or to systems in nondegenerate excited states. We examine the new method in combination with either generalized valence bond perfect pairing (GVB) or complete active space self-consistent field (CASSCF) description of the interacting monomers. For model multireference dimers in ground states (H2···H2, Be···Be, He···H2), exchange-dispersion energies are reproduced accurately. For the interaction between the excited hydrogen molecule and the helium atom we found unacceptably large errors which is attributed to the neglect of diagonal double excitations in the employed approximation to the linear response function.
Collapse
Affiliation(s)
- Michał Hapka
- Institute of Physics , Lodz University of Technology , ul. Wolczanska 219 , 90-924 Lodz , Poland.,Faculty of Chemistry , University of Warsaw , ul. L. Pasteura 1 , 02-093 Warsaw , Poland
| | - Michał Przybytek
- Faculty of Chemistry , University of Warsaw , ul. L. Pasteura 1 , 02-093 Warsaw , Poland
| | - Katarzyna Pernal
- Institute of Physics , Lodz University of Technology , ul. Wolczanska 219 , 90-924 Lodz , Poland
| |
Collapse
|
18
|
Pastorczak E, Hapka M, Veis L, Pernal K. Capturing the Dynamic Correlation for Arbitrary Spin-Symmetry CASSCF Reference with Adiabatic Connection Approaches: Insights into the Electronic Structure of the Tetramethyleneethane Diradical. J Phys Chem Lett 2019; 10:4668-4674. [PMID: 31356083 DOI: 10.1021/acs.jpclett.9b01582] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
The recently proposed approach to multireference dynamic correlation energy based on the adiabatic connection (AC) is extended to an arbitrary spin symmetry of the reference state. We show that both the spin-free AC approach and its computationally inexpensive approximation, AC0, when combined with a complete active space wave function, constitute viable alternatives to the perturbation-based and density-functional-based multiconfiguration methods. In particular, the AC0 approach, thanks to its favorable scaling with the system size and the size of the active space, allows for treating larger systems than its perturbation-based counterparts while maintaining comparable accuracy. We show the method's robustness on illustrative chemical systems, including the elusive tetramethyleneethane (TME) diradical, potential energy surfaces of which present a challenge to most computational approaches. For the latter system, AC0 outperforms other methods, staying in close agreement with the full configuration interaction quantum Monte Carlo benchmark. A careful analysis of the contributions to the correlation energy of TME's lowest singlet and triplet states reveals the subtle interplay of the dynamic and static correlation as the key to understanding the shape of the diradical's potential energy surfaces.
Collapse
Affiliation(s)
- Ewa Pastorczak
- Institute of Physics , Lodz University of Technology , ul. Wolczanska 219 , 90-924 Lodz , Poland
| | - Michał Hapka
- Faculty of Chemistry , University of Warsaw , ul. L. Pasteura 1 , 02-093 Warsaw , Poland
| | - 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
| | - Katarzyna Pernal
- Institute of Physics , Lodz University of Technology , ul. Wolczanska 219 , 90-924 Lodz , Poland
| |
Collapse
|
19
|
Abstract
We propose a method for calculating a second-order dispersion energy for weakly interacting multireference systems in arbitrary electronic states. It is based on response properties obtained from extended random phase approximation equations. The introduced formalism is general and requires only one- and two-particle reduced density matrices of monomers. We combine the new method with either generalized valence bond perfect pairing (GVB) or complete active space (CAS) self-consistent field description of the interacting systems. In addition to a general scheme, three approximations, leading to significant reduction of the computational cost, are developed by exploiting Dyall partitioning of the monomer Hamiltonians. For model multireference systems (H2···H2 and Be···Be) the method is accurate, unlike its single-reference-based counterpart. Neither GVB nor CAS description of single-reference monomers improves the dispersion energy with respect to the Hartree-Fock-based results.
Collapse
Affiliation(s)
- Michał Hapka
- Institute of Physics , Lodz University of Technology , ul. Wolczanska 219 , 90-924 Lodz , Poland.,Faculty of Chemistry , University of Warsaw , ul. L. Pasteura 1 , 02-093 Warsaw , Poland
| | - Michał Przybytek
- Faculty of Chemistry , University of Warsaw , ul. L. Pasteura 1 , 02-093 Warsaw , Poland
| | - Katarzyna Pernal
- Institute of Physics , Lodz University of Technology , ul. Wolczanska 219 , 90-924 Lodz , Poland
| |
Collapse
|
20
|
Guńka PA, Hapka M, Hanfland M, Dranka M, Chałasiński G, Zachara J. How and Why Does Helium Permeate Nonporous Arsenolite Under High Pressure? Chemphyschem 2018; 19:857-864. [PMID: 29341365 DOI: 10.1002/cphc.201701156] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2017] [Revised: 12/06/2017] [Indexed: 11/06/2022]
Abstract
Investigations into the helium permeation of arsenolite, the cubic, molecular arsenic(III) oxide polymorph As4 O6 , were carried out to understand how and why arsenolite helium clathrate As4 O6 ⋅2 He is formed. High-pressure synchrotron X-ray diffraction experiments on arsenolite single crystals revealed that the permeation of helium into nonporous arsenolite depends on the time for which the crystal is subjected to high pressure and on the crystal history. The single crystal was totally transformed into As4 O6 ⋅2 He within 45 h under 5 GPa. After release of the pressure, arsenolite was recovered and a repeated increase in pressure up to 3 GPa led to practically instant As4 O6 ⋅2 He formation. However, when a pristine arsenolite single crystal was quickly subjected to a pressure of 13 GPa, no helium permeation was observed at all. No neon permeation was observed in analogous experiments. Quantum mechanical computations indicate that there are no specific attractive interactions between He atoms and As4 O6 molecules at the distances observed in the As4 O6 ⋅2 He crystal structure. Detailed analysis of As4 O6 molecular structure changes has shown that the introduction of He into the arsenolite crystal lattice significantly reduces molecular deformations by decreasing the anisotropy of stress exerted on the As4 O6 molecules. This effect and the pΔV term, rather than any specific As⋅⋅⋅He binding, are the driving forces for the formation As4 O6 ⋅2 He.
Collapse
Affiliation(s)
- Piotr A Guńka
- Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664, Warszawa, Poland
| | - Michał Hapka
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093, Warszawa, Poland
| | - Michael Hanfland
- European Synchrotron Radiation Facility, 71 avenue des Martyrs, CS 40220, 38043, Grenoble Cedex 9, France
| | - Maciej Dranka
- Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664, Warszawa, Poland
| | - Grzegorz Chałasiński
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093, Warszawa, Poland
| | - Janusz Zachara
- Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664, Warszawa, Poland
| |
Collapse
|
21
|
Pastorczak E, Shen J, Hapka M, Piecuch P, Pernal K. Intricacies of van der Waals Interactions in Systems with Elongated Bonds Revealed by Electron-Groups Embedding and High-Level Coupled-Cluster Approaches. J Chem Theory Comput 2017; 13:5404-5419. [DOI: 10.1021/acs.jctc.7b00797] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Ewa Pastorczak
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Jun Shen
- Department
of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Michał Hapka
- Faculty
of Chemistry, University of Warsaw, ul. Pasteura 1, 02-093 Warsaw, Poland
| | - Piotr Piecuch
- Department
of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
- Department
of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, United States
| | - Katarzyna Pernal
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| |
Collapse
|
22
|
Hapka M, Rajchel Ł, Modrzejewski M, Schäffer R, Chałasiński G, Szczęśniak MM. The nature of three-body interactions in DFT: Exchange and polarization effects. J Chem Phys 2017; 147:084106. [PMID: 28863509 DOI: 10.1063/1.4986291] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose a physically motivated decomposition of density functional theory (DFT) 3-body nonadditive interaction energies into the exchange and density-deformation (polarization) components. The exchange component represents the effect of the Pauli exclusion in the wave function of the trimer and is found to be challenging for density functional approximations (DFAs). The remaining density-deformation nonadditivity is less dependent upon the DFAs. Numerical demonstration is carried out for rare gas atom trimers, Ar2-HX (X = F, Cl) complexes, and small hydrogen-bonded and van der Waals molecular systems. None of the tested semilocal, hybrid, and range-separated DFAs properly accounts for the nonadditive exchange in dispersion-bonded trimers. By contrast, for hydrogen-bonded systems, range-separated DFAs achieve a qualitative agreement to within 20% of the reference exchange energy. A reliable performance for all systems is obtained only when the monomers interact through the Hartree-Fock potential in the dispersion-free Pauli blockade scheme. Additionally, we identify the nonadditive second-order exchange-dispersion energy as an important but overlooked contribution in force-field-like dispersion corrections. Our results suggest that range-separated functionals do not include this component, although semilocal and global hybrid DFAs appear to imitate it in the short range.
Collapse
Affiliation(s)
- Michał Hapka
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Łukasz Rajchel
- Faculty of Chemistry, University of Duisburg-Essen, Universitätsstraße 5, 45117 Essen, Germany
| | - Marcin Modrzejewski
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Rainer Schäffer
- Faculty of Chemistry, University of Duisburg-Essen, Universitätsstraße 5, 45117 Essen, Germany
| | - Grzegorz Chałasiński
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | | |
Collapse
|
23
|
Kłos J, Hapka M, Chałasiński G, Halvick P, Stoecklin T. Theoretical study of the buffer-gas cooling and trapping of CrH(X6Σ+) by3He atoms. J Chem Phys 2016; 145:214305. [DOI: 10.1063/1.4968529] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jacek Kłos
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742-2021, USA
| | - Michał Hapka
- Faculty of Chemistry, University of Warsaw, ul. Pasteura 1, 02-093 Warsaw, Poland
| | - Grzegorz Chałasiński
- Faculty of Chemistry, University of Warsaw, ul. Pasteura 1, 02-093 Warsaw, Poland
| | - Philippe Halvick
- Institut des Sciences Moléculaires, Université de Bordeaux, CNRS-UMR 5255, 33405 Talence, France
| | - Thierry Stoecklin
- Institut des Sciences Moléculaires, Université de Bordeaux, CNRS-UMR 5255, 33405 Talence, France
| |
Collapse
|
24
|
Jachymski K, Hapka M, Jankunas J, Osterwalder A. Experimental and Theoretical Studies of Low‐Energy Penning Ionization of NH
3
, CH
3
F, and CHF
3. Chemphyschem 2016; 17:3776-3782. [DOI: 10.1002/cphc.201600608] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2016] [Revised: 08/02/2016] [Indexed: 11/10/2022]
Affiliation(s)
- Krzysztof Jachymski
- Faculty of Physics University of Warsaw Pasteura 5 02-093 Warsaw Poland
- Institute for Theoretical Physics III and Center for Integrated Quantum Science and Technology University of Stuttgart Pfaffenwaldring 57 70550 Stuttgart Germany
| | - Michał Hapka
- Faculty of Chemistry University of Warsaw Pasteura 1 02-093 Warsaw Poland
| | - Justin Jankunas
- Institute for Chemical Sciences and Engineering Ecole Polytechnique Fédérale de Lausanne (EPFL) 1015 Lausanne Switzerland
- Deceased
| | - Andreas Osterwalder
- Institute for Chemical Sciences and Engineering Ecole Polytechnique Fédérale de Lausanne (EPFL) 1015 Lausanne Switzerland
| |
Collapse
|
25
|
Affiliation(s)
- Sirous Yourdkhani
- Department
of Chemistry, Tarbiat Modares University, P.O. Box 14115-175, Tehran, Iran
- Faculty
of Chemistry, University of Warsaw, ul. Pasteura 1, 02-093 Warsaw, Poland
| | - Michał Chojecki
- Faculty
of Chemistry, University of Warsaw, ul. Pasteura 1, 02-093 Warsaw, Poland
| | - Michał Hapka
- Faculty
of Chemistry, University of Warsaw, ul. Pasteura 1, 02-093 Warsaw, Poland
| | - Tatiana Korona
- Faculty
of Chemistry, University of Warsaw, ul. Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
26
|
Jankunas J, Jachymski K, Hapka M, Osterwalder A. Communication: Importance of rotationally inelastic processes in low-energy Penning ionization of CHF3. J Chem Phys 2016; 144:221102. [DOI: 10.1063/1.4953908] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Justin Jankunas
- Institute for Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Krzysztof Jachymski
- Faculty of Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw, Poland
- Institute for Theoretical Physics III and Center for Integrated Quantum Science and Technology, University of Stuttgart, Pfaffenwaldring 57, 70550 Stuttgart, Germany
| | - Michał Hapka
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Andreas Osterwalder
- Institute for Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| |
Collapse
|
27
|
Affiliation(s)
- Justin Jankunas
- Institute for Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | | | - Michał Hapka
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Andreas Osterwalder
- Institute for Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| |
Collapse
|
28
|
Hapka M, Dranka M, Orłowska K, Chałasiński G, Szczęśniak MM, Zachara J. Noncovalent interactions determine the conformation of aurophilic complexes with 2-mercapto-4-methyl-5-thiazoleacetic acid ligands. Dalton Trans 2015; 44:13641-50. [DOI: 10.1039/c5dt01627d] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The synthesis of three new gold(i) salts is reported. The intermolecular interactions stabilizing the aurophilic dimers are examined by several first principles interpretative techniques.
Collapse
Affiliation(s)
| | - Maciej Dranka
- Warsaw University of Technology
- Faculty of Chemistry
- Department of Inorganic Chemistry and Solid State Technology
- 00-664 Warsaw
- Poland
| | - Katarzyna Orłowska
- Warsaw University of Technology
- Faculty of Chemistry
- Department of Inorganic Chemistry and Solid State Technology
- 00-664 Warsaw
- Poland
| | | | | | - Janusz Zachara
- Warsaw University of Technology
- Faculty of Chemistry
- Department of Inorganic Chemistry and Solid State Technology
- 00-664 Warsaw
- Poland
| |
Collapse
|
29
|
Hapka M, Rajchel Ł, Modrzejewski M, Chałasiński G, Szczęśniak MM. Tuned range-separated hybrid functionals in the symmetry-adapted perturbation theory. J Chem Phys 2014; 141:134120. [DOI: 10.1063/1.4896608] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michał Hapka
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Pasteura 1, Poland
| | - Łukasz Rajchel
- Faculty of Chemistry, University of Duisburg-Essen, Universitätsstraße 5, 45117 Essen, Germany
- Interdisciplinary Centre for Mathematical and Computational Modelling, University of Warsaw, 00-838 Warsaw, Prosta 69, Poland
| | - Marcin Modrzejewski
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Pasteura 1, Poland
| | - Grzegorz Chałasiński
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Pasteura 1, Poland
- Department of Chemistry, Oakland University, Rochester, Michigan 48309-4477, USA
| | | |
Collapse
|
30
|
Jankunas J, Bertsche B, Jachymski K, Hapka M, Osterwalder A. Dynamics of gas phase Ne* + NH3 and Ne* + ND3 Penning ionisation at low temperatures. J Chem Phys 2014; 140:244302. [PMID: 24985633 DOI: 10.1063/1.4883517] [Citation(s) in RCA: 75] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Justin Jankunas
- Institute for Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Benjamin Bertsche
- Institute for Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | | | - Michał Hapka
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Andreas Osterwalder
- Institute for Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| |
Collapse
|
31
|
Koppen JV, Hapka M, Modrzejewski M, Szczęśniak MM, Chałasiński G. Density functional theory approach to gold-ligand interactions: Separating true effects from artifacts. J Chem Phys 2014; 140:244313. [DOI: 10.1063/1.4885137] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jessica V. Koppen
- Department of Chemistry, Oakland University, Rochester, Michigan 48309, USA
| | - Michał Hapka
- Faculty of Chemistry, Warsaw University, Pasteura 1, 02-093 Warszawa, Poland
| | - Marcin Modrzejewski
- Faculty of Chemistry, Warsaw University, Pasteura 1, 02-093 Warszawa, Poland
| | | | - Grzegorz Chałasiński
- Department of Chemistry, Oakland University, Rochester, Michigan 48309, USA
- Faculty of Chemistry, Warsaw University, Pasteura 1, 02-093 Warszawa, Poland
| |
Collapse
|
32
|
Hapka M, Kłos J, Korona T, Chałasiński G. Theoretical Studies of Potential Energy Surface and Bound States of the Strongly Bound He(1S)–BeO (1Σ+) Complex. J Phys Chem A 2013; 117:6657-63. [DOI: 10.1021/jp404467b] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Michał Hapka
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Jacek Kłos
- Department of Chemistry and
Biochemistry, University of Maryland, College
Park, Maryland 20742-2021, United States
| | - Tatiana Korona
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Grzegorz Chałasiński
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
- Department of Chemistry, Oakland University, Rochester, Michigan 48309-4477,
United States
| |
Collapse
|
33
|
Hapka M, Chałasiński G, Kłos J, Żuchowski PS. First-principle interaction potentials for metastable He(3S) and Ne(3P) with closed-shell molecules: Application to Penning-ionizing systems. J Chem Phys 2013; 139:014307. [DOI: 10.1063/1.4812182] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
34
|
Hapka M, Żuchowski PS, Szczęśniak MM, Chałasiński G. Symmetry-adapted perturbation theory based on unrestricted Kohn-Sham orbitals for high-spin open-shell van der Waals complexes. J Chem Phys 2013; 137:164104. [PMID: 23126692 DOI: 10.1063/1.4758455] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Two open-shell formulations of the symmetry-adapted perturbation theory are presented. They are based on the spin-unrestricted Kohn-Sham (SAPT(UKS)) and unrestricted Hartree-Fock (SAPT(UHF)) descriptions of the monomers, respectively. The key reason behind development of SAPT(UKS) is that it is more compatible with density functional theory (DFT) compared to the previous formulation of open-shell SAPT based on spin-restricted Kohn-Sham method of Żuchowski et al. [J. Chem. Phys. 129, 084101 (2008)]. The performance of SAPT(UKS) and SAPT(UHF) is tested for the following open-shell van der Waals complexes: He···NH, H(2)O···HO(2), He···OH, Ar···OH, Ar···NO. The results show an excellent agreement between SAPT(UKS) and SAPT(ROKS). Furthermore, for the first time SAPT based on DFT is shown to be suitable for the treatment of interactions involving Π-state radicals (He···OH, Ar···OH, Ar···NO). In the interactions of transition metal dimers ((3)Σ(u)(+))Au(2) and ((13)Σ(g)(+))Cr(2) we show that SAPT is incompatible with the use of effective core potentials. The interaction energies of both systems expressed instead as supermolecular UHF interaction plus dispersion from SAPT(UKS) result in reasonably accurate potential curves.
Collapse
Affiliation(s)
- Michał Hapka
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Pasteura 1, Poland.
| | | | | | | |
Collapse
|
35
|
Koppen JV, Hapka M, Szczęśniak MM, Chałasiński G. Optical absorption spectra of gold clusters Aun(n = 4, 6, 8,12, 20) from long-range corrected functionals with optimal tuning. J Chem Phys 2012; 137:114302. [DOI: 10.1063/1.4752433] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
|
36
|
Rajchel Ł, Żuchowski PS, Hapka M, Modrzejewski M, Szczęśniak MM, Chałasiński G. A density functional theory approach to noncovalent interactions via interacting monomer densities. Phys Chem Chem Phys 2010; 12:14686-92. [DOI: 10.1039/c0cp00626b] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|