-
MiMiC: A High-Performance Framework for Multiscale Molecular Dynamics Simulations
Authors:
Andrej Antalík,
Andrea Levy,
Sonata Kvedaravičiūtė,
Sophia K. Johnson,
David Carrasco-Busturia,
Bharath Raghavan,
François Mouvet,
Angela Acocella,
Sambit Das,
Vikram Gavini,
Davide Mandelli,
Emiliano Ippoliti,
Simone Meloni,
Paolo Carloni,
Ursula Rothlisberger,
Jógvan Magnus Haugaard Olsen
Abstract:
MiMiC is a framework for performing multiscale simulations, where individual subsystems are handled at different resolutions and/or levels of theory by loosely coupled external programs. To make it highly efficient and flexible, we adopt an interoperable approach based on a multiple-program multiple-data paradigm, serving as an intermediary responsible for fast data exchange and interactions betwe…
▽ More
MiMiC is a framework for performing multiscale simulations, where individual subsystems are handled at different resolutions and/or levels of theory by loosely coupled external programs. To make it highly efficient and flexible, we adopt an interoperable approach based on a multiple-program multiple-data paradigm, serving as an intermediary responsible for fast data exchange and interactions between the subsystems. The main goal of MiMiC is to avoid interfering with the underlying parallelization of the external programs, including the operability on hybrid architectures (e.g., CPU/GPU), and keep their setup and execution as close as possible to the original. At the moment, MiMiC offers an efficient implementation of electrostatic embedding QM/MM that has demonstrated unprecedented parallel scaling in simulations of large biomolecules using CPMD and GROMACS as QM and MM engines, respectively. However, as it is designed for high flexibility with general multiscale models in mind, it can be straightforwardly extended beyond QM/MM. In this article, we illustrate the software design and the features of the framework, which make it a compelling choice for multiscale simulations in the upcoming era of exascale high-performance computing.
△ Less
Submitted 27 March, 2024;
originally announced March 2024.
-
SC1MC-2022: A database of transition metal complexes for training ML models to predict one-site entropies and mutual information
Authors:
Pavlo Golub,
Pavel Beran,
Andrej Antalik,
Jiri Brabec
Abstract:
We introduce a new version of the database SC1MC (SC1MC-2022), obtained by extension of the recent SC1MC-2020, which includes artificial mono transition metal complexes. The database involves reference data used as inputs for training of machine learning models, one- and two-site entropies, and mutual information obtained at the DMRG level for canonical and split-localised orbitals. The purpose of…
▽ More
We introduce a new version of the database SC1MC (SC1MC-2022), obtained by extension of the recent SC1MC-2020, which includes artificial mono transition metal complexes. The database involves reference data used as inputs for training of machine learning models, one- and two-site entropies, and mutual information obtained at the DMRG level for canonical and split-localised orbitals. The purpose of this database is to obtain as much as possible information about the electronic correlation structure, which could be exploited by machine learning models to estimate these important information without a significant computational cost for any similar type of systems with some degree of transferability.
△ Less
Submitted 27 January, 2023; v1 submitted 15 January, 2021;
originally announced January 2021.
-
Automatic selection of active spaces for strongly correlated systems using machine learning algorithms
Authors:
Pavlo Golub,
Andrej Antalik,
Libor Veis,
Jiri Brabec
Abstract:
The active-space quantum chemical methods could provide very accurate description of strongly correlated electronic systems, which is of tremendous value for natural sciences. The proper choice of the active space is crucial, but a non-trivial task. In this article, we present the neural network (NN) based approach for automatic selection of active spaces, focused on transition metal systems. The…
▽ More
The active-space quantum chemical methods could provide very accurate description of strongly correlated electronic systems, which is of tremendous value for natural sciences. The proper choice of the active space is crucial, but a non-trivial task. In this article, we present the neural network (NN) based approach for automatic selection of active spaces, focused on transition metal systems. The training set has been formed from artificial systems composed from one transition metal and various ligands, on which we have performed DMRG and calculated single-site entropy. On the selected set of systems, ranging from small benchmark molecules up to larger challenging systems involving two metallic centers, we demonstrate that our ML models could correctly predict the importance of orbitals with the high accuracy. Also, the ML models show a high degree of transferability on systems much larger than any complex used in training procedures.
△ Less
Submitted 30 November, 2020;
originally announced November 2020.
-
Ground State of the Fe(II)-porphyrin Model System Corresponds to the Quintet State: A DFT and DMRG-based Tailored CC Study
Authors:
Andrej Antalík,
Dana Nachtigallová,
Rabindranath Lo,
Mikuláš Matoušek,
Jakub Lang,
Örs Legeza,
Jiří Pittner,
Pavel Hobza,
Libor Veis
Abstract:
Fe(II)-porphyrins (FeP) play an important role in many reactions relevant to material science and biological processes, due to their closely lying spin states. However, this small energetic separation also makes it challenging to establish the correct spin state ordering. Although the prevalent opinion is that these systems posses the triplet ground state, the recent experiment on Fe(II)-phthalocy…
▽ More
Fe(II)-porphyrins (FeP) play an important role in many reactions relevant to material science and biological processes, due to their closely lying spin states. However, this small energetic separation also makes it challenging to establish the correct spin state ordering. Although the prevalent opinion is that these systems posses the triplet ground state, the recent experiment on Fe(II)-phthalocyanine under conditions matching those of an isolated molecule points toward the quintet ground state. We present a thorough study of FeP model by means of the density functional theory and density matrix renormalization group based tailored coupled clusters, in which we address all previously discussed correlation effects. We examine the importance of geometrical parameters, the Fe-N distances in particular, and conclude that the system possesses the quintet ground state, which is in our calculations well-separated from the triplet state.
△ Less
Submitted 12 June, 2020; v1 submitted 14 January, 2020;
originally announced January 2020.
-
Near-linear Scaling in DMRG-based Tailored Coupled Clusters: An Implementation of DLPNO-TCCSD and DLPNO-TCCSD(T)
Authors:
Jakub Lang,
Andrej Antalík,
Libor Veis,
Jan Brandejs,
Jiří Brabec,
Örs Legeza,
Jiří Pittner
Abstract:
We present a new implementation of DMRG-based tailored coupled clusters method (TCCSD), which employs the domain-based local pair natural orbital approach (DLPNO-TCCSD). Compared to the previous LPNO version of the method, the new implementation is more accurate, offers more favorable scaling and provides more consistent behavior across the variety of systems. On top of the singles and doubles, we…
▽ More
We present a new implementation of DMRG-based tailored coupled clusters method (TCCSD), which employs the domain-based local pair natural orbital approach (DLPNO-TCCSD). Compared to the previous LPNO version of the method, the new implementation is more accurate, offers more favorable scaling and provides more consistent behavior across the variety of systems. On top of the singles and doubles, we include the perturbative triples correction (T), which is able to retrieve even more dynamic correlation. The methods were tested on three systems: tetramethyleneethane, oxo-Mn(Salen) and Iron(II)-porphyrin model. The first two were revisited to assess the performance with respect to LPNO-TCCSD. For oxo-Mn(Salen), we retrieved between 99.8-99.9% of the total canonical correlation energy which is the improvement of 0.2% over the LPNO version in less than 63% of the total LPNO runtime. Similar results were obtained for Iron(II)-porphyrin. When the perturbative triples correction was employed, irrespective of the active space size or system, the obtained energy differences between two spin states were within the chemical accuracy of 1 kcal/mol using the default DLPNO settings.
△ Less
Submitted 14 April, 2020; v1 submitted 29 July, 2019;
originally announced July 2019.
-
Towards the Efficient Local Tailored Coupled Cluster Approximation and the Peculiar Case of Oxo-Mn(Salen)
Authors:
Andrej Antalík,
Libor Veis,
Jiří Brabec,
Örs Legeza,
Jiří Pittner
Abstract:
We introduce a new implementation of the coupled cluster method tailored by matrix product states wave functions (DMRG-TCCSD), which employs the local pair natural orbital approach (LPNO). By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/o…
▽ More
We introduce a new implementation of the coupled cluster method tailored by matrix product states wave functions (DMRG-TCCSD), which employs the local pair natural orbital approach (LPNO). By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/or with larger basis sets. We assessed the accuracy of the approximation using two systems: tetramethyleneethane (TME) and oxo-Mn(Salen). Using the default cut-off parameters, we were able to recover over 99.7% and 99.8% of canonical correlation energy for the triplet and singlet state of TME respectively. In case of oxo-Mn(Salen), we found out that the amount of retrieved canonical correlation energy depends on the size of the active space (CAS) - we retrieved over 99.6% for the larger 27 orbital CAS and over 99.8% for the smaller 22 orbital CAS. The use of LPNO-TCCSD allowed us to perform these calculations up to quadruple-$ζ$ basis set amounting to 1178 basis functions. Moreover, we examined dependance of the ground state of oxo-Mn(Salen) on CAS composition. We found out that the inclusion of 4d$_{xy}$ orbital plays an important role in stabilizing the singlet state at the DMRG-CASSCF level via double-shell effect. However, by including dynamic correlation the ground state was found to be triplet regardless of the size of the basis set or composition of CAS, which is in agreement with previous findings by canonical DMRG-TCCSD in smaller basis.
△ Less
Submitted 7 May, 2019;
originally announced May 2019.
-
Numerical and Theoretical Aspects of the DMRG-TCC Method Exemplified by the Nitrogen Dimer
Authors:
Fabian M. Faulstich,
Mihály Máté,
Andre Laestadius,
Mihály András Csirik,
Libor Veis,
Andrej Antalik,
Jiří Brabec,
Reinhold Schneider,
Jiří Pittner,
Simen Kvaal,
Örs Legeza
Abstract:
In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate chemical properties of the used method, such as energy size extensivity and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use…
▽ More
In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate chemical properties of the used method, such as energy size extensivity and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use of a so-called CAS-ext gap describing the basis splitting between the complete active space and the external part. Moreover, the behavior of the energy error as a function of the optimal basis splitting is discussed. We show numerical investigations on the robustness with respect to the bond dimensions of the single orbital entropy and the mutual information, which are quantities that are used to choose the complete active space. Furthermore, we extend the mathematical analysis with a numerical study on the complete active space dependence of the error.
△ Less
Submitted 20 September, 2018;
originally announced September 2018.
-
Full configuration interaction quantum Monte Carlo benchmark and multireference coupled cluster studies of tetramethyleneethane diradical
Authors:
Libor Veis,
Andrej Antalík,
Örs Legeza,
Ali Alavi,
Jiří Pittner
Abstract:
We have performed a FCI-quality benchmark calculation for the tetramethyleneethane molecule in cc-pVTZ basis set employing a subset of CASPT2(6,6) natural orbitals for the FCIQMC calculation. The results are in an excellent agreement with the previous large scale diffusion Monte Carlo calculations by Pozun et al. and available experimental results. Our computations verified that there is a maximum…
▽ More
We have performed a FCI-quality benchmark calculation for the tetramethyleneethane molecule in cc-pVTZ basis set employing a subset of CASPT2(6,6) natural orbitals for the FCIQMC calculation. The results are in an excellent agreement with the previous large scale diffusion Monte Carlo calculations by Pozun et al. and available experimental results. Our computations verified that there is a maximum on PES of the ground singlet state ($^1\text{A}$) $45^{\circ}$ torsional angle and the corresponding vertical singlet-triplet energy gap is $0.01$ eV. We have employed this benchmark for assessment of the accuracy of MkCCSDT and DMRG-tailored CCSD (TCCSD) methods. MR MkCCSDT with CAS(2,2) model space, though giving good values for the singlet-triplet energy gap, is not able to properly describe the shape of the multireference singlet PES. Similarly, DMRG(24,25) is not able to correctly capture the shape of the singlet surface, due to the missing dynamic correlation. On the other hand, the DMRG-tailored CCSD method describes the shape of the ground singlet state with an excellent accuracy, but for the correct ordering requires computation of the zero-spin-projection component of the triplet state ($^3\text{B}_1$).
△ Less
Submitted 3 January, 2018;
originally announced January 2018.
-
Coupled cluster method with single and double excitations tailored by matrix product state wave functions
Authors:
Libor Veis,
Andrej Antalík,
Jiří Brabec,
Frank Neese,
Örs Legeza,
Jiří Pittner
Abstract:
In the last decade, the quantum chemical version of the density matrix renormalization group (DMRG) method has established itself as the method of choice for calculations of strongly correlated molecular systems. Despite its favourable scaling, it is in practice not suitable for computations of dynamic correlation. We present a novel method for accurate "post-DMRG" treatment of dynamic correlation…
▽ More
In the last decade, the quantum chemical version of the density matrix renormalization group (DMRG) method has established itself as the method of choice for calculations of strongly correlated molecular systems. Despite its favourable scaling, it is in practice not suitable for computations of dynamic correlation. We present a novel method for accurate "post-DMRG" treatment of dynamic correlation based on the tailored coupled cluster (CC) theory in which the DMRG method is responsible for the proper description of non-dynamic correlation, whereas dynamic correlation is incorporated through the framework of the CC theory. We illustrate the potential of this method on prominent multireference systems, in particular N$_2$, Cr$_2$ molecules and also oxo-Mn(Salen) for which we have performed the first "post-DMRG" computations in order to shed light on the energy ordering of the lowest spin states.
△ Less
Submitted 14 November, 2016; v1 submitted 20 June, 2016;
originally announced June 2016.