ONE-DIMENSIONAL SOLITON DYNAMICS IN T H E P R E S E N C E OF A PINNING POTENTIAL By Benoit Leduc B . Sc. (Mathematiques-Physique) Universite de Montreal, 1993 A T H E S I S S U B M I T T E D IN P A R T I A L F U L F I L L M E N T O F T H E REQUIREMENTS FOR T H E DEGREE OF M A S T E R OF SCIENCE in T H E F A C U L T Y O F G R A D U A T E STUDIES D E P A R T M E N T O F PHYSICS We accept this thesis as conforming to the required standard T H E U N I V E R S I T Y O F BRITISH C O L U M B I A December 1995 © Benoit Leduc, 1995 In presenting this degree at the thesis in University of partial fulfilment of of department this thesis for or by his or scholarly purposes may be her representatives. permission. Department The University of British Columbia Vancouver, Canada for an advanced Library shall make it agree that permission for extensive It publication of this thesis for financial gain shall not DE-6 (2/88) requirements British Columbia, I agree that the freely available for reference and study. I further copying the is granted by the understood that head of copying my or be allowed without my written Abstract The dissipative effect of (what could be modeled as) a pinning center on the lowamplitude motion of a quasi-one-dimensional ferromagnetic soliton (domain wall) at temperature T = 0 is investigated. T h e method of collective coordinates is used to eliminate the problem of the zero-frequency eigenmode and the question is, subsequently, cast into a Caldeira-Leggett form to calculate the parametrized spectral density function. It is shown that this function equals zero below a specific frequency w i t h the consequence that magnons, i n the presence of a domain wall, have no dissipative effect at T = 0 on a slow moving wall trapped by a soft pinning center. ii Table of Contents Abstract ii Table of Contents iii List of Figures v Acknowledgements vi Dedication vii Epigraphs viii 1 2 Introduction 1 1.1 T h e quantum measurement paradox 3 1.2 The microscopic quantum tunnelling 8 1.3 T h e Caldeira-Leggett formalism 10 1.4 The domain wall and its dynamics 14 1.5 T h e Sine-Gordon equation 22 1.6 Sources of dissipation . 29 A theoretical model for ferromagnetic "wires w i t h pinning centers 2.1 A model for some types of easy-axis ferromagnetic 31 quasi-one-dimensional magnets 31 2.2 Quantization about classical solutions 38 2.3 Toy-model for pinning centers in magnetic wires 46 iii 3 4 5 Soliton bremsstrahlung and the m e t h o d of collective coordinates 49 3.1 Classical treatment of the soliton bremsstrahlung 49 3.2 T h e method of collective coordinates 59 3.3 Quantization of the bremsstrahlung problem 78 Oscillations of a p i n n e d soliton 83 4.1 Classical treatment 83 4.2 Q u a n t u m treatment 91 4.3 M a p p i n g to the Caldeira-Leggett formalism 96 Discussion 103 5.1 Interpretation of the spectral density function J(u) 103 5.2 W h a t could be done next? 104 5.3 Conclusion 105 Bibliography 107 Appendices A Detailed derivation of equation (2.10) and (2.13) 110 B Inverse of the velocity matrix D 117 C Lifetime of an environment-coupled harmonic oscillator 130 iv List of Figures 1.1 P r o b a b i l i t y densities of low-lying quantum energy levels for a finite square (a) and harmonic (b) potential well 8 1.2 A metastable potential 9 1.3 A 180° domain wall, also called Bloch wall 14 1.4 T h e coordinate system 15 1.5 T h e physical pinning potential in the absence (a) and presence (b) of a magnetic force H 21 1.6 A soliton and anti-soliton at rest 23 1.7 A mechanical model for a soliton obeying a Sine-Gordon equation 1.8 T h e propagation of a kink 25 1.9 T w o kinks of same helicity 28 2.1 A general potential V(x) 38 3.1 T h e impurity wavefunction "moving" towards a stationary soliton . . . . 51 3.2 T h e integration path 56 4.1 T h e physical potential provided by the pinning center 86 4.2 Q u a n t u m and classical oscillations of a soliton 92 4.3 T h e normalized spectral density function J(u) 101 C.l A shifted harmonic potential in the presence of a constant force 128 C.2 T h e integration path Y x . . . . : R 24 for the partial evaluation of the second order perturbation energy term 135 v Acknowledgements First of a l l , I wish to thank my supervisor, D r . P h i l i p C . E . Stamp, for suggesting the topic for this thesis and guiding me through to its completion. To various degrees, I am also indebted to D r . D o u g B o n n , D r . Ian Affleck, D r . G o r d o n Semenoff, and M a r t i n Dube for helpful discussions and comments. It is a pleasure for me to express my deepest gratitude to my mother for providing me with a good education, for meeting my needs from the most basic to the most emotional w i t h a motherly love, and for her unfailing support i n regard to my studies and, most especially, my recovery from a past and recent affliction. T h i s small section could hardly contain the overflow of thankfulness I have for D r . B . C . E u , without whom I would simply never have persevered in this noble field of physics. In this regard, I also wish to mention D r . L . J . Lewis and D r . G . R . B r o w n . Lastly, I wish to thank w i t h a l l my heart O u r Blessed V i r g i n M a r y for her constant protection and intercession on my behalf w i t h her Son, O u r L o r d Jesus Christ; and I humbly pray that she may remember me at the hour of my death. Also, I a m most grateful for the undeserved help provided me by St. Francis of Assisi, St. Augustine, St. M o n i c a , Ste. Therese de Lisieux, St. M a r y Magdalene, my guardian angel, and my late maternal grandmother whom I loved very much. Benoit Leduc Feast of Saint John of the Cross, December 14th, 1995. vi T o Jeannie and D r . E u A p a r t from me you can do nothing. John 15:5 If any one shall not be ashamed to assert that, except for matter, nothing exists; let h i m be anathema. T h e V a t i c a n Council, Session 3, C a n o n 2, A p r i l 24th, 1870. viii Chapter 1 Introduction Ever since Descartes hinted at a world view based on a materialistic and mechanical reductionistic m o d e l , natural philosophers of the modern era have attempted to explain 1 complex phenomena as being nothing more than the sum of more fundamental selfexistent ones occurring at a microscopic level regulated by some laws, inferred or simply given axiomatically. In contrast w i t h most Scholastics of the Middle-Ages, for w h o m , Aristotelian as they were, four general kinds of causes were recognized — namely: 1) the material cause (the scholastic causa materialis), which provided the passive receptacle on which the causes act; 2) the formal cause (causa formalis), which contributed the essence, idea, or quality of the t h i n g concerned; 3) the motive force of efficient cause efficiens) (causa which was the external compulsion that bodies had to obey; and 4) the final cause (causa finalis) as the goal to which everything strove and which everything served — modern science would be concerned only w i t h the efficient causes as the ontological pattern of being and becoming. A l t h o u g h the sciences in the nineteenth and twentieth centuries came to recognize this restriction as inappropriate for the faithful description and explanation of natural phenomena — by introducing, for instance, the notions of A s is well known, Descartes's secret philosophy, expounded in Le Monde (1633), which he did not dare publish, was almost entirely materialistic; it was only his public philosophy that was rationalistic to the extent to which reason did not conflict with religious dogma [1, 2]. Descartes considered the living as machines which could entirely be described by matter and motions [2, 3]. Motion, as confirmed by Newton's laws, was for him - as for Spinoza, Bruno, and Leibniz - self-caused and self-moving. Indeed, Spinoza went even further by asserting categorically that substance was not only self-subsistent but also self-existent, self-moving, and self-caused. J 1 Chapter 1. Introduction 2 reciprocal interaction, fields, and statistical determined events — very few scientists d i d , 2 and still to this day, do support the existence of organizational levels whose behaviors are different from and inherently irreducible to the superposition of their constitutive elements and their accompanying laws. Materialistic reductionism is appealing to reason owing to its simplicity and its exclusion of supernumerary transcendental principles. Physicists were among the first to wholeheartedly sell their souls over to this model p e r m i t t i n g them to have control over phenomena which, until then, were clouded i n an awe-inspiring mist of mystery. Poets would therefore be relegated to singing the odes of love and leaving Nature to those new mechanists. A t the end of the nineteenth century, various results, such as the blackbody ultraviolet catastrophe, based on the established mathematical models of classical physics of the time, were i n outright conflict w i t h the experiments which they proposed to explain. These eventually paved the way for a new theory, essentially microscopic: quantum mechanics. A s this latter model dealt w i t h more fundamental units of matter, many physicists, imbued as they were w i t h reductionistic ideas, suggested that the classical world should entirely be explained w i t h i n its framework. However, that was without reckoning w i t h the powerful personality of Niels B o h r and his idealistic Copenhagen interpretation of quantum mechanics which quickly quelled the ambition of these reductionists. Although Newton did introduce those last two notions in Principia with his third law of motion and theory of gravitation, they were hardly ever recognized or understood by the scientific community until much later. (However, one must note that his notion of a "field" was different from that in which it will be understood in the nineteenth century.) The fact that the agent can be influenced by the so-called passive object on which it acts, which then acts on the latter in a new way so that the cycle repeats itself, was hardly present in the scientific thought of the seventeenth century. [1] 2 Chapter 1. Introduction 3 T h e q u a n t u m measurement paradox 1.1 T h i s interpretation, maintained by B o h r throughout his life, can be subsumed in four major points: [4] 1. Microscopic entities (such as electrons and atoms) are not even to be thought of as possessing properties in the absence of specification of the macroscopic experimental arrangement. 2. Macroscopic experimental arrangements, and the results of experiments, are to be expressed in classical, realistic terms. 3 . There exists what B o h r repeatedly refers to as an "unanalyzable link" between the microsystem and the macroscopic measuring apparatus. 4. T h e principle of complementarity: different experimental arrangements exclude one another, and the measurement of one property may therefore be "complementary" to the measurement of another. It can be inferred from his numerous papers [5, 6] that B o h r himself d i d not rule out the possibility of obtaining macroscopic quantum phenomena in which the principle of complementarity is retained and merely extended to a new level, the macroscopic one. 3 S t i l l , very few physicists hoped to see quantum phenomena at the macroscopic level owing to an argument developed by B o h r . F i r s t , he stressed the fact that quantum effects are important only when the action S is on the order of the quantum of action h. Second, i f E is the energy of the macroscopic system in question and r i ?a u c l d the order of magnitude of its classical period of motion, then in a typical experiment we have S^ETCI SO that the condition S implies E£hu i. c T h i r d , since for a macroscopic system the relevant classical frequencies are certainly not greater than 1 0 s , this i n 1 6 _ 1 Namely, that macrophysical objects have objective existence and intrinsic properties in one set of circumstances such as when used as a measuring apparatus, and have properties relative to the observer in another set of circumstances. 3 Chapter 1. Introduction 4 turn implies that the characteristic energy scale for the macroscopic motion can be no greater t h a n a few e V (the order of ionization energy of a single atom). B o h r concluded that it was simply unlikely that any macroscopic variable should be associated w i t h so tiny an energy. Since the time of Bohr, however, physics w i t h the help of technology has brought to our attention cases where the motion of a macroscopic microscopic variable can be controlled by a energy on the order of the thermal energy of an atom at room temperature or even lower (e.g., a macroscopic variable such as the current or trapped flux in a bulk superconducting ring). Meanwhile, the reductionistic compulsion resurfaced w i t h the publication of von Neumann's book [7] on the axiomatization of quantum mechanics and his interpretation of measurements, and its popularized adaptation by L o n d o n and Bauer [8] of the problem of the quantum measurement as exemplified by the Schrodinger's C a t Paradox. Let us briefly expose this paradox. For simplicity, consider an ensemble of microscopic systems which have only two available microscopically different states ipi and ip2- I n order to mea- sure which of these states is realized for a particular system of the ensemble, we couple it to a macroscopic measuring apparatus assumed to be entirely described by quantum mechanics. T h i s measuring device is such that, upon coupling, the state of the sys- tem will induce final state \&i of the apparatus, while ip2 will induce the macroscopically distinguishable final state of the apparatus $2- Schematically, we have ^1*0 -» ^ 1 * 1 and ip ^o 2 -r ^ 2 ^ 2 where ty corresponds to the initial state of the apparatus. 0 Now, suppose we prepare an ensemble for the microsystems which is described by the linear superposition state aipi + b T/> , a,b ^ 0. W h e n coupled to a measuring apparatus, 2 what is its description after the interaction? In order to find this out, we consider the Chapter 1. Introduction 5 quantum evolution of the "universe" consisting of the apparatus plus the microsystem which, by making good use of the linearity of quantum mechanics, is given by ( a + aipx * i 6^2)^0 + 6^2*2. (1.1) We come to the astounding conclusion that we obtain a linear superposition of states of the "universe" corresponding to macroscopically distinct behavior of the measuring apparatus. Various solutions have been proposed i n the past to account for the lack of such macroscopic linear superpositions. A m o n g the staunch supporters of the universality of quantum mechanics i n the physical world (who seem to form a majority in the physics community), it is argued that it w i l l be impossible i n real life to discriminate between the experimental predictions made by equation (1.1) and those made by a classical "mixture" description, in which the universe is simply assigned a probability pi = | Ci | of being in 2 the macroscopic state ipi^i, i = 1,2 (ci = a a n d c 2 = b). In the technical language of the quantum mechanical density m a t r i x p , it is claimed that its off-diagonal elements in a representation corresponding to the macroscopically distinguishable states i are unobservable in any realistic experiment so that the correct p, which has elements Pij i = c c ii may be safely replaced by the one corresponding to the classical mixture, namely d 0 = I 12 r T h i s reduction is claimed to be effected by the dissipative interaction that the so-called universe (which is never totally isolated) has w i t h its environment, so that possible quant u m macroscopic phenomena are washed out; using their jargon, the environment would decohere the wavefunction (1.1) w i t h macroscopic distinct states. However, as was so Chapter 1. Introduction 6 well pointed out by Leggett [9], the conceptual problem is not whether at the macroscopic level the object behaves as if it were i n a definite macroscopic state, but whether it is i n such a state. In other words, can we simply assume that the linear superposition (1.1) ever exist at the classical level? It seems evident from the foregoing discussion that a crucial test for the pretension of quantum mechanics as the all-encompassing physical law, extending its dominion right up to the classical world, is the observation of a macroscopic quantum coherence ( M Q C ) which presupposes the existence of such macroscopic states (1.1). U n t i l now, no such conclusive test has been shown. O n the other hand, quantum tunnelling of a macroscopic variable has been observed! [10, 31, 33, 34] Before discussing at length this promising vindication of quantum mechanics in the classical regime, we return to the decoherence effect of the environment discussed above. It was thought (on the implicit assumption of the universality of quantum mechanics) that this decoherence effect could be more or less severe so as to permit us to expect that, under some specific circumstances, a remnant of quantum macroscopic behavior could be made manifest. To that end, it seemed imperative to construct a quantum theoretical model which would take into account the dissipative effect of the environment; the Caldeira-Leggett formalism [12], to be introduced below, is a response to such a need. B y a greater understanding of the environment's effect on a system, one could certainly devise ways i n which the disastrous decoherence would be brought low enough to salvage a possible quantum character to a macroscopic object. T h i s study w i l l bear on one macroscopic quantum effect: the macroscopic quantum tunnelling, abbreviated as M Q T , i n the presence of dissipation. A s Leggett [11] rightly observed, i n most discussions of macroscopic quantum phenomena, the "system" in question is microscopic, and it is the the interaction w i t h a macroscopic system, the measuring apparatus, which destroys the coherence of its wave function. B y contrast, below, we Chapter 1. Introduction 7 consider a "system" which is itself macroscopic, the coherence of whose wavefunction may be destroyed by its dissipative interaction w i t h a host of microscopic degrees of freedom representing the "environment." Let us briefly review the elementary features of microscopic quantum tunnelling in the following section. Chapter 1.2 1. Introduction 8 T h e microscopic q u a n t u m tunnelling In figure (1.1(a)), we show the first energy levels of a particle in a finite square potential well of w i d t h a and arbitrary depth U, along w i t h their respective probability curves given by the modulus squared of the wavefunctions. (a) (b) Figure 1.1: P r o b a b i l i t y densities of low-lying quantum energy levels for a finite square (a) and harmonic (b) potential well. In (a), the probability densities for the lowest three energy levels of a square well of height U and width a are shown with a non-zero probability of finding a particle with energy less than U outside the well. The same is shown i n (b) for a harmonic potential well; the vertical bars indicate the classical turning points and the dashed curves represent the classical probability densities [14]. Only the lowest three probability densities are explicitly shown. Classically, as long as the energy E of the particle is less t h a n U, it is impossible for the particle to be present outside the potential well, the barriers on both sides preventing it from crossing beyond the classical turning points x = 0 and x = o. However, quantum mechanics allows the wavefunctions to leak through by a small amount into the classically forbidden zone w i t h i n which they adopt an exponentially decaying form. Therefore, by the orthodox interpretation of quantum mechanics, there should be a non-zero probability of finding the particle in that region. Figure (1.1(b)) shows the same for a quadratic potential well. Chapter 1. Introduction 9 T h i s is the origin of the a-decay i n nuclei; although nucleons are totally confined classically, some nucleons will make it through the potential barrier thanks to their wavefunctions extending beyond it. It is important to note that what distinguishes M Q T from the ar-particle decay and other microscopic tunnelling effects is that the classically accessible regions which are separated by the barrier are effectively macroscopically distinct. Now, consider a metastable potential well V(Q) as shown i n figure (1.2). T h e stability V(Q) 1 \ V F i g u r e 1.2: width B. A metastable potential. \B \ Q This metastable potential has a barrier of height U and Q frequency at the bottom of this well located at Q = 0 is given by 1 2 M for a particle of mass M. 1/2 dV dQ 21 ( ( L 2 ) T h e low-lying quantum levels associated w i t h this well w i t h energy less than U w i l l , also, see their wavefunctions leaking through beyond the point 0 Q = B. B y construction of V(Q), we see that, once the particle has tunnelled through the barrier to the right, it w i l l simply roll off and move away from point B w i t h a very small likelihood of tunnelling right back i n . Let the Lagrangian of this system be given by L(Q,Q) = -MQ* l V(Q) (1.3) Chapter 1. Introduction 10 B y the method of Wentzel-Kramers-Brillouin ( W K B ) , one can calculate the probability per unit time that the particle escape from the well. T h e result i s P = Ae' with B / h 4 A = CLo {B/2-Kh) / 1 2 0 (1.4) where B = 2 j * yJ2MV{Q) dQ (1.5) and C is a dimensionless constant of order unity depending on the shape of the potential V(Q) [12]. T h i s is the standard result for a particle, be it microscopic or macroscopic, subjected to a conservative potential V(Q). However, we noted above that the difficulty to observe quantum tunnelling at the macroscopic level had to do w i t h the dissipative effect of the environment on the macroscopic object. Therefore, one should modify the simple Lagrangian (1.3) i n order to take this into account. T h i s is discussed in the next section. 1.3 T h e Caldeira-Leggett formalism T h i s treatment w i l l be rather brief; we suggest the reader consult the proper references for more details [12,13]. In the Caldeira-Leggett paper [12], the authors sought to answer the questions of how and to what degree dissipation can affect the tunnelling probability of a particle whose quasi-classical equation of motion is given by MQ + rjQ + ^ = F (t), ext • (1.6) w i t h r] being the phenomenological friction coefficient which may be frequency- and amplitude-dependent [13], as compared w i t h MQ 4 + ^ The symbol = stands for "by definition." = F (t), ext (1.7) Chapter 1. Introduction 11 the quasi-classical equation of motion derived from the dissipation-free Lagrangian (1.3) (complemented w i t h the required term QF (t) ext which accounts for the forced motion). We stress the fact that both potentials V(Q) i n equations (1.6) and (1.7) are exactly the same regardless of possible potential-renormalization effects [13]. 5 B y the term quasi- classical, we mean that the above equations of motion apply for the expectation value Q = (Q) where we assume that it is legitimate to replace ( ) dQ by ~ ' d V { Q ) dQ J Such a substitution is valid whenever the potential V(Q) Q=(Q) varies appreciably only over macroscopic scales and the wavefunction associated w i t h the system has negligible uncertainties or fluctuations about its expectation values of the coordinate Q , {Q), and its conjugate momentum P , ( P ) . Caldeira and Leggett showed that the most general Lagrangian we need ever be concerned w i t h i n modeling the interaction between a system and its environment is L = \Mt? - J V(Q) + QF (t) + \ £ ext 3 m j {x^-rfx)) (1.8) ^i 2m under the. express assumption of three provisos, namely, 1) any one degree of freedom of the environment, X j , is sufficiently weakly perturbed so as to neglect nonlinear effects [15], 2) the interaction Lagrangian coupling the system to its environment contains terms either (i) linear in the system coordinate Q and its conjugate momentum P or (ii) quadratic in Q and P , but not containing the environment's coordinates Xj or their conjugate momenta, and 3) condition of time-reversal invariance. W e note that M , V(Q), m,j, LUJ, and T h e dissipative interaction with the environment may cause the potential V(Q) appearing in (1.3) to be different from the one in (1.6). For purposes of comparison, we suppose, then, that the potential V(Q) in (1.3) already includes this possible interaction and enquire about the sheer effect of the term 7] Q on the tunnelling rate. 5 Chapter 1. Introduction 12 the environment coordinates Xj i n (1.8) may carry renormalization factors arising from the interaction itself. T h e stability frequencies Uj are those associated w i t h the small oscillations of the environment coordinates Xj consequent to their small perturbation as required under proviso 1 above. In other words, we can represent a system's environment, insofar as its effect on the system is concerned, by a bath of harmonic oscillators. T h i s idea is not new, already Feynman and Vernon in 1963 [15] came to the same conclusion. In reference [12], only proviso 1 is explicitly retained plus some other restrictions so that the general Lagrangian takes the form L = \MQ - 2 V(Q) + QF (t) + 1 £ ext mj (x* - u>? x*) (1.9) Here, we note that the interaction is more general than that i n (1.8) w i t h a coupling linear i n the environment coordinates Xj only. T w o cases are to distinguished: the quasilinear and the the strictly linear dissipation mechanisms. T h e former applies when, for every fixed typical amplitude Q, it is possible to find a frequency co(Q) small enough such that equation (1.6) holds, and the latter when, for any typical amplitude Q, it is possible to find such a frequency u independent of Q. It is then shown in reference [12] that, in order that the Lagrangian (1.9) be consistent w i t h the quasi-classical damped equation of motion for the system's coordinate Q i n the quasi-linear case, we should have the condition — 77 1 2 y l^-(Q) m u) d Q 1 2 , 8{u- ). For the more specific case of strict linearity, we should have Fj(Q) condition (1.10) reduces to v = (1.10) Uj 3 f E - ^ I Q I 2 ^ - ^ - ) , — CjQ so that Chapter 1. Introduction 13 or, cast in a different form, to J(UJ) = rju (1.11) with z . rrij ujj denoted as the spectral density function. T h e condition of strict linearity, therefore, requires that the coupling be linear in both the system's and environment coordinates, Q and{xj}. Equipped w i t h these theoretical tools and making good use of the instanton technique, Caldeira and Leggett can calculate the quantum tunnelling rate at temperature T = 0 of a system coupled to a bath of harmonic oscillators and deduce an expression of the form (1.4) w i t h the various parameters carrying an integration of the environment's effect. However, i n this study, we are not directly concerned w i t h the calculation of the quantum tunnelling rate. Instead, we seek to calculate a piece of information which gives an idea of the dissipation involved in a ferromagnetic system at T = 0 and is used for the evaluation of the M Q T rate itself: the spectral density function J(w). It is time, now, to introduce our macroscopic object that w i l l be associated w i t h the macroscopic coordinate Q; to this, the next section is devoted. Chapter 1.4 1. Introduction 14 T h e domain wall and its dynamics In this study, our macroscopic objects w i l l be of a magnetic nature: domain walls in magnetic wires [16, 17]. A domain wall is the transition region of finite length separating magnetic domains w i t h different orientations of the magnetization. Here, we w i l l focus on a wall separating domains w i t h a 180° reversal of magnetization called a B l o c h wall; a sketch of which is shown i n figure (1.3) [18]. Figure 1.3: A 180° domain wall, also called B l o c h wall. Consider the coordinate system for a 180° domain wall shown i n figure (1.4) w i t h the magnetization vector pointing i n the ^-direction for z —» —oo and in the negative ^-direction for z —» +oo. We suppose that in any x-y plane, a l l magnetization vectors have the same relative orientation so that the magnetization varies only along the z-axis. T h e isotropic exchange energy for nearest-neighbour interaction is given by (1.13) where Sj is the t o t a l magnetization vector, or spin, of the i ih netic substances. ion and J > 0 for ferromag- B y construction of our domain wall, it is more useful to consider the energy per unit area e . ex If we assume for simplicity that we deal w i t h a simple hexagonal Chapter 1. Introduction 15 M n z M F i g u r e 1.4: T h e coordinate system. The angles 6 and <j> are, respectively, the polar and azimuthal coordinates of the magnetization vector M . For a Bloch wall, the polar angle has value 7r/2 for each magnetization vector along the z-axis and the azimuthal angle takes the value of zero at z -> - o o and of TT at z —¥ +oo. lattice whose primitive vectors are of equal length a, then J s *s n (1.14) % w i t h the index n running i n one dimension only, the z-axis. Expression (1.14) can be rewritten as JS 2 ex = coso; e (1.15) n where 5 is the magnitude of the magnetization, i.e., the total spin quantum number of each ion, and a n the angle subtending'the magnetization vectors S a n d S „ n + 1 . For small a , we can expand c o s a to quadratic order to obtain n n £-P.T. S Yl l V3a J a 92 + constant. (1.16) n In the continuum l i m i t , the angle between neighbouring magnetization vectors is given by 'df dz, Chapter 1. Introduction 16 w i t h 4> being the azimuthal angle because the spins are totally confined to x-y planes. T h e n , equation (1.16) transforms as JS /•+«> 2 (d^ a i-oo ^ \dz) Z 2 ' (1-1^) where we have discarded the arbitrary constant and scale z i n units of a. In chapter 2, we w i l l consider a model for ferromagnetic wires w i t h an anisotropy energy t e r m K y e ra 2 n which makes the x-axis the easy axis. A g a i n , i n the continuum l i m i t , this anisotropy t e r m per unit area for our previous simple hexagonal lattice becomes KS r ° 2 £ a = +0 / dz sin <f>(z) 2 (1.18) for spin vectors confined to x-y planes so that 6 = n/2 for each of them. T h e total energy per unit area is thus In order to determine the stable wall configuration, we must minimize the energy expression e by the usual variational procedure which yields the following Euler equation dz 2 ' Ja 2 smn(z) (1.20) where n = 2 (j). Expression (1.20) is the static Sine-Gordon equation. W e w i l l have more to say about the Sine-Gordon (SG) equation i n the next section; for the time being, suffice it to say that the non-trivial solutions to (1.20) are r){z) m 2 = A = 4 t a n " (e ( )^j 1 K J a ±m z with (1.21) Chapter 1. Introduction 17 T h e ± signs i n (1.21) correspond to the two possible helicities of the rotation of the magnetization vector S w i t h i n the wall (cf. sec. 4). Also, it is supposed i n (1.21) that the spin vector S points exactly i n the positive y-direction at z = 0. 6 For simplicity, we shall choose the helicity corresponding to the positive sign. In terms of the original azimuthal angle </>, we have from (1.21) that sin^(z) = sech(mz). (1.22) T h e energy of the wall per unit area, equation (1.19), can also be rewritten i n terms of the angle n as e = ^£ £> {l ( I ) 2 + ™ 2 ( i - « • " > } • < L 2 3 > so that upon using (1.21) we find 4 JS m 2 y/3a A J S 2 V K J v^a (1.24) 2 D o t i n g [19] i n 1948 was the first to discover that a domain wall exhibits an inertia, despite the lack of any mass displacement. T h e mass of a domain wall has its origin i n the angular momenta of the spins forming the wall. We recall that a spin w i l l adopt a precession motion about the direction of the applied magnetic field i n the same fashion that a spinning top would when placed in a gravitational field. Hence, when we apply a magnetic field in the x-direction, each spin of the wall w i l l start precessing about the x-axis resulting i n a rotation out of the x-y plane. T h i s rotation w i l l simply induce the appearance of magnetic free poles giving rise to a small demagnetization field in the z-direction, H z = 1.25 Ho Should the spin vector point in the positive y-direction at z = z , the exponential in (1.21) would have as its argument ± m ( z — z ) instead. 6 0 0 Chapter where M 1. z Introduction 18 is the z-component of the induced magnetization and fj, the permeability of 0 the vacuum. T h i s demagnetization field then, i n t u r n , acts on each spin so as to induce a precession motion i n the x-y plane which results in the displacement of the wall i n the z-direction. T h e rotational velocity of this precession motion about the z-axis is given by <j> = gH (1.26) z where g is the gyromagnetic constant such that g = g //{, w i t h //& being the B o h r magneton and g the gyromagnetic ratio. B u t , dt-~ dz dt - V dz [ where v is the translational velocity of the wall as it moves along the z-axis. n Using equations (1.26) and (1.27), we find M. = --^£-V ."• 9 (1.28) oz Thus, the m o v i n g wall w i l l have the additional surface energy 1 r+°° *k = —z / • dz Z J—oo = ^r 2 fi 0 i z •„» M MH Z • l Z . ' (i.29) = i-oo z (^T~) = 2g 2 7-oo \dz) y 1 F r o m (1.22), we find that m sech (mz), 2 2 so that substituting the latter into (1.29), e& becomes Ho™ e k i = , (1.30) M , Chapter 1. Introduction 19 where can be interpreted as a v i r t u a l mass per unit area for the wall. W h a t is the energy supplied by the applied magnetic field H x which makes the wall move? For low fields and velocities, assume that the wall has moved by a distance Q, in units of the lattice constant a, from its original center's position at z = 0. T h e work done per unit area by the field H x W is given symbolically by = A Energy = —H x • j dz { W a l l configuration &tz = Q — W a l l configuration a t z = 0 }. O n each spin vector, we have the following change in direction F r o m (1.22), we deduce for a wall that cos 0(2) = -tgh(mz). (1.32) T h u s , when the wall's center has moved to position z = Q , we have cos 4>f(z) = w where 4>f(z) and 4>i{ ) z a r e — tgh (m(z — Q)) cos(f>i(z) + mQ sech (mz), (1.33) 2 > respectively, the azimuthal angle for the i n i t i a l and configurations, and the last approximation is valid for small displacements. final Since, for Chapter 1. Introduction 20 each spin, the work done by the field is given by -H S j(z) x x + H S (z) x Xyi = -H S (cos <f>f(z) - f» —H S mQ x x cos<fo(z)) sech (mz), 2 we obtain by integration +O0 / = dz sech 2 (mz) -oo -2H SQ (1.34) X as the energy per unit area supplied by the magnetic field H x i n order to push the wall by a small distance Q. Therefore, AEnergy 2H x Q S A P can be thought of as the pressure P exerted on the wall by the applied field. In general, this pressure w i l l cause the wall to swell and a radius of curvature is therefore established. In quasi-one-dimensional ferromagnetic systems, we assume that the area is small enough to neglect this curvature effect. A s mentioned previously, D o r i n g showed that the domain wall has dynamic properties analogous to a particle of mass M. Indeed, one can write down a phenomenological equation of motion for a wall as MQ(t) + r,Q(t) + = P , ext (1.35) where Q(t) is the domain wall's center coordinate, M the v i r t u a l mass of the wall given in (1.31), r] the friction coefficient, V(Q) a conservative surface potential, and, finally, P ext an external forcing pressure term. Note that equation (1.35) is given i n terms of quantities Chapter 1. Introduction 21 per unit area. In general, the magnetic material will contain impurities, defects, or voids; all of which contribute to a non-uniform background for the domain wall in addition to the magnetic anisotropy. T h i s means that the wall will see its surface energy vary as it moves from one position to the other. In this study, we w i l l show in chapter 4 that pinning centers can be thought of, qualitatively, as giving rise to a potential of the form V(Q) = 2o!sech(m<2) + e (1.36) where e is the surface energy of the wall (1.24) and a a negative small coupling constant. A d d i n g the potential energy arising from the application of the external magnetic field H , equation (1.34), we obtain, by discarding the constant term e, x V(Q) = 2asech(mQ) - 2H SQ. (1.37) X B o t h potentials are shown i n figure (1.5). i v * Q 1— B (a) (b) Figure 1.5: T h e physical pinning potential i n the absence (a) and presence (b) of a magnetic force H . x A pinning impurity in a quasi-one-dimensional ferromagnet generates the above potentials for a Bloch domain wall whose center is parametrized by coordinate Q. In (a), the potential is of a hyperbolic secant form. The application of a magnetic field H x in (b) effectively tips potential (a) asymmetrically and gives rise to a metastable potential with a barrier of height U and width B. 0 Compare with fig. (1.2). Chapter 1. Introduction 22 In this case of an applied magnetic field H , x MQ(t) + nQ(t) equation (1.35) reads + = 2SH . (1.38) X One can immediately notice the parallel between the equation of motion for the wall's coordinate Q(t) (1.35) and figure (1.5(b)) and the quasi-classical equation of motion (1.6) and figure (1.2) in section 2. In other words, Q(t) becomes the coordinate for the macroscopic object we are interested i n : the domain w a l l . T h e phenomenological friction coefficient n carries the dissipative effect of the environment on the motion of a domain wall. One would then be curious to know how this dissipative interaction affects, at temperature T = 0, the tunnelling of a domain wall under a barrier of height U 0 as shown i n figure (1.5(b)). T h i s is exactly the question the Caldeira-Leggett formalism in section 3 was tailored to answer. T h i s study w i l l be mainly concerned w i t h one source of dissipation, namely, the magnons. We w i l l discuss in the last section the other sources of dissipation likely to be present at T = 0 and T ^ 0. B u t first, we w i l l t u r n to the Sine-Gordon equation and some of its properties. 1.5 T h e S i n e - G o r d o n equation The Sine-Gordon equation [20, 21] a non-linear wave-equation of the form d ib(x,t) \ <,d ib(xA) ~ <Z 2 2 d 2 ) 2 dx . ,, ; + "o s m V M ) = 9 n 0 . (1.39) where c is a characteristic velocity and to a characteristic frequency. T h i s is the equation 0 0 of motion obtained from the Lagrangian density £ = \ i ^ 2 - c ^ } ,2 0 - u; (1 - cos i>). 2 (1.40) We note that i n the case of very small amplitude | ip | - C 1, one can replace sin-0 by tp so that the Sine-Gordon equation reduces to the ordinary K l e i n - G o r d o n equation d ib 2 d ib 2 Chapter 1. Introduction 23 It has been shown by Faddeev and Takhtadzhian [22] that the classical Sine-Gordon system represents a completely integrable H a m i l t o n i a n system whose spectrum is totally exhausted by free oscillation modes (or magnons), solitons, and breathers. T h e last two represent bounded large-amplitude solutions of the S G equation (1.39) for which the soliton has the form r (x,t) = 4tgh- {exp[±^ (x-^)]}, (1.41) (l - (1.42) 1 s 7 where 7 = ^ ) \v\< C o , and the ± signs refer respectively to a soliton and anti-soliton Of opposite helicities. In the literature, one also calls a soliton a kink. A soliton and an anti-soliton i n their rest frames are shown i n figure (1.6). -10 10 0 20 30 40 X Figure 1.6: A soliton and anti-soliton at rest. These are bounded since ip] —> 0 (mod 2-K) as | x \—> oo. Furthermore, these solutions are Lorentz invariant, as are the Lagrangian density and its resulting equation of motion ((1.39), (1.40)) where c plays the role of a l i m i t i n g speed. 7 0 M o v i n g i n the rest frame of a soliton w i t h speed v, one obtains il> {x,t) = a 7 4 tgh- j e ^ / ^ J 1 . In magnetic systems, this limiting value is simply the Walker velocity. (1.43) Chapter 1. Introduction 24 Before finding the excitations of the field ip(x,t) in the presence of a soliton tp (x,t), s it might be useful to digress for a brief moment and offer a mechanical analogy whose equation of motion is of the Sine-Gordon form. A mechanical model is shown below in figure (1.7). , i i , i , 7, ,fi.Oi.in mm- -g) llilllllllllllllllllllllllllllllll Mg. 3. - Mechanical model of the SUE: a) spring 0.2 diameter, b) solder, c) brass, d) tap and thread, e) piano wire, /) nail, g) and h) ball bearings, i) wooden base. Figure 1.7: A mechanical model for a soliton obeying a Sine-Gordon equation [21]. Its motion is governed by the difference differential equation M #4>i dt 2 K [<p i+1 - 2(f>i + & _ i ] - T sin 4 , where fa is the angle between the direction of the gravitational field and the i M (1.44) t h pendulum, the moment of inertia of a single pendulum, K the constant torque of a section of spring between two pendula, and T s i n ^ j the gravitational restoring torque of the i ih pendulum. For waves which vary slightly over the distance Ax separating two pendula, one can t u r n the difference equation (1.44) into a scaled continuum partial differential Chapter 1. Introduction 25 equation d (j> 2 dx 2 + s'mcj) = where distance is measured i n units of 7„ = yK/T 0, (1.45) Ax and time i n units of t = 0 yM/T. Expression (1.45) is simply the normalized S G equation (1.39). A strobe photograph of a soliton travelling from right to left is shown i n figure (1.8). F i g u r e 1.8: T h e propagation of a kink. A single kink is moving from right to left for the mechanical model shown in figure (1.7). The strobe photographs are taken at time intervals of 0.6 second. The kink slows down because of friction effects which are not accounted for in the equation of motion (1.45). Notice that the kink becomes wider as its velocity decreases, which demonstrates the "Lorentz" contraction effect we expect from the "Lorentz" invariant continuum equation of motion (1.46). The limiting velocity is 50 crh/s for this mechanical system as the length and time scales, 7 and 0 r , are 5 cm and 0.1 s respectively [21]. (cf. text) 0 Now, back to our Sine-Gordon Lagrangian. It is well known that, in condensed-matter physics, many properties of materials are explained simply by the notion that s m a l l amplitude waves (phonons, magnons, etc . . . ) are the fundamental entities that enter the description of their thermodynamics and response to external probes. One, therefore, Chapter 1. Introduction 26 would like to find the excitations of the field ib about a stationary soliton solution ib in s the form of small amplitudes 4>. Thus, consider ib{x,t) = ib {x) + <P{x,t) (1.46) s with | <f>(x,t) | < 1. We seek to find a solution for (j)(x,t) such that ip(x,t) (1.47) obeys the Sine-Gordon equa- tion (1.39). Substituting (1.46) i n (1.39) and making use of the expansion sm(t/;(x,t)) = s'm(ip (x)) + cos(tp (x)) <j>{x, t), s s (1-48) which is sufficient on account of (1.47), we obtain '<j> - c\<b" + u> (l - 2 2sech (^V)) 2 0 = 0. (1.49) In deriving (1.49), we used the identity cos(ibs(x)) = 1 — 2 s e c h ( — xj 2 \ c 0 J and equation (1.39) for the static solution ib {x). a In order to diagonalize equation (1.49), we seek its eigenfunctions by assuming a solution (j)(x, t) w i t h harmonic time-dependence <f>(x,t) = f(x)e-^, (1.50) and substitute in (1.49) to obtain the following equation for , —c: 0 d f(x) f(x): 2 dx 2 + u; (l 2 0 2sech ^x)) 2 f(x) = co f(x). 2 T h e above equation has the form of a Schrodinger's equation w i t h potential V(x) = u (l 0 - 2sech (^V)). 2 (1.51) Chapter 1. Introduction T h e solutions f(x) 27 are combinations of hypergeometric functions [23] which generate exactly one bound state w i t h "2 = o and eigenfunction f (x) = 2— sech (—x) , Cg \C / and oscillatory states w i t h a continuum of eigenvalues (1.52) b 0 co? = ell 2 + u? (1.53) 0 and corresponding eigenfunctions T h e zero-frequency eigenmode fb(x) corresponds roughly, as we w i l l see i n more details in chapter 3, to the translational excitation of the soliton in the limit of small displacements; whereas fi(x) are the small-amplitude oscillatory excitation waves for a magnetic system (the magnons) but in the P R E S E N C E of the stationary soliton and N O T built around the ground-state of a ferromagnetic material, in which a l l spins are aligned pointing in the direction of the applied field. For this latter case, we recover the usual magnons w i t h which we are a l l familiar. T h e different magnons, represented by fi(x), correspond to precessional modes in the domain wall, or magnetic soliton, configuration. Because the functions {fb(x), fi(x)} are eigenmodes of the self-adjoint operator they form a complete basis-set of orthogonal functions w i t h (fb, fb) = 8 ^ , (1.55) (fiJm) = 5(1-m), (1.56) (fiJb) = 0, (1.57) Chapter 1. Introduction 28 where > </,<?> = J—oo dx f{x)* g(x). A s a final note on the Sine-Gordon equation, we mention that there exist solutions t/)(x, t) which, under some circumstances, can be thought of as the linear superposition of various solitons and anti-solitons, and are classified by a topological index N defined as N = -L<ty(oo,t) - V(-oo,<)}. (1.58) For instance, two N = 1 solitons are shown i n figure (1.9) so that the t o t a l index NT is s given by 2N 8 = 2. q x Figure 1.9: T w o kinks of same helicity. Solitons fa and fa have the same topological index N = 1 so that the total function fa(x) = fa(x) + fa(x - q) has a total index of N = 2 N = 2. s t s In the mechanical analogy of the pendula, this corresponds to two successive revolutions of the pendula separated by a distance Q = Q\ + Qi so that the total angular displacement adds up to 47T. Chapter 1.6 1. Introduction 29 Sources of dissipation In addition to the dissipative effect of magnons on the quantum motion of a domain wall, this latter can also interact w i t h phonons via magnetoelastic coupling, photons, other impurities and defects from a distance, nuclear spins via the hyperfine coupling, and free electrons. F o r magnetic insulators, we may disregard this last dissipative interaction. A t temperature T = 0 , only the nuclear spin interaction should survive w i t h an appreciable effect in the absence of impurities. A s pointed out by Stamp [24], magnons and phonons should have no dissipative effect on the wall motion at T = 0 i n the absence of such defects. Besides, several workers, such as Chudnosky and Stamp [24, 26], have shown that the coupling to photons and other impurities gives negligible dissipation on the motion of the domain wall. One is thus left w i t h the interactions of phonons, magnons, and nuclear spins as the major sources of dissipation i n the presence of impurities. T h e case of phonons was first addressed by G a r g and K i m [27, 28] who showed that the dissipation might be weak for some single-domain a superohmic spectral density J{u) ~ ferromagnetic materials on account of to . 3 A t this point, we may digress slightly to acquaint the reader w i t h the various types of spectral density functions. For to —>• 0, the behavior of J (to) can be classified as sub-Ohmic, O h m i c , or super-Ohmic according as the exponent s in J(w) ~ to s (1.59) is s < 1, s = 1, o r s > 1. T h i s classification is coarse and needs refinement, but does serve our purpose for the time being. T h e reader may consult Leggett et al. [29] for more details. A s a rule of thumb, the dissipation can be considered small for the super-Ohmic case and large for the sub-Ohmic case. 8 T h e normal O h m i c case, s = 1, strikes a middle position between these two extremes and corresponds to the usual decomposition of the 8 Actually, the sub-Ohmic case is pathological and leads to serious problems of divergence. Chapter 1. Introduction 30 spectral density as J(CJ) = r/u), where n is the phenomenological friction coefficient in expression (1.6). Now, the application of phonon-related dissipation in the presence of impurities on a domain wall motion may require more care and, according to the few studies which have been made so far [24, 26], may give rise to an O h m i c dissipation. T h e dissipation caused by a spin bath, such as nuclear spins and magnetic impurities, was studied by Stamp and Prokof'ev [30] where a model different from the Caldeira-Leggett type must be considered w i t h the consequence that such a bath can have disastrous dissipative effects on macroscopic quantum phenomena. T h e interaction w i t h magnons in the presence of pinning centers is the object of this study. To be specific, our goal is to resolve the dissipative effect of magnons at T = 0 on the quantum motion of a domain wall bound to a pinning center i n the l i m i t of a one-dimensional ferromagnetic system such as a magnetic wire. W e w i l l have the domain wall oscillate back and forth at the bottom of the potential well V(Q) in (1.36), generated by the pinning defect, and extract the spectral density J(co) which gives a qualitative measure of the dissipation to be expected from magnons. T h i s whole calculation is motivated by the fact that some experiments might provide the evidence for domain wall tunnelling; the reader is urged to consult the works of Uehara and B a r b a r a et al. [31, 32, 33, 34] T h i s thesis w i l l be divided as follows: i n chapter 2, we justify the use of solitons and their attendant Sine-Gordon equation for the study of domain walls i n a 1-D ferromagnet; in chapter 3, we review in detail the effect of a soliton colliding w i t h some form of potential barrier and ,thus, prepare the stage for the manner w i t h which we w i l l deal w i t h a soliton bound to a pinning center in chapter 4; finally, i n chapter 5, a discussion of the results found in the previous chapter is given and possible refinements are proposed. Chapter 2 A theoretical model for ferromagnetic wires w i t h pinning centers In the first and t h i r d sections of this chapter, we give a general theoretical model for ferromagnetic wires w i t h impurities which display a soliton-like structure, the B l o c h wall. In the second section, we introduce the scheme of quantization about classical solutions which w i l l be further developed in the next chapter. 2.1 A model for some types of easy-axis ferromagnetic quasi-one-dimensional magnets We use a one dimensional spin-chain having the following H a m i l t o n i a n H = -J £ S •S n n+1 +K, £ n where K y <C K . n (S*) 2 + K Y, (SZ) 2 V (2-1) n T h i s H a m i l t o n i a n corresponds to a ferromagnet w i t h isotropic exchange z supplemented by a form of uniaxial anisotropy in which the symmetry in the x-y plane is broken by the weak anisotropy constant K y i n the y-direction. We w i l l show below that such a H a m i l t o n i a n i n the continuum and classical limit admits soliton solutions for which our spins are virtually confined to the x-y plane, this confinement being the result of the relatively high K . z The soliton solution for the time-independent case is the classical B l o c h W a l l configuration separating domains w i t h a 180° reversal of magnetization, (cf. figure (1.3)) We note that equation (2.1) is one-dimensional whereas magnetic wires are obviously 3-dimensional objects. However, for domain walls, we can always assume the behavior 31 Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 32 w i t h i n the two-dimensional cross-section of the wire as uniform for the following reason. In (1.21), for a planar domain wall w i t h the same anisotropy K S 2 y s i n 4> i n the y2 direction as above, m~ gives a measure of the w i d t h of the wall i n units of the lattice l constant a: We assume, then, that the cross-section of a magnetic wire is much less than a K /J so 2 y that we can, w i t h confidence, consider the magnetization as uniform across any crosssectional plane of the wire. Indeed, ignoring the demagnetization field, this rules out any wall configuration i n directions other than the z-axis and any curvature of the domain wall. In order to derive the equations of motion from the H a m i l t o n i a n (2.1) i n the continuum and classical l i m i t (i.e., loosely speaking, treating the magnetization S as a field and letting h go to zero), two methods are at our disposal : either we apply the continuum approximation to the equations of motion on the discrete lattice, or we work w i t h the continuum limit of the H a m i l t o n i a n (2.1) using the Poisson brackets derived from the quantum commutation relations. We shall choose the first method [35]. In the Heisenberg picture of quantum mechanics, the operators carry the full burden of the dynamical evolution of the system's observables. W e define the following decomposition of the magnetization vector S : S n = y/S(S + l) ( s i n 0 „ ( i ) cos0 (<), s i n 0 ( t ) s i n 0 ( i ) , cos 0 (t)) n n n n = (S*, S » , S ) (2.2) z n where {0 , 0 „ } are the polar and azimuthal angles of the magnetization vector for the n th n atom i n the spin chain. We obviously have the usual operations on an eigenvector \S] M ) n such as 15, M ) = hM \S, n n M ) = h ^S(S n + 1) cos/9 15, M ) (the angles {6 ,<j> } are n n n n quantized by virtue of the respective quantization of the magnetization components) as Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 33 well as the general commutation relationships for angular momentum observables such ] — ih as [ S^, S% 5 n > m . T h e H a m i l t o n i a n (2.1) can be rewritten as H=-JY,\Sn- S z n+l + l(S: • S- + S~ +l +1 • St)} + K Y,{S ? + K E ( ^ ) , (2.3) z z 2 n y the time evolution of the S^(t) being given by ij ~ ~_ ~_ ~_ - - - - - W r i t i n g the above using our decomposition of the magnetization vector, we obtain: < J 5 ^ + 1 > j s i n ^ s m ^ c * ^ - A - i ) -sm0 m + 1 sin0 e^™+ m 1 ~ (2.4) - sinti sin0 _ e~ (^ ~ ^ - 0 i l + sin § + 1) K y cos6 m sin t ^ e * ^ " * - m + 1 — 2 S ( S + 1) K s i n 6 cos$ m y n 1 m m ~ 1 \ s i n 9 sin 0 . m m B y the Ehrenfest theorem, we know that the time evolution of the expectation value of an observable is given by d(Cl(t)) _ *dt -i h <[«(*),#]>• (2-5) Since we are interested by the classical limit, we may always consider a quantum state | ^ ) which has the characteristics we would ascribe to a classical one, chiefly, that of being able to specify conjugate dynamical coordinates within the experimental uncertainty Such a state can always be constructed i n principle by the use of gaussian-like wavefunctions Chapter 2. A theoretical model for ferromagnetic wires with pinning centers which saturate the Heisenberg uncertainty inequality. For instance, for a free particle, the wavefunction / 1 \ V = [^) Mx) ip x 4 -(x-Xp) 2 0 e n e 2A* (2.6) w i l l give uncertainties A X = A and A F = ft/A, both of which can be made small i n the classical regime by choosing A appropriately and consonant w i t h the experimentally detectable range. For instance, by choosing A = 1 0 - 1 3 cm,which is the size of a proton, we obtain A P = 1 C T g cm/sec which translates for a particle of mass 1 g to a mere 14 uncertainty i n speed Av of about 1 0 ~ 14 cm/sec\ Hence, i n the classical scale, such a state can always be said to have well defined values of X and P, namely, XQ and po, since the fluctuations around these values are truly negligible. In the same manner, we assume that we have a state \$f )c such that we can ascribe to it angular momentum w i t h precise s components i n each directions i n the sense understood above. T h e n , the angles <$>} w i l l have definite values w i t h small fluctuations which, i n any case, w i l l be negligible as we let h —» 0 to m i m i c the classical behavior. Therefore, w i t h impunity, we w i l l consider such a state and sandwich (2.5) w i t h i t leading t o our particular form of the Ehrenfest theorem : 1 h^/S(S+ l)j t h * 2 J S cos e m = ^ + 1) { s i n 0 s i n 0 _ * t o » m - sin0 ra sin6 e ^ + i m + 1 (2.7) lC m m 1 ~ ^ - sin0 s i n 0 _ m + sin 9 2 z i e ^™ ~ sin 0 e~^ + m m+1 + ih \JS(S + 1) K cos 9 m l ^-i) _ m <M j — 2h S(S + 1) K sin 9 cos (j> sin 9 sin <j> . 2 m y m m m m We want a solution where the magnetization vector is nearly confined to the x-y plane; This state is l i m ® I '0s)c,m where | i> )c,m is the nearly classical state constructed above for the lattice point m. Also, note that we rule out possible linear combinations of such general state vectors since we know, by experience, that such combinations for the classical realm have never been observed. lr s 34 Chapter 2. A theoretical model for ferromagnetic wires with pinning so we expand the functions of 9 around the value $ ' 9, Hence, we o b t a i n m the deviations Q* = 0' — 9 . m m m m cos0 ~ m 9 M centers 35 = TV/2 and consider, instead of 2 (2.8) s m and sin0 ~ m 1 - (9 ) . s (2-9) 2 m Taking the discrete result over to the continuum case involves the use of expansions about the "field points" 6(z,i) and (f>(z,t) and the assumption that those fields are weakly varying over z and t. The details are left in appendix A and we merely state the result reminding the reader that we disregard any small variations of order 2 or higher. n ^ ^ y S J ^ A 2 h = Sfs + l j * 4 ^ W + ift ^S(S + 1) K 2 . ( 2 9 s z - 2h S{S + l)K 2 y 1 0 ) cos <b sin </>. T h e classical limit is obtained by having % —¥ 0 and S —> oo such that the product HS remains finite and gives the classical magnetization s in the l i m i t i n g process. Thus, by taking the l i m i t on b o t h sides of equation (2.10), we obtain d9 (z,t) ±-L-L dt s d (j>(z,t) a gJ V ' dz • - 2 2 2 = T ; n „ 2K . cos 0 sin <b, ± v ± 2 / n „ . (2.11 y J the second term disappearing since l i m h JS(S + 1) = 2 S-tOO .lint h S = s l i m ft = 0 . 2 S—HX) S - t o o T h e value a is the lattice spacing of our near-dense spin chain. Let n = 2 4>. T h e n , the above equation reduces to 89 (z,t) _ s dt a Js 2 2 d n(z,t) 2 dz 2 K y s sin 77(2,*). (2.12) We could also have made this expansion for the operators cos# and sin# resulting, respectively, in 6" and 1 - (6 ) where cos# is the operator measuring the deviation from the x-y plane. 2 TO s m 2 m m TO Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 36 pp In like manner, we can proceed w i t h the time evolution of operator S*(t) and obtain the following equation in the continuum and classical limit: „ J . ' . * ^ 2K..».. + (2,3) In this derivation, which is left in appendix A , it is important to note that we explicitly made use of the fact that K y <C K. z Differentiating equation (2.13) with respect to time, discarding the resulting term dtdz 2 as being negligible, and using equation (2.12), we get Rearranging the terms and dividing through by a s K 2 2 z d r]{z,t) —^-2 1 d r){z,t) ~2 2 2 where c = as {2K J) l , l 2 z 2 dt = 2 m J leads to . smr}(z,t) (2.15) a characteristic velocity, and m — (2K /a J) l 2 playing the role l 2 y of a mass. T h e above equation is none other than the Sine-Gordon equation we met i n chapter 1! Hence, from the requirement that the magnetization vector be confined to the x-y plane, we have shown that the field-like classical limit of our ferromagnetic Heisenberg spin chain admits a soliton solution. Incidentally, the Poisson bracket derived from the quantum commutator relation yJS(S + \) [cos<?„(<), (f> (t)] = -ih5 5 >, m nm u or using the small fluctuations approximation about the x-y plane ^/S(S + l)[e (t)J (t')] s n is s {On{t), <t> {t')} = -5 8 > m nm m = -ih5 6 ,, nm tt which transforms in the continuum case to tt {se (z,t), s <t>(z',t')} = -5 S ,. nm tt (2.16) Chapter 2. A theoretical Therefore, s6(z,t) model for ferromagnetic wires with pinning centers 37 becomes the conjugate momentum to the angular position <j)(z,t) so d6 that, upon using equation (2.13) and discarding again the small term J s - ^ — we have oz 2 s yg£-=2K..r[z,t). 3 (2.17) W e note that the Sine-Gordon equation is Lorentz invariant w i t h respect to the characteristic velocity c. A s shown above, this characteristic velocity is proportional to the lattice constant a, which is consistent w i t h a reduction of intermediate lattice sites between any two points thereby increasing the response time of the spin chain. For the time-independent case, equation (2.15) reduces to the fixed B l o c h W a l l solution d r\ 2 m sin r?, (2.18) = 4K se (z,t)- (2.19) 2 and 9 equals zero by virtue of S s z confirming that the magnetization is t r u l y confined to the x-y plane. W i t h the soliton moving, 6 (z, t) no longer equals zero and thus provides a magnetization i n the z-direction s to sustain this motion (as explained i n chapter 1). One might wonder why we bothered deriving the classical l i m i t of the Heisenberg spin chain (2.5). T h e reason for doing so is to provide us w i t h the platform for the quantization of the magnetization field by a method which treats quantum fluctuations about a classical solution. One may seriously objects to our starting w i t h a perfectly quantum H a m i l t o n i a n (equation (2.5)), deriving the classical l i m i t , and then t r y i n g to quantize the system a l l over again! "Should you not have stuck w i t h the quantum mechanical H a m i l t o n i a n that you had and worked w i t h it directly t a k i n g good care w i t h the transition towards a quantum field?," such an astute reader may retort. Very well, but the Recall that the transformation from the Poisson brackets of dynamical quantities A(t), B(t) to the quantum commutators is given by the prescription {A,B}PB -> —(i/h)[A, B]. 3 Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 38 Heisenberg equations of motion derived from this H a m i l t o n i a n are intractable and so must be supplemented by some classical intuition where we suppose that the quantum phenomenon we are interested in manifests itself in the form of small quantum fluctuations (which may have nonetheless important observable effects such as the a-decay) about a stable configuration in classical phase space (such as a classical particle to sit quietly at the b o t t o m of a harmonic potential well). After, we hope, dispelling some misgivings, it seems appropriate to introduce this quantization scheme about a classical solution. T h e exposition of the subject follows closely that given in a book by Rajamaran [36], and we refer the reader to this book for more details. 2.2 Quantization about classical solutions In order to give a qualitative flavor of the method and help to build up our intuition, we consider a non-relativistic unit-mass particle under the influence of a potential V(x) in one dimension as shown in figure (2.1). V I a F i g u r e 2.1: A general potential V(x). I : b 1 c x The point a can be thought as the absolute minimum of the potential whereas point b constitutes a relative minimum. The point b, although an extremum of the potential function, is obviously an unstable one. Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 39 In the classical regime, the particle obeys Newton's law of motion given by ^ - i r = — -7—• dt dx 2 Q u a n t u m mechanically, the particle is assigned a state vector \tp) i n a one-particle Hilbert space whose projection in the rr-basis is the familiar wavefunction ip(x). T h e energy eigenstates are given by solving the time-independent Schrodinger equation: Hipn = \(P 2 + V[x))^ = £ n i n (2.20) represented in the rc-basis w i t h p — —ift-^-. W e w i l l now draw some relationships between dx the classical and quantum cases. Classically, the stable time-independent solutions for the particle w i t h the potential V(x) shown are x = a and x = c, the first case being an absolute m i n i m u m while the second a relative one (the point x = b is obviously an unstable^point). T h e absolute m i n i m u m corresponds to the the classical ground-state of the particle w i t h energy Eg = V(a). 1 In quantum mechanics, such a tranquil state is forbidden by the Heisenberg uncertainty principle; small fluctuations exist about this classical ground state. For energy states close enough to Eg , the potential V(x) 1 may not differ much from V(a) so that a truncated Taylor expansion about x = a is legitimate V(x) « V(a) + ^uj (x-a) 2 + ^X (x-a) 2 3 3 + ^ A (x - a)\ 4 (2.21) where to corresponds to the curvature of V(x) at x — a and the linear term is absent as 2 V(a) is an extremum. Now, for those state vectors \ip) that satisfy the condition A M (X - a) \t/}) < r r u (<ip\ ( X - a) 2 2 (2.22) for r = 3, 4, and higher, the anharmonic terms i n (2.15) w i l l be negligible so as to leave us w i t h a quadratic potential for a l l intents and purposes. Furthermore, we exactly know how to treat a harmonic oscillator quantum mechanically; hence, for those energy levels Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 40 close to Eg and states obeying inequality (2.22), we have the approximate energies E n (2.23) = Note that the first term corresponds to the classical energy, the second, the first order quantum correction with oo being the classical stability frequency, and the last, higher order corrections which may be treated by standard perturbation techniques. Already, i n this simple equation, the right- hand side of which contains classically obtainable values, we can appreciate the many connections between the quantum and classical regimes (recall that even \ ' ° are classical values, they represent successive derivatives of the potential V{x) at the classical position x = a). E q u a t i o n (2.23) is valid only up to some quantum number n at which point it is unwise any longer to consider E n as a low-lying state. Note that inequality (2.22) implies state vectors whose wavefunctions are localized around x = a, once again the classical solution. Now, i f we return to our classical particle, we can easily determine by inspection of V(x) another static solution: the particle sits at the b o t t o m of the well located at x = c, the local m i n i m u m . T h e energy of the particle is then Ef x that at the absolute m i n i m u m . T h e energy Ef corresponds to an excited state above the x classical ground state energy Ef = V(c) which is higher than about which one can also construct a set of quantum energy levels using an expansion of V(x) around x = c instead. Those energy levels will be higher than any of those built around Eg i n (2.23). R u n n i n g quickly through the 1 previous procedure, we obtain : V(x) « V(c) + \ (u'f (x - cf + i A' (x - cf 3 + i A4 (re - c ) 4 (2.24) for the state vectors. \ip) such that A; <^| ( x - c) |^> « r (u/) w 2 ( x - cf w (2.25) Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 41 holds for a l l r > 3. Here, co' is the classical stability frequency at x = c and X' the r corresponding derivatives of V(x) at the same point. Once again, those state vectors have wavefunctions strongly localized about the local m i n i m u m x = c. T h e energy levels built around this excited classical solution is then E, ex n = V(c) + {n + \)Ku)' + 0(X' ). (2.26) r T h e same limitations apply for this expression as for the previous one. Now, here lies the full force of our analogy: for a classical field <f>(x,t) there also exist static solutions, say <f>i(x,t) and (f>2(x,t), for which the energy of (j>2(x,t) may, for instance, be higher than that of <j>\(x, t). T h i s latter solution may therefore be thought of as the classical ground state of the system in the same way as x = a was for the discrete case. Thus, one can also build approximate quantum energy levels about those respective classical field solutions by following the same method outlined above but generalized to fields; namely, by " T a y l o r expanding" the potential functional V[(j>) and keeping only the quadratic part. T h e important feature to bear in m i n d is that the energetically higher classical field solution faix,!) may very well be a non-perturbative solution i n the sense that as A —»• 0 (A being a parameter of the H a m i l t o n i a n which is present i n <fo(x,t)) this solution fails to make sense . . . it simply blows up! In other words, the solution (j)2(x,t) cannot be continuously deformed to the classical ground state solution <j>\(x,t) and corresponds, in a more sophisticated jargon, to a completely different topological sector of the H a m i l t o n i a n spectrum. Before giving the details concerning the transition to fields, the perspicacious reader may wonder what happens i f co equals zero as inequality (2.22) can never be fulfilled. In this case, the curvature is zero, so that the wavefunction is less localized that it would be in the presence of a harmonic-like potential well, and our procedure completely breaks down insofar as a series of harmonic-like low-lying energy states cannot be found (for those vanishing curvatures cases only). One might want to Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 42 construct states keeping only the more significant stable anharmonic term (for instance, A4 x ); however, we assume that most m i n i m a we are interested in w i l l be of a harmonic 4 type. A more interesting case arises when V{x) is independent of x in some region; classically, the particle can be at rest at any of those points w i t h i n this region for which all derivatives co, {X } of the potential equal zero. In other words, there is absolutely no constraint to r localize the wavefunction, and consequently, it w i l l tend to spread w i t h i n this region and adopt some form of wavefunction resembling e^ x (the momentum eigenfunction). T h e extreme l i m i t would be a potential independent of the coordinate X\ over its entire range even though it may be constrained along the other coordinates X{ (we consider a m u l t i dimensional potential V^(x)). T h i s independence of V(x) for one dimension corresponds to a translational symmetry of the system along this direction. Therefore, our procedure needs be modified for those coordinates absent i n the potential V(x); we w i l l deal w i t h this particular problem as the case arises later in our treatment of soliton dynamics. [Indeed, even for fields, there exist symmetries for which the corresponding frequencies co vanish (these are called zero-frequency or translational modes).] In the meantime, in order to give a brief idea of how we w i l l accomplish this, note that for a quantum particle i n a constant potential V{x), where p n are the conserved its energy levels are given by E = V + | (Pn) 2 n momenta; the conservation of which arises precisely from the translational symmetry of the problem. Hence, for zero-frequency modes, one might expect to deal w i t h conserved momenta conjugate to those coordinates displaying the various symmetries (in the Lagrangian formalism, it is useful to think about these as the cyclic coordinates). We are now ready to move on to the case of fields. Consider for simplicity a scalar field (j)(x,t), vector fields being easy to generalize to. T h e dynamics of the field is governed Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 43 by a Lagrangian functional m i = J dx. = (d<j>(x,ty 2 V - ( V 0 ) - u{4>) 2 dt T[4>] - V[J>}, (2.27) (2.28) w i t h kinetic energy 'd<j>{x,t) dt and potential energy V[<b] = fdx {^(V0) 2 --£7(0} (2.29) T h e classical solutions for this Lagrangian is obviously given by the Euler-Lagrange equations of motion d <b{x,t) 8V[(j>] 2 (2.30) df 2 where <50(x,i) is understood as a functional derivative. If we enquire about the static solutions, these are simply given by the generalization of the extremum condition for the functional potential V[<b\. H(x,t) =0 ( 2 ' 3 1 ) F r o m equation (2.31) (and the second order version), one may find the field m i n i m a <fo(x) which can be arranged i n ascending order of energy so that c6 (x) is considered as the 0 absolute m i n i m u m or ground-state of the system (also referred to as the v i r t u a l vacuum). Consider one such m i n i m u m 0 ( x ) . n One may then expand the functional V[<j>] about such a m i n i m u m i n a functional Taylor expansion dU 2 V[4>] = VfoJ + j d x \ L x ) -V + 2 d<b 2 (2.32) where 77(x) = 0(x) - 0 ( x ) , the dot representing higher terms i n ?y(x). T h e second term n is the second order or harmonic variation i n 77(x) i n which the term between brackets Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 44 corresponds to 5 V[tff] 2 Scf) 2 0(x) = 0„(x) d V(x) (note that integration by parts has been used) i n an obvious generalization of ~ dx 2 2 for the one-dimensional case discussed above. T h e linear term is absent by the very fact that we expand about a m i n i m u m so that 8V[<f>] = 5(f> 0. W h a t we want to do next is to diagonalize the operator d U\ 2 -V 1 + d<j) 2 l^(x)=(/>„(x) so as to express the harmonic term i n equation (2.32) i n terms of normal modes /i(x); these being found by solving the eigenvalue equations dU 2 -V 2 + d < /<(x) = uo /,(x) (2.33) 2 ^ tf>(x)=^(x)j These modes constitute an orthonormal basis in terms of which one may write the variation n(x) about the m i n i m u m </>(x) as n r?(x) = J2 C i / i W , i (2.34) so that the harmonic term i n (2.32) becomes simply a sum over the amplitudes of the modes times their corresponding energies : (1/2) ^ (ci) uo . 2 2 i Setting ??(x,i) = 0(x,t) - </>(x) n = X>(i)/ (x) i t and assuming that the fluctuation 7?(x, t) remains small enough over time to neglect the higher order terms i n (2.32) represented by the dots, one may write the Lagrangian L as L = I Et^W] 2 - (v[<Pn] + \ [ci{t)} u ) 2 2 + 0ft (x,<)), s (2.35) Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 45 where the correction (9(?y (x, t)) can always be treated by the usual perturbation tech3 niques. T h e most significant terms i n (2.35) correspond then to harmonic oscillations which (as the one-dimensional case discussed above) we perfectly know how to quantize. Therefore, one constructs a set of low-lying energy states around the classical static field solutions (/> (x) in the same manner as before, n E. n = V[(j) ] + n hJ2( i n + 1 / ) 2 u ) i higher order terms, + (2.36) owing to the quantization of the normal modes coefficients Cj. Once again, we remind the reader that the higher order terms can be treated perturbatively and that 0 (x) n may turn out to be a non-perturbative result i n the sense outlined above. T h e classical expression (2.35) is valid as long as the fluctuation ^(x, t) around a given m i n i m u m </> (x) m is negligible. B u t , it may happen, as before, that a symmetry exists for the Lagrangian L so that V[(j)\ has a m i n i m u m </>(x) which is translational invariant, i.e. one for which n </> (x+a), for any vector a, are also m i n i m a w i t h the same energy y[0 (x)]. In such a case, n n we w i l l obtain eigenfunctions /i(x) w i t h eigenvalues tOi = 0 for which the procedure given above also breaks down. These translational modes are the exact equivalent of those 4 symmetries encountered before, for which V(x) was independent of some coordinates if one is w i l l i n g to view <^ ( ) as a point i n the field space 4>{~x) and consider 0 ( x + a) x n n as a ray or line in this space. A s promised, these anomalies w i l l be treated i n due course. Note from equation (2.36) that the energy levels are built around a classical static solution <f> (x) w i t h the leading quantum correction given by the harmonic oscillators. F o r n the absolute m i n i m u m of V[4>], 0 (x) which is space-independent on account of the o term ( V 0 ) i n V[4>] and which may, incidentally, be degenerate, this corresponds to the 2 2 vacuum state and its familiar associated quanta of the field. The word translational may convey the wrong impression that symmetries may exist only for "rectilinear dimensions", which is not true. One may also have angular-like symmetries. Perhaps, one should adopt the more appropriate term zero-frequency mode. 4 Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 46 In our case, however, we w i l l build quantum energy levels around a classical static solution of the Euler-Lagrange equations of motion which corresponds to a soliton w i t h higher energy than the ground-state solution </>(x) = 0. 0 After this introduction to the quantization procedure about classical solutions, we discuss in the next section the toy-model we w i l l be considering in this thesis. 2.3 T o y - m o d e l for pinning centers in magnetic wires T h e mathematical model we w i l l consider is the one-dimensional Sine-Gordon H a m i l t o dip nian density plus a perturbation term of the form OL — V(x) where a is a small coupling constant: U = H + H 0 with p U = \n n = o ^ V ( 4 0 p 2 + I (2.37) + u, (l-cos^), (2-38) 2 0 (2.39) c and u) being the characteristic velocity and frequency of the system. T h e coupling 0 0 w i t h the gradient of the field itself has been chosen so as to eliminate possible ambiguity w i t h regards to the different values of the S G soliton ib (x) at x = ± o o (ip (—oo) = s 3 0 and ^ ( + o o ) = 2TT for N = 1). T h e field ib(x, t) is dimensionless since it corresponds to angular values and the above H a m i l t o n i a n density can be made totally so by introducing the following scaling t = co t, x 0 = (—)x, c 0 and a j. A oc „ „ T h e reason why we choose this H a m i l t o n i a n density is clear from the discussion at the beginning of this chapter: we found t h a t the spin chain modeling a ferromagnetic wire admits in the continuum limit a Sine-Gordon soliton solution as exemplified by equation Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 47 (2.15). T h e correspondence w i t h equation (2.38) is achieved by setting m in equation (2.15) to u> /c . F r o m now on, we w i l l deal exclusively w i t h the dimensionless H a m i l t o n i a n 0 0 density and drop the asterisks on the space-time coordinates and the coupling constant for convenience. As mentioned in the introduction, we want to investigate the effect of pinning centers on the motion of solitons. Thus, how do we model those pinning potentials? We know by experience that those pinning potentials are very strong w i t h i n a very narrow range. Qualitatively, such potentials can be approximated by a Dirac-delta functions 5(x). In any case, the real pinning potential V(x) can always be decomposed as a linear superposition of such delta functions, V(x) so that, since we w i l l consider a linear approximation about the soliton solution, the response of the system to such a potential can be inferred from that of a single delta function. One may have some difficulties to make sense of the coupling of 5(x) to the gradient — which results i n a a force-like term of the form 5'(x) in the Euler-Lagrange equations of motion. After some reasoning, the reader may convince himself that such an idealized force 5'(x) (which carries only the qualitative aspect of what the real force might be) w i l l indeed constrain a moving soliton between two points once we recall our mechanical analogy w i t h weights and springs discussed in the introduction for which the wave propagating is of an angular sort and not of the usual linear one, which we have come to think about whenever the term wave is evoked. Indeed, such a reasoning will be confirmed later by the derivation of a hyperbolic secant potential given in terms of the soliton's center coordinate. Chapter 2. A theoretical model for ferromagnetic wires with pinning centers 48 A s a preliminary to the study of solitons bound to impurities, we w i l l examine i n the next chapter the response of a soliton to a potential V(x) = Q(x — x ) in equation (2.39). 0 T h i s w i l l introduce us to the notion of emission of radiation consequent to a soliton h i t t i n g such a step potential and acquaint us to the way w i t h which we w i l l similarly c tackle the pinning potential case. T h i s emission of radiation is an example of excitation of the quantum field built around the soliton solution. * Chapter 3 Soliton bremsstrahlung and the m e t h o d of collective coordinates T h e treatment in the first section of this chapter is largely based on a paper by Eilenberger published in 1977 [37]. T h i s w i l l be divided in two sections; the first one treats the subject at a classical level where we seek small-amplitude solutions about the soliton solution, the second and t h i r d go beyond Eilenberger and attempts the j u m p to the quantized version. 3.1 Classical treatment of the soliton bremsstrahlung T h e H a m i l t o n i a n density is given by equation (2.43) and (2.44) w i t h V(x) = — Q(x—x ). 0 A n integration by parts for allows us to rewrite it as with K = l i m _ . a ; ) + 0 0 %j)(x,t) which we consider as being equal to 2n, the l i m i t i n g angu- lar value for a soliton w i t h topological index N equal to one (cf. introduction). A s a constant appearing i n a H a m i l t o n i a n is benign insofar as the dynamics of the system is concerned, we may disregard it totally. Hence, the perturbation H a m i l t o n i a n density can be considered as being equal to +a ip(x,t) 5(x — x ) where —a 5(x — x ) — F(x) 0 0 plays the role of a Dirac-delta force acting on the system. Since the Euler-Lagrange equations of motion for H 0 are linearized about the soliton solution (cf. introduction), the linear response to a general force F(x) can be determined by superposition. 49 Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 50 We consider a fast soliton moving w i t h speed v to the right w i t h the step potential positioned at x , the location of the impurity or defect. Intuitively, the soliton begins to 0 feel the step potential only when its center moves in its vicinity. A t this point, a force acts on the soliton which has it slow down until it has gone beyond the potential's immediate range. Thus, we expect the soliton to emerge from its interaction w i t h the potential w i t h a speed v' which is less than the i n i t i a l speed v. However, since we suppose a very fast soliton, v' should not be so different from v in the presence of a soft potential owing to the smallness of the coupling constant a. B y conservation of energy, the impurity should in principle also gain a small speed; however, we assume that the restoring forces in the material are strong enough to counter this motion so as to leave the impurity stationary. These forces in the material are not accounted for i n the above H a m i l t o n i a n as we isolate the relevant part of the universe and "ignore" the surrounding whose effect, in any case, manifests itself i n the renormalized masses of the various parts of the system under consideration. T h e important point to note, here, is that radiation can also be emitted as a result of the collision i n the same way as we expect heat to be given off when a fast ball grazes the surface of a large massive block. A still better analogy is the case of a charged particle which, being decelerated by some force, gives off radiation (for instance, an electron colliding w i t h a nucleus surrounded by its electrical field). One may then substitute the charged particle for our soliton and the decelerating force for that provided by the step potential; this is the origin of the term "bremsstrahlung" to describe this process of a soliton collision. Solving for the equation of motion derived from our H a m i l t o n i a n density % , we obtain ip - ib" + smib = -a 8(x - x ). (3.1) 0 Linearizing about the small-amplitude (b(x, t) as outlined i n the introduction leads to <j>(x,t) - (f>"{x t) + (1 - 2sech [ (:r - vt)]) <j>{x,t) = -a 2 t 7 5(x - x ). 0 (3.2) Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 51 For t —¥ — oo, the i n i t i a l condition where the soliton is far away to the left of the step potential, the hyperbolic secant term in (3.2) v i r t u a l l y vanishes to leave us w i t h 4>(x,t) — cf)"(x,t) + (j)(x,t) = —a5(x — x ). (3.3) 0 We suppose that the behavior of <f> for this i n i t i a l condition is relatively independent of time so that a Fourier transform of equation (3.3) can be performed w i t h respect to the position only, (k 2 + l)(j>(k) = -a, where (f)(k) is the Fourier transform of <f)(x) in the vicinity of "t ~ —oo." One can easily solve this equation, transform back to the rc-space, and obtain *o(x) = -|e-l*—I (3.4) as the i n i t i a l small-amplitude field generated by the sole presence of the step potential. T h i s can be referred to as the "dressing" of the i m p u r i t y located at x . 0 We expect, then, to recover the same dressing at t —> —oo w i t h the soliton moving w i t h a slightly different speed to the right of the impurity. It is convenient to transform to a different inertial frame i n which the soliton is stationary. In this reference frame, we have instead the i m p u r i t y or step potential coming from the right w i t h speed v (cf. fig. (3.1)). V(T -T) C Figure 3.1: T h e i m p u r i t y wavefunction "moving" towards a stationary soliton. In the rest frame of the soliton, whose center is located at coordinate z = 0, it is the impurity center, located at — ( T - r) where r is the soliton's proper time, that moves towards it with speed v. The impurity dressing assumes an even decaying exponential form at large distances from the soliton's center. Z V 0 Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 52 After performing the required Lorentz transformation [38] x — 7 ( z + vr), t — 7 ( r + vz), v^r =x , 0 (3.5) 0 the last equality m a k i n g the i m p u r i t y collide w i t h the stationary soliton located at z = 0 at time r ~ r , we obtain 0 $(z,t) - 5(z + V(T - r ) ) , (t>"{z,t) + (1 - 2 s e c h [ ( z - vt)]) <t>{z,t) = 2 7 0 (3.6) 171 the partial derivatives being w i t h respect to the new coordinates (z, r ) . Since we consider a soliton w i t h topological index N equal to one, i.e., a kink, 7 > 0 and we may drop the absolute sign i n equation (3.6). E x p a n d i n g 4>(z,r) i n the orthonormal basis {fb(z), f (z)} that we found earlier, k 4>{*,T) = °-^-h{z) + Jdk 1 c (r) f (z) k (3.7) k where (c^(r))* = c_fe(r) ensuring thereby that <j>(z,r) remains a real field, we get \ {[c (r) + c ( r ) (1 - 2sech z)] f (z) - 2 b 6 b c b + [dk {[c (r) + c (r) (1 - 2sech z)] f (z) 2 fc J k k f»(z)} - c f' \z)} k k = - - 5(z + v(r - r ) ) . (3.8) 7 0 In order to obtain the equations of motion for the normal modes' coefficients, we simply m u l t i p l y the above equation by the complex conjugate of a n o r m a l mode f {z), integrate p the resulting equation w i t h respect to z, and use the orthonormality relations. D o i n g so w i t h the zero-frequency mode, /(,(z), knowing the following relations +00 d z [ ( l - 2 s e c h z ) / ( z ) - n'(z)f (z)} = 0 dz[2sech zf (z)f (z) = 2 2 b / -00 / 2 j—00 k b + f' \z)f {z)] k b 0, because of the scaling we performed on our Hamiltonian, the normal modes (1.52), (1.54) become f ( x ) = 2 sech(z), f ( x ) = -j=^-e (tgh(z) - ik) with eigenvalues wf = k? + 1, k = (LJ /C )1, and lkz b (h,h)=s, k (f ,f ) k P = 0 6(k-p)™" 0 Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 53 leads to 2a Cft(r) = sech(-y(r - r ) ) 0 7 = ~ / 6 ( » ( r - r 0 ) ) . (3.9) One can infer a good deal of physics from equation (3.9). T h e mode fb{z) provides the translational motion of the soliton for small displacements since then * «*) + ^ ^ f (3.10) Thus, equation (3.9) can be viewed (at least for small displacements) as some form of Newton's equation of motion where the n o r m a l i z e d soliton's center coordinates c ( r ) / 8 2 D responds to a force i n the immediate vicinity of the i m p u r i t y (realized when r —> r ) . 0 T h i s is even more so by the fact that v is large since we deal w i t h very fast solitons (in the lab frame). Furthermore, as the coupling constant a makes the force small, c (r) will only acquire a small speed by the time the force becomes v i r t u a l l y unimportant. This b smallness of c (r) b guarantees the approximation of a moving soliton ipsil (z + Q,(r)/8)) by "0*( + h ( ) / 8 ) w i t h 7 ^ 1 . Hence, this confirms our earlier expectation of a soliton z c r emerging from a collision w i t h the i m p u r i t y w i t h a speed slightly less than initially (in the moving reference frame, this translates into a soliton acquiring a small velocity i n the negative direction). Note, however, that this interpretation breaks down once c (r)/8 b becomes important since then the approximate equation (3.10) is no longer valid; this w i l l occur, though, only after the soliton leaves the immediate vicinity of the i m p u r i t y at which point the force is vanishingly small and the soliton cruises w i t h a nearly constant slow speed (in the moving reference frame). I t is normalized with respect to the the mass of the soliton which has the value of eight in dimensionless units. 2 Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 54 Now, tackling the other modes f (z), we obtain w i t h the help of the following relations p fk{z) = - f (z)[k +2 sech z], 2 2 k f£i*) +00 / = dz f (z)f;(z) f (z)(tgh (z)-sech (z)), 2 2 b (sech (z) + tgh (z)) = 0, 2 b 2 -00 the expression c (r) + OJ C (T) = 2 p P -SLf;( ( -T )). V T (3.11) 0 T h i s equation is rather difficult to solve unless we resort t o a trick proposed by Eilenberger. We decomposed the small field amplitudes Cp(r) as C {T) = c {r) + c (r) P p (3.12) op with W h a t does this decomposition correspond to ? One may find out what the terms c ( r ) op generate i n the (z, r) space: +00 / = dk c (r) f (z) ok k -oo H S T / - . * ^ ( * ^ ) + « (tgh , tgh r + - r.)]) H - t g h 0 t g h [ w ( r (-a) [ 2 \ v\ z + v ( r T n ) - r )]} 0 [V(T-T )}) \ {j+tgh z) (7+tgh 0 (72 + 1) (l+tghz)(l+tgh[t;(T-r )])l 0 (7 2 7 Note that for | v(r - T ) | -> +oo, (j> {z,r) 0 0 - (a/2) e ^ I - z+ u ( + l) r _ r J- ° ) I which is exactly the Lorentz-transformed expression of the impurity dressing (3.4) and corresponds Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 55 to the expected situation of the i m p u r i t y being v i r t u a l l y unaffected a long time after its collision w i t h the soliton. So, m a k i n g use of decomposition (3.12) and substituting equation (3.13) in (3.11), we subtract the i m p u r i t y dressing contribution and obtain + LO2C(T) = + uo c {r) £(r) c [r) v p p ~ W i k ( A .2 a e + 2jj*^ ikv(r-r ) i 7 ) 3 0 l e [ + cosh «(r-r ) 2 + ± 1 dr cosh v(r - r ) v ( V { 0 2 } ] ) l 0 T h i s equation is easier to handle i n frequency-space; so taking the Fourier transform on both sides and using the integral dx e sech (arc) ikx 2 = k cosech(--) we get ~ iane ° 7(0; — kv) lu,T CP(W) = ~(LO + IC) -LO 2 oo ^ 2 2 p p (<*V^) + + k) 2 s\nh[{Tt/2)((u /v) k k , + k)\ ' , 1 } where we added an infinitesimal imaginary part to the frequency denominator to place the poles i n the lower half of the complex plane. Transforming back i n the r-space and setting to as a complex variable u, we have c( )= ~ l i m / du e H « ( 27T(7 + A: ) e ^ o i Z T P { } m 2 r 7 2 u " r w {(u + it) e 2 r H where T R k (u-kv)((u/v) + k) - LO ) sinh[(7r/2)((w/u) + k)} 2 is the path of integration i n the complex plane and the subscript R reminds us that we w i l l let the radius of a semicircle tend to infinity i n order to evaluate 5p(r). For r < T , we need to close the contour of integration i n the upper 0 half-plane (cf. fig. (3.2)); since there is no pole, c (r) = 0, which is consistent w i t h what we expect p before the soliton appreciably begins to feel the impurity at r w T . Therefore, our choice 0 for p u t t i n g the poles below the real line turns out to satisfy the i n i t i a l conditions of the problem. For T > T , we close the contour i n the lower half-plane, make use of the Residue 0 theorem and the fact that jv (LO /V - p ) 2 2 (7 + k ). 2 2 2 1 v7 Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 56 (a) F i g u r e 3.2: T h e integration path. The poles, designated by x on the graphs, are barely located in the lower half-plane because of an infinitesimal imaginary part added to the frequencies — ui and ui . p p The integration path is made up of the line from —R to R on the real axis and a semicircle CR. of radius R in the upper half-plane in (a) and in the lower i n (b) to accommodate the initial conditions of the problem. This complex integration is necessary to evaluate the inverse Fourier transformation of the Fourier coefficients in equation (3.14). For that purpose, we let R tend to infinity. C (OJ) P to find CT = P -iu} (r 2u jv iu) {r - T ) - T) p OJ7T 0 0 ^sinh[(7r/2)(w/fc/i; + k)} 2 0 v smh[(n/2)(-u /v k + k)} (3.15) So, subtracting off the intrinsic i m p u r i t y dressing and its effect during the collision, equation (3.15) gives us wave-packets r+oo f+oo dz C T (J){Z,T) = K f (z) k J—oo which move w i t h a dispersion relation u> = y/1 + k . T o second order i n a, we can calcu2 k late the energy of the system i n the small-amplitude approximation using the following decomposition of % 0 ri 0 — = 'Hose — % s + 'Hose + 2 ^ * + ^) with + + (1-cos + ^ cos(&), 2 Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 57 where the field is written as = *b (z) + </>(Z,T), ip{z,r) s ib (z) being the static soliton located at z = 0 and (J)(Z,T) given by equation (3.7). s Integrating with respect to the coordinate z leads to +oo / Hose dz H s = M = 0 8, the rest mass of the soliton, -oo ^ dz Hose +00 / -oo -oo r+oo +oo / -oo = 2 {^f 1 dz Hi = 0, dzU = aib (x ) p since we assume | (j) \ + J dk I cl(r) | +u,l | c ( r ) | } (3.16) 2 2 2 and s 0 + acj)(x ,r) « 0 ib (x ) s 0 1. After a sufficiently long time, the two wave packets, moving i n opposite directions in the soliton rest frame, w i l l become entirely detached from each other and, presumably, from the impurity dressing as well on account of the impurity's great speed. Therefore, we can ignore the mixed terms i n | Cfc(r) | and isolate the terms pertaining to the wavepackets 2 alone an ° ^ k T ) ~ - k{r i(7UJ e - r) 0 2UIJV sinh[(7r/2)(( /t;) + A;)] ( ? ) W/k with <J = ± 1 ; the sign referring respectively to a right- or a left-moving wavepacket and the superscript r standing for "radiation" which is how we coin this particular excitation of the small-amplitude field. Plugging equation (3.17) into (3.16), we obtain for the energy of this radiation Chapter 3. Soliton bremsstrahlung Ur{r > r ) 0 = and the method of collective Y, d <r=l J "1 k I c*,* T '2jv' Hence, c r ka 58 (3-18) -°° r+oo = coordinates 1 dk J-oc ^ UT -) smh [(Tv/2)({aoj /v) 2 k dkcosech [(7r/2)((oj /v) 2 1 k + k)] + k)] can be interpreted as the amplitude of the excited normal mode k w i t h equation (3.18) being the usual sum of individual mode energies. A general step potential V(x) can always be constructed using the decomposition +oo dx V(x ) 0 / 0 Q(x - x ), 0 -oo from which we obtain the corresponding force +O0 dx 0 / F{x ) 0 5(x — x ) 0 -oo with F(x ) 0 = -V(x ). 0 It is then a simple matter to calculate the radiation amplitudes by linear superposition. Naively, one would assume that equation (3.17) is the correct result at the classical level; one for which the amplitudes c\ simply carry over to the quantum world as the amplitudes of the corresponding quantum oscillators. However, this classical result is incorrect, and the reason has to do w i t h the warning we issued when we talked about the quantization scheme about a classical solution: the presence of a translational mode, namely in our case, fb(x) w i t h its eigenvalue ui = 0. In classical or quantum mechanics, b the decomposition of the harmonic H a m i l t o n i a n density (3.16) is insufficient since the translational mode has no confinement simply means that, in the f (x) b to approximate it as a harmonic mode. This "direction", the anharmonic terms in a functional Taylor Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 59 expansion of a potential V (cf. equ. (2.32)) are likely to be as important as the harmonic term, whether the state is classical or quantal. A s promised, we give i n the next section a method to handle those zero-frequency modes. W i t h this method, we w i l l see that equation (3.16) is nonetheless valid only after t a k i n g into account some facts specific to our physical situation of a fast soliton's colliding w i t h a soft impurity. 3.2 T h e m e t h o d of collective coordinates A s pointed out i n chapter 2, eigenfunctions of the operator w i t h eigenvalues ui equal to zero may exist for an expansion about some m i n i m u m c6 (x). 0 n T h i s arises from continuous symmetries displayed by a Lagrangian and i n particular from the translational invariance of m i n i m a 0„(x) for which V[0 (x + a)] = 0 (x), V a . n n A m i n i m u m 0j(x) which is independent of x w i l l not produce a zero-frequency mode for the reason that <j>j{x. + a) = 0,-(x) and, therefore, can be considered as a point i n "field space". T h e zero-frequency eigenfunctions / (x) correspond to the infinitesimal 0 translational generators of the m i n i m a </>(x) i n the sense that n 5(f> (x) = n = 0 (x + <fa) - 0 (x) n n (*a)/ (x) 0 In our case, this is the function fb(x) = 2 sechrc which is exactly the derivative dib (x)/dx s of the soliton solution. In chapter 2, we gave an idea of how symmetries i n a classical Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 60 Lagrangian translate over to the quantum world: their respective conjugate momenta operators turned out to be conserved. Before stepping into relativistic field space, let us deal w i t h a simple example in non-relativistic one-particle quantum mechanics. T h i s w i l l provide the reader w i t h the basic physics behind the collective coordinates method applied to fields. Consider a particle in two dimensions w i t h the following Lagrangian L = \ {x\ + xl) - V(y/x\ + x%). (3.19) The potential enjoys an obvious continuous symmetry: it is rotationally invariant. Let x ^ 0 be a m i n i m u m of V. This is not the only classical solution since V ( x ) such that r p | x P |= | x r | w i l l also be a m i n i m u m . If we naively apply our method of quantization about the classical solution x , we obtain r 2 #n n 1; 2 = V(xr) + £ fa + 1 /2) ftUJ { + higher orders (3.20) 1=1 where are the eigenvalues of / dV 2 [dxidxjj \ x = x One finds that one eigenvalue, say ojj, equals zero on account of the rotational symmetry of the potential V. T h i s zero eigenvalue implies that no confinement exists for the wavefunction i n the angular direction. Therefore, it fails to be localized on the circle of radius | x r |= R although, radially, it is so owing to the nearly quadratic potential for each of those points x p i n the radial direction. This shows that the quantum stationary states of the system are not localized around any specific classical solution, w i t h the consequence that the approximate equation for their energies (3.20) is incorrect. One simple way to deal w i t h this problem is to change from cartesian to polar coordinates, which change is naturally dictated by the symmetry of the potential. The Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 61 Lagrangian (3.19) transforms then to 1 - r 2 + - r9 2 2 - V(r) (3.21) One immediately notices the absence of the angular coordinate which makes it ipso facto a cyclic one, so that its conjugate momentum I = r 9 becomes a conserved quantity. 2 Therefore, one may readily reduce the problem t o a one-dimensional case involving the radial coordinate r and the momentum conjugate to 9, I, considered as a given constant. Thus, we can write the Lagrangian as U = r £ + 2 \f + 2 ™ + V (r;l) (3.22) eff where V / / ( r ; /) becomes an effective potential acting radially on the particle w i t h a given e angular momentum I. T h e term Z /(2 r ) simply corresponds to a fictitious centrifugal 2 2 potential. N o w , we can apply the usual procedure to this Lagrangian by finding first, a static solution r(t) = R[ given by the extrema of the effective potential V ff(r; I), e V f(r-l)\ ef dr r=R t dV dr (3.23) r=Ri and second, by building a series of low-lying quantum energy levels about this solution in the harmonic oscillator approximation, Ej n « V (r;l) = V{Ri) +. \ {n + l/2)hLO 2 eff + - (n +1/2) hco whose eigenvalues Co are given by d V eff 2 LO P dr 2 \r=R t + ^ (3.24) Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 62 Note the difference between equations (3.20) and (3.24) where the indices i n the former, {ni,ri2}, have been replaced in the latter by {n, 1} w i t h / being a number associated w i t h the coordinate responsible for the symmetry of the Lagrangian. T h i s difference is a manifestation of our changing from one set of coordinates, the cartesian, to another, the polar. In the field case, we w i l l encounter the same difference resulting from a change of coordinates. A more i l l u m i n a t i n g way of looking at this reduction insofar as it bears on its generalization to fields is to view the construction of quantum levels in equation (3.24) as being built around a classical time-independent solution instead of a static one. Indeed, the equations of motion derived from (3.21) are r• = P + ^ dr r = 0 and • d(r B ) 2 dt If we set r (t) d — {r i(t),0 i(t)) c c 2 0. = (i?;, ut) w i t h u = l/R 2 a constant, the second equation is automatically solved while the first reduces to equation (3.23). Therefore, the approximate energy equation (3.24) can be thought of as an expansion about a time-dependent rotating solution which samples uniformly a l l points on the circle R and which ought to give a better expansion than that about a specific point x on this circle (which rotational r symmetry does not favour in any way). After this short excursion in the discrete non-relativistic realm, we w i l l deal first, w i t h relativistic but classical fields and then, move on to the quantized format. In order to keep things simple, we w i l l confine ourselves to a scalar field in one d i mension w i t h one non-degenerate static solution. T h e generalization to fields in multiple dimensions is straightforward. Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 63 Consider the Lagrangian +oo 1 dx - (0 (x,t) - (<b'(x,t)) ) 2 / 2 Z -oo - U{cb), (3.25) which is the one-dimensional version of equation (2.27), and a static m i n i m u m which we designate by (j) (x). We recall that this m i n i m u m is a solution to p 5V[d>]\ 0, 5(b namely, it satisfies ^L MM _ d dx 1 o. = d(p (3.26) p However, by translational invariance, (j> (x) is not the only solution; indeed, 4> {x + X) p p for any X is also a m i n i m u m . T h i s invariance signals that we w i l l run into problems i f we attempt a harmonic expansion about a specific solution 4> (x + X), p arbitrarily, just as we did about the specific solution x Lagrangian. r where X is fixed for the rotational symmetric W e exhort the reader to be vigilant because we w i l l use the exact same methods that we used before in order to get out of the fix! T h e naive expansion would be to express <f>(x) about the classical static solution <f> (x+X) in the orthonormal basis {f (x+X), p b f (x+X)} which consists of the eigenmodes n (cf. equ. (2.34)) that diagonalize the functional potential V[4>] w i t h f (x b + X) being a zero-frequency mode (i.e., one whose eigenvalue equals zero), 00 <j>{x,t) = <f> ( + ) x + Cb(t)f {x) x p + £ b c (t) n f (x n + X). (3.27) n=0 T h e approximate energy equation is then given by E = V[4> ] + ft 53 ( i + V ) i n ni n 2 w + higher order terms. (3.28) We have just recapitulated what we already mentioned i n section 2.2. Due to the zerofrequency mode fb(x) we know that this expression is incorrect. reader to bear closely with us. Now, we entreat the Chapter 3. Soliton bremsstrahlung and the method of collective coordinates In function space, the basis {fb(x), f {x)} 64 acts as a reference system i n the same way as n the basis {1,3} does for a 2-dimensional space w i t h decomposition x ( t ) = x\(t)i Hence, c {t) and c (t) act as the coordinates b n +xzit)]. of the system w i t h respect to a reference frame whose origin is at the "point" <f> (x+X) in field space. In order to fix our problem i n p the non-relativistic case, we recall that we made a change of coordinates from Xi(t), to the polar representation r(t),6(t); X2(t) which representation naturally lends itself to on account of the rotational symmetry the Lagrangian enjoyed. Furthermore, we pointed out that the expansion was made about a time-dependent classical solution instead; for which the symmetry-related coordinate 6(t) somehow sampled uniformly a l l the points connected by the rotational symmetry. T h e reader may have guessed already that one would have to change the set of coordinates {c (t), c (t)} to one including the coordinate b n responsible for the symmetry of the Lagrangian. W e know what this coordinate is: this is the value X i n <j> {x + X). n Therefore, our new coordinates would be something like {X(t), q (t)} where q (t) should be associated w i t h the usual basis functions f (x) n n n as the wavefunction is perfectly confined i n these "directions" (as it was i n the radial direction for the previous one-particle case). T h e proper expansion is then as follows: 4>(x,t) = (f> [x-X{t)) + £ n q (t) f (x-X(t)), n n (3.29) n=0 where f (x n — X(t)) is as appropriate as f (x) n due to the translational invariance of the Lagrangian L (equ. (3.25)). Notice two things: the zero-frequency mode fb(x), which caused so much trouble, has T O T A L L Y disappeared; a n d the expansion is made about a time-dependent solution <j) (t)(x — X(t)) which samples all m i n i m a related to each other n by the translational symmetry. This sampling w i l l be uniform as we discover later that the momentum conjugate to the coordinate X(t) is a conserved quantity i n the same way the angular momentum r 9 previously was. Expression (3.29) is what we refer to as the 2 method of collective coordinates [36, 41, 42, 43, 44] where X(t) is called the collective Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 65 coordinate associated w i t h the corresponding translational symmetry. For convenience, let us denote q (t) = X{t). A s mentioned earlier, the potential term b V{4>] is translation invariant and, therefore, does not depend on X(t). Hence, using expansion (3.29), we obtain r+oo V[t] = = ~ (-1 |-(0') + 2 f_™dx U(t)} V(M ,{q }) 0 V[4> ] n ll l + higher order terms w i J ! + n (3.30) in the small-amplitude harmonic approximation w i t h V[0„] and {w } being the same as n in equation (3.28). We note that 2 ± M 0 corresponds to the "rest mass" of the field for this particular classical excitation, (j) {x). n T h e n the kinetic energy term r+oo +oo _ 1\ / dX - (<j>f Z 2 -00 becomes T = 1 f ^{Qb(t)D q (t) bb r+oo + J dk q (t) D b b bk / dk q {t) D •oo ° ' ^ Ut) Dij <ij{t), hj=b,—oo k = 1 2 + k r+oo /-+oo +oo A q (t) kb q {t) + / b J— oo J—OO 0 +oo / 53 +oo Oj k kp q (t) p J (3.31) where i=6,—oo dk dp q (t) D = a 6 + 5Z °« i=:—oo Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 66 for any indexed quantity a;, and Dij = D, (3.32) jti +oo c -oo dx U' {x) I. / r+oo +oo / = k p / dkq (t)f' (x)\ J—oo k ^ rfa: j ^ ( r r ) + / -oo 6{k+p), dpq (t) p r+°° / dp q (t) f' (x)j p p f (x) k dx f' {x) f {x), p (3.34) k ./—oo (3.35) n are orthonormal to each other by k Also, we have (D - )* n (3.33) J ' ^| in which we made use of the fact that <f)' (x) and f (x) virtue of <f>' (x) = fb(x). k /-foo C +oo D, + n -\ 2 b = D . k E q u a t i o n (3.35) is the simple btk generalization to continuous case of the Kronecker function 5 k p in the discrete case. (If the minus sign bothers, recall that expansion (3.29) is made up of complex eigenfunctions whose coefficients are related by (q (t))* = q~ (t) to ensure a real field.) k In equations k appear w i t h their arguments set to x k (3.33) and (3.34), 4> and the eigenfunctions f n instead of x + X(t) since under the integral sign it does not matter. The field Lagrangian written in terms of variables {q (t), q (t)} is b 1 L 2 +00 / E m i,j=b,-oo where V{M ,{q }) 0 k k DiddQk}) 4 2 V(M , qj(t) ~ V(M , 0 {q }) k (3.36) k b b 6 {q }). 0 Note i n the above equation that the coordinate q (t) = X(t) its speed, q (t) = X(t), + 0{<j) ) nowhere appears, though does. In order to derive the H a m i l t o n i a n , we proceed as usual by finding the canonical momenta 7Tj conjugate to q^t) dL ^ = = qiV<) u °°' E E>ijQj(t) j=b,-oo + Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 67 so that formally the inversion reads as +00 / E Qi(t) = (S)r>i(*)> (3-37) j=b,—oo and so obtains +00 / i=6,—oo = | +oo g E / TO T i + ^(Mo, 1 fe}). (3.38) i,j=6,—oo In appendix B , we derive the inverse of the m a t r i x D and, therefore, merely state the results: = a n d = - D b ' - k ^ - for * , p # 6 p where V = determinant of Z> = -^(c* + M ) 0 2 with (3.39) +O0 a / /-+O0 / -oo \J—oo It is very i l l u m i n a t i n g to see to what 7r , the momentum conjugate to X(t), 6 T h e total momentum for a given field if) is defined as J-oo In our case, this t o t a l field momentum is dt dx corresponds. Chapter 3. Soliton bremsstrahlung p = and the method of collective coordinates - T d x ^ ^ - dt dx J-oo +oo / ^ ( / f+oo dx { + \ . dk q (t) f' (x)j k X(t)- k / / dx U' (x) n -oo + / \ k dk q (t) / = dx {<t>' (x) + / n J — oo dk D q {t) dk D q (t) b>k -oo +oo b ) f e k f+OO / k / f' (x)} k X(t) k f+OO -oo +oo dk q {t) k / +O0 / f (x)j . 2 dk q {t) f (x)) J—oo dk q (t) k + £J \ f+OO > r+00 x +O0 P 68 dk q (t) f' (x)} k k f (x) k J—oo k k -oo n- Hence, ir is none other than the total field momentum of the system; which is quite b reassuring as we know that its conjugate coordinate X ( t ) can be thought of as the "centerof-mass" of the field <b(x, t) giving an idea of its location in space. F o r example, for the Sine-Gordon Lagrangian, <f> (x — X(t)) = ib (x — X{t)) where X(t) locates the center of n s the soliton (even though the Lorentz factor 7, which acts as a scaling factor, is absent). Furthermore, qb(t) or X(t) is a cyclic coordinate i n H a m i l t o n i a n (3.38) and, therefore, we expect P to be a conserved quantity in the same manner as the angular momentum I = r^O was for the discrete case. O n e can surely appreciate, by now, the many analogies we can draw between the discrete and field cases. U s i n g 7T(, = P, we can rewrite the classical H = ~ +g f_l dk + field-relativistic \l H a m i l t o n i a n as jdkd [5{k+ )- >- £ ^} Db P P + V(M {q (t)}). 01 n k h 7T k TTp (3.40) Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 69 Now, we must quantize this H a m i l t o n i a n . We w i l l use a method developed by Christ and Lee [41] which makes use of the usual canonical quantization scheme based on the H a m i l t o n i a n . T h i s latter consists i n promoting the field 4>(x,t) and its conjugate field momentum n(x, t) into operators obeying the equal-time commutation rule (h = 1) [fx,t),Tt(y,t')] = i5{x-y)5 ,. (3.41) tt T h e quantum H a m i l t o n i a n is simply derived from the classical one by elevating the fields to the rank of operators H 1 l = - r+oo r+°° / „ „ dx{t {x,t) 2 + V((f>{x,t))} (3.42) Z J—oo with and Formally, the commutator (3.41) is fulfilled by defining 5 ir(x,t) 6<f>(x, t) in the sense of a functional derivative operator on the appropriate Hilbert space. However, what we want to do is to perform a change of variables from 0(x, t) to {X(t), qk(t)} which can be executed unambiguously. In order to see things i n a more i l l u m i n a t i n g way, the reader can always think of <f>{x,t) as (f> (t), an operator w i t h a continuous index x , that x acts on "functions" ip(y) w i t h y = (y )\ i n other words, f u n c t i o n a l . Thus, n (t) is x x correspondingly viewed as Tt (t) which translates as the operator —i5/8(j) (t) x x on this space of "wavefunctionals." Hence, in effect, we want to achieve the change of variables {k(t)} -> {X{t)Ak{t)} (3.43) Chapter 3. Soliton bremsstrahlung and the method of collective coordinates Under this transformation , how do our operators Tv (t) = —i5/5(f) (t) x x 70 look like ? W e stress the fact that the only difficulty i n the quantum H a m i l t o n i a n (3.42) under this change of variables lies i n the momentum operator Ti{x,t) since the potential operator term V((j)(x, t)) can be handled i n the same manner as before, V[<f>(x,t)] = / dxV{<j>(x,t)) = V(M ,{q (t)}) 0 n J—oo w i t h the simple prescription of elevating all the q (t) n i n V(M , {q (t)}) 0 n into operators. 3 This is obviously due to the quantum mechanical decomposition of the field operator <f>{x,i) as +oo = (t> {x-X{t)) + n dk q {t) f (x k / - X(t)) k (3.44) -oo p{x) which lends itself to the same treatment of V[(f)(x, t)} on account of the integration over the x coordinate. G o i n g back to the transformation of % {t) under the change of variables (3.43), we, once x more, t u r n to classical concepts to give us a hint of how to handle this problem in the continuous case; this time, the classical concept is the standard rules of calculus. Were we to have only a finite number of variables ( c v i , a ) instead of the continuum-indexed N ( 0 ) , the L a p l a c i a n j dx5 /84> 2 X x would be given by ^2 N »=i v v c i T h i s L a p l a c i a n would act on wavefunctions of the form ip(osi,a ). N B y the standard rules of calculus, a change of variables {o^} —> {/3i} would have the L a p l a c i a n take the form Pragmatically, this involves covering the coordinates q (t) with a nice little hat n Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 71 where ^ Sa 8a k Pi 5 *=l k Pj d and B = Det(B)ij. It is now easy to see the parallel w i t h the continuum case; the following equivalences can be made N < E > <t>x i 1=1 > J +00 I N A E <—> Qi > dx <—• E i,j=l i,j=b,-oo ' so that the generalization of equations (3.46) a n d (3.45) amounts to Bij = J+°° dx B^ = J+°° dx Sfc 5_$x Sqi 5qj ' S(j)(x) 5<f>(x) oo +oo / -oo dx 7i(x,t) = r+oo / Sqi 5qj dx TT (t) x J—oo g2 +oo I dx oo +00 ^2 dx —00 One can easily convince himself that the m a t r i x B corresponds exactly to the matrix D we encountered before. Hence, w r i t i n g the full H a m i l t o n i a n i n the coordinate basis {qb(t),q (t)}, k we obtain Chapter 3. Soliton bremsstrahlung and the method of collective coordinates \C^^m ^m 72 - {Drl) 1 1 5 , r- (3 47) 5 \ £> {im^^sw) + s w ) ^ ^ j x \ t ) + j V(M ,{q (t)}), 0 k where we chose to unmask q (t) to reveal its true nature, X(t), and affixed to H the b subscript c to remind us that such an operator is represented i n a basis. Setting 5 5 ~ JW) l = U t ) 4 p eim a n d * = k ( t h consistent with the commutation rules [qk(t),it (t')] = iS(k+p)5 , p [X(t),P(t')] u = i5 , tt (3.48) [X(t),n (t')} = 0 = [&(*), P(*0] k (which corresponds exactly to the usual canonical quantization scheme applied on coordinates {X(t), we can finally write down the representation-free quantum q (t),P(t),Tt(t)}), k H a m i l t o n i a n as H = 1 ./IA, 1 1 /•+«> 1 . , D+ U £ > + 2 \ II ^ dkdp k j^+P) ~ ^ j k + V'(M ,{q (t)}), 0 where f+OO V = -D b>b 1 + / dk D+ D k 2 btk (3.50) k Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 73 with f+oo +OO / (D~\ b (D-\ k (D- )^ 1 ^ = -I, = ^ / dx<j>' {x) \ dk q {t) n f' (x)j k k (3.51) = = r+oo r+OO 5 = { k (3.52) + p ) ( 3 . 5 3 ) Note that a l l the various operators A j and © are made so by virtue of {q (t)} of which k they are function (cf. equs. (3.32) - (3.35)). Also, as before, the operator X(t) is completely absent from equation (3.47); only its conjugate momentum P comes along for the ride. Finally, we bring the reader's attention to the fact that the above H a m i l t o n i a n is exact i n the sense that no approximation whatsoever has been made i n its derivation. H a d we carelessly quantized the H a m i l t o n i a n (3.40) by using the commutation relations (3.48), we would have been stuck w i t h the ordering ambiguity of the second and third terms since D j and 7r do not commute. O u r expedient to treating 7r as differential i% fc fc operators has just naturally taken care of this ordering. B y virtue of the absence of X{t) in the H a m i l t o n i a n (3.49), we obtain and thus, P, the total field momentum operator, is conserved. Therefore, we can choose simultaneous eigenstates of H and P. Let us consider the eigenvalue P = 0. In this case, the H a m i l t o n i a n reduces to H (P=0) = -j l + jdkdp-j=t ^8(k+p) k V'(M ,{q (t)}). 0 k - ^X^j 4i* v (3-54) Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 74 T h e above is an exact result ( V is just expressed in terms of M and {qk(t)} 7 through 0 the decomposition (3.45); in principle, this can be an intractable expression except in the harmonic approximation). Now, we can start our expansion and keep only the leading terms since | p(x) | - C | (t> {x — X(t)) |; this translates in keeping low powers of qk(t) and n bearing in m i n d that +oo dx <j) (x - X(t)) n / -oo {qk(t)}- is huge compared w i t h any expressions including only Let us deal first w i t h the determinant operator V. T h e leading term is — M 0 £ = ~ ( & + M) 0 = ~ V = iJM 2 M o - 2 - ^ f & ~ as ~ o- (3-55) M Therefore, I ZQl + - = 0 (3.56) and 1 TM K ia + A : - + - (3 57) Substituting equations (3.55)- (3.57) into (3.54) and dropping terms of order 1 / M , (q ) k 0 3 and higher, we find H( =o) P ~ \j + J dkdp 7T+ S(k+p)it+ 2Mr II dkdp 0 + V (M ,{q (t)}) h 0 (3.58) k ^ S A * > *? - d i j / / d k d p ^ s { k + p ) { & ) 2 where Vh(M , 0 {q (t)}) k = M w i t h u)k being the eigenvalue of the mode 0 + ^Jdk fk(x). q+(t) q (t) k uo\ (3.59) Chapter 3. Soliton bremsstrahlung and the method of collective If we discard terms involving order of 1/M and 7r q (t) q (t) 0 k t k p p leading quantum correction to the classical result M 1 {P=0) « MI 0 K n {t), k p we obtain the — V[(f> (x)]: 0 H 75 n or higher and treat them m as negligible w i t h respect to terms of order M , q (t) q (t),and 0 coordinates n r+°o + - / dk{7v£it + L0 q+(t)q (t)}. (3.60) 2 k k The above H a m i l t o n i a n in the rest frame of the t o t a l field is the result we tried so much to achieve where the zero-frequency to no longer appears. T h i s expression, therefore, is enb tirely correct insofar as the small-amplitude and harmonic approximations are concerned; the wavefunction is t r u l y confined in the "directions" of the eigenfunctions f (x). k For eigenvalues P different from zero, we need not go back to equation (3.49) but simply invoke the Lorentz invariance so that Hp = \/p + 2 H (3.61) P=0 A c c o r d i n g to H a m i l t o n i a n (3.60), the approximate eigenstates i n the rest frame ( P = 0) w i l l have energies Kest,{n ) k ~ o M + E (»* + V 2 ) Wfc, \ (3.62) k Z where the asterisk reminds the reader that we deal w i t h quantum "excited" states. For field momenta P ^ O , once again, Lorentz invariance w i l l give us P = 7 E; e s U n k } V (3.63) and T = 7 E* (3.64) rest>{nk) as the t o t a l field momentum and energy, where 7= , : 1 V i - v 2 (3.65) v ' Chapter 3. Soliton bremsstrahlung and and the method of collective coordinates 76 V is the "speed" of the total field as determined by P V -7r=w = - (3 66) Practically, the above relations lead us to a situation where one would be unable to isolate the different variables, even though the relations remain self-consistent. In principle, for P^O, one can always move t o an inertial frame where P — 0, determine E * E S I | n f c p and finally boost the result i n order to get the answer i n the original inertial frame i n which P^O. However, for slow moving classical solutions 4> (x — X(t)), one can extract a H a m i l t o n i a n n in leading orders where the field operators P stands by itself. Indeed, we can make the following a p p r o x i m a t i o n 4 T = "T^rest + — 1 — rp* ^rest 2 I Brest , E *rest V 2 2 E* o ,T?*2 2 (E* E = rest E + P , r m + P ) 2 S T 2 (3.67) 2 E S T by use of equation ( 3 . 6 6 ) . Furthermore, E* 1 ^rest (E; + P ) 2 1 [E* 2 E S T * R E S T [(M + (P /E; )} 2 E S T + 6) + ( P / E ; 2 0 E S T where E* 4 E S T = M 0 + 5 We have remove {n } from the expression E* k eat with 5 < M 0 for simplicity. . )} ( 3 - 6 8 ) Chapter 3. Soliton bremsstrahlung and the method of collective Now, we can assume for low velocity X(t) P 77 that '« 2 coordinates MX, 2 (3.69) 2 so that P MIX 2 1 E* M +5 rest 0 Mo X6 2 2 MX = 2 lv± 0 Mo 2 MX + X 8. 2 (3.70) 2 0 Substituting equation (3.70) into (3.68) leads to 1 E* {Ef + P ) 2 est , M + <5 + M X 0 + X 5 2 2 0 ~ 1 ~ M + 5+ MX 2 0 0 1 5 M + MX (M + 2 0 0 0 MX) 2 2 0 1 M + MX 2 0 0 h -^- - 1 (3 71) W i t h this approximation, the total energy T becomes P 2 -X ) 2 2M ( 1 0 P PX ^r*«t + 7777" ~ ^2M T7^2M 2M 2 2 2 I r e s t 00 ( - ) 3 v 7 2 y 0 We shall suppose that the last term in the above equation is negligible compared to the other terms, because of the smallness of X , 2 is reasonable. and verify a posteriori i f such an assumption Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 78 Hence, we assume, for small P ^ 0, the following H a m i l t o n i a n in the small-amplitude and harmonic expansion: P 2M 1 r+°° 2 /-oo 2 R S p = + 0 + + »I- (3-73) M The above H a m i l t o n i a n concludes our program for quantization of fields about a classical solution. It w i l l be considered thereafter as the unperturbed H a m i l t o n i a n H 0 small perturbation Hamiltonians H are grafted. p 3.3 to which Quantization of the bremsstrahlung problem For a classical field i/)(x,t), the approximations we made to derive the low P - v a l u e d quantum H a m i l t o n i a n (3.73) from the exact one (3.49) can also be made to the classical H a m i l t o n i a n (3.40) to obtain Hp = 2M + 0 2 Loo ^ * ' 11 k 2W k for low values of the total field momentum P w i t h iTk(t) = 1 Ck |2> ' (3 74) c (t). k The proper expansion i n the moving reference frame attached to the i m p u r i t y is t/>(*,T) = il) {z - X{T)) a + (j)(z,r) (3.75) where +oo / dk c {t) f (z-X(r)). k (3.76) k -oo In section 3.1, we recall that we used instead the expansion ip(z,r) = ip (z) + (f)(z,r) s with = - ^ M z ) + dk c (t) k f {z), k Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 79 which proves a problem for the correct evaluation of the classical or quantum energy of the system. W i t h expansion (3.75), the perturbation H a m i l t o n i a n becomes Hp a f+°° — dz ib(z,r) 5(z+ v(r - = 7 = r )) 0 J—oo -MV(T-TO)-X(T)).+-7 where the term aib(x,t) f °dkc (t)J (v(T-r )-X(T)) 7 J-oo (3.77) +0 k k 0 8(x) making up the perturbation H a m i l t o n i a n density has been transformed to (a/7) ib(z,r) S(z + v(r — r )) under the Lorentz transformation (3.5). 0 U s i n g approximation (3.69), P « M X(T),.valid for slow moving solitons, the full H a m i l - 0 tonian reads as H = M + 0 (X(T)) + - 2 / dk {| 7r I +u + -IP (V(T-T )-X{T)) s 2 2 fe Ic + - 0 7 (3.78) I } 2 k f dkc (t)f (v(r-T )-X(r)). +CO k 7 k 0 J-00 It is a simple matter to derive the equations of motion for X(t) and c (t) from the above k H a m i l t o n i a n ; these are M X(r) 0 - T ) - X(T)) = - f (v(r 7 b 0 + - / 7 dk c (r) k f' (v(r k - T ) - X(T)) 0 (3.79) J-00 and C (T) + UJ C (T) = 2 P P P / ; ( ( _ ) - X(T)) . W T (3.80) TO These equations are different from the corresponding ones found previously (equs. (3.9) and (3.11)) by the presence of the additional terms —X(r) and (a/7) J dkc (r) k f' (v(r k — r ) — X(r)) 0 i n the functions' arguments in (3.79). Those additional terms destroy the 5 simplicity we had previously. It seems that the only way to solve those mixed equations 5 Expression (3.79), apart from the difference pointed above, is for small X(T) approximately the same as equation (3.9) once we make the identification Cfc(r)/M = —X(T), 0 where M A = 8, so that Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 80 is by self-consistent iteration. However, there is a lucky break. T h i s is provided by the fact that the soliton moves very slowly after its impact w i t h the i m p u r i t y so that we can expect I -%) I < v in such a way that V{T -T)-X{T) « V(T -T) 0 (3.81) 0 for all time r . W i t h this approximation, our equations of motion become M X(T) c p ( r ) + ^ c p rv = 0 ( r ) rv -f ( (r-r )) 7 b v + - f;(v(r = r+oo / 7 0 dk c (r) f' (v(r k - r )) k 0 (3.82) ./-oo - T0)) (3.83) T h i s last expression has the exact same form as equation (3.11) and, therefore, solving for Cp(r), gives the same answer as for c (r) i n section 3.1. Hence, we have not laboured p in vain in deriving that section's results. We know that those Cp(r) contain a to linear order so as to make the second t e r m in (3.82) proportional to cv . Since we are interested 2 only in effects linear in a because of its smallness, we can neglect this second term and recover essentially the previous expression (3.9). Thus, we have M X(T) « 0 - f (v(r 7 b - To)). ' (3.84) In the quantum mechanical domain, we make use of expansion + / *I>(Z,T) = il> {z-X{T)) a dkq (r)f (z-X(r)) k k J—oo so that the perturbed H a m i l t o n i a n H p H p = —I 7 = - 7 is dz I/J(Z,T) 8(z +• V(T - T ) ) 0 J-oo MV(T - T) 0 - X(T)) + - 7 / ^-oo dk q (T) k f ( (T k V - To) - X(T)) . (3.85) Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 81 The H a m i l t o n i a n of the problem is therefore H =H +H 0 P w i t h H = Hp as given by equation (3.73). 0 We can then determine the equations of motion for the operators X{T) and q {j) i n k the usual way by m a k i n g good use of the commutation relations (3.48). F o r X(T), we have ^ = -i[*(r),£(r)] 2M.0 dX(r) V dr ' ~ P , (3.86) x M 0 and — [P,1>MT-T )-X{T))] - 0 %-r )-I(r)) - ft 0 i dP - H X(T)) + •/-«> dkq (T)[PJ (v(T-T )-X(r))} k 0 ° d k J-oo 7 - r„) - / f 1 dx {fb(v(r - + dk U T ) ^ { ax v { r _ T o ) k 0 _ x { r ) ) HX q (r) f' (v(r k - k - X(r))| , r) 0 (3.87) so that substituting the former expression i n the latter yields ° ^ T M = ~\f (v(r-r )-X(r)) b + 0 dk q (r) f' (v(r k k - r) - X(r))|. 0 (3-88) For q {T), p we have ^ ( T dr ) = = -i[q (r),H(r)] p 7T, V (3.89) Chapter 3. Soliton bremsstrahlung and the method of collective coordinates 82 and — i [7T , H p ~dr~ [ T T P , q+(r) q (r) ] - = -i u = u? q (r) 2 p p — - p - p ^ f;(v(r i [TT , g_ (r) ] f- (v(r p p p - r) 0 X(r)) - r ) - X(r)), (3.90) 0 so that + < «P(T) 2 dT = - - ~ r ) - X(T)) . 7 0 (3.91) We note that the form of these operator equations is exactly the same as that of their classical counterparts, equations (3.79) and (3.80). In a like manner, we can make the same approximations on X(T) and the second term in equation (3.82) to obtain Mo—TT ar 1 d q (r) 2 « -fb{v(r-T )) = - - f;(v(r 7 (3.92 0 2 p dr +^g (r) P 7 0 - r)). (3.93) We remind the reader that the above equations and their respective classical analogues are valid in the l i m i t of a slow-moving soliton as viewed in the frame we designated as the soliton rest frame (whose term assumes its true meaning only for r —> —oo). In conclusion, by virtue of Ehrenfest theorem, the above shows that the evolution of the expectation value of the Heisenberg state | ip), which corresponds as closely as possible to our initial conditions of a soliton at rest and an i m p u r i t y dressing moving w i t h speed u a t r - 4 —oo, w i l l be governed for all intents and purposes by the classical equations of motion (3.9) and (3.11). Therefore, we can expect that, by comparing w i t h the classical energy expression (3.74), the average energy for the quantum excitation of the field to be also given by the classical bremsstrahlung energy H (r r > r ) , equation (3.18). In order 0 to prove this assertion, we had to take a lengthy detour i n the direction of the collective coordinates land; a trip, nonetheless, that w i l l prove most profitable in the next chapter. Chapter 4 Oscillations of a pinned soliton T h e previous two chapters cleared the ground for the main problem we seek to solve and which constitutes the raison d'etre of this study. T h e perturbation H a m i l t o n i a n density we are interested in is given by equation (2.39) w i t h V(x) = 5{x) (4.1) for the i m p u r i t y (or pinning center) located at the origin. A s in the last chapter, we divide the discussion between a classical and a quantum treatment. 4.1 Classical treatment T h e equation of motion flowing from the H a m i l t o n i a n density % = 7i + H 0 •0 - -0" + s i n ^ = a6'(x). p is (4.2) A s expected, we use the expansion of the field ij)(x,t) as follows: iP{x,t) = i> {x-X(t)) + (f>(x,t) s (4.3) where +O0 / -oo dkc (t)f (x-X(t)). k k W i t h expansion (4.3), the perturbation H a m i l t o n i a n yields 83 (4.4) Chapter 4. H a p = Oscillations of a pinned dx —^—J-oo soliton 84 Six) OX +CO / dk c (t) f' (-X(t)). k (4.5) k -00 If we assume from the outset that we deal w i t h slow solitons, then, upon using the unperturbed H a m i l t o n i a n Hp (3.74) and the approximation (3.69), the total H a m i l t o n i a n becomes A/f H = M . + ^X (t) 2 0 f+OO 1 + i / k +oo / | +u>l | c 2 dk{\* |} 2 k dk c {t) f' (-X(t)). k (4.6) k -00 T h e equations of motion for X ( i ) and c (t) from the above H a m i l t o n i a n are found to be p +O0 / c (t) p + uJ (t) 2 pCp = dk c (t) f' \-X(t)) k (4.7) k -oo -af;'(-X(t)) (4.8) in the l i m i t of a slow soliton. Equations (4.7) and (4.8) couple the coordinates X(t) resort to a self-consistent iterative method. consider first solving for X(t) and c (t). p We must therefore For the start of our iteration process, we by assuming c (t) = 0 i n (4.7). Under some assumptions to be defined later, the resulting solution X(t), upon substitution into equation (4.8), k w i l l generate solutions Cp(t) containing a to linear order which in turn, upon feeding back into (4.7) for the second iteration step, would make the term a J quadratic in a. dkc (t)f (—X(t)) k k Since we consider, as a first approximation, equations of motion only linear i n the small coupling constant a, we may as well ignore the higher order term in equation (5.7) and obtain instead Chapter 4. Oscillations of a pinned soliton M X{t) 85 = -af (X{t)). 0 (4.9) b Note that the Caldeira-Leggett formalism, which is taken up below, w i l l permit us to tackle the full equation (5.7) in the l i m i t of low amplitudes for X{t). Once again, (5.7) gives some form of Newton's law of motion for the coordinate of the soliton X(t) w i t h a force term —af (X). center-of-mass Indeed, this is exactly what b we would expect from a potential V ( X ) felt by a stationary soliton determined as follows: Consider the terms V(x;X) = I ( # ) + (1-cos ik) + <*^8(x) 2 given by the H a m i l t o n i a n density H for a stationary soliton ip (x — X) whose center is s located at x = X, a distance | X | from the impurity. We may always consider V(x;X) as the total potential density term as a function of the parameter X , the soliton's center. T h e corresponding potential is therefore r+oo +oo / dx V(x; X) -oo I-I {l^s( - )) dX +oo / x X 2 1 dx{-mx-X)f -oo + (l-cos^.(x-X)) + + a d M X X Q x ) 5(x) (l-cos^(a:-X))} Z ^ + a r+oo J-oo ^ { x - X ) dx ^ - 5(x) f / OX +O0 / = V(X) M + 0 dx f (x - X) 5(x) b -00 af {-X) b = M + 2asech(X). (4.10) 0 T h i s potential is drawn in figure (4.1). Notice that the force term in equation (4.9), —af (X), b gradient of V(X), —dV(X)/dX, is exactly the negative of the reinforcing even further the identification of (4.9) as a Chapter 4. Oscillations of a pinned soliton 86 V(X) x M -2lal 0 Figure 4.1: T h e physical potential provided by the pinning center. The potential for the domain wall, whose center is parametrized by coordinate X, is of hyperbolic secant form with M being 0 the rest mass of the wall and a a negative small coupling constant. Newton's law of motion for the soliton's coordinate X: Hence, insofar as quantum fluctuations are set apart, the effect of the pinning center is to provide a potential well of hyperbolic secant form to a moving soliton which may become trapped if it does not possess sufficient kinetic energy In order to solve equation (4.9), we suppose that we have a small amplitude X(t) the m i n i m u m of the potential V(X). In this case, we may expand the potential about V(X) and keep only the quadratic order so that dV(X) M X{t) 0 dX 8X(t) 2aX, X(t) (4.11) (4.12) where we used the mass of the soliton M = 8. 0 A simple solution to equation (4.12) is given by X(t) = X 0 sin(wTj) (4.13) Chapter 4. Oscillations of a pinned soliton 87 w i t h to = —a/4 (recall that a is negative). O w i n g to the fact that a is small, coX , 2 0 which gives the m a x i m a l speed of the oscillating soliton, w i l l have this one move slowly. E q u i p p e d w i t h this approximation for X(t), we can hope to solve for the coefficients Cp{t) by keeping terms only up to 0(X ). Let us see what we get for / * '(X(t)) 0 i n this soliton's low-amplitude approximation. r \x) p = -i f; + =/=p— P -ipX e = —f=—{p 2 + sech X - 2 zptghX}. (4.14) V 2 7 T LOp For small X, we have e vX 2 -i x P tghX sech X 2 w « 1 - ipX X, 1 - X. - 2 — , P (4.15) (4.16) 2 Note that the approximation for the exponential e^P X is valid only for sufficiently low values of p. If we look at our problem from the physics point of view, we would expect only the low-lying quantum energy states to be excited by the slow motion of the soliton at the bottom of the hyperbolic secant potential. Therefore, only states w i t h low momenta should make an important contribution to the quantum fluctuation. W e remind the reader that this whole calculation is done at temperature T « 0 so that the higher energy states w i l l hardly be populated at a l l . Therefore, we consider the approximation (4.15) as valid under those circumstances. Hence, w i t h the above approximations, we obtain /;'(-*(*)) * 1 + _ x (t) 2 + V 2 7 T LOp + i x(t)} P Chapter 4. Oscillations of a pinned soliton UJ 88 ip(p + y/2n \/27r • to order 0(X(t)), + 2)X(t) 2 p ( 4 1 7 ) UJ P so that equation (4.8), to this order, becomes c (t) v + = - O p ""•(•* + »>*(«). - VZ7T (4,8) V Z n Up Note that for a stationary soliton positioned right at the pinning center, equation (4.18) becomes, w i t h c (t) = 0 = X(t) —» 0, and X(t) p 2 T h i s value simply corresponds to the screening of the soliton caused by the presence of the impurity; in other words, this is the i m p u r i t y "dressing" in the presence of the soliton. Therefore, the screening is given by r+oo +O0 / dk c° k f (x) k -oo ptfAj*!. \J2lX J - o o Ul (420) k We can subtract off this i m p u r i t y contribution from equation (4.18) by setting c (t) P = c° + dp(t). (4.21) p Thus, we obtain + _ " P ( ^ ' > * ( « ) , ( 4 . 2 2 ) V Z7T Up or ~~Cp{t) + co Cp(t) 2 = - a ^ y/2n P + X sm(cot) 2) 0 ojp (4.23) Chapter 4. Oscillations of a pinned soliton 89 where we used equation (4.13) for X(t). Solving for Cp{t) is easy and the result is ~/s Cj,{t) = ae- " lu ^ i,,t + Be*""* t iaX p(uj + l) s'm(ojt) ^==-^ , _ >, , 4.24 2 0 V with A and B being two constants of integration. In order to fix those constants, we need the initial conditions. These w i l l be determined by the physical situation at time t =..— 7T/(2OJ) of a soliton at its m a x i m a l amplitude, x = —X , w i t h no velocity, about to be pulled in towards the m i n i m u m of the potential 0 V(X). For this case, obviously, 6p{-n/(2u)) = 0, (4.25) and from equation (4.18) 2 / , .pcpi-n/^)) f 0 \\ = au --j= = ^ iaX p(uj + l) 2 + p m y 0 + 1 (4-26) ) 3 W i t h initial conditions (4.25) a n d (4.26), we find f e = A ^2 -i7TWp/(2w) 7-1—5T-2 > {u) ~u )u 2 2 2 5 = 2 with Tp A iaX p{u)l 0 +1) P u l l i n g everything together, we obtain C PW = ° + Cp(^) C P T with A e "" p / " ( - ) 4 2 7 Chapter 4. Oscillations of a pinned soliton 90 We need not worry about the vanishing of the denominator (co — co ) since i t w i l l never 2 2 occur as to > 1 whereas to = —a/4 <C 1. 2 2 Therefore, r+oo too (f>(x,t) = / dk c (t) f (x~X(t)) k -00 k represents the classical small-amplitude fluctuation generated by the presence of a pinning center apart from the screening effect. T h e energy fluctuation is given by the classical expression i n equation (3.74), E } L = H S P ( M - + —^y . 0 "I r+oo = dk{\c (t)\ \ c (t) \ } 2 . 2 k (4.29) aX /-+A k (u + l) dk + ——2- / dk - j - ^ a — / 2 = + .co 2 k Z J — OO 2 2 2 47T 7 - A 47T J-A + 2 2 W | (wf - 2a; sina;t cos 2 \co . . . \ — + (u} + co ) sincot 4.30 {Co\ ^ A 4 2 A 2 k CO ) 2 2 {KCO k — \2co + * V \ ' CO . ) ) /7VC0 ' M l k sin —— + co t) \ \ 2co / J J , k co k where a momentum cut-off A has been introduced i n keeping w i t h our previous discussion of a limited number of modes involved i n the l i m i t of T —>• 0 a n d a slowly oscillating soliton. T h e first term in (4.30) corresponds to the i m p u r i t y screening effect alone. Also, for a slow soliton, we have P 2M _ 2 M (X) 2 0 0 and -too dx "rip r+oo p / -00 +oo / dkc (t)f' (-X(t)) k * + £ / > (a 2 ? r ~ J A [ k -oo 1+fc2( ^ 5 1)2 w (w* - w ) 2 k *»*l***+-ico-+•*«) [ w j* V 2a; fc / (4.31) Chapter 4. Oscillations of a pinned soliton 91 as the energy of the perturbation proper. Classically, we can speak of the eigenmodes fk(x — X(t)) being subjected to a time-dependent potential V(X(t)) i n <f>(x,t) of equation 4.4 as = 2a sech (X(t)), which destroys the simple harmonic motion that they would otherwise enjoy (in the l i m i t of a slow soliton), as can be inferred from equation (4.18) w i t h a = 0 and X(t) « 0. T h e above discussion was entirely classical; the next section deals w i t h the quantum treatment. 4.2 Q u a n t u m treatment In the last chapter, we expanded at length on the transition from the classical to the quantum world. We recall that such a transition needed a special treatment on account of the zero-frequency eigenmode fb(x) associated w i t h the translation invariance of the classical static solution tp (x — X). s However, here, i n the presence of the perturbation term a - % ) , we found that this led to a potential well 4.10 for which, classically, we expect the soliton to sit at the b o t t o m motionless, (cf. fig. (4.2 (a))) Therefore, the classical solution is no longer ip {x—X) s for any X, but rather ip (x) which s carries N O translational invariance. T h i s implies that we could attempt to diagonalize the operator ' ~2 (d U d<f> 2 a (4.32) 2 with U {4>) = a 1 - cos0 + d(f> a-^5(x) w i t h no risk of obtaining a zero-frequency mode. The classical picture (figure 1.2a) of a soliton at rest at the b o t t o m of the potential well is forbidden in quantum mechanics and Chapter 4. Oscillations of a pinned soliton 92 Figure 4.2: Q u a n t u m and classical oscillations of a soliton. The little filled dot is the virtual particle representing a magnetic soliton, the domain wall. In (a), the soliton is expected classically to be at the bottom of the potential well V(x); an effaced state that is strictly forbidden in quantum mechanics which, therefore, imparts the soliton in (b) with quantum oscillations about the minimum. In case of high energy states, which may come about in many different ways, the oscillations in (c) may assume a macroscopic character as befit classical states. we expect, therefore, to have small quantum oscillations about the m i n i m u m . Furthermore, as the soliton is very massive, these quantum fluctuations are believed to be minute. Since the potential well 2 a s e c h ( x ) is somewhat shallow, because of the small coupling constant a, we anticipate the eigenmodes {fb(x), fk(x)} of operator(4.32) to resemble our previous eigenmodes {f (x), fk(x)} obtained with the potential U((j>) = 1 — cos</>. b Therefore, it is reasonable to expect the coordinate c (t) associated w i t h fb{x) to reflect b the position of the soliton i n the same manner as X(t) d i d . In addition to this quantum effect for coordinate X (or i f you prefer c ), there also w i l l be quantum fluctuations for b the coordinate c i n the "directions" of the eigenmodes fk(x). k A l l this is very fine insofar as the exact eigenmodes {fb(x), fk(x)} flowing from the operator (4.32) can easily be found. Unfortunately, this might not be the case. However, once again, because of the smallness of a, we can think that expansion (4.3) w i t h eigenmodes fk(x — X(t)) w i l l not be a bad choice; the price to pay being the sad consequence of obtaining equations of motion m i x i n g the coordinates X(t) and Ck{t), as can easily be seen by inspection of equation (4.7) and (4.8), away from the simple harmonic form they Chapter 4. Oscillations of a pinned soliton 93 enjoy w i t h the proper coordinates {cb(t),k(t)}} Perhaps, this a better trade-off as the diagonalization of operator (4.32) can prove most difficult. There is another way i n which the original expansion (4.3) is preferable. T h e soliton may very well perform classical oscillations in the sense that the amplitude is large enough to be measured w i t h arbitrarily experimental precision (cf. fig. (4.2 (c))). Such a motion could be initiated by the application of a force for a very short interval of time even at T = 0. For a magnetic soliton, the Bloch wall, this force is obviously magnetic. In our case, the amplitude X Q for the harmonic motion (4.13) could be large compared with that associated w i t h the quantum fluctuations (cf. fig.( 4.2 (b))), and yet, small enough, on the length scale as determined in section 2.2, to develop the real potential sech (X) up to quadratic order only. In that case, expansion 4.3 can be regarded as one about a classical time-dependent solution ip (x — X(t)) w i t h X(t) s obeying the equation of motion (4.9). In keeping w i t h our special quantization scheme of collective coordinates i n chapter 3, we consider, therefore, the unperturbed H a m i l t o n i a n p2 H 0 = 2M 0 '+00 + dk t^k —oo + ^c c + f c } + Mj, (4.33) w i t h the perturbation term (4.34) # M ) = Mx~X{t)) + / dkc (t)f' (-X(t)). k k J—oo T h e equations of motion for X(t) 1 and c can be found as follows: k This is after all what diagonalization is all about ! (4.35) Chapter 4. Oscillations of a pinned dX{t) soliton 94 -i[X(t),H(t)} dt 2M [X(t),P*] 0 M 0 dP dt = -i[P,H] = -ia[PJ (X)} - b ia r+oo dkc (t)[P,f' (-X(t))] k k J—oo —a +O0 MOT dX dk r c,u,.n,dfL(-x) (t) /J-oo -OO dX +00 k w so that d X(t) 2 and dc (t) k dt = = d7r {t) k dt = -i[c (t),H(t)} k i[ck{t)>n n] k k -i[7T (t),H(t)} k = -i [% ,u) cl c ] - = -colc (t) 2 k k k k - ia[it ,ck k f'_ {-X{t))] k afZ'{-X(t)), so that + uolc {t) = -af*'(-X(t)). k (4.37) Chapter 4. Oscillations of a pinned soliton 95 W e note immediately that the operator equations for X(t) and c (t) correspond exactly to the form for the equations of motion of X(t) k ((4.36) and (4.37)) and Cp(t) found above, equations (4.7) and (4.8). Under the same assumptions of low amplitudes for X(t), so that and of omission of terms higher than linear order in a , we can approximate (4.36) by Thus, in the l i m i t of a slow soliton and small coupling constant a, the time evolution of the operators X(t) and c (t) m i m i c that of their classical counterparts X(t) k and c (t). As k discussed at the end of section 3.3, we expect, then, the average energy of the quantum fluctuation to be given by the classical result for the small-amplitude approximation, equation (4.30). Formally, the quantum fluctuation energy is given by the harmonicoscillator form of H a m i l t o n i a n (4.33) (minus the classical energy M ) , 0 X M 2 k z M x 0 2 2 . +00 r+oo / dk {n + l/2)u>k, k (4.39) -oo where X is the average velocity of the oscillating soliton. Expressions (4.30) and (4.39) 2 do not include the energy of the perturbation proper, equation (4.31). It is important to point out that some of those energies may contain infinities which should be handled with the greatest care to obtain a meaningful physical result. Here, for T = 0, the introduction of momentum cut-offs and rejection of the zero-point energy would seem to be sufficient for our purpose. In order to carry out the perturbation to higher orders, one would have to deal w i t h the usual techniques of normal ordering, counter terms, and the like, to handle the troublesome divergences. 2 Recall that we use the convention h = 1. Chapter 4. Oscillations of a pinned soliton 96 T h e next section deals w i t h the mapping of our problem onto the Caldeira-Leggett formalism and the calculation of the spectral density function J(u>), the key result of this study. 4.3 M a p p i n g to the Caldeira-Leggett formalism In the introduction, we mentioned that a general Lagrangian, which expresses the physical situation of a system coupled to its environment and whose quasi-classical equation of motion is given by MQ + nQ + = F (t), (4.40) ext is L = '-MQ 2 - V(Q) + QF (t) + J E K ' S ext 2 - rrijrfx)) (4.41) i For the case of strictly linear dissipation, we have i(Q) F = QQ, a linear coupling in the system's coordinate Q . It is a simple exercise to derive the equations of motion for the harmonic oscillators Xj in the strictly linear dissipation case: rrijXj + rrijOJ^Xj = —QCj y (4.42) or for unit-mass oscillators, x p + ulxp = -Q{t)C , p (4.43) showing the time-dependence of Q explicitly and changing the index for later convenience. Chapter 4. Oscillations of a pinned soliton 97 Recall the equation of motion we obtained for the coordinates Cp(t) in the classical treatment of the pinning potential i n the l i m i t of a slow oscillating soliton i n section 4.1 (cf. equation (4.18)) W ) + .1 cAt) = " V27T A p a r t from the screening contribution —auo /y/2ir, p "* p P + ). i ) (4.44) X{t VZTTUJp expressions (4.43) and (4.44) would be somewhat equivalent should the identification of Q(t) with the soliton's center coordinate X(t), Xp(t) with Cp(t), and C with — (iap p (a»p + l))/(\/2~7rOJ ) be made . One's first reac3 P tion would be to regard as suspicious the identification of a whole system parametrized by the coordinates {X(t), c (t)}, whose interaction with the surroundings is not accounted p for by the H a m i l t o n i a n (2.37), with the Caldeira-Leggett universe consisting, respectively, of the system and the environment coordinates, Q(t) and {x (t)}. p However, we must admit that the partition of the universe is always arbitrary and consonant with the kinds of questions one may wish to address. Furthermore, we showed in section 4.1 that, under some specific circumstances, the soliton's coordinate might be considered as a classical one subjected to a harmonic forced motion. It seems therefore reasonable to entertain the idea of the soliton's constituting the system (parametrized by X(t)) cou- pled to a bath of harmonic oscillators (parametrized by {cp(t)}) which may be viewed as the "environment." In our specific case, those oscillators correspond to magnons, the elementary quantum excitations of the magnetization field about the soliton solution. Above, we promised to deal w i t h equation 4.7 without o m i t t i n g the second term +O0 / dk {t) Ck m-x(t)). -oo Indeed, it is the presence of this term which makes it possible to think of the coordinate 3 The minus sign arises from the fact that we deal with complex coefficients c . p Chapter X(t), 4. Oscillations of a pinned soliton 98 the soliton's center, as being subjected to viscous forces reflected by the phenomeno- logical friction parameter n i n equation (4.40). In order to see this more clearly, let us write down the Lagrangian derived from the H a m i l t o n i a n 4.6 i n the small-amplitude approximation, L = M ^ o " ( V W X + !£ k {I I <*(*) | } I " (4.45) 2 2 k where V(X{t)) Kk{t) = af (X(t)) = c (t), + M b = 2asech(X(<)) + 0 M, 0 k and a discrete summation has been chosen instead of the integral sign for convenience. In the l i m i t of low amplitudes for X(t), approximation (4.17) is valid and so L transforms to L = M ° ~ " E V { X { t ) ) I " "l I k(t) | } + k k(t)f° - c k c 2 2 (4.46) c (t)c x(t) E k k where and C = _ H * M ± I 2 . !i ± t V Z7T T h e terms (1/2) { ^ | | c (i) | fe ( 4 . 4 7 ) V27T Wfc 2 + 2ac (t)f } k k lend themselves to completion of a square by the addition of a constant to the Lagrangian; which is immaterial insofar as energy differences and equations of motion are concerned. 4 U s i n g the fact that c* = C- , k k we find that 4 ? (£) This simply amounts to a canonical transformation for which the coordinates suffer a translation by specific constants [40]. 4 Chapter 4. Oscillations of a pinned soliton 99 does the trick so that L = M ^ 0 + \Y,{\ - V(X(t)) k 4 + o? Y, k - ]T q (t) C X(t) k k 9k{t) I Ck I ^{t) 2 |} (4.48) 2 X(t) k u where Qk(t) = c (t) + Zli k = + C k { t ) . V 2 7 T 10 U) k The last term in (4.48) can always be combined with V(X) to give an effective potential which takes directly into account the screening effect of the impurity V (X) - = V(X(t)) eff - i i ) yiix k (4.49) k 5 —X(t). co (4.50) k Thus, we obtain L = - V (X(t)) + | E eff (I ?*(<) I - fc I 2 I} w (4-51) 2 which fits a Lagrangian of the Caldeira-Leggett type (4.41) in the strictly linear case, with m k = 1 and no potential renormalization term Fi(X) E * "t ' 2 The equation of motion for the coordinate q (t) is given by p q (t) + u q (t) 2 P / P \ P Q (t) + u q (t) p c x(t), p ( « + 1) -X(t), V ^vr a>p W 9 2 P = P = . J (4.52) which is exactly the equation for dp(t) (4.23) that we found previously in the limit of a slow soliton and low momentum p, as could have easily been inferred from a comparison of equations (4.21) and (4.49). In the limit of low amplitudes, V(x) w UJX so that V (X) simply corresponds to a displaced harmonic potential centered around a value of X ^ 0 because of the constant force (a /y/2n) £ c /u . 5 2 EFF 2 f e k k Chapter 4. Oscillations of a pinned soliton 100 Similarly, should we derive the equation of motion for X(t) from Lagrangian (4.51), we would get the original equation 4.7 in which Ck(t) has been decomposed as (4.49) and f (—X(t)) k derived from the approximation for f' (—X(t)) in (4.17). k Thus, having cast our problem into a Caldeira-Leggett form, we can proceed with the application of various results associated w i t h this formalism; in particular, the calculation of the spectral density given by 2 y In our particular case, we have C p 5{u-u)j). rrij ujj (4.53) assuming the form in (4.47) with m p — 1 and the summation replaced by an integral on account of the continuous variable p. Hence, we derive J(u) = dp 5(v - u u) p n '•+ , p {ul + l ) C? f+°° 2 = 00 T L Y Ji * i P I = 2 u) p 2 2 2 ( 5 u - v) u 2 8(y - Ui 1 2 + l) p n 2 5{y- 2 {u - l ) / 2 2 K - I ) K + I) , — -jr n fifY(u / + J(u) u.v up du) O? o o — ' du a 2 {v + l ) 2 Up) 2 1/2 v > 1 (4.54) V < 1 . Let v = 1 + x. For v —> 1, we have x —> 0 so that J(i/) EE J ( x ) ^ 2V5a x / . 2 1 2 (4.55) E q u a t i o n (4.54) is the major result we sought to derive in this study. A graph of the normalized spectral function is shown in figure (4.3) along with its l i m i t i n g behavior in Chapter 4. Oscillations of a pinned soliton 101 the region of v « 1 as reflected by (4.55). It is very important to bear in m i n d that this spectral density is for the case T = 0. T h i s function gives the key to discussions of dissipation and its effect on the motion of a magnetic soliton, a B l o c h wall, at temperature T « 0. T h i s dissipation, we recall, is generated by the pinning potential which, indirectly, couples the magnetic harmonic oscillators (the magnons) to the quasi-classical motion of the B l o c h wall. J(v) frequency v Figure 4.3: T h e normalized spectral density function J(v). This is the normalized spectral density function at T = 0 for the interaction between magnons and an oscillating Bloch wall at low amplitudes. Those magnons are not the typical spin-waves, but rather the elementary quantum excitations built in the presence of a domain wall. The function given above is normalized in the sense that the factor a /2 2 in the real spectral density function J{v) (equation (4.54)), where a is a small coupling constant parameter, has been dropped in J(v). In the inset, the behavior of the function near the frequency v = 1 is shown magnified. In passing, we mention that, should the dissipative mechanism be quasi-linear in the sense given by Caldeira-Leggett [12], then the phenomenological friction coefficient n would have the following relation j(Q) dQ dF V 2 y ^n~u7) 5{v-u)j), (4.56) Chapter 4. Oscillations of a pinned soliton 102 using the notation of expression (4.41). L o o k i n g at our specific Lagrangian (4.45), where Q «->• X, Xj '<->• Cfc, and m = 1, we note that F (X) k £ J—oo K n o w i n g that f {X) k n(u) = -{k r°°dfc = 4 [./-oo so that k (4.57) tO + 2 sech x)f (X), 2 we obtain 2 k 4 8(* ~ "*) W = a f' (-X) k + 4 s e c h X t g h X f ™ dk 4 2 + 5 { ~ y 7-oo fe + S p ech>* ( t t to k *'(4^-*') 8{u - to ) k J-00 2 \Ji 8{v-u ) 2a 2 4 ( 2 v 1)1/2 V (l/2 _ sech X 2 _ 1)2 2 ^ ' ^ a; ^ u ; - 1 + k k + r + 4 s e c h X t g h X [ °°dto * k to + 3 2 X 4 s e c h . 4 gh X 2 X t ^ y > 1 u2^(^-l)(3^ 2 + l) + sech^X- - 0 (4.58) v<\. T h i s calculation of n only assumes a slow moving soliton in contrast w i t h that of J{y) (4.54) which, in addition, also assumed low amplitudes for X(t) c (t). p In the l i m i t of low value for X(t), f 77(1/) -> { a 2 2 0 (v 2 - l) and low momenta for expression (4.58) takes the form 2 r , 2 , 2 2 l " (4.59) i/ < 1. T h e key point to note is that the friction coefficient n still vanishes for u < 1 as J(i/) does. Hence, the conclusion to be drawn i n the next chapter w i l l remain unaffected. Chapter 5 Discussion 5.1 Interpretation of the spectral density function J(u) T h e interpretation of our specific form of J(y) is very interesting. For driving frequencies v less than one, we expect the oscillating soliton to experience no friction as J{u) = 0 in this domain. In other words, the magnons do not interact w i t h the soliton at such low frequencies. A s the natural frequency of the system can be thought of as oo = (the stable frequency of the potential V(X) = 2 a sech (X)) y/—a/2 and a is already a very small coupling constant, we think likely that the driving frequency v, which would cause resonance, w i l l fall w i t h i n the dissipation-free domain. In practice, as v —> oo, we should expect a resonance peak at u & oo w i t h virtually no linewidth ! Indeed, i n appendix C , we show that Lagrangian (4.51) admits the lifetime r n energy level E n for the harmonic oscillator X(t) of an coupled to a bath of harmonic oscillator {qk} to be expressed as " J M = where oo is the natural frequency of the oscillator X(t). 0 (5-1) ^ ~ For our specific case, oo = oo < 1 0 so that, in principle, the lifetime of each energy level of the soliton, parametrized by X(t), is infinite. In other words, whatever small oscillatory motion the soliton may possess, it will not be damped by the magnon environment at T = 0. In section 1.5, we discussed the various sources of dissipation beside that of magnons. We have suggested that, at T = 0, only the dissipative effect of a spin bath and possibly 103 Chapter 5. Discussion 104 of phonons could be considered important on the motion of a domain wall. O u r result for the spectral density (4.54) arising from magnons at T = 0 vindicates this suggestion by ruling out any dissipation whatsoever for walls oscillating at a sufficiently low frequency. Therefore, the quantum tunnelling of domain walls at zero temperature under a pinning potential barrier can only be possibly suppressed by nuclear spins, magnetic impurities or phonons in the presence of defects. 1 Hence, magnons, at T — 0 to the very least, should not be singled out as the main cause for the suppression of quantum tunnelling for states that are macroscopically distinct from each other, such as the case is for a domain wall pinned at one location and depinned at another on a macroscopic scale. 5.2 W h a t could be done next? The first obvious extension of this study would be to investigate the dissipation for the system at non-zero temperature. In so doing, one might want to adopt a treatment of the zero-frequency modes that would be amenable to the path-integral formalism since the latter offers a straightforward generalization to quantum field theory at temperature T ^ O and may, i n some cases, be more gentle on the well-known ordering ambiguity i n going from classical to quantum formalism. There exist adaptations of the collective coordinates method to the functional-integral formalism; the reader may refer, for instance, to Rajamaran and Weinberg [46] and Lee [41]. A t T 7^ 0, we may have to worry about an additional degree of freedom besides the direction of the magnetization vector M , namely, its magnitude M which may very well differ form its saturation value M 0 expected at T « 0. Also, one could m i n d taking into account the scattering of solitons among themselves; a behavior which is unlikely to be present at T « 0. Furthermore, one could consider higher-order kinks or breather modes, 2 1 2 The effect of nuclear spins is important despite the presence of defects. Even though the scattering amplitude is time-independent, there is a higher likelihood of collision Chapter 5. Discussion 105 all of which translates into more massive, and therefore more energetic, solitons which a non-zero temperature could "produce" and support. In any event, unless one considers relatively high temperature T in the sense that ~ ' / E k T e is non-negligible for energies E' associated w i t h the above phenomena, these would constitute higher-order corrections to the magnons' effect on the motion of a simple soliton. 5.3 Conclusion We recall that the goal of this study was to give a qualitative idea of the dissipa- tive effect of magnons on the macroscopic quantum tunneling of a B l o c h wall pinned a defect. W e have effectively shown above that magnons do N O T interact dissipatively w i t h such a domain wall which oscillates at a low frequency at temperature T = 0. In order to demonstrate this, we needed a method, the collective coordinates, to eliminate the zero-frequency eigenmode fb(x) which rendered invalid, in the classical or quantum description, the approximation of the low-lying energy levels of a domain wall by a set of harmonic oscillators which included this mode. Yet, the eigenmode fb(x) was shown to play a crucial role in the establishment of a Newton-like law governing the center-ofmass motion of a domain wall i n the presence of impurities. Also, we cast our problem in the Caldeira-Leggett formalism in order to apply the usual techniques related to the suppression of macroscopic quantum tunnelling by a microscopic environment. One obvious extension of this study was shown to be the case at non-zero temperature which corresponds to the real situation in which experiments are conducted. Unless the dissipation arising from a spin bath and from phonons i n the presence of impurities are resolved in detail, the relative importance of magnons on the quantum tunnelling of a domain wall w i l l be difficult to assess. For the same reasons, a comparison between the theoretical predictions and experimental results cannot be fully realized. events at higher temperature. This new feature considerably complicates the statistical behavior of the quasi-one-dimensional system; preferably, we strive to obtain magnetic wires with one domain wall only. Chapter 5. Discussion 106 In closing, we mention that the existence of quantum tunnelling at a level as macroscopic as a domain wall provides a strong argument for the validity of quantum mechanics in the classical world. A l t h o u g h a more convincing test in this regards would be the observation of M Q C , a greater understanding of the effect of dissipation on quantum phenomena could surely indicate us on how "macroscopic" a scale we should ever hope to observe a quantum tunnelling event. Should experimentalists repeatedly fail to observe the quant u m tunnelling of a macroscopic object which theory allows, after all factors have been exhaustively taken into account; we might witness either the breakdown of quantum mechanics at this object's macroscopic level or the inadequacy of our theoretical description of quantum dissipative processes. In any case, this would certainly rekindle the con- troversy around the quantum measurement problem and, perhaps, prompt physicists to entertain the idea that entirely new irreducible levels of reality may emerge at certain stages of complexity of a system. However, reductionism have served us so well for the understanding of Nature that it truly deserves our vote of confidence for her investigation. T h e author have nevertheless the solid impression that Nature, somehow, will never allow to be so easily contained w i t h i n the confines of a reductionistic receptacle; her true beauty lying i n her mysteries. Bibliography M . Bunge, Causality, T h e Place of the Causal P r i n c i p l e i n M o d e r n Science (The W o r l d P u b l i s h i n g Company, Cleveland, 1963). R . Descartes, i n Oeuvres Philosophiques, edited by F . A l q u i e (Garnier Freres, Paris, 1963), tome I, p. 315; i b i d . , pp. 479-480. Encyclopaedia corpus 19. Universalis (Encyclopaedia Universalis France, Paris, 1989), A . J . Leggett, i n The Lesson of Quantum Theory, edited by J . de Boer, E . D a l , and O . Ulfbeck (Elsevier Science Publishers, A m s t e r d a m . 1986), p. 35. M . Jammer, The Philosophy of Quantum (John W i l e y & Sons, New York, 1974). Mechanics N . B o h r , i n Albert Einstein: Philosopher-Scientist, edited by P. A . Schilpp (Library of L i v i n g Philosophers, Evanston, Illinois,1949), pp. 221-222. J . von Neumann, Mathematische Grundlagen der Quantenmechanik (Springer, B e r l i n , 1932). [translated i n English by R . T . Beyer i n Mathematical foundations of quantum mechanics (Princeton University Press, Princeton, N . J . , 1955)] F . L o n d o n and E . Bauer, i n Actualites scientifiques et industrielles: Exposes de physique generale, edited by P. Langevin (Hermann, Paris, 1939), 775. A . J . Leggett, i n Quantum implications: essays in honour of David Bohm, edited by B . J . Hiley and F . D a v i d Peat (Routledge & K e g a n P a u l , L o n d o n , 1987), p. 85. J . M . M a r t i n i s , M . H . Devoret, and J . Clarke, Phys. Rev. Lett., 55, 1543 (1985); S. Washburn et al, Phys. Rev. Lett., 54, 2712 (1985); D . B . Schwartz et al, Phys. R e v . Lett., 55, 1547 (1985). A . J . Leggett, i n Quantum Tunnelling in Condensed Media, edited by Y . K a g a n and A . J . Leggett (Elsevier Science Publishers, A m s t e r d a m . 1992), p. 1. A . O . Caldeira, A . J . Leggett, A n n . Phys. ( U S A ) , 149, 374 (1984). A . J . Leggett, Phys. R e v . B , 30, 1208 (1984). 107 Bibliography 108 [14] R . A . Serway, C . J . Moses, and C . A . Moyer, Modern Physics, Sunburst Series (Saunders College, Philadelphia, 1989). Saunder Golden [15] R . P. Feynman and F . L . Vernon, Jr., A n n . Phys. ( U S A ) , 24, 118 (1963). [16] C . K i t t e l and J . K . G a i t , in Solid-State Physics, edited by F . Seitz and D . Turnbull (Academic Press, Inc., New York, 1956), v o l . 3, p. 437. [17] A . P . Malozemoff and J . C S l o n c z e w s k i , Magnetic Domain (Academic Press, Inc., New. York, 1979). [18] C . K i t t e l , Introduction 3rd edition. Walls in Bubble Materials to Solid State Physics (John W i l e y & Sons, New York, 1966), [19] W . D o r i n g , Z. Naturforsch., 3a, 373 (1948). [20] J . Rubinstein, J . M a t h . Phys., 11, 258 (1970). [21] A . Barone et ai, R i v . Nuovo Cimento, 1, 227 (1971). [22] L . D . Faddeev and L . A . Takhtadzhyan, Teor. M a t . F i z . (Russia), 21, 160 (1974). [Theor. M a t h . Phys., 21, 1046 (1974)] [23] P. M . Morse and H . Feshbach, Methods of Theoretical Physics, Part II, International Series in Pure and A p p l i e d Physics ( M c G r a w - H i l l , New York, 1953), p. 1651. [24] P. C . E . Stamp, Phys. Rev. Lett., 66, 2802 (1991). [25] P. C . E . Stamp, E . M . Chudnovsky, and B . Barbara, Int. J . M o d . Phys. B , 6, 1355 (1992). [26] O . Iglesias, E . M . Chudnovsky, and P. C . E . Stamp, Phys. Rev. B , 46, 5392 (1992). [27] A . G a r g and G . - H . K i m , Phys. Rev. Lett., 63, 2512 (1989). [28] A . G a r g and G . - H . K i m , J . A p p l . Phys., 67, 5669 (1990). [29] A . J . Leggett et ai, Rev. M o d . Phys., 59, 1 (1987). [30] N . V . P r o k o f ev and P. C . E . Stamp, unpublished. [31] B . B a r b a r a et al, J . A p p l . Phys., 73, 6703 (1993); M . Uehara and B . Barbara, J . Physique, 47, 235 (1986). [32] M . Uehara et al, Phys. Lett. A , 114, 23 (1986). [33] E . M . Chudnovsky, J . A p p l . Phys., 73, 6697 (1993). Bibliography 109 [34] K . H o n g and N . Giordano, unpublished. [35] H . J . Mikeska, J . Phys. C : Solid State Phys., 11, L29 (1978). [36] R . Rajamaran, Solitons and Instantons, A n Introduction to Solitons and Instantons in Q u a n t u m F i e l d Theory (North-Holland, Amsterdam, 1982). [37] G . Eilenberger, Z . Phys. B , 27, 199 (1977). [38] M . B . Fogel et al, Phys. Rev. B , 15, 1578 (1977). [39] R . P. Feynman, Statistical Mechanics: A Set of Lectures, Frontiers in Physics, Lecture Note Series (Addison-Wesley, Redwood City, California, 1972). [40] H . Goldstein, Classical Mechanics 2nd edition. (Addison-Wesley, Reading, Massachusetts, 1980), [41] N . H . Christ and T . D . Lee, Phys. Rev. D , 12, 1606 (1975). [42] M . Creutz, Phys. Rev. D , 12, 3126 (1975). [43] J . F . Janak, Phys. Rev., 134, A411 (1964). [44] J . M . W i n t e r , Phys. Rev., 124, 452 (1961). [45] D . Esteve, M . H . Devoret, and J . M . M a r t i n i s , Phys. Rev. B , 34, 158 (1986). [46] T . D . Lee, Non-Topological Solitons ( C o l u m b i a University Preprint C o . , 1976), p. 2271, Lectures given at the Symposium on Frontier Problems in H i g h Energy Physics, Pisa; R . Rajamaran and E . J . Weinberg, Phys. Rev. D , 11, 2950 (1975). Appendix A Detailed derivation of equation (2.10) a n d (2.13) In this appendix, we provide the details for the derivation of equation (2.13) i n chapter 2. In equation (2.8), we have the following terms sin0 m sin0 , (A.2) sin# m sin# _i, (A.3) m + 1 m (f> +i - <Pm, and m 4>m ~ (A.4) 4>m-l- (A.5) O n account of the assumption that the direction of the magnetization vector varies weakly over the lattice spacing a, we would like to express the above equations i n terms of quantities defined at the m t h site only. Starting w i t h (A.3), we have Pm+1 - <Pm ~ ~fa i a where the derivative is taken at the point z = ma, i.e., the coordinate of the m t h ion. F o r (A.4), we have to take into account the second derivative at the point z = ma as shown below: Let cj) m - A0 _x m A0 <£ _i = = A 0 m — = A 0 m H m T O _i. Then, Az, Az — where Az = ma—(m — l)a, a (A.6) z=ma • 110 Appendix A. Detailed derivation of equation (2.10) and (2.13) 111 A n d so, we obtain A A A<^m-1 a A0 _i m (m — 1) a in fi — A(ft a A(A0) ma — (m — 1) a Az 2 m ma — (m — \)a z=ma , using equation (A. 5), 0m(<). - ^ m - i W -> A ~Q^ ~ Q^- A ( A- 7) In like manner, proceeding with ( A . l ) , we get sin e m+l « sin 0 m + cos 9 (6 m - 9 ), m+1 m so that sm6 (t) m smO (t) m+l « sin 6> (£) + sin0 (r;) cos9 (t) -)• sin 0(z,7j) + sin 6(z,t) cos 6>(z, t) — ^ ' ^ <9z 2 m m m (O +i(t) - 2 1 - (9 (z,t)) s 2 + 0 (<)) m m - 0 (z,t)^^-, (A.8) s a az where we used approximations (2.8) and (2.9) applied to the continuous case, cos0(z,7j) « r(z,0 sinc?(z,t) « 1 - (A.9) (0*(z,/;)) , (A.10) 2 and discarded variations i n (A.7) of order two and higher. For s i n # _ i i n (A.4), we m obtain in the continuum limit sin0 _iOO m -> sin 6(z,t) - a cos 9{z,t) + a cos 9{z,t) OZ so that , 2 _ 0Z Appendix A. Detailed derivation of equation (2.10) and (2.13) s'm0 (t) m sin0 _i(£) m — > sin 6(z,t) d9(z,t) cos9(z,t) dz - a sin6(z,t) 2 112 d 6{z,t) + a s'm9(z,t) cos 9(z,t) dz d9(z, t) 1 - (9 (z,t)) + 9 (z,t) ( - a + a dz dz 2 2 2 « s 2 s 2 where we neglected (9 ) . s 3 Also, using (A.8) and (A.9), we obviously have s'm'9 (t) c o s 0 ( i ) sin6> (t) s i n 0 ( t ) m m m cos 0(z, rj) sin (b(z, t) (1 - 2 (9 (z, t)) ). s m 2 Substituting all those relations into (2.8) yields in the continuum limit hy/S(S + l)h us(S 2 dt + i) (1 - (9 ) ) s 2 H 0z a e + a9 dz*)- a ' e dz L u dz dz e - e dz ) a + 2 i a e d z - i a e d z dz e^ dz a dz ' ^ dz ) i ; 2 dz > - a 2h S(S + l)K cos4>sm(b(\-2(9 ) ) 2 s ^ 2 e a dz a dz ) 2 + itfy/S(S + 1) K 9 2 S y Z (A.ll) ift J 2 5 ( 5 + 1) {(I-(9 ) ) s 2 '.{ d(f> d cf>\ -2h S(S + l)K cos(bsm(j>(l-2(9 ) ) 2 / d<f 2 2 s + ih y/S(S + 1) K 9 2 2 y s z = n a J5(5+l)|(l-2(^) )0' 2 2 2 (A.12) 2h•S(S + 1) K cos 0 sin cf> (1 - 2(0*) ) + ih ^JS(S + l)K 9 , 2 y 2 s z Appendix A. Detailed derivation of equation (2.10) and (2.13) 113 where we discarded in (A.11) the last two terms w i t h i n the parentheses for they involve variations of order two and higher, which which we neglect. If we also neglect the term (9 ) s w i t h respect to one in expression (A. 12), we obtain 2 (2.10) 7 ^ ( 5 + 1) — « h Ja S(S 2 + l)-^ 2 - 2h S(S + l)K 2 cost sin<f) y + ih y sJs~+V)K 9 . 2 f (2.10) s z Let us now calculate the time evolution of the operator S+(t). § - -11^*1 + K i {S+ S* z = iJ{—S+S^+i m — S^-iS^ + K z ^ ( 5 + l ) | s i n 0 m e ^ iJ S(S + 1) + S° S+} + K m + S^S^-i i {5+ + S S* + S + S+} z m } z y m m (A.13) 'S'm+l ^ m ) + K {S^ y + S^J , = j - sin 9 + e^ l m cos 9 m cos0 + K y sin 9 m sin0 _! ^ ~ l m + i K S(S Z - m+1 m m l e e^ l e^ + m m + 1 (A. 14) m cos0 j 1 m cos0 m + cos0 m sin</> c o s # m + cos# m w m m m + sin0 + 1) { s i n 0 e ^ S(S + 1) { s i n # cos 9 ^ m sin0 e^™} m sin# m sin</> }. Sandwiching equation (A.14) w i t h a classical state vector as outlined i n chapter 2, we obtain: m Appendix ft^5(5 A. Detailed derivation +l)|(sin0 e^) m of equation (2.10) and (2.13) 114 = ih JS(S5 ( 5 + + 1) ^ |- e ^ 2 sin0 m m (cos0 i + cos0 _i) m+ + cos9 (e ^ (A.15) m s i n 0 _ ! + e^m+i i m m + 2 z ft iv~ 5 ( 5 + 1) sin 0 m + 2 / i A'j, 5 ( 5 + 1) sin 0 cos 0 2 z 2 m s in0 cos 0 m + 1 e^ l( m )} m sin 0 . m m Concentrating on the first term on the right-hand side of (A.15), denoted as X, we make use of the following expansions cos (9 i TO+ cos0 _i m ss cos 0 m - sin 0 m (0 - 0 ) m + 1 m —> cos0(z,t) — s'm9(z,t) 90 —a, oz « cos0 - 0 ) -» cos0(z,i) - m - sin0 (0 _i m m (A.16) m sin0(z,t) ( _ a ^ + a 3^2)' to obtain i n the continuum limit X = i f t J 5 ( 5 + 1) j - e ^ ' ^ ) s i n 0 ( z , t ) 2 ^2cos0(z,<) - s i n 0 ( , M ) a 2 cos 0(z, t) + e y + = ih JS(S 2 e i ( ^ , t ) + a f ) f ^cos0(z,t)] + l ) e ^ { - 2 sin 0 cos 0 + s i n 0 a 2 + COS0 a 2 2? ^(-a|^ + a f | ) ^ i a ^ \ . 1 9« / + e 3 « 1 sin 0 2 v 2 a Appendix A. Detailed derivation + l e \ dz of equation dz ) 2 (2.10) and (2.13) 115 [ -a cos 9 — + a cos 0 — + e \ dz dz J 2 m a — cos 9 dz dz 1 Once more, we discard variations of order two or higher that would result by expanding the exponentials and obtain X « ih JS(S+ j-2 l)e ^ 2 i( cos 9 tt JS(S {e a ^ 2 + l)e^ 2 s sin 0a 0 2 sine? cos0 + (2 + a 2 + a 2 2 9 2 ( /\ . „ sine? - ^ 2 (A.17) + a d9 — 2 o cos0 2 n 0} where we used (A.8) and (A.9) and neglected 9 . d (j) We assume that a term like 9 -^-r constituted of two small quantities is negligible so dz that (A.17) reduces to Z S 2 s 1 ifiVS^ + l j e * - ' ^ 2 ^ . oz (A.18) 1 In the continuum l i m i t , the second and t h i r d terms on the right-hand side of equation (A.15) simply transform to 2ih K S(S + l)e ^ ^ 2 i z + 2k 2 K S(S y z sin 9(z,t) cos 9(z,t) + 1) s i n 0 ( z , t ) cos 9{z,t) sin (j>(z,t) « 2ih K S(S + l)e ^ ) 2 i z + 2hK 2 y We mentioned i n chapter 2 that K y t 9 (z,t) s (A.19) S(S + 1) 9 (z, t) sin <f>(z, t). 8 is much smaller than K ; using this fact and knowing z that sin</>(z, t) have the same order of magnitude as e^ ^\ z it is legitimate, therefore, to approximate (A.19) with 2 ih K S(S + 1) e ^ ( ' *) 0,(z, t). 2 z z (A.20) Appendix A. Detailed derivation of equation (2.10) and (2.13) 116 T h e left-hand side of equation (A.15) yields then the following i n the continuum limit: hy S(S / + l)|(sin0 m e^) -4 hy S(S + 1) ^ = nJ5(5 + l ) ( c o s ^ , t ) ^ + i s i n ^ . ^ e ^ ' * ) ^ / (sin0(M).e^ ( M ) ) ihy/SiS + VeW*'*)^, where we neglect terms like # (A.21) and 0 — because they are composed of two small 2 s quantities. Therefore, using (A.20) and (A.22), equation (A.15) i n the continuum l i m i t reads as ih y/S{S + l) ^ 5(5 + 1) a H « ih J « ft 5(5+ 2 + 2 ift X , 5(5 + 1) 2 2 IJa 0 t) 2 7^5(5 + 1 ) ^ 2 1) + 27v" 0 | z s (A.22) Taking the classical limit (ft —>• 0, 5 —>• oo) on both sides of (A.24) leads to equation (2.13) where s is the value of the classical spin, i.e., the magnitude of the magnetization. Appendix B Inverse of the velocity matrix D T h i s appendix is not meant to be rigorous, but gives a heuristic derivation of the inverse of the m a t r i x D introduced i n chapter 3. One should realize that m a t r i x D is one indexed by a continuum, the variable k, and therefore, it may not obvious, or, for that matter, even true, to generalize the method we adopt in this appendix (valid for the usual m a t r i x w i t h a discrete index) to a continuum. A discrete indexed version of m a t r i x D is the following <lbb &b,-n 0-b,n 1 o (B.l) ^0,6 V a n, o y, i b where a j = a^,, (a _ )* = a ^ , a n d a it h) \ fe kjP = 5„. k W e let i, j denote the total index p set including the index "6", whereas we reserve k,p for the numeral indices proper. T h e special case n = 1 is shown below where, for simplicity, we denote A by A only; x ^ bb a 06,-1 0>bfi b,l a G-1,6 0 0 1 «0,6 0 1 0 1 0 0 117 ^ Appendix B. Inverse of the velocity matrix ^ a D a -i bb b 118 a bi0 'a 0 ) 1 0 0 1 b,0 0 1 0 \ a>b,i 1 0 0 a ^ (B.2) T h e real m a t r i x D in chapter 3 is obtained from two successive l i m i t i n g transformations: first, the index k running from — n to n has n —> oo to obtain an infinite countable index set, and second, is made a continuum to obtain an uncountable index set. We w i l l not show the generalization of our method for either of these two l i m i t i n g transformations. We know that the inverse of a matrix C is given by C' = 1 (B.3) DetC where Cn = DetM {j and Mij is the m a t r i x obtained from C by deleting the i th Let us compute the inverse of m a t r i x A n t h column. w i t h the help our special case, m a t r i x A . F i r s t , let us count the numeral indices as k -n -n + 1 so that for (B.2) we have row and j k' <-> 1 <-> 2 Appendix B. Inverse of the velocity matrix D 119 k k' -1 <-> 1 0 <-> 2 1 O 3 and obtain respectively °66 0-6,1 ••• Ofi.n+l 1 0-6,1 : O (BA) 1 b,n+l a ^ 0&,2n+l o 1 and O-bb A 0-6,1 «6,2 06,3 06,1 0 0 1 0 0 1 0 1 0 0 6 ) 2 \ o Now, denoting A n M by C for convenience, we have O C bb = £>e* O (B.5) Appendix B. Inverse of the velocity matrix D A l s o , C ,k> = (-l) +(*'+ ) DetM , , 1 with 1 b b k ( a ,_ 6 M,, = M b k 120 = b:k 0 0 ... 0 1 ^ n 0 0 ... 0 0 a ,_* 6 \ 1 0 ... 0 0 ) a, b n where the resulting m a t r i x has its first column made up of a l l the a btP zero following only the row starting w i t h the specific value a ^ . b> k w i t h a string of It is, then, a simple matter to see that DetM , = btk ( - l f + 1 a , (-1) byP w i t h p' being the value corresponding to p = — k. Thus, = C, b k Now, for C k:P ( - l ) ^ ' ) ( - l ) ^ a = °>b,p' > = ab k v • (B.6) of numeral indices k, this involves a bit more bookkeeping. For this reason, we w i l l explicitly deal w i t h the specific case of n = 1, i.e. the m a t r i x A i n (B.2), in the hope of finding a general pattern. We have A, k p = A ,, , k p = (-lf ' +p +2 DetM , . k y We need calculate only the six different sub-matrices M / > derived from A. k yP Appendix B. Inverse of the velocity matrix D ( a ao bb Dei M n = Det b o 0 0 6 ) 3 ^ 3 1 db,2 V 121 0 J (a , ) ( a ) b 3 M K ) > 2 or 3 Det M _ i _ i -Ki) w i t h the original indices fc. 2 ^ a Dei M x2 = Dei a aA a 0 0 B 6 j 2 \ «6,3 «6,2 a&,3 , Dei M _ i j 0 = a 6 ; 0 a ^ 6 3 0 0 1 0 or M ^ ^66 flft.l ^6,2 ^ DeiMi j 3 = Dei -o D e i Af.i,! = G6,2 0 1 V 06,3 1 0 + a i a bb + a ,-i a \ + -a 66 by b -06b + | o 2 = Det 06,1 V 06,3 (a ) 2 b> | M ^ abb o DetM ,2 + a6,2 «6,2 , 6 i 3 2 M 0 1 bfi + (a a bi3 1 0 i ) 0 ^ ) 2 or Appendix B. Inverse of the velocity matrix D Det M , 0 0 122 = -a + 0-6,1 o = -a b + a _i a i + a ia _i bb b b o ,i Det = \ b ) 66 06,1 06,2 06,1 0 0 M 1 0 a De*M 3 > 3 = ( a a Det a ,i 06,1 Oj,,2 \ bb = { b,2 a lfl / 06,-1 O£,,o b DetM \ or b 0 b 2 a ,i 0^,2 , 73e<M ,i or b 6 a 2 3 b 6 ) — —abb + 2 DetM , + a ,3 a ,i , 6 j 3 = —a = -(o ) , 6)1 0• 0 0 1 ai bt or 2 M = / -{a -if. K We note that the determinant of m a t r i x A is given byDei A = -a bb + a ia Det A = -a bb + a _ a = ~a so that Det M _ a ) 1 and Det M 6 ) b a M b i 2 a 4- a^o a + 2 |a b>0 | , + (a 0 ) 0 could be rewritten as 2 + a 6>2 bb b ) 0 ) + a 6 ) 3 M 2 6 ) 3 + a M a M , a _ or b x Appendix B. Inverse of the velocity matrix Det M _ i i ; DetM 123 = Det A = D e i A - | a ,i | a _i — a b M 2 6 = 0fi D Det A - (a , ) b 2 0 Thus, we obtain for A ^ = -(a ) , 2 M -4-1,0 = — A_i,i = Det A - a _i a A , 0 = Det A - (a ) , Ao,i = Ai,! = 0 O(, a )0 M5 b 6 ; 1 , 2 M — o ^ o , -(a _i) . 2 6 ) We immediately notice that the above formulae can be described by a general expression like — o~k,- Det A — a -h A ,p k P • b Therefore, from (B.3), we obtain, for our specific example, the inverse A' 1 -1 Det A ' a>b,-k {A ~ > 1)b k A itoi' _ ( A-l\ { = ) h ' p ~ n d b,-k b,-p a r >~ dk a p ~ a Det A • Hence, the "non-trivial" generalization to a continuous index set yields where Appendix B. Inverse of the velocity matrix D 124 Db-k V ' 5(k+p) Db-k Db,-p V which is the result we sought to show. O n the other hand, the determinant of D (D) can easily be determined for an infinite countable index set. Indeed, let D bb J^b,-n D, • • hb b,-n 6,0 n b 1 o D„ Du D ,n 6,0 (B.7) o 1 o where is a (2n + l ) x (2n + l ) m a t r i x whose only non-zero o V ) entries are numbers one on the diagonal running from the top right to the b o t t o m left corners. Appendix B. Inverse of the velocity matrix D 125 One can easily show that DetDr, A , - n A,-n+l A4-1) Y, ( - l ) + ( i + 1 + 0 0 0 0 1 0 0 1 0 0 (B.8) A , t Det n ) D _i 0 0 D xy 1 0 ... 0 0 0 y bi y h M 0 0 0 > where the m a t r i x w i t h i n the summation possesses only one row w i t h a complete string of zeros following the first element, namely, the row that starts w i t h element Hence, V n = -D + b b £ {-V D ,i ( - V - (i+l+n) { b l + n ) D ,_ b { (-1), i=—n since by crossing out the row containing the element we are left w i t h a (2n) x (2n) m a t r i x of the same form as (B.5) and whose determinant is —1. Thus, we obtain, making good use of (D b -D + bb J2 ( - 1 ) 2 |A,i A* = D , b>i | 2 i=—n n -D bb + Y I D ,k | , (we changed the index from i to k) 2 b k=—n and taking n to infinity V = lim V„ = ~D bb + £ | D (B.9) bJl k——oo We suppose that for a continuous index set, we have the obvious generalization +00 / -00 dk I D btk (B.10) Appendix B. Inverse of the velocity matrix D 126 Now, we would like to prove expression (3.39) for the determinant T>. A c c o r d i n g to (B.10), (3.36), and (3.38), we have V = L + Io where +O0 ( 1 r+oo (B.ll) / dx U' (x) + / n -oo +O0 / l f+OO ( dk dk q (t) f' (x)\ k dpq (t) U ) (B.12) -\ / p \J—oo dx f {x) f {x)\ p k J—oo J +OO 1 f+OO dvq* (t) / dy f;(y) v -oo r (y)\ k ) . J—oo The term I\ can be expanded as I = II + Ib ith x , k J—oo f+oo / -oo (B.13) w +O0 f+OO / dx (<l>' (x)) -oo f d We recall that (j) (x) = f (x) w i t h {f (x), b b +oo / f (x)} c f+OO r+oo 2 n -oo = (B.14) 1 dp q (t) f (x)j p p (B.15) forming an orthonormal basis where k dx 4> {x) = k k f'^jy^ {/_«, dk q (t) x dk q (t)fax)f' (x), J—oo 1 f+OO k n dx J—oo f+OO = f+OO 2/ 2 n 7 l 2 J—oo M 0 dx f (x) 2 b . (B.16) Therefore, one can express +00 / dk q (t) f' {x) k = c (t)f (x) k b b + / -oo f+OO drcr(t)f (x) r J—oo (B.17) U s i n g decomposition (B.17), we have for l\ l\ = - M c (t) + J J drds Cr(t)c (t) = - c {t)M + / 2 0 b 2 0 s J—oo dr\c {t) r | 5(r + s) (B.18) Appendix B. Inverse of the velocity matrix D 127 and, from 7 , 2 r+oo +oo / r+oo r+oo dp q (t) p r+oo r+oo / dx f'(x) f (x) = k dx f (x) k J — oo -oo r+oo r+oo / dp q (t) f'Jx) p J — oo J—oo +0O dr c (t) 5(r + k) / -oo c- (t), r (B.19) k and, similarly, r+oo +oo / r+oo r+oo dvq* (t) dy f';(y)r (y) v -oo r+oo r+oo = k / J—oo r+oo r+oo dyft{y) J—oo dv q* (t) f' *(y) v v J—oo +O0 / -oo drc* {t) S(r + k) d (t). r (B.20) k Using (B.19) and (B.20), 7 takes the simple form 2 r+oo +00 / dk c- {t) c*_ {t) k k -oo +oo dk\c (t)\ . (B.21) 2 / k -oo Therefore, adding l\ and I i n their simplified expressions (B.18) and (B.21), we obtain 2 r+oo 7 i + h = -j^ = -c {t)M ( k = 0 -c (t)M . 2 b 0 - / J—oo dp q (t) f (x)j k r+oo n ) r+OO dx { j f ^ dk q (t) f' (x)^j_^ 2 b -\ ( r+OO p p + I 2 r+oo dr\c (t)\ 2 r + / J—oo dr\cr(t)\ 2 (B.22) Appendix B. Inverse of the velocity matrix D 128 But ++ 0000 / / dx (f>' (x) n -00 f+00 / dk q (t) f' (x) k k \J—oo +00 / / dx f (x) ct(t) fb(x) b -00 + dr c (t) f (x) r \ r+00 dx f (x) -00 -,2 dr c (t) / dx f ( r J—oo [Cb(t) J r+oo + / 2 b b J—oo Mf 0 ct(t)M , 2 (B.23) so that equation (B.22) can be expressed as -cl(t) M I\ + h 0 1 (<*(*) M ) 2 1 / f+00 dx <j>' {x) ^ dk q (t) r r+oo = [j^ n f' (x) k k Hence, the determinant T> becomes V + I\ A° + h r+oo +OO / f+OO dx (<t>' (x)) - 2/ 2 n -00 I r f+00 — JM dx dk q (t) k \ 1 [ y ^ dx <p' (x) ( y ^ dk (t) n + 2 J ~dx + 0 qk n M~ L/-00 ^ = -{M + 2a + — a } 2 0 k n 0 ' 2 r (x)j </>' (x) (j^dk 0 f' (x) <f>' (x) J—oo { f+00 + V f+OO J—oo r J—OO +00 / N r+oo n ( X ) where k q (t) k V-00 ^ q k { t ) f k { X ) ) Appendix B. Inverse of the velocity matrix A a D f °° 129 ( f °° + = J \ + dx <f>'Jx) ( d k q (t) k Therefore, v = "k^ ' = which is exactly equation (3.39). + (« + ° 2M a M )\ c ° + M } f' (x)\ k Appendix C L i f e t i m e of a n environment-coupled h a r m o n i c oscillator In this appendix, we show that the H a m i l t o n i a n , derived from the Lagrangian (4.51) in the l i m i t of low amplitudes for the domain wall, makes the lifetime of a domain wall's energy level E take the form n nJM = (CI) where J ( ^ ) is the spectral density function and ui the natural frequency of the oscillating 0 domain wall parametrized by the coordinate X(t). For low amplitudes of the domain wall, V(X) in (4.50) can be approximated by a harmonic potential so that V (X) eff w i t h UJ 2 « V(X) = co X 2 2 0 - -2= £ VZ7T K reflecting the curvature at the bottom of the well ^X, (C.2) U) K V(X). A s explained in the footnote of page 99 in chapter 4, expression (C.2) can be thought of as a displaced harmonic potential well centered around a value of X different from zero owing to the presence of the constant force term E ipk/^k)Let us denote k the new center by X = X' (cf. fig. ( C . l ) ) . T h e energy eigenstates for the H a m i l t o n i a n Hx = ^ X \ t ) (a /-\/27r) 2 + V{X) (C.3) are the same as those for the ordinary harmonic oscillator except for their eigenvalues being all offset by the same constant term V(X') [Yi^k/^k)} - = -(a /Snul) 2 4 k In the l i m i t of low amplitudes for a domain wall, we have the total H a m i l t o n i a n H = H DW + \ £ {| q (t) k k | 2 + u 2 | q {t) | } z 130 2 k + X(t) £ k q (t) C , k k (C.4) Appendix C. Lifetime of an environment-coupled harmonic oscillator 131 Figure C . l : A shifted harmonic potential in the presence of a constant force. This figure shows how a harmonic potential is shifted by the presence of a constant force f . The new minimum for c the shifted potential V(X) where H = Hx DW = wx - fX 2 is now given by X' = f /u>. c c w i t h the subscript DW reminding us that the classical object is a domain wall. This H a m i l t o n i a n is composed of the usual harmonic oscillator term for the domain wall, HDW] a (infinite) series of harmonic oscillator H a m i l t o n i a n s 1 representing the environ- ment; and a coupling H a m i l t o n i a n for the last term. We suppose that H 0 = H + DW W 1 {\ q (t) k | + 2 LOI | q (t) k | } 2 (C.5) k represents the unperturbed H a m i l t o n i a n , whose energy eigenstates are denoted by I I)DW I n) b where +00 I n) = b (g) I n) k fe=—00 = I ri0) 0 I ni) ® I n_i) 0 ••• , T h e fact that we have complex coordinates q (t) does not invalidate our procedure, although this may appear odd. These coordinates are not totally independent of each other so as to ensure that we deal after all with a real field. This is analogous to going from a real basis {cos(/ca:),sin(fca;)} with real coefficients to a complex basis {e^ ] with complex coefficients. : k x Appendix C. Lifetime of an environment-coupled harmonic oscillator 132 and = x Hi Y k q (t)c k k represents the perturbation H a m i l t o n i a n . A c c o r d i n g to the standard perturbation theory, the second order correction to the energy level E , is given by n {nk} A£ , A 2 P n { n f c } i( V" = | m) —— {n \ H \DW N ^ — b x | m) DW | 2 b (C.6) m,- where {n } represents the quantum number for each oscillator q (t)" of the bath. u k k Now, (n \ D W {n \ Hi | m) b DW | m) = b (n \ D W = ((n \ D W {n \ | X Y ^X b Qk Ckj I rn) | m) DW b Y k {n |t q \ rn) m) ) C DW k b k . = {(n \ D W X Y fc | m) ) DW k {n C k | qk \ m ) 8 ,m , k nk k where w i t h a hat over a particular Kronecker symbol signifying that this symbol is to be omitted from the expression in which it appears. For the harmonic oscillator eigenstates, we know that ( n\ (n k | q k | m) k = j I— 1 / 2 {(n k \ a k \ m) + {n k k | a + | m )} k using the usual lowering and raising o p e r a t o r s , ^ a n d a , + for the = A; oscillator, th (v^fe ( k n \2u7 ) k \m -l) ^ v ^ f c ^ ^ - i k + Vm + Vm k k + 1 (n \m k + 1 6 , n k m k + i.} k + 1)} Appendix C. Lifetime of an environment-coupled harmonic oscillator 133 Hence, we obtain | (n \ D W < (n \ (n \ #1 | m) b D W | m) DW | = 2 b X | }n <m ( S k x j«n \ D X W (^"J I m) y Dw = 52 Cfc , 77 + ym (m p 1 ( v + l)£„ f e P m m * 5 { A ^ ^ - i I (n ' * » 5 n p low , m p n p m p + 1 }c) n p ) m t > | 2 + \/(m + \Jrn {m p + ^/mp + l < 5 , I ^BW I x ^ - i k k + l)(m„ + l)5„ , fc f c + l ) £ n , m + i 5„ , p p p fc mA: m f c + i 5 ,, n j m p + i _i j Thus, we have for the summation over the indices { m } , r ^ I (n | D V V (n | H t x (^n,{n } Im) | m) \ \ (n \ 2 D V V b = -E'm,{m }) r mi,ro m ,... [(-^n I (n \ Im) r 2l = D 3 D W X W — W ^ X \ m) x | m) DW | 2 b ^ m ) + Z)r ( ^ n ~ -Em )] | D V V E D b r m I m) r x 2 E) mi,m2,m3,... [(-^n I (n \ (n \ H + 5TJ (-En r -^m )] r r I x (C.7) 2 DW ^ k,p({ r}) a k,p m i , r o 2 , m . . . [ C ^ n m -^m) + zDr (-^n 3l -^m )] r r where Gfc, ({™r}) P = C f c J y ^ m * * n , m - l <^n ,m -l ^ , m C; ^ t + \/(m \jm [m v k + l)5 n p i m p _! 5 n f c ! fc m + l <W,m ^ n , m f c fc p t p p f c 4 ,m p (C.8) p + l)(mp + 1) <y ,m* + l ^n ,m +l 5 „ , nfc p + y/m (m k p + 1) 5 p n p > m p + 1 p fc mjl c5 , np + mp ^ . m f c - l 4*,™* 4 , m | p p Appendix C. Lifetime of an environment-coupled harmonic For the " r " oscillator, the Kronecker symbol <5„ t h r:TOr oscillator 134 is present for each product 8 8 n k t t n k n p t m t unless k = p = r , in which case we have — (o~n ,m ) $n ,m 5n ,m k k k k 0 ($n_i,m_i) 0 (^m,mi) ' ' " (^n ,m ) ' ' " (4i 3,mi 3) 2 k k 2 W h a t this means is that, when we perform our summation 52 m • 2 52 = > o y m t n e mi,m2,ms,... r terms for which k = p i n (C.8) may possibly survive since all other terms for which k / p vanish because would require n = m for a l l r in order to be non-zero, but in 8 m 5n ,m nk> k p r p such a case the Kronecker symbols 5 , n t and 5 n m p _ i 5rn;,m +i i fc na m t r _i 6 _i, ^ n p > m p <5 r l f e > m f c + 1 7ip;7Tlp+ i, <5 , np o^^-i, mp+1 fc,p({ v}) would all vanish causing afc ({m }) itself to have the rj !P r same fate; and i f any of those latter Kronecker symbols were non-zero, then it would mean that for some r , n ^ m making 8 r r n k t m k 8 n and, therefore, a ({m }) m kjP both zero. r Hence, we obtain 1 / (n \ m y D (n \ Hi | m) / -p rp W b \&n,{n ) r r ~ | m) \ _ 2 DW \ b ~ \ \ N \DW ^V ^m,{m }) . . I m) 2 D W | X r ^ ^ ^ Jfc,p mi,ra ,ms,... [C^n E) 2 = | (n \ D W ^ I 12 Ck | 2 u ) k 2 ft ( 52 1™ *( ^ , m , - l ) 2 | {n k | 2 r k k X E )] 1lr mr x ktk 2 \DW r a ({m }) mi,m ,m3,... [(-^n 5n ,m 0~n ,m k D V V ({m }) + J2 (E j2 k = X | m) m ) P |m ) D W + ^ E) m r (i? ^m )] n r r | X + (m + 1) (8 ,m k r nk k + lf 3n ,m 0~n ,k k [(£„ - E ) + 5J (£„ - £ m r m k k TOr , Appendix C. Lifetime of an environment-coupled = ^ |C | h | 2u; n since the term containing 5 m -i o~n ,m +i llk>mk k k nk | - E i) nk n nk m 1lk - E ^) nk x J, always vanishes no matter what o~n ,m 5 ,m k | (E - E ) + (E nk+ k 135 | x DW k \ (E - E ) + (E fe oscillator X | m) D W n +1 2 k k | (n \ harmonic k m i , 77Z2,777.3, • • • are. Therefore, making use of the fact that E = {n + 1/2) huo nk k k for each oscillator k, expression (C.6) yields - ? M " m ~~2^r~ ? 1 w - (»«,)' I, ' <c 9) w here Mim = I AE = E = ( IDW ^ I n — E n ft {iO ft w n n m - m) | , and 2 DW m uo ) m . T h e second order correction to the energy level where a l l oscillators of the bath would be at their ground states (n* = 0 for all &), were it not for the perturbation, is given by 2&E {n =0} 7h k = E M «™ E 00 fc=—00 77-7^ I r i 12 2wji T77T 1 (w n m — ojk) Appendix C. Lifetime of an environment-coupled harmonic oscillator 136 However, by definition, we have [13, 45] z U) k k for a bath of unit-mass oscillators. Substituting ( C . l l ) into (C.10) yields ^ M = 2^ / m Jo nm 2AE n{n 0 } K foo dv J{y) (iO - nm . (C.12) v) ' In order to compute this integral, we add an infinitesimal imaginary part to the denominator and split the integral as follows 2A£ ,{ n B t = 0 } = £ ~ ~ ~ m = K since M n n o} = = | (n \ l 0 / + d 0 J lim / m<n nA;= n u , (W n m 7T _ lim fO° It) - (U) , dv — e->-0+ Jo {y [V ~M • \ + • \ V + —— {V / v , - n m dv e-M)+ Jo IT M n) i ^ > m 2AS { l e > m^n X | n) — m - - Km lim / T It)) J { » ) r°° nm e^0+ 7i + h , D W n Jo it)) aV - [V J(y) — — - (0J nm - —r, le)) (C.13) | = 0. 2 D V V In the classical regime, the macroscopic coordinate X{t) obeys a general phenomenological equation of motion of the form where K is a linear operator subject to the requirement of causality [45]. F o r our case, we have V(X) = V(X) i n expression (C.3). T h e Fourier transform of K, K(v), is analytical i n the lower half of the complex plane, K(0) = 0, and we have the following 2 2 A s a consequence of causality. Appendix C. Lifetime of an environment-coupled harmonic oscillator 137 relation [13, 45]: so that K{u) = i J(u) Im K{u) = J(v). for v € IR (C.14) Now, in order to solve , h = f°° M y E ^ where u) (~ nm) V- l 27T N i n l *->Q+ J{y) , / d u 7 7 ^TV ( kI- K J-OO r a - «€)) > 0, we transform u into a complex variable z and choose the integration nm path TR by closing the contour in the lower half-plane, where the pole oj nm — ie is located, and let R —>• oo (cf. fig. (C.2)) to obtain V- ( - M ) f / J{\z\) \ „ > -— l i m < l i m / ——;—^—-——•} = > £? 2TT *+o*\R->ooJr (\z\-(u> -ie))j ^ (-M n m n R ?m n n m 2TT ) . l i m 2inJ{\ uj —ie \) " U N nm M by the Residue theorem. We suppose that J ( | z |) is sufficiently well-behaved so that the integrand in the above integral vanishes on the semicircle Cn as R —> oo. Thus, ii = E ( n m ) i limi 7(| co - M nm m<n = h = e ^ (-M E M n m n m ^ 0 )i/(w AT'K m ) - it \) + B m ), (C.16) m<n according to (C.14). For the t e r m 72, we do not have to resort to a complex integral on account of the fact that for m > n , oj nm < 0. Indeed, we have Appendix C. Lifetime of an environment-coupled Figure C.2: T h e integration path T R harmonic oscillator 138 for the partial evaluation of the second order per- t u r b a t i o n energy term. The pole, denoted by x in the figure, is barely located in the lower half-plane on account of an infinitesimal imaginary part added to the harmonic oscillator's frequency difference ui . nm The path T R is made up of the line extending from - R to R on the real line and a semi-circle C R of radius R in the lower half-plane. For the proper evaluation of our expression of interest, we let R tend to infinity. m>n - 2s TV 6-S-0+ JO \V + I w n m I + L E ) m>n m>n •n.^ E n ("Mnm) r TT JO J d v HH=) r 4, 7T Jo ^ (" + I Wnm |) I-^M I w n T O (C . 17) | 4-ie) using expression (C.14) once again. O f the terms I\ and 7 in (C.13), only Ii is complex-valued. According to the standard 2 interpretation [39], one relates the imaginary part of 2&E ^ =Q} w i t h the lifetime of the n energy level E n nk of the particle (parametrized by the coordinate X(t)) for which, were it not for the perturbation, all the bath's oscillators would normally be in their ground states; a situation which is likely to be realized at very low temperatures and consistent Appendix C. Lifetime of an environment-coupled harmonic oscillator 139 w i t h our requirement that each oscillator be weakly perturbed. T h e relation is as follows 3 Im( A£ , 2 where r n n { n f c } ) = Im(/!) = - \ - is the lifetime of energy level E (CAS) for the system's harmonic oscillator. n We showed previously that Ji with M nm = | (n \ D W X \ m) M nm 2 = (m(5 , _i) nm K*(to ) nm | . Now, DW M Y = n 2 m + ( m + 1) (5 )} 2 n>m+1 for a harmonic oscillator so that h = h nm K*{oj ) M n<m = n A"*(w = -inJ(vi) n>n _i), (C.19) using expression (C.14). Substituting equation (C.19) into (C.18), we obtain K \ n J{OJ ,n-l) n n J(hco ) 0 1 1 = ~ — = ^ 2 T , (C.20) N In reality, Im(2 AE„,{n }) relates to the lifetime of the energy level E' which corresponds to the unperturbed energy level E ^ y, with the harmonic oscillator X(t) having quantum number n and the bath's oscillators, quantum numbers n , in the sense that, by slowly turning off the coupling between the particle and the bath, E' ^ y would approach the value of E ^t y. We implicitly suppose that the interaction does not create additional quantum states for a complete description of the universe (system plus bath). Obviously, this is always a matter of semantics when we assign the full thrust of the perturbation effect to the system's particle alone, as we did above. 3 fc n n nk k n Hk n nk Appendix C. Lifetime of an environment-coupled harmonic oscillator 140 since, for a harmonic oscillator, a l l energy levels are separated by the same difference hto 0 where to is the oscillator's natural frequency. Using the convention h~ = 1, we obtain 0 expression (5.1) n J(io ) 0 = - —, which is the result we wished to demonstrate in this appendix.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- One-dimensional soliton dynamics in the presence of...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
One-dimensional soliton dynamics in the presence of a pinning potential Leduc, Benoit 1995
pdf
Page Metadata
Item Metadata
Title | One-dimensional soliton dynamics in the presence of a pinning potential |
Creator |
Leduc, Benoit |
Date Issued | 1995 |
Description | The dissipative effect of (what could be modeled as) a pinning center on the lowamplitude motion of a quasi-one-dimensional ferromagnetic soliton (domain wall) at temperature T = 0 is investigated. The method of collective coordinates is used to eliminate the problem of the zero-frequency eigenmode and the question is, subsequently, cast into a Caldeira-Leggett form to calculate the parametrized spectral density function. It is shown that this function equals zero below a specific frequency with the consequence that magnons, in the presence of a domain wall, have no dissipative effect at T = 0 on a slow moving wall trapped by a soft pinning center. |
Extent | 5528110 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Date Available | 2009-02-06 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0086976 |
URI | http://hdl.handle.net/2429/4200 |
Degree |
Master of Science - MSc |
Program |
Physics |
Affiliation |
Science, Faculty of Physics and Astronomy, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 1996-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-ubc_1996-0064.pdf [ 5.27MB ]
- Metadata
- JSON: 831-1.0086976.json
- JSON-LD: 831-1.0086976-ld.json
- RDF/XML (Pretty): 831-1.0086976-rdf.xml
- RDF/JSON: 831-1.0086976-rdf.json
- Turtle: 831-1.0086976-turtle.txt
- N-Triples: 831-1.0086976-rdf-ntriples.txt
- Original Record: 831-1.0086976-source.json
- Full Text
- 831-1.0086976-fulltext.txt
- Citation
- 831-1.0086976.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0086976/manifest