DOCSLIB.ORG

Soler Model: Stability, Bi-Frequency Solitons, and SU(1,1) Andrew Comech, Texas A&M University and I.I.T.P

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Soler Model: Stability, Bi-Frequency Solitons, and SU(1,1) Andrew Comech, Texas A&M University and I.I.T.P Soler model: stability, bi-frequency solitons, and SU(1,1) Andrew Comech, Texas A&M University and I.I.T.P. (Moscow) Nabile Boussa¨ıd, Universite´ Bourgogne – Franche-Comte,´ Besanc¸on, France July 2016 Einstein: E2 = p2 + m2 Schrodinger:¨ (i∂ )2ψ = ( i )2ψ + m2ψ t − ∇ 2 E = m2 + p2 m + p ≈ 2m p [Schrodinger¨ 26]: i∂ ψ = 1 ( i )2ψ t 2m − ∇ [Dirac28]: E = p2 + m2 = α p + βm, · p i∂ ψ = ( iα + βm) ψ, ψ(x,t) C4, x R3 t − ·∇ ∈ ∈ Dm | {z } 2 2 Dirac matrices α (1≤j≤3) and β: self-adjoint, ( iα + βm) = ∆+ m j − ·∇ − 0 σj I2 0 Standard choice: αj = (Pauli matrices), β = σj 0 0 I2 − 1 Dirac–Maxwell system ψ¯ := ψ∗β γµ(i∂ eA )ψ = mψ, ✷Aµ = eψγ¯ µψ, ∂ Aµ = 0. µ − µ µ [Esteban et al.96, Abenda98]: solitary waves ψ(x,t) = φ(x)e−iωt, ω ( m,m) ∈ − Dirac–Maxwell and Dirac–Coulomb systems 40 30 −iωt Q(φωe ) 20 10 0 E(φ e−iωt) ω -10 -20 -1 -0.5 0 0.5 1 m<ω m − ≤ [Comech15]: solitary waves with ω = m for Dirac–Coulomb; φ(x) e−const√|x|. ∼ γ0 = β, γj = βαj , 1 j 3 ≤ ≤ 2 Self-interacting spinors ψ¯ := ψ∗β Models of self-interacting spinor field: [Ivanenko38, Finkelstein et al.51, Finkelstein et al.56, Heisenberg57] [...] Soler model [Ivanenko38, Soler70], scalar self-interaction: µ k+1 LSoler = ψ¯(iγ ∂ m)ψ + ψψ¯ k > 0 µ − | | in (1+1)D: massive Gross-Neveu model [Gross & Neveu74] Massive Thirring model [Thirring58], vector self-interaction: k+1 µ µ 2 LMTM = ψ¯(iγ ∂ m)ψ + ψγ¯ ψ ψγ¯ ψ k > 0 µ − µ J µ(x) = ψ¯(x)γµψ(x): “charge-current density” 3 Soler model: NLD with scalar self-interaction i∂ ψ = ( iα + mβ) ψ (ψψ¯ )kβψ, ψ(x,t) CN , x Rn t − ·∇ − ∈ ∈ Dm | {z } [Soler70, Cazenave & Vazquez´ 86]: existence of solitary waves in R3, • ψ(x,t) = φ (x)e−iωt, ω (0,m), φ H1(R3) ω ∈ ω ∈ Attempts at stability: [Bogolubsky79, Alvarez & Soler86, Strauss & Vazquez´ 86] ... • Numerics suggest that (all?) solitary waves in 1D cubic Soler model are stable: • [Alvarez & Carreras81, Alvarez & Soler83, Cooper et al.10, Berkolaiko & Comech12], [Mertens et al.12, Shao et al.14, Cuevas-Maraver et al.14] Assuming linear stability, one tries to prove asymptotic stability • [Pelinovsky & Stefanov12][Boussaid & Cuccagna12][Comech, Phan, Stefanov14] 4 Nonrelativistic limit of NLD: ω . m [Ounaies00, Guan08] v(x) Solitary wave: ψ(x,t) = e−iωt; v, u C2 x R3 u(x) ∈ ∈ 0 σ 1 0 iψ˙ = i ·∇ + m (ψψ¯ )k ψ, − σ 0 − 0 1 n ·∇ − o v u v ω iσ + m v 2k u ≈ − ·∇ v −| | u − If ω . m, u v 1: 2mu iσ v, v satisfies NLS: | |≪| | ≪ ≈ − ·∇ 1 (m ω)v = ∆v v 2kv − − −2m −| | Scaling: v(x) = ǫ1/kΦ(ǫx), ǫ = √m ω, − 1 Φ = ∆Φ Φ 2kΦ, Φ(x) R, x Rn. − −2m −| | ∈ ∈ 5 NLD: linearization at a solitary wave Given φ (x)e−iωt, ω ( m,m), consider ψ(x,t) = φ (x)+ r(x,t) e−iωt ω ∈ − ω Linearized eqn on r(x,t), i∂ r = D r ωr + ... t m − (m + ω)i Re r 0 Dm ω + ... Re r ∂t = − Im r Dm + ω + ... 0 Im r − Aω | {z } (m − ω)i ∂tR = AωR; σ(Aω) iR (“spectral stability”) ??? ⊂ ✲ σess(Aω) σ(D ω) m − r0 m ω m ω − − − 6 σ(Aω) Linear instability of NLD 2mi (m + ω)i i∂ ψ = D ψ (ψψ¯ )kβψ, x Rn t m − ∈ −iωt ψ(x,t) = φω(x)e (m − ω)i Theorem 1 ([Comech, Guan, Gustafson12]). t tt ✲ ✲ −λω λω If k > 2/n, so that NLSk is linearly unstable, then for ω . m, λ σ (A ), ∃ ± ω ∈ d ω Re λω > 0, λω 0 ω−→→m Unstable for n = 1, k > 2 (above quintic) −2mi Unstable for n = 2, k > 1 (above cubic) ω < m ω m Unstable for n = 3, k > 2/3 = Proof: Rescale; use Schur reduction and Rayleigh-Schrodinger¨ perturbation theory. ✷ 7 σ(Aω) Bifurcations from σess (m + ω0)i λω Let 0 ω0 m λ 0 ω r r r ✈r r ≤ ≤ r r r r Theorem 2 ([Boussaid & Comech16]). (m − ω0)i Assume: λω σp(Aω), Re λω = 0, ∈ 6 t ✲ R λω λω0 i ω−→→ω0 ∈ Then: r r r r r r ✈r r r λ m + ω0, | ω0 | ≤ λ σ (A ) 0 i(m + ω0) ω0 ∈ p ω0 ∪ { } ∪ { } LAP methods of [Agmon75], [Berthier & Georgescu87]: u 1 c (D z)u 2 , z C [ m,m], s> 1/2, z m > δ k kH−s ≤ z,s,δk m − kLs ∈ \ − | ± | s u 1 := (1 + x ) u) 1 k kHs k | | kH 8 σ(Aω) Linear stability of NLD 2mi (m + ω)i 2ωi t i∂ ψ = D ψ (ψψ¯ )kβψ, x Rn t m − ∈ −iωt ψ(x,t) = φω(x)e (m − ω)i t ✲ ✲ t Theorem 3 ([Comech11]). If n 3, then λ = 2ωi σ (A ). ≤ ± ∈ p ω −2ωi t Theorem 4 ([Boussaid & Comech16]). n 3, k . 2/n: −2mi ≤ −iωt φωe are linearly stable for ω . m ω < m ω = m Including the “charge-critical” case k = 2/n 9 Eigenvalues λ = 2ωi in the Soler model ± λ = 2ωi is bad for proving asymptotic stability! ± [Boussaid & Cuccagna12], [Comech, Phan, Stefanov14] Asymptotic stability of solitons in Soler model for “radially symmetric” case: 1. perturbations orthogonal to translations; 2. perturbations orthogonal to 2ωi eigenvectors ± 10 SU(1,1) invariance of the Soler model and the Soler charge Theorem 5 ([Boussaid & Comech16]). 1. Soler model hamiltonian is invariant under transformations ψ(x,t) (a ibγ2C)ψ(x,t), a, b C, a 2 b 2 = 1 −→ − ∈ | | −| | a b 2. These transformations form the group SU(1, 1) U(1) ⊃ ¯b a¯ 3. Conserved charges: Q = ψ∗ψdx > 0 Λ = R ψt( iγ2)ψ dx C − ∈ R 4. Bi-frequency solutions: if φ(x)e−iωt is a solitary wave, so is ψ(x,t) = ae−iωtφ(x) beiωtiγ2C(φ(x)) − Note: if b 1, then ψ(x,t) e−iωt φ(x) be2iωtiγ2C(φ(x)) | | ≪ ≈ − 11 Q(ω), 10 10 8 8 E(ω) 6 6 4 4 2 2 0 0 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 1 1 0.8 0.8 0.6 0.6 0.4 0.4 λ iR+ ∈ 0.2 0.2 0 0 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 Figure 1: (1+1)D Soler model. LEFT: k = 2 (“charge-critical”); RIGHT: k = 3. TOP ROW: charge and energy of the solitary waves as functions of ω ∈ (0, 1). BOTTOM ROW: Spectrum on the upper half of the imaginary axis. Note the exact eigenvalue λ = 2ωi. 12 1 0.5 ) λ 0 Im( −0.5 −1 0 0.5 1 ω Figure 2: (2+1)D cubic “charge-critical” Soler model [Cuevas-Maraver et al.16] 13 1 0.5 ) λ 0 Im( −0.5 −1 0 0.5 1 ω Figure 3: (3+1)D cubic “charge-supercritical” Soler model [Cuevas-Maraver et al.16] 14 Theorem 6 ([Berkolaiko, Comech, Sukhtayev15]). Both Q′(ω) = 0 and E(ω) = 0 indicate collision of eigenvalues at λ = 0 15 −iωt 10 Q(φωe ) 5 −iωt E(φωe ) 0 -5 -1 -0.5 0 0.5 1 1 0.8 0.6 λ iR+ 0.4 ∈ 0.2 0 -1 -0.5 0 0.5 1 m<ω<m − Figure 4: “Quadratic” massive Thirring model, k = 1/2. TOP: energy and charge as functions of ω ∈ (−m, m). BOTTOM: The spectrum of the linearization at a solitary wave. Dotted eigenvalue collides with its opposite at the origin when ω⋆ ≈ −0.6276m, where E(ω⋆) = 0 15 References [Abenda98] S. Abenda, Solitary waves for Maxwell-Dirac and Coulomb-Dirac models, Ann. Inst. H. Poincare´ Phys. Theor.,´ 68 (1998), pp. 229–244. [Agmon75] S. Agmon, Spectral properties of Schrodinger¨ operators and scattering theory, Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4), 2 (1975), pp. 151–218. [Alvarez & Carreras81] A. Alvarez & B. Carreras, Interaction dynamics for the solitary waves of a nonlinear Dirac model, Phys. Lett. A, 86 (1981), pp. 327–332. [Alvarez & Soler83] A. Alvarez & M. Soler, Energetic stability criterion for a nonlinear spinorial model, Phys. Rev. Lett., 50 (1983), pp. 1230–1233. [Alvarez & Soler86] A. Alvarez & M. Soler, Stability of the minimum solitary wave of a nonlinear spinorial model, Phys. Rev. D, 34 (1986), pp. 644–645. [Berkolaiko & Comech12] G. Berkolaiko & A. Comech, On spectral stability of solitary waves of nonlinear Dirac equation in 1D, Math. Model. Nat. Phenom., 7 (2012), pp. 13–31. 16 6 Q(ω) 5 4 3 2 E(ω) 1 0 1 -1 -0.5 0 0.5 1 0.8 0.6 0.4 λ iR+ ∈ 0.2 0 -1 -0.5 0 0.5 1 m<ω<m − Figure 5: Massive Thirring model with k = 1. TOP: energy and charge as functions of ω ∈ (−1, 1). BOTTOM: The spectrum of the linearization at a solitary wave on the upper half of the imaginary axis. No nonzero eigenvalues in completely integrable model. 17 5 Q(ω) 4 3 2 1 E(ω) 0 -11 -1 -0.5 0 0.5 1 0.8 0.6 0.4 λ iR+ ∈ 0.2 0 -1 -0.5 0 0.5 1 Figure 6: Massive Thirring model with k = 2, “charge critical”. TOP: energy and charge as functions of ω ∈ (−1, 1).
Load more

Recommended publications
  • What Is the Dirac Equation?
    What is the Dirac equation? M. Burak Erdo˘gan ∗ William R. Green y Ebru Toprak z July 15, 2021 at all times, and hence the model needs to be first or- der in time, [Tha92]. In addition, it should conserve In the early part of the 20th century huge advances the L2 norm of solutions. Dirac combined the quan- were made in theoretical physics that have led to tum mechanical notions of energy and momentum vast mathematical developments and exciting open operators E = i~@t, p = −i~rx with the relativis- 2 2 2 problems. Einstein's development of relativistic the- tic Pythagorean energy relation E = (cp) + (E0) 2 ory in the first decade was followed by Schr¨odinger's where E0 = mc is the rest energy. quantum mechanical theory in 1925. Einstein's the- Inserting the energy and momentum operators into ory could be used to describe bodies moving at great the energy relation leads to a Klein{Gordon equation speeds, while Schr¨odinger'stheory described the evo- 2 2 2 4 lution of very small particles. Both models break −~ tt = (−~ ∆x + m c ) : down when attempting to describe the evolution of The Klein{Gordon equation is second order, and does small particles moving at great speeds. In 1927, Paul not have an L2-conservation law. To remedy these Dirac sought to reconcile these theories and intro- shortcomings, Dirac sought to develop an operator1 duced the Dirac equation to describe relativitistic quantum mechanics. 2 Dm = −ic~α1@x1 − ic~α2@x2 − ic~α3@x3 + mc β Dirac's formulation of a hyperbolic system of par- tial differential equations has provided fundamental which could formally act as a square root of the Klein- 2 2 2 models and insights in a variety of fields from parti- Gordon operator, that is, satisfy Dm = −c ~ ∆ + 2 4 cle physics and quantum field theory to more recent m c .
    [Show full text]
  • Transverse Instability of Line Solitary Waves in Massive Dirac Equations
    J Nonlinear Sci (2016) 26:365–403 DOI 10.1007/s00332-015-9278-1 Transverse Instability of Line Solitary Waves in Massive Dirac Equations Dmitry Pelinovsky1 · Yusuke Shimabukuro1 Received: 7 November 2014 / Accepted: 23 October 2015 / Published online: 6 November 2015 © Springer Science+Business Media New York 2015 Abstract Working in the context of localized modes in periodic potentials, we con- sider two systems of the massive Dirac equations in two spatial dimensions. The first system, a generalized massive Thirring model, is derived for the periodic stripe poten- tials. The second one, a generalized massive Gross–Neveu equation, is derived for the hexagonal potentials. In both cases, we prove analytically that the line solitary waves are spectrally unstable with respect to periodic transverse perturbations of large periods. The spectral instability is induced by the spatial translation for the general- ized massive Thirring model and by the gauge rotation for the generalized massive Gross–Neveu model. We also observe numerically that the spectral instability holds for the transverse perturbations of any period in the generalized massive Thirring model and exhibits a finite threshold on the period of the transverse perturbations in the generalized massive Gross–Neveu model. Keywords Massive Thirring model · Massive Gross-Neveu equation · Line solitary waves · Transverse instability · Lyapunov-Schmidt reduction method · Chebyshev interpolation method Mathematics Subject Classification 35Q41 · 37K40 · 37K45 Communicated by Michael I. Weinstein. B Dmitry Pelinovsky [email protected] 1 Department of Mathematics and Statistics, McMaster University, Hamilton, ON L8S 4K1, Canada 123 366 J Nonlinear Sci (2016) 26:365–403 1 Introduction Starting with pioneer contributions of V.E.
    [Show full text]
  • AN OVERVIEW on LINEAR and NONLINEAR DIRAC EQUATIONS 1. Introduction and Basic Properties of the Dirac Operator. in Relativistic
    DISCRETE AND CONTINUOUS Website: http://AIMsciences.org DYNAMICAL SYSTEMS Volume 8, Number 2, April 2002 pp. 381–397 AN OVERVIEW ON LINEAR AND NONLINEAR DIRAC EQUATIONS MARIA J. ESTEBAN AND ERIC SER´ E´ Abstract. In this paper we review the variational treatment of linear and nonlinear eigenvalue problems involving the Dirac operator. These problems arise when searching for bound states and bound energies of electrons sub- mitted to external or self-consistent interactions in a relativistic framework. The corresponding energy functionals are totally indefinite and that creates difficulties to define variational principles related to those equations. Here we describe recent works dealing with this kind of problems. In particular we describe the solutions found to linear Dirac equations with external potential, the Maxwell-Dirac equations, some nonlinear Dirac equations with local non- linear terms and the Dirac-Fock equations arising in the description of atoms and molecules. 1. Introduction and basic properties of the dirac operator. In relativistic quantum mechanics [11], the state of a free electron is represented by a wave func- tion Ψ(t, x) with Ψ(t, .) ∈ L2(IR3, CI 4) for any t. This wave satisfies the free Dirac equation: 3 X 2 i ∂t Ψ = H0Ψ, with H0 = −ic~ αk∂k + mc β, (1.1) k=1 where c denotes the speed of light, m > 0, the mass of the electron, and ~ is Planck’s constant. In the Dirac equation, α1, α2, α3 and β are 4×4 complex matrices, whose standard form (in 2 × 2 blocks) is I 0 0 σk β = , αk = (k = 1, 2, 3), 0 −I σk 0 with 0 1 0 −i 1 0 σ = , σ = , σ = .
    [Show full text]
  • Tight-Binding Model and Ab Initio Calculation of Silicene with Strong
    Tight-binding model and ab initio calculation of silicene with strong spin-orbit coupling in low-energy limit Chen-Huan Wu ∗ Key Laboratory of Atomic & Molecular Physics and Functional Materials of Gansu Province, College of Physics and Electronic Engineering, Northwest Normal University, Lanzhou 730070, China April 6, 2018 We carry out both the tight-binding model and the ab initio to study the layered silicene, the spin, valley, sublattice degrees of freedom are taken into consider and the effects of electric field, magnetic field, and even the light in finite frequency together with its interesting optical propertice, which are all closely related to the spin-orbit coupling and Rashba coupling and lead to the tunable phase transitions (between the nontrivial topological phase which has nonzero Chern number or nonzero spin Chern number and the trivial phase). An exotic quantum anomalous Hall insulator phase are found which has nonzero spin Chern number and nonzero valley Chern number and as a giant-application-potential spintronic and valleytronics in the two- ternimate devices based on the monolayer silicene for the information transmission. In fact, the gap-out action can be understanded by analyse the Dirac mass as well as the Zeeman splitting or the external-field-induced symmetry-broken, and the changes of gap has a general nonmonotone-variation characteristic under both the effects of electron filed-induced Rashba-coupling R2(E ) and the electron field-induced band gap, and the band inversion related to the spin-orbit⊥ coupling which absorbs both the spin and orbital angular of momentum may happen during this process.
    [Show full text]
  • Variational Methods in Relativistic Quantum Mechanics
    VARIATIONAL METHODS IN RELATIVISTIC QUANTUM MECHANICS MARIA J. ESTEBAN, MATHIEU LEWIN, AND ERIC SER´ E´ Abstract. This review is devoted to the study of stationary solutions of lin- ear and nonlinear equations from relativistic quantum mechanics, involving the Dirac operator. The solutions are found as critical points of an energy func- tional. Contrary to the Laplacian appearing in the equations of nonrelativistic quantum mechanics, the Dirac operator has a negative continuous spectrum which is not bounded from below. This has two main consequences. First, the energy functional is strongly indefinite. Second, the Euler-Lagrange equations are linear or nonlinear eigenvalue problems with eigenvalues lying in a spectral gap (between the negative and positive continuous spectra). Moreover, since we work in the space domain R3, the Palais-Smale condition is not satisfied. For these reasons, the problems discussed in this review pose a challenge in the Calculus of Variations. The existence proofs involve sophisticated tools from nonlinear analysis and have required new variational methods which are now applied to other problems. In the first part, we consider the fixed eigenvalue problem for models of a free self-interacting relativistic particle. They allow to describe the localized state of a spin-1/2 particle (a fermion) which propagates without changing its shape. This includes the Soler models, and the Maxwell-Dirac or Klein- Gordon-Dirac equations. The second part is devoted to the presentation of min-max principles al- lowing to characterize and compute the eigenvalues of linear Dirac operators with an external potential, in the gap of their essential spectrum. Many con- sequences of these min-max characterizations are presented, among them a new kind of Hardy-like inequalities and a stable algorithm to compute the eigenvalues.
    [Show full text]
  • LAN-DISSERTATION-2019.Pdf (967.2Kb)
    ON THE SPECTRAL STABILITY OF SOLITARY WAVES A Dissertation by RUOMENG LAN Submitted to the Office of Graduate and Professional Studies of Texas A&M University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Chair of Committee, Andrew Comech Co-Chair of Committee, Alexei Poltoraski Committee Members, Artem Abanov Gregory Berkolaiko Head of Department, Sarah Witherspoon December 2019 Major Subject: Mathematics Copyright 2019 Ruomeng Lan ABSTRACT We study the spectral stability of the solitary wave solutions to the nonlinear Dirac equations. We focus on two types of nonlinearity: the Soler type and the Coulomb type. For the Soler model, we apply the Evans function technique to explore the point spectrum of the linearized operator at a solitary wave solution to the 2D and 3D cases. For the toy Coulomb model, the solitary wave solutions are no longer SU(1; 1) symmetric. We show numerically that there are no eigenvalues near 2!i in the nonrelativistic limit (! . m) and the spectral stability persists in spite of the absence of SU(1; 1) symmetry. ii DEDICATION To my wife Yao and my son Xiao. iii ACKNOWLEDGMENTS I would like to express my deepest gratitude to my advisor, Dr. Andrew Comech, for intro- ducing me to the nonlinear Dirac equation and the spectral theory. I appreciate his patience and encouragement. I cannot complete the dissertation without his guidance and help. I would like to thank the committee member, Dr. Alexei Poltoratski, Dr. Gregory Berkolaiko and Dr. Arterm Abanov for offering their time and support. I thank Dr. Ciprian Foias and Dr.
    [Show full text]
  • Stability and Well-Posedness in Integrable Evolution Equations
    STABILITY AND WELL-POSEDNESS IN INTEGRABLE NONLINEAR EVOLUTION EQUATIONS STABILITY AND WELL-POSEDNESS IN INTEGRABLE EVOLUTION EQUATIONS By YUSUKE SHIMABUKURO, B.Sc., M.Sc. A Thesis Submitted to the Department of Mathematics and Statistics and the School of Graduate Studies of McMaster University in Partial Fulfilment of the Requirement for the Degree of Doctor of Philosophy McMaster University c Copyright by Yusuke Shimabukuro DOCTOR OF PHILOSOPHY (2016) McMaster University Department of Mathematics and Statistics Hamilton, Ontario, Canada TITLE: Stability and well-posedness in integrable evolution equations AUTHOR: Yusuke Shimabukuro B.Sc. (University of Arizona) M.Sc. (University of Arizona) SUPERVISOR: Dr. Dmitry Pelinovsky NUMBER OF PAGES: IX, 153 II Abstract This dissertation is concerned with analysis of orbital stability of solitary waves and well-posedness of the Cauchy problem in the integrable evolution equations. The analysis is developed by using tools from integrable systems, such as higher- order conserved quantities, B¨acklund transformation, and inverse scattering trans- form. The main results are obtained for the massive Thirring model, which is an integrable nonlinear Dirac equation, and for the derivative NLS equation. Both equations are related with the same Kaup-Newell spectral problem. Our studies rely on the spectral properties of the Kaup-Newell spectral problem, which convey key information about solution behavior of the nonlinear evolution equations. III Acknowledgements First and foremost I would like to express my deepest gratitude to my advisor, Professor Dmitry Pelinovsky. It has been truly an exciting journey under his guidance. His vast knowledge, ideas, and insights have inspired and helped me to complete this dissertation.
    [Show full text]
  • Solitary Waves in the Nonlinear Dirac Equation
    Solitary waves in the Nonlinear Dirac Equation Jesus´ Cuevas-Maraver, Nabile Boussa¨ıd, Andrew Comech, Ruomeng Lan, Panayotis G. Kevrekidis, and Avadh Saxena Abstract In the present work, we consider the existence, stability, and dynamics of solitary waves in the nonlinear Dirac equation. We start by introducing the Soler model of self-interacting spinors, and discuss its localized waveforms in one, two, and three spatial dimensions and the equations they satisfy. We present the associ- ated explicit solutions in one dimension and numerically obtain their analogues in higher dimensions. The stability is subsequently discussed from a theoretical per- spective and then complemented with numerical computations. Finally, the dynam- ics of the solutions is explored and compared to its non-relativistic analogue, which is the nonlinear Schrodinger¨ equation. Jesus´ Cuevas-Maraver Grupo de F´ısica No Lineal, Universidad de Sevilla, Departamento de F´ısica Aplicada I, Escuela Politecnica´ Superior. C/ Virgen de Africa,´ 7, 41011-Sevilla, Spain, Instituto de Matematicas´ de la Universidad de Sevilla (IMUS). Edificio Celestino Mutis. Avda. Reina Mercedes s/n, 41012-Sevilla, Spain e-mail: [email protected] Nabile Boussa¨ıd Universite´ de Franche-Comte,´ 25030 Besanc¸on CEDEX, France Andrew Comech St. Petersburg State University, St. Petersburg 199178, Russia Department of Mathematics, Texas A&M University, College Station, TX 77843-3368, USA IITP, Moscow 127994, Russia Roumeng Lan Department of Mathematics, Texas A&M University, College Station, TX 77843-3368, USA Panayotis G. Kevrekidis Department of Mathematics and Statistics, University of Massachusetts, Amherst, MA 01003- 4515, USA Avadh Saxena Center for Nonlinear Studies and Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 1 2 J.
    [Show full text]
  • Global Questions for Evolution Equations Landau-Lifshitz Flow and Dirac Equation
    Global questions for evolution equations Landau-Lifshitz flow and Dirac equation by Meijiao Guan M.Sc., Huazhong University of Science and Technology, 2003 Ph.D., The University of British Columbia, 2009 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in The Faculty of Graduate Studies (Mathematics) THE UNIVERSITY OF BRITISH COLUMBIA (Vancouver) July 2009 © Meijiao Guan 2009 Abstract This thesis concerns the stationary solutions and their stability for some evolution equations from physics. For these equations, the basic questions regarding the solutions concern existence, uniqueness, stability and singular ity formation. In this thesis, we consider two different classes of equations: the Landau-Lifshitz equations, and nonlinear Dirac equations. There are two different definitions of stationary solutions. For the Landau-Lifshitz equation, the stationary solution is time-independent, while for the Dirac equation, the stationary solution, also called solitary wave solution or ground state solution, is a solution which propagates without changing its shape. The class of Landau-Lifshitz equations (including harmonic map heat flow and Schrödinger map equations) arises in the study of ferromagnets (and anti-ferromagnets), liquid crystals, and is also very natural from a geometric standpoint. Harmonic maps are the stationary solutions to these equations. My thesis concerns the problems of singularity formation vs. global regu larity and long time asymptotics when the target space is a 2-sphere. We consider maps with some symmetry. I show that for m-equivariant maps with energy close to the harmonic map energy, the solutions to Landau Lifshitz equations are global in time and converge to a specific family of harmonic maps for big m, while for m = 1, a finite time blow up solution is constructed for harmonic map heat flow.
    [Show full text]
  • Solitary Waves of a PT-Symmetric Nonlinear Dirac Equation Jesús Cuevas–Maraver University of Seville
    University of Massachusetts Amherst ScholarWorks@UMass Amherst Mathematics and Statistics Department Faculty Mathematics and Statistics Publication Series 2015 Solitary Waves of a PT-Symmetric Nonlinear Dirac Equation Jesús Cuevas–Maraver University of Seville Panayotis G. Kevrekidis University of Massachusetts Amherst Avadh Saxena Los Alamos National Laboratory Fred Cooper Santa Fe Institute Avinash Khare Savitribai Phule Pune University See next page for additional authors Follow this and additional works at: https://scholarworks.umass.edu/math_faculty_pubs Part of the Mathematics Commons Recommended Citation Cuevas–Maraver, Jesús; Kevrekidis, Panayotis G.; Saxena, Avadh; Cooper, Fred; Khare, Avinash; Comech, Andrew; and Bender, Carl M., "Solitary Waves of a PT-Symmetric Nonlinear Dirac Equation" (2015). Mathematics and Statistics Department Faculty Publication Series. 1244. Retrieved from https://scholarworks.umass.edu/math_faculty_pubs/1244 This Article is brought to you for free and open access by the Mathematics and Statistics at ScholarWorks@UMass Amherst. It has been accepted for inclusion in Mathematics and Statistics Department Faculty Publication Series by an authorized administrator of ScholarWorks@UMass Amherst. For more information, please contact [email protected]. Authors Jesús Cuevas–Maraver, Panayotis G. Kevrekidis, Avadh Saxena, Fred Cooper, Avinash Khare, Andrew Comech, and Carl M. Bender This article is available at ScholarWorks@UMass Amherst: https://scholarworks.umass.edu/math_faculty_pubs/1244 1 Solitary waves of a PT -symmetric Nonlinear Dirac equation Jesus´ Cuevas–Maraver, Panayotis G. Kevrekidis, Avadh Saxena, Fred Cooper, Avinash Khare, Andrew Comech and Carl M. Bender (Invited Paper) Abstract—In the present work, we consider a prototypical analogy of the paraxial approximation of Maxwell’s equations example of a PT -symmetric Dirac model.
    [Show full text]
  • Numerical Methods and Comparison for the Nonlinear Dirac Equation in the Nonrelativistic Limit Regime
    Numerical methods and comparison for the nonlinear Dirac equation in the nonrelativistic limit regime BAO WeiZhua, JIA XiaoWeia, YIN Jiab aDepartment of Mathematics, National University of Singapore, Singapore 119076, Singapore bNUS Graduate School for Integrative Sciences and Engineering (NGS), National University of Singapore, Singapore 117456, Singapore Abstract We present and analyze several numerical methods for the discretization of the nonlinear Dirac equation in the nonrelativistic limit regime, involving a small dimensionless parameter 0 <ε 1 which is inversely proportional to the speed of light. In this limit regime, the solution is highly oscillatory≪ in time, i.e. there are propagating waves with wavelength O(ε2) and O(1) in time and space, respectively. We begin with four frequently used finite difference time domain (FDTD) methods and establish rigorously error estimates for the FDTD methods, which depend explicitly on the mesh size h and time step τ as well as the small parameter 0 < ε 1. Based on the error bounds, in order to obtain ‘correct’ numerical solutions in the nonrelativistic limit≤ regime, i.e. 0 <ε 1, the FDTD methods share the same ε-scalability: τ = O(ε3) and h = O(√ε). Then we propose and analyze≪ two numerical methods for the discretization of the nonlinear Dirac equation by using the Fourier spectral discretization for spatial derivatives combined with the exponential wave integrator and time-splitting technique for temporal derivatives, respectively. Rigorous error bounds for the two numerical methods show that their ε-scalability is improved to τ = O(ε2) and h = O(1) when 0 < ε 1 compared with the FDTD methods.
    [Show full text]
  • SPECTRAL STABILITY of BI-FREQUENCY SOLITARY WAVES in SOLER and DIRAC–KLEIN–GORDON MODELS 1. Introduction. the Soler Model [
    COMMUNICATIONS ON doi:10.3934/cpaa.2018065 PURE AND APPLIED ANALYSIS Volume 17, Number 4, July 2018 pp. 1331–1347 SPECTRAL STABILITY OF BI-FREQUENCY SOLITARY WAVES IN SOLER AND DIRAC–KLEIN–GORDON MODELS Dedicated to Professor Vladimir Georgiev on the occasion of his sixtieth birthday NABILE BOUSSA¨ID Universite´ Bourgogne Franche-Comte´ 16 route de Gray, 25030 Besanc¸on CEDEX, France ANDREW COMECH∗ Texas A&M University, College Station, TX 77843, USA IITP, Moscow 127051, Russia St. Petersburg State University, St. Petersburg 199178, Russia ABSTRACT. We construct bi-frequency solitary waves of the nonlinear Dirac equation with the scalar self-interaction, known as the Soler model (with an arbitrary nonlinearity and in arbitrary dimension) and the Dirac–Klein–Gordon with Yukawa self-interaction. These solitary waves provide a natural implementation of qubit and qudit states in the the- ory of quantum computing. We show the relation of ±2!i eigenvalues of the linearization at a solitary wave, Bogo- liubov SU(1; 1) symmetry, and the existence of bi-frequency solitary waves. We show that the spectral stability of these waves reduces to spectral stability of usual (one-frequency) solitary waves. 1. Introduction. The Soler model [18, 29] is the nonlinear Dirac equation with the mini- mal scalar self-coupling, ∗ n N i@t = Dm − f( β )β ; x 2 R ; (t; x) 2 C ; (1) where f is a continuous real-valued function with f(0) = 0. Above, ¯ = ∗β, with ∗ the hermitian conjugate. This is one of the main models of the nonlinear Dirac equation, alongside with its own one-dimensional analogue, the Gross–Neveu model [17, 21], and with the massive Thirring model [30].
    [Show full text]