arxiv:cond-mat/ v1 7 Nov 1997

Transcription

1 Stepwise quantum decay of self-localized solitons V.Hizhnyakov and D.Nevedrov Institute of Theoretical Physics, University of Tartu, Tähe, EE Tartu, Estonia. Institute of Physics, Riia 1, EE Tartu, Estonia. arxiv:cond-mat/971167v1 7 Nov 1997 The two-phonon decay of self-localized soliton in a one-dimensional monatomic anharmonic lattice caused by cubic anharmonicity is considered. It is shown that the decay takes place with emission of phonon bursts. The average rate of emission of phonons is of the order of vibrational quantum per vibrational period. Characteristic time of the relaxation is determined by the quantum anharmonicity parameter; this time may vary from a few (quantum lattices, large anharmonicity to thousands (ordinary lattices, small anharmonicity of vibrational periods. PACS numbers: 63.Kr, 63.Pw, 63.Ry The dynamics of strong nonlinear excitations in polymers and quasi-one-dimensional biomolecular chains is an active research field. Theoretical studies of anharmonic perfect lattices have shown the existence of localized vibrations (self-localized solitons, SLSs with frequencies above the phonon band or in the gap of the phonon spectrum (see 1 1 and references therein. SLSs are soliton-like excitations in discrete lattices, and are thus closely related to ordinary solitons. The existence of stationary and moving SLSs was derived in the frame of classical mechanics. Until now there was a little discussion about the influence of quantum (and thermal fluctuations on stability of the SLS. An exception is the paperofovchinnikov 1 whoarguesthatthedecayoftheslscausedbythesefluctuationsdiminisheswith the increase of the mode amplitude. This statement, however, is based on a perturbational consideration which can not be applied to the description of evolution of vibrations with a strong amplitude. Recently we showed that the perturbational treatment of the effect of quantum and thermal fluctuations on a local mode associated with a defect atom in a lattice fails for the case when the amplitude of vibration is large. In fact, the two-phonon damping of the local mode, caused by cubic anharmonicity, behaves dramatically with the change of the amplitude: at definite critical amplitudes relaxation jumps take place being accompanied by a generation of phonon bursts. This effect indicates that the quantum and thermal fluctuations may dominate in the dynamics of strong vibrations with the energy of the mode being within a specific range. The strong field of the local vibration causes the transformation of phonon operators and the increase of the number of phonons in time 11,1. This mechanism of phonon generation by a local vibration is analogous to the mechanism of black hole radiation In this communication we extend the theory to the SLS in a monatomic one-dimensional lattice (chain with cubic and hard quartic anharmonicity. This case is of special interest since the monatomic chain with both anharmonicities is the simplest model for the investigation of the quantum relaxation of SLS in a perfect lattice. Hard quartic anharmonicity is a prerequisite of existence of the SLS, while the cubic anharmonicity stands for the two-phonon decay of the SLS. We note that SLSs in this model are stable in classical limit, since all harmonics of the mode are out of resonance with lattice phonons 3. The odd SLS is examined and analytical nonperturbative solution of the problem and results of numerical calculations are presented. The potential energy operator in a monatomic one-dimensional lattice, which includes linear, and the first two nonlinear terms and takes account of the nearest neighbour interaction, has the following form: ˆV = n r= K r r (Ûn+1 Ûn r, (1 where Ûr is the operator of the longitudinal displacement of the n-th atom from its equilibrium position, K r are harmonic (r = and anharmonic (cubic: r = 3, quartic: r = springs. The operators Ûn satisfy the following equations of motion: Û n t = r= K r [(Ûn+1 Ûn r 1 (Ûn Ûn 1 r 1], ( 1

2 K r = K r /M, M is the mass of an atom. We suppose that an SLS with the frequency ω l < ω D is excited at the time t = at the sites n = and its neighbours (ω D is the maximum harmonic frequency, ω D = K. Anharmonic interactions are supposed to be weak, satisfying the condition a ( A, where a ( h/(ω D M is the amplitude of zero-point vibrations, A K /K is the amplitude of the self-localized vibration (from physical reasons the value of A should not exceed the value of the lattice constant d, i.e. K > K /d. The condition a ( A means that the characteristic energy of the SLS, being of the order of K /K is much larger than the characteristic vibrational quantum hω D : M 1/ K = K /K = K3/ 1. (3 hω D hk Note that the reversed dimensionless parameter K 1 characterizes the degree of the quantum anharmonicity of the lattice: it increases as the anharmonic term K increases and as the mass of the atoms M decreases. In quantum crystals He and Ne K is of the order of 1, while in ordinary crystals (characterized by small amplitude of zero-point vibrations as compared to the lattice constant it can reach many hundreds. In the problem under investigation K plays an essential role determining the time-scale of the energy relaxation (see below. In order to take account of the SLS decay we introduce the operators Ûn in the form Û n = u n (t+ξ n + ˆq n, ( where the classical displacements u n are supposed to be nearly periodic functions: u n (t A n (tcosω l t, A n (t and ξ n (t are, slowly changing with time, amplitudes and shifts satisfying the equations ω l A n = K (Ān+1 Ān+ K 3 (Ān+1 ξ n+1 Ān ξ n K ( Ā 3 n+1 = K ( ξ n+1 ξ n + K 3 ( Ā n+1 + K ( 3 Ā n+1 ξ n+1 +Ān+1ξn+1 Ā3 n Ān ξ n, (5 + ξ 3 n+1 3Ā n ξ n + ξ n+1 Ā n ξ n ξ n +, (6 wherewekeepjustthe cosωttermneglecting cos3ω l t, cos5ω l t, andintroducethe newvariables: Ā n = A n A n 1, ξ n = ξ n ξ n 1. Variations of the A n and ξ n in time, due to quantum fluctuations, are described by means of time-dependent operators ˆq n (t, which satisfy the equations d ˆq n dt = W n+1(t(ˆq n+1 ˆq n W n (t(ˆq n ˆq n 1. (7 Here W n (t = K +k n +w n cosω l t, where k n = K 3 ξn +3 K (Ā n /+ ξ n and w n = K 3 Ā n +3 K Ā n ξn determine the change of springs caused by the SLS. The terms K 3 (ˆq n ˆq n 1, K Ā n (ˆq n ˆq n 1, K ξn (ˆq n ˆq n 1 and K (ˆq n ˆq n 1 3 are supposed to be small and they were neglected. It is important to take these terms into account when ω l > ω D, i.e. when the two-phonon decay under consideration is forbidden by the energy conservation law. The phonon hamiltonian is defined as Ĥph(t = ˆT(t+Ŵ(t, where ˆT = 1 ˆq n n and Ŵ = Ŵ+Ŵ1+ Ŵ t are the operatorsof kinetic and potential energy, Ŵ = K / n (ˆq n ˆq n 1, Ŵ 1 = 1/ n k n(ˆq n ˆq n 1, Ŵ t = cosω l t n w n(ˆq n ˆq n 1 ; h = 1. Ŵ 1 describes the stationary perturbation of the phonon subsystem by the SLS, while Ŵt takes account of the oscillatory time-dependence of the springs induced by the SLS. This is the time-dependence of Ĥ(t which causes the phonon emission by the SLS (phonons are generated by pairs 11,1. Introducing properly chosen configurational coordinates ˆX i = n S inˆq n and ˆQ m = n s mnˆq n, one can diagonalize Ŵ1 and Ŵt: Ŵ 1 = 1 i η i ˆX i, Ŵ t = cosω l t m v m ˆQ m. Both Ŵ1 and Ŵt have one zero eigenvalue; all other η i are positive, while v m are sign-alternating, being symmetric with respect to

3 the sign change. The mode with zero eigenvalue is totally symmetric and therefore does not enter into ˆq n ˆq n 1. The Hamiltonian Ĥph is diagonalized by introducing time-dependent operators as follows. First the time-independent Hamiltonian Ĥ = Ĥ( ph +Ŵ is diagonalized by applying the standard methods 18,19 of local dynamics (Ĥ( ph is the Hamiltonian of the chain in harmonic approximation. Then the Hamiltonian Ĥph(t is diagonalized by the method presented in 11,1,. The time-dependent phonon Hamiltonian in the diagonal representation allows one to find energy of phonons E ph (t and the rate of phonon generation. The latter equals to 11,1 Here E l n de ph (t dt = de l(t dt ( ω l M A n + r K r ξr n r is the energy of the SLS, = I(t. (8 I = ωd ω l ω D dω m πω l v m 8 ρ m (ωρ m (ω l ω(1+n(ω 1 v m G m (ω (9 denotes the intensity of the emission of phonons, ρ m = π I(G m(ω, the temperature factor is included by the multiplier (1+n(ω, Gm (ω = G m (ωg m (ω ω l, G m (ω = i dte iωt ǫt [ˆQ m (t, ˆQ m (], ǫ, (1 is the Green s function of the mode m, ˆQm (t e iĥt ˆQm e iĥt. Use of the diagonal representation for the operator Ŵ1 allows one to easily express the Green s functions G m via the one-cite Green s function of the perfect lattice 1 G ( = i/(ω ω D ω as follows G m = G ( σ imσ i ms ins i n cos(κ(n n [ ( nn ii 1 ηi G ( ( 1 ηi G ( i i 1 δ ii η i G ( i ]. (11 Here G ( i = G ( nn S ins in cos(k(n n, σ mi = n s mns in, cosκ = 1 ω /ωd, ω < ωd. In the limit of small v m one can neglect the term vm G m (ω in the denominator of formula (9. In this limit (9 coincides with the corresponding formula of the perturbation theory 17. For realistic K 3 and the considered values of ω l between 1.3ω D and ω D the term vm G m is not small. At some ω l = ω k and ω both 1 vmr( G m (ω and I( G m (ω turn to zero and the perturbation theory fails. Near such ω k I ω l ω k 1 t t k 1/, i.e a sharp burst of phonons is generated causing a relaxation jump (t k corresponds to the time moment when ω l = ω k. Note that vm is not a linear function of the energy of the SLS. Therefore the relaxation of the SLS is always nonexponential, including the case when perturbation theory is applied. In the latter case E l (t (t l t α, whereα ,t l isthelifetime oftheslswhichisfinite inthisapproximation. AsanexamplethequantumdecayofanoddSLSisexaminedbelow. ThepropertiesoftheSLSdepend essentially on the dimensionless parameter δ = K3 /(K K. For realistic one-well pair-potentials this parameter has the value between 1 and 3/9. For Ar-Ar and Kr-Kr potentials δ = 1.37, for K-Br potential δ = For such values of δ the odd SLS is well localized if ω l > 1.3ωD. In this case SLS is localized on three central atoms: A 1 = A 1 A /, A n <.5 A, n ; ξ n = ξ n ξ, n 1, ξ = 3,. This allows one to take account of only A, A 1 and ξ 1 = ξ 1 ξ. The parameters A and ξ of the classical SLS problem in this approximation are determined by the equations (5 (n = and (6 (n = 1: ω l = 3 K +6 K 3 ξ K A +9 K ξ, (1 = K ξ + K 3 (A +ξ + K ( 7 ξ 8 A +ξ. (13 3

4 Introducing dimensionless amplitude a = K A / K and shift z = ξ K / K, one gets where z is determined by the real solutions of the equation ω l = a = z(z +zδ +1 δ z, (1 z(5z zδ + 7 δ 1 δ z (15 (here and below the units ω D = 1 are used. This is a cubic equation for z which for considered ω l and δ has only one solution with a >. This solution determines the parameters of the SLS in the classical limit. In the approximation considered, the SLS causes changes of only in central and next to central springs: k = k 1 =(16zδ+z +7a /3, k =k 1 =3a /3, w = w 1 = a(δ +3z/8, w = w 1 = aδ/8. In this case Ŵ1 and Ŵt can be diagonalized analytically. The eigenvalues η i and v m and the components of the eigenvectors are the following: η 1 =, ( η,3 = k + 3k ( ± (k 3k 1/ +k, ( η,5 = k + k ( ± (k k 1/ k, v 3 =, v 1,,,5 = ± 1 ( 8w +6w 1 +w 1w ± (16w +8w1w +16ww w1 +1w1w (++ corresponds to v 1, + to v, + to v and to v 5 ; S kn S kn = S, S1n = S k,1 = 1, S k = 1 η k /β, n kn S k3 = +η k /β, Sk = 1 η k /β, Sk5 = 1, k =,3; S k3 =, Sk = 1+η k /β, Sk5 = 1, k =,5. Energy of the SLS equals s mn 5 s mn =, s i1 = s in, n s mn n= s i = 1 v i (a i w 1 /w +1/w 1, s i3 = 1, s i = 1+a i, s i5 = 1+v i /w, a i = v i (v i +w 1 +w /(w 1 w. E l hω D Kε l ; ε l = 3ωl a +z ( 1 + z 3 + z ; (16 one sees that parameter K determines the quantum scale of the SLS energy. Note that the intensity of phonon emission does not depend on K (the unit of I in (9 is hω D. Consequently, parameter K determines the time-scale of energy relaxation. We performed calculations of the frequency dependence of the intensity of phonon emission by the SLS and of the time dependence of the intensity of phonon emission and the SLS energy (at the temperature T =. The calculation procedure was the following. First we fixed the initial SLS frequency (usually at ω l = 1.98ω D and calculated z, a, E l, G(ω, I. The second step was to take a lower value of ω l and respectively calculate new values of parameters. We repeated this procedure until the frequency of the SLS got value 1.3. The time step was calculated using the relation t = E l /I.

5 In Fig. 1 the dependence of the intensity of phonon emission on the value of ω l is shown. Indeed, the intensity I has peaks ω l ω k 1 at critical frequencies ω k (this was confirmed numerically. This means that the SLS is unstable in the vicinity of ω k. The time dependences of the intensity of phonon emission and of the mode energy E l are given in Fig.. In all cases the lowest considered frequency ω l = 1.3ω D corresponds to a small energy as compared to the initial energy for ω l ω D. Consequently, the main part of the energy of the SLS with initial frequency ω l ω D is lost during the time considered (note that an extrapolation of E l (t to larger t also gives finite lifetime of SLS. Depending on the value of δ and on the initial energy, this time can be either longer or shorter than that given by the perturbation theory. The relaxation law is essentially nonexponential: sharp bursts of phonons are generated at critical values of time t k when ω l approaches ω k causing relaxation jumps. The time dependence of the intensity I near the t k is I t t k 1/. The peaks mentioned are caused by the poles in the integrand function (9; they correspond to the emission of quasimonochromatic phonons and give sharp lines in the spectrum of generated phonons J(ω (the integrand function in (9 summed over time in Fig. 3. The SLS is not fully localized on three central atoms: it always has wings. Equations ( and (5 give the following solutions for the amplitudes A n and the shifts ξ n in the wings: A n ( 1 n A e µ n, ξ n sign(n K 3 A K e µ n, n 1 (17 where µ = arch ( ω l /ω D. It can be shown that taking into account these amplitudes and shifts leads to relaxation jumps of SLSs with ω l close to ω D (ω D ω l <.1. These relaxation jumps are weaker than those caused by the central atoms. In the case δ < 3/ there may exist large-size SLS with ω l close to ω D. The energy of such an SLS E l 8 (ω l 1K /(K (3 δ exceeds the unit ( hω D only if δ approaches 3/ from below. The problem of quantum decay of this large size SLS is analogous to that of solitons in continuous media and it will be considered elsewhere. In conclusion, we examined the quantum relaxation of the odd SLS in a one-dimensional anharmonic monatomic lattice, caused by the zero-point fluctuations of the phonons and gave an analytical nonperturbative solution of the problem. We found that the fluctuations can cause generation of very sharp bursts of quasimonochromatic phonons and very fast, step-wise jumps of the SLS energy. The rate of energy loss is on average of the order of the phonon quantum hω D per period of vibrations π/ω D. The full relaxation time is determined by the quantum anharmonicity parameter K; in chains with potentials corresponding to the quantum crystals, this time is of the order of a period of vibrations, while in lattices with potentials of ordinary crystals it may reach thousands of periods. Supported by the Estonian Science Foundation Grant No A.A.Ovchinnikov, Zh. Eksp. Teor. Fiz., 57, 63 (1969; see also A.A.Ovchinnikov and N.S.Erihman, Usp. Fiz. Nauk 138, 89 (198 [Sov. Phys. Usp. 5, 738 (198]. A.M.Kosevich, A.S.Kovalev, Zh. Eksp. Teor. Fiz. 67, 1793 (197 [Sov. Phys. JETP, 891 (197]. 3 A.S.Dolgov, FTT 8, 161 (1986 [Sov. Phys. Solid State 8, 97 (1986]. A.J.Sievers, S.Takeno, Phys. Rev. Lett. 61, 97 ( J.B.Page, Phys. Rev. B 1, 7835 ( V.M.Burlakov, S.A.Kiselev, Zh. Eksp. Teor. Fiz. 99, 156 ( G.S.Zavt, M.Wagner, and A.Lütze, Phys. Rev. E 7, 18 ( S.A.Kiselev, S.R.Bickham, and A.J.Sievers, Phys. Rev. B 5, 9135 ( K.W.Sandusky and J.B.Page, Phys. Rev. B, 5, 866 ( R.Dusi, G.Viliani, and M.Wagner, Phys. Rev. B 5, 989 ( V.Hizhnyakov, Proceedings of the XIIth Symposium on the Jahn-Teller Effect [Proc. Estonian Acad. Sci. Phys. Math., 36 (1995]. 1 V.Hizhnyakov, Phys.Rev. B 53, ( V.Hizhnyakov, D.Nevedrov, Proc. Estonian Acad. Sci. Phys. Math., 376 ( N.O.Birel and P.C.Davies, Quantum Fields in Curved Space (Cambridge University Press, Cambridge, A.A.Grib et al. Vacuum Quantum Effects in a Strong Field (Energoizdat, Moscow, 1988 (in Russian. 5

6 16 S.W.Hawking, Nature (London, 3, 3 (197; Commun. Math. Phys. 3, 199 ( P.Klemens, Phys. Rev. 1, 3 ( A.A.Maradudin et. al, Theory of Lattice Dynamics in Harmonic Approximation (Academic, New York, A.A.Maradudin, Theoretical and Experimental Aspects of the Effects of Point Defects and Disorder of the Vibrations of Crystals (Academic, New York V.Hizhnyakov, J. Phys. C: Solid State Physics, 673 ( E.N.Economou, Green s Functions in Quantum Physics (Springer-Verlag, Berlin,

7 Figure Captions Figure 1. Frequency dependence of the intensity of emission of phonons by the odd SLS.The peaks correspond to the regions where the SLS is unstable withrespect to quantum fluctuations; δ = (K 3 /K K 1/. Figure. Time dependence of the SLS energy E l (thick dashed line and intensity of phonon emission by the SLS I (thick solid line; decay of the SLS energy calculated with perturbation theory (thin solid line. Sharppeaksin the I(tdependence describeemissionofphononburstsbythesls. E l = ε l K /K, K = K /(K hω D. I is measured in units of hω D. Figure 3. Spectrum of generated phonons J(ω (in a.u.. 7

8 ε l δ = δ = δ =.8 6 t D ω /κ I

9 I δ =. δ = 1.37 δ = ωl D

10 1 δ =. J 1 δ = δ = ω/ωd