This work reports numerical experiments intended to clarify the internal
equilibration process in large molecules, following vibrational excitation. A model
of an amorphous and oxygenated hydrocarbon macromolecule (∼ 500
atoms)–simulating interstellar dust—is built up by means of a chemical simulation
code. Its structure is optimized, and its normal modes determined. About 4.5 eV
of potential energy is then deposited locally by perturbing one of the C–H
peripheral bonds, thus simulating the capture of a free H atom by a dangling C
bond. The ensuing relaxation of the system is followed for up to 300 ps, using a
molecular mechanics code. When steady state is reached, spectra and time
correlation functions of kinetic energy and bond length fluctuations indicate
that most normal modes have been activated, but the motion remains
quasi-periodic or regular. By contrast, when the molecule is violently excited
or embedded in a thermal bath (modelled by Langevin dynamics), the
same markers clearly depict chaotic motions. Thus it appears that even
such a large system of oscillators is unable to provide the equivalent of
a thermal bath to any one of these, barring strong resonances between
some of them. This conclusion is of consequence for the interpretation of
astronomical infrared spectra.
Collateral numerical experiments show that (a)
relaxation times increase as perturbation energy decreases by spreading
through the system; (b) energy deposited in the highest-frequency modes
does not relax preferentially into the lowest-frequency modes but follows
specific paths determined by near resonances and coupling strength; (c)
energy deposited in the lowest-frequency modes is able to flow up the
whole frequency ladder, albeit less easily than in the opposite direction.