A mathematical model is introduced for weakly nonlinear wave phenomena in molecular systems like DNA and protein molecules that includes thermal effects: exchange of heat energy with the surrounding aqueous medium.
The resulting equation is a stochastic discrete nonlinear Schrödinger equation with focusing cubic nonlinearity and "Thermal'' terms modeling heat input and loss: PDSDNLS.
New numerical methods are introduced to handle the unusual combination of a conservative equation, stochastic, and fully nonlinear terms.
Some analysis is given of accuracy needs, and the special issues of time step adjustment in stochastic realizations.
Numerical studies are presented of the effects of thermalization on solitons, including damping induced
self-trapping of wave energy, a discrete counterpart of single-point blowup.