(1) Combining Monte Carlo method, Molecular Dynamics method and the statistical perturbation method, the force fields on the free energy hypersurface were calculated for several interesting chemically reacting systems in solution.
(2) The numerical algorythm of stochastic partial differential equation that the head investigator completed in 2001, was coded to develop a program to solve a simultaneous equations of motion for a quantum-classical mixed reacting system.
(3) We have clarified the microscopic reaction dynamics in solution, especially for Menshutkin reaction, and explained some typical characteristics by calculating the free energy profile along the intrinsic reaction coordinate, which was executed by the free energy gradient method.
(4) At the end of the research term, several results with respect to the free energy profiles for several chemical reactions in solution, were published in a few important meetings, e.g., the annual meetings of Japan Chemical Society. Further, they were collected and published as several research articles as described in the following section, i.e., 11 REFERENCE.