A stroboscopic averaging technique for highly-oscillatory Schrödinger equation
|
|
- Shauna Evans
- 5 years ago
- Views:
Transcription
1 A stroboscopic averaging technique for highly-oscillatory Schrödinger equation F. Castella IRMAR and INRIA-Rennes Joint work with P. Chartier, F. Méhats and A. Murua Saint-Malo Workshop, January 26-28, 2011
2 Goal of the present study Numerical schemes and asymptotic study for nonlinear Schrödinger equations of the form i t ψ ε (t, x) = H ( x ε ψε (t, x) + F ψ ε (t, x) 2) ψ ε (t, x), where the Hamiltonian H x is the harmonic oscillator (x R d ) H x = 1 2 x x 2. and F is a given nonlinear function. Context: Bose-Einstein condensates H x /ε represents strong confinement in x. nonlinear term describes bosons interactions.
3 Goal of the present study (2) References (in that context): Bao, Markowich, Schmeiser, Weishäupl (M3AS, 2005), Ben Abdallah, Méhats (JMPA, 2005), Ben Abdallah, Méhats, Schmeiser, Weishäupl (SIAM, 2005). See also, in the context of fluid mechanics: Grenier (JMPA, 1997), Schochet (JDE, 1994), Métivier, Schochet (JDE 2003), and, for other Schrödinger-like or Vlasov-like equations: Bidégaray, Castella, Degond (DCDS 2004), Castella, Degond, Goudon (JSP 2006).
4 Rough analysis (and wrong approach) Linear part: Eigenvalues of H x : {E n = n n d + d 2 ; n = (n 1,..., n d ) N d }, Eigenfunctions: explicitely known χ n (x) s (n N d ), Hermite polynomial Gaussian. Project ψ ε (t, x) = n ψ ε n (t)χ n(x), and write, say if F( ψ 2 )ψ = ψ 2 ψ, i t ψn ε (t) = E n ε ψε n (t) + A n,p,q,r ψp ε (t) ψq ε (t)ψε r (t) p,q,r ( ) with A n,p,q,r = ψp ε (t) ψq ε (t)ψε r (t). R d
5 Rough analysis (and wrong approach) (2) Filter out the oscillatory term and introduce ( φ ε n (t, x) = exp +it E ) n ψ ε (t, x). ε It satisfies i t φ ε n(t) = p,q,r A n,p,q,r exp ( it E ) n + E p E q E r φ ε ε p(t) φ ε q(t)φ ε r (t). System of the form ( ) t t u ε = G ε, uε, with G periodic (E n s are half-integers) and u ε = (φ ε 0,φε 1,...) (infinite, nonlinearly coupled system).
6 Rough analysis (and wrong approach) (3) Questions How to numerically average out the ODE ( ) t t u ε = G ε, uε? Answer inspired by Chartier, Murua, Sanz-Serna (ODE s). How to control the norms of u ε, and that of the nonlinear term G (i.e., typically φ ε and the sums p,q,r : it is not clear that the Sobolev smoothness n nα φ n 2 < + implies n nα p,q,r 2 < +, because of the An,p,q,r s)? Answer inspired by Ben Abdallah, Castella, Méhats (functional analytic framework).
7 Right approach and analytic framework Avoid projecting on the χ n s and write directly (Ben Abdallah, Castella, Méhats) ( φ ε (t, x) = exp +it H ) x ψ ε (t, x). ε ( ) t It satisfies i t φ ε (t, x) = G ε,φε (t, x), ( F(z) = F( z 2 ) z) where, for u = u(x) given, we have ( ) ( t G ε, u = exp +it H ) ( x F exp ε ( it H x ε ) ) u(x)
8 Right approach and analytic framework (2) Bounds on the nonlinear terms G(t/ε,φ ε ) or F(ψ ε )? When F 0, we have (take a large N N) i t ψ ε = H x ε ψε while = t H N x ψ ε 2 = 2Re L 2 H N x ( i H ) x ε ψε, Hx L N ψε = 0, 2 (H x commutes with Hx N and H x is self-adjoint) = Hx N ψ ε 2 const. L 2 t N x ψ ε 2 L 2 = O ( ) 1 ε = no uniform in ε bound on N x ψε 2 L 2.
9 Right approach and analytic framework (3) This is due to t N x ψ ε 2 = 2Re = O ( 1 ε N x = 2Re L 2 ). N x ( i x 2 ε ψε ( i H x ε ψε ), x N ψ ε ), x N ψ ε Only a uniform bound on H N x ψε L 2 (and not on N x ψε L 2) can be obtained in the general case. The good norm is u B N = u L 2 + H N x u L 2, instead of the Sobolev scale u H N = u L 2 + N x u L 2. L 2 L 2
10 Right approach and analytic framework (4) Note that : ( Hx N L φε = H N 2 x exp +it H ) x L ψ ε ε 2 ( = exp +it H ) x L ε HN x ψ ε (commutation) 2 = Hx N ψ ε L (self-adjointness) 2 Remaining question: we know ( u H N := u L 2 + x Nu L 2) F(ψ ε ) const. ( ψε H N), H N yet is it true that ( u B N := u L 2 + Hx Nu L 2) F(ψ ε ) const. ( ψε B N B N)?
11 Right approach and analytic framework (5) Answer (Ben Abdallah, Castella, Méhats, see also Helffer, Nier, and Bony, Chemin requires Weyl-Hörmander calculus for general Hamiltonians): Yes because of the following Theorem (H x = x + x 2 ) For any N, we have the equivalence of norms u L 2 + 2N x u L 2 + x 2N u L 2 u L 2 + ( 2 x + x 2 ) N u L 2 or, in other words u L 2 + 2N x u L 2 + x 2N u L 2 u B N.
12 Right approach and analytic framework (6) By standard nonlinear analysis (Gronwall + the above nonlinear estimates on F(u) B N + standard fixed point in the original PDE), this immediately implies the existence of T 0 > 0 such that ψ ε (t, x) B N, φ ε (t, x) B N const. (0 t T 0 ), and all nonlinear terms are well defined and uniformly bounded in the space B N there remains to average out, in the space B N, the equation ( ) t i t φ ε (t, x) = G ε,φε (t, x)
13 Highly oscillatory ODE s Problems of the form ẏ = dy dt = f(y, t/ε), y(0) = y 0, where : ε is a small parameter (the inverse of the frequency). f is a smooth vector-field, 2π-periodic in t/ε. Assumptions: a clear and explicit separation of the time-scales a periodic (and not only quasi-periodic at this stage) dependence in the fast time
14 Highly oscillatory ODE s (2) Appears in various forms: Standard form in KAM theory Two-scale expansions in quantum mechanics (Wentzel-Kramers-Brillouin, 1926) Magnus expansions (A. Iserles, 2002) Modulated Fourier expansions (D. Cohen, E. Hairer and Chr. Lubich, 2003) Theory has been gradually improved: Krylov and Bogoliubov (1934) : basic idea Bogoliubov and Mitropolski (1958) : rigorous statement for second order approximation and general scheme Perko (1969) : almost complete theory with error estimates for the periodic and quasi-periodic cases (see also the book Sanders, Verhulst and Murdock, 2007)
15 Highly oscillatory ODE s (3) Theorem (Perko, SIAM 1969 : Averaging of high-order) Under a smoothness assumption on f and a non-resonance condition on ω, given the system { y = εf(y,θ) R n, y(0) = y 0, θ = ω T d, θ(0) = θ 0, there exists a transformation from R n T d to R n such that y = U(Y,θ) = Y + εu 1 (Y,θ) ε k u k (Y,θ) Y = εf 1 (Y) ε k F k (Y), Y(0) = ξ. and y(τ) U(Y(τ),θ 0 + τω) Cε k for τ C/ε.
16 Highly oscillatory ODE s (4) The functions u i (and thus F i ) are not unique except F 1 F 1 (Y) = f(y,θ), j 1 [ 1 F k f ( ) j (Y,θ) = k! y k u i1,..., u ik u ] k Y F j k, k=1 i i k =j 1 F j (Y) = 1 (2π) d Fj (Y,θ)dθ, T d ω u j θ (Y,θ) = F j (Y,θ) F j (Y). For d = 1, one can impose u j (Y, 0) = 0: this is stroboscopic averaging, in the sense that U(Y, 2kπ) = Y.
17 An example For f(y,τ) = 1 + cos(τ)y, we have Φ τ τ 0 (y 0 ) = εe ε sin(τ) τ so that τ 0 e ε sin(s) ds + e ε(sin(τ) sin(τ 0)) y 0, Φ τ 0+2π τ 0 (y 0 ) = 2πενe ε sin(τ 0) + y 0, ν = 1 2π ( We see that Φ τ 0+2π τ 0 2π 0 e ε sin(τ) dτ. ) k (y0 ) coincides with the solution at points τ = 2kπ of the differential equation { Y (τ) = ενe ε sin(τ 0) Y(0) = y 0, which happens to be the modified equation we are looking for.
18 An example (2) Example 2.4 Exact solutions versus their averaged versions y and Y Time We search for a differential equation whose flow coincides at times that are multiple of 2π with the flow of the highly-oscillatory equation.
19 B-series expansions Mode-coloured B-series with time-dependent coefficients Consider the rescaled system expanded in Fourier coefficients { y (τ) = ε k Z eikτ f k (y), y(0) = y 0 θ (τ) = 1, θ(0) = θ 0 where the f k s are the Fourier coefficients of f. Chartier, Murua and Sanz-Serna [FOCM 2010] consider B-series of the form where y 0 + u T ε u α u(τ) σ u F u (y 0 ) T is a set of mode-decorated trees. F u (y 0 ) are associated derivatives of the f k s. σ u is a normalisation coefficient. α u (τ) are coefficients defining the series.
20 B-series expansions (2) Tree Order σ u F u k 1 1 f k l k = [ k ] l 2 1 f l f k k l = [ k, l ] m δ k,l f m(f k, f l ) m k k l [ l ] m = m 3 1 f mf l f k Table: Mode-decorated trees and their associated coefficients
21 B-series expansions (3) Nota : The general procedure is described as in Chartier, Murua, Sanz-Serna, FOCM 2010 Construction of averaged vector field F from εf(y, τ): 1 Compute the B-series expansion of the exact solution y(τ) = y 0 + u T ε u α u(τ) σ u F u (y 0 ) 2 Notice that α u (τ) = P u (τ,θ) θ=τ with P u (τ,θ) = a u kl τ l e ikθ 3 Define averaged (or frozen) solution by ᾱ u (τ) = P u (τ, 0) Y(τ) = y 0 + u T ε u ᾱu(τ) σ u F u (y 0 ) 4 Define averaged vector field εf(y,τ) such that { Y (τ) = εf(y(τ),τ) Y(0) = y 0 5 Notice that F is autonomous and Hamiltonian as soon as f is Hamiltonian.
22 The averaged vector field is autonomous The vector field corresponding to Y(τ) = Φ τ τ 0 (y 0 ) writes F(Y,τ) = ( Φ τ τ 0 τ ( Φτ τ0 ) 1)(Y) = u T where β(τ) = ᾱ 1 (τ) ᾱ (τ). Since Y(τ j+k ) = ε β u u (τ) F u (Y) σ u ( Φ τ 0+2π τ 0 ) j+k (y0 ): ᾱ u (τ j + τ k ) = ᾱ u (τ k ) ᾱ u (τ j ) = ᾱ u (τ j ) ᾱ u (τ k ), j, k = 0, ±1, ±2,... and by an interpolation argument k, τ, ᾱ u (τ + τ k ) = ᾱ u (τ k ) ᾱ u (τ). Differentiating and taking the value at τ = 0 gives k, ᾱ (τ k ) = ᾱ(τ k ) ᾱ (0), so that by interpolation once again, τ, ᾱ (τ) = ᾱ(τ) ᾱ (0) i.e. τ, ᾱ 1 (τ) ᾱ (τ) = ᾱ (0).
23 The averaged vector field is Hamiltonian If α is symplectic, the B-series ᾱ remains symplectic Since α u (τ) = P u (τ,τ) with P u (τ,θ) := a u lk τ l e ikθ we see that the symplecticity relation τ, α u v (τ) + α v u (τ) = α u (τ)α v (τ) implies τ, θ, P u v (τ,θ) + P v u (τ,θ) = P u (τ,θ)p v (τ,θ) and hence τ, ᾱ u v (τ) + ᾱ v u (τ) = ᾱ u (τ)ᾱ v (τ)
24 The numerical scheme The idea is to solve the averaged equation at two levels, in the spirit of Heterogeneous Multiscale Methods: Approximate the averaged vector field F by central differences of the form F(Y) 1 4π (Φ2π 0 (Y) Φ 2π 0 (Y)) with a numerical method with constant stepsizes (micro-steps). Solve the averaged equation by a numerical method with possibly variable stepsizes (macro-steps) Nota : in the formulation ẏ = f(t/ε, y), this gives rise to micro-steps δt ε alternating with macro-steps ε t 1.
25 p Stroboscopic averaging for nonlinear Schrödinger equations The numerical scheme (2) Example (Van Der Pol oscillator) 2.5 Van Der Pol oscillator { q = p ṗ = q + ε(1 q 2 )p q which, after the change of variables writes q = cos(t)x + sin(t)y p = sin(t)x + cos(t)y j ẋ = sin (t) ε `1 (cos (t) x + sin (t) y) 2 ( sin (t) x + cos (t) y) ẏ = cos (t) ε `1 2 (cos (t) x + sin (t) y) ( sin (t) x + cos (t) y)
26 Stroboscopic averaging for nonlinear Schro dinger equations The numerical scheme (3) Van Der Pol oscillator p q Figure: Exact solution (red) and numerical solution obtained by SAM (blue)
27 The numerical scheme : references References P. Chartier, J.M. Sanz-Serna and A. Murua, Higher-order averaging, formal series and numerical integration I: B-series, FOCM, M.P. Calvo, P. Chartier, J.M. Sanz-Serna and A. Murua, A stroboscopic numerical method for highly oscillatory problems, submitted. P. Chartier, J.M. Sanz-Serna and A. Murua, Higher-order averaging, formal series and numerical integration II: the multi-frequency case, in preparation.
28 Application to the considered nonlinear Schrödinger equations Theorem The Perko formulae apply in our context, when conveniently rephrased in the B N spaces: one can define functions u j, Y j as in Perko, and averaging can be performed at any order, to provided φ ε K (t, x) s such that φ ε K (t, x) = φε (t, x) + O(ε K ) in B N, whenever 0 t T 0. Numerical counterpart The above described numerical scheme, when conveniently adapted, provides the following results.
29 Application to the considered nonlinear Schrödinger equations (2) In the case of i t ψ ε = 1 ε x + ψ ε 2 ψ ε x [0, 1] with periodic boundary conditions, t [0, 2/2π] error eps=4.4e 2 eps=1.1e 2 eps=2.8e eps=6.9e 4 eps=1.7e 4 eps=8.6e time step Figure: direct Strang splitting (order 2 if ε 1)
30 Application to the considered nonlinear Schrödinger equations (3) error eps=4.4e 2 eps=1.1e 2 eps=2.8e 3 eps=6.9e 4 eps=1.7e 4 eps=8.6e time step Figure: direct order 4 splitting
31 Application to the considered nonlinear Schrödinger equations (4) error eps=4.4e 2 eps=1.1e 2 eps=2.8e 3 eps=6.9e 4 eps=1.7e 4 eps=8.6e time step Figure: direct order 6 splitting
32 Application to the considered nonlinear Schrödinger equations (5) eps=2.2e 2 eps=5.5e 3 eps=1.4e 3 eps=3.5e 4 eps=8.6e error micro step Figure: Stroboscopic averaging constant numerical cost, various values of ε
33 Application to the considered nonlinear Schrödinger equations (6) energy time Figure: Energy conservation.
34 Application to the considered nonlinear Schrödinger equations (7) u p p=1 p=2 p=3 p=4 p=5 p=6 p=7 p= time Figure: Evolution of modes 1 to 8 in the cubic NLS: modes 1 and 5 feed mode 7 (and only this mode).
Quasi-stroboscopic averaging: from B-series to numerical methods?
Outline Quasi-stroboscopic averaging: from B-series to numerical methods? Philippe Chartier INRIA, ENS Cachan Bruz, IRMAR, University of Rennes I Joint work with Ander Murua and Jesus-Maria Sanz-Serna
More informationQuasi-stroboscopic averaging: The non-autonomous case
Outline Quasi-stroboscopic averaging: The non-autonomous case Philippe Chartier INRIA and ENS Cachan Bruz Joint work with Ander Murua and Jesus-Maria Sanz-Serna Geometric Numerical Integration, Oberwolfach
More informationUniformly accurate averaging numerical schemes for oscillatory evolution equations
Uniformly accurate averaging numerical schemes for oscillatory evolution equations Philippe Chartier University of Rennes, INRIA Joint work with M. Lemou (University of Rennes-CNRS), F. Méhats (University
More informationNumerical stroboscopic averaging for ODEs and DAEs
Numerical stroboscopic averaging for ODEs and DAEs M. P. Calvo Universidad de Valladolid, Spain Joint work with Ph. Chartier, A. Murua, J. M. Sanz-Serna 1 I. HIGHLY OSCILLATORY PROBLEMS 2 Consider the
More informationThree approaches for the design of adaptive time-splitting methods
Three approaches for the design of adaptive time-splitting methods Mechthild Thalhammer Leopold Franzens Universität Innsbruck, Austria Workshop on Geometric Integration and Computational Mechanics Organisers:
More informationStroboscopic averaging for the nonlinear Schrödinger equation
Stroboscopic averaging for the nonlinear Schrödinger equation F. Castella, Ph. Chartier, F. Méhats and A. Murua July 25, 213 Abstract In this paper, we are concerned with an averaging procedure, namely
More informationKAM for NLS with harmonic potential
Université de Nantes 3rd Meeting of the GDR Quantum Dynamics MAPMO, Orléans, 2-4 February 2011. (Joint work with Benoît Grébert) Introduction The equation : We consider the nonlinear Schrödinger equation
More informationStroboscopic Averaging for the Nonlinear Schrödinger Equation
Stroboscopic Averaging for the Nonlinear Schrödinger Equation François Castella, Philippe Chartier, Florian Méhats, Ander Murua To cite this version: François Castella, Philippe Chartier, Florian Méhats,
More informationAVERAGING OF NONLINEAR SCHRÖDINGER EQUATIONS WITH STRONG MAGNETIC CONFINEMENT
AVERAGING OF NONLINEAR SCHRÖDINGER EQUATIONS WITH STRONG MAGNETIC CONFINEMENT RUPERT L. FRANK, FLORIAN MÉHATS, AND CHRISTOF SPARBER Abstract. We consider the dynamics of nonlinear Schrödinger equations
More informationModeling, Analysis and Simulation for Degenerate Dipolar Quantum Gas
Modeling, Analysis and Simulation for Degenerate Dipolar Quantum Gas Weizhu Bao Department of Mathematics National University of Singapore Email: matbaowz@nus.edu.sg URL: http://www.math.nus.edu.sg/~bao
More informationPeriodic oscillations in the Gross-Pitaevskii equation with a parabolic potential
Periodic oscillations in the Gross-Pitaevskii equation with a parabolic potential Dmitry Pelinovsky 1 and Panos Kevrekidis 2 1 Department of Mathematics, McMaster University, Hamilton, Ontario, Canada
More informationZ. Zhou On the classical limit of a time-dependent self-consistent field system: analysis. computation
On the classical limit of a time-dependent self-consistent field system: analysis and computation Zhennan Zhou 1 joint work with Prof. Shi Jin and Prof. Christof Sparber. 1 Department of Mathematics Duke
More informationON THE GROSS-PITAEVSKII EQUATION WITH STRONGLY ANISOTROPIC CONFINEMENT: FORMAL ASYMPTOTICS AND NUMERICAL EXPERIMENTS
Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company ON THE GROSS-PITAEVSKII EQUATION WITH STRONGLY ANISOTROPIC CONFINEMENT: FORMAL ASYMPTOTICS AND NUMERICAL EXPERIMENTS
More informationHigh-order time-splitting methods for Schrödinger equations
High-order time-splitting methods for Schrödinger equations and M. Thalhammer Department of Mathematics, University of Innsbruck Conference on Scientific Computing 2009 Geneva, Switzerland Nonlinear Schrödinger
More informationNormal form for the non linear Schrödinger equation
Normal form for the non linear Schrödinger equation joint work with Claudio Procesi and Nguyen Bich Van Universita di Roma La Sapienza S. Etienne de Tinee 4-9 Feb. 2013 Nonlinear Schrödinger equation Consider
More informationDegree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m.
Degree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m. Candidates should submit answers to a maximum of four
More informationA Minimal Uncertainty Product for One Dimensional Semiclassical Wave Packets
A Minimal Uncertainty Product for One Dimensional Semiclassical Wave Packets George A. Hagedorn Happy 60 th birthday, Mr. Fritz! Abstract. Although real, normalized Gaussian wave packets minimize the product
More informationExponential integrators
Exponential integrators Marlis Hochbruck Heinrich-Heine University Düsseldorf, Germany Alexander Ostermann University Innsbruck, Austria Helsinki, May 25 p.1 Outline Time dependent Schrödinger equations
More informationMagnetic wells in dimension three
Magnetic wells in dimension three Yuri A. Kordyukov joint with Bernard Helffer & Nicolas Raymond & San Vũ Ngọc Magnetic Fields and Semiclassical Analysis Rennes, May 21, 2015 Yuri A. Kordyukov (Ufa) Magnetic
More informationPhilippe Chartier 1, Florian Méhats 2, Mechthild Thalhammer 3 and Yong Zhang Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher CONVERGENCE OF MULTI-REVOLUTION COMPOSITION TIME-SPLITTING METHODS FOR HIGHLY OSCILLATORY
More informationGroup Method. December 16, Oberwolfach workshop Dynamics of Patterns
CWI, Amsterdam heijster@cwi.nl December 6, 28 Oberwolfach workshop Dynamics of Patterns Joint work: A. Doelman (CWI/UvA), T.J. Kaper (BU), K. Promislow (MSU) Outline 2 3 4 Interactions of localized structures
More informationStochastic nonlinear Schrödinger equations and modulation of solitary waves
Stochastic nonlinear Schrödinger equations and modulation of solitary waves A. de Bouard CMAP, Ecole Polytechnique, France joint work with R. Fukuizumi (Sendai, Japan) Deterministic and stochastic front
More informationThe Density Matrix for the Ground State of 1-d Impenetrable Bosons in a Harmonic Trap
The Density Matrix for the Ground State of 1-d Impenetrable Bosons in a Harmonic Trap Institute of Fundamental Sciences Massey University New Zealand 29 August 2017 A. A. Kapaev Memorial Workshop Michigan
More informationExponential integrators for oscillatory second-order differential equations
Exponential integrators for oscillatory second-order differential equations Marlis Hochbruck and Volker Grimm Mathematisches Institut Heinrich Heine Universität Düsseldorf Cambridge, March 2007 Outline
More informationAnton ARNOLD. with N. Ben Abdallah (Toulouse), J. Geier (Vienna), C. Negulescu (Marseille) TU Vienna Institute for Analysis and Scientific Computing
TECHNISCHE UNIVERSITÄT WIEN Asymptotically correct finite difference schemes for highly oscillatory ODEs Anton ARNOLD with N. Ben Abdallah (Toulouse, J. Geier (Vienna, C. Negulescu (Marseille TU Vienna
More informationSmilansky-Solomyak model with a δ -interaction
Smilansky-Solomyak model with a δ -interaction Jiří Lipovský University of Hradec Králové, Faculty of Science jiri.lipovsky@uhk.cz joint work with P. Exner Ostrava, April 6, 2018 Jiří Lipovský Smilansky-Solomyak
More informationApplications of the periodic unfolding method to multi-scale problems
Applications of the periodic unfolding method to multi-scale problems Doina Cioranescu Université Paris VI Santiago de Compostela December 14, 2009 Periodic unfolding method and multi-scale problems 1/56
More informationWell-Posedness and Adiabatic Limit for Quantum Zakharov System
Well-Posedness and Adiabatic Limit for Quantum Zakharov System Yung-Fu Fang (joint work with Tsai-Jung Chen, Jun-Ichi Segata, Hsi-Wei Shih, Kuan-Hsiang Wang, Tsung-fang Wu) Department of Mathematics National
More informationBose-Einstein condensation and limit theorems. Kay Kirkpatrick, UIUC
Bose-Einstein condensation and limit theorems Kay Kirkpatrick, UIUC 2015 Bose-Einstein condensation: from many quantum particles to a quantum superparticle Kay Kirkpatrick, UIUC/MSRI TexAMP 2015 The big
More informationDerivation of quantum hydrodynamic equations with Fermi-Dirac and Bose-Einstein statistics
Derivation of quantum hydrodynamic equations with Fermi-Dirac and Bose-Einstein statistics Luigi Barletti (Università di Firenze) Carlo Cintolesi (Università di Trieste) 6th MMKT Porto Ercole, june 9th
More informationLIMITATIONS OF THE METHOD OF MULTIPLE-TIME-SCALES
LIMITATIONS OF THE METHOD OF MULTIPLE-TIME-SCALES Peter B. Kahn Department of Physics and Astronomy State University of New York Stony Brook, NY 794 and Yair Zarmi Department for Energy & Environmental
More informationKAM for quasi-linear KdV
KAM for quasi-linear KdV Massimiliano Berti ST Etienne de Tinée, 06-02-2014 KdV t u + u xxx 3 x u 2 + N 4 (x, u, u x, u xx, u xxx ) = 0, x T Quasi-linear Hamiltonian perturbation N 4 := x {( u f )(x, u,
More information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More informationWORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS. Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid
WORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid 1 I. OVERVIEW 2 The importance of splitting integrators
More informationWeakly non-linear completely resonant hamiltonian PDEs and the problem of weak turbulence
Sergei Kuksin Weakly non-linear completely resonant hamiltonian PDEs and the problem of weak turbulence (Toronto, 10 January, 2014 ) 1 1 Introduction: weak turbulence (WT) (One of) origines: Rudolf Peierls,
More informationFrom the N-body problem to the cubic NLS equation
From the N-body problem to the cubic NLS equation François Golse Université Paris 7 & Laboratoire J.-L. Lions golse@math.jussieu.fr Los Alamos CNLS, January 26th, 2005 Formal derivation by N.N. Bogolyubov
More informationExponential multistep methods of Adams-type
Exponential multistep methods of Adams-type Marlis Hochbruck and Alexander Ostermann KARLSRUHE INSTITUTE OF TECHNOLOGY (KIT) 0 KIT University of the State of Baden-Wuerttemberg and National Laboratory
More informationMetastability for the Ginzburg Landau equation with space time white noise
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz Metastability for the Ginzburg Landau equation with space time white noise Barbara Gentz University of Bielefeld, Germany
More informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
More informationLecture 4: Birkhoff normal forms
Lecture 4: Birkhoff normal forms Walter Craig Department of Mathematics & Statistics Waves in Flows - Lecture 4 Prague Summer School 2018 Czech Academy of Science August 31 2018 Outline Two ODEs Water
More informationSemi-Classical Dynamics Using Hagedorn Wavepackets
Semi-Classical Dynamics Using Hagedorn Wavepackets Leila Taghizadeh June 6, 2013 Leila Taghizadeh () Semi-Classical Dynamics Using Hagedorn Wavepackets June 6, 2013 1 / 56 Outline 1 The Schrödinger Equation
More informationHigher order asymptotic analysis of the nonlinear Klein- Gordon equation in the non-relativistic limit regime
Higher order asymptotic analysis of the nonlinear Klein- Gordon equation in the non-relativistic limit regime Yong Lu and Zhifei Zhang Abstract In this paper, we study the asymptotic behavior of the Klein-Gordon
More informationSimple Harmonic Oscillator
Classical harmonic oscillator Linear force acting on a particle (Hooke s law): F =!kx From Newton s law: F = ma = m d x dt =!kx " d x dt + # x = 0, # = k / m Position and momentum solutions oscillate in
More informationA Time Splitting for the Semiclassical Schrödinger Equation
A Time Splitting for the Semiclassical Schrödinger Equation V. Gradinaru Seminar for Applied Mathematics, ETH Zürich, CH 809, Zürich, Switzerland vasile.gradinaru@sam.math.ethz.ch and G. A. Hagedorn Department
More informationLocal smoothing and Strichartz estimates for manifolds with degenerate hyperbolic trapping
Local smoothing and Strichartz estimates for manifolds with degenerate hyperbolic trapping H. Christianson partly joint work with J. Wunsch (Northwestern) Department of Mathematics University of North
More informationNonlinear Modulational Instability of Dispersive PDE Models
Nonlinear Modulational Instability of Dispersive PDE Models Jiayin Jin, Shasha Liao, and Zhiwu Lin Georgia Tech ICERM workshop on water waves, 4/28/2017 Jiayin Jin, Shasha Liao, and Zhiwu Lin Georgia Tech
More informationDIMENSION REDUCTION FOR DIPOLAR BOSE-EINSTEIN CONDENSATES IN THE STRONG INTERACTION REGIME. Weizhu Bao. Loïc Le Treust and Florian Méhats
Kinetic and Related Models doi:1.3934/krm.21722 c American Institute of Mathematical Sciences Volume 1, Number 3, September 217 pp. 553 571 DIMENSION REDUCTION FOR DIPOLAR BOSE-EINSTEIN CONDENSATES IN
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
More informationNumerical aspects of the nonlinear Schro dinger equation. Christophe Besse
Numerical aspects of the nonlinear Schro dinger equation by Christophe Besse Laboratoire Paul Painleve, Universite Lille 1, CNRS, Inria Simpaf Team Lille Outline of the talk 1 Motivation 2 Numerical methods
More informationPY 351 Modern Physics - Lecture notes, 3
PY 351 Modern Physics - Lecture notes, 3 Copyright by Claudio Rebbi, Boston University, October 2016. These notes cannot be duplicated and distributed without explicit permission of the author. Time dependence
More informationFavourable time integration methods for. Non-autonomous evolution equations
Favourable time integration methods for non-autonomous evolution equations International Conference on Scientific Computation and Differential Equations (SciCADE) Session on Numerical integration of non-autonomous
More informationDimer with gain and loss: Integrability and PT -symmetry restoration
Dimer with gain and loss: Integrability and PT -symmetry restoration I Barashenkov, University of Cape Town Joint work with: Dima Pelinovsky (McMaster University, Ontario) Philippe Dubard (University of
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationStochastic homogenization 1
Stochastic homogenization 1 Tuomo Kuusi University of Oulu August 13-17, 2018 Jyväskylä Summer School 1 Course material: S. Armstrong & T. Kuusi & J.-C. Mourrat : Quantitative stochastic homogenization
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More informationNote on semiclassical analysis for nonlinear Schrödinger equations. Satoshi Masaki
Note on semiclassical analysis for nonlinear Schrödinger equations Satoshi Masaki Preface This article is a de facto doctoral thesis of the author. This consists of the summary of the author s research
More informationPeriodic solutions for a class of nonlinear partial differential equations in higher dimension
Periodic solutions for a class of nonlinear partial differential equations in higher dimension Guido Gentile 1 and Michela Procesi 2 1 Dipartimento di Matematica, Università di Roma Tre, Roma, I-00146,
More informationequations in the semiclassical regime
Adaptive time splitting for nonlinear Schrödinger equations in the semiclassical regime Local error representation and a posteriori error estimator Workshop, Numerical Analysis of Evolution Equations Vill,
More informationarxiv:math/ v1 [math.ap] 28 Oct 2005
arxiv:math/050643v [math.ap] 28 Oct 2005 A remark on asymptotic completeness for the critical nonlinear Klein-Gordon equation Hans Lindblad and Avy Soffer University of California at San Diego and Rutgers
More informationEstimating Accuracy in Classical Molecular Simulation
Estimating Accuracy in Classical Molecular Simulation University of Illinois Urbana-Champaign Department of Computer Science Institute for Mathematics and its Applications July 2007 Acknowledgements Ben
More informationLecture 3 Dynamics 29
Lecture 3 Dynamics 29 30 LECTURE 3. DYNAMICS 3.1 Introduction Having described the states and the observables of a quantum system, we shall now introduce the rules that determine their time evolution.
More informationEffective dynamics of many-body quantum systems
Effective dynamics of many-body quantum systems László Erdős University of Munich Grenoble, May 30, 2006 A l occassion de soixantiéme anniversaire de Yves Colin de Verdiére Joint with B. Schlein and H.-T.
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationHarmonic Oscillator. Robert B. Griffiths Version of 5 December Notation 1. 3 Position and Momentum Representations of Number Eigenstates 2
qmd5 Harmonic Oscillator Robert B. Griffiths Version of 5 December 0 Contents Notation Eigenstates of the Number Operator N 3 Position and Momentum Representations of Number Eigenstates 4 Coherent States
More informationBlow-up on manifolds with symmetry for the nonlinear Schröding
Blow-up on manifolds with symmetry for the nonlinear Schrödinger equation March, 27 2013 Université de Nice Euclidean L 2 -critical theory Consider the one dimensional equation i t u + u = u 4 u, t > 0,
More informationSymmetry of the Dielectric Tensor
Symmetry of the Dielectric Tensor Curtis R. Menyuk June 11, 2010 In this note, I derive the symmetry of the dielectric tensor in two ways. The derivations are taken from Landau and Lifshitz s Statistical
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationAngular momentum. Quantum mechanics. Orbital angular momentum
Angular momentum 1 Orbital angular momentum Consider a particle described by the Cartesian coordinates (x, y, z r and their conjugate momenta (p x, p y, p z p. The classical definition of the orbital angular
More informationHomogenization method for elastic materials
Outline 1 Introduction 2 Description of the geometry 3 Problem setting 4 Implementation & Results 5 Bandgaps 6 Conclusion Introduction Introduction study of the homogenization method applied on elastic
More informationanalysis for transport equations and applications
Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg
More informationThe Mathematics of Nonlinear Optics
CHAPTER 3 The Mathematics of Nonlinear Optics Guy Métivier Université Bordeaux 1, IMB UMR 5251, 351 Cours de la Libération, 33405 Talence Cedex, France E-mail: guy.metivier@math.u-bordeaux1.fr Contents
More informationQuantum mechanics. Chapter The quantum mechanical formalism
Chapter 5 Quantum mechanics 5.1 The quantum mechanical formalism The realisation, by Heisenberg, that the position and momentum of a quantum mechanical particle cannot be measured simultaneously renders
More informationharmonic oscillator in quantum mechanics
Physics 400 Spring 016 harmonic oscillator in quantum mechanics lecture notes, spring semester 017 http://www.phys.uconn.edu/ rozman/ourses/p400_17s/ Last modified: May 19, 017 Dimensionless Schrödinger
More informationH = ( H(x) m,n. Ω = T d T x = x + ω (d frequency shift) Ω = T 2 T x = (x 1 + x 2, x 2 + ω) (skewshift)
Chapter One Introduction We will consider infinite matrices indexed by Z (or Z b ) associated to a dynamical system in the sense that satisfies H = ( H(x) m,n )m,n Z H(x) m+1,n+1 = H(T x) m,n where x Ω,
More informationSplittings, commutators and the linear Schrödinger equation
UNIVERSITY OF CAMBRIDGE CENTRE FOR MATHEMATICAL SCIENCES DEPARTMENT OF APPLIED MATHEMATICS & THEORETICAL PHYSICS Splittings, commutators and the linear Schrödinger equation Philipp Bader, Arieh Iserles,
More informationMulti-revolution composition methods for highly oscillatory differential equations
Multi-revolution composition methods for highly oscillatory differential equations Philippe Chartier 1, Joseba Makazaga 2, Ander Murua 2, and Gilles Vilmart 3 March 1, 2013 Abstract We introduce a new
More informationThe Dirichlet problem for non-divergence parabolic equations with discontinuous in time coefficients in a wedge
The Dirichlet problem for non-divergence parabolic equations with discontinuous in time coefficients in a wedge Vladimir Kozlov (Linköping University, Sweden) 2010 joint work with A.Nazarov Lu t u a ij
More informationModified Equations for Variational Integrators
Modified Equations for Variational Integrators Mats Vermeeren Technische Universität Berlin Groningen December 18, 2018 Mats Vermeeren (TU Berlin) Modified equations for variational integrators December
More informationProblems and Multiple Choice Questions
Problems and Multiple Choice Questions 1. A momentum operator in one dimension is 2. A position operator in 3 dimensions is 3. A kinetic energy operator in 1 dimension is 4. If two operator commute, a)
More informationLong-term behavior of the numerical solution of the cubic non-linear Scrödinger equation using Strang splitting method
Long-term behavior of the numerical solution of the cubic non-linear Scrödinger equation using Strang splitting method Marianna Khanamiryan a, Olavi Nevanlinna a, Tiina Vesanen a, a Department of Mathematics
More informationHomogenization of micro-resonances and localization of waves.
Homogenization of micro-resonances and localization of waves. Valery Smyshlyaev University College London, UK July 13, 2012 (joint work with Ilia Kamotski UCL, and Shane Cooper Bath/ Cardiff) Valery Smyshlyaev
More informationHilbert Space Problems
Hilbert Space Problems Prescribed books for problems. ) Hilbert Spaces, Wavelets, Generalized Functions and Modern Quantum Mechanics by Willi-Hans Steeb Kluwer Academic Publishers, 998 ISBN -7923-523-9
More information3 Generation and diffusion of vorticity
Version date: March 22, 21 1 3 Generation and diffusion of vorticity 3.1 The vorticity equation We start from Navier Stokes: u t + u u = 1 ρ p + ν 2 u 1) where we have not included a term describing a
More informationBrief review of Quantum Mechanics (QM)
Brief review of Quantum Mechanics (QM) Note: This is a collection of several formulae and facts that we will use throughout the course. It is by no means a complete discussion of QM, nor will I attempt
More informationAsymptotic Analysis of the Approximate Control for Parabolic Equations with Periodic Interface
Asymptotic Analysis of the Approximate Control for Parabolic Equations with Periodic Interface Patrizia Donato Université de Rouen International Workshop on Calculus of Variations and its Applications
More informationBackward error analysis
Backward error analysis Brynjulf Owren July 28, 2015 Introduction. The main source for these notes is the monograph by Hairer, Lubich and Wanner [2]. Consider ODEs in R d of the form ẏ = f(y), y(0) = y
More informationThe Schrödinger equation with spatial white noise potential
The Schrödinger equation with spatial white noise potential Arnaud Debussche IRMAR, ENS Rennes, UBL, CNRS Hendrik Weber University of Warwick Abstract We consider the linear and nonlinear Schrödinger equation
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationAnalytic invariants associated with a parabolic fixed point in C 2
Ergod. Th. & Dynam. Sys. (2008), 28, 1815 1848 doi:10.1017/s0143385707001046 c 2008 Cambridge University Press Printed in the United Kingdom Analytic invariants associated with a parabolic fixed point
More informationTWO-SCALE ASYMPTOTIC EXPANSION : BLOCH-FLOQUET THEORY IN PERIODIC STRUCTURES
TWO-SCALE ASYMPTOTIC EXPANSION : BLOCH-FLOQUET THEORY IN PERIODIC STRUCTURES JACK ARBUNICH 1. Setting of Periodic Structures Our aim is to study an application of Bloch-Floquet Theory in the multiscale
More informationCollection of formulae Quantum mechanics. Basic Formulas Division of Material Science Hans Weber. Operators
Basic Formulas 17-1-1 Division of Material Science Hans Weer The de Broglie wave length λ = h p The Schrödinger equation Hψr,t = i h t ψr,t Stationary states Hψr,t = Eψr,t Collection of formulae Quantum
More informationNosé-Hoover Thermostats
Nosé- Nosé- Texas A & M, UT Austin, St. Olaf College October 29, 2013 Physical Model Nosé- q Figure: Simple Oscillator p Physical Model Continued Heat Bath T p Nosé- q Figure: Simple Oscillator in a heat
More informationExponentially small splitting of separatrices of the pendulum: two different examples. Marcel Guardia, Carme Olivé, Tere M-Seara
Exponentially small splitting of separatrices of the pendulum: two different examples Marcel Guardia, Carme Olivé, Tere M-Seara 1 A fast periodic perturbation of the pendulum We consider a non-autonomous
More informationA stroboscopic averaging algorithm for highly oscillatory delay problems
IMA Journal of Numerical Analysis (2018) Page 1 of 21 doi:10.1093/imanum/drn000 A stroboscopic averaging algorithm for highly oscillatory delay problems J. M. SANZ-SERNA Departamento de Matemáticas, Universidad
More information2:2:1 Resonance in the Quasiperiodic Mathieu Equation
Nonlinear Dynamics 31: 367 374, 003. 003 Kluwer Academic Publishers. Printed in the Netherlands. ::1 Resonance in the Quasiperiodic Mathieu Equation RICHARD RAND Department of Theoretical and Applied Mechanics,
More information