High-order time-splitting methods for Schrödinger equations
|
|
- Loreen Reynolds
- 6 years ago
- Views:
Transcription
1 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
2 Nonlinear Schrödinger equations Bose Einstein condensation In our laboratories temperatures are measured in micro- or nanokelvin... In this ultracold world... atoms move at a snail s pace... and behave like matter waves. Interesting and fascinating new states of quantum matter are formed and investigated in our experiments. Grimm et al. Practical realisation. Observation of Bose Einstein condensation in physical experiments. Theoretical model. Mathematical description by systems of nonlinear Schrödinger equations. Numerical simulation. Favourable discretisations rely on time-splitting pseudospectral methods. See recent work by Bao, Cancès, Dion, Du, Jaksch, Markowich, Pérez-García, Shen, Tang, Tiwari, Shukla, Vazquez, Zhang etc.
3 Objectives Nonlinear Schrödinger equations Gross Pitaevskii equation (GPE). Nonlinear Schrödinger equation of the form i t ψ(x, t) = ( 2 2 m + U(x) + g ψ(x, t) 2) ψ(x, t). Numerical solution. Aims. Space discretisation based on pseudospectral methods. Time discretisation based on exponential operator splitting methods. Convergence analysis for linear Schrödinger equations. Comparison of high-order splitting methods regarding accuracy, efficiency and conservation of geometric properties (particle number, energy).
4 Outline Nonlinear Schrödinger equations 1 Nonlinear Schrödinger equations Gross Pitaevskii equation Abstract formulation 2 Pseudospectral methods Hermite vs. Fourier Numerical solution of nonlinearity 3 Exponential operator splitting Convergence analysis of linear Schrödinger equations Long-term integration of GPE
5 Gross Pitaevskii equation Gross Pitaevskii equation Abstract formulation Gross Pitaevskii equation (GPE). Normalised nonlinear Schrödinger equation i t ψ(x, t) = ( V (x) + ϑ ψ(x, t) 2) ψ(x, t). describes wave function of Bose Einstein condensate. Harmonic trap. Consider physically relevant case of harmonic potential V (x) = 1 2 V H(x) = 1 2 dj=1 γ 2 j x 2 j Geometric properties. Conservation of particle number ψ(, t) L t x 2 1 and energy E ( ψ(, t) ) ( ( = V ϑ ψ(, t) 2) ) ψ(, t) ψ(, t) L 2 3
6 Gross Pitaevskii equation Gross Pitaevskii equation Abstract formulation Gross Pitaevskii equation (GPE). Nonlinear Schrödinger equation i t ψ(x, t) = ( V (x) + ϑ ψ(x, t) 2) ψ(x, t). describes the wave function ψ : R d [0, ) C of a Bose Einstein condensate. Groundstate solution. Special solution of the form ψ(x, t) = e iµt φ(x) that minimizes energy functional. Minimisation approach. Compute ground state by direct minimisation of energy functional, see Bao and Tang (2003) and Caliari et al. (2009). Use ground state as reliable reference solution for time integration.
7 Evolutionary Schrödinger equations Gross Pitaevskii equation Abstract formulation Abstract formulation. Interpret time-dependent nonlinear Schrödinger equation as abstract evolution equation u (t) = ( A + B(u(t)) ) u(t). Numerical solution. Split right hand side into two parts Hermite Fourier A = i 1 ( ) 2 VH A = i 1 2 B(u) = i ϑ u 2 B(u) = i ( 1 2 V H + ϑ u 2 ) Numerical method relies on solutions of evolution equations v (t) = A v(t) v(0) = v 0 given w (t) = B(w(t)) w(t) w(0) = w 0 given
8 Hermite spectral method Pseudospectral methods Hermite vs. Fourier Numerical solution of nonlinearity Spectral decomposition of harmonic oscillator. Properties of the harmonic oscillator: Hermite functions (H µ ) form ON-basis of L 2 (R d ) v = v µ H µ, v µ = (v H µ ) L 2. µ N d Hermite functions are eigenfunctions of i A = 1 ( ) 2 + VH ( ) VH Hµ = λ µ H µ, λ µ = d j=1 γ j ( µ j). Exact solution of v (t) = A v(t), v(0) = v µ H µ is given by v(t) = µ N d e itλµ v µ H µ.
9 Hermite pseudospectral method Pseudospectral methods Hermite vs. Fourier Numerical solution of nonlinearity Numerical approximation. Truncate infinite sum and use the first M Hermite functions v(t) = v µ (t)h µ. µ M 1 Computation of spectral coefficients. Approximation of Hermite spectral coefficients by Gauss Hermite quadrature formula v µ = v(x) H µ (x) dx ω µ e ξ k 2 v(ξ k ) H µ (ξ k ). R d k Implementation. Hermite transform and inverse Hermite transform in 2D can be realised by two matrix matrix multiplications. Complexity of Hermite transform in 2D is O(M 3 ).
10 Fourier spectral method Pseudospectral methods Hermite vs. Fourier Numerical solution of nonlinearity Spectral decomposition of Laplacian. Properties of the Laplacian: Fourier basis functions (F µ ) form ON-basis of L 2 (( a, a) d ) v = µ Z d v µ F µ, v µ = (v F µ ) L 2. Fourier basis functions are eigenfunctions of i A = F µ = λ µ F µ, λ µ = 1 2 a 2 π 2 d µ 2 j. j=1 Exact solution of v (t) = A v(t), v(0) = v µ F µ is given by v(t) = µ Z d e itλµ v µ F µ.
11 Fourier pseudospectral method Pseudospectral methods Hermite vs. Fourier Numerical solution of nonlinearity Numerical approximation. Truncate infinite sum and use the first M Fourier functions v(t) = M 2 v µ (t)h µ. µ M 2 1 Computation of spectral coefficients. Approximation of Fourier spectral coefficients by trapezoidal rule with equidistant grid points v µ = v(x) F µ (x) dx ω v(ξ k ) F µ (ξ k ). ( a,a) d Implementation. Use Fast Fourier Transformation (FFT). Complexity of FFT in 2D is O(log(M) M 2 ).
12 Numerical experiment (CPU-time) Computation time of spectral methods. Pseudospectral methods Hermite vs. Fourier Numerical solution of nonlinearity 10 2 Hermite 1D Fourier 1D Hermite 2D Fourier 2D CPU seconds CPU seconds degree of freedom degree of freedom Figure: Computation time of the Hermite and Fourier spectral methods in one (left picture) and two (right picture) space dimensions using M = 2 i, 4 i 8, basis functions in each space direction.
13 Numerical experiment (spatial error) Pseudospectral methods Hermite vs. Fourier Numerical solution of nonlinearity error ϑ = 1 ϑ = 10 ϑ = 100 ϑ = degree of freedom error ϑ = 1 ϑ = 10 ϑ = 100 ϑ = degree of freedom Figure: Spatial error of the Hermite (left picture) and Fourier (right picture) spectral method for different values of the coupling constant ϑ.
14 Pseudospectral methods Hermite vs. Fourier Numerical solution of nonlinearity Numerical approximation of nonlinear part B Exact solution. Exact solution of the nonlinear part ( w (t) = i ɛv + ϑ w(t) 2) w(t), ɛ {0, 1}, is given by ( ( w(t) )(x) = e it ɛv (x)+ϑ w 0 (x) 2) w 0 (x). Numerical approximation. Evaluate pointwise multiplication only on the Gauss-Hermite quadrature nodes ξ k. Implementation. Numerical approximation in 2D can be realised by a pointwise matrix-matrix multiplication. Complexity in 2D is O(M 2 ).
15 Exponential operator splitting methods Problem class. Linear evolution equation Exponential operator splitting Convergence analysis of linear Schrödinger equations Long-term integration of GPE u (t) = A u(t) + B u(t), t 0, u(0) given. Compute numerical approximation u n u(t n ) at time t n = h n. Example method. Strang or symmetric Lie Trotter splitting u n+1 = e h 2 A e hb e h 2 A u n, see Trotter (1959) and Strang (1968). Method class. Higher-order exponential operator splitting methods of the form s u n+1 = e bjha e ajhb u n. j=1 with coefficients a j, b j (1 j s).
16 Exponential operator splitting Convergence analysis of linear Schrödinger equations Long-term integration of GPE Splitting methods for nonlinear Schrödinger equations Problem class. Nonlinear evolution equation of the form u (t) = A u(t) + B(u(t)) u(t), t 0, u(0) given. Example method. Strang or symmetric Lie Trotter splitting u n+1 = e h 2 A e hb(e h 2 A u n) e h 2 A u n. Method class. Higher-order exponential operator splitting methods of the form U n,1 = e a 1hA u n 1, U n,j = e a jha e b j 1hB(U n,j 1 ) U n,j 1, 2 j s, u n = e bshb(un,s) U n,s, with coefficients a j, b j (1 j s).
17 Exponential operator splitting Convergence analysis of linear Schrödinger equations Long-term integration of GPE High-order exponential operator splitting methods method order #Compositions McLachlan McLachlan p = 2 s = 3 Strang Strang p = 2 s = 2 BM4-1 Blanes & Moan PRKS 6 p = 4 s = 7 BM4-2 Blanes & Moan SRKN b 6 p = 4 s = 7 M4 McLachlan p = 4 s = 6 S4 Suzuki p = 4 s = 6 Y4 Yoshida p = 4 s = 4 BM6-1 Blanes & Moan PRKS 10 p = 6 s = 11 BM6-2 Blanes & Moan SRKN b 11 p = 6 s = 12 BM6-3 Blanes & Moan SRKN a 14 p = 6 s = 15 KL6 Kahan & Li p = 6 s = 10 S6 Suzuki p = 6 s = 26 Y6 Yoshida p = 6 s = 8 Table: Splitting methods of order p involving s compositions.
18 Exponential operator splitting Convergence analysis of linear Schrödinger equations Long-term integration of GPE Schrödinger equation involving unbounded potential Problem class. Linear Schrödinger equation of the form i t ψ(x, t) = ( V (x) + W (x)) ψ(x, t) involving unbounded polynomial potential W (x) = α µ x µ, α µ R. µ N d. Consider abstract evolution equation and split right hand side into A v(t) = i 1 2 ( V H(x)) v(t) and B w(t) = i W (x) w(t). Aim. Convergence result for exponential operator splitting Hermite pseudospectral methods of arbitrary order.
19 Exponential operator splitting Convergence analysis of linear Schrödinger equations Long-term integration of GPE Convergence result for linear Schrödinger equation Theorem (Neuhauser and Thalhammer (2009)) Suppose that the coefficients of the splitting method satisfy the classical order conditions for p 1. Then, provided that the exact solution u(t) is sufficiently regular, the following estimate holds u n u(t n ) L 2 C u 0 u(0) L 2 + C h p, 0 nh T Y KL S BM6 1 BM BM Strang McLachlan Y S M BM4 1 BM
20 Exponential operator splitting Convergence analysis of linear Schrödinger equations Long-term integration of GPE Numerical experiment (long-term integration) Numerical experiment (Caliari et al. (2009)). Illustrates the temporal order of Hermite and Fourier pseudospectral splitting methods for the GPE i t ψ(x, t) = ( V (x) + ψ(x, t) 2) ψ(x, t) involving the potential V (x) = 1 2 (x x 2 2). Choose the ground state solution of the GPE 0.7 involving the potential V (x) = x1 2 + x as initial value. Compute a reference solution at T = 400 ( basis functions, 2 17 time steps) t x
21 Exponential operator splitting Convergence analysis of linear Schrödinger equations Long-term integration of GPE Numerical experiment (long-term integration) Numerical experiment. Illustrates the temporal order of Hermite and Fourier pseudospectral splitting methods for the GPE i t ψ(x, t) = ( V (x) + ψ(x, t) 2) ψ(x, t). Compute numerical approximations by time-splitting spectral methods with 2 i, 6 i 8, basis functions in each space direction and 2 i, 6 i 15, time steps. Compute numerical approximations by standard explicit Runge Kutta methods of order four. Compute error in discrete L 2 -norm. Measure the particle number and energy conservation.
22 Exponential operator splitting Convergence analysis of linear Schrödinger equations Long-term integration of GPE Numerical experiment (long-term integration) tol. method d.o.f. #transf. pn E < 10 2 Hermite < 10 2 Fourier < 10 2 Hermite < 10 2 Fourier < 10 2 Hermite < 10 2 Fourier < 10 2 Hermite rk < 10 2 Fourier rk < 10 2 Hermite ode < 10 2 Fourier ode
23 Exponential operator splitting Convergence analysis of linear Schrödinger equations Long-term integration of GPE Numerical experiment (long-term integration) tol. method d.o.f. #transf. pn E < 10 4 Hermite < 10 4 Fourier < 10 4 Hermite < 10 4 Fourier < 10 4 Hermite rk < 10 4 Fourier rk < 10 4 Hermite ode < 10 4 Fourier ode
24 Observations Nonlinear Schrödinger equations Exponential operator splitting Convergence analysis of linear Schrödinger equations Long-term integration of GPE The number of required Hermite basis functions is smaller than the number of Fourier basis functions. The number of required Hermite transforms is (in general) smaller than the required Fourier transforms. The splitting methods outperform the explicit Runge Kutta methods. The splitting methods of order four and six are superior to the second order Strang splitting.
25 Conclusion and future work Exponential operator splitting Convergence analysis of linear Schrödinger equations Long-term integration of GPE Contents. High accuracy discretisations by time-splitting pseudospectral methods. Convergence analysis for linear evolutionary Schrödinger equations. Comparison of time-splitting methods regarding accuracy, efficiency, and geometric properties. Future work. Provide convergence analysis of high-order time-splitting pseudospectral methods for nonlinear Schrödinger equations. Study other basis functions, well adapted for a certain potential.
Three 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 informationRapporto di ricerca Research report
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report September 2009 76/2009 Spectral methods for dissipative nonlinear Schrödinger equations Laura M. Morato
More informationNumerical simulation of the dynamics of the rotating dipolar Bose-Einstein condensates
Numerical simulation of the dynamics of the rotating dipolar Bose-Einstein condensates Qinglin TANG Inria, University of Lorraine Joint work with: Prof. Weizhu BAO, Daniel MARAHRENS, Yanzhi ZHANG and Yong
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 informationQinglin TANG. INRIA & IECL, University of Lorraine, France
Numerical methods on simulating dynamics of the nonlinear Schrödinger equation with rotation and/or nonlocal interactions Qinglin TANG INRIA & IECL, University of Lorraine, France Collaborators: Xavier
More informationA simple Dufort-Frankel type scheme for the Gross-Pitaevskii equation of Bose-Einstein condensates on different geometries
A simple Dufort-Frankel type scheme for the Gross-Pitaevskii equation of Bose-Einstein condensates on different geometries Ming-Chih Lai Chung-Yin Huang Te-Sheng Lin Department of Applied Mathematics,
More informationBifurcations of Multi-Vortex Configurations in Rotating Bose Einstein Condensates
Milan J. Math. Vol. 85 (2017) 331 367 DOI 10.1007/s00032-017-0275-8 Published online November 17, 2017 2017 Springer International Publishing AG, part of Springer Nature Milan Journal of Mathematics Bifurcations
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 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 informationHigh-order actions and their applications to honor our friend and
High-order actions and their applications to honor our friend and collaborator Siu A. Chin Universitat Politecnica de Catalunya, Barcelona, Spain Lecture: Magnus Expansion and Suzuki s Method Outline of
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 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 informationA uniformly convergent numerical method for singularly perturbed nonlinear eigenvalue problems
A uniformly convergent numerical method for singularly perturbed nonlinear eigenvalue problems Weizhu Bao and Ming-Huang Chai Abstract In this paper we propose a uniformly convergent numerical method for
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 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 informationResearch Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
Applied Mathematics Volume 013, Article ID 591636, 5 pages http://dx.doi.org/10.1155/013/591636 Research Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationA Uniformly Convergent Numerical Method for Singularly Perturbed Nonlinear Eigenvalue Problems
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol 4, No, pp 35-6 Commun Comput Phys July 8 A Uniformly Convergent Numerical Method for Singularly Perturbed Nonlinear Eigenvalue Problems Weizhu Bao, and Ming-Huang
More informationREVIEW: The Matching Method Algorithm
Lecture 26: Numerov Algorithm for Solving the Time-Independent Schrödinger Equation 1 REVIEW: The Matching Method Algorithm Need for a more general method The shooting method for solving the time-independent
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 informationEnergy-preserving Pseudo-spectral Methods for Klein-Gordon-Schrödinger Equations
.. Energy-preserving Pseudo-spectral Methods for Klein-Gordon-Schrödinger Equations Jingjing Zhang Henan Polytechnic University December 9, 211 Outline.1 Background.2 Aim, Methodology And Motivation.3
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationPAPER 84 QUANTUM FLUIDS
MATHEMATICAL TRIPOS Part III Wednesday 6 June 2007 9.00 to 11.00 PAPER 84 QUANTUM FLUIDS Attempt TWO questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY REQUIREMENTS
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 informationLecture I: numerical computation of the dynamics of Schrödinger equations. Extension to the GPE for BECs
/67 Lecture I: numerical computation of the dynamics of Schrödinger equations. Extension to the GPE for BECs Xavier Antoine Institut Elie Cartan de Lorraine (IECL) Université de Lorraine INRIA-CORIDA Team
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More informationNumerical Methods for the Dynamics of the Nonlinear Schrodinger / Gross-Pitaevskii Equations. Weizhu Bao.
Numerical Methods for the Dynamics of the Nonlinear Schrodinger / Gross-Pitaevskii Equations Weizhu Bao Department of Mathematics National University of Singapore Email: matbaowz@nus.edu.sg URL: http://www.math.nus.edu.sg/~bao
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 informationSolitons and vortices in Bose-Einstein condensates with finite-range interaction
Solitons and vortices in Bose-Einstein condensates with finite-range interaction Luca Salasnich Dipartimento di Fisica e Astronomia Galileo Galilei and CNISM, Università di Padova INO-CNR, Research Unit
More informationc 2017 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 39, No. 5, pp. A198 A1950 c 017 Society for Industrial and Applied Mathematics HERMITE SPECTRAL METHODS FOR FRACTIONAL PDEs IN UNBOUNDED DOMAINS ZHIPING MAO AND JIE SHEN Abstract.
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationA stroboscopic averaging technique for highly-oscillatory Schrödinger equation
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,
More informationApplied Numerical Mathematics
Applied Numerical Mathematics 131 (2018) 1 15 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum The spectral-galerkin approximation of nonlinear eigenvalue
More informationNo-hair and uniqueness results for analogue black holes
No-hair and uniqueness results for analogue black holes LPT Orsay, France April 25, 2016 [FM, Renaud Parentani, and Robin Zegers, PRD93 065039] Outline Introduction 1 Introduction 2 3 Introduction Hawking
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
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 informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationContents. Xavier ANTOINE Romain DUBOSCQ. 1 Introduction 2
Robust and efficient preconditioned Krylov spectral solvers for computing the ground states and dynamics of fast rotating and strongly interacting Bose-Einstein condensates Xavier ANTOINE Romain DUBOSCQ
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 informationQuantized Vortex Stability and Dynamics in Superfluidity and Superconductivity
Quantized Vortex Stability and Dynamics in Superfluidity and Superconductivity Weizhu Bao Department of Mathematics National University of Singapore Email: matbaowz@nus.edu.sg URL: http://www.math.nus.edu.sg/~bao
More informationEnergy Preserving Numerical Integration Methods
Energy Preserving Numerical Integration Methods Department of Mathematics and Statistics, La Trobe University Supported by Australian Research Council Professorial Fellowship Collaborators: Celledoni (Trondheim)
More informationNumerical simulation of Bose-Einstein Condensates based on Gross-Pitaevskii Equations
1/38 Numerical simulation of Bose-Einstein Condensates based on Gross-Pitaevskii Equations Xavier ANTOINE Institut Elie Cartan de Lorraine (IECL) & Inria-SPHINX Team Université de Lorraine, France Funded
More informationLecture 2: Weak Interactions and BEC
Lecture 2: Weak Interactions and BEC Previous lecture: Ideal gas model gives a fair intuition for occurrence of BEC but is unphysical (infinite compressibility, shape of condensate...) Order parameter
More informationTHE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS
THE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS MARIANNA KHANAMIRYAN Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road,
More informationLecture 6 Quantum Mechanical Systems and Measurements
Lecture 6 Quantum Mechanical Systems and Measurements Today s Program: 1. Simple Harmonic Oscillator (SHO). Principle of spectral decomposition. 3. Predicting the results of measurements, fourth postulate
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationPolynomial Approximation: The Fourier System
Polynomial Approximation: The Fourier System Charles B. I. Chilaka CASA Seminar 17th October, 2007 Outline 1 Introduction and problem formulation 2 The continuous Fourier expansion 3 The discrete Fourier
More informationSoliton-like Solutions to NLS on Compact Manifolds
Soliton-like Solutions to NLS on Compact Manifolds H. Christianson joint works with J. Marzuola (UNC), and with P. Albin (Jussieu and UIUC), J. Marzuola, and L. Thomann (Nantes) Department of Mathematics
More informationASC Report No. 14/2017. W. Auzinger, I. Brezinova, H. Hofsta tter, O. Koch, M. Quell
ASC Report No. 14/2017 Practical splitting methods for the adaptive integration of nonlinear evolution equations. Part II: Comparison of local error estimation and step-selection strategies for nonlinear
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 informationLecture 5: Harmonic oscillator, Morse Oscillator, 1D Rigid Rotor
Lecture 5: Harmonic oscillator, Morse Oscillator, 1D Rigid Rotor It turns out that the boundary condition of the wavefunction going to zero at infinity is sufficient to quantize the value of energy that
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 informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationSTRANG SPLITTING METHODS FOR A QUASILINEAR SCHRÖDINGER EQUATION - CONVERGENCE, INSTABILITY AND DYNAMICS 1. INTRODUCTION.
STRANG SPLITTING METHODS FOR A QUASILINEAR SCHRÖDINGER EQUATION - CONVERGENCE, INSTABILITY AND DYNAMICS JIANFENG LU AND JEREMY L. MARZUOLA ABSTRACT. We study the Strang splitting scheme for quasilinear
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 informationQuasi-statically driven Otto engine operating simple few-body quantum systems in a 1D infinite square well
October 9, 017 Master thesis Quasi-statically driven Otto engine operating simple few-body quantum systems in a 1D infinite square well Division of Mathematical Physics, Department of Physics, Lund University
More informationMATH3383. Quantum Mechanics. Appendix D: Hermite Equation; Orthogonal Polynomials
MATH3383. Quantum Mechanics. Appendix D: Hermite Equation; Orthogonal Polynomials. Hermite Equation In the study of the eigenvalue problem of the Hamiltonian for the quantum harmonic oscillator we have
More informationNotes on Quantum Mechanics
Notes on Quantum Mechanics Kevin S. Huang Contents 1 The Wave Function 1 1.1 The Schrodinger Equation............................ 1 1. Probability.................................... 1.3 Normalization...................................
More informationGround State Patterns of Spin-1 Bose-Einstein condensation via Γ-convergence Theory
Ground State Patterns of Spin-1 Bose-Einstein condensation via Γ-convergence Theory Tien-Tsan Shieh joint work with I-Liang Chern and Chiu-Fen Chou National Center of Theoretical Science December 19, 2015
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 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 informationPhysics 137A Quantum Mechanics Fall 2012 Midterm II - Solutions
Physics 37A Quantum Mechanics Fall 0 Midterm II - Solutions These are the solutions to the exam given to Lecture Problem [5 points] Consider a particle with mass m charge q in a simple harmonic oscillator
More informationA Novel Approach with Time-Splitting Spectral Technique for the Coupled Schrödinger Boussinesq Equations Involving Riesz Fractional Derivative
Commun. Theor. Phys. 68 2017 301 308 Vol. 68, No. 3, September 1, 2017 A Novel Approach with Time-Splitting Spectral Technique for the Coupled Schrödinger Boussinesq Equations Involving Riesz Fractional
More informationVortices in Bose-Einstein condensates. Ionut Danaila
Vortices in Bose-Einstein condensates 3D numerical simulations Ionut Danaila Laboratoire Jacques Louis Lions Université Pierre et Marie Curie (Paris 6) http://www.ann.jussieu.fr/ danaila October 16, 2008
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
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 informationFRACTIONAL ERROR ESTIMATES OF SPLITTING SCHEMES FOR THE NONLINEAR SCHRÖDINGER EQUATION
FRACTIONAL ERROR ESTIMATES OF SPLITTING SCHEMES FOR THE NONLINEAR SCHRÖDINGER EQUATION JOHANNES EILINGHOFF, ROLAND SCHNAUBELT, AND KATHARINA SCHRATZ Abstract. We investigate the Lie and the Strang splitting
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationMATHEMATICAL THEORY AND NUMERICAL METHODS FOR GROSS-PITAEVSKII EQUATIONS AND APPLICATIONS CAI YONGYONG
MATHEMATICAL THEORY AND NUMERICAL METHODS FOR GROSS-PITAEVSKII EQUATIONS AND APPLICATIONS CAI YONGYONG (M.Sc., Peking University) A THESIS SUBMITTED FOR THE DEGREE OF DOCTOR OF PHILOSOPHY DEPARTMENT OF
More informationCOMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS
COMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS MIKKO SALO AND JENN-NAN WANG Abstract. This work is motivated by the inverse conductivity problem of identifying an embedded object in
More informationSplitting and composition methods for the time dependent Schrödinger equation
Splitting and composition methods for the time dependent Schrödinger equation S. Blanes, joint work with F. Casas and A. Murua Instituto de Matemática Multidisciplinar Universidad Politécnica de Valencia,
More informationStability and instability of solitons in inhomogeneous media
Stability and instability of solitons in inhomogeneous media Yonatan Sivan, Tel Aviv University, Israel now at Purdue University, USA G. Fibich, Tel Aviv University, Israel M. Weinstein, Columbia University,
More informationNumerical Study of Oscillatory Regimes in the KP equation
Numerical Study of Oscillatory Regimes in the KP equation C. Klein, MPI for Mathematics in the Sciences, Leipzig, with C. Sparber, P. Markowich, Vienna, math-ph/"#"$"%& C. Sparber (generalized KP), personal-homepages.mis.mpg.de/klein/
More informationhttp://www.springer.com/3-540-30725-7 Erratum Spectral Methods Fundamentals in Single Domains C. Canuto M.Y. Hussaini A. Quarteroni T.A. Zang Springer-Verlag Berlin Heidelberg 2006 Due to a technical error
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationA butterfly algorithm. for the geometric nonuniform FFT. John Strain Mathematics Department UC Berkeley February 2013
A butterfly algorithm for the geometric nonuniform FFT John Strain Mathematics Department UC Berkeley February 2013 1 FAST FOURIER TRANSFORMS Standard pointwise FFT evaluates f(k) = n j=1 e 2πikj/n f j
More informationLecture 4.6: Some special orthogonal functions
Lecture 4.6: Some special orthogonal functions Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4340, Advanced Engineering Mathematics
More informationAppendix C: Recapitulation of Numerical schemes
Appendix C: Recapitulation of Numerical schemes August 31, 2009) SUMMARY: Certain numerical schemes of general use are regrouped here in order to facilitate implementations of simple models C1 The tridiagonal
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationSplitting Methods for Non-autonomous Hamiltonian Equations
Journal of Computational Physics 170, 205 230 2001 doi:10.1006/jcph.2001.6733, available online at http://www.idealibrary.com on Splitting Methods for Non-autonomous Hamiltonian Equations S. Blanes and
More informationSymplectic time-average propagators for the Schödinger equation with a time-dependent Hamiltonian
Symplectic time-average propagators for the Schödinger equation Symplectic time-average propagators for the Schödinger equation with a time-dependent Hamiltonian Sergio Blanes,, a Fernando Casas, 2, b
More informationExponential integrators for semilinear parabolic problems
Exponential integrators for semilinear parabolic problems Marlis Hochbruck Heinrich-Heine University Düsseldorf Germany Innsbruck, October 2004 p. Outline Exponential integrators general class of methods
More informationNumerical methods for computing vortex states in rotating Bose-Einstein condensates. Ionut Danaila
Numerical methods for computing vortex states in rotating Bose-Einstein condensates Ionut Danaila Laboratoire de mathématiques Raphaël Salem Université de Rouen www.univ-rouen.fr/lmrs/persopage/danaila
More information1 Distributions (due January 22, 2009)
Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =
More informationNUMERICAL ANALYSIS 2 - FINAL EXAM Summer Term 2006 Matrikelnummer:
Prof. Dr. O. Junge, T. März Scientific Computing - M3 Center for Mathematical Sciences Munich University of Technology Name: NUMERICAL ANALYSIS 2 - FINAL EXAM Summer Term 2006 Matrikelnummer: I agree to
More informationA UNIFORMLY CONVERGENT NUMERICAL SCHEME FOR SINGULARLY PERTURBED NON-LINEAR EIGENVALUE PROBLEM UNDER CONSTRAINTS CHAI MING HUANG
A UNIFORMLY CONVERGENT NUMERICAL SCHEME FOR SINGULARLY PERTURBED NON-LINEAR EIGENVALUE PROBLEM UNDER CONSTRAINTS CHAI MING HUANG (B.Sc.(Hons.), NUS) A THESIS SUBMITTED FOR THE DEGREE OF MASTER OF SCIENCE
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationOne-dimensional Schrödinger equation
Chapter 1 One-dimensional Schrödinger equation In this chapter we will start from the harmonic oscillator to introduce a general numerical methodology to solve the one-dimensional, time-independent Schrödinger
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
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 informationNonstandard Finite Difference Time Domain (NSFDTD) Method for Solving the Schrödinger Equation
Pramana J. Phys. (2018)?: #### DOI??/?? c Indian Academy of Sciences Nonstandard Finite Difference Time Domain (NSFDTD) Method for Solving the Schrödinger Equation I WAYAN SUDIARTA mitru et al. [11] have
More informationComputational Methods for Optoelectronics
IPT 526 Computational Methods for Optoelectronics Ray-Kuang Lee Department of Electrical Engineering and Institute of Photonics Technologies National Tsing-Hua University, Hsinchu, Taiwan e-mail: rklee@ee.nthu.edu.tw
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationChemistry 432 Problem Set 4 Spring 2018 Solutions
Chemistry 4 Problem Set 4 Spring 18 Solutions 1. V I II III a b c A one-dimensional particle of mass m is confined to move under the influence of the potential x a V V (x) = a < x b b x c elsewhere and
More informationExamination paper for TMA4125 Matematikk 4N
Department of Mathematical Sciences Examination paper for TMA45 Matematikk 4N Academic contact during examination: Anne Kværnø a, Louis-Philippe Thibault b Phone: a 9 66 38 4, b 9 3 0 95 Examination date:
More informationMatter-Wave Soliton Molecules
Matter-Wave Soliton Molecules Usama Al Khawaja UAE University 6 Jan. 01 First International Winter School on Quantum Gases Algiers, January 1-31, 01 Outline Two solitons exact solution: new form Center-of-mass
More informationA Vector-Space Approach for Stochastic Finite Element Analysis
A Vector-Space Approach for Stochastic Finite Element Analysis S Adhikari 1 1 Swansea University, UK CST2010: Valencia, Spain Adhikari (Swansea) Vector-Space Approach for SFEM 14-17 September, 2010 1 /
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More information