## Rare transition event with self-consistent theory of large-amplitude collective motion

Kyosuke Tsumura, Yoshitaka Maeda, Hiroyuki Watanabe Analysis Technology Center, Research &amp; Development Management Headquarters,

Fujifilm Corporation, Kanagawa 250-0193, Japan

## Abstract

A numerical simulation method, based on Dang et al. 's self-consistent theory of large-amplitude collective motion, for rare transition events is presented. The method provides a one-dimensional pathway without knowledge of the final configuration, which includes a dynamical effect caused by not only a potential but also kinetic term. Although it is difficult to apply the molecular dynamics simulation to a narrow-gate potential, the method presented is applicable to the case. A toy model with a high-energy barrier and/or the narrow gate shows that while the Dang et al. treatment is unstable for a changing of model parameters, our method stable for it.

Keywords: PACS: 82.20.-w, 87.10.Tf

Rare transition event, Collective path

## 1. Introduction

Transition events with long timescales are of great importance in many fields of molecular science. In particular, examples in biology range from conformational changes in proteins associated with ligand binding [1] to allosteric transitions that occur throughout the global protein domain [2, 3]. Their typical timescales are the order of 10 -5 -1 sec [4]. Such rare transition events are due to the separation between the initial and final configurations in phase space. The separation is caused by the existence of a high energy barrier or a narrow gate in a potential energy surface; schematic examples of such separations are shown in Fig. 1.

Email address: kyosuke.tsumura@fujifilm.com (Kyosuke Tsumura)

Figure 1: Schematic examples of the separations in the potential energy surface. (a) A high energy barrier or (b) a narrow gate separate the initial and final configurations represented by the circle and cross, respectively. The thin lines denote the contours of the potential energy surfaces.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0364v3-figure-379e4e36-5fa0-4b2b-8064-aca8701fc719.jpg)

Numerical simulation is an essential tool to elucidate the rare transition event and obtain the final configuration. There exist two kinds of approaches to the problem. The first one is to find a pathway determined only by the geometrical feature of the potential energy together with some criteria. The second one is to obtain a trajectory which respects both the potential energy and the kinetic energy. In this work, we employ the latter one because the rare transition events are generated by the potential and kinetic effects.

A standard method for the latter one is molecular dynamics (MD). However, the computational time requirement of conventional MD makes this method insufficient for simulation of these rare transition events. A number of efficient methods that use artificial forces to overcome high energy barriers have been proposed [5, 6, 7, 8]. Most of the previous methods, however, are not efficient for rare transition events caused by a narrow gate. This is because the transition rate through the narrow gate is determined by a frequency factor and has a weak (at most power-law) dependence on the magnitude of the artificial forces as shown in chemical kinetics. Both the high energy barrier and the narrow gate are of significance for rare transition events especially in biology [4].

In this paper, we describe the development of a method to simulate rare transition events caused by not only the high energy barrier but also the narrow gate; we are interested in transitions between local minima in a potential energy surface. Typical thermally fluctuating motions are confined in space near the initial configuration; the motions have complex multi-

dimensional structures. On the other hand, the motion in the rare transition event moves along a simple sharp curve for the following reasons. For the narrow gate, coordinates of the system are limited to specific values as a result of the narrowness of the potential energy surface. For the high energy barrier, the coordinates must obtain large momenta to overcome the barrier, and the momenta restrict the direction of motion to a specific one. From this discussion, the method to be developed is one that can separate one relevant coordinate, which parametrizes the curve, from other coordinates orthogonal to the relevant one. We call the coordinate a collective coordinate and the curve a collective path . Many construction methods for the collective path have been proposed over the past five or more decades [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19]. This approach includes both the potential energy and the kinetic energy as we shall see in Appendix A.

The self-consistent theory of large-amplitude collective motion developed by Dang et al. [15, 16, 17, 18] is one of the powerful methods, which is a pioneering work in applications to molecular systems. Central to their theory are decoupling conditions given by four equations whose solution gives the collective path. Because exact decoupling is not anticipated for realistic systems, that is, the four equations cannot have a simultaneous solution in general, the problem becomes that of defining a constructive procedure for the approximate collective path. In Refs. [15, 16, 17, 18], two approximations based on a subset of the four equations have been proposed.

However, which approximations in the theory of Dang et al. give a transition event between local minima in a potential energy surface depends highly on the details of the potential parameters, as will be demonstrated in Sec. 2.3. Owing to the lack of the robustness in the previous works, defining a unique approximation to the decoupling conditions and a method for evaluating the quality of it is a central issue for development of the collective path theory.

A natural improvement of the approximations proposed by Dang et al. is a construction of the collective path as an optimal solution to all the four decoupling condition equations rather than a subset of them. To manipulate this approximation and evaluate the quality of the resultant collective path, we introduce a unique scalar quantity that measures the deviations from the equalities in the four decoupling condition equations. We then minimize this quantity to formulate the proper construction method of the collective path. The resultant method is based on a first-order differential equation whose solubility algorithm requires calculation of the first and second derivatives of

the potential energy at each step for a stepwise construction of the collective path. We show this method can robustly construct the collective path even if the initial and final configurations are separated by a high energy barrier or narrow gate as shown in Fig. 1.

This paper is organized as follows: In Sec. 2, we introduce the decoupling conditions in the theory of Dang et al. [15, 16, 17, 18] and demonstrate that the two approximations cannot work as a robust construction method of the collective path. A detailed derivation of the decoupling conditions is presented in Appendix A. In Sec. 3, we define the scalar quantity and then formulate the proper method for the construction of the collective path by minimizing this quantity. In Sec. 4, we use two simple Hamiltonians to demonstrate the ability of our method to construct the collective path. We devote Sec. 5 to a summary and concluding remarks.

## 2. Preliminary

In this section, we present a brief account of the self-consistent theory of large-amplitude collective motion developed by Dang et al. [15, 16, 17, 18]. First, we present an explicit form of the decoupling conditions, and prove that there exists no simultaneous solution to them in the generic case. Then, we introduce the two approximations to the decoupling conditions that Dang et al. [15, 16, 17, 18] have proposed. Furthermore, we demonstrate that whether the two approximations can detect the collective path connecting local minima depends highly on the shape of a potential energy surface. We note that the proof of the absence of the solution and the demonstration of the lack of the robustness property are original arguments of this work.

## 2.1. Decoupling conditions

Because we are interested in the applications to molecular systems, we start with a Hamiltonian defined by

$$H = \sum _ { i = 1 } ^ { N } \, \frac { 1 } { 2 \, m _ { i } } \, p _ { i } ^ { 2 } + V ( x _ { 1 }, \, \cdots, \, x _ { N } ),$$

where N is the total number of degrees of freedom, x i , p i , and m i denote the coordinate, conjugate momentum, and mass of the i -th degree of freedom, respectively, and V ( x , 1 · · · , x N ) is the potential energy.

We introduce the decoupling conditions that define the collective path in the Hamiltonian (1). A detailed derivation of the decoupling conditions is presented in Appendix A. When we write the collective path as

$$x _ { i } = x _ { i } ( q ),$$

where q denotes the collective coordinate, the decoupling conditions to determine the explicit q -dependence of x i [15, 16, 17, 18] are

$$\frac { \mathrm d x _ { i } } { \mathrm d q } = \varphi _ { i } / \sqrt { m _ { i } },$$

where ϕ i are given as a solution to the following equations:

$$\sum _ { j = 1 } ^ { N } \left ( V _ {, i j } - \lambda \, \delta _ { i j } \right ) \varphi _ { j } \ = \ 0, \quad & & ( 4 ) \\ \nu, \quad & & \nu \searrow$$

$$V _ {, i } \ = \ \gamma \, \varphi _ { i },$$

$$\sum _ { i = 1 } ^ { N } \varphi _ { i } \, \varphi _ { i } \ = \ 1.$$

$$V _ {, i } \ \equiv \ \frac { 1 } { \sqrt { m _ { i } } } \, \frac { \partial V ( \boldsymbol x ) } { \partial x _ { i } }, \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$V _ {, i j } \ \equiv \ \frac { \cdot } { \sqrt { m _ { i } } \sqrt { m _ { j } } } \frac { \partial ^ { 2 } V ( \boldsymbol x ) } { \partial x _ { i } \, \partial x _ { j } },$$

with x = ( x , 1 · · · , x N ). We note that Eqs. (3)-(6) show that ϕ / i √ m i is a tangent vector to the collective path x = x ( q ) at the collective coordinate q , ϕ i should be an eigenvector of the matrix V ,ij , be proportional to the vector V ,i , and be normalized and that λ and γ denote an eigenvalue of the eigenvector ϕ i and a proportionality factor, respectively. By solving the decoupling conditions given by Eq. (3) together with Eqs. (4)-(6), we can construct the collective path x = x ( q ). We use the obtained x ( q ) to reduce the Hamiltonian (1) into a Hamiltonian that governs the dynamics of the collective coordinate q [15, 16, 17, 18].

where

## 2.2. Proof for absence of solution to decoupling conditions

We emphasize that for a generic potential energy function V ( x ), there exists no ϕ i that satisfies Eqs. (3)-(6), except for the local minima or saddle points defined by V ,i = 0. To explain this fact in an explicit manner, first we identify ( x , ϕ , λ, γ i i ) with would-be independent variables and consider an infinitesimal expansion of these variables around ( x , ϕ , λ, γ i i ), which satisfies Eqs. (4)-(6). Substituting x ′ i = x i + d x i , ϕ ′ i = ϕ i + d ϕ i , λ ′ = λ + d λ , and γ ′ = γ +d γ into Eqs. (4)-(6), we have equations governing (d x , i d ϕ , i d λ, d γ ), which are necessary for ( x , ϕ , λ , γ ′ i ′ i ′ ′ ) to remain the solution to Eqs. (4)-(6). Then, by dividing the equations by d q and incorporating Eq. (3) with them, we reduce the decoupling conditions (3)-(6) into a compact set of equations given by

$$\frac { \mathrm d } { \mathrm d q } ( x _ { i }, \varphi _ { i }, \lambda, \gamma ) = ( X _ { i } / \sqrt { m _ { i } }, \sum _ { a = 1 } ^ { N - 1 } \varphi _ { \perp i } ^ { a } \Phi _ { a }, \Lambda, \Gamma ),$$

where ϕ a ⊥ i with a = 1 , · · · , N -1 denote a complete set of the vectors orthogonal to ϕ i with the normalizations

$$\sum _ { i = 1 } ^ { N } \varphi _ { \perp i } ^ { a } \varphi _ { \perp i } ^ { b } \ = \ \delta ^ { a b },$$

and Y ≡ t ( X , i Φ a , Λ Γ) are defined as the solution to ,

$$A \, Y - b = 0,$$

with where

$$V _ {, i j k } \equiv \frac { 1 } { \sqrt { m _ { i } } \sqrt { m _ { j } } \sqrt { m _ { k } } } \frac { \partial ^ { 3 } V ( \boldsymbol x ) } { \partial x _ { i } \, \partial x _ { j } \, \partial x _ { k } }.$$

Here, we have used the Einstein summation convention for the dummy indices. We note that the matrix A is a (3 N ) × (2 N + 1) matrix. Because

$$A \ \equiv \ \left ( \begin{array} { c c c } V _ {, i j k } \varphi _ { k } & ( V _ {, i j } - \lambda \, \delta _ { i j } ) \, \varphi _ { \perp j } ^ { a } & - \varphi _ { i } & 0 \\ V _ {, i j } & - \gamma \, \varphi _ { \perp i } ^ { a } & 0 & - \varphi _ { i } \\ \delta _ { i j } & 0 & 0 & 0 \end{array} \right ), \quad \ \ ( 1 2 ) \\ \cdot \quad t / \sim \ \sim \ \, \end{array}$$

$$b \ \equiv \ ^ { t } ( 0, 0, \varphi _ { i } ),$$

the number of the equations is obviously larger than that of the variables, Eq. (11) has no solution with respect to Y = ( t X , i Φ a , Λ Γ) in general, and , hence the collective path x ( q ) also cannot be constructed as a solution to Eq. (9).

/negationslash

We prove that Eq. (11) has a definite solution of Y = ( t X , i Φ a , Λ Γ) , when x i agrees with x eq i , satisfying V ,i ( x eq ) = 0. We write Y obtained at x eq i as Y eq = ( t X eq i , Φ eq a , Λ eq , Γ eq ). To obtain an explicit form of Y eq , first we construct ϕ i , λ , γ , and ϕ a ⊥ i which are necessary to build the matrix elements of A and b in Eq. (11). Substituting V ,i ( x eq ) = 0 into Eq. (5), we have γ = 0. Furthermore, because Eq. (4) requires ϕ i and λ to be eigenvectors and eigenvalues of V ,ij ( x eq ), respectively, we have ϕ i = φ 0 i and λ = λ 0 . Here, we have introduced φ µ i and λ µ with µ = 0 , · · · , N -1 as the eigenvectors and eigenvalues of V ,ij ( x eq ), respectively. We suppose that λ 0 = λ a for any a = 1 , · · · , N -1. Using the eigenvectors φ a i orthogonal to φ 0 i , we have ϕ a ⊥ i = φ a i . Then, we use ( x , ϕ , λ, γ, ϕ i i a ⊥ i ) = ( x eq i , φ , λ , 0 i 0 0 , φ a i ) to construct A and b . After substituting these into Eq. (11), we come to the solution that

$$X _ { i } ^ { \text{eq} } \ = \ \phi _ { i } ^ { 0 },$$

$$\Phi _ { a } ^ { \text{eq} } & \ = \ \frac { 1 } { \lambda ^ { 0 } - \lambda ^ { a } } \, V _ {, i j k } ( x ^ { \text{eq} } ) \, \phi _ { i } ^ { 0 } \, \phi _ { j } ^ { 0 } \, \phi _ { k } ^ { a }, \\ \Lambda ^ { \text{eq} } & \ = \ V _ {, i j k } ( x ^ { \text{eq} } ) \, \phi _ { i } ^ { 0 } \, \phi _ { i } ^ { 0 },$$

$$\Lambda ^ { \text{eq} } \ = \ V _ {, i j k } ( \boldsymbol x ^ { \text{eq} } ) \, \phi _ { i } ^ { 0 } \, \phi _ { j } ^ { 0 } \, \phi _ { k } ^ { 0 },$$

$$\Gamma ^ { \text{eq} } \ = \ \lambda ^ { 0 }.$$

Indeed, Y = ( t X , i Φ a , Λ Γ) has been constructed as a definite solution to , Eq. (11) at a stationary point.

## 2.3. Approximations to decoupling conditions

It is a significant task to obtain an approximate solution to Eqs. (3)-(6). In Refs. [17, 18], Dang et al. proposed two methods. One respects Eqs. (4)-(6), but not (3), which is an approximation 1 (Approx1) method. The other takes into account Eqs. (3), (4), and (6), but not (5), which we call an approximation 2 (Approx2) method.

We examine whether the Approx1 method or the Approx2 method can construct the collective path connecting local minima, by using the following two-dimensional Hamiltonian ( N = 2):

$$H = \frac { 1 } { 2 } \, ( p _ { 1 } ^ { 2 } + p _ { 2 } ^ { 2 } ) + V _ { \text{MB} } ( x _ { 1 }, \, x _ { 2 } ).$$

Figure 2: The potential energy surfaces made by the M¨ller-Brown potential [20], i.e., u (a) the modified one V ( a ) MB and (b) the original one V ( b ) MB , together with the collective paths extracted by the Approx1 and Approx2 methods. The circles denote the local minima and the squares the saddle points. The solid (dashed) line represents the collective path calculated by the Approx1 (Approx2) method, which starts from one of the local minima. The thin lines denote the contours of the potential energy surfaces.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0364v3-figure-1b927465-fcfe-4d2f-912f-6c302d575e64.jpg)

Here, V MB ( x , x 1 2 ) is the M¨ller-Brown potential [20], u

$$V _ { \text{MB} } ( x _ { 1 }, x _ { 2 } ) = \sum _ { \alpha = 1 } ^ { 4 } C ^ { \alpha } \, \exp \left [ \frac { 1 } { 2 } \, \sum _ { i, j = 1 } ^ { 2 } \left ( x _ { i } - \bar { x } _ { i } ^ { \alpha } \right ) M _ { i j } ^ { \alpha } \left ( x _ { j } - \bar { x } _ { j } ^ { \alpha } \right ) \right ], \quad \ ( 2 0 )$$

where ( C , x α ¯ α i , M α ij ) with α = 1 , · · · , 4 and i, j = 1 , 2 denote parameters. The parameters are set equal to [20]

$$C ^ { \alpha } \ & = \ ( - 2 0 0, \, - 1 0 0, \, - 1 7 0, \, 1 5 ), \quad \ ( 2 1 ) \\ \bar { \tau } ^ { \alpha } \ & = \ ( 1 \, \ 0 \, \ = 0 \, 5 \, \ = 1 \, \L )$$

$$\bar { x } _ { 2 } ^ { \alpha } \ = \ ( 0, \, 0. 5, \, 1. 5, \, 1 ),$$

$$\bar { x } _ { 1 } ^ { \alpha } \ & = \ ( 1, \, 0, \, - 0. 5, \, - 1 ), \\ \bar { r } _ { \mu } ^ { \alpha } \ & = \ ( 0 \, 0 \, 5 \, 1 \, 5 \, 1 )$$

$$M _ { 1 1 } ^ { \alpha } \ & = \ ( - 2, \, - 2, \, - 1 3, \, 1. 4 ), \\ M _ { \alpha } ^ { \alpha } \ & = \ M _ { \alpha } ^ { \alpha } = \left ( 0 \ 0 \ 1 1 \ 0 \ 6 \right )$$

$$M _ { 1 2 } ^ { \alpha } \ = \ M _ { 2 1 } ^ { \alpha } = ( 0, \, 0, \, 1 1, \, 0. 6 ),$$

$$M _ { 2 2 } ^ { \alpha } \ = \ ( - 2 0, \, - 2 0, \, - 1 3, \, 1. 4 ).$$

The potential energy produced by this parameter set is denoted as V ( b ) MB and the set modified by altering C 3 from -170 to -120 is denoted V ( a ) MB .

In Fig. 2, we summarize four collective paths that are extracted from the Hamiltonians with V ( a ) MB and V ( b ) MB by the Approx1 and Approx2 methods. In

Figure 3: All the paths derived by the Approx1 method for (a) V ( a ) MB and (b) V ( b ) MB .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0364v3-figure-1dd1a620-3672-49bf-8720-911a3c4ffa1f.jpg)

Refs. [17, 18], the Approx1 method was thought to be a better method than the Approx2 method. For V ( a ) MB , indeed, the Approx1 method can produce the collective path between local minima, while the Approx2 method is not the case. For V ( b ) MB , however, the approximation that extracts the transition event is not the Approx1 method but the Approx2 method. Figure 3 depicts all paths based on the Approx1 method. We note that, in contrast to V ( a ) MB , V ( b ) MB has no single path connecting the local minima. It turns out that whether the Approx1 method or the Approx2 method can extract the transition event sensitively depends on the shape of the potential energy surface, and hence both the Approx1 and Approx2 methods have not been a robust construction method of the collective path.

## 3. Improvement of approximations to decoupling conditions in self-consistent theory of large-amplitude collective motion

In this section, we develop the numerical simulation method that can robustly construct the collective path connecting local minima, which is an improvement of the Approx1 and Approx2 methods. A key point for the development of the method is the construction of Y = ( t X , i Φ a , Λ Γ) as an , optimal solution to Eq. (11).

By constructing the optimal solution Y and substituting Y into Eq. (9), we can obtain a first-order differential equation to determine ( x , ϕ , λ, γ i i ) as a solution to an initial-value problem.

We shall construct a procedure to obtain the optimal solution to Eq. (11).

To this end, first we define a vector ∆, which represents a deviation from the equality in Eq. (11), that is,

$$\Delta \equiv A Y - b.$$

Here, we treat the problem to minimize the magnitude of ∆ to derive an equation for the construction of the optimal Y .

## 3.1. Measure of deviation from exact decoupling

We introduce three new vectors d i , e i , and f i to represent ∆ as a form of the modification into ϕ i in the elements of A and b . Indeed, by parameterizing ∆ as

$$\Delta = \left ( \begin{array} { c } G _ { i j } \left ( d _ { j } - \varphi _ { j } \right ) \\ \Gamma \left ( e _ { i } - \varphi _ { i } \right ) \\ f _ { i } - \varphi _ { i } \end{array} \right ),$$

with

$$G _ { i j } \equiv \Lambda \, \delta _ { i j } - V _ {, i j k } \, X _ { k },$$

we can convert Eq. (27) into the same form as Eq. (11):

$$\tilde { A } Y - \tilde { b } = 0,$$

where

$$\begin{array} {$$

$$\begin{array} {$$

Here, ˜ and ˜ are obtained by the replacement of A b ϕ i in the first and second rows of A and the third row of b to d i , e i , and f i , respectively. Because the condition d i = e i = f i = ϕ i reproduces the exact decoupling conditions, we can reach the following natural definition of the scalar quantity that measures the deviation from exact decoupling:

$$Q = \sum _ { i = 1 } ^ { N } \, \left ( w _ {$$

where w d , w e , and w f denote weight parameters to determine which decoupling conditions (3)-(5) are respected. Because we are interested in the case where all of the decoupling conditions (3)-(5) are respected equivalently, we must set

$$( w _ { d }, w _ { e }, w _ { f } ) = ( 1$$

By definition, the set of ( w , w , w d e f ) = (1 , 1 , 0) corresponds to the Approx1 method, while the set of ( w , w , w d e f ) = (1 , 0 , 1) leads to the Approx2 method.

To make clear a geometrical meaning of Q , we rewrite Q as

$$Q = { ^ { t } } \Delta \, C \, \Delta = { ^ { t$$

This representation means that Q is identical to the norm, i.e., the square of the length of the vector ∆ = AY -b measured on the metric C given by

$$C \equiv \left ( \begin{array} { c c c } \left [ ^ { t } ( G ^ { - 1 } ) \, G ^ { - 1 } \right ] _ { i j } & 0 & 0 \\ 0 & \Gamma ^ { - 2 } \, \delta _ { i j } & 0 \\ 0 & 0 & \delta _ { i j } \end{array} \right ).$$

It is obvious that C has a dependence on ( X , i Λ Γ) that is not yet deter-, mined.

As shown in Eqs. (15)-(18), Y = ( t X , i Φ a , Λ Γ) can be obtained without , reference to C at the local minima and saddle points x eq i . If we take x eq i as a starting point to detect the collective path, we can carry out a straightforward construction of C by substituting ( X , i Λ Γ) in Eqs. (15), (17), and (18) into , Eq. (36). Here, we refer to C constructed at the starting point x eq i as C eq .

In this work, we use this metric C eq as a substitute for C in the course of the detection of the collective path;

$$C = C ^ { \text{eq} }.$$

Thanks to the metric C being constant, Q is a function bilinear with respect to Y , which is easily minimized. It is left to future work to carry out the nonlinear optimization of Q with respect to Y by taking into account the Y -dependence in C .

## 3.2. Basic equation of proper construction method of collective path

By applying the minimization (stationary) condition ∂Q/∂Y = 0 to Q with C = C eq , we arrive at

$$( ^ { t } A \, C ^ { \text{eq} } \, A ) \, Y - ^ { t } A \, C ^ { \text{eq} } \, b = 0.$$

We note that in contrast to Eq. (11), Eq. (38) can have a definite solution in general because t AC eq A is a (2 N +1) × (2 N +1) matrix and the number of the equations agrees with that of the variables. We stress that Eq. (9) combined with Y that is a solution to Eq. (38) is the main equation for the construction of the collective path, which is one of the results in this work.

A few remarks are in order here: (i) The basic equation of our method given by the set of Eqs. (9) and (38) is a first-order differential equation with respect to ( x , ϕ , λ, γ i i ), which can be solved with a single iteration for the stepwise construction of the collective path by repeating the calculation of the elements of the matrix A . In fact, we can construct A without an iteration procedure, because the elements of A are given by only the first, second, and third derivatives of the potential energy with respect to the coordinate variables, which can be calculated at a given structure of the system. Furthermore, we note that the elements associated with the third derivatives of V ( x ) in the matrix A in Eq. (38) can be calculated in the same computational cost as that in the case of the second derivative terms because the third derivative terms enter into A as a form of V ,ijk ϕ k . Indeed, we can construct V ,ijk ϕ k from the two second derivatives on the basis of the formula

$$V _ {, i j k } ( x ) \, \varphi _ { k } \sim \frac { V _ {, i j } ( x + n \, \Delta x ) - V _ {, i j } ( x - n \, \Delta x ) } { 2 \, \Delta x },$$

with n ≡ ( ϕ / 1 √ m , 1 · · · , ϕ N / √ m N ). Here, ∆ x denotes an infinitesimal value. We note that the calculations of the two second derivatives can be run in parallel. Thanks to Eq. (39), the resultant method requires only the numerical calculation of the first and second derivatives in each step of the single iteration. This property reduces the computational costs considerably. (ii) As mentioned in the final paragraph in Sec. 2.2, we suppose that the initial value of ( x , ϕ , λ, γ i i ) is set equal to ( x eq i , φ , λ , 0 i 0 0). Here, we present an account of the eigenvalue λ 0 and eigenvector φ 0 i . In this work, we adopt the smallest eigenvalue as λ 0 . This is because we are interested in the rare transition event of molecular systems where the amplitude of the collective

Figure 4: The collective paths constructed by the method based on the set of Eqs. (9) and (38). The potential energy surfaces are the same as those in Fig. 2, and the solid line represents the collective path.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0364v3-figure-6096dbfe-7fc7-4cb1-84c4-2e8ad9f756ca.jpg)

coordinate is large compared to those of the other coordinates, which can be translated to the requirement to use the smallest eigenvalues of V ,ij . We note that the direction of the construction of the collective path from the initial configuration depends on the sign of φ 0 i , because φ 0 i gives the step from the initial structure to the next one.

## 4. Applications of self-consistent theory of large-amplitude collective motion to simple Hamiltonians

In this section, first we examine our method's ability to construct the collective path connecting local minima. Then, we elucidate the ability of the method to extract the rare transition event by using potential energy surfaces with either a high energy barrier or a narrow gate to separate the two local minima as shown in Fig. 1.

## 4.1. M¨ uller-Brown potential

To examine the robustness property of the numerical simulation method based on the set of Eqs. (9) and (38), we utilize the Hamiltonians with the M¨ uller-Brown potential [20], which are the same as those presented in Fig. 2.

In Fig. 4, we show the resultant collective paths. Both of them are shown as trajectories that connect almost all the local minima and saddle points. As that it is not the case in the collective paths obtained with the Approx1

and Approx2 methods, shown in Fig. 2, we conclude that our method is a robust construction method of the collective path describing the transition between local minima, in contrast to the Approx1 and Approx2 methods.

We now make one remark here. The collective paths depicted in Fig. 4 pass close to but not exactly through the saddle points. This is a critical problem for the studies of minimum energy path [21, 22]. However, the path studied in this paper is the one including not only the potential energy but also the kinetic energy. In this case, that the path passes exactly through the saddle points is not necessarily required. A further examination of this point is a future problem.

## 4.2. Barrier and gate potential

We consider a Hamiltonian given by

$$H = \frac { 1 } { 2 } \, ( p _ { 1 } ^ { 2 } + p _ { 2 } ^ { 2 } ) + V _ { \text{BG} } ( x _ { 1 }, \, x _ { 2 } ).$$

Here, V BG ( x , x 1 2 ) denotes the barrier and gate potential defined by

$$V _ { \text{BG} } ( x _ { 1 }, x _ { 2 } ) = - \frac { 1 } { 2 } \, a \, x _ { 1 } ^ { 2 } + \frac { 1 } { 4 } \, b \, x _ { 1 } ^ { 4 } + ( x _ { 2 } - c \, x _ { 1 } ) ^ { 2 } \left [ \alpha \, x _ { 1 } ^ { 2 } + \beta \, e ^ { - \gamma x _ { 1 } ^ { 2 } } + \delta \right ], \ \ ( 4 1 )$$

with ( a, b, c, α, β, γ, δ ) being parameters. The parameters ( a, b, c ) determine locations and energy values of stationary points. Indeed, this potential energy surface has two local minima and one saddle point located at ( ± √ a/b, ± c √ a/b ) and (0 , 0), respectively, and the activation energy in the transition between the two local minima is given by a / 2 (4 b ). The parameters ( α, β, γ, δ ) determine the narrowness of the pathway connecting the local minima and the saddle point.

To show the ability of the method developed to simulate the rare transition event, we introduce three potential energy surfaces: a potential energy surface containing neither a high energy barrier nor a narrow gate V ( a ) BG , a surface with a high energy barrier V ( b ) BG , and a surface with a narrow gate V ( c ) BG . The parameters ( a, b, c, α, β, γ, δ ) which create V ( a ) BG , V ( b ) BG , and V ( c ) BG read (80 , 80 , 3 , 0 0045 . , 10 , 10 , 1), (240 , 240 , 3 , 0 , 0 , 10 , 20), and (80 , 80 , 3 , 0 31 . , 690 , 10 , 9 7), . respectively. All the potential energies have the same locations of the local minima as ( ± 1 , ± 3).

In Fig. 5, we show the collective paths corresponding to V ( a ) BG , V ( b ) BG , and V ( c ) BG together with trajectories obtained by the conventional MD method.

Figure 5: The collective paths obtained by the method developed in this work and the trajectories obtained by the MD method on the V ( a ) BG , V ( b ) BG , and V ( c ) BG potential energy surfaces. The solid lines in (a-1) and (a-2) ((b-1) and (b-2) or (c-1) and (c-2)) show the collective path and the MD trajectory derived from V ( a ) BG ( V ( b ) BG or V ( c ) BG ), respectively. In (a-1) and (a-2) ((b-1) and (b-2) or (c-1) and (c-2)), the thin lines denote the contours of V ( a ) BG ( V ( b ) BG or V ( c ) BG ), the circles the local minima, and the squares the saddle point.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0364v3-figure-b7abfc14-b2ac-4e77-b989-befd00081551.jpg)

We find that the collective paths provide pathways connecting the two local minima even if a high energy barrier or narrow gate exists between the two local minima. Conversely, the MD method applied to V ( a ) BG can provide a trajectory near the collective path, but the trajectories obtained for the V ( b ) BG and V ( c ) BG cases are seen to be trapped near one local minimum and cannot reach the other local minimum. We note that the qualitatively same results are obtained with less symmetric potentials, e.g., V BG ( x , x 1 2 )+ x 3 1 ; in Appendix B, the results for the potentials are discussed. We stress that in contrast to the MD method, our method works well as a simulation technique for the rare transition event.

## 5. Summary and concluding remarks

Basic notion adopted in this work is that simulation of rare transition events can reduce to construction of a collective path given as a simple sharp curve along which only one coordinate increases if the initial coordinates are adequately chosen. To detect the collective path, we have developed a method to separate a collective coordinate, which parametrizes the collective path, from the other coordinates orthogonal to the collective coordinate. The method has been formulated by constructing the collective path as an optimal solution to the decoupling conditions in the self-consistent theory of large-amplitude collective motion proposed by Dang et al. [15, 16, 17, 18]. A basic equation of the method is given by the set of Eqs. (9) and (38). We have found that the equation is a first-order differential equation that determines the collective path as a solution to an initial-value problem from an arbitrary point in the collective path with the use of the first and second derivatives of the potential energy with respect to the atomic coordinates. By using the M¨ller-Brown potential [20], we have demonstrated that the u method, which is an improvement of the theory of Dang et al. , can robustly construct the collective path connecting the local minima in the potential energy surface. We have used a simple potential energy surface to show that our method uses only the initial configuration to construct a pathway to the final configuration, even if the configurations are separated by a high energy barrier and/or a narrow gate.

We consider that our study can stimulate discussions about the rare transition events in various systems of interest. In fact, the method presented in this work can be applied to simulate the rare transition events of generic

molecular systems whose potential energy surfaces are well-defined. For example, by using a potential energy for biomolecules, e.g., the generalized amber force field [23], we will demonstrate [24] that our method can be applied to describe allosteric transitions [2, 3] which are typical rare transition events for protein systems. Furthermore, if we utilize potential energies calculated numerically using computer programs developed on the basis of quantum chemistry, we can describe a chemical reaction as a rare transition event. In fact, we will show [25] that combining this method with the density functional theory [26] allows simulation of an intermolecular proton transfer reaction, i.e., molecular dynamics with quantum tunneling. We will also show [27] that our method in combination with time-dependent density functional theory [28] can be applied to simulate photochemical reactions, i.e., the dynamics of molecular systems in electronic excited states.

Finally, we remark again the following study that should be reported in the near future. We apply the non-linear optimization to Eq. (35) to determine the collective path in a full consistent manner and compare it with the collective path obtained in the approximation (37).

## Acknowledgments

The authors are grateful to Dr. Furuya for his encouragement and support.

## Appendix A. Detailed derivation of decoupling conditions

In this section, we present a review of the derivation of the decoupling conditions in the self-consistent theory of large-amplitude collective motion developed by Dang et al. [15, 16, 17, 18]. We focus on the collective path parametrized by one collective coordinate, although Dang et al. [15, 16, 17, 18] have presented a construction method of a multi-dimensional collective path, i.e., the so-called collective surface. Basically, we obey the notation used in Refs. [17, 18] in this section. In the final stage, we convert the derived decoupling conditions into those introduced in Sec. 2.1.

We treat an N -body system whose Hamiltonian is

$$H = \frac { 1 } { 2 } \, \pi _ { \alpha } \, B ^ { \alpha \beta } \, \pi _ { \beta } + V ( \xi ^ { 1 }, \, \cdots, \, \xi ^ { N } ).$$

Here, α and β denote indices that vary from 1 to N ξ , α and π α are a coordinate and its conjugate momentum of the α -th degree of freedom, respectively,

B αβ is a reciprocal mass tensor, and V ( ξ 1 , · · · , ξ N ) is a potential energy. We suppose that B αβ is diagonalized as

$$B ^ { \alpha \beta } = \text{diag} ( 1 / m _ { 1 }, \cdots, 1 / m _ { N } ),$$

with m α being the mass of the α -th degree of freedom.

We use a canonical transformation to obtain conditions that define the collective path. In the canonical transformation, we consider mappings of the form

$$\xi ^ { \alpha } = g ^ { \alpha } ( q ^ { 1 }, \, \cdots, \, q ^ { N } ) \equiv g ^ { \alpha } ( q ),$$

whose inverse relations are given by

$$q ^ { \mu } = f ^ { \mu } ( \xi ^ { 1 }, \cdots, \xi ^ { N } ) \equiv f ^ { \mu } ( \xi ).$$

Using a chain rule relation, we find that g α ( q ) and f µ ( ξ ) satisfy the following relations:

$$\delta ^ { \alpha } _ { \ \beta } \ = \ \partial \xi ^ { \alpha } / \partial \xi ^ { \beta } = g ^ { \alpha } _ {, \mu } ( q ) \, f ^ { \mu } _ {, \beta } ( \xi ),$$

$$\delta ^ { \dot { \mu } } _ { \, \nu } \ = \ \partial q ^ { \mu } / \partial q ^ { \nu } = f ^ { \dot { \mu } } _ {, \, \alpha } ( \xi ) \, g ^ { \dot { \alpha } } _ {, \, \nu } ( q ),$$

where the comma indicates a partial derivative; F , α ( ξ ) ≡ ∂F ξ /∂ξ ( ) α and G , µ ( q ) ≡ ∂G q /∂q ( ) µ . The indices α, β, · · · represent those of the initial coordinates, and µ, ν, · · · those of the final coordinates, although each set has the same range 1 , · · · , N . The momentum p µ conjugate to q µ is given by

$$p _ { \mu } = g _ {, \mu } ^ { \alpha } ( q ) \, \pi _ { \alpha },$$

with inverse

$$\pi _ { \alpha } = f _ {, \alpha } ^ { \mu } ( \xi ) \, p _ { \mu }.$$

The Hamiltonian as described with the canonical coordinates ( q µ , p µ ) is

$$H = \frac { 1 } { 2 } \, p _ { \mu } \, \bar { B } ^ { \mu \nu } ( q ) \, p _ { \nu } + \bar { V } ( q ).$$

Here, ¯ B µν ( q ) is an inverse matrix of ¯ B -1 µν ( q ) given by

$$\bar { B } _ { \mu \nu } ^ { - 1 } ( q ) \equ$$

and ¯ ( V q ) is the potential energy written by q µ as

$$\bar { V } ( q ) \equiv V ( \xi ^ { 1 } = g ^$$

Equations of motion derived from the Hamiltonian (A.9) are

$$\dot { q } ^ { \mu } \ = \ \bar { B } ^ { \mu$$

$$\begin{array} {$$

In the new canonical coordinates, we identify the collective path as follows: We divide the set of ( q µ , p µ ) into two subsets, ( q 1 , p 1 ) and ( q a , p a ) with a = 2 , · · · , N , and suppose this division to ensure that if at time t = 0 both q a = 0 and p a = 0, then q a ( ) = 0 and t p a ( ) = 0. t The time evolution of q 1 and p 1 in the new canonical coordinates can be converted into that of ξ α and π α which takes place on a two-dimensional subspace

$$\begin{array} { r c l } \xi ^ { \alpha } & = & \hat$$

$$\pi _ { \alpha } \ = \ \hat { \pi } _ { \alpha } \equiv f _ {, \alpha } ^ { \mu } ( \hat { \xi } ) \, \hat { p } _ { \mu },$$

with

$$\hat { q } ^ { \mu } \ & \equiv \ ( q ^ { 1 }, \, 0, \, \cdots, \, 0 ), \\ \hat { n } \ & = \ ( n, \, 0 \, \dots \, \, \mathfrak { m } ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$\hat { p } _ { \mu } \ \equiv \ ( p _ { 1 }, \, 0, \, \cdots, \, 0 ).$$

We note that q 1 is the collective coordinate and ˆ ξ α is identical to the collective path parametrized by q 1 .

The requirements ( ˙ q a , ˙ p a ) = (0 , 0) can be compatible with the requirements ( q a , p a ) = (0 , 0) only if the equations

$$0 \ = \ \bar { B } ^ { a 1 } ( \hat { q } ) \, p _ { 1 },$$

$$0 \ = \ - \frac { 1 } { 2 } \, p _ { 1 } \, \bar { B } _ {, a } ^ { 1 1 } ( \hat { q } ) \, p _ { 1 } - \bar { V } _ {, a } ( \hat { q } ), \quad \ \ ( A. 1 9 )$$

are satisfied, as one sees from Eqs. (A.12) and (A.13). Equations (A.18) and (A.19) are equivalent to three conditions

$$\bar { B } ^ { a 1 } ( \hat { q } ) \ = \ 0,$$

$$\bar { V } _ {, a } ( \hat { q } ) \ = \ 0,$$

$$\bar { B } ^ { 1 1 } _ {, \, a } ( \hat { q } ) \ = \ 0.$$

Here, we have assumed that p 1 is not a constant of the motion. We emphasize that these equations determine an explicit form of ˆ . ξ α In fact, if g α , 1 (ˆ) q is obtained as a solution to Eqs. (A.20)-(A.22), we can construct ˆ ξ α = ˆ ( ξ α q 1 ) by solving a differential equation

$$\frac { \mathrm d } { \mathrm d q ^ { 1 } } \hat { \xi } ^ { \alpha } = g _ {, 1 } ^ { \alpha } ( \hat { q } ),$$

which has been derived from the definition (A.14).

From here on, we convert Eqs. (A.20)-(A.22) into convenient forms to obtain an explicit form of g α , 1 (ˆ). q For this purpose, first we introduce two scalar quantities

$$X ^ { ( 0 ) } ( \xi ) \ & \equiv \ V ( \xi ) \\ & = \bar { \ V } ( q ) \equiv \bar { X } ^ { ( 0 ) } _ { \hat { \ } } ( q ),$$

$$X ^ { ( 1 ) } ( \xi ) \ & \equiv \ \stackrel { \theta \, \stackrel { \theta \, \cdot \, } { \, } } { \, } \stackrel { \theta \, \stackrel { \theta \, \cdot \, } { \, } } { \, } \stackrel { \theta \, \stackrel { \theta \, \cdot \, } { \, } \, } \, \stackrel { \theta \, \stackrel { \theta \, \cdot \, } { \, } \, } \, } { \, } \, \stackrel { \theta \, \stackrel { \theta \, \cdot \, } { \, } \, } \, \stackrel { \theta \, \stackrel { \theta \, \cdot \, } { \, } \, } \, } \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \cdot \, } \, } \, } { \, } \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \cdot \, } \, } \, } \, } \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \colon \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \ch | \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \textrel \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \. \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \, \stackrel { \theta \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \hat { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar {'\bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \ } { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { } { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \ \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \ | \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \ { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { - \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \ - \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { { { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \bar { \ }$$

$$\overset { \cdot \, \mathbb { \upsilon } } { \bar { V } } ( q ) \equiv \bar { X } ^ { ( 0 ) } ( q ), & & ( \text{A.24} ) \\ V _ {, \alpha } ( \xi ) \, B ^ { \alpha \beta } \, V _ {, \beta } ( \xi ) \\ \bar { V } _ {, \mu } ( q ) \, \bar { B } ^ { \mu \nu } ( q ) \, \bar { V } _ {, \nu } ( q ) \equiv \bar { X } ^ { ( 1 ) } ( q ). & & ( \text{A.25} ) \\ \overset { \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \colon \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdots } { \dots \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cd. \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdelta \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdets \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cd!.$$

These quantities have an important property given by

$$\bar { X } _ {, a } ^ { ( \sigma ) } ( \hat { q } ) = 0, \ \sigma = 0, \, 1.$$

Here, we present a proof for the above equations. The proof for the case of σ = 0 is trivial as follows:

$$\bar { X } _ {, a } ^ { ( 0 ) } ( \hat { q } ) = \bar { V } _ {, a } ( \hat { q } ) = 0, & & ( \text{A.27} )$$

which is just Eq. (A.21). In the case of σ = 1, we simply compute

$$\begin{array} { r c l } \kappa & \ell & \ell & \ell & \ell & \ell \\ \bar { X } _ {, a } ^ { ( 1 ) } ( \hat { q } ) & = & 2 & \bar { X } _ {, \mu a } ^ { ( 0 ) } ( \hat { q } ) & \bar { X } _ {, \nu } ^ { ( 0 ) } ( \hat { q } ) & \bar { B } ^ { \mu \nu } ( \hat { q } ) + \bar { X } _ {, \mu } ^ { ( 0 ) } ( \hat { q } ) & \bar { X } _ {, \nu } ^ { ( 0 ) } ( \hat { q } ) & \bar { B } ^ { \mu \nu } _ {, a } ( \hat { q } ) \\ & = & 2 & \bar { X } _ {, b a } ^ { ( 0 ) } ( \hat { q } ) & \bar { X } _ {, 1 } ^ { ( 0 ) } ( \hat { q } ) & \bar { B } ^ { b 1 } ( \hat { q } ) + 2 & \bar { X } _ {, 1 a } ^ { ( 0 ) } ( \hat { q } ) & \bar { X } _ {, 1 } ^ { ( 0 ) } ( \hat { q } ) & \bar { B } ^ { 1 1 } ( \hat { q } ) \\ & & & + \bar { X } _ {, 1 } ^ { ( 0 ) } ( \hat { q } ) & \bar { X } _ {, 1 } ^ { ( 0 ) } ( \hat { q } ) & \bar { B } _ {, a } ^ { 1 1 } ( \hat { q } ) \\ & = & 0. & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & && & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &
 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \end{array} \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & - & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &   & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  \end{array} \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & & & &  & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & & & & \ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &. & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \\ & & & & & \ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \ & & & & & \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \\ & & & & & \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &$$

In passing to the second equality, we have used ¯ X (0) , a (ˆ) = 0 proved in Eq. q (A.27). We have then used Eqs. (A.20) and (A.22) in the first and third terms, respectively, to make these vanish, whereas the second term vanishes because ¯ X (0) , 1 a (ˆ) = 0 which can be derived from q ¯ X (0) , a (ˆ) = 0. q

Then, we rewrite Eq. (A.26) as

$$B ^ { \alpha \beta } \, X _ {, \beta } ^ { ( \sigma ) } ( \hat { \xi } ) = \bar { X } _ {, 1 } ^ { ( \sigma ) } ( \hat { q } ) \, g _ {, 1 } ^ { \alpha } ( \hat { q } ) \, \bar { B } ^ { 1 1 } ( \hat { q } ), \ \sigma = 0, \, 1. \quad \text{(A.29)}$$

Here, we have used relations

$$X _ {, \alpha } ^ { ( \sigma ) } ( \hat { \xi } ) = \bar { X } _ {, \mu } ^ { ( \sigma ) } ( \hat { q } ) \, f _ {, \alpha } ^ { \mu } ( \hat { \xi } ) = \bar { X } _ {, 1 } ^ { ( \sigma ) } ( \hat { q } ) \, f _ {, \alpha } ^ { 1 } ( \hat { \xi } ),$$

$$B ^ { \alpha \beta } \, f ^ { 1 } _ {, \beta } ( \hat { \xi } ) = g ^ { \alpha } _ {, \mu } ( \hat { q } ) \, \bar { B } ^ { \mu 1 } ( \hat { q } ) = g ^ { \alpha } _ {, 1 } ( \hat { q } ) \, \bar { B } ^ { 1 1 } ( \hat { q } ),$$

which have been derived from the chain rule and Eqs. (A.5), (A.6), (A.10), and (A.24)-(A.26). We note that Eq. (A.29) is an alternative form of Eqs. (A.20)-(A.22) because Eq. (A.29) can be derived only if all of Eqs. (A.20)(A.22) are used.

Substituting the explicit forms

$$X _ {, \alpha } ^ { ( 0 ) } ( \hat { \xi } ) \ = \ V _ {, \alpha } ( \hat { \xi } ),$$

$$X _ {, \alpha } ^ { ( 1 ) } ( \hat { \xi } ) \ = \ 2 \, V _ {, \alpha \beta } ( \hat { \xi } ) \, B ^ { \beta \gamma } \, V _ {, \gamma } ( \hat { \xi } ),$$

into Eq. (A.29), we have

$$B ^ { \alpha \beta } \, V _ {, \beta } ( \hat { \xi } ) = \gamma ( \hat { q } ) \, g _ {, 1 } ^ { \alpha } ( \hat { q } ),$$

with

$$\left [ B ^ { \alpha \gamma } \, V _ {, \gamma \beta } ( \hat { \xi } ) \, - \lambda ( \hat { q } ) \delta ^ { \alpha } _ { \ \beta } \right ] g ^ { \beta } _ {, 1 } ( \hat { q } ) = 0,$$

$$\gamma ( \hat { q } ) \ \equiv \ \bar { X } _ {, 1 } ^ { ( 0 ) } ( \hat { q } ) \, \bar { B } ^ { 1 1 } ( \hat { q } ), \\ \bar { \nu } ^ { ( 1 ) \, / \, \hat { \varepsilon } \, }$$

$$\lambda ( \hat { q } ) \ \equiv \ \frac { \stackrel {, \cdot \cdot \cdot \cdot } { \bar { X } _ {, 1 } ^ { ( 1 ) } ( \hat { q } ) } } { 2 \, \bar { X } _ {, 1 } ^ { ( 0 ) } ( \hat { q } ) }. \quad \quad \quad$$

Owing to a scale indefiniteness in the definition of q 1 , there exists an ambiguity for the overall amplitude of g α , 1 (ˆ). q Without loss of generality, we can adopt the normalization

$$g _ {, 1 } ^ { \alpha } ( \hat { q } ) \, B _ { \alpha \beta } ^ { - 1 } \, g _ {, 1 } ^ { \beta } ( \hat { q } ) = 1.$$

We note that Eqs. (A.23), (A.34), (A.35), and (A.38) are identical to the decoupling conditions. By setting

$$\xi ^ { \alpha = i } \ = \ x _ { i },$$

$$\begin{array} {$$

$$\ q ^ { 1 } \ = \ q,$$

we find that Eqs. (A.23), (A.34), (A.35), and (A.38) correspond to Eqs. (3), (5), (4), and (6), respectively.

## Appendix B. Barrier and gate potential with less symmetry

A Hamiltonian we consider is given by

$$H = \frac { 1 } { 2 } \, ( p _ { 1 } ^ { 2 }$$

Here, V ′ BG ( x , x 1 2 ) denotes the barrier and gate potential with less symmetry defined by

$$V _ { \text{BG} } ^ { \prime } ( x _ { 1 }, \$$

where V BG ( x , x 1 2 ) has been given in Eq. (41). We note that the term associated with /epsilon1 breaks the symmetry with respect to x 1 → -x 1 . This potential energy surface has two local minima and one saddle point located at ( /epsilon1/ 2 b ± √ /epsilon1 2 / 4 b 2 + a/b, c/epsilon1/ 2 b ± c √ /epsilon1 2 / 4 b 2 + a/b ) ≡ ( x eq 1 ± , x eq 2 ± ) and (0 , 0), respectively, and the activation energies in the transitions from ( x eq 1+ , x eq 2+ ) to ( x eq 1 -, x eq 2 -) and from ( x eq 1 -, x eq 2 -) to ( x eq 1+ , x eq 2+ ) are given by a / 2 4 b + a/epsilon1 2 / 4 b 2 + /epsilon1 4 / 24 b 3 ± ( a/epsilon1/ 3 b + /epsilon1 3 / 12 b 2 ) √ /epsilon1 2 / 4 b 2 + a/b , respectively.

In Fig. B.6, we show the collective paths for V ′ ( a ) BG , V ′ ( b ) BG , and V ′ ( c ) BG . It reveals that the collective paths provide pathways connecting the two local minima. We stress that these results are qualitatively the same as those of the symmetric potentials V ( a ) BG , V ( b ) BG , and V ( c ) BG examined in Sec. 4.2.

We introduce three potential energy surfaces: a potential energy surface containing neither a high energy barrier nor a narrow gate V ′ ( a ) BG , a surface with a high energy barrier V ′ ( b ) BG , and a surface with a narrow gate V ′ ( c ) BG . The parameters ( a, b, c, α, β, γ, δ, /epsilon1 ) which create V ′ ( a ) BG , V ′ ( b ) BG , and V ′ ( c ) BG read (80 , 80 , 3 , 0 0045 . , 10 , 10 , 1 , 20), (240 , 240 , 3 , 0 , 0 , 10 , 20 , 60), and (80 , 80 , 3 , 0 31 . , 690 , 10 , 9 7 . , 20), respectively. All the potential energies have the same locations of the local minima as ((1 ± √ 65) / 8 , 3(1 ± √ 65) / 8).

## References

- [1] P. G. Debrunner and H. Frauenfelder, Ann. Rev. Phys. Chem. 33 (1982), 283.
- [2] Chou KC, Biophys. Chem. 20 (1984), 61.
- [3] O. Miyashita, J. N. Onuchic, and P. G. Wolynes, Proc. Natl. Acad. Sci. USA, 100 (2003), 12570.

Figure B.6: The collective paths obtained by the method developed in this work on the V ′ ( a ) BG , V ′ ( b ) BG , and V ′ ( c ) BG potential energy surfaces. The solid lines in (a) ((b) or (c)) show the collective path derived from V ′ ( a ) BG ( V ′ ( b ) BG or V ′ ( c ) BG ), respectively. In (a) ((b) or (c)), the thin lines denote the contours of V ′ ( a ) BG ( V ′ ( b ) BG or V ′ ( c ) BG ), the circles the local minima, and the squares the saddle point.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0364v3-figure-ab1e109b-7b4d-4fad-bbd4-a30e119440cc.jpg)

- [4] J. A. McCammon, Rep. Prog. Phys. 47 (1984), 1.
- [5] N. Nakajima, J. Phys. Chem. B 101 (1997), 817.
- [6] Y. Sugita and Y. Okamoto, Chemical Phys. Lett. 314 (1999), 141.
- [7] A. Laio and M. Parrinello, Proc. Natl. Acad. Sci. USA 99 (2002), 12562.
- [8] Y. Fukunishi, Y. Mikami, and H. Nakamura, J. Phys. Chem. B 107 (2003), 13201.
- [9] S. Tomonaga, Prog. Theor. Phys. 13 (1955), 467.
- [10] S. Tomonaga, Prog. Theor. Phys. 13 (1955), 482.
- [11] T. Marumori, T. Maskawa, F. Sakata, and A. Kuriyama, Prog. Theor. Phys. 64 (1980), 1294.
- [12] M. Matsuo and K. Matsuyanagi, Prog. Theor. Phys. 74 (1985), 288.
- [13] M. Matsuo, Prog. Theor. Phys. 76 (1986), 372.
- [14] M. Matsuo, T. Nakatsukasa and K. Matsuyanagi, Prog. Theor. Phys. 103 (2000), 959.
- [15] G. Do Dang, A. Bulgac, and A. Klein, Phys. Rev. C 36 (1987), 2661.
- [16] N. R. Walet, A. Klein, and G. Do Dang, J. Chem. Phys. 91 (1989), 2848.
- [17] A. Klein, N. R. Walet, and G. Do Dang, Ann. Phys. 208 (1991), 90.
- [18] G. D. Dang, A. Klein and N. R. Walet, Phys. Rep. 335 (2000), 93.
- [19] H. Aoyama and H. Kikuchi, Nucl. Phys. B 369 (1992), 219.
- [20] K. M¨ uller and L. D. Brown, Theor. Chim. Acta. 53 (1979), 75.
- [21] K. Fukui, S. Kato, and H. Fujimoto, J. Am. Chem. Soc. 97 (1975), 1.
- [22] G. Henkelman and H. J´nsson, J. Chem. Phys. o 113 (2000), 9978.
- [23] J. Wang, R. M. Wolf, J. W. Caldwell, P. A. Kollman, and D. A. Case, J. Comput. Chem. 25 (2004), 1157.

- [24] K. Tsumura et al. , in preparation.
- [25] K. Tsumura et al. , in preparation.
- [26] W. Kohn and L. J. Sham, Phys. Rev. 140 (1965), A1133.
- [27] K. Tsumura et al. , in preparation.
- [28] E. Runge and E. K. U. Gross, Phys. Rev. Lett. 52 (1984), 997.
Kyosuke Tsumura Yoshitaka Maeda Hiroyuki Watanabe
Rare transition event with self-consistent theory of large-amplitude collective motion
A numerical simulation method, based on Dang et al.'s self-consistent theory
of large-amplitude collective motion, for rare transition events is presented.
The method provides a one-dimensional pathway without knowledge of the final
configuration, which includes a dynamical effect caused by not only a potential
but also kinetic term. Although it is difficult to apply the molecular dynamics
simulation to a narrow-gate potential, the method presented is applicable to
the case. A toy model with a high-energy barrier and/or the narrow gate shows
that while the Dang et al. treatment is unstable for a changing of model
parameters, our method stable for it.