sbml4md: A computational platform for System-Bath Modeling via Molecular Dynamics powered by Machine Learning

Abstract

We introduce sbml4md, a newly developed algorithm implemented as a software package to extract parameters of multimode anharmonic Brownian (MAB) models from molecular dynamics (MD) trajectories for simulating nonlinear vibrational spectra of intramolecular modes of molecular liquids. By leveraging machine learning (ML) techniques to capture vibrational anharmonicity, intermolecular couplings, and bath correlation functions for each mode, sbml4md obviates empirical fitting and enables the modeling of environments with spatial and temporal heterogeneity. This work provides a set of parameters specifically tailored for the Hierarchical Equations of Motion (HEOM) framework, enabling numerically "exact" simulations of nonlinear vibrational spectra. Building upon our previous implementation for intramolecular vibrational modes [Park, Jo, and Tanimura, J. Chem. Phys. 163, 214104 (2025)], the present code enhances optimization efficiency by explicitly accounting for intermolecular vibrational contributions. This extension enables sbml4md to broaden the applicability of HEOM-based dynamical modeling by seamlessly integrating classical MD approaches, thereby providing a flexible and scalable framework for simulating both linear and nonlinear spectra under realistic conditions with minimal empirical input. The accompanying ML code, written in Python, is provided as supporting material.

Special Topic, Ultrafast Spectroscopy and Dynamics

References in corpus (11)

• Perspective: Numerically "exact" approach to open quantum dynamics: The hierarchical equations of motion (HEOM)
• Analysis of 2D THz-Raman spectroscopy using a non-Markovian Brownian oscillator model with nonlinear system-bath interactions
• Critical role of quantum dynamical effects in the Raman spectroscopy of liquid water
• Linear, third- and fifth-order nonlinear spectroscopy of a charge transfer system coupled to an underdamped vibration
• Modeling intermolecular and intramolecular modes of liquid water using multiple heat baths: Machine learning approach
• Discretized hierarchical equations of motion in mixed Liouville--Wigner space for two-dimensional vibrational spectroscopies of liquid water
• Full molecular dynamics simulations of liquid water and carbon tetrachloride for two-dimensional Raman spectroscopy in the frequency domain
• Simulating two-dimensional correlation spectroscopies with third-order infrared and fifth-order infrared--Raman processes of liquid water
• Quantum hierarchical Fokker-Planck equations with U(1) gauge fields: Application to the Aharonov-Bohm ring
• A Multimode Classical Hierarchical Fokker-Planck Equations Approach to Molecular Vibrations: Simulating Two-Dimensional Spectra
• System-Bath Modeling in Vibrational Spectroscopy via Molecular Dynamics: A Machine Learning Framework for Hierarchical Equations of Motion (HEOM)