Double-well potentials are widely used to model phenomena in physics and chemistry. The system is assumed to be formed in a metastable state, its dynamical evolution providing the clues for the interpretation of the experimental data. Quantum mechanics predicts coherent oscillations of probability between wells if the double-well potential is nearly symmetric and irreversible exponential decay if the final well has an infinite width. For very asymmetric double-well potentials, these two extreme behaviors are expected to coexist. The purpose of the present paper is to investigate this coexistence and its evolution as a function of the width of (or density of states in) the second well. In this sense, increasing the density of states can be regarded as a mechanism for coherence breakdown. The dynamical evolution of the metastable state can be simulated by solving numerically the time-dependent Schroedinger equation (TDSE). This approach is general, intuitive, and gives access to time scales and dynamical effects. It also allows the inclusion of phenomenological dissipation.
Double-well potentials are widely used to model phenomena in physics and chemistry. The system is assumed to be formed in a metastable state, its dynamical evolution providing the clues for the interpretation of the experimental data. Quantum mechanics predicts coherent oscillations of probability between wells if the double-well potential is nearly symmetric and irreversible exponential decay if the final well has an infinite width. For very asymmetric double-well potentials, these two extreme behaviors are expected to coexist. The purpose of the present paper is to investigate this coexistence and its evolution as a function of the width of (or density of states in) the second well. In this sense, increasing the density of states can be regarded as a mechanism for coherence breakdown. The dynamical evolution of the metastable state can be simulated by solving numerically the time-dependent Schroedinger equation (TDSE). This approach is general, intuitive, and gives access to time scales and dynamical effects. It also allows the inclusion of phenomenological dissipation.
For a model 1d asymmetric double-well potential we calculated so-called survival probability (i.e. the probability for a particle initially localised in one well to remain there). We use a semiclassical (WKB) solution of Schroedinger equation. It is shown that behaviour essentially depends on transition probability, and on dimensionless parameter which is a ratio of characteristic frequencies for low energy non-linear in-well oscillations and inter wells tunnelling. For the potential describing a finite motion (double-well) one has always a regular behaviour. For the small value of the parameter there is well defined resonance pairs of levels and the survival probability has coherent oscillations related to resonance splitting. However for the large value of the parameter no oscillations at all for the survival probability, and there is almost an exponential decay with the characteristic time determined by Fermi golden rule. In this case one may not restrict oneself to only resonance pair levels. The number of perturbed by tunnelling levels grows proportionally to the value of this parameter (by other words instead of isolated pairs there appear the resonance regions containing the sets of strongly coupled levels). In the region of intermediate values of the parameter one has a crossover between both limiting cases, namely the exponential decay with subsequent long period recurrent behaviour.
Abstract The out-of-equilibrium quantum dynamics of an interacting Bose gas trapped in a one-dimensional asymmetric double-well potential is studied by solving the many-body Schrödinger equation numerically accurately. We examine how the gradual loss of symmetry of the confining trap affects the macroscopic quantum tunneling dynamics of the system between the two wells. In an asymmetric double well, the two wells are not equivalent anymore, say, the left well is deeper than the right one. Accordingly, we analyze the dynamics by initially preparing the condensate in both the left and the right wells. The dynamics of the system is characterized by the time evolution of a few physical quantities of increasing many-body complexity, namely, the survival probability, depletion and fragmentation, and the many-particle position and momentum variances. In particular, we have examined the frequencies and amplitudes of the oscillations of the survival probabilities, the time scale for the development of fragmentation and its degree, and the growth and oscillatory behavior of the many-particle position and momentum variances. There is an overall suppression of the oscillations of the survival probabilities in an asymmetric double well, except for resonant values of asymmetry for which the one-body ground state energy in the right well matches with one of the one-body excited states in the left well, thereby resulting in resonantly enhanced tunneling from the right well ground state. Overall, depending on whether the condensate is initially prepared in the left or right well, the repulsive inter-atomic interactions affect the survival probabilities differently. For a sufficiently strong repulsive interaction, the system is found to become fragmented. The degree of fragmentation depends both on the asymmetry of the trap and the initial well in which the condensate is prepared in a non-trivial manner. Furthermore, we show that the phenomenon of resonantly enhanced tunneling can be accompanied by a large degree of fragmentation (depletion) for the strong (weak) interaction. The many-particle position and momentum variances follow the density oscillations of the system in the asymmetric double well and bears prominent signatures of the degree of depletion or fragmentation, depending on the strength of the interactions. These quantities further exhibit a fine structure signifying a breathing-mode oscillation. Finally, a universality of fragmentation for systems made of different numbers of particles but the same interaction parameter is also found and its dependence on the asymmetry is investigated. The phenomenon is robust despite the asymmetry of the junction and admits a macroscopically-large fragmented condensate characterized by a diverging many-particle position variance. This is as far as one can get from the dynamics of the density in the junction.
The time evolution of quantum mechanical states in a square well of infinite depth with a Dirac δ function at its center has been examined for cases where the initial state was localized in one of the wells. Let well A denote the well in which the initial state is localized, and let P(t) denote the integrated probability density Ψ*Ψ in well A. For opaque barriers the time-dependent system is adequately described by a two-state model in which only the pair of stationary states of even and odd parity are considered whose wave functions in well A, apart from an arbitrary phase factor, are largely identical with the wave function of the initial state. For P(t) a harmonic oscillation is observed whose frequency νt is well approximated by the well-known formula for the tunneling frequency νt≂ΔE/h, where ΔE represents the energy separation among the pair of states in the model. For the present model of a symmetric double well it has been shown that for highly transparent barriers a three-state model can describe the time evolution adequately. The three stationary states involved in this model are a state of odd parity whose wave function in well A is largely the same as the wave function of the initial state and a pair of stationary states of even parity which on the energy scale are immediately above and below the first state. In this three-state model the function P(t) is a superposition of two sinusoidal functions with nearly identical amplitudes and frequencies plus a constant. As a consequence, the amplitude of P(t) changes harmonically. In the present model a δ function has been used as a barrier in order to minimize the mathematical detail involved in the time-dependent treatment. It is to be expected that the beating in P(t) can be observed also in the time evolution of a state localized in one of the wells of a symmetric double-minimum potential with a more realistic ‘‘low’’ barrier if the density of energy levels near the energies of the levels to be considered in the corresponding three-state model varies slowly with energy.
Metastability separates apparent equilibration from true convergence: a dynamical system can remain near an intermediate state long before reaching its stationary state. In quantum dynamics used for optimization, such a plateau can delay progress toward a target concentrated near the global minimizers of a nonconvex potential. In this paper, we initiate the study of the general mechanisms by which coherent evolution and dissipation shape metastability, using quantum Langevin dynamics (QLD) as a concrete example [CLW+25], where the objective function is encoded in the potential and coherent transfer between wells interacts with tunable dissipation. To quantify this interaction, we solve a two-state model precisely and derive a spectral-separation criterion in terms of tunneling frequency, dissipation strength, and damping imbalance. Extending the analysis to multi-state models, we show that sufficiently weak tunneling between wells can separate slow population transfer from fast relaxation within each well. In addition, we establish conditions under which changing the structure of dissipation shortens the time needed to approach a prescribed target state. Numerical experiments on double-well and multi-well potentials demonstrate these improvements while keeping the target and accuracy fixed and accounting for the time spent on control.
The motion of the center of a wave packet localized in one of the wells of an asymmetric double-well potential is studied. The asymmetric potential is obtained from a solvable symmetric potential using the Gel'fand–Levitan formalism, by keeping the energy eigenvalues unchanged, but changing the normalization of the ground-state wave function. While in a symmetric double well the wave packet tunnels back and forth between the two wells, for an asymmetric well, in general, tunneling does not take place. This is independent of the initial location of the wave packet, i.e., whether it is centered at the minimum of the shallower or of the deeper well. However, when such a system is coupled to a dissipative force, then the tunneling becomes possible. For instance it is shown that if, as a model of dissipative coupling, one chooses Gisin's nonlinear evolution equation, then the center of the wave packet ends up in the deeper well (decay of false vacuum). This result depends on the particular model of dissipation, for instance, an optical potential model yields a different result.
For a model one-dimensional asymmetric double-well potential we calculated the so-called survival probability (i.e., the probability for a particle initially localized in one well to remain there). We use a semiclassical (WKB) solution of the Schrödinger equation. It is shown that behavior essentially depends on transition probability, and on a dimensionless parameter Lambda that is a ratio of characteristic frequencies for low-energy nonlinear in-well oscillations and interwell tunneling. For the potential describing a finite motion (double-well) one has always a regular behavior. For Lambda<<1, there are well defined resonance pairs of levels and the survival probability has coherent oscillations related to resonance splitting. However, for Lambda>>1 there are no oscillations at all for the survival probability, and there is almost an exponential decay with the characteristic time determined by Fermi golden rule. In this case, one may not restrict himself to only resonance pair levels. The number of levels perturbed by tunneling grows proportionally to square root of [Lambda] (in other words, instead of isolated pairs there appear the resonance regions containing the sets of strongly coupled levels). In the region of intermediate values of Lambda one has a crossover between both limiting cases, namely, the exponential decay with subsequent long period recurrent behavior.