-
High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation
Authors:
Jan P. Unsleber,
Hongbin Liu,
Leopold Talirz,
Thomas Weymuth,
Maximilian Mörchen,
Adam Grofe,
Dave Wecker,
Christopher J. Stein,
Ajay Panyala,
Bo Peng,
Karol Kowalski,
Matthias Troyer,
Markus Reiher
Abstract:
Quantum chemical calculations on atomistic systems have evolved into a standard approach to study molecular matter. These calculations often involve a significant amount of manual input and expertise although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for e…
▽ More
Quantum chemical calculations on atomistic systems have evolved into a standard approach to study molecular matter. These calculations often involve a significant amount of manual input and expertise although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput lectronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject hem to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all omponents of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow at the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.
△ Less
Submitted 13 April, 2023; v1 submitted 26 November, 2022;
originally announced November 2022.
-
Triple excitations in Green's function coupled cluster solver for studies of strongly correlated systems in the framework of self-energy embedding theory
Authors:
Avijit Shee,
Chia-Nan Yeh,
Bo Peng,
Karol Kowalski,
Dominika Zgid
Abstract:
Embedding theories became important approaches used for accurate calculations of both molecules and solids. In these theories, a small chosen subset of orbitals is treated with an accurate method, called an impurity solver, capable of describing higher correlation effects. Ideally, such a chosen fragment should contain multiple orbitals responsible for the chemical and physical behavior of the com…
▽ More
Embedding theories became important approaches used for accurate calculations of both molecules and solids. In these theories, a small chosen subset of orbitals is treated with an accurate method, called an impurity solver, capable of describing higher correlation effects. Ideally, such a chosen fragment should contain multiple orbitals responsible for the chemical and physical behavior of the compound. Handing a large number of chosen orbitals presents a very significant challenge for the current generation of solvers used in the physics and chemistry community. Here, we develop a Green's function coupled cluster singles doubles and triples (GFCCSDT) solver that can be used for a quantitative description in both molecules and solids. This solver allows us to treat orbital spaces that are inaccessible to other accurate solvers. At the same time, GFCCSDT maintains high accuracy of the resulting self-energy. Moreover, in conjunction with the GFCCSD solver, it allows us to test the systematic convergence of computational studies. Developing the CC family of solvers paves the road to fully systematic Green's function embedding calculations in solids. In this paper, we focus on the investigation of GFCCSDT self-energies for a strongly correlated problem of SrMnO$_3$ solid. Subsequently, we apply this solver to solid MnO showing that an approximate variant of GFCCSDT is capable of yielding a high accuracy orbital resolved spectral function.
△ Less
Submitted 22 November, 2022;
originally announced November 2022.
-
Equation of motion coupled-cluster approach for intrinsic losses in x-ray spectra
Authors:
John J. Rehr,
Fernando D. Vila,
Joshua J. Kas,
Nitzan Y. Hirshberg,
Karol Kowalski,
Bo Peng
Abstract:
We present an equation of motion coupled cluster approach for calculating and understanding intrinsic inelastic losses in core level x-ray absorption spectra (XAS). The method is based on a factorization of the transition amplitude in the time-domain, which leads to a convolution of an effective one-body spectrum and the core-hole spectral function. The spectral function characterizes these losses…
▽ More
We present an equation of motion coupled cluster approach for calculating and understanding intrinsic inelastic losses in core level x-ray absorption spectra (XAS). The method is based on a factorization of the transition amplitude in the time-domain, which leads to a convolution of an effective one-body spectrum and the core-hole spectral function. The spectral function characterizes these losses in terms of shake-up excitations and satellites, and is calculated using a cumulant representation of the core-hole Green's function that includes non-linear corrections. The one-body spectrum also includes orthogonality corrections that enhance the XAS at the edge.
△ Less
Submitted 16 February, 2020; v1 submitted 13 February, 2020;
originally announced February 2020.
-
Massively parallel quantum chemical density matrix renormalization group method
Authors:
Jiří Brabec,
Jan Brandejs,
Karol Kowalski,
Sotiris Xantheas,
Örs Legeza,
Libor Veis
Abstract:
We present, to the best of our knowlegde, the first attempt to exploit the supercomputer platform for quantum chemical density matrix renormalization group (QC-DMRG) calculations. We have developed the parallel scheme based on the in-house MPI global memory library, which combines operator and symmetry sector parallelisms, and tested its performance on three different molecules, all typical candid…
▽ More
We present, to the best of our knowlegde, the first attempt to exploit the supercomputer platform for quantum chemical density matrix renormalization group (QC-DMRG) calculations. We have developed the parallel scheme based on the in-house MPI global memory library, which combines operator and symmetry sector parallelisms, and tested its performance on three different molecules, all typical candidates for QC-DMRG calculations. In case of the largest calculation, which is the nitrogenase FeMo cofactor cluster with the active space comprising 113 electrons in 76 orbitals and bond dimension equal to 6000, our parallel approach scales up to approximately 2000 CPU cores.
△ Less
Submitted 19 June, 2020; v1 submitted 14 January, 2020;
originally announced January 2020.
-
Green's Function Coupled Cluster Simulation of the Near-valence Ionizations of DNA-fragments
Authors:
Bo Peng,
Karol Kowalski,
Ajay Panyala,
Sriram Krishnamoorthy
Abstract:
Accurate description of the ionization process in DNA is crucial to the understanding of the DNA damage under exposure to ionizing radiation, and the exploration of the potential application of DNA strands in nano-electronics. In this work, by employing our recently developed Green's function coupled-cluster (GFCC) library on supercomputing facilities, we have studied the spectral functions of sev…
▽ More
Accurate description of the ionization process in DNA is crucial to the understanding of the DNA damage under exposure to ionizing radiation, and the exploration of the potential application of DNA strands in nano-electronics. In this work, by employing our recently developed Green's function coupled-cluster (GFCC) library on supercomputing facilities, we have studied the spectral functions of several guanine$-$cytosine (G$-$C) base pair structures ([G$-$C]$_n$, $n=1-3$) for the first time in a relatively broad near-valence regime ([-25.0,-5.0] eV) in the coupled-cluster with singles and doubles (CCSD) level. Our focus is to give a preliminary many-body coupled-cluster understanding and guideline of the vertical ionization energy (VIE), spectral profile, and ionization feature changes of these systems as the system size expands in this near-valence regime. The results show that, as the system size expands, even though the lowest VIEs keep decreasing, the changes of spectral function profile and the relative peak positions get unexpectedly smaller. Further analysis of the ionized states associated with the most intensive peak in the spectral functions reveals non-negligible $|2h,1p\rangle$'s in the ionized wave functions of the considered G$-$C base pair systems. The leading $|2h,1p\rangle$'s associated with the main ionizations from the cytosine part of the G$-$C base pairs feature a transition from the intra-base-pair cytosine $π\rightarrowπ^\ast$ excitation to the inter-base-pair electron excitation as the size of G$-$C base pairs expands, which also indicates the minimum quantum region in the many-body calculations of DNA systems.
△ Less
Submitted 11 December, 2019; v1 submitted 11 November, 2019;
originally announced November 2019.
-
Green's function coupled cluster formulations utilizing extended inner excitations
Authors:
Bo Peng,
Karol Kowalski
Abstract:
In this paper we analyze new approximations of the Green's function coupled cluster (GFCC) method where locations of poles are improved by extending the excitation level of inner auxiliary operators. These new GFCC approximations can be categorized as GFCC-i($n,m$) method, where the excitation level of the inner auxiliary operators ($m$) used to describe the ionization potentials and electron affi…
▽ More
In this paper we analyze new approximations of the Green's function coupled cluster (GFCC) method where locations of poles are improved by extending the excitation level of inner auxiliary operators. These new GFCC approximations can be categorized as GFCC-i($n,m$) method, where the excitation level of the inner auxiliary operators ($m$) used to describe the ionization potentials and electron affinities effects in the $N$$-$1 and $N$+1 particle spaces is higher than the excitation level ($n$) used to correlate the ground-state coupled cluster wave function for the $N$-electron system. Furthermore, we reveal the so-called "$n$+1" rule in this category (or the GFCC-i($n$,$n$+1) method), which states that in order to maintain size-extensivity of the Green's function matrix elements, the excitation level of inner auxiliary operators $X_p(ω)$ and $Y_q(ω)$ cannot exceed $n$+1. We also discuss the role of the moments of coupled cluster equations that in a natural way assures these properties. Our implementation in the present study is focused on the first approximation in this GFCC category, i.e. the GFCC-i(2,3) method. As our first practice, we use the GFCC-i(2,3) method to compute the spectral functions for the N$_2$ and CO molecules in the inner and outer valence regimes. In comparison with the GFCCSD results, the computed spectral functions from the GFCC-i(2,3) method exhibit better agreement with the experimental results and other theoretical results, particularly in terms of providing higher resolution of satellite peaks and more accurate relative positions of these satellite peaks with respect to the main peak positions.
△ Less
Submitted 16 November, 2018; v1 submitted 18 September, 2018;
originally announced September 2018.
-
Relativistic ideal Fermi gas at zero temperature and preferred frame
Authors:
K. Kowalski,
J. Rembielinski,
K. A. Smolinski
Abstract:
We discuss the limit T->0 of the relativistic ideal Fermi gas of luxons (particles moving with the speed of light) and tachyons (hypothetical particles faster than light) based on observations of our recent paper: K. Kowalski, J. Rembielinski and K.A. Smolinski, Phys. Rev. D, 76, 045018 (2007). For bradyons this limit is in fact the nonrelativistic one and therefore it is not studied herein.
We discuss the limit T->0 of the relativistic ideal Fermi gas of luxons (particles moving with the speed of light) and tachyons (hypothetical particles faster than light) based on observations of our recent paper: K. Kowalski, J. Rembielinski and K.A. Smolinski, Phys. Rev. D, 76, 045018 (2007). For bradyons this limit is in fact the nonrelativistic one and therefore it is not studied herein.
△ Less
Submitted 17 December, 2007;
originally announced December 2007.
-
Lorentz covariant statistical mechanics and thermodynamics of the relativistic ideal gas and preferred frame
Authors:
K. Kowalski,
J. Rembielinski,
K. A. Smolinski
Abstract:
The Lorentz covariant classical and quantum statistical mechanics and thermodynamics of an ideal relativistic gas of bradyons (particles slower than light), luxons (particles moving with the speed of light) and tachyons (hypothetical particles faster than light) is discussed. The Lorentz covariant formulation is based on the preferred frame approach which among others enables consistent, free of…
▽ More
The Lorentz covariant classical and quantum statistical mechanics and thermodynamics of an ideal relativistic gas of bradyons (particles slower than light), luxons (particles moving with the speed of light) and tachyons (hypothetical particles faster than light) is discussed. The Lorentz covariant formulation is based on the preferred frame approach which among others enables consistent, free of paradoxes description of tachyons. The thermodynamic functions within the covariant approach are obtained both in classical and quantum case.
△ Less
Submitted 18 December, 2007; v1 submitted 17 December, 2007;
originally announced December 2007.