Force Field Development with Density-Based Energy Decomposition Analysis

2019 ◽  
Vol 15 (9) ◽  
pp. 5001-5013 ◽  
Author(s):  
Akshaya K. Das ◽  
Lars Urban ◽  
Itai Leven ◽  
Matthias Loipersberger ◽  
Abdulrahman Aldossary ◽  
...  

Author(s):  
Yudong Qiu ◽  
Daniel Smith ◽  
Chaya Stern ◽  
mudong feng ◽  
Lee-Ping Wang

<div>The parameterization of torsional / dihedral angle potential energy terms is a crucial part of developing molecular mechanics force fields.</div><div>Quantum mechanical (QM) methods are often used to provide samples of the potential energy surface (PES) for fitting the empirical parameters in these force field terms.</div><div>To ensure that the sampled molecular configurations are thermodynamically feasible, constrained QM geometry optimizations are typically carried out, which relax the orthogonal degrees of freedom while fixing the target torsion angle(s) on a grid of values.</div><div>However, the quality of results and computational cost are affected by various factors on a non-trivial PES, such as dependence on the chosen scan direction and the lack of efficient approaches to integrate results started from multiple initial guesses.</div><div>In this paper we propose a systematic and versatile workflow called \textit{TorsionDrive} to generate energy-minimized structures on a grid of torsion constraints by means of a recursive wavefront propagation algorithm, which resolves the deficiencies of conventional scanning approaches and generates higher quality QM data for force field development.</div><div>The capabilities of our method are presented for multi-dimensional scans and multiple initial guess structures, and an integration with the MolSSI QCArchive distributed computing ecosystem is described.</div><div>The method is implemented in an open-source software package that is compatible with many QM software packages and energy minimization codes.</div>


2015 ◽  
Vol 11 ◽  
pp. 2727-2736 ◽  
Author(s):  
Diego M Andrada ◽  
Nicole Holzmann ◽  
Thomas Hamadi ◽  
Gernot Frenking

Fifteen cyclic and acylic carbenes have been calculated with density functional theory at the BP86/def2-TZVPP level. The strength of the internal X→p(π) π-donation of heteroatoms and carbon which are bonded to the C(II) atom is estimated with the help of NBO calculations and with an energy decomposition analysis. The investigated molecules include N-heterocyclic carbenes (NHCs), the cyclic alkyl(amino)carbene (cAAC), mesoionic carbenes and ylide-stabilized carbenes. The bonding analysis suggests that the carbene centre in cAAC and in diamidocarbene have the weakest X→p(π) π-donation while mesoionic carbenes possess the strongest π-donation.


2021 ◽  
Vol 23 (36) ◽  
pp. 20533-20540
Author(s):  
Gustavo Cárdenas ◽  
Álvaro Pérez-Barcia ◽  
Marcos Mandado ◽  
Juan J. Nogueira

The interactions that control the permeation of cisplatin through a DOPC bilayer are unveiled by a QM/MM EDA scheme.


Sign in / Sign up

Export Citation Format

Share Document