# Next Generation Methods for Advanced Condensed Phase Simulations in Q-Chem

> **NIH NIH R44** · Q-CHEM, INC. · 2020 · $504,923

## Abstract

PROJECT SUMMARY
Next Generation Methods for Advanced Condensed Phase Simulations in Q-Chem
 Biophysical systems exist in the condensed phase, and that is the environment in which their
properties should be computer-modeled. The correct theory to describe the electrons is using ab
initio (AI) quantum mechanics (QM), whilst nuclear motion requires molecular dynamics (MD). The
combination, AIMD, is thus the appropriate tool for biophysical simulations. While use of AIMD
is vastly more expensive than MD with empirical potentials, it is nonetheless the standard to aspire
to. AIMD enables correct treatment of bond-breaking for reactive processes, as well an accurate
description of the non-bonded interactions that determine solvation and conformational preferences.
 This Phase II proposal has the objective of bringing a production level AIMD code to the Q-Chem
software package. The key justiﬁcation for the proposed work, and the potential value of the resulting
product is that it will bring together capabilities that are not found jointly in any other AIMD code.
The valuable synergy between the density functional theory implementation for periodic boundary
conditions (DFT-PBC), and advanced algorithms for efﬁciently and accurate propagating the MD is
the core innovation of this project.
 With regard to DFT-PBC (the ﬁrst speciﬁc aim), the focus is on implementing high precision, high
efﬁciency algorithms for the critical components of DFT with advanced functionals. Our code will
support the latest meta-generalized gradient approximations (mGGAs), with inclusion of non-local
van der Waals density functionals, that are not available in DFT-PBC codes to date. We will addition-
ally provide support for range-separated exact exchange, with high efﬁciency. These capabilities will
come with energies and gradients. Our software framework can also permit all-electron calculations
as needed e.g. for NMR properties that depend on the electron density at the nucleus. Our modular
code will support efﬁcient on-node parallelism.
 To propagate MD efﬁciently and stably (the second speciﬁc aim), we employ two innovative statis-
tical mechanics (SM) algorithms that have been proven in conventional MD, but are not yet available
in any production AIMD code. First, we are extending the inertial extended Lagrangian self-consistent
ﬁeld (iEL/SCF) method to work robustly and efﬁciently with AIMD, building upon promising Phase
I results, by combining it with a stochastic-isokinetic integration (SII) scheme to enable a single but
larger MD time step. Second, we will explore the combination of iEL/SCF-SII with a multiple time-
stepping method in which will explore whether different components of the QM force can be updated
on different timescales in the AIMD.
 In ﬁnal Aim 3 we test the combined DFT-PBC and iEL/SCF-SII capabilities on biophysical appli-
cations including zwitterionic glycine and valine peptides in aqueous solution and molecular crystals.
1

## Key facts

- **NIH application ID:** 10011528
- **Project number:** 2R44GM128480-02A1
- **Recipient organization:** Q-CHEM, INC.
- **Principal Investigator:** Evgeny Epifanovsky
- **Activity code:** R44 (R01, R21, SBIR, etc.)
- **Funding institute:** NIH
- **Fiscal year:** 2020
- **Award amount:** $504,923
- **Award type:** 2
- **Project period:** 2018-09-05 → 2022-03-31

## Primary source

NIH RePORTER: https://reporter.nih.gov/project-details/10011528

## Citation

> US National Institutes of Health, RePORTER application 10011528, Next Generation Methods for Advanced Condensed Phase Simulations in Q-Chem (2R44GM128480-02A1). Retrieved via AI Analytics 2026-05-22 from https://api.ai-analytics.org/grant/nih/10011528. Licensed CC0.

---

*[NIH grants dataset](/datasets/nih-grants) · CC0 1.0*
