Systematically Constructing QM regions for Calculating QM/MM Absorption Spectra
5 affiliations
- 1Department of Chemistry and The PULSE Institute
- 2Stanford University
- 3SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025
- 4Department of Chemistry
- 5University of North Carolina at Chapel Hill
Abstract
Hybrid simulations utilizing both quantum mechanics and molecular mechanics (QM/MM) are widely used to model complex problems, such as the electronic excitation of chromophores embedded in a protein environment. Its success requires a careful choice of the portion of the system to include in the QM region to obtain accurate results. Chemical intuition and/or proximity arguments are often used to select the QM region for excited state calculations aiming to compute absorption spectra, but such choices do not guarantee accuracy. We compare four schemes for selecting the minimal set of protein residues that must be treated quantum mechanically to reproduce a reference spectrum: distance analysis, charge shift analysis, point charge variation analysis, and subtractive analysis. Charge shift analysis is the most robust: in both systems studied it identifies the smallest QM region that reproduces the reference spectral shape and absorption maximum, at moderate computational cost, using a threshold that transfers between systems. Because the resulting regions are both accurate and minimal, converged excited state calculations become affordable at scale, laying a foundation for large, internally consistent datasets of excited state properties in extended systems. We release an open-source implementation that automates the full workflow, including residue selection and determination of the QM/MM boundary, so a converged QM region can be built for a new system without manual intervention.
Related
- Optical absorption spectra from the GW-Bethe-Salpeter framework for thousands of atoms
- Middle-range-energy structures in Ti L 2 , 3 − edge x-ray absorption spectra of titanium oxides
- Probing multi-state dark matter via optical absorption lines in DESI spectra
- Absorption and scattering spectra of massive scalar waves in charged regular black hole spacetimes
- From GROMACS Molecular Dynamics Simulations to Electronic Absorption Spectra: A Tutorial for Small Molecules in Organic Solvents.
BibTeX
@misc{zhu2026systematically,
title = {{Systematically Constructing QM regions for Calculating QM/MM Absorption Spectra}},
author = {Mingning Zhu and Vinícius Wilian D. Cruzeiro and Elisa Pieri and Todd J. Martínez},
howpublished = {ChemRxiv},
year = {2026},
month = oct,
doi = {10.26434/chemrxiv.15010595/v1},
url = {https://chemrxiv.org/doi/full/10.26434/chemrxiv.15010595/v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Systematically Constructing QM regions for Calculating QM/MM Absorption Spectra
AU - Mingning Zhu
AU - Vinícius Wilian D. Cruzeiro
AU - Elisa Pieri
AU - Todd J. Martínez
PY - 2026
DA - 2026/10/16
PB - ChemRxiv
DO - 10.26434/chemrxiv.15010595/v1
UR - https://chemrxiv.org/doi/full/10.26434/chemrxiv.15010595/v1
AB - Hybrid simulations utilizing both quantum mechanics and molecular mechanics (QM/MM) are widely used to model complex problems, such as the electronic excitation of chromophores embedded in a protein environment. Its success requires a careful choice of the portion of the system to include in the QM region to obtain accurate results. Chemical intuition and/or proximity arguments are often used to select the QM region for excited state calculations aiming to compute absorption spectra, but such choices do not guarantee accuracy. We compare four schemes for selecting the minimal set of protein residues that must be treated quantum mechanically to reproduce a reference spectrum: distance analysis, charge shift analysis, point charge variation analysis, and subtractive analysis. Charge shift analysis is the most robust: in both systems studied it identifies the smallest QM region that reproduces the reference spectral shape and absorption maximum, at moderate computational cost, using a threshold that transfers between systems. Because the resulting regions are both accurate and minimal, converged excited state calculations become affordable at scale, laying a foundation for large, internally consistent datasets of excited state properties in extended systems. We release an open-source implementation that automates the full workflow, including residue selection and determination of the QM/MM boundary, so a converged QM region can be built for a new system without manual intervention.
ER -