Susceptibility of quasiclassical Brownian motion in harmonic nonlinear potentials

This work sets the exact equations for the quasiclassical response function and susceptibility of a Brownian particle immersed in a bath of quantum harmonic oscillators driven by nonlinear harmonic potentials. A delta force perturbation gives rise to a response whose susceptibility is the combination of a linear term, owner of the harmonic oscillator, plus a nonlinear one involving an integral equation. It is provided a recursion method to find its solutions based on functional equations in the Banach space. The ODE for the response function is a highly nonlinear damped nonautonomous Duffing equation for which the aforementioned method is used to get its solution.