Recent Advances in First-Principles Based Molecular Dynamics.

First-principles molecular dynamics (FPMD) and its quantum mechanical-molecular mechanical (QM/MM) extensions are powerful tools to follow the real-time dynamics of a broad variety of systems in their ground as well as electronically excited states. The continued advances in computational power have enabled simulations of QM regions of larger sizes for more extended time scales. In addition, development of the parallel algorithms has boosted the performance of QM/MM methods even on existing computer architectures. In the case of density functional-based FPMD, systems of several hundreds to thousands of atoms can now be customarily simulated for tens to hundreds of picoseconds. In spite of this progress, the time scale limitations remain severe, especially when high-rung exchange-correlation functionals or high-level wave function based quantum mechanical methods are used. To ameliorate this, a large number of enhanced sampling methods have been introduced but most of the approaches that have been developed to increase the efficiency of FPMD based simulations sacrifice the real-time dynamics in favor of enhancing sampling. Here, we present some recent advances in boosting the efficiency of FPMD based simulations while keeping the full dynamic information. These include a highly efficient recent implementation of FPMD-based QM/MM simulations that not only enables fully flexible combinations of different electronic structure methods and force fields via a highly efficient communication library, it also fully exploits parallelism for both quantum and classical descriptions. The second type of acceleration methods we discuss is a large family of specially devised multiple-time-step algorithms that make use of suitable breakups of the total nuclear forces into fast components that can be calculated via lower level methods and slowly varying correction forces evaluated with a high-level method at long time intervals. The computational gain of this scheme mostly depends on the cost difference between the two methods and advantageous combinations can yield large speedups without compromising the accuracy of the high-level method. And finally, the third class of FPMD acceleration methods presented here are machine learning models to accelerated FPMD and their powerful combinations with multiple-time-step techniques. The combination of all the approaches enables substantial speedups of FPMD simulations of several orders of magnitude while fully preserving the real-time dynamics and accuracy.

A universal co-solvent dilution strategy enables facile and cost-effective fabrication of perovskite photovoltaics.

Cost management and toxic waste generation are two key issues that must be addressed before the commercialization of perovskite optoelectronic devices. We report a groundbreaking strategy for eco-friendly and cost-effective fabrication of highly efficient perovskite solar cells. This strategy involves the usage of a high volatility co-solvent, which dilutes perovskite precursors to a lower concentration (<0.5 M) while retaining similar film quality and device performance as a high concentration (>1.4 M) solution. More than 70% of toxic waste and material cost can be reduced. Mechanistic insights reveal ultra-rapid evaporation of the co-solvent together with beneficial alteration of the precursor colloidal chemistry upon dilution with co-solvent, which in-situ studies and theoretical simulations confirm. The co-solvent tuned precursor colloidal properties also contribute to the enhancement of the stability of precursor solution, which extends its processing window thus minimizing the waste. This strategy is universally successful across different perovskite compositions, and scales from small devices to large-scale modules using industrial spin-coating, potentially easing the lab-to-fab translation of perovskite technologies.

Wavefunction-Based Electrostatic-Embedding QM/MM Using CFOUR through MiMiC.

We present an interface of the wavefunction-based quantum chemical software CFOUR to the multiscale modeling framework MiMiC. Electrostatic embedding of the quantum mechanical (QM) part is achieved by analytic evaluation of one-electron integrals in CFOUR, while the rest of the QM/molecular mechanical (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 self-consistent field 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 ab initio MD (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 implementations were tested and validated for a large number of wavefunction approaches, including Hartree-Fock and post-Hartree-Fock methods like Møller-Plesset, coupled-cluster, and complete active space self-consistent field.

Accuracy of Molecular Simulation-Based Predictions of koff Values: A Metadynamics Study.

The koff values of ligands unbinding to proteins are key parameters for drug discovery. Their predictions based on molecular simulation may under- or overestimate experiment in a system- and/or technique-dependent way. Here we use an established method-infrequent metadynamics, based on the AMBER force field-to compute the koff of the ligand iperoxo (in clinical use) targeting the muscarinic receptor M2. The ligand charges are calculated by either (i) the Amber standard procedure or (ii) B3LYP-DFT. The calculations using (i) turn out not to provide a reasonable estimation of the transition-state free energy. Those using (ii) differ from experiment by 2 orders of magnitude. On the basis of B3LYP DFT QM/MM simulations, we suggest that the observed discrepancy in (ii) arises, at least in part, from the lack of electronic polarization and/or charge transfer in biomolecular force fields. These issues might be present in other systems, such as DNA-protein complexes.

Molecular Basis of CLC Antiporter Inhibition by Fluoride.

CLC channels and transporters conduct or transport various kinds of anions, with the exception of fluoride, which acts as an effective inhibitor. Here, we performed sub-nanosecond DFT-based QM/MM simulations of the E. coli anion/proton exchanger ClC-ec1 and observed that fluoride binds incoming protons within the selectivity filter, with excess protons shared with the gating glutamate E148. Depending on E148 conformation, the competition for the proton can involve either a direct F-/E148 interaction or the modulation of water molecules bridging the two anions. The direct interaction locks E148 in a conformation that does not allow for proton transport, and thus inhibits protein function.

Extreme Scalability of DFT-Based QM/MM MD Simulations Using MiMiC.

We present a highly scalable DFT-based QM/MM implementation developed within MiMiC, a recently introduced multiscale modeling framework that uses a loose-coupling strategy in conjunction with a multiple-program multiple-data (MPMD) approach. The computation of electrostatic QM/MM interactions is parallelized exploiting both distributed- and shared-memory strategies. Here, we use the efficient CPMD and GROMACS programs as QM and MM engines, respectively. The scalability is demonstrated through large-scale benchmark simulations of realistic biomolecular systems employing non-hybrid and hybrid GGA exchange-correlation functionals. We show that the loose-coupling strategy adopted in MiMiC, with its inherent high flexibility, does not carry any significant computational overhead compared to a tight-coupling scheme. Furthermore, we demonstrate that the adopted parallelization strategy enables scaling up to 13,000 CPU cores with efficiency above 70%, thus making DFT-based QM/MM MD simulations using hybrid functionals at the nanosecond scale accessible.

MiMiC: A Novel Framework for Multiscale Modeling in Computational Chemistry.

We present a flexible and efficient framework for multiscale modeling in computational chemistry (MiMiC). It is based on a multiple-program multiple-data (MPMD) model with loosely coupled programs. Fast data exchange between programs is achieved through the use of MPI intercommunicators. This allows exploiting the existing parallelization strategies used by the coupled programs while maintaining a high degree of flexibility. MiMiC has been used in a new electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) implementation coupling the highly efficient CPMD and GROMACS programs, but it can also be extended to use other programs. The framework can also be utilized to extend the partitioning of the system into several domains that can be treated using different models, such as models based on wave function or density functional theory as well as coarse-graining and continuum models. The new QM/MM implementation treats long-range electrostatic QM-MM interactions through the multipoles of the QM subsystem which substantially reduces the computational cost without loss of accuracy compared to an exact treatment. This enables QM/MM molecular dynamics (MD) simulations of very large systems.