Related Experiment Video
Updated: May 22, 2026

05:57
Synthesizing Amino Acids Modified with Reactive Carbonyls in Silico to Assess Structural Effects Using Molecular Dynamics Simulations
Published on: April 26, 2024
Comment on "A minimal implementation of the AMBER-GAUSSIAN interface for ab initio QM/MM-MD simulation"
Benjamin P Roberts1, Gustavo M Seabra, Adrian E Roitberg
1Quantum Theory Project, Department of Physics, University of Florida, Gainesville, Florida 32611-8435, USA.
Journal of Computational Chemistry
|May 10, 2012
Summary
The Program for User Package Interfacing and Linking (PUPIL) system
Area of Science:
- Computational Chemistry
- Molecular Modeling
- Quantum Mechanics
Background:
- The Program for User Package Interfacing and Linking (PUPIL) is used for multiscale simulations.
- A recent critique questioned the use of PUPIL for linking AMBER and GAUSSIAN.
- Accurate force calculations are crucial in QM/MM simulations.
Purpose of the Study:
- To address a recent critique of the PUPIL system.
- To clarify the implementation of specific computational methods within PUPIL.
- To discuss the broader context of open-source software in computational chemistry.
Main Methods:
- Commenting on a published critique.
- Reviewing the implementation details of the PUPIL system.
- Discussing the calculation of forces in QM/MM simulations.
Main Results:
- The force calculation method criticized was implemented in PUPIL version 1.3 (released December 2009).
- Specific characterizations in the critique are considered doubtful.
- The discussion is framed within the context of open-source computational chemistry tools.
Conclusions:
- The PUPIL system's capabilities, particularly regarding QM/MM force calculations, are clarified.
- The critique may not reflect the most current or complete understanding of PUPIL's features.
- Highlights the importance of accurate awareness of open-source software functionalities.
Related Concept Videos
The Quantum-Mechanical Model of an Atom
Shortly after de Broglie published his ideas that the electron in a hydrogen atom could be better thought of as being a circular standing wave instead of a particle moving in quantized circular orbits, Erwin Schrödinger extended de Broglie’s work by deriving what is now known as the Schrödinger equation. When Schrödinger applied his equation to hydrogen-like atoms, he was able to reproduce Bohr’s expression for the energy and, thus, the Rydberg formula governing hydrogen spectra. Schrödinger...
Hybridization of Atomic Orbitals II
sp3d and sp3d 2 Hybridization
Atomic Absorption Spectroscopy: Atomization Methods
Atomic Absorption Spectroscopy (AAS) atomizes samples through flame atomization or electrothermal atomization. Flame atomization typically involves a nebulizer and spray chamber assembly to combine the sample with a fuel–oxidant mixture, creating a fine aerosol mist that enters a burner. Typically, the fuel and oxidant are combined in an approximately stoichiometric ratio. However, for atoms that are easily oxidized, a fuel-rich mixture may be more advantageous. Only about 5% of the aerosol...
Hybridization of Atomic Orbitals I
The mathematical expression known as the wave function, ψ, contains information about each orbital and the wavelike properties of electrons in an isolated atom. When atoms are bound together in a molecule, the wave functions combine to produce new mathematical descriptions that have different shapes. This process of combining the wave functions for atomic orbitals is called hybridization and is mathematically accomplished by the linear combination of atomic orbitals. The new orbitals that...
Molecular Orbital Theory II
Molecular Orbital Energy Diagrams
NMR Spectroscopy Of Amines
In proton NMR spectroscopy, primary amines and secondary amines showcase their N–H protons as a broad signal in the chemical shift range between δ 0.5 and 5 ppm. The exact position in this range depends on several factors, including sample concentration, hydrogen bonding, and the type of solvent used. Since amine protons undergo fast proton exchange in solution, the protons are labile and therefore do not participate in any splitting with adjacent protons. Thus, the observed peak is broad and...

