Related Experiment Video
Updated: Aug 13, 2026

Isotopic Effect in Double Proton Transfer Process of Porphycene Investigated by Enhanced QM/MM Method
Published on: July 19, 2019
Lennard-Jones parameters for the combined QM/MM method using the B3LYP/6-31G*/AMBER potential
Marek Freindorf1, Yihan Shao, Thomas R Furlani
1Center for Computational Research, The State University of New York at Buffalo, Buffalo, NY 14260, USA. mfrein@buffalo.edu
Abstract:
A combined DFT quantum mechanical and AMBER molecular mechanical potential (QM/MM) is presented for use in molecular modeling and molecular simulations of large biological systems. In our approach we evaluate Lennard-Jones parameters describing the interaction between the quantum mechanical (QM) part of a system, which is described at the B3LYP/6-31+G* level of theory, and the molecular mechanical (MM) part of the system, described by the AMBER force field. The Lennard-Jones parameters for this potential are obtained by calculating hydrogen bond energies and hydrogen bond geometries for a large set of bimolecular systems, in which one hydrogen bond monomer is described quantum mechanically and the other is treated molecular mechanically. We have investigated more than 100 different bimolecular systems, finding very good agreement between hydrogen bond energies and geometries obtained from the combined QM/MM calculations and results obtained at the QM level of theory, especially with respect to geometry. Therefore, based on the Lennard-Jones parameters obtained in our study, we anticipate that the B3LYP/6-31+G*/AMBER potential will be a precise tool to explore intermolecular interactions inside a protein environment.
Related Concept Videos
Electronic Structure of Atoms
An atom comprises protons and neutrons, which are contained inside the dense, central core called the nucleus, with electrons present around the nucleus. Taking into account the wave–particle duality of electrons and the uncertainty in position around the nucleus, quantum mechanics provides a more accurate model for the atomic structure. It describes atomic orbitals as the regions around the nucleus where electrons of discrete energy exist, characterized by four quantum numbers: n, l, ml, and...
Molecular Geometry and Dipole Moments
¹H NMR: Complex Splitting
Splitting diagrams or splitting tree diagrams are routinely used to depict such complex couplings. While drawing splitting diagrams, the splitting with the larger coupling constant is usually applied first.
π Molecular Orbitals of 1,3-Butadiene
The simplest conjugated diene is 1,3-butadiene: a four-carbon system where each carbon is sp2-hybridized and has an unhybridized p orbital that contains an unpaired electron. According to molecular orbital theory, atomic orbitals combine to form molecular orbitals such that the number...
Spin–Spin Coupling Constant: Overview
Qualitatively, any spin plus-half nucleus polarizes the spins of its electrons to the minus-half state. Consequently, the paired electron in the hydrogen–carbon bond must have a...
Van der Waals Equation
First, the attractive forces between molecules, which are stronger at higher densities and reduce the pressure, are considered by adding to the pressure a term equal to the square of the molar density multiplied by a positive coefficient a. Second, the volume...

