Monte Carlo integration relies on sampling from a certain distribution of random numbers. Importance sampling is a category of Monte Carlo approach which is widely used to compute multi-dimensional integral. In this approach, one first identifies a sampling distribution which resembles the target integrand. One then proceeds to sample the numbers from this distribution, and evaluates the modified integrand at these numbers. It is connected to umbrella sampling, and is popular among computational physicists. In importance sampling, the following approximation is made
where is a normalized sampling distribution. We use Monte Carlo technique to evaluate the classical flux-side correlation function, . is used to compute classical rates, especially those that proceed at high temperatures or involve heavy reacting particles. is defined as
where potential energy , and represents the shift in coupling.
We want to evaluate the by using the method of Monte Carlo integration. For this, we first need to identify efficient sampling distribution. The normalized Gaussians from which we wish to sample , and are defined as,
We need not sample , as this coordinate will be fixed by the properties of delta function. We therefore identify . Therefore, after using the properties of delta function to integrate out -degree of freedom, we can simply rewrite the integral as,
is the Heaviside step function, and is obtained by propagating the particle classically using velocity Verlet algorithm. The reaction coordinate is taken to be -mode.
is computed for two different choices of dividing surface parameter . In both cases, at long time, the value of are identical. This suggests that classical rate is independent of the choice of dividing surface. The TST approximation to the classical rate is obtained in the short-time limit of , which obviously depends on the choice of .
Note: I prepared this problem for Advanced Kinetics exercise at ETH Zurich in Spring 2020.