Consistent and Generalizable Effective Model Hamiltonian Framework for Studying Nonadiabatic Dynamics in the Condensed Phase

Simulating nonadiabatic dynamics in complex, condensed-phase systems presents a formidable computational challenge, demanding the development of effective model Hamiltonians that capture the essential physics of electronic and nuclear interactions. This Review charts the evolution of such models, from the foundational two-state spin-boson model to multistate Frenkel exciton models, highlighting the limitations of traditional approaches, particularly the isolated bath assumption, which neglects crucial environmental correlations. We focus on the recently developed multistate harmonic (MSH) model, a general and consistent framework for mapping information from all-atom simulations onto an effective Hamiltonian. The MSH model overcomes the shortcomings of previous models by systematically satisfying all pairwise reorganization energy constraints for a multistate system. This is achieved through a novel extension of the nuclear coordinate space, which provides a physically grounded and geometrically intuitive representation of a globally shared, correlated bath. We detail the construction of the MSH Hamiltonian and its equivalent multistate reaction coordinate (MRC) representation and discuss its applications in conjunction with various dynamical methods, including perturbative quantum master equations and semiclassical nonadiabatic dynamics approaches, as well as rate constant and time-dependent rate. The MSH/MRC models not only provide a robust platform for predictive simulations of charge and energy transfer in the condensed phase but also serve as an invaluable tool for benchmarking the accuracy of approximate quantum dynamics methods.