Graph-based analysis of kinetics on multidimensional potential-energy surfaces

Mean first passage times reconstruct the slowest relaxations in potential energy landscapes of nanoclusters

Relaxation modes are the collective modes in which all probability deviations from equilibrium states decay with the same relaxation rates. In contrast, a first passage time is the required time for arriving for the first time from one state to another. In this paper, we discuss how and why the slowest relaxation rates of relaxation modes are reconstructed from the first passage times. As an illustrative model, we use a continuous-time Markov state model of vacancy diffusion in KCl nanoclusters. Using this model, we reveal that all characteristics of the relaxations in KCl nanoclusters come from the fact that they are hybrids of two kinetically different regions of the fast surface and slow bulk diffusions. The origin of the different diffusivities turns out to come from the heterogeneity of the activation energies on the potential energy landscapes. We also develop a stationary population method to compute the mean first passage times as mean times required for pair annihilations of particle-hole pairs, which enables us to obtain the symmetric results of relaxation rates under the exchange of the sinks and the sources. With this symmetric method, we finally show why the slowest relaxation times can be reconstructed from the mean first passage times.

Slowest kinetic modes revealed by metabasin renormalization

Understanding the slowest relaxations of complex systems, such as relaxation of glass-forming materials, diffusion in nanoclusters, and folding of biomolecules, is important for physics, chemistry, and biology. For a kinetic system, the relaxation modes are determined by diagonalizing its transition rate matrix. However, for realistic systems of interest, numerical diagonalization, as well as extracting physical understanding from the diagonalization results, is difficult due to the high dimensionality. Here, we develop an alternative and generally applicable method of extracting the long-time scale relaxation dynamics by combining the metabasin analysis of Okushima et al. [Phys. Rev. E 80, 036112 (2009)PLEEE81539-375510.1103/PhysRevE.80.036112] and a Jacobi method. We test the method on an illustrative model of a four-funnel model, for which we obtain a renormalized kinematic equation of much lower dimension sufficient for determining slow relaxation modes precisely. The method is successfully applied to the vacancy transport problem in ionic nanoparticles [Niiyama et al., Chem. Phys. Lett. 654, 52 (2016)CHPLBC0009-261410.1016/j.cplett.2016.04.088], allowing a clear physical interpretation that the final relaxation consists of two successive, characteristic processes.