Transfer matrices and solution of the problem of stochastic dynamics of aerosol clusters by Monte Carlo method

A Monte Carlo algorithm based on the use of transfer matrices is developed to describe the stochastic dynamics of the rotational--translational motion of aerosol clusters taking into account fluctuations in the molecular fluxes of the gas medium. In the general case, the cluster is immersed into a rarefied gas medium, the temperatures of its surfaces may differ from the temperature of the surrounding gas, for example, due to absorption of visible and infrared radiation. The motion of the cluster is described based on Langevin motion equations. The algorithm allows one to calculate parameters of the probability distribution of a six-dimensional vector consisting of components of the momentum and angular momentum vectors transmitted to the cluster by molecular flows. The numerical method allows one to apply preliminary analytical averaging modulo velocities of molecules for both the average values of the components of momentum and angular momentum and their correlation characteristics, which significantly reduces the calculation time.