getting back to the spirit of the xix@xmath0 century  when purely mechanical models were ubiquitous even for understanding systems involving electromagnetic fields  this paper discusses a simple model of an oscillator coupled to a string that presents a host of interesting features in a very accessible way .
it will be presented in detail in    [ sec : model ] ( see figure  [ fig : cordosc ] ) . even though the required technical background is on an upper undergraduate / early graduate level ,
if one remains within a classical ( _ i.e. _ non - quantum ) context , this model encapsulates many relevant features of physical interest which can be found in more realistic models , for instance when considering the interaction of an atom ( the oscillator ) with light ( the string ) .
first (   [ sec : scattering ] ) , when the string is infinite , it allows one to gain insight into ( or to discover for the first time in academic studies ) the scattering of waves by a dynamical system .
it provides the opportunity to introduce some of the ingredients of scattering theory : when dealing with a one - dimensional mechanical system , the notion of asymptotic states ( the so - called `` in '' and `` out '' states ) is made as simple as possible and the @xmath1 matrix takes a particular simple form . unlike the transmission of light in a fabry - perot interferometer , the resonant scattering appears here explicitly in connection with the familiar resonance phenomenon of a forced oscillator and we can also understand how the motion of the oscillator affects in return the external excitation .
second (   [ sec : radiation ] ) , from the point of view of the oscillator , this model constitutes the most elementary example of radiation .
it shows concretely how an interaction not only induces a shift in the natural frequency of the oscillator ( this can already be seen when only two degrees of freedom are coupled ) but also that the coupling to a large number of degrees of freedom ( the string being seen as a large collection of oscillators ) induces a friction term for the oscillator , although no dissipation exists in the system as a whole .
indeed , as far as i know , there are very few places where this model is discussed and always in the specialised literature with a dalembert string @xcite .
some of its variants @xcite are introduced precisely for studying dissipation at a quantum level .
as we will see , a klein - gordon string allows to keep one discrete mode ( a bounded state ) without dissolving it in the continuous spectrum of the string and therefore allows to mimic the interaction of a field with a stable particle , not just a metastable one .
this model is particularly relevant to see how backreaction works : for instance , for a dalembert string , the oscillator is governed by an abraham - lorentz - dirac - like equation . besides
, we will be able to illustrate precisely the deep connection between the resonance and the poles of the @xmath1 matrix ( a major feature in high energy particle physics and in condensed matter physics ) .    before we give some guidelines for further developments in the conclusion (   [ sec : conclusion ] ) , we will complete our classical study in
[ sec : hamdiag ] by the detailed diagonalization of the hamiltonian and the discussion of the completeness of the basis of modes that are used to described the dynamics
. this will be the occasion to sketch the finite size effects if one wants to use this model for describing an atom placed in a cavity .
this study also prepares the ground for the quantization which will be proposed in a future paper .
model of a harmonic oscillator coupled to a string via a massless spring .
the long - dashed lines refer to the equilibrium positions of the string and the oscillator .
the inset shows the three forces that cancel at the massless attachment point  @xmath2 ( @xmath3 ) : the left and right tensions of the string and the elastic force of the coupling spring .
, width=340 ]
a very thin homogeneous string of linear mass  @xmath4 is considered to have only transverse displacements in one direction . at equilibrium
it forms a straight line along the @xmath5-axis with the uniform tension  @xmath6 ( we will never take into account the effects of gravity ) .
the string is coupled to a system of mass  @xmath7 with one degree of freedom connected to a fixed support with a massless spring of stiffness  @xmath8 ( see figure  [ fig : cordosc ] ) .
the coupling is modelised by a second massless spring of stiffness  @xmath9 attached to the string at the massless point  @xmath2 at  @xmath3 .
all the vibrations will be considered within the harmonic approximation of small amplitudes .
we will denote by  @xmath10 the transverse displacement at time  @xmath11 of the string element located at  @xmath5 at equilibrium .
the displacement of  @xmath7 with respect to its equilibrium position will be denoted by  @xmath12 .
the equation of motion of the oscillator is    @xmath13    where @xmath14 is the displacement of  @xmath2 and @xmath15 is the harmonic frequency of the free oscillator . for  @xmath16
, @xmath17 fulfills the one - dimensional dalembert equation @xmath18 where the wave velocity on the string is given by  @xmath19 .
the derivation of   is a major step in every introductory course on waves ( * ? ? ?
* , for instance ) .    a mechanical model that leads to a klein - gordon equation : the string is attached to an elastic support that can be seen as a tight collection of massless springs uniformly distributed per unit length along the string . for another mechanical example ,
see @xcite.,width=321 ]    less common is perhaps the refinement of sticking the string to a `` mattress '' ( see figure  [ fig : cordematelas1 ] ) made of  @xmath20 massless springs per unit length along the  @xmath5-axis , each of them having a stiffness  @xmath21 .
the restoring force per unit length due to the mattress , @xmath22 , turns the dalembert equation into a klein - gordon equation @xmath23 where @xmath24 .
this equation governs also the propagation of the electromagnetic field in a rectangular waveguide @xcite and is also the relativistic equation of a quantum particle of mass  @xmath25 , @xmath26 being the velocity of light in the vacuum .    since the attachment point  @xmath2 is massless , the sum of the three forces applied to it vanishes as depicted in the inset of figure  [ fig : cordosc ] . to first order in the slope  @xmath27 ,
the transverse component of the right tension applied to  @xmath2 is given by  @xmath28 the limit when  @xmath29 keeping  @xmath30 .
on the other side , the transverse component of the left tension is  @xmath31 .
the restoring force of the coupling string corresponds to the opposite of the left hand side of , namely  @xmath32 .
therefore the coupling introduces a discontinuity of the slope of the string at  @xmath3 : @xmath33      the fundamental units will be chosen to be @xmath34 for the times , @xmath35 for the lengths and @xmath7 for the masses .
the model is therefore uniquely determined in terms of dimensionless quantities defined by the appropriate rescaling : @xmath36 , @xmath37 , @xmath38 , @xmath39 , @xmath40 , @xmath41 , etc .
these conventions correspond to @xmath42 .
there is no need for any quantization as long as the effective planck constant  @xmath43 . for simplifying the notations , in the following we will drop the `` effective '' subscript and work directly with  @xmath44 , @xmath45 , @xmath46 and @xmath47 . introducing the shifted that  @xmath7 actually feels . ]
frequency @xmath48 the equations governing the dynamics of the system become    [ subeq : motiona ] @xmath49 @xmath50 and @xmath51 with the help of the dirac distribution  @xmath52 , equations   and   can be combined in one equation , valid for all  @xmath5 , @xmath53 indeed , is recovered after integrating between  @xmath54 and  @xmath55 when  @xmath56 , since  @xmath17 and its time derivatives are continuous ( and finite ) everywhere .
if one wants to prepare the ground for some perturbative treatment of some non - linear corrections , if one wants to quantize the model and/or to couple it to a thermal bath , one possible starting point is the hamiltonian of the system expressed in term of some canonical variables .
the continuous part of the system ( the string ) corresponds to a hamiltonian density @xcite involving a pair of canonically conjugate fields  @xmath57 whereas the oscillator is described in terms of  @xmath58 .
the poisson bracket has also a mixed structure of continuous and discrete variables : for any two  @xmath59 , @xmath60 that are functions of  @xmath58 and functionals of  @xmath61 , @xmath62 the hamiltonian of the whole system generates the evolution of any function(al )  @xmath63 via  @xmath64 . for our system
we have @xmath65{\mathrm{d}}x\;;\ ] ] the corresponding hamilton equations yield directly to .
the interaction term is completely different from the one used in a recent one - dimensional pedagogical model  @xcite .
the latter emphasizes the recoil effect of the field on the mass  @xmath7 ( unlike in the present work , the external ( translational ) degrees of freedom for @xmath7 are considered ) whereas we are more interested in the resonance effects .
as for any vibrating system , a normal mode is defined to be a particular collective motion where all the degrees of freedom oscillate with the same frequency . in the absence of coupling  ( @xmath66 ) , for a given frequency  @xmath67 , one can choose two independent normal modes ( the free modes ) on the string given by @xmath68 the wave number is obtained from the dispersion relation of the klein - gordon equation : @xmath69 when  @xmath70 , @xmath71 is purely imaginary with a positive imaginary part and does not correspond to any traveling wave .
when @xmath72 , @xmath71 is real positive and the two modes are two monochromatic counter - propagating waves .
a real field  @xmath17 obeying the klein - gordon equation   carries an linear energy density @xmath73 $ ] and a linear density of energy current @xmath74 . the conservation of energy takes the form  @xmath75 for  @xmath76 . for a monochromatic traveling wave of complex amplitude  @xmath77 , @xmath78
, an elementary calculation shows that the average current over one period is given by @xmath79 whereas for an evanescent wave @xmath80      when a coupling is present ( @xmath81 ) , the two free waves   are not solutions of anymore .
a monochromatic traveling wave will be partially reflected ( resp . transmitted ) by the oscillator with a reflection ( resp .
transmission ) coefficient  @xmath82 ( resp .
@xmath83 ) that is a complex function of  @xmath71 or @xmath84 .
the linearity of equations guarantees that the frequency will be unchanged by scattering . indeed , to describe such a scattering process , a relevant choice for the two modes of frequency  @xmath84 is to look for in - states , defined ( for real positive  @xmath71 ) to be of the form ( see figure  [ fig : inoutt ] a ) ) :    [ subeq : inmodek ] @xmath85\scriptscriptstyle \rightarrow}}}(x , t)=\frac{1}{\sqrt{2\pi } }    \begin{cases}{\mathrm{e}^{{{\mathrm{i}}}(kx-\omega t)}}+\rho\ , { \mathrm{e}^{{{\mathrm{i}}}(-kx-\omega t ) } }                 & \mathrm{for}\ x\leqslant0\;;\\            \tau\ , { \mathrm{e}^{{{\mathrm{i}}}(kx-\omega t ) } }            & \mathrm{for}\ x\geqslant0\ ; ;   \end{cases}\ ] ] and , since the scattering is symmetric with respect to  @xmath86 , @xmath87\scriptscriptstyle   \leftarrow}}}(x , t){\stackrel{\mathrm{def}}{=}}\xi^{\mathrm{in}}_{{\substack{k \\[-.7ex]\scriptscriptstyle \rightarrow}}}(-x , t)=\frac{1}{\sqrt{2\pi } }    \begin{cases }   \tau\ , { \mathrm{e}^{{{\mathrm{i}}}(-kx-\omega t ) } }                 & \mathrm{for}\ x\leqslant0\;;\\     { \mathrm{e}^{{{\mathrm{i}}}(-kx-\omega t)}}+\rho\ , { \mathrm{e}^{{{\mathrm{i}}}(kx-\omega t ) } }                 & \mathrm{for}\ x\geqslant0\ ; ;   \end{cases}\ ] ]     several choices for the normal modes with frequency  @xmath84
. the wavy arrows symbolise monochromatic traveling waves whose complex amplitude is specified below .
a ) the in - state given by  ; b ) the in - state given by  ; c ) the out - state given by   and d )  shows the general coefficients that will be linked by the @xmath1 or @xmath88 matrix defined respectively by   and  . ]    with , for both modes , the same amplitude  @xmath89 for the oscillator that can be interpreted , using the language of linear response theory , as a susceptibility : @xmath90 these modes form a set of waves that is suited for constructing localised wave - packets that look like free wave - packets when  @xmath91 ( _ i.e. _ far away from the oscillator ) .
for instance , when considering a localised wave - packet traveling to the right , we shall have is concentrated around  @xmath92 , the wave - packets travel with the group velocity  @xmath93 . when  @xmath91 , the dominant contributions to the integrals in form one wave - packet located in  @xmath94 and traveling to the right .
] : @xmath95\scriptscriptstyle \rightarrow}}}(x , t)\,{\mathrm{d}}k \xrightarrow[t\to-\infty ] { } \int\tilde\varphi(k)\,\xi^{\mathrm{fr}}_k(x , t)\,{\mathrm{d}}k= \frac{1}{\sqrt{2\pi}}\int\tilde\varphi(k)\,{\mathrm{e}^{{{\mathrm{i}}}(kx-\omega t)}}{\mathrm{d}}k\;.\ ] ] the continuity of  @xmath17 at  @xmath3 implies that @xmath96 conservation of energy implies that the average energetic current is conserved in the stationary regime . from , we must have @xmath97 moreover , the equations are real and invariant under the time reversal , therefore if  @xmath10 is a solution , so is its complex conjugate  @xmath98 and  @xmath99 .
if we consider the solution @xmath100\scriptscriptstyle   \leftarrow}}}(x ,- t)\big)^*$ ] , then we obtain the mode depicted in  [ fig : inoutt ] c ) obtained from  [ fig : inoutt ] b ) by reversing the orientation of the arrows and by conjugating the amplitudes .
this procedure defines the out - modes that behave like free modes for remote future times when packed in localised superpositions .
we will have ( see figure  [ fig : inoutt ] c ) )    [ subeq : outmodek ] @xmath101\scriptscriptstyle \rightarrow}}}(x , t){\stackrel{\mathrm{def}}{=}}\left(\xi^{\mathrm{in}}_{{\substack{k \\[-.7ex]\scriptscriptstyle   \leftarrow}}}(x ,- t)\right)^ *   = \frac{1}{\sqrt{2\pi } }    \begin{cases }   \tau^*\ , { \mathrm{e}^{{{\mathrm{i}}}(kx-\omega t ) } }                 & \mathrm{for}\ x\leqslant0\;;\\     { \mathrm{e}^{{{\mathrm{i}}}(kx-\omega t)}}+\rho^*\ , { \mathrm{e}^{{{\mathrm{i}}}(-kx-\omega t ) } }                 & \mathrm{for}\ x\geqslant0\ ; ,   \end{cases}\ ] ] and @xmath102\scriptscriptstyle   \leftarrow}}}(x , t){\stackrel{\mathrm{def}}{=}}\left(\xi^{\mathrm{in}}_{{\substack{k \\[-.7ex]\scriptscriptstyle \rightarrow}}}(x ,- t)\right)^ *   = \frac{1}{\sqrt{2\pi } }   \begin{cases}{\mathrm{e}^{{{\mathrm{i}}}(-kx-\omega t)}}+\rho^*\ , { \mathrm{e}^{{{\mathrm{i}}}(kx-\omega t ) } }                 & \mathrm{for}\ x\leqslant0\;;\\            \tau^*\ , { \mathrm{e}^{{{\mathrm{i}}}(-kx-\omega t ) } }            & \mathrm{for}\ x\geqslant0\ ; ;   \end{cases}\ ] ]    with @xmath103 if we interpret  @xmath104\scriptscriptstyle \rightarrow}}}$ ] as a superposition of in - modes coming from both sides that conspire to product no wave traveling to the right for  @xmath94 , we get @xmath105    more mathematically , @xmath104\scriptscriptstyle \rightarrow}}}(x , t)$ ] can be seen as the continuation of  @xmath106\scriptscriptstyle   \leftarrow}}}(x , t)$ ] to the domain of negative  @xmath71 s . indeed , we have @xmath104\scriptscriptstyle \rightarrow}}}(x , t)=\xi^{\mathrm{in}}_{{\substack{-k \\[-.7ex ] \leftarrow}}}(x , t)$ ] provided that we define @xmath107      the in - modes and the out - modes are two possible bases for describing a scattered wave - packet
. these bases can be obtained one from each other by linear transformations ; the linearity of the equations of our model implies that they connect waves with the same frequency only , which is a major simplification .
a typical scattering experiment consists in preparing one wave - packet traveling towards the scatterer ( the oscillator ) .
long before the diffusion , this ingoing wave - packet is a simple superposition of in - modes .
long after the diffusion , we get two outgoing wave - packets that are naturally described in term of out - modes .
the passage from the in - basis to the out - basis is described in term of the scattering matrix  @xmath1 that encapsulates all the information about the possible scattering processes , to get a complete basis one may add to the in - states   the bounded modes when existing . ] .
it is made of @xmath108 blocks  @xmath109 defined by @xmath110\scriptscriptstyle \rightarrow}}}\\[1ex]\xi^{\mathrm{out}}_{{\substack{k \\[-.7ex]\scriptscriptstyle   \leftarrow}}}\end{pmatrix } = s(\omega ) \begin{pmatrix}\xi^{\mathrm{in}}_{{\substack{k \\[-.7ex]\scriptscriptstyle \rightarrow}}}\\[1ex]\xi^{\mathrm{in}}_{{\substack{k \\[-.7ex]\scriptscriptstyle   \leftarrow}}}\end{pmatrix}.\ ] ] the decomposition of each out - mode in term of the two in - modes for , say , @xmath111 , leads to @xmath112                 { \displaystyle}-\frac{\rho}{\rho+\tau}&{\displaystyle}\frac{\tau}{\rho+\tau } \end{pmatrix } = \begin{pmatrix}\tau^ * & \rho^*\\ \rho^*&\tau^ * \end{pmatrix}.\ ] ] in other words , for the general monochromatic wave @xmath113 the @xmath1 matrix connects linearly the coefficients : @xmath114 in the absence of scattering ( @xmath115 , @xmath116 ) , @xmath1 simply reduces to the identity .
the unitarity of  @xmath1 , which can be checked on  , can be seen as a direct consequence of the conservation of energy since , from  , the norm of the two vectors involved in   is preserved : @xmath117 .    if one wants to calculate the diffusion by several scatterers , it is more convenient to introduce the transfer matrix  @xmath88 whose @xmath108 blocks are defined to connect the left coefficients to the right coefficients , @xmath118 then we have @xmath119                 { \displaystyle}-\frac{\rho}{\tau}&{\displaystyle}\frac{1}{\tau } \end{pmatrix}\ ] ] whose determinant is one .
the addition of one scatterer on the string corresponds to a multiplication by a  @xmath88 matrix .
the definitions and the general properties presented in     [ subsec : inoutdef ] and [ subsec : stdef ] are valid for any non - dissipative punctual scatterer .
as far as our model is concerned , inserting the expression   with the oscillation   in equations yields to a linear system that can be solved straightforwardly : @xmath120 @xmath121 with @xmath122 and @xmath123    even though the coefficients  @xmath82 and  @xmath83 were first defined for traveling waves , _
i.e. _ for  @xmath72 , the above expressions can be continued for  @xmath70 .
we will understand the physical interpretation of this procedure when we will study radiation in    [ sec : radiation ] .
as long as  @xmath72 , , and hold .
we can check that the ultra - violet limit  @xmath124 is equivalent to the limit of weak coupling where the oscillator becomes transparent :  @xmath125 and  @xmath126 .
another case where the coupling is inefficient is when  @xmath2 and  @xmath127 both oscillate in phase with  @xmath128 since the coupling spring remains unstretched .
graphs of @xmath129 ( dashed line ) and  @xmath130 ( solid line ) given by   and   when  @xmath131 .
one can observe a resonance scattering spike for  @xmath132 and an antiresonance for  @xmath133 where no scattering occurs .
the quality factor is a ) @xmath134 and b )  @xmath135 .
let us mention that when  @xmath136 , the antiresonance has vanished and the local maximum  @xmath137 is too soft to be called a resonance . when @xmath138 , no scattering resonance occurs and  @xmath139 decreases monotonically from 1 at  @xmath140 to 0 when  @xmath141 . ]    more interesting is the resonant scattering that occurs , provided that @xmath142 , when the ingoing wave that forces the oscillator has precisely the same frequency as the shifted frequency of the latter , that is  @xmath132 .
then , the scattering is the most efficient since no transmission occur ( @xmath143 , @xmath144 ) .
the resonance spike can be seen in figure  [ fig : rhotau ] and its quality factor can be evaluated from its width  @xmath145 when  @xmath146 : for a small coupling , @xmath147 and therefore , the smaller  @xmath9 , the better the quality of the resonance .
the general idea that damping and therefore irreversibility emerge because of the interaction with a large number of degrees of freedom can be illustrated explicitly on our model .
if we choose initial conditions such that the spring is at rest at @xmath148 , the entire energy being contained in the oscillator , for instance @xmath149 the energy transfer to the string will damp the oscillations of  @xmath7 and the latter may completely lose its energy far before the energy can get back from the string if its boundary is far away from the oscillator ( for a string of length  @xmath150 , poincar recurrence time is of order @xmath151 if there were no dispersion ) .
the equation of motion of the oscillator is particularly simple when the string is non - dispersive ( @xmath152 ) and therefore we will start by studying this case .
however , we will also consider the case of the klein - gordon string because when  @xmath153 , we will see that there exists a stable mode of the oscillator at a frequency  @xmath154 whose dissipation is blocked because  @xmath155 .
its vibration does not decay because at this frequency , only evanescent waves can exist on the string , which do not carry away energy current on average  ( see  ) .      with the initial conditions  ,
the general form of the radiated waves on the string will be  @xmath156 : each of the two wave - packets travels away from  @xmath3 without distortion when  @xmath152 .
equation   becomes @xmath157 and the elimination of  @xmath158 from   and yields to @xmath159 the two terms in the right hand side are dissipative forces .
the first one has the familiar taste of the viscous resistive force whereas the second has the flavour of the schott term  @xmath160 ( * ? ? ?
* eq . ( 2.7b ) ) in the abraham - lorentz - dirac equation which governs the dynamics of an electric charge  @xmath161 that takes into account the electromagnetic self - force of the charge .
the major difference is the sign of the coefficient in front of the third - derivative . unlike the schott term ,
the negative sign in   prevents the spurious exponentially accelerating solutions .
it can be clearly seen that the irreversibility due to dissipation comes straightforwardly from the choice of initial conditions that break the time - reversal symmetry under which the original equations are invariant .
many models of an oscillator coupled to one - dimensional waves are recovered in the limit of strong coupling  @xmath162 ( see the references given in the introduction , for instance when the mass is directly attached on the string ) . in that case , only the viscous force remains in   and we immediately get the well - known damped oscillator @xmath163 .    looking for exponential solutions  @xmath164 yields to the characteristic equation of  : @xmath165
there are three solutions , one real @xmath166 and two complex  @xmath167 all having a strictly negative real part .
perturbatively in  @xmath9 , we have    @xmath168    for generic initial conditions , including  , where the string is at rest the energy of the oscillator will exponentially decay like  @xmath169 at the rate @xmath170 in the language of particle physics , the stable non - interacting particle ( the mode of the free oscillator ) has been destabilised into a metastable particle of lifetime  @xmath171 because of its interactions .
when  @xmath172 , one can not get a differential equation for  @xmath12 but must keep working with its temporal fourier transform @xmath173 together with a superposition of purely radiated waves of the form  @xmath174 . inserting them in , @xmath175 must satisfy . ] @xmath176\tilde{x}(\omega)=0\;.\ ] ] therefore  @xmath175 vanishes everywhere but at the frequencies that cancel the brackets .
these are precisely the poles of  @xmath83 and therefore of  @xmath177 given by   and  . indeed , for pure radiative modes , the ingoing waves vanish ( @xmath178 ) and therefore the matrix element  @xmath179 must go to infinity in order to keep  @xmath180 finite ( see figure  [ fig : inoutt ] d ) and equations   and  ) . letting  @xmath181
, we look for the solutions of the cubic equation @xmath182 perturbatively in  @xmath9 , those are    @xmath183    when  @xmath184 , we recover  @xmath185 and  @xmath186 .
the physical frequencies will be the three square roots  @xmath187,@xmath188 whose imaginary part is not positive : the typical decay rate of energy will be given by the nearest root  @xmath189 to the real axis : @xmath190 . as long as  @xmath191 ,
all the three frequencies have strictly negative real part .
the decay rate is given by @xmath192 it is a very general feature that the poles of the  @xmath1 matrix are associated with resonances and , more precisely , that their imaginary part provide the decay rates , which are proportional to the inverse of the quality factor of the resonances .
when  @xmath193 , @xmath194 is negative , one residual oscillation persists at frequency  @xmath195 . for  @xmath155 ,
no transfer of energy is allowed ; only evanescent waves are created and those do not carry any average energy current . unlike the scattering states , this non - decaying mode is spatially localised .
more generally , any bounded mode has a purely real frequency  @xmath196 that must be less than  @xmath197 since  @xmath198 in order to fulfil  . from  , @xmath199 is therefore purely imaginary .
moreover , in order to cancel the bracket in  , @xmath200 must lie in between  @xmath201 and  @xmath202 . for simplicity ,
let us introduce the auxiliary parameter @xmath203 and the real positive variable  @xmath204 ; a stable mode will exist if the cubic equation @xmath205 has a positive real solution .
this can be achieved for  @xmath206 only _ i.e. _ in a regime where  @xmath207 . for @xmath206 , the product of the roots of the left hand side of  , @xmath208 is  @xmath209 . if two roots are complex conjugated , therefore the third one is necessarily positive .
if the three roots are real , either only one is positive or all three of them are .
the latter case must be ruled out since the sum  @xmath210 is strictly negative .
we have therefore proved that a sufficient and necessary condition for a stable mode to exist is that  @xmath153 .
its frequency is given by   @xmath211 where @xmath212 is the unique positive real solution of  .
perturbatively in  @xmath9 , we have @xmath213 and the corresponding bounded mode is given by    [ subeq : bmode ] @xmath214 @xmath215    the choice of the normalization , @xmath216 will be justified below ( equation [ eq : xibnorm ] ) .
what makes the model completely tractable is of course that it remains linear .
however , the direct diagonalization of the quadratic hamiltonian   remains particularly difficult . in that case , the trick is to solve the equations of motion to determine the normal modes first  this is precisely what we have done in the previous paragraphs  and then write the hamiltonian in its diagonal form , @xmath217 in term of some ( real ) canonical coordinate  @xmath218 or ( complex ) normal coordinates  @xmath219 that are associated with modes labelled by the discrete and/or continuous index  @xmath220 to recall which indices are running and what is their range  @xmath221 if the latter does matter . ] .
we have @xmath222 @xmath223 and , for each pair  @xmath224 , @xmath225 @xmath226 where  @xmath52 stands for the kronecker symbol or the dirac distribution .
the second step consists in determining the canonical transformation that expresses @xmath227 in terms of some a priori known normal coordinates , namely some free normal modes  @xmath228 . in our case
this transformation is linear and will be transposed directly into the quantum theory by replacing the complex number  @xmath227 ( resp .
@xmath229 ) by the creation ( resp .
annihilation ) operator  @xmath230 ( resp .
its hermitian conjugate  @xmath231 ) of the  @xmath220th one - particle eigenstate whose energy is  @xmath232 .
as we have seen , all the scattering states are twice degenerate , in the sense that each normal frequency  @xmath84 is associated with two independent states labelled by  @xmath71 and  @xmath233 .
these modes both diagonalize the hamiltonian  .
an infinite number of pairs of eigenvectors can be chosen to constitute a basis , among them , the in and out - states , which are particularly relevant as soon as we get into a quantum field theory the normal coordinates  @xmath234 and  @xmath235 constructed from the scattering modes , once quantized into  @xmath236 and  @xmath237 , allow the interpretation of the quantum states in term of asymptotic ( quasi-)particles ; more precisely the linear transformations from the free  @xmath238 provide the explicit connection between the non - interacting states ( the fock space for bare particles including the free vacuum ) and the interacting states ( the fock space for dressed particles including the interacting vacuum ) . ] .
but it order to get   properly one must check that the set of modes is actually complete  _ i.e. _ that any kind of motion of our system can be described as a linear superposition of modes  and orthonormalized correctly in order to deal with canonical complex coordinates .
fourier analysis assures that the free states   constitute a complete set for describing the waves on the string . when interacting with the oscillator , if  @xmath153 one bounded state exists that must be added to the in - modes ( or to the out - modes ) to get a genuine basis .
then , including the normal coordinates  @xmath239 of the bounded mode if there is any , reads @xmath240 we chose the convention that , when  @xmath241 , @xmath234 ( resp .
@xmath235 ) is constructed from @xmath106\scriptscriptstyle \rightarrow}}}$ ] ( resp .
@xmath104\scriptscriptstyle \rightarrow}}}$ ] ) while @xmath242 ( resp .
@xmath243 ) is constructed from @xmath106\scriptscriptstyle   \leftarrow}}}$ ] ( resp .
@xmath104\scriptscriptstyle   \leftarrow}}}$ ] ) .      having a complete set of modes
does not guarantee that they are orthogonal .
indeed , it may happen that two eigenvectors having a common eigenfrequency are not .
for instance , one must check in one way or another that the modes   are orthogonal and properly normalized . if we denote by  @xmath244 a _ classical _ state represented by the displacement  @xmath12 of the oscillator and the wave  @xmath10 on the string , the scalar product between two states  @xmath245 and  @xmath246 is @xmath247 it is shown in the appendix that if @xmath248 ( resp .
@xmath249 ) stands for the mode @xmath106\scriptscriptstyle \rightarrow}}}(x , t)$ ] ( resp .
@xmath106\scriptscriptstyle   \leftarrow}}}(x , t)$ ] ) both with  @xmath250 , then we have , for any ( positive and/or negative ) real pair  @xmath251 , @xmath252 it is easy to see that we chose the normalization in order to get @xmath253 the free states @xmath254 represented by  @xmath255 and   are clearly orthonormalized ,  @xmath256 , and form a complete basis if we add the state that allow to describe the motion of the oscillator , namely @xmath257 represented by  @xmath258 and  @xmath259 .
the potential in   is a real definite positive symmetric quadratic form and therefore can be diagonalized in an orthogonal basis of _ real _ vectors .
the natural choice of retaining the real or the imaginary part of  @xmath260 actually provides two real modes but that are not orthogonal .
a way to assure that we deal with an orthogonal basis , is to pick up a symmetry , say the parity  @xmath86 , of the hamiltonian and classify the eigenmodes accordingly . the bounded state , if there is any , remains even .
the in and out modes are not symmetric under space inversion but it is straightforward to obtain eigenmodes that are also eigenvectors of parity . for any complex factors
@xmath261 , the combinations @xmath262 are symmetric / antisymmetric eigenvectors at any time with the eigenvalues  @xmath263 . after some algebraic manipulations using the expressions   of  @xmath82 in term of  @xmath264 given by  , the symmetric and antisymmetric modes are represented , for  @xmath241 , by    @xmath265\scriptscriptstyle \rightarrow}}}^{\mathrm{in}}(x
, t)+\xi_{{\substack{k \\[-.7ex]\scriptscriptstyle   \leftarrow}}}^{\mathrm{in}}(x , t)\big)&= & c_+\,{\mathrm{e}^{-{{\mathrm{i}}}\eta}}\,\sqrt{\frac{2}{\pi}}\cos(k|x|-\eta)\,{\mathrm{e}^{-{{\mathrm{i}}}\omega t}}\ ; ; \\    c_-\big(\xi_{{\substack{k \\[-.7ex]\scriptscriptstyle \rightarrow}}}^{\mathrm{in}}(x , t)-\xi_{{\substack{k \\[-.7ex]\scriptscriptstyle   \leftarrow}}}^{\mathrm{in}}(x , t)\big)&= & { { \mathrm{i}}}c_-\sqrt{\frac{2}{\pi}}\,\sin(kx)\,{\mathrm{e}^{-{{\mathrm{i}}}\omega t}}\;.\end{aligned}\ ] ]    the choice  @xmath266 and  @xmath267 leads to the real normalized symmetric modes , defined for  @xmath241 by @xmath268 with the unitary matrix @xmath269                 { \displaystyle}-{{\mathrm{i}}}&{\displaystyle}{{\mathrm{i}}}\end{pmatrix}.\ ] ]    constructing the real ( anti)symmetric states from the out - modes at any time leads to the same @xmath270 since the common eigenspace to @xmath271 and to the parity is of dimension one . to sum up , for  @xmath241 , @xmath270
is represented by    [ eq : xi+ ] ^+_k(x , t)&=^-t ; & x^+(t)&=^-t ; + ^-_k(x , t)&=(kx)^-t ; & x^-(t)&=0 ,    and  @xmath272 is , at any time , an orthonormalized eigenbasis of symmetric or antisymmetric eigenvectors of the hamiltonian with eigenvalues given respectively by and  .
the explicit canonical linear transformation that connects the free canonical variables to the interacting ones ( the in or out modes via the real symmetric ones ) is beyond the scope of this article and will be given and extensively interpreted in a future paper where we will quantize our model . as explained above ( see   [ subec : hamiltonian ] and also the note [ fn : qft ] ) , this is really interesting and beyond a purely academic exercise only if one wants to switch to quantum theory and/or statistical physics .
the quantum linear transformation between  @xmath236 and  @xmath238 appears to be a generalised bogoliubov transformation and our model provides an explicit construction of quasi - particles in terms of free particles .
graphical resolution of equations   that provide the even frequency spectrum when the oscillator is attached to the middle of string of finite size  @xmath150 .
the frequencies , represented by small up triangles on an horizontal line at the bottom of the figure , are centered on the abscissae of the intersections of the graph @xmath273 ( thick solid line ) with the graphs @xmath274 ( thin solid lines ) for  @xmath72 . as in figure
[ fig : rhotau ] a ) , the quality factor of the resonance is  @xmath134 .
the down triangles indicate the spectrum of vibration of the klein - gordon string alone of finite size  @xmath150 obtained for @xmath275 , that is for  @xmath276 with  @xmath20 a positive integer . ]
however , the real symmetric modes that have been founded in the previous section remain interesting at the less advanced level of the present article because they are the natural modes to work with when the finite size  @xmath150 of the string becomes relevant . indeed , when @xmath151 is not too large compared to the typical time @xmath171 characterising the radiations of the oscillator , the discrete character of the spectrum of the non - interacting string can be `` feeled '' by the oscillator .
when boundary conditions are imposed , say  @xmath277 , the discrete ( even ) spectrum of the whole system is modified by the presence of the oscillator and , from given by the  @xmath278 that fulfill the equations @xmath279 that can be solved graphically ( figure  [ fig : cavite ] ) .
the frequency  @xmath280 of the free oscillations of the mass inserts in the spectrum of the string .
the even spectrum will differ from the non - interacting case when  @xmath281 is significantly different from one . for resonances with high quality
, it will not affect the frequencies that are away from the resonant frequency .
what one gets here , for  @xmath152 is an elementary model of an atom in a ( perfect ) electrodynamics cavity of size  @xmath150 ( some imperfections can be taken into account if we relax the dirichlet boundary conditions and put partially reflectives `` mirrors '' on the string ) .
the field may or may not be quantized and , not to speak of lasers , we obtain a sort of primer for the widespread physics of quantum electrodynamics cavities that have been realised in laboratory to test successfully some fundamental concepts in quantum physics @xcite .
the purely mechanical model for infinite  @xmath9 ( the mass is directly attached on the string ) has been carefully studied with experiments in @xcite .
in addition to a more detailed study of the finite size effects , another natural development of the present work would be to deal with multiple scatterers .
for instance , when there are two identical scatterers with  @xmath153 , we expect that the degeneracies of the two bounded modes is broken and that the splitting between the symmetric and the antisymmetric bounded modes decreases exponentially with the separation of the oscillators . starting with initial conditions where only one oscillator has some energy , the beating between the two oscillators is an example of tunnelling due to the presence of evanescent waves connecting the two oscillators
.    even before we quantize the whole system , our model may be interesting to keep the field classical whereas only the oscillator is quantized .
it would provide an illustration of say , the fermi golden rule within the context of time - dependent perturbation theory @xcite .
however , we have proven that this golden rule transpires in our classical model since the transition rate   to the continuum of the modes is proportionnal to the square of the coupling strength  @xmath9 to first order in the perturbation .
as it has been demonstrated , this model captures many fondamental phenomena that are important in many areas of physics and offers wide possibilities for pedagogical use .
above all , i hope it will help the reader to discover and/or to understand them more deeply .
many thanks to domique delande and benot grmaud for their hospitality at the laboratoire kastler - brossel ( paris ) and to stam nicolis ( laboratoire de mathmatiques et de physique thorique , tours ) for his careful reading of the manuscript .
the construction of the real symmetric modes presented in   [ subsec : realsymmetric ] leads to an orthogonal basis  @xmath272 .
when it exists ( @xmath153 ) , the bounded state @xmath282 has a norm unity .
that they are eigenvectors for different eigenvalues of  @xmath271 or of the parity guarantees that  @xmath283 and this appendix proves that the proportionality factor is indeed unity .
rewriting   with the help of   and with the heaviside step function  @xmath284 , @xmath285\scriptscriptstyle \rightarrow}}}(x , t)=\frac{1}{\sqrt{2\pi } } \left[{\mathrm{e}^{{{\mathrm{i}}}kx } } + \rho(k)\,{\mathrm{e}^{-{{\mathrm{i}}}kx}}\theta(-x)+ \rho(k)\,{\mathrm{e}^{{{\mathrm{i}}}kx}}\theta(x)\right]{\mathrm{e}^{-{{\mathrm{i}}}\omega t}}\;,\ ] ] we have , with  @xmath286 and @xmath287 , @xmath288\scriptscriptstyle \rightarrow}}}^*(x , t)\,\xi^{}_{{\substack{k_2 \\[-.5ex]\scriptscriptstyle \rightarrow}}}(x , t)\,{\mathrm{d}}x=\delta(k_1-k_2 ) + \frac{{{\mathrm{i}}}}{2\pi}\left [   \rho^{}_2\left(\frac{1}{k_2+k_1+{{\mathrm{i}}}0^+}+\frac{1}{k_2-k_1+{{\mathrm{i}}}0^+}\right)\right.\\ + \rho^*_1\left(\frac{1}{k_2-k_1+{{\mathrm{i}}}0^+}-\frac{1}{k_2+k_1-{{\mathrm{i}}}0^+}\right ) \left.+\frac{2\rho^*_1\rho^{}_2}{k_2-k_1+{{\mathrm{i}}}0^+ }                   \right].\end{gathered}\ ] ] we have used the identity ( @xmath289 stands for the limit  @xmath290 keeping  @xmath291 ) @xmath292 valid for any real  @xmath71 .
the other identity @xmath293 allows to convert in terms of the @xmath52 distribution and of the cauchy principal value  @xmath294 : @xmath295\scriptscriptstyle \rightarrow}}}^*(x , t)\,\xi^{}_{{\substack{k_2 \\[-.5ex]\scriptscriptstyle \rightarrow}}}(x , t)\,{\mathrm{d}}x=\delta(k_1-k_2)\\ + \rho^*_1\delta(k_1+k_2 ) + \frac{{{\mathrm{i}}}}{2\pi}\left [ ( \rho^{}_2-\rho^*_1)\frac{\wp}{k_2+k_1}+(\rho^*_1+\rho_2 + 2\rho^*_1\rho_2)\frac{\wp}{k_2-k_1}\right].\end{gathered}\ ] ] in fact , the coefficients of the principal values both vanish when the respective denominators cancel ( we use , and , then  @xmath296 follows ) and we can drop the symbol  @xmath294 .
the  @xmath297 can also be forgotten for , to constitute the basis , we retain only strictly positive values of  @xmath298 and  @xmath299 . a little bit of algebraic jugglery with and allows to check that @xmath300 = -\frac{\kappa^2}{2\pi}\,\frac{\tau_1^*\tau^{}_2}{(\omega_\kappa^2-\omega_1 ^ 2)(\omega_\kappa^2-\omega_2 ^ 2 ) } = -\chi^*(\omega_1)\chi(\omega_2)\ ] ] with @xmath301 , @xmath302 ( @xmath303 ) .
then we have proved that , for  @xmath304 and  @xmath305 , @xmath306\scriptscriptstyle \rightarrow}}}^*(x , t)\,\xi^{}_{{\substack{k_2 \\[-.5ex]\scriptscriptstyle \rightarrow}}}(x , t)\,{\mathrm{d}}x=\delta(k_1-k_2)\;.\ ] ] the space inversion of this identity leads immediately to  @xmath307 .
at last , with the same techniques we can obtain @xmath308\scriptscriptstyle \rightarrow}}}^*(x , t)\,\xi^{}_{{\substack{k_2 \\[-.5ex]\scriptscriptstyle \leftarrow}}}(x , t)\,{\mathrm{d}}x=(1+\rho^*_1)\delta(k_1+k_2)\ + \\ \frac{{{\mathrm{i}}}}{2\pi}\left[\frac{2\tau^{}_2\tau^*_1-\tau^{}_2-\tau^*_1}{k_2-k_1}+\frac{\tau^{}_2-\tau^*_1}{k_2+k_1 } \right].\end{gathered}\ ] ] as above , for  @xmath298 and @xmath299 both strictly positive , @xmath297 vanishes whereas the second term on the right hand side is precisely  @xmath309 .
therefore @xmath310\scriptscriptstyle \rightarrow}}}^*(x , t)\,\xi^{}_{{\substack{k_2 \\[-.5ex]\scriptscriptstyle \leftarrow}}}(x , t)\,{\mathrm{d}}x=0\;.\ ] ] the complex conjugation and the time reversal  @xmath311 of the above relations allows to show that the out - modes are also orthonormalised .
the orthonormalization of  @xmath312 follows from and . to sum up
, we have obtained the following scalar products , for any  @xmath304 and  @xmath305 : @xmath313 @xmath314 @xmath315 and , when  @xmath153 for a unique bounded state to exist , @xmath316 @xmath317 three eigenbases for the hamiltonian have been chosen : @xmath318 , @xmath319 and @xmath272 . the passage from one to the other
is done with unitary matrices made of independent @xmath108 unitary blocks of  @xmath109 or  @xmath320 given by   and   respectively .