-
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 in which loosely coupled external programs describe individual subsystems at different resolutions and levels of theory. To make it highly efficient and flexible, we adopt an interoperable approach based on a multiple-program multiple-data (MPMD) paradigm, serving as an intermediary responsible for fast data exchange and interactions betwe…
▽ More
MiMiC is a framework for performing multiscale simulations in which loosely coupled external programs describe individual subsystems at different resolutions and levels of theory. To make it highly efficient and flexible, we adopt an interoperable approach based on a multiple-program multiple-data (MPMD) 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 26 September, 2024; v1 submitted 27 March, 2024;
originally announced March 2024.
-
Multiscale Biomolecular Simulations in the Exascale Era
Authors:
David Carrasco-Busturia,
Emiliano Ippoliti,
Simone Meloni,
Ursula Rothlisberger,
Jógvan Magnus Haugaard Olsen
Abstract:
The complexity of biological systems and processes, spanning molecular to macroscopic scales, necessitates the use of multiscale simulations to get a comprehensive understanding. Quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations are crucial for capturing processes beyond the reach of classical MD simulations. The advent of exascale computing offers unprecedented opp…
▽ More
The complexity of biological systems and processes, spanning molecular to macroscopic scales, necessitates the use of multiscale simulations to get a comprehensive understanding. Quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations are crucial for capturing processes beyond the reach of classical MD simulations. The advent of exascale computing offers unprecedented opportunities for scientific exploration, not least within life sciences, where simulations are essential to unravel intricate molecular mechanisms underlying biological processes. However, leveraging the immense computational power of exascale computing requires innovative algorithms and software designs. In this context, we discuss the current status and future prospects of multiscale biomolecular simulations on exascale supercomputers with a focus on QM/MM MD. We highlight our own efforts in developing a versatile and high-performance multiscale simulation framework with the aim of efficient utilization of state-of-the-art supercomputers. We showcase its application in uncovering complex biological mechanisms and its potential for leveraging exascale computing.
△ Less
Submitted 3 April, 2024; v1 submitted 3 March, 2024;
originally announced March 2024.
-
Wavefunction-based electrostatic-embedding QM/MM using CFOUR through MiMiC
Authors:
Till Kirsch,
Jógvan Magnus Haugaard Olsen,
Viacheslav Bolnykh,
Simone Meloni,
Emiliano Ippoliti,
Ursula Rothlisberger,
Michele Cascella,
Jürgen Gauss
Abstract:
We present an interface of the wavefunction-based quantum-chemical software CFOUR to the multiscale modeling framework MiMiC. Electrostatic embedding of the quantummechanical (QM) part is achieved by analytic evaluation of one-electron integrals in CFOUR, while the rest of the QM/MM operations are treated according to the previous MiMiC-based QM/MM implementation. Long-range electrostatic interact…
▽ More
We present an interface of the wavefunction-based quantum-chemical software CFOUR to the multiscale modeling framework MiMiC. Electrostatic embedding of the quantummechanical (QM) part is achieved by analytic evaluation of one-electron integrals in CFOUR, while the rest of the QM/MM operations are treated according to the previous MiMiC-based QM/MM implementation. Long-range electrostatic interactions are treated by a multipole expansion of the potential from the QM electron density to reduce the computational cost without loss of accuracy. Testing on model water/water systems, we verified that the CFOUR interface to MiMiC is robust, guaranteeing fast convergence of the SCF cycles and optimal conservation of the energy during the integration of the equations of motion. Finally, we verified that the CFOUR interface to MiMiC is compatible with the use of a QM/QM multiple time-step algorithm, which effectively reduces the cost of AIMD or QM/MM-MD simulations using higher level wavefunction-based approaches compared to cheaper density-functional theory based ones. The new wavefunction-based AIMD and QM/MM-MD implementation was tested and validated for a large number of wavefunction approaches, including Hartree-Fock and post-Hartree-Fock methods like Moller-Plesset, coupled cluster, and complete active space self-consistent field.
△ Less
Submitted 24 August, 2021;
originally announced August 2021.
-
The DIRAC code for relativistic molecular calculations
Authors:
Trond Saue,
Radovan Bast,
Andre Severo Pereira Gomes,
Hans Jørgen Aagaard Jensen,
Lucas Visscher,
Ignacio Agustın Aucar,
Roberto Di Remigio,
Kenneth G. Dyall,
Ephraim Eliav,
Elke Faßhauer,
Timo Fleig,
Loıc Halbert,
Erik Donovan Hedegård,
Benjamin Helmich-Paris,
Miroslav Iliaš,
Christoph R. Jacob,
Stefan Knecht,
Jon K Laerdahl,
Marta L. Vidal,
Malaya K Nayak,
Małgorzata Olejniczak,
Jógvan Magnus Haugaard Olsen,
Markus Pernpointner,
Bruno Senjean,
Avijit Shee
, et al. (2 additional authors not shown)
Abstract:
DIRAC is a freely distributed general-purpose program system for 1-, 2- and 4-component relativistic molecular calculations at the level of Hartree--Fock, Kohn--Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, coupled cluster and electron propagator theory. At the self-consistent-field level a highly original scheme, bas…
▽ More
DIRAC is a freely distributed general-purpose program system for 1-, 2- and 4-component relativistic molecular calculations at the level of Hartree--Fock, Kohn--Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, coupled cluster and electron propagator 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 the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding, and frozen density embedding models. DIRAC was one of the earliest codes for relativistic molecular calculations and remains a reference in its field.
△ Less
Submitted 14 February, 2020;
originally announced February 2020.
-
Response properties of embedded molecules through the polarizable embedding model
Authors:
Casper Steinmann,
Peter Reinholdt,
Morten Steen Nørby,
Jacob Kongsted,
Jógvan Magnus Haugaard Olsen
Abstract:
The polarizable embedding (PE) model is a fragment-based quantum-classical approach aimed at accurate inclusion of environment effects in quantum-mechanical response property calculations. The aim of this tutorial is to give insight into the practical use of the PE model. Starting from a set of molecular structures and until you arrive at the final property, there are many crucial details to consi…
▽ More
The polarizable embedding (PE) model is a fragment-based quantum-classical approach aimed at accurate inclusion of environment effects in quantum-mechanical response property calculations. The aim of this tutorial is to give insight into the practical use of the PE model. Starting from a set of molecular structures and until you arrive at the final property, there are many crucial details to consider in order to obtain trustworthy results in an efficient manner. To lower the threshold for new users wanting to explore the use of the PE model, we describe and discuss important aspects related to its practical use. This includes directions on how to generate input files and how to run a calculation.
△ Less
Submitted 10 April, 2018;
originally announced April 2018.
-
Modeling the absorption spectrum of the permanganate ion in vacuum and in aqueous solution
Authors:
Jógvan Magnus Haugaard Olsen,
Erik Donovan Hedegård
Abstract:
The absorption spectrum of the MnO$_{4}$$^{-}$ ion has been a test-bed for quantum-chemical methods over the last decades. Its correct description requires highly-correlated multiconfigurational methods, which are incompatible with the inclusion of finite-temperature and solvent effects due to their high computational demands. Therefore, implicit solvent models are usually employed. Here we show t…
▽ More
The absorption spectrum of the MnO$_{4}$$^{-}$ ion has been a test-bed for quantum-chemical methods over the last decades. Its correct description requires highly-correlated multiconfigurational methods, which are incompatible with the inclusion of finite-temperature and solvent effects due to their high computational demands. Therefore, implicit solvent models are usually employed. Here we show that implicit solvent models are not sufficiently accurate to model the solvent shift of MnO$_{4}$$^{-}$, and we analyze the origins of their failure. We obtain the correct solvent shift for MnO$_{4}$$^{-}$ in aqueous solution by employing the polarizable embedding (PE) model combined with a range-separated complete active space short-range density functional theory method (CAS-srDFT). Finite-temperature effects are taken into account by averaging over structures obtained from ab initio molecular dynamics simulations. The explicit treatment of finite-temperature and solvent effects facilitates the interpretation of the bands in the low-energy region of the MnO$_{4}$$^{-}$ absorption spectrum, whose assignment has been elusive.
△ Less
Submitted 22 February, 2017; v1 submitted 20 January, 2017;
originally announced January 2017.