Quantum-tunneling transitions and exact statistical mechanics of bistable systems with parametrized Dikandé–Kofané double-well potentials

We consider a one-dimensional system of interacting particles (which can be atoms, molecules, ions, etc.), in which particles are subjected to a bistable potential the double-well shape of which is tunable via a shape deformability parameter. Our objective is to examine the impact of shape deformability on the order of transition in quantum tunneling in the bistable system, and on the possible existence of exact solutions to the transfer-integral operator associated with the partition function of the system. The bistable potential is represented by a class composed of three families of parametrized double-well potentials, whose minima and barrier height can be tuned distinctly. It is found that the extra degree of freedom, introduced by the shape deformability parameter, favors a first-order transition in quantum tunneling, in addition to the second-order transition predicted with the ϕ4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi ^4$$\end{document} model. This first-order transition in quantum tunneling, which is consistent with Chudnovsky’s conjecture of the influence of the shape of the potential barrier on the order of thermally assisted transitions in bistable systems, is shown to occur at a critical value of the shape-deformability parameter which is the same for the three families of parametrized double-well potentials. Concerning the statistical mechanics of the system, the associate partition function is mapped onto a spectral problem by means of the transfer-integral formalism. The condition that the partition function can be exactly integrable, is determined by a criterion enabling exact eigenvalues and eigenfunctions for the transfer-integral operator. Analytical expressions of some of these exact eigenvalues and eigenfunctions are given, and the corresponding ground-state wavefunctions are used to compute the probability density which is relevant for calculations of thermodynamic quantities such as the correlation functions and the correlation lengths.


Introduction
In low-dimensional systems, structural phase transitions are triggered by symmetry breakings that induce phase instabilities under specific conditions [1][2][3][4][5][6][7][8]. In general these instabilities are governed by changes in equilibrium properties of the systems, and in the displacive regime [3][4][5][6] they can be accompanied by defects in their topological structures including disclinations, discommensurations, domain walls, large-amplitude charge-density waves (phasons) and instantons to name these few ones. Characteristic properties of these structural defects, which are usually referred to as solitons or solitary waves [9,10], have attracted a great deal of attention over the past years given their widespread applications in a broad range of physical contexts [9,10]. a e-mail: fernand.naha.nz@gmail.com b e-mail: dikande.alain@ubuea.cm (corresponding author) c e-mail: esmkam@yahoo.com In the specific context of condensed matter Physics, the role of instantons in tunneling processes observed in some of these systems at low temperatures is now well established [11][12][13][14]. As a matter of fact, in multistate systems, the transitions between two states separated by an energy barrier usually occur either in the classical regime, where they are driven by thermal activations, or in the quantum regime driven by tunneling processes. At high temperatures, classical thermal activations govern the transition which occurs as a hopping over the barrier. However, as the temperature approaches zero, quantum fluctuations gain in importance and the transition is driven by quantum tunnelings through the barrier. However in the later case, the system dynamics can be described by classical configurations which are either instantons or vacuum, that dominate the thermal rate [11][12][13][14][15]. As the temperature increases, thermally-induced crossovers become more and more relevant and at some critical temperature a phase transition in quantum tunneling can take place [13][14][15]. The phenomenon of phase transition in quantum tunneling has been extensively investigated over the past years [13][14][15]; thus, it is known that some physical systems can exhibit not only a smooth second-order transition in quantum tunneling at a critical temperature T (2) 0 , but a transition of the first order can also occur [13,[16][17][18][19][20][21]. Instructively, this first-order transition in quantum tunneling is closely related to the variation of the instanton action with respect to the instanton energy [13,[16][17][18][19][20][21].
Another relevant issue in the study of phase transitions in the presence of topological defects (i.e., kink solitons), is the statistical mechanics in the relevant temperature regions where phase transitions are expected. Addressing this issue [22][23][24][25][26] a phenomenological theory was introduced, that rests on the functional-integral operator formalism and assumes that at low temperature, the statistical mechanics of one-dimensional systems with kinks and phonons can be mapped onto an eigenvalue problem associated with the transfer-integral operator. For classic models such as the sine-Gordon and φ 4 models, the eigenvalue problem of the transfer-integral operator does not admit exact solutions [23,24], nevertheless their lowest eigenstates, obtained via approximate methods such as the WKB approximation, can permit to construct acceptable statistical mechanical quantities [23]. Since in the theory of critical phenomena phase instabilities are predicted to occur at temperatures where the free energy, and/or its derivatives representing thermodynamic functions, are nonanalytic, an appropriate formulation of the partition function for a given system therefore represents a key prerequisite for an accurate description of phase instabilities that might take place in the system.
While theoretical approaches used in the studies of phase transitions in quantum tunneling, and of the statistical mechanics of one-dimensional systems, are usually well accepted, the mathematical models describing these physical systems are sometimes weak. For instance the sine-Gordon model was introduced to describe the formation of solitons in one-dimensional systems with periodic energy landscapes, and to investigate soliton contributions to the statistical mechanics of these specific systems [9,23]. However, due to its rigid profile characterized by fixed barrier height and extrema positions as well as a constant period of the substrate potential, the sine-Gordon model could not be exploited for contexts where shape profiles of the periodic substrate could vary in response to processes intrinsic to the systems [9]. To overcome the weakness related to its rigid shape profile, a parametrized version of the sine-Gordon model was proposed, i.e., the socalled Remoissenet-Peyrard potential [27][28][29]. On the other hand the φ 4 model was introduced by Landau in his mean-field theory of second-order phase transitions [25,26]. Like the sine-Gordon model, the rigid shape of this double-well model (fixed barrier height and minima positions) restricts its scope to a very limited number of real bistable systems. This limitation motivated the proposal of bistable models with deformable doublewell (DW) shapes.
Besides the famous double-Morse model [30], introduced in the study of hydrogen-bonded systems such as organic ferroelectrics [31] or DNA conformational transitions [32], some bistable models with parametric DW potentials are found in the literature [30,33,[33][34][35][36][37]. In a series of papers [38][39][40], three different families of parametrized DW potentials were proposed by Dikandé and Kofané (hereafter referred to as DK potentials) that had in common the φ 4 model as a specific limit. The three models differ in their respective distinct parametrization, and thus offer three possible different deformable features for the DW energy profile: in one family the barrier height is constant but positions of the potential minima can be varied [38], in the second family the parametrization favors continuous variation of the barrier height without changing the two degenerate minima [39], and in the third family both the barrier height and minima of the DW potential can be tuned simultaneously [40]. Given that the three families tend to the φ 4 model when their common deformability parameter tends to zero, investigating phase transitions in quantum tunneling and the classical statistical mechanics of the three families of models, is a relevant exercise susceptible to unveil phenomena that the φ 4 model cannot account for. For instance, the φ 4 model predicts the transition in quantum tunneling to be only of second order [15], and we mentioned the fact that the statistical mechanics of this model is non exact within the framework of the transfer-integral formalism [22].
Concerning the possibility of a first-order transition in quantum tunneling in bistable systems, this issue was addressed in a relevant work by Chudnovsky [18], who postulated that in general this would depend on the shape of the potential barrier. Exploiting Chudnovsky's conjecture, Zhou et al [15] examined the problem by considering a parametrized double-well potential with a variable barrier height but fixed positions of potential minima [39]. They concluded that the increase of the barrier height with the deformability parameter μ, is actually the main factor favoring a first-order transition in quantum tunneling in this particular model of bistable system. However in a recent work [41], we carried out a study on the same problem using a version of the parametrized double-well potential with fixed barrier height and fixed potential minima, but with variable shape of the barrier top. We obtained that despite the fixed barrier height and fixed minima positions, a first-order transition in quantum tunneling could still be observed in this other model of bistable system. Motivated by the two recent studies [39,41], we have undertaken to investigate conditions for the occurrence of a first-order transition in quantum tunneling for the three distinct members of the so-called Dikandé-Kofané (DK) potentials [38][39][40]. The present work will address this issue. We shall also investigate values of the shape deformability parameter, for which exact wavefunctions and eigenvalues of the transfer-integral operator associated with the three distinct parametrized DK potentials, can be found as done recently with the new mem-ber of this family introduced in ref. [41], characterized by a fixed barrier height and fixed minima positions, but a tunable flatness of the barrier top.

Three families of parametric double-well potentials
One-dimensional bistable systems can be represented as a linear chain of particles (atoms, molecules or ions) interacting via two-body forces, with the individual particles lying in one-body potentials with a DW profile. In some physical contexts, shapes of the DW potential can be tuned. For instance particles in the chain can interact with their neighbors via bonds whose lengths are constantly varying (e.g., bond stretching or compression processes), due to variations of pressure, isotopic substitutions or chemical reactions that inevitably affect equilibrium properties of particles along the chain.
In the present study we are interested in three specific families of DW potentials with deformable shape profiles [38][39][40], which can be expressed most generally as [40]: where a(μ) > 0 and α(μ) are real functions of the shape deformability parameter μ. V (u, μ) is a class of parametric DW potentials whose shapes can be tuned differently, but which admit the φ 4 potential [23] as a common asymptotic limit. The first family of parametrized DW potentials associated with V (u, μ), has the following parameters: a 0 = constant. This first model represents a DW potential for which positions of the two degenerate minima can be varied, leaving unchanged the barrier height [38]. Some shape profiles of this first family of DW potentials are sketched in Fig. 1a. An increase of μ maintains the potential barrier fixed at a 0 , whereas positions of the two potential minima; are continuously shifted near to each other as μ is varied. As the figure shows this shift induces a narrowing of potential wells, and as a sequel a reduction of the barrier width. A narrowing of potential wells will enhance confinement of particles inside the wells, thus favoring tunneling processes given that the stiffness of the two walls culminating at the barrier top is expected to increase with the confinement. For the second family, characteristic parameters are given by: This second model corresponds to a DW potential with fixed degenerate minima, but a tunable barrier height [39]. Figure 1b illustrates profiles of the second model of DW potential, for some values of the deformability parameter μ. The two minima are indeed always fixed at u = ± 1, but when μ is increased the barrier height continuously decreases with the barrier top gradually flattened.
The third model has the following parameter values: This third model mimics a DW potential whose degenerate minima and barrier height can be simultaneously varied [40]. In Fig. 1c it is seen that an increase of μ causes the increase of the barrier height, as well as a shift away of positions of the two potential minima. This third model is clearly most prone to tunneling processes when μ is large, given that by increasing the deformability parameter the barrier height is increased while the distance from the equilibrium positions of a particle inside the wells to the position of the potential barrier, is relatively large. Therefore it will require more and more energy for a particle to move from inside the wells to the barrier top, to hope from one well to the other unless it proceeds via tunneling.
Remark that When μ tends to zero, the three families of parametrized DW models reduce to the standard φ 4 model, i.e., [23]: In the next section, we investigate phase transitions in quantum tunneling mediated by periodic instantons, for the three families of parametrized DW models. We shall sometimes refer to the three families of parametrized DW potentials as DK potentials [38][39][40].

Periodons and transitions in quantum tunneling
Consider the following Euclidean action for onedimensional quantum models (in dimensionless unity) [13][14][15]: where u represents a scalar field in one time and zero space dimension. The variable τ = it is the imaginary time, and V (u, μ) is a nonlinear potential energy of the general form Eq. (1). The integral in Eq. (7) is taken over the period τ p of the trajectory in Euclidean space. In statistical mechanics, the period τ p is related to temperature For simplicity and without loss of generality, in this section we shall assume ( , k B ) ≡ 1.
To an exponential accuracy, the decay rate in the semiclassical limit is given by: where S min is the minimum effective Euclidean action obtained by minimizing Eq. (7), along a trajectory u c (τ ) of the classical pseudoparticle. This trajectory satisfies the equation: where the integration constant E ≥ 0 can be regarded as the energy of the pseudoparticle of mass unity. When E = 0, The trajectory u c (τ ) is a regular vacuum instanton described by the single-kink soliton: where: is the kink width. However when E > 0, using periodic boundary conditions and identifying τ p as the period of motion, we find that the instanton solution to Eq. (9) should now read: In this later expression, the quantity sn(τ, κ) is one of Jacobi elliptic functions [43], its modulus κ is given by: .
(13) The two parameters C 1 and C 2 in the instanton solution Eq. (12) are defined as: The instanton trajectory (12) possesses real periods 4mK(κ), with m an integer and K(κ) the quarter period defined in terms of the complete elliptic integral of the first kind [43]. Hence, the trajectory Eq. (12) is a periodic instanton, with a period τ p corresponding to m = 1 and given by: We refer to the periodic instanton Eq. (12) as a periodon. From Eq. (12) we can see that the periodon exists only in the range 0 ≤ E ≤ a(μ), and that the expression of the vacuum instanton, i.e., the single-kink soliton Eq. (10), will be recovered when E = 0 (i.e., κ = 1 [43]). Using the energy-integral equation (9), the classical action S p of the periodon can be computed and is found to be: In this formula E(κ) and Π(C 1 , κ) are complete elliptic integrals of the second and third kinds, respectively. We note that when E = a(μ) the trajectory reduces to a trivial configuration u c (τ ) = 0, which is periodic with an arbitrary period. This configuration, which we shall call sphaleron [44], has an action coinciding with the thermodynamic action S 0 given by [44]: For the sphaleron configuration, the escape rate is the Boltzmann formula representing a pure thermal activation, i.e.: From the above discussions, we can conclude that a periodon (i.e., the periodic instanton) interpolates between the sphaleron and the vacuum instanton (i.e., the single-kink instanton). Therefore the actual escape rate should be of the form of Eq. (8), where the minimum action S min (E) is given by: When dealing with periodic problems in classical statistical mechanics, it is known that the derivative of the action with respect to the energy is equivalent to the oscillation time τ at this energy, where τ is proportional to the inverse temperature T . Hence, with the action corresponding to the motion in the potentials we find: Taking Eq. (19) together with Eqs. (14) and (15), we can determine the influence of the shape deformability parameter on the temperature dependence of S min , using the dependence of the period of periodon on the energy E. In the following analysis these equations will enable us to point out the existence of a critical value of the shape deformability parameter, for which a firstorder transition from quantum to thermal regimes can occur for the three families of DKDW models. Chudnovsky [18] postulated two criteria for the occurrence of transitions in quantum tunneling. The first states that the transition is of first order if τ p (E) decreases to a minimum and then rises again, when E increases from the potential bottom (i.e., minima)   to the barrier top. On the other hand, if the instanton period τ p (E) is monotonically decreasing when E increases, the transition is of second order. The second criterion states that if at some critical temperature the first derivative of S min (T ) is discontinuous, and an abrupt change is observed in the temperature dependence of the action, then the transition in quantum tunneling is a first-order transition in temperature. Chudnovsky's first criterion is equivalent to solving the equation dτ p (E)/dE = 0, which is tractable when seeking nontrivial solutions and considering the temperature dependence of the period given by (14). In a numerical procedure we have chosen a 0 = 1/2, and computed the energy E 1 of the minimum τ p (E) by varying the deformability parameter μ from the value μ = 3 downwards. Results for the three models are listed in Tables 1, 2 and 3.
The emerging trend from the three tables is E 1 tending to the maximum energy E 0 of the pseudoparticle energy, when μ approaches 3/2. Hence μ c = 3/2 is a critical value of the shape deformability parameter, values of μ smaller than μ c will correspond to for- bidden energy configurations (i.e., nonanalytical values of E). In Chudnovsky picture, a first-order quantumto-classical transition will occur when μ > μ c . When μ ≈ μ c we obtain E 1 ≈ E 0 . Moreover E 1 → 0 as μ is increased above μ c , so [42] increasing the shape deformability parameter above the critical value μ c results in a sharper first-order transition in quantum tunneling. Figures 2, 3 and 4 show the energy dependence of the period τ p of the periodon for the three models, considering two values of the shape deformability parameter μ, respectively, below and above the critical value μ c .
For the three models the figures suggest a monotonic decrease of the period with increasing energy for μ = 1 (a value lower than μ c ). However, for μ = 3, we notice an increase of the period after a critical value of the energy E. This increase reflects favorable conditions for a first-order transition in quantum tunneling.
Characteristic features emerging from Figs. 2, 3 and 4 are also reflected in the corresponding plots of the thermodynamic action and the periodon action, shown in Figs. 5, 6 and 7 still for μ = 1 and μ = 3. Here too the three models display similar behaviors, indeed a smooth change from S p to S 0 is observed when the temperature increases for μ = 1 as we can see in Figs. 5a, 6a and 7a. Clearly such change is, from the standpoint of Chudnovsky criteria, the signature of a second-order transition. Figures 5b, 6b and 7b are characterized by an abrupt change in the temperature dependence of the minimum action, so they satisfy Chudnovsky's second criterion [18] for a first-order transition from quantum to thermal regimes. To provide an analytical proof of a possible finite critical value of μ for a first-order transition in quantum tunneling, we follow Ref. [45]. We postulate that if the period τ p (E → a(μ)) of the periodon close to the barrier top can be found, then the relation τ p (E → a(μ))−τ s < 0 (or ω 2 −ω 2 s > 0) is a sufficient condition for the system to exhibit a first-order transition. Note that τ s denotes the period of small oscillations around the sphaleron, ω s is the corresponding oscillation frequencies and ω is the frequency of oscillations around the energy barrier of height a(μ). For the three models, this relation yields [45]: where u sph = 0 corresponds to the position of the sphaleron solution. Since all the three families of DKDW potentials are symmetric we should have V (0, μ) = 0, and the above criterion (i.e., Eq. (20)) reduces to V (0, μ) < 0. Using the general expression of V (u, μ) given by formula (1) we finally obtain: The last inequality reduces to μ 2 > 3/2, valid for the three families of DKDW potentials.

Exact statistical mechanics
The aim of the transfer-integral formalism [22][23][24] is to express the partition function associated with the statistical mechanics of low-dimensional systems, into a spectral problem whose solutions lead to thermodynamic quantities such as the free energy, the internal energy, the entropy, correlation functions, the correlation length and so on [23]. This eigenvalue problem is represented by a linear Schrödinger-like equation, which in some contexts do not admit exact solutions but nevertheless can be treated within the framework of approximation methods such as the WKB method [23]. The φ 4 is a typical non-integrable model, for which the WKB approximation has made it possible to obtain almost exact low-lying eigenmodes [23]. Using the φ 4 as generic model a generalized potential of the form V (x) = [(η + ξ)x 2 − (η − ξ)] 2 , for which the Schrödinger equation can admit exact solutions for specific values of the characteristic parameters η and ε, was proposed [46]. Quite remarkably this generalized model can be connected to some existing double sinh-Gordon models [47,48], that turns out to be quasi-exactly solv-able (QES) under appropriate conditions. Note that by quasi-exact solvability we mean the possibility to carry out exact partial diagonalization of the partition function.
Consider a one-dimensional chain of N interacting atoms of identical mass m (which we set to unity), placed in the field of a substrate potential V (u, μ), where u stands for the displacement field of atoms to which we associate a conjugate π. The (dimensionless) Hamiltonian of such system can be written: where V (u, μ) is given by Eq. (1). As found in the previous section, the Hamiltonian (22) admits single and periodic-kink soliton solutions given by (10) and (12), which we called, respectively, vacuum instanton and periodon. Here we focus on the thermodynamics of the DKDW models in the presence of vacuum instanton. Proceeding, our main concern is the calculation of the partition function which can be factorized as: In this factor the kinetic part is [23]: with β = (k B T ) −1 . The configurational part can be written [23]: where j are eigenvalues of the transfer-integral operator, involving the strain (i.e., interatomic) and on-site potential (i.e., one-body) energy terms in the Hamiltonian (22). In the displacive regime [22,23,49] the eigenvalues j are those of the Schrödinger-like equation: Let us adopt the variable change z = α(μ) u and introduce the new parameters: With these new quantities Eq. (26) becomes: In Ref. [36], solving the eigenvalue problem for a quantum-mechanical system with a hyperbolic potential closely similar to Eq. (28), Konwent established that provided specific constraints are imposed between characteristic parameters of the spectral equation, exact eigen solutions could exist. Following the same consideration let us introduce a constraint between the temperature T and the deformability parameter μ, through: where q is assumed to be a positive integer. Substituting Eq. (29) into Eq. (28), we obtain [36]: Equation (30) in this representation the function r(z) should be a polynomial and should be a linear combination of either cosh m p z or sinhm p z, where z 0 and m p ( p = 0, 1, 2, ...) are constants to be determined. Expressions for the (unnormalized) ground states and the associated eigen energies, for the first four values of q, are obtained as: q = 1: ψ 0 (u) = 6 1 + 2μ 2 exp − 3 cosh(2α(μ)u) 2(1 + 2μ 2 ) Solutions corresponding to eigenstates above the ground state, are given in the Appendix for the same four selected values of q. Equation (29) actually represents the condition for quasi-exact solvability of the system, and can be exploited to analyze the ratio a(μ)/β −1 of the energy barrier to the thermodynamic energy β −1 . In Fig. 8 we plotted the ratio a(μ)/β −1 as a function of μ, for four different values of q. One sees that for the three families of DKDW models, the solvability requires β −1 → ∞ in the limit μ → 0. This is in agreement with known results of the transfer-integral formalism for the φ 4 model, i.e., the spectral problem for this model does not admits exact solutions at finite temperatures [23].
For small values of μ, the solvability condition suggests that k B T >> a(μ) for the three models. As μ increases, temperatures obtained from the lowest to the largest values of q steadily decrease. But while going far below the energy barrier in the first model as it is noticeable in Fig. 8a, in Fig. 8b, c we observe that the temperatures in model 2 and model 3 instead tend to a particular limit. We found analytically that this limit is coincidentally the same for both models and a(μ)/β −1 → q/ √ 2 as μ → ∞. Therefore the solvability condition will hold at four temperatures below the symmetry breaking temperature in model 1, whereas only three temperatures below the energy barrier will meet the condition in both model 2 and model 3 given that the temperature at q = 1, in these two last models, is always greater than the energy barrier despite its decrease with an increase of μ.
According to the transfer-integral formalism, the free energy is intimately related to the ground-state energy 0 [23,49]. The influence of the shape deformability parameter μ on the ground-state energies is represented in Fig. 9. The three models are expected to display infinitely high ground-state energies in the limit μ → 0, in agreement with the non-integrability of the φ 4 model within the framework of the transfer-integral formalism.
As μ is increased, the ground-state energies in model 1 and model 2 drastically decrease to a critical limit as seen in Fig. 9a, b, respectively. As for model 3, Fig. 9c shows that the ground-state energies corresponding to different q decrease to a minimum, then steadily rise back to infinity. Let us also examine the influence of μ on the competition between the ground-state energy and the barrier height a(μ), taking into consideration the variation of q and hence of temperature. In Fig. 10 we plot 0 /a(μ) as a function of μ, for four different values of q.
The motivation for studying this ratio resides in the important fact that it does actually not depends on the choice of a(μ) and α(μ), so the observations done on the relative position of the energy levels will be valid for the three models. In Fig. 10 we observe that the ground-state energy is infinitely higher than the barrier height for μ → 0, and drops with an increase of μ  Fig. 10 Plot of the ratio of the ground-state energy 0 to the barrier height a(μ), versus the shape deformability parameter μ for q = 1 (solid line), q = 2 (dashed line), q = 3 (dash-dotted line), q = 4 (dotted line). Here a0 = 1, and the horizontal line marks the region where the groundstate energy is located exactly on top of the barrier irrespective of the temperature. Moreover, increasing q for a fixed value of μ will yield a lower energy level. This means that the choice of q, combined with an increase of μ, has a great impact on the position of the groundstate energy level with respect to the energy barrier. For instance, by choosing q = 1, the energy level critically drops but will always be located above the energy barrier. A similar drop is observed when choosing higher values of q. However, when μ grows higher than a critical value namely μ s = 3/2, μ s = 3/4 and μ s ≈ 0.717 for q = 2, q = 3 and q = 4, respectively, the energy level falls below the energy barrier and tends to a finite value and the ratio ( 0 /a(μ))(μ → ∞) = (3/4) for q = 2, ( 0 /a(μ))(μ → ∞) = (5/9) for q = 3 and ( 0 /a(μ))(μ → ∞) = (7/16) for q = 4.
Knowing the exact ground-state energies and the corresponding eigenfunctions, it is motivating to illustrate their use in the evaluation of some thermodynamic functions. In this goal, we consider a relevant quantity that plays the key role in the calculation of correlation functions [23,49], namely the probability density. Since the wavefunctions corresponding to eigenstates of the transfer-integral operator are all real, the probability density associated with the classical field u will be nothing else but the square of the normalized groundstate wavefunction. Sketched in Fig. 11 are ground-state wavefunctions for the three families of DKDW models, computed for two different temperatures and considering several values of the shape deformability parameter μ. To understand the physics in the curves of Fig.  11 as concerns their variations with the deformability parameter μ, it is useful to recall that in Fig. 10 we have seen that for small values of μ the ground-state energies were in the high-temperature region, as they were found to be greater than the height of the potential barrier for the three models. Consequently in this regime a pseudoparticle will possess sufficient energy to escape from one well to another over the barrier. Thus, the whole configurational space becomes probable and the full state space is covered with power-law tails, with a maximum probability density at the barrier peak. As μ increases, the probability density increases in amplitude for the three models, but while its width decreases for the first and second model the width of the probability density is instead increasing for the third model. An increase of μ in the case of model 1 and model 2 causes a decrease in the width and height of the potential barrier. In the two models the shape deformability will also strengthen the steepness of the reflective walls, and thus restricts the attainable space to the region covered by just the valleys and the energy barrier. In model 3, the potential barrier increases with its width gradually increasing as the degenerate minima move far from each other. Increasing the deformability parameter μ consequently broadens the attainable space covered by the energy barrier. The ground-state energies being high in the three models at large μ, the energy barrier thus remains the most probable. To end this section it is relevant to stress that when analyzing properties of the probability density in the present specific context, it is relevant to always keep in mind the importance of the choice of the value of q. Taking q = 1 for instance, the ground-state energy is always above the energy barrier irrespective of μ. Indeed even if an increase of μ lowers the temperature the probability density will always be single picked, with the peak position coinciding with the position of the potential barrier (top graphs in Fig. 11). For values of q greater than one, an increase of μ above μ s lowers the ground-state level below the energy barrier. The transitions from one well to another thus become less probable as the temperature decreases below the transition temperature. In this latter case the pseudoparticle stays much longer confined in the potential wells, such that the valleys become more probable. This trend is illustrated in Fig.  11 by the probability density exhibiting a double-peak shape, with each peak located in the vicinity of one of the two degenerate potential minima.

Conclusion
The decay of metastable states over an energy barrier is a physical process inherent to a large number of systems, and particularly those undergoing phase transitions [51][52][53]. It is observed in the motion of atoms trapped in the field of force of a multi-well potential [52], the motion of kinks and dislocations across the Peierls-Nabarro barrier [54,55], the excitations of Cooper pairs in a Josephson loop [56], the dynamics of macromolecules in DNA strands and other molecular chains across hydrogen bridges and so on. In a pioneer study, Kramers [51] suggested that the lifetime of a classical particle in a metastable state, separated from a stable equilibrium state by an energy barrier, would obey Arrhenius law. Later on, it was established that Kramers theory is valid main at high temperatures, where thermal activations are likely to overcome high potential barriers. At low enough temperatures thermal fluctuations are not enough to favor a jump over the potential barrier, and only by quantum tunneling the particle will cross the barrier. It is known [42] that there exists a critical temperature T c at which there is a transition from classical to quantum decay regimes. The transition can be of first order with a discontinuous first derivative of the Euclidean action, or of second-order with only a second derivative of the action developing a jump. This transition, also referred to as transition in quantum tunneling, has recently been investigated for some physical systems [13][14][15][16][17][18][19][20][21]. For systems described by the φ 4 model, the transition in quantum tunneling is strictly of second order [13][14][15]42].
In this study, we determined conditions for the occurrence of first-order transitions in quantum tunneling in bistable systems, using a class of three double-well potential models whose shapes can be tuned by varying a deformability parameter. The first-order transition supplements the second-order transition in quantum tunneling predicted with the φ 4 model, which turns out to be a specific limit of the three parametrized double-well models. The first-order transition occurs at a finite critical value of the shape deformability parameter, which is the same for all three models. We have also examined conditions for exact integrability of the partition function associated with the statistical mechanics of the three models, within the framework of the transfer-integral operator formalism. In this formalism, the partition function is mapped onto a linear Schrödinger equation whose eigenvalues and eigenfunctions contribute to the formulation of relevant thermodynamic quantities such as the free energy, the correlation length and the correlation functions. We established that when the thermodynamic temperature and the deformability parameter are connected by a specific relation, exact solutions to this spectral problem can exist. We derived the eigenvalues and eigenfunctions of some of these exact solutions analytically, and discussed the physical implications of their dependence on the shape deformability parameter.
Bistable processes are among the most frequent physical phenomena observed in nature [57], they usually emerge in form of a transition between two states of equivalent energy across an potential barrier and are present in structural phase transitions, dynamical properties of anharmonic electromechanical oscillators displaying hysteresis features [58], organized cellular structures in tissues and embryonic cells [59,60], glycolitic oscillations in suspensions of yeast cells in unison [59,60], the assembly of pacemaker cells in the sino-atrial node and so on [58]. Most commonly the bistability in these systems are represented by the φ 4 model, a universal model which unfortunately suffers from a weakness due to its rigid double-well profile. Fofr instance this model predicts that bistable systems can undergo only a second-order transition in quantum tunneling, whereas experiments have established that several physical systems with double-well energy landscapes could also undergo a first-order transition in quantum tunneling. Moreover, the transfer-integral formalism for the φ 4 generates an eigenvalue problem which is not exactly solvable. The present study emphasizes the need for parametrizing the standard φ 4 potential, to account for several physical processes observed experimentally but which the φ 4 has not been able to predict. We note that the influence of the double-well shape profile of the bistable potential, on the order of transitions in quantum tunneling in multistate systems, has already been discussed in some previous works [42]. Our study provides a more comprehensive insight onto this problem, given that the shapes of the three families of double-well potentials can be tuned distinctively, and consequently represent a very large variety of bistable physical systems.