# Full text of "A model of interband radiative transitions"

## See other formats

A Model of Interband Radiative Transition J. Dittrich,"'^ P. Exner,"'^ and M. Hirokawa'^ a) Department of Theoretical Physics, Nuclear Physics Institute, Academy of Sciences, 25068 Rez, Czech Republic b) Doppler Institute, Czech Technical University, Bfehovd 7, 11519 Prague, Czech Republic c) Department of Mathematics, Faculty of Science, Okayama University, 3-1-1 Tsushima-naka, Okayama 700-8530, Japan dittr ich@uj f . cas . cz , exner @uj f . cas . cz , hirokawa@math.okayama-u.ac.jp We consider a simple model which is a caricature of a crystal inter- acting with a radiation field. The model has two bands of continuous spectrum and the particle can pass from the upper one to the lower by radiating a photon, the coupling between the excited and deexcited states being of a Friedrichs type. Under suitable regularity and an- alyticity assumptions we find the continued resolvent and show that for weak enough coupling it has a curve-type singularity in the lower halfplane which is a deformation of the upper-band spectral cut. We then find a formula for the decay amplitude and show that for a fixed energy it is approximately exponential at intermediate times, while the tail has a power-like behaviour. 1 I Introduction A rigorous description of decay and resonance processes in quantum theory has a long history starting from the Friedrichs model presented in and discussed later in numerous papers - see, e.g. [0, 0. A systematic study of the problem started in the seventies. J. Rowland and H. Baumgartel with collaborators - see [Q, 0, 0, and the papers quoted there - used operator methods to establish the existence of resonance poles and to prove the Fermi rule for various systems with perturbed embedded eigenvalues. At the same time the seminal paper ^ by J. Aguilar and J.-M. Combes initiated the development of complex-scaling methods which are nowadays a very efficient tool to study resonances of Schrodinger operators. In the eighties many papers dealing with quantum-field decay models appeared. A phenomenological models based on the Langevin equation were investigated in 0, [1^, [|ll|| and ||12[. Moreover a generalization of them was given by A. Aral |T^ within the Hamiltonian formalism. In a last few years the long-time behavior of canonical correlation functions for general Hamiltonians was investigated in |T^ by applying the results of |]13| and 15| via a quantum Langevin equation. From the point of view of virtual transitions, the long-time behavior of a correlation function was studied in | 16| . It is also worth of noticing that, revisiting the decay problem. Bach, Frohlich, and Sigal have developed a new manner to analyze the resonace problems for a class of models in quantum electrodynamics fl^, |[T8[ . In most of these models the unstable states come from perturbation of eigenvalues, either embedded in the continuous spectrum or isolated as in the case of Stark effect. Much less attention has been paid to the situation when the states which should decay belong to the continuous spectrum of the unperturbed Hamiltonian. An archetypal example of such a situation is a crystal in which an electron can radiate a photon and pass to a lower spectral band. A natural model in this case would be a Schrodinger operator with a periodic potential coupled to a quantized field. This is not easy, however. To start with a simpler case, we discuss in this paper a model of Friedrichs type with transitions between two bands of the absolutely continuous spectra which can be regarded as a one-photon approximation of the more realistic description. While perturbed embedded eigenvalues typically give rise to resonance poles in the analytically continued resolvent, we are going to show that in the mentioned model the cut-like singularity corresponding to the "excited" 2 spectral band gets deformed to the lower complex halfplane. Recall that a similar behavior has been observed in a completely different type of systems which involve a perturbation of a band spectrum, namely for scattering in finitely periodic systems [1^ . Here we have a situation with a finite number of resonances which accumulate, however, along curves in the lower half- plane which are close to the spectral bands of the infinite system when the interaction is weak. Let us describe briefly the contents of the paper. After formulating the model in the next section we shall compute in Section III the projection of the Hamiltonian resolvent onto the subspace of excited states corresponding to the upper spectral band of the "crystal". Under natural regularity as- sumptions we prove the mentioned claim about the change of the spectral singularity caused by a decay with the radiation of a "photon" . Then we turn to the time evolution of the undecayed state and show that its projection onto the upper-band subspace is - at least for a weak enough coupling - realized as multiplication by a function which we evaluate explic- itly. The rest of the paper is devoted to properties of this decay amplitude. We show that in the weak-coupling case the latter is dominated at interme- diate times by an exponential function. Hence the population of the excited spectral band changes in the course of the evolution: the wavefunction com- ponents supported in the regions where the deformed singularity is closer to the real axis survive longer. On the other hand, similarly to the usual decay theory, the deexcitation process cannot be purely exponential; we show that the decay amplitude has a power-like tail at long times. II Description of the model The "crystal part" of our model is assumed to have the simplest nontrivial spectrum consisting of a pair of disjoint absolutely continuous bands Iq = and h = [d'^e!""^] with -oo < < d""^ < d"^ < < oo. Using the spectral representation we can assume without loss of general- ity that the crystal state space is L^(JiU/o, w{x) dx) with the Hamiltonian acting as multiphcation by the variable x; the weight function w is positive a.e., Lebesgue integrable, and satisfies w{x) dx = 1 . h'Jio 3 As we have said the "field part" is represented by the vacuum and one-photon (or phonon) states, which coexist with the upper and lower band of the "crys- tal" , respectively. The photon vacuum is by assumption a single state of zero energy, while the single-photon states belong to the space L'^{[i^, oo), uj{z) dz), u >0, on which the free Hamiltonian Hp acts as a multiplication by the vari- able z. The weight function cu is again Lebesgue integrable, non-negative a.e., and satisfies poo / uj{z) dz — 1 . J u Putting the two components together we get the total state space of our model in the form H = Ho©Hi ■.^L'{h,w^{x)dx)®[L'{Io,w^{y)dy)®L'{K,uj{z)dz)], (2.1) where K — [u, oo) and Wa '■— wl'Ia, Oi — 0,1. The free Hamiltonian acts as g J \{H^®I + I® Hp)g which means with the arguments x & Ii, y E Iq, and z & K. Next we have to choose the interaction part of the Hamiltonian. Being inspired by the Priedrichs model we require (i) the interaction includes necessarily a single photon emission/absorption, or in other words, the projections of Hint on L'^{Ii,wi{x) dx) and its orthogonal complement in Ti are zero, (a) the interaction is "minimal" in the sense that the action of i^int can be written in terms of multiplication by a " formfactor" , integration, and possibly a change of variables. It follows from (i) that i?int = kL with an interaction constant k and an "off-diagonal" operator L, where : Hj — > Hi, i.e. 4 Furthermore, in accordance with (ii) the operator Lio should be chosen in the form {L,of){y,z) = X{y,z)f{u{y,z)), (2.4) where X : Iq x K —* C and u : Iq x K Ii are functions containing the dynamical information about the system. This choice in turn restricts Loi because the full Hamiltonian (with a real coupling constant k) must be symmmetric, which means / f{x){Loig){x)wi{x)dx= X{y, z)f{u(y, z))g{y, z)wo{y)uj{z) dy dz Jh J JloxK (2.5) for all / and g from the operator domain. Suppose now that there are functions u, v such that {y,z) i— > {u{y, z),v{y, z)) : Iq x K ^ Ii x K is a. bijective diffeomorphism which can be used as a substitution at the r.h.s. of ( |2.5| ) leading to {Loig){x)wi{x) = / X{y,z)g{y,z K D{y,z) D{u,v) wo{y)u{z)dt, (2.6) the variables y, z being expressed as the inverse of a; = u{y, z) and t = v{y, z) at the r.h.s. Remarks 2.1 (a) For the sake of simplicity, assume that u depends on a single variable mapping Jq onto Ji. This will reduce the dependence of the transition between a pair of states in Ji and Jq, respectively, on the photonic component of the system, (b) In the same vein we could suppose that X{y,z) = X^{y)XK{z) (2.7) which will turn Hq + i^int - up to the isomorphism between Ii and /q - into a direct integral of Friedrichs-type Hamiltonians. However, we choose a nontrivial setup and do not require that the dependence of the interaction strength on the energies of the excited state and the photon contained in the function A factorizes. In other words, we will keep a general A : Jq x — s> C After this heuristic discussion, let us define the Hamiltonian which we shall consider in the following. We suppose that (al) u : Iq —>■ Ii is a bijective C^-diffeomorphism, 5 then the interaction term i^int acts according to ( |2.3| ) with {Liof){y,z) := X{y,z)f{u{y)) , iLo,g)ix) := I ^ [ Kn-Kx),z)9{u~\x),z)u;iz)dz (2.8) \u'{u ^{x))\wi{x) Jk with X E Ii, y E Iq, and z E K. The second expression makes sense because the two factors in the denominator are positive a.e. by assumption. The operator L defined in this way is formally symmetric and unbounded in general. To get a self-adjoint Hamiltonian we add a boundedness assumption. Specifically, we assume that (a2) X is Lebesgue measurable in Jq x i^' and there are positive C, Ci such that / \X{y,z)Mz)dz<C, wo{y)<Ci\u'{y)\w,{u{y)) holds for every y E Iq] the last inequality means that the Radon-Nikodym derivative appearing as the first factor in Lqiq is bounded. Proposition 2.2 Under the assumptions (al) and (a2), H\^t is hounded and symmetric. Consequently, H = H{k) = Ho + ifint =Ho + kL is self-adjoint on the domain of Hq. Proof: It remains to verify the boundedness of Hi^t which amounts to check- ing that the operators Lio : Ho Hi and Lqi : Hi ^ Ho are bounded. This is easily seen from the following estimates: ll^io/llLx = // \X{y,z)\''\f{u{y))\'wo{y)u{z)dydz J JIqxK < C [ \fix)\' ^f",l'^l dx<CiC\\f\\l 6 and I -^015' I -fi h u)o{u \x)) X{u ^{x),z)g{u ^{x), z)u!{z) dz K X \{u-^{x),t)g{u-^{x),t)u{t) dt Wi{x) dx lo \'^'(y)\M'^(y)) X KxK X{y,z)X{y,t) g{y,z)g{y,t) uj{z)uj{t) dz dt\ dy \\{y,z)\^uj{z) dz K \g{y,t)\Mt)dt K wo{y) dy < C^C\\g\\l,^, where we have used the Pubini theorem in combination with the Schwarz inequahty for the scalar product in L'^(K x K,u;{z)u;{t) dzdt). | Before proceeding further let us make a comment on the assumptions, part physical and part technical, which we will have to make in the following. Since the present model is rather a motivation study for a more realistic one, we do not strive for the maximal possible generality. On the other hand, we do not want to impose many unnecessary restrictions which would correspond to a fully specific system such as the one given below. Example 2.3 Let -Ej(-), j = 0, 1 be the lowest two dispersion curves of a one-dimensional crystal. Since we are discussing a caricature model, we neglect the multiplicity of the eigenvalues. In other words, we consider just a half of the Brillouin zone and regard Ej as maps [0, vr] Ij with Eq strictly increasing and Ei strictly decreasing. Moreover, both are restrictions to [0, vr] of real-analytic functions with the first derivatives vanishing at the endpoints of the interval but nonzero in its interior. To rewrite the band projections of the crystal Hamiltonian in our formal- ism, we employ the operators Uj : L^([0, vr]) L'^{Ij,Wj{y) dy) defined by {Ujf){y) := f {Ej^ {%)))] the definition makes sense since the inverse functions Ej^ exist by assumption. The operators Uj are unitary provided we put Wj{y)^\E'{E-\y))\- (2.9) 7 These functions are C°° in (0, vr) with singularities at the endpoints but the latter are integrable. In particular, if E'-{'d) 7^ at "i? = 0, vr we have Wj{y) = O V^^^) there. One of the basic ingredients is, of course, the function u. Since the system of the crystal plus the radiation field is invariant w.r.t. the discrete group of translations on a multiple of the lattice constant, it is natural in the present example to suppose that the interaction does not couple states whose quasimomentum support in the upper and lower bands are disjoint. This is achieved if we choose uiy) = E, {E^\y)) ; (2.10) it is easy to see that it is a C°° function and . WM p.ii) has finite limits at ^q^"* assuming that Eq and Ei have the first non-vanishing derivative at resp. tt of the same order. On the other hand we think of the radiation field as of the electromagnetic field in the rotating wave approximation. In this case we put the threshold energy z/ = and 00 (z) = X[o,i/niax](-2^) where u-^^x is a possible ultraviolet cut-off. Under these model assumptions (al) is satisfied automatically and the same is true for the second part of (a2); it follows from ( |2.9| ) and ( |2.11| ) that it is valid for any Ci > 1. The only remaining restriction is thus the boundedness condition /q'^"'"' \X{y, z)\'^ dz < C for the formfactor. Ill The resolvent As usual the spectral information is contained in the resolvent of the Hamil- tonian. Under our assumptions, we can find it explicitly by solving the equation 9 J \ 91 for ( in the resolvent set, in particular, for all C £ C\R. It is straightforward to check that fix) = r{x,C)fi{x)-Kr{x,() '''^''''^''^^ \u'{u ^{x))\wi{x) 8 ^^"^ '^^^'""^ gi{u-\x),z)u;{z)dz, (3.1) K where y + z-C y + z-C y + z-C Let P be the projection onto the subspace Hq = L'^{Ii,Wi dx) of "undecayed" states in 7i, By (|3.1| ), the reduced resolvent acts then as muhiphcation by the function r, P{H-Cr'P = r{-X)P- (3.2) For the sake of brevity we introduce the following notation, v(y,z) := \\(y,z)\^^(z), (3.3) .""'"''ir" , ■ (3.4) \u'[u ^[x))\wi[x) GivX) ■= I ^^^rf., (3.5) so the function r can be written as r{x, C) = {x - C - K^Qix) G{u-\x), C)Y' (3.6) for $5C ^ 0. Remark 3.1 In the particular case of Example ^]3| it follows from (2^) and (|2.11| ) that q{x) = 1, and moreover, v{y,z) = \X{y, z)\'^X[o,un,..]iz). To reveal the analytic properties of r{x, ■) let us begin with those of G{y, ■). 9 Lemma 3.2 Let v{y,-) have a locally bounded derivative in (z/, oo). Then for any y & Iq and a real ( > y + i/ there exists finite principal value of the integral ny,0:-v[^^i.. (3.7) Moreover, for any k & {0,( — y — ^) , ^f-'-' v(y,z) , /■■:-»+' „(9,z)-t.fe,C-!;) Ju y + z-C Jc dz r ^(ml,, (3.8) Jc-y+k y + z-C + where all the three integrals are Lebesgue convergent. Proof: Choose any A; e (0, C — 1/ — z^)- As the integrals r-'^ivA.,, and r Ju y + z-C Jc- ,^^y+ky + z-C exist due to the assumption (a2) it is sufficient to check the convergence of fC-y+k ( \ h{y,O^V 7' V ^i^. (3.9) Jc-y-k y + z-Q We employ the identity v{y, z) — v{y, C ~ y) + [v{y, z) — v{y, ^ — y)] together with the estimate \v{y, z) - v{y, C-y)\< ci\y + ^ - CI with a finite ci independent of z. We see that finite Jc-y-k y + z-C exists and it is sufficient to check V L , , ^ which is easily seen to exist JC-y-k y+z-c •' and to be equal to zero. | As usual in similar situations to proceed one needs some analyticity as- sumption about the formfactor. In the present case we suppose that 10 (a3) for all y E Iq the function v{y, ■) can be holomorphically extended to an open set fly^y D {u, oo); we denote the extension again as v{y, ■). Let us further assume that there is an open set in C such that (^i"^ + z/, oo) C C Hyeio {y + ^v,y) ■ Notice that the hypothesis of the previous lemma is satisfied under (a3). Now we can make the following claim. Lemma 3.3 Let y E Iq and ^ E {y + u, oo). Then lim g{y, C) = I{y, ± ^T^viy, ^ - y) . Proof: Let us write again Q{y, Q defined by ( |3.5| ) as a sum of three integrals over the intervals {1^,$, —y — k), H — y — k,^ — y + k) and — y + k, 00) with 0<k<^ — y — u. The first and the third integral can be interchanged with limit by dominated convergence. The set fl^^y is open and contains (i/, 00), hence there is fci > such that any -i? G C satisfying {"^ — C,+y\ < ki belongs to Qv,y Let us consider only ( satisfying |C — ^| < ^1 (so that ( — y E ^v,y) in the second integral and denote (i = then we employ the identity v{y, z) = v{y, (i-y) + [v{y, z) - v{y, Ci - y)] and observe that \viy,z) - v{y,Ci -y)\< ci{y,^,k,ki) \y + z - Ci\ ■ The contribution from the difference can be thus also handled by dominated convergence. In view of ( ^.8[ ) we get lim g{y,0 = Hy^O + v{y,^-y) lim dz and the result follows by an easy calculation. Lemma 3.4 Define the functions Qq : Iq x Q G{yX) Qn{y,C) = { iiyX) + i^v{y,C~y) G{yX) + "^iiTviyX- y) 5>C > (QiyX) - 2ivrt;(?/, C - y) I{yX)~nxv{yX-y) G{y,0 ••• ^(<o Under our assumptions (al)-(a3), the functions Qn{y,-) (ind G^{y, holomorphic in Q \ {—00 , y + u] for any fixed y E Iq. C and : IqxQ s^c > ^Cby (3.10) (3.11) are 11 Proof: By Lemma ^.21 and assumption (a3), is a finite function. Notice tliat C — y E y for ( E and y E Iq. According to Lemma |3.3| , tlie function Qn{y, ■) is continuous in {( E fl\Q( > 0} \ {—oo,y + i^] - see, e.g., Thm 146 in Ref. Alternatively, the continuity of I{y, •) in (y + z/, oo) can be established directly from the dominated convergence used in the proof of Lemma p^. Similarly, the continuity in G f2 1 < 0}\(— oo, y + u] is seen and thus Qn{y, ■) is continuous in \ (— oo,?/ + z/]. As it is holomorphic in {( E Q\'^C > 0}U{C G fi| < 0} it is also holomorphic in oo, y+i'] due to a corollary (dubbed the edge-of-wedge theorem) of the Morera's theorem (stating that the continuous function is holomorphic iff the integrals over all rectangles with the edges parallel to the axes are zero - see, e.g., p^ , Thm 168] or [El, Thm 10.17]). As to Q^{y, ■), we can prove our statement in the same way as for ^^(y, ■). | Now we are in position to show what happens with the upper spectral band under influence of the perturbation. Let us formulate some further assumptions before. (a4) The functions g{x)QQ{u^^{x), () and ^(j;) ^^"^"^^^^)''^) continuous in the set {(x,C) E h x Q\( ^ {-oo,u~'^{x) + z/]}. (a5) For all x E h, X > u^^{x) + z/ . Remarks 3.5 (a) In the particular case of Example |2.3| the factor g{x) = 1 can be dropped in (a4) and the assumption (a5) is satisfied, (b) While most assumptions we make are of a technical nature, (a5) is a physical hypothesis saying that in no part of the excited spectral band the decay is prevented by energy conservation. It is satisfied, of course, if u = 0. Let us denote ilx = ^ \ ( — oo, U^^{x) + l'] . Theorem 3.6 Assume (al)-(a5). Then the following statements hold. (a) There exist A > 0, 5 > and a unique function ( : Ii x [—6,6) C satisfying C(x, k) E {x - A,x + A) + i{-A, A) C fi^r , (3.12) X - ({x, k) - K^g{x) gfi{u'^{x), C(x, k)) = . (3.13) 12 The function ( is continuous in Ii x (—5, 5) and Ci^-,') ^ C°°(— 5, 5). (b) The resolvent has no singularity in the upper complex half-plane, in par- ticular ^Cix,K)<0 (3.14) holds for all x & Ii, k E (—5, 6). Moreover, if g{x)v{u^^{x), X — u^'^{x)) ^ (3.15) for all X from some compact I' G Ii, then there exists a 5 G (0, 6i] such that QC{x,k)<0 (3.16) holds for all < |k| < 6i and x G Proof: (a) Let us denote D+{x, C) := X - C - K^gix) gn{u'^{x), () ■ (3.17) The functions and are continuous in {{x, k,, ()\x G /i, k G M, C £ ^x} by assumption and -D+(x, ■, ■) G C°°(M x Q^) by Lemma |3]^. Furthermore, D+{x, 0,x) = and dD+{x,0,x) ^ _^ / Q By the imphcit function theorem - see, e.g. |2T|, Thm 211] - to any Xq G Ii there exist (i^„ > 0, > and A^,;, > such that for all x G (xq — c^^.^, xq + 4Jn/iand/tG {-S^^, S^^) there is just oneCxo(a;, «) G (xq- A^^,, xo + A^J + i(— Aa;Q, Aa;^) C fix (recall (a5)) satisfying D+{x, nXxoi^: i^)) = 0; the relation (|3.13| ). The function (^o is continuous in ((xq — (i^o? + ^xo) H Ji) x (-^i^o'i^i^o) andCxo(a;,-) G C'°°(-5a;o, (^xo) for any fixed x G (xo-4o,a;o + 4o)- We put d'xo = min(Axo,4J . As Ii is compact by assumption, the open covering of Ji defined in this way has a finite subcovering, i.e. there exist a finite number of points Xj G /i, j = 1, . . . , n, such that h C U^=,ir, ; we employ here the notation Kj = {xj - , Xj + ) , Jj = Kj + i{-d'^^ , , 13 for J = 1, ... , n. Let us pick a point Xj^ G Kj n K^. for given j, A; = 1, . . . , n; then there is < 5jk < min((5-r^., S^^) such that Cx, (a;jfc, n) e n 4 , Cx^Xjk, e n Jk (3.18) for |k| < 6jk. Moreover, (x { ) for all X G KjClKk and \k\ < 6jk', otherwise the uniqueness of (xj and would be violated near at least one of the points sup{x e KjnKk\x> Xjk, Cxjiy, = Cxkiv, i^) iorxjk <y <x,\k\ < Sjk} , mf{x e Kj n Kk\x < Xjk, Cxjiy, /«) = Cxkiv, K)ioix <y < Xjk, < Sjk} ■ Choosing a number 5' > with 6' < mini<j<„ Sx^ and 6' < mmK^nKkj^iD ^jk, we conclude that there exists a unique ( : Ii x {—S', S') C such that ({x, k) G Jj for x e Kj and D^{x,k,C,{x,k)) — 0. The function C, is continuous in h X 5') and C{x, ■) e C°°{{-5\ 5')) for any fixed x e h. Put hj{x) = min(a; — Xj + d'^.,Xj + d'^^ — x) . The function hj : /i — > M defined in this way is continuous and x e Kj if and only if hj{x) > 0. Then h(x) := max hAx) l<j<n specifies a positive continuous function h on Ii. Let us denote D = min h{x) > , A = min I D, min d'^ . ) > . As ( is uniformly continuous on compact subsets of /i x {—5', 6') there exists < 5 <5' such that C{x, K)e{x-A,x + A)+ i{-A, A) for X & Ii and |k| < 5; hence the existence of the numbers 5, A and the function ( is demonstrated. Finally, to check the uniqueness of C let us assume that C is another function satisfying C{x, k) e {x - A,x + A) + i{-A, A) , D+{x, k, C{x, k)) = 14 for X G Ji, K G {—S,S). Suppose that x G Ji and \k\ < S are given. There exists an index j = 1, . . . ,n such that hj{x) = h{x) > D , X G Kj , ({x, k) G Jj . As the inequahties X — Xj + d'^. > D , Xj + d'^. — X > D , —D < y — x < D hold, where y := ^({x, k), we have also y-Xj + > , Xj + - y>0 and y G Kj. Furthermore, |53C(x, < A < d'^^ < Ar^.. Then ({x,k,) = ({x, k) and the uniqueness is proven. (b) Assume first that 3C > 0, then ^g^{u-^{x)X) = ^G{u-\x),C) > by Eqs. ( pT5| ) and ( p.lO| ), so the r.h.s of Eq. ( p. 17] ) has negative imaginary part. Consequently, there are no solutions ({x, k) with positive imaginary parts, in other words ( p.l4| ) holds. We have checked here only that the open upper half-plane is a part of the resolvent set for the Hamiltonian. In the lower half-plane, the function D+( meromorphic continuation of r{x, ■) and may have singularities. Suppose now that ( p.l5| ) holds. The expression g{x) ^^"^^ ^^).C(a^.K)) jg continuous in (x, k) G /i x {—6, 6). It follows that M := max '2'2 (x,«)G/'x[-f,|] < oo . dC Differentiating the equation defining ({x, k) with respect to we get (3.19) In combination with the previous inequality we conclude that dCjx, k) ^ g{x)gQ{u-\x)X{x, k)) 15 defining M' is continuous in (x, k) E I'x ^- 00 for M = 0. Furtliermore, the assumption ( p.l5|) together with (|3.10|) imphes '^^^^j- < 0, hence there is < 5i < min(|, M~^ ) such that 3 < for (x, k) el' X (0, 61) and ( CT) holds. | Remarks 3.7 (a) Putting k = in ( p. 19 ) we obtain 9C(x,0) (3.21) where right-hand side is given by Lemma |3.4 This relation can be regarded as an analogue of the Fermi golden rule in the present situation. (b) Notice that for the factorization (2]7) the term \\q{u^^ {x))\'^ factorizes from Qn{u^^{x), Q and C,{x, k) = x holds whenever \q{u^^{x)) = 0. IV Decay of excited states In accordance with the physical motivation, we are interested in transitions from a given state supported in /i into those in /q. To find the time profile of the de-excitation probability it is sufficient to know the reduced evolution operator PU{t)P = Pe^*^*P. Suppose that the initial state is of the form \l/r for some ipo G L'^{Ii,Wi{x) dx) with ||\l/o|P = //^ I'lpoi^)]'^ "Wiix) dx = 1. Its time evolution is given by the Stone formula, ^0= lini 7^ / [{H~^-zr])-'-{H-^ + zr]r']e-'^'^od^, according to [^, Thm VIII. 5], and the projection P can be interchanged with the limit and the integral being a bounded operator. This yields the reduced evolution operator. 16 where (4.1) and r is given by ( |3.6| ). The integral and the hmit refer to functions with values in TCo = L'^{Ii,wi{x) dx); they are known to be convergent as the Hamiltonian H is self- adjoint. Let us now look for conditions under which the interchange of the limit and the integral in ( |4.1| ) is possible. To this end, we need more assumptions. (a6) v{y,z) < C2 and dv{y,z) dz < C3 holds for some positive constants C2, C3 and all y G Iq, z G K. (a?) vijj.v) = for all y G Jq. (a8) There exists a zero-measure set N G Ii and a number ui > di := sup [x — u^^{x) — u] > such that q{x)v{u~^{x),^) > for all X G /i \ and ^ G (i^, + Ui). Lemma 4.1 Assume (al)-(a7). Then there exists a number C 4 such that \g{y,^±tr])\<C, (4.2) holds for all y e Iq, ^ eR, and eR. Proof: Recall the definition Using the first part of (a6), we get \r]\v{y,z) dz {y + z - + 7]"^ v{y, z) dz . \'^Q{y,i^i'n)\ {y + z-^f + r]- dz < C. 2 / 9 I 9 00 ^ + ^ dz = ttCo . 17 We fix q; > and distinguish several cases. (i) ^-y>'^ + a. Then where by (a2) we have I/+J/— ^ J —a a zv{y,z - y + 2-2 _|_ jj2 dz zv{y,z-y + ^) ' v+y-i J a / z'^ -\- rf' Using the mean value theorem, zv{y, z-y + ^) dz < a v{y, z) dz < C a J 2 : ^ z dz -a Z + T v{y, ^-y) + zd2v{y, i^{y, ^, z))] dz with '&{y, ^, z) between ^ — y and ^ — y + z. The integral of the first term is zero due to the antisymmetry in z while the second term can be estimated by (a6) giving | J2I < 2C3Q; and (ii) v<^ — y<v-\-a. Then Ug{y,i±irii) + / + zv{y,z-y + ^) ^ z^ + r}^ where i-v-y ^ -^—^v{y,z-y + i) dz < 2C3(C -i^-y)< 2C3a follows by the same procedure as for the integral J2 in case (i) and a Z"^ + v{y, z-y + ^)dz 18 due to (a2). In the remaining integral, \v{y,z-y + 0\ <C3iz-y + ^-iy) by (a6) and (a7). Denoting for a while A — ^ — i/ — y,we have now z pot I 2 I 2 '^(^' z-y + i)dz z^ + r)' <C3 + 7]^ , , f a A a — A-\ In— — \r]\ arctan— — arctan— 2 A'^ + r]'^ \ \r]\ \r] H i^—r^ 7, 2 + taking into account that < A < a in the last inequality. Let us estimate the maximum of function in the mentioned interval of A. Clearly /(O) = /(a) = and f{A) > for < A < a. Hence / has a maximum at some point Aq e (0, a) satisfying 42 From the last equation, /(A) ^" <Ao<a As a result. and — — ^ z-y + ^)dz < 2C3Q; |3?6;(|/,e±ir^)| <4C3Q; + Ca;-^ (^mj V — a<^ — y<u. Then |K6;(t/,e±ir7)| = 2; a / Z^+T]^ 19 v{y,z-y + C) dz Here the second integral is bounded hj Ca ^ and the first one we estimate similarly as in the case (ii). Denoting here B — u + y — ^ e (0, a], we obtain —^v{y,z-y + dz and \^g{y,^±t7])\<C3a + Ca-\ (iv) ^ — y < u — a . Then \miy,^±tv)\ = ;v{y,z-y + i) dz iv+y-e, + "n^ Summing up the discussion, we have found that in all the cases the inequality \m{y,i±ir})\ <4C3a + Ca-i holds. Minimizing the right-hand side with respect to a > 0, we get \m{y,i±iri)\<^Vccl and |6;(y,e±i77)|<v/l6CC3 + 7r2C|, what we set out to prove. | (4.3) Theorem 4.2 Assume (al)-(a8). Then there exists ^2 > such that for all < k| < 82 and t e M tl^it.x) ^U{t,x)ilJo{x) holds for almost every x E h, where POO U{t,x) = / W{x,C)e-'^'d^, K.'^g{x)v{u~^{x),^ — u~^{x)) (4.4) (4.5) W{x,0 = [x—^ — K'^g{x)I{u ^{x),^)] + TT'^K'^g{xyv{u ^{x),^—u ^{x)y (4.6) 20 Proof: Let S, A and ({x, k) = Ci( ) be as in Theorem 13.0. We first verify that C2{x, > for x & Ii\N and n / I I / A' • I X /'^i / ^ U < |k| < O2 := nun | 0, O1O4 V '-^iL'' where d := min \x — v — u ^{x\\ > by assumption (^aJj. It is sufficient to show that ({x,k,) is not real as we know that (2{x, > 0. By assumption (a2) and Lemma [4.1| , for real ^ and there is no solution in {—oo,x — K^CiCi) U {x + ti^CiCijOo). li ({x, K,) = ^ then the imaginary part of the Eq. (|3.13|) reads 'Tcg{x)v{u ^{x),^ — u ^(x)) = Thus there are no real solutions in (z/ + z/ + z^i + D (x — (i,z/ + z/i + u~^{x)) by assumption (a8). For the considered values of k the intervals without real solutions ^ cover the whole real axis. To any natural number n there exists an open set C M of Lebesgue measure smaller then - such that C A^„+i C A^„. Let us denote = Ii \ Nn- Let if be an arbitrary vector from TYq and = <^Xii,- The scalar product 1 {ipn,ip{t,-)) = lim - / ipn{x){^r{x,^ + ir]))e '^^iIjo{x)wi{x) dx V^0+ IT Jj^xR (4.7) with r{x,C, + if]) given by ( p.6|) . Fubini theorem can be used here as \g{u-\x),i + ir])\ < - , \r{x,i + i7])\ < - , rj rj \Qr{x,^ + ir])\ = C(r^) as ^^±00 ioT 1] > using (a2), ( ^.51 ) and ( |3.6D only. Let us next choose < Ai < A, < r/i < A, denote 21 and consider further only < \k,\ < S2, < rj < rji. Now \K'Qix)giu'\x),^ + ir])\ < ^Ai by Lemma [4.1| . We divide the integration range Ji x M of ([4.7|) into the parts where I^C ~ 2^1 ^ ^1 and I^C — 2^1 ^ Ai, respectively, and construct the integrable majorant allowing us to use the dominated convergence in (^4.7| ). For 1^ — x| > Ai, clearly ' ^ '^ '^'-(ie-x|-iA,r- A? • Let us define function ^ : M ^ R as (recall that h = ( - ^{-^ + ^r\m + f) ••• e<d"^-Ai 9(0= \ ••• d'^-Ai<e<d^^ + Ai (4.9) Then (x,^) 1-^ g{^)\ip{x) \ \iIjo{x) \ wi{x) is the sought majorant. For 1^ — a;| < Ai we can consider only a; G as (pn{x) = elsewhere. By Theorem |3.(j| , m-K,™ := min \D+{x, k,^ + iri)\ > where is defined in ( 3.17 ) and the minimum is taken over the considered set of variables x E I'^, C, E [x — Ai, x + Ai], i] G [0, rji] and a fixed value of < \k\ < 62 (notice that ^ + irj E fix due to our choice of Ai, rji and the inclusion in ( |3.12|) ). The majorant can be now chosen as Interchanging the limit with the integral in ( [4.7|) , using Lemma |3]^ and realizing that the integrand limit vanishes for ^ < + u~^{x), we obtain (<^„,7/>(t,-))/; = (<^„,W(t,-)^o)/; (4.10) with U given by ([4.5|)-(^l6D and scalar products in the space L'^^I'n^ Wi{x) dx). Here U(t,-)4'o G L^{r^,Wi{x) dx) as U is bounded in R x which can be seen using the majorant constructed above. As ipn ^ L'^{I'n,Wi{x) dx) may be arbitrary, Eq. ( [4. 4]) follows for a.e. x G J^. Now we see (|4.4| ) for a.e. a; G /i in the limit n —>■ 00. | 22 V Exponential decay at intermediate times Recall that decays of unstable quantum systems are nonexponential at very short and very long times, however, they are usually exponential in a very good approximation over a wide range of intermediate times. Our aim here is to show that the present models exhibits a similar behaviour in the sense that the function U{-,x) appearing in the restricted time evolution operator (|4.4| ) can be approximated by an exponential for a.e. fixed a; G /i. The way to prove that is inspired by 0. We employ the fact that the continued resolvent is for any fixed x a meromorphic function and show that for a sufficiently weak coupling the time evolution is dominated by the contribution from the residue term in ( |4.5| ). In addition to the hypotheses made above, let us assume that there exist a constant C5 such that (a9) < C5 holds for all y G Jq and z & K. Lemma 5.1 For any a > — z/, x G Ii, and ^ > u ^(x) + u the following estimates hold: +C3 In < Co a + V ^ — u ^{x) — V a + u + C,i^-u-\x)-u) dgniu-\x),0 Proof: Let us estimate dgiu~\x),c) dc °° v(u-Hx),z) iu-\x) + z-cy ■ dz (5.1) (5.2) (5.3) for ^ = ^ + i?7, > 0, and get the result on the real axis by taking the limit ?7 — > 0^ using Lemma |3.4| . We rewrite the derivative as dgn{u-\x),0 dc z"^ — rf + 2iriz _^ v{u \x),z + i-u ^{x))dz (5.4) 23 and denote for a moment V . l = i u -1 (a;) + a ; (5.5) by assumption we have < /3 < 7. In the expression for the imaginary part of (|5.4| ) we separate the integrals over {—(3,(3) and {(3, 00). In the second integral the limit 77 — >■ gives zero as can be seen easily by the dominated convergence. In the integral over {—(3,(3), we insert the Taylor expansion v{u ^{x),z + ^ u -1, x)) v{u ^{x),i-u ^{x)) dv{u ^{x),^ — u ^{x)) ld'^v{u ^{x),^ — dz 2 dz^ u X + (5.6) where 9 in the error term lies between and z. The contribution of the z° term to the integral vanishes because it gives rise to an odd function. The contribution of the second term is bounded by ttCs in the limit r/ ^ as it follows from assumption (a6) and an explicit calculation. The term again does not contribute in view of assumption (a9) and an explicit calculation. In this way, inequality ( |5.2| ) is proved. As for the real part of Eq. (|5.4|) , we proceed similarly. Inserting the expansion (|5.6|) into the integral over {—(3,(3) we obtain from the z^ term --nu {u'\x),^-u-\x)) < 2C. where the assumptions (a6) and (a7) were used in the last inequality. The term with z does not contribute and the term with z"^ is estimated by 0^(3 in the limit rj —>■ 0. The integral over (/3, 7) in ( p. 41 ) can be handled by means of (a6) and (a7). v{u^^{x), z + ^ — u~^{x)) dz {z + i- ^2 _|_ ^2^2 z'^ — rf dv{u~^{x),6i) [z^ _|_ j^2p ^z + ^- dz u ^(x u \x dz = C-i /3 In 1 + u) dz a + z/\ a + u P 7 In -1, x] a + u 24 where we have employed u < 6i < z+C,—u ^(x) and the inequality In(H-x) < 1 + I lna:|. Finally, we have ^v{u ^{x),z + ^-u ^{x))dz 7 ■ - — ■■ °° dz . C2 7 putting all these estimates together, we arrive at (|5.1|). | <C, -< a + u Lemma 5.2 There is S3 > such that for allO < \k\ < 83 and almost every X e Ii, the function W{x, •) defined by formula ( |^.6| j for ^ > z/ + u'^{x) and extended by zero to the rest of the real axis, W{x, = for C, < i' + u^^{x), is absolutely continuous in any compact subinterval of M. Proof: From the proof of Theorem 4.2 we know that W(x,^) = — lim '^r(x, ^ + irj) and therefore W{x,0 = -^rnix,0, (5.7) 71 Mx,0 ■■= 7 2 (\r ( -K ^ ,. =[D+{x,K,0]-' (5.8) X - ^ - K^g{x)gn{u ^{x},^) for ^ > z/ + u^^{x). Let S and A be the numbers from Theorem |3.6| . For K < AC I , Dj^ has no zeros if |^ — x| > A (see Lemmas |3.4| and [43 and assumption (a2)). On the other hand, for |^ — x| < A and < < ^2 real zeros can exist for at most zero-measure set of x which we neglect (see the proof of Theorem [4.2| ). Apart of it W{x^ ■) has a continuous derivative in [v + u ^(x), 00) and therefore it is absolutely continuous in any compact subinterval. Let us denote d := min \x — v — u^^ix)] > 3 where the positivity follows from assumption (a5). Then |D_|_(x, k, ^)| > for u + u'^{x) < ^ < u + u'^{x) + |, < 3^^, and ^^^^^ is bounded by an expression similar to the r.h.s. of ( ^.1|) in the considered interval of ^. Due to the integrability of | ln(^ — v — n~^(x))| and the estimate W{^.i)\<^-^{i-^-u-\x)) 25 (see assumptions (a6)-(a7) and ( ^.10| )) W{x,-) is absolutely continuous in [z/ + 'u~^(x), p + u~^{x) + |]. Consequently, it is absolutely continuous in any compact subinterval of M. Choosing 53 = min (5,^2, A d C1C4 V 3C1C4 we get the desired result. | Lemma 5.3 There exists ^4 > such that \Q{x)Qn{u''^{.x)X{x,K))\ < 00 dgn{u-^{x)X{x, k)) Ml := max x£/i,|k|<(54 Mo max x£/i,|k|<(54 Q{X} < 00 - (5.9) (5.10) and for all \k\ < := min(54, (2M2) 2), x E Ii, we have |Ci(x, k) -x\< 2MiK^ , < (2(3;, K,) < 2Mik2 \C{x,k) -x\< 2MiK^ d({x, k) < 2Mi where ({x, k) = Ci{x, k) — iC2{x, k) is the function from Theorem ^76 . Proof: By Theorem |3.6| ( is uniformly continuous in Ji x [— |, |]. Hence there is < ^4 < I such that for \k\ < 5^ and all x E h we have \C{x, k) — x\ < d. Then Ci{x, k) > u + u^^^x) and the functions in the r.h.s of Eqs. (|5^ ), (|5.10| ) are continuous. Consequently, Mi, M2 are finite. For |k| < ^5 we now have dC{x, k) Q{.x)gn{u i(x),C(a;, k)) 1 + k^q{x dGn{u~''- (x) ,({x,k)) < Ml 1 - k'^M2 < 2Mi and the sought estimates on ({x, k) — x follow. Lemma 5.4 Let a be a number such that < a < d := inf (x — u ^(x) — u a < dist {Ii,C\n) and let us denote eC\\^ -x\ <a} , {{x,^) E h X C \ X E Ji, \^-x\<a}. Then 26 (i) c {{x,^) ehxCl^ en\{-oo,u-\x) + (a) if X E Ii and 'd G N^^^ then (x, ^9) G N^, (in) Na is closed m M x C, (iv) the numbers max \p(x)Qq(u ^ (x) , ■&)\ < oo , max g[x) dgn{u-\x),^) < oo (5.11) (5.12) are finite, (v) there exists an a' > a such that for any x E Ii, d & Na',x, md \k,\ < Sei^a) := min ^^5, j (see Lemma \5.3[ ) we have gn{u'\x),^) =gn{u-\x),ax,^)) dQn{u'^{x),(ix, k)) + (^ - C(x, k)) + T{x, ^)(^ - C(x, K)Y (5.13) where T{x^ ■) is a function holomorphic in the interior of N^'^x one? 8M3(a) \q{x)J^{x, dJ^{x,^) g[x) < < m3[a) , =: 7714(0;) (5.14) (5.15) holds for i) in the interior of Na/2,x ■ Proof: The claims (i)-(iii) trivially follow from the definitions, the claim (iv) follows from the assumption (a4) and the claims (i), (iii). Under our assump- tions there exists a' > a satisfying all the assumptions of the lemma. Then for any x G Ji, the function Qn{u~^{x),-) is holomorphic in the interior of Nai,x^ the function JF defined by Eq. (|5.13|) exists and jF(x, ■) is holomorphic in the interior of Na'^x- For \n\ < 5Q{a) now (^(x, n) is in the interior of Na^x and for all d in the interior of A^^v -r we have ^(x)jF(x, 'd) Qix)gn{u-'^{x),z) 1 f 0{x)Qn{u-'^{x), z) dz , dz . 27 If \k\ < (Jfifa) and G N^^^, then a \z-C{x,k)\ > \z-x\-\x-C{x,K)\>a-2MiK^>-, (5.16) (5.17) a \z — > — a;| — |x — -i^l > — by Lemma [5.3| , and the inequahties ( |5.14| ), ( |5.15| ) follow immediately. | Theorem 5.5 Assume (al)-(a9). Then there exist finite constants 5' > and Cg > such that for all\n\ < 5' and t > we have U(t, x) — A{x, K)e' < (5.18) for a.e. x G /i where ({x, k) = Ci{x, k) — iC,2{x-, k) is the singularity location ( with C,i real, C2 > - cf. Theorem \3.(^ and 1,2/ ^^gn{u \x)X{x,k)) 1 + K g{x) dC Proof: If K = we have ({x, 0) = x by (^) and U{t, x) = e"*^* (see dg^)) so the theorem holds with any Cq. Let us further suppose that /? 7^ 0. By Theorem ^]6| and assumption (a8), (2{x,k) > for a.e. x G /i if \k\ < 62- Let us exclude the remaining zero- measure set of x's from our considerations. Then the integral where /oo pH e-'«V(x,Orf^= lim / e~'^'V{x,OdC: A{x, k) V{x,0 53 (5.19) (5.20) vr C{x,k)-C exists in the generalized sense (|5.19|) . While the Lebesgue integral does not exist due to the behavior at large |^|, the existence of generalized integral is well known and will be in fact seen from our calculations below. We shall estimate the difference between U{t,x) in Eq. ([4.5|) and the integral ( ^.19| ). Let us recall from the proof of Theorem [4.2| that , (5.21) W{x, = lim - 5> r(x, ^ + ir]) = -<^ 2 / \r t \ c\ 28 where the last equahty should be used for ^ > u + u~^{x) only. Combining this with (|4.6| ), assumptions (a2), (a6) and Lemma |4.1| we arrive at the estimate for > x + H?CiC4^. Due to Lemma p]2|we can integrate by parts for < ^3, Let us choose an a > satisfying the assumptions of Lemma and consider only the values of the coupling constants such that < |k| < min (^6s, 6^(0), \y^^J ■ (5-23) To calculate ^W{x,^) let us denote for a while D+ = D+(x, = X - ^ - K^g{x)gn{u'^{x),^) , D, = ^D+ = x-^-K^gix)^gniu~\x),0, D2 = '^D+ = -K^g{x)Qgn{u~\x),0, f)J~) f) D[ = ^ = -l-n'g{x)^^gniu-\x),0, D', = ^ = -K'g{x)^-g^{u-\x),0- ^W{x,0 = -Tv-'lD+r" [{Dl - Dl)D'^ - 2D,D[D2] . If now 1^ — x| > I the assumption (a2) together with Lemmas ^]T| and [O (where we denote the constant as C3) imply a Then \D', > -/s:"CiC4 > - > 0, - 21"- I' < \^ - x\ + K^Cid < - x\ , < C's + K^CiCsl ln(e - u'\x) - z/)| + K^CiC^i^ - u-\x) - v) 29 From here we get -^Wix,0\d^<C7K'. where the exphcit value of the constant C7 can be expressed from the above estimates if necessary. What is important is that C7 can be chosen indepen- dent of K in the considered range. Let us consider the term V{x,^) now. We have the bounds 1 1 1 - K^M2 by Lemma 5.3 so < \A{x,k)\ < 2 holds for (5.24) (5.25) \k\ < min((54, (2M2)-^/2) . Denoting for a while Ai = ^A{x, k), A2 = QA{x, k), we have l^il < \A{x,i)\ < 2, IA2I < k'^M2\A{x,k)\^ < Ak'^M2 and d_ _ 1 A2U - Ci(x, K)f - (2(3:, Kf] - 2ArU^, K){i - Ci(x, k)) If 1^ - x| > f and IkI < min ( 5\ a 8M1 (5.26) we have \i — Ci{x^ k)\ > j by Lemma |0, and therefore i with a «;-independent finite constant Cg which can be given explicitly if nec- essary. 30 Let us now turn to^ G (x — |,a; + f). Using the expansion (|5.13|) , W{x,i)-V{x,i) 3 d ^\W{x,i)-V{x,0] = -K^'^iAix.K TT X g{x)'-Z^ - n^A{x,K)Q{xf:F{x,if [1 + K)Q{x){i - C(x, K))T{x, 0? Using (|5.24| ), ( |5.14| ), ( |5.15| ) together with Lemma and assuming that 1 \k\ < min I ^^{a), I we obtain where ^[Wix,0-Vix,0] 16 2V Ml ' 2^ms{a)a^ (5.27) Cg := — (mJa) + 26Q(aYm3(aY) . TT ^ ^ Putting all the estimates together, we get roo o J e-^^'—[W{x,0-V{x,0]d^ where and Ce = Cv + Ca + Cga k| < S' 5' being the minimum of 5, 82 and the r.h.s. in ( ^.23|) , (|5.25|) -( pl27D . Evalu- ating the generalized integral e-'^'V{x,C)di = A{x, K)e by closing the integration contour in the lower half-plane for t > the in- equality ( p.l8| ) is obtained. | 31 fi:-r/2Tie-" "^^^ = Cg , ^ = 1, 2. (5.30) The theorem is apparently useless for very short and very large times when the error estimate 0{k,H~^) is much larger then the amplitude value ~ exp{—(2{x,K)t). On the other hand, we get a nontrivial bound for the times when ^ « e-f^^-''^)* (5.28) where we take into account that A{x, n) ~ 1. Let us write C2(x, k) = K'^ri2{x,K), ?72(x, k) = TTg{x)v{u~^{x), X — u~^{x)) + 0{k'^) (5.29) for small coupling k. In the subsequent formulas we do not write the argu- ments of T]2, however, its x-dependence should be kept in mind in general. The relation ( |5.28| ) is valid for Ti <^ t <^ T2 where Ti, T2 are two solutions of the equation If K^Ti is small we can approximate the equation by replacing the exponential with one obtaining On the other hand, if ^^772 ^2 ^ 1 we do not enlarge the range (Ti,T2) by dropping the linear factor in ( |5.30|) . Then we obtain T2^-^ln(C6A2). The r.h.s. here is an decreasing function of 772 in the interval {0,Cq^k~^). By ( ^.29|) and assumptions (a2), (a6) we have < r/2 < 7rCiC2 in the k*^ approximation. Restricting ourselves then to the coupling constant values with \k\ < {tiCiC2Cq)-^'\ we can safely use 1 ln(7rCiC2C6/€^ 7rCiC2K2 Hence we see that the announced approximately exponential behaviour of U{-,x) holds in the weak-coupling regime over wide time range, roughly _ 1 _i 3 speaking from "C^ "C^n to Cf ^Ca'^K'^. 32 VI Long time behavior The fact that the bound given by Theorem |5.5| becomes useless at very large times is not coincidental, because the decay rate is indeed slower there. To illustrate this claim, for instance, let x G /i be such that by Lemma |3.4|, Theorem |3.6| (b). Theorem |4.2| and assumptions (a6)-(a8), we have W{x,^) is finite and continuous w.r.t. ^ G [z/ + u^^ (x) , oo) (6.1) for < \k,\ < 62, where 62 is the number from Theorem |4.2| . This holds for a.e. X G /i. By (a6) and ([4. 61), we get K'^g{x)C2 [x~^-K^g{x)I{u-^{x),Of Since lim^^oo -^(^^^(a^), = by (|37^ ), we get P^- (6-2) Thus, g{^) := X[u+u-^{x),oo){C)W{x,^) is in L^(M), but its support is not the whole M, and l^f (-.01 < . , no- POO U{t,x) = / g{i)e-'^'di by ( [4.5| ). Applying now Corollary C2], we find that for a.e. a; G /i and \k\ < 62 U{t,x) does not decay exponentially as |t| 00. (6.3) To learn more about the long-time asymptotic behavior of U{t, x), we adopt the conditions (al0)-(al3) below, and employ the results of [|13| and ||I5| in the same way as in Thm 3.2(ii)]. Given i/ > and 6 G (0, 7r/2), we define g by D^,e ■={CeC\^C>v,-e < argC < 0} . (6.4) If < z/', we have therefore ^u'fi C D,,,. (6.5) Let us denote ^{v) :=n,e/o(^-?/)- (6.6) Notice that VL{v) C ^y^io^v,y by (a3). We shall assume: 33 (alO) There exists Oq E (0,7r/4) such that Du,eo *^ ^{v). (all) v{y,^) > holds for each y E Iq and ^ > z/. (al2) Given y G Jq, there exists Cy > and > such that \viyX)\<Cy\cry holds for any ( G Q{v). Notice that for v which is continuous by (a3), the assumption (all) implies, in particular, that for each x G Ji and a, /5 G (z/, oo) we have I^x,a,f3 ■= ^inf ^ V {u~^ (x),^) > . (6.7) For fixed x E Ii and ( E —oo, {x) + p], we have defined -D+(x, k, () by ( |3.17| ). In a similar way, we define three other functions, D_{x, kX), W{xX), and ^^(C) by D4x,kX) := x-C-K^Q{x)g''{u-\x)X), (6.8) ^ '^^ ■ D4x,kX)D.{x,kX) ' ^ ^ ^ (x, C + (x)) D_ {x, K,,( + (x)) (6.10) in the last case ( E Q{v) \ (— oo, u]. Then, for a.e. x E Ii and k G M with < |k| < ^2, can be regarded as measurable with g.j. E {{u, oo) , d^) (6-11) by ( |6.1| ) and ( |6.2| ), and we can write the time evolution as follows, POO Uit,x) = e-^'^-'(^)* / g^iOe-'^'d^, (6.12) J u by (pD and (g^ Next we need several lemmas. The first of them follows from (aS), (alO), and Lemma 3.4: 34 Lemma 6.1 QxiC) is meromorphic in I^ufio f^i" every x G /i and k G M. Lemma 6.2 For every x G /i with g{x) 7^ 0, ^ G M with ^ > v, and k G with < |k| < 82, (6.13) Proof: Let ^' = ^ + u^^{x). By Lemma |3.4| we have DAx,K,i') := lim (6.14) = X~^' -K'Q{x)^l{u-\x),i')+l'Kv{u-\x),0^, D4x,K,^') lim D_(x, K, — is) (6.15) = x-e'-/s:2f?(a;){/(n-i(x),e') -^vrt;(n-Ma;),0} for ^' > M^^(x) + z/ (^ > z/), which imphes that D+(x,K,0^-(a;,'^,0 (6.16) = [x-^'- k^q{x)I {u-^ {x) , 0] ' + ^^t^^Q{xfv {u-^ (x) , e)' . Then hm£^o+ W{x, ^' - ie) = W{x, follows from ( |4^ ) giving ( |6.13| ). | Lemma 6.3 For every x G Ii, with q{x) 7^ 0, all sufficiently small e > 0, every a,(3 & ('^i 00) wi/i a < and every k G R wi/i < < 62, there exists a constant Cx,a,i3 > independent of e such that sup \gx{^- is)\ < Cx,a,f3- (6.17) Proof: Set S, p,<i •" {( eC\p<^( <q, -z/tan^o < < 0}. Fix e' G with < < 1 arbitrarily, v {u~^ (x) , ■) is uniformly continuous in Sa,i3 by (aS) and ('aiOj since Sa,f3 C Dj^gg. So there exists a constant ei = ei{x,e') > such that <6' {u ^ (x) , ^ - ie) - (m ^ (x) , ^) for a < ^ < /5 and < £ < £1 and we have {u-'{x),^-ie) < {l + e')v{u~'{x),^) V {u ^ (x) , ^) (6.18) 35 for a < ^ < /3 and < e < ei. Since D±{x,k,,-) is holomorphic in f2 \ ( — oo, u^^{x) + u] by Lemma [3l^ , D±{ uniformly continuous in Sa,f3+ u^^{x). In view of ( p.l4| ) and ( |6.15| ) there exists 62 = e2{x,e') > such that \D^ix, K, - te) - D^{x, K, O I < s'\D±ix, k, ^ I for ^' = ^ + u~^{x) with a < ^ < /5 and < e < 82- Hence we have {1-e') D±{x,K,C + u-^{x)) < D±{x,K,^ + u~^{x)-i£) (6.19) a a <^ < P andO < e < 62. Using further Q, (|6l0|) , ( |6l6D , ( |6l8D , and ( |6.19| ), we get {l + e')^''g{x)\v{u-'{x),0\ < < (1 - E'y\D^ (x, K, e + n-1 (x)) D_ (x, /£, e + u-^ (x)) I (l-£')'vr2fi:4^(x)2|t;(M-i(x),0P (1 + ^0 (1 -£:')^Aix,«,/3vr2K2^(x) ioT a < < P and < £ < £0 = min{£:i,e2}5 which imphes the desired result. I Lemma 6.4 For every x G h, all sufficiently large \(\ with ( G Dy^g^,, and every /? G M satisfying < \k\ < 62, \9. (C) I < C 10 |C|2+<?a wt/i a constant Cio > independent 0/ G D Proof: In this proof, we set ?/ = m ^(x), ^' = ^ + m ^(x), and let ^ > z/ > 0. Since L)_(x, K, ^' - ze) = X - - ie) - K^g{x)g{y, ^ - ie) for every e > 0, we get \D_ (x, K, e + (x) - P > + A,,. (0)' 36 where Set Ae,. (0 := K'Q{x)^g{y, i' - le) + ^-^i x] — X . Bx ■■= K^^i(x)C4 + \u'^{x)\ + |a;| > . Then we get |y4e_3;(^)| < by Lemma ^TTl Since we now take ^ > 0, we get for every C_ with < C_ < 1, Bl l-C •2 • Thus there exists C_ with < C_ < 1 and = ^-(x) > independent of £ > such that for every > As for -D+(x, k, ^' — ie), we have K, ^' — ie:) = -D-(x, k, ^' — ie) — q{x)t[v {y, ^ — is) for any e > 0. Moreover, by (al2) we get (6.20) 1^2 _|_ ^2 jSy/S ^(Jh for ^ — ie E VL{y). Thus there exists .^^ > independent of e > such that if ^ > then C_ Together we get (x, ^ + m"^ (x) — ie) > 4:'kk?'q{x) [x, K,^ + u (x) — ie) for ^ > max l} =: by ( |6.20D ; notice that ,^_|_ is independent of e > 0. On the other hand, we get ^ 2 - 2 ^ 37 for ^ > Now we set C+ := C_/2; then < C+ < 1 and (x, K, ^ + M ^ {x) — is) (6.21) for every ^ > Put ^ = 3f?C and — ?7 = so that r/ > 0. Then, having ^ > r/, we get 2^^ _ (^^2 ^ ^ ^2 _ ^2 > q_ Ylence by (|;20D and (|;2g) we obtain C± ^ ^ |/^±(x,K,C + ^-M^))l for C = ^ — ZT] with ,^ > max 77). If C ^ '^yfio with > max we have Using then (|6.10|) , (al2), and ( |6.22| ), we arrive at for sufficiently large \(\ with ( G Tiufio- I Next we set for any x E Ii dl = X - [y + u ^ (x)) - k^q{x) ^v{u-^{x),z) dz . z — V (6.22) (6.23) Remark 6.5 Recall that by (a5) is positive for sufficiently small Let us finally state the last assumption: (al3) Given x G /i, there are constants Ai, .j. 7^ and p^.^^ > such that v{u-^x) ,C + u) _ 71 — ^u.x ■ We set X — V — u (x) g{x)I {u ^{x),i' + u ^ (x)) {x — V ~ u ^ {x)) I q{x) v{u~^{x),z) dz z — V -1 (6.24) By {a5), this quantity satisfies kI^> and d^ ^ for /t^ 7^ k^^. 38 Lemma 6.6 Assume (al)-(a3), (al) and (alO). Then lim D±{x,K,^ - ir] + u^^{x)) = dl- (6.25) Proof: By ([XT7| ), and Lemma K,^ — if] + u^^{x)) = X — ^ — 'U^"'^(x) + i?7 —K^g{x)[Q{u^^{x),^ — ir] + u^^{x)) + 2iTTv{u^^{x), ^ — ir])] , D_{x, — ir] + u~^{x)) = X — C, — u^-^i^x) + ir] —n'^g{x)Q{u^^{x),C, — if] + u^-^i^x)) for Tj > 0, ^—ir] G Q which is the sufficient range of variables as i>+u~^{x) e fl by (alO). Under assumption (alO), there exists A > such that v{u~^{x), •) is holomorphic in the set {( E C\\( — < 2A}. Taking into account (a7) then hm v{u~'^{x),^ — ir]) = v^u'-^^x), u) = . Let us write ^/ -1/ w • , -1/ XX r^"^ v{u-\x),z) -v{u-\x),^ -ir]) Q[u {x),t,-ir] + u [x)) = / : dz Ju Z-^ + lT] v+A dz f°° v{u ^{x),z) + v{u-\x),i-ir^) / 7^. + / ' dz. (6.26) For the ffist and third integral, dominated convergence theorem can be used giving (recall (a7)) °^ v{u-\x),z) z — V as the limit of their sum as ^ — > z/+, ?7 ^ 0+. The second integral '^+^ dz 1, (z/ + A-n2 + r/2 = — In z — ^ + ir] 2 (^ — z/)2 + ■ f V V \ + r arctan h arctan tt V u+A-^ i-v J for 1/ < ^ < z/ + A, ?7 > 0. As \v{u^^{x),^ ^ ^^7)1 ^ c^/{^^^liY^^rf for ^ — ir] in a neighborhood of u with a suitable constant c, the limit of the second term in ( |6.26D is zero. Now ( |6.25| ) is seen. | 39 Lemma 6.7 Let < |k| < ^2, x G /i and ^ 0. Then the function has no poles in {C^ G Dj,_eg | |(^ — z/| < Sq] with a constant > ^^^^^ '^he limit Wy^^ := hm — — (6.27) 2 (x, K, Z/ + (x)) D_ (x, K, Z/ + (x)) dl Proof: The poles of QxiC) come only from the zeroes of D± (x, ft, C + "^"^ If < ^ then lim g^iC) = (6.28) by ( |6.10|) and Lemma 3^. By Lemma 3J-, (^3; is meromorphic in Dj^^p so its only possible singularities there are isolated poles; they also do not accumu- late at V due to ( |6.28| ). Thus gx{C) has no poles in a small neighborhood of C, = V in D,^^e(j. By (al3), we therefore have Now we can formulate the main theorem of this section: Theorem 6.8 Assume (al)-(a7), (al0)-(al3). Then for every x G /i and K G M satisfying g{x) > 0, df, 0, < \k\ < 62 we have the following asymptotic behaviour: U{t, x) ~ ^^i^,xe~^[''+" Hx)]t^-iAp.,^+i)/2Y ^^^^^ + 1) ^-(p.. +1) t—rOO where F is the gamma function. Proof: It is sufficient to apply [13, Theorem 2.1(b)] to ( |6.12| ) with the help of Lemmas |6.1[ - |6.7| and we obtain the desired result. | Acknowledgment The research has been partially supported by GA ASCR and Czech Ministry of Education under the contracts A1048801 and ME170. M.H. was supported by Grant-in- Aid 11740109 for Encouragement of Young Scientists from Japan Society for the Promotion of Science. 40 References [1] K.O. Friedrichs: "On the perturbation of continuous spectra" , Commun. (Pure and) Appl. Math. 1, 361-406 (1948). [2] H. Araki, Y. Munakata, M. Kawaguchi, T. Goto: "Quantum field theory of unstable particles", Progr. Theor. Phys. 17, 419-442 (1957). [3] J. Dittrich, P. Exner: "A non-relativistic model of two-particle decay I-IV", Czech. J. Phys. B37, 503-515, 1028-1034 (1987); B38, 591-610 (1988) ; B39, 121-138 (1989). [4] J. Rowland: "Puiseux series for resonances at an embedded eigenvalue". Pacific J. Math. 55, 157-176 (1974). [5] J. Rowland: "The Livsic matrix in perturbation theory" , J. Math. Anal. Appl. 50, 415-437 (1975). [6] R. Baumgartel, M. Demuth: "Perturbation of unstable eigenvalues of finite multiphcity", J. Funct. Anal. 22, 187-203 (1976). [7] R. Baumgartel, M. Demuth, M. WoUenberg: "On the equality of res- onances (poles of the scattering amplitude) and virtual poles". Math. Nachr. 86, 167-174 (1978). [8] J. Aguilar, J.-M. Combes: "A class of analytic perturbations for one- body Schrodinger Ramiltonians" , Commun. Math. Phys. 22, 269-279 (1971). [9] R. Grabert, U. Weiss, P. Talkner: "Quantum theory of the damped harmonic oscillator", Z. Phys. B55, 87-94 (1984). [10] C. Aslangul, N. Pitter, D. Saint- James: "Time behavior of the cor- relation functions in a simple dissipative model", J. Stat. Phys. 40, 167-189 (1985). [11] R. Jung, G.-L. Ingold, R. Grabert: "Long-time tails in quantum Brow- nian motion", Phys. Rev. A32, 2510-2512 (1985). [12] E. Braun, P.A. Mello: "The correlation function for a quantum oscillator in a low-temperature heat bath", Physica A143, 547-567 (1987). [13] A. Aral: "Long-time behavior of two-point functions of a quantum har- monic oscillator interacting with bosons", J. Math. Phys. 30, 1277-1288 (1989) . [14] M. Rirokawa: "An inverse problem in quantum field theory and canon- ical correlation functions. An application of a solvable model called the rotating wave approximation", J. Math. Soc. Japan 51, 337-369 (1999). 41 [15] A. Arai: "Spectral analysis of a quantum harmonic oscillator coupled to infinitely many scalar bosons", J. Math. Anal. Appl. 140, 270-288 (1989). [16] E. Karpov, I. Prigogine, T. Petrosky, G. Pronko: "Friedrichs model with virtual transitions. Exact solution and indirect spectroscopy" , J. Math. Phys. 41, 118-131 (2000). [17] V. Bach, J. Frohlich, I. M. Sigal: "Renormalization group analysis of spectral problems in quantum field theory". Adv. Math. 137, 205-298 (1998). [18] V. Bach, J. Frohhch, I. M. Sigal: "Return to equihbrium", J. Math. Phys. 41, 3895-4060 (2000). [19] F. Barra, P. Gaspard: "Scattering in periodic systems: from resonances to band structure", J. Phys. A32, 3357-3375 (1999). [20] M. Reed, B. Simon: Methods of Modern Mathematical Physics. I. Func- tional Analysis., Academic Press, New York, 1972. [21] V. Jarmk: Differential Calculus II, Academia, Prague, 1956 (in Czech). [22] I. Cerny: Foundations of Analysis in the Complex Domain, Academia, Prague, 1967 (in Czech). [23] W. Rudin: Real and Complex Analysis, McGraw-Hill, New York, 1974. 42