Faster simulations of molecular interactions could accelerate drug and material design
Researchers have combined quantum computing techniques with artificial intelligence to simulate how molecules interact at ultrafast timescales—while cutting computational costs by orders of magnitude. The breakthrough could enable faster, cheaper development of new pharmaceuticals, materials, and chemical processes by making detailed molecular simulations practical for industry use.
Originaltitel: Two-dimensional infrared spectroscopy of solute–solvent complexes from linear-scaling DFT and machine learning
Two-dimensional infrared (2DIR) spectroscopy captures vibrational correlations on femtosecond timescales, offering direct insights into hydrogen(H)-bonding dynamics and other ultrafast molecular processes. However, interpreting these spectra requires simulations that accurately describe solute-solvent interactions over realistic timescales and system sizes. While classical approaches using empirical frequency maps are common, ab initio molecular dynamics (AIMD) offers a more rigorous alternative by treating dynamics and vibrational frequencies on a consistent theoretical level. The primary drawback of AIMD is its high computational cost, which typically limits simulations to short trajectories. Here, we introduce a hybrid strategy that combines linear-scaling density functional theory (LS-DFT) with a machine-learned (ML) interatomic potential. We use short LS-DFT simulations to generate reference energies, forces, and electron-density-derived dipole moments, which then serve as training data for a DeepMD model. The resulting ML potential allows nanosecond-scale dynamics at a fraction of the ab initio cost. We demonstrate this approach for N-methylacetamide in methanol, a model system known to form distinct H-bonded subpopulations. A key advantage of our method is that it bypasses the need for empirical frequency maps. Instead, molecular dipoles are learned directly from the electron density, and instantaneous vibrational frequencies are calculated from stable numerical Hessians. The resulting linear and 2DIR spectra show excellent agreement with experiment, accurately reproducing the characteristic doublet structure of the amide I band. This framework provides a practical and accurate route to simulating vibrational spectra at the AIMD level of theory for a wide range of IR-active solutes and H-bonded complexes.