Fragmentation kinetics of a Morse oscillator chain under tension

The bond dissociation kinetics of tethered atomic (Morse potential) chains under tensile stress is studied. Both RRKM (fully anharmonic, Monte Carlo) and RRK (harmonic appproximation) theory are applied to predict bond dissociation rate constants as a function of energy and tensile force. For chains with N ⩾ 3 atoms a hybrid statistical theory is used involving a harmonic approximation for motion in the transition state for bond dissociation. For chains with N = 2-5 atoms, while the RRK approximation significantly overestimates the dissociation rate constant, the fully anharmonic RRKM rate is quite close to simulation results. For the N = 2 chain, a novel approach to the extraction of decay rate constants based on the classical spectral theorem is implemented. Good agreement between the RRKM and dynamical rate constants is obtained for N = 2 despite the fact that the reactant phase space contains a significant fraction of relatively short-lived trajectories.