Numerical aspects of the nonlinear Schro dinger equation. Christophe Besse
|
|
- Madeleine Hicks
- 5 years ago
- Views:
Transcription
1 Numerical aspects of the nonlinear Schro dinger equation by Christophe Besse Laboratoire Paul Painleve, Universite Lille 1, CNRS, Inria Simpaf Team Lille
2 Outline of the talk 1 Motivation 2 Numerical methods for NLS 3 Absorbing boundary conditions for the Schrödinger equations 4 Some numerical experiments.
3 Motivation Nonlinear Schrödinger Equation (NLS) iε tψ = ε2 2 ψ + V (x)ψ + β ψ 2 ψ, (x, t) R d x [0; T ], u(x, 0) = ψ 0 (x), x R d x. ψ(x, t) : complex-valued wave function V (x) : real-valued external potential β: interaction constant (= 0: linear, +1: repulsive interaction, 1: attractive interaction) ε: scaled Planck constant Numerical approximations ε = 1: standard; 0 < ε 1 and β = ±1 semi-classical regimes efficient, convergent and accurate numerical schemes compute efficient boundary conditions to deal with bounded computational domains.
4 Application of NLS In quantum physics Interaction between particles with quantum effect Bose-Einstein condensation (BEC): bosons at low temperature Superfluidity In plasma physics wave interaction between electron and ion In semiconductor industry In quantum chemistry chemical interaction based on the first principle In nonlinear optics & atom laser In fluid mechanics (water waves) NLS is the prototype for many dispersive systems Examples : Davey-Stewartson systems (water waves) { i tu + λ 2 x u + 2 y u = ν u 2 u + u xψ, α 2 x ψ + 2 y ψ = χ x( u 2 ).
5 Properties of NLS Time reversible Time transverse invariant (gauge invariant) V (x) V (x) + α ψ ψe itα/ε ψ unchanged Mass (wave energy) conservation N ψ (t) := ψ(x, t) 2 dx N ψ (0) := R d ψ(x, 0) 2 dx := R d ψ 0 (x) 2 dx, R d t 0 Energy ( or Hamiltonian) conservation ε 2 E ψ (t) := R d 2 ψ(x, t) 2 + V (x) ψ(x, t) 2 + β 2 ψ(x, t) 4 dx E ψ (0), t 0 Dispersion relation without external potential ψ(x, t) = ae i(k x ωt/ε) (plane wave solution) ω = ε2 2 k 2 + β a 2
6 Numerical difficulties Explicit vs Implicit (or computational cost) Spatial/temporal accuracy Stability / Convergence Keep the properties of NLS in the discretized level Time reversible & time transverse invariant Mass & energy conservation Dispersion conservation Resolution in the semiclassical regime: 0 < ε 1 Take BC into account
7 Outline of the talk 1 Motivation 2 Numerical methods for NLS 3 Absorbing boundary conditions for the Schrödinger equations 4 Some numerical experiments.
8 Numerical methods for NLS Classical schemes : Crank-Nicolson type : Delfour-Fortin-Payre, Dúran-Sanz Serna,... Runge-Kutta type : Akrivis-Dougalis-Karakashian,... etc... High order schemes, convergent, preserving of invariants, BUT semi-implicit : need of a nonlinear step (Newton) Example : Crank-Nicolson scheme (CNFD) iε ψn+1 ψ n = ( ε22 δt + β ψn ψ n 2 ) ψ + V (, t n+1/2 n+1 + ψ n ) 2 2 ψ 0 (x) = ψ 0 (x), x R d. = 0, where u n is the approximate value of u at time t n = nδt. Example : Dúran-Sanz Serna scheme (SSFD) ( iε ψn+1 ψ n = ε2 δt 2 + β ψ n+1 + ψ n ) V (, t n+1/2 ψ n+1 + ψ n ) 2 ψ 0 (x) = ψ 0 (x), x R d. = 0,
9 Numerical methods for NLS Classical schemes : Crank-Nicolson type : Delfour-Fortin-Payre, Dúran-Sanz Serna,... Runge-Kutta type : Akrivis-Dougalis-Karakashian,... etc... High order schemes, convergent, preserving of invariants, BUT semi-implicit : need of a nonlinear step (Newton) Example : Crank-Nicolson scheme (CNFD) iε ψn+1 ψ n = ( ε22 δt + β ψn ψ n 2 ) ψ + V (, t n+1/2 n+1 + ψ n ) 2 2 ψ 0 (x) = ψ 0 (x), x R d. = 0, where u n is the approximate value of u at time t n = nδt. Example : Dúran-Sanz Serna scheme (SSFD) ( iε ψn+1 ψ n = ε2 δt 2 + β ψ n+1 + ψ n ) V (, t n+1/2 ψ n+1 + ψ n ) 2 ψ 0 (x) = ψ 0 (x), x R d. = 0, Possible cures : to decouple linear and nonlinear parts of (NLS)
10 Numerical methods for NLS Relaxation scheme 1 Relaxation scheme: to regard (NLS) as a Schrödinger-Poisson system where the Poisson equation for the potential Υ is replaced by the explicit formula Υ = ψ 2, that is, { Υ = ψ 2, (x, t) R d R +, iε tψ = ε2 2 ψ + V (x)ψ + βψυ, (x, t) Rd R +,
11 Numerical methods for NLS Relaxation scheme 1 Relaxation scheme: to regard (NLS) as a Schrödinger-Poisson system where the Poisson equation for the potential Υ is replaced by the explicit formula Υ = ψ 2, that is, { Υ = ψ 2, (x, t) R d R +, iε tψ = ε2 2 ψ + V (x)ψ + βψυ, (x, t) Rd R +, Relaxation scheme (ReFD) CB 04 Υ n+1/2 + Υ n 1/2 = ψ n 2, x R d, 2 iε ψn+1 ψ n δt = ( ε22 ) ψ n+1 + V (x) + + ψ n βυn+1/2 = 0, x R d, 2
12 Numerical methods for NLS Relaxation scheme 1 Relaxation scheme: to regard (NLS) as a Schrödinger-Poisson system where the Poisson equation for the potential Υ is replaced by the explicit formula Υ = ψ 2, that is, { Υ = ψ 2, (x, t) R d R +, iε tψ = ε2 2 ψ + V (x)ψ + βψυ, (x, t) Rd R +, Relaxation scheme (ReFD) CB 04 Υ n+1/2 + Υ n 1/2 = ψ n 2, x R d, 2 iε ψn+1 ψ n δt = ( ε22 ) ψ n+1 + V (x) + + ψ n βυn+1/2 = 0, x R d, 2 Advantages convergent scheme, efficient, preserves the invariants, easily adaptable Drawbacks Υ u 2 at the discrete level but Υ(x, t) = t 0 s u 2 (x, s)ds make proof of convergence harder, second order, Υ n > 0 n?.
13 Numerical methods for NLS Splitting scheme 2 Splitting scheme : Splitting methods are based on a decomposition of the flow of (NLS). We introduce S t the flow of (NLS) : ψ(t, ) = S t ψ 0 and the two flows X t and Y t solutions to tv = i ε 2 v : v(, y) = Xt v(, 0) iε tw = V (x)w + β w 2 w = 0 : w(, y) = Y t w(, 0) Solve nonlinear ODE analytically since t( w(x, t) 2 ) = 0 w(x, t) = w(x, 0). w(x, t) = e it[v (x)+β w(x,0) 2 ]/ε w(x, 0). Idea : to approximate the flow S t by combining the two flows X t and Y t Order 1 (Lie) Z t L = Xt Y t or Z t L = Y t X t Order 2 (Strang) Z t S = Xt/2 Y t X t/2 or Z t S = Y t/2 X t Y t/2 Higher order See Descombes, Thalhammer
14 Numerical methods for NLS Splitting scheme Convergence : Theorem (CB, Bidégary, Descombes 02) u 0 H 2, T > 0, C and h 0 such that h (0, h 0 ], n such that nh T ( ) n ZL h u0 S nh u 0 Ch u0 H 2. ( ) Moreover, if u 0 H 4 n, then ZS h u0 S nh u 0 Ch 2 u 0 H 4.
15 Numerical methods for NLS Splitting scheme Convergence : Theorem (CB, Bidégary, Descombes 02) u 0 H 2, T > 0, C and h 0 such that h (0, h 0 ], n such that nh T ( ) n ZL h u0 S nh u 0 Ch u0 H 2. ( ) Moreover, if u 0 H 4 n, then ZS h u0 S nh u 0 Ch 2 u 0 H 4. Advantages convergent scheme, efficient, easily adaptable, very well adapted in semi-classical regimes Drawbacks energy is not preserved, periodic boundary conditions if Fourier for the space variable.
16 Numerical methods for NLS Method TSSP CNFD SSFD ReFD Time Resersible Inv. by Gauge Change Mass Conservation Energy Conservation Dispersion conservation Unconditional L 2 Stability Convergence Explicit Scheme Time accuracy 2 th or 4 th 2 th 2 th 2 th
17 Outline of the talk 1 Motivation 2 Numerical methods for NLS 3 Absorbing boundary conditions for the Schrödinger equations 4 Some numerical experiments.
18 ABCs for Schrödinger equations NLS in R d i tψ + ψ + V (x, t, ψ) ψ = 0, (x, t) R d [0; T ] (S) lim ψ(x, t) = 0, t [0; T ] x + ψ(x, 0) = ψ 0 (x), x R d where V = V (x) + f(ψ). Motivation : the numerical resolution of Schrödinger equations leads to introduce a finite computational domain Ω. = introduction of a regular fictive boundary Γ Natural questions : 1 which are the satisfactory conditions to set on Γ [0, T ] such that the approximate solution coincides with the restriction to Ω of the true solution? 2 If these conditions exist, how do we handle them numerically? 3 Does the energy of the approximate solution decay in Ω (stability of the numerical schemes)?
19 ABCs for Schrödinger equations NLS in R d i tψ + ψ + V (x, t, ψ) ψ = 0, (x, t) R d [0; T ] (S) lim ψ(x, t) = 0, t [0; T ] x + ψ(x, 0) = ψ 0 (x), x R d where V = V (x) + f(ψ). Motivation : the numerical resolution of Schrödinger equations leads to introduce a finite computational domain Ω. = introduction of a regular fictive boundary Γ Natural questions : 1 which are the satisfactory conditions to set on Γ [0, T ] such that the approximate solution coincides with the restriction to Ω of the true solution? 2 If these conditions exist, how do we handle them numerically? 3 Does the energy of the approximate solution decay in Ω (stability of the numerical schemes)?
20 Domain truncation Problem : Mesh an unbounded domain (here in 1D) R d [0; T ] t ψ 0 x
21 Domain truncation Problem : Mesh an unbounded domain (here in 1D) R d [0; T ] t T Σ T Ω T ψ x l 0 x r x Truncation R [0; T ] Ω T := ]x l, x r[ [0; T ] Introduction of a fictitious boundary Σ T := Σ [0; T ] with Σ := Ω = {x l, x r} BC on the boundaryσ T
22 Domain truncation Problem : Mesh an unbounded domain (here in 1D) R d [0; T ] t T Σ T Ω T f(ψ, n ψ) = 0 ψ x l 0 x r x Truncation R [0; T ] Ω T := ]x l, x r[ [0; T ] Introduction of a fictitious boundary Σ T := Σ [0; T ] with Σ := Ω = {x l, x r} BC on the boundaryσ T The boundary condition on x l must represent the effect of the potential on ], x l ]. Expression of the boundary condtion with the help of the Dirichlet-to-Neumann map: nψ + iλ + ψ = 0, on Σ T
23 What happen if we do not take BC into account t Potential V (x) = x Initial datum: gaussian function ψ 0 (x) = e x2 +10ix Domain: Ω = [ 5; 15] u(x,t) x 2 Exact Solution Evolution of ψ w.r.t time t = t x x
24 What happen if we do not take BC into account t Potential V (x) = x Initial datum: gaussian function ψ 0 (x) = e x2 +10ix Domain: Ω = [ 5; 15] u(x,t) x Evolution of ψ w.r.t time t = 0.2 Approximated numerical solution Homogeneous Dirichlet Boundary Conditions: ψ Σ = 0 Parasistic reflexion t x x
25 The 1D Eq. without potential L L 1 Splitting between interior and exterior problems Interior Problem (i t + 2 x )v = 0, x Ω, t > 0, xv = xw, x Σ, t > 0, v(x, 0) = ψ 0(x), x Ω. Exterior problem (i t + 2 x )w = 0, x Ω, t > 0, w(x, t) = v(x, t), x = x l,r, t > 0, lim w(x, t) = 0, t > 0, x + w(x, 0) = 0, x Ω. left exterior problem interior problem right exterior problem output: Neumann data wx(x,t) input: Dirichlet data v(x,t) v (x,t) xl xr 2 Laplace transform w.r.t time (t τ) on Ω r: ODE in x 3 Argument: w L 2 (Ω r) to select the outgoing wave 4 Inverse Laplace transform 5 Exact boundary condition on x r: xv(x, t) x=xr = e iπ/4 1/2 t v(x r, t) TBC on Σ : where 1/2 t f(x, t) = 1 π t t 0 f(x, s) x=xr t s ds nv + e iπ/4 1/2 t v = 0, on Σ [0; T ]. fractional derivative operator of order 1/2
26 The 1D Eq. without potential The problem (S) is thus transformed in (S app) The Schrödinger Eq. in Ω i tψ + ψ = 0, (x, t) Ω T (S app) nψ + e iπ/4 1/2 t ψ = 0, on Σ T, ψ(x, 0) = ψ 0 (x), x Ω T Remark: 2 x + i t = ( n + i i t)( n i i t)
27 Extension to simple cases (1) V = 0 Expression of the exact Dirichlet-to-Neumann operator nψ + e iπ/4 1/2 t ψ = 0, on Σ T. (ABC 0 ) V = V l,r constant at the exterior of Ω ( ) nψ + e iπ/4 e itv l,r 1/2 t e itv l,r ψ = 0, on Σ T. V = V (t) : Gauge change t We set: v(x, t) = ψ(x, t)e iv(t) with V(t) = V (s) ds. 0 Then i tψ = (i tv V (t) v) e iv(t), where v is solution to the equation without potential Exact condition ABC 0 for v : nv + e iπ/4 1/2 t v = 0, on Σ T ( ) We come back to ψ : nψ + e iπ/4 e iv(t) 1/2 t e iv(t) ψ = 0, on Σ T
28 Non constant potentials If V = V (x, t), then the Laplace transform can not be used anymore. One have to introduce a new tool: pseudodifferential calculus.
29 Non constant potentials If V = V (x, t), then the Laplace transform can not be used anymore. One have to introduce a new tool: pseudodifferential calculus. Pseudodifferential Operators in 1D A pseudodifferential operator P (x, t, t) is described by its total symbol p(x, t, τ) in the Fourier space (τ is the covariable of t) P (x, t, t) u(x, t) = Ft 1 ( ) p(x, t, τ) û(x, τ) = p(x, t, τ) û(x, τ)e itτ dτ R Notations: P = Op(p), p(x, t, τ) = σ(p (x, t, t))
30 Non constant potentials Examples The fractional operators 1/2 t and I α/2 t 1/2 t f(t) = 1 t f(s) t ds π 0 t s I α/2 1 t t f(t) = (t s) α/2 1 f(s) ds Γ(α/2) 0 Nonlocal w.r.t time convolution operator Operator t 1/2 t I 1/2 t I t Symbol iτ e iπ/4 τ e iπ/4 τ 1 iτ Class OP S 1 OP S 1/2 OP S 1/2 OP S 1
31 Non constant potentials Case V = 0 TBC: nψ + e iπ/4 1/2 t ψ = 0, on Σ T. Case constant V = V TBC: nψ + e iπ/4 e itv 1/2 ( t e itv ψ ) = 0, on Σ T.
32 Non constant potentials Case V = 0 TBC: nψ + e iπ/4 1/2 t ψ = 0, on Σ T. nψ i Op ( τ ) ψ = 0, on Σ T. Case constant V = V TBC: nψ i e itv Op ( τ ) ( e itv ψ ) = 0, on Σ T.
33 Non constant potentials Case V = 0 TBC: nψ + e iπ/4 1/2 t ψ = 0, on Σ T. nψ i Op ( τ ) ψ = 0, on Σ T. Case constant V = V TBC: nψ i e itv Op ( τ ) ( e itv ψ ) = 0, on Σ T. ( τ ) nψ i Op + V (ψ) = 0, on Σ T. Lemma If a is a symbol belonging to S m independent of t, and V = V (x), then ( ) Op (a(τ V (x))) ψ = e itv (x) Op (a(τ)) e itv (x) ψ
34 Non constant potentials Case V = 0 TBC: nψ + e iπ/4 1/2 t ψ = 0, on Σ T. nψ i Op ( τ ) ψ = 0, on Σ T. Case constant V = V TBC: nψ i e itv Op ( τ ) ( e itv ψ ) = 0, on Σ T. ( τ ) nψ i Op + V (ψ) = 0, on Σ T. Lemma If a is a symbol belonging to S m independent of t, and V = V (x), then ( ) Op (a(τ V (x))) ψ = e itv (x) Op (a(τ)) e itv (x) ψ Case V = V (t) : Gauge change Antoine, Besse et Descombes, 2006 nψ i e iv(t) Op ( τ ) ( ) e iv(t) ψ = 0, on Σ T.
35 With potential V = V (x, t) + f( ψ ): two strategies Use of pseudodifferential calculus No More Exact! : Artifical Boundary Condition 1) Gauge change t v(x, t) = e iv(x,t) ψ(x, t), with V(x, t) = V (x, s, ψ) ds. 0 Involve operators e iv(x,t) Op ( τ ) ( ) e iv(x,t) ψ 2) Direct method No Gauge change Involve operators ( ) Op τ + V (x, t, ψ) (ψ) Equivalent strategies for V = V (x), non equivalent for V = V (x, t) In both cases, we still have to approximate Boundary Conditions, at different order M. Antoine, Besse et Klein, J. Comput. Phys., 2009
36 With potential V = V (x, t) + f( ψ ): two strategies Nonlinearities and general repulsive potentials x xv > 0 for x Ω V = f(x, ψ) V = f(x, ψ) and V(x, t) = t 0 f(x, ψ(x, s))ds Absorbing boundary conditions (ABC) for M = 4: ( ) ABC1 4 : nψ + e iπ/4 e iv 1/2 t e iv ψ i nv ( ) 4 eiv I t e iv ψ = 0 ABC 4 2 : nψ i i t + V ψ nv (i t + V ) 1 ψ = 0
37 Extension to higher dimensions Take into account the geometry: convex set with general boundary, smooth, with curvature κ. Ω n Generalized coordinates system of the boundary with respect to normal variable r and curvilinear abscissa s = 2 r + κ r r + h 1 s ( h 1 s ) κ r = h 1 κ : curvature of a parallel surface Σ r to Σ h(r, s) = 1 + rκ L = 2 r + κ r r + i t + h 1 s ( h 1 s ) + V τ s M (s) r n M (r, s) Schrödinger equation with variable coefficients: pseudo-differential calculus t t dual variable τ x r y s dual variable ξ Σ Σ r
38 Nonlocality Nonlocal both in space and time. Localizing in space: two approaches, valid for both strategies Taylor approach: Taylor expansion of the symbols for τ ξ 2 Thereby: τ ξ 2 + V = τ (1 + ξ2 τ V ) τ τ ξ 2 + V τ (1 + ξ2 2τ V ) = τ ξ2 2τ 2 = Localizing in space only 1 τ + V 2 1 τ Padé approximation approach: Op ( τ ξ 2 + V ) i t + Σ + V mod OP S 1 formal approximation of by Padé approximants = Localizing both in space AND time
39 Conclusion Two possible approaches for each strategy, so 4 families of ABC. Taylor approach Padé approach Gauge change ABC M 1,T ABC M 1,P Direct method ABC M 2,T ABC M 2,P ( 1/2 t, I 1/2 t, I t Op τ ξ 2) ( ) v = e iv u Op τ ξ 2 + V or
40 ABC: Taylor approach Gauge change ABC 2 1,T nψ + e iπ/4 e iv 1/2 t ( ) e iv ψ + κ 2 ψ ( κ ABC 3 1,T e iπ/4 e iv Σ 2 + i sv s + 1 ) 2 (i 2 s V ( sv)2 ) ) ( ) I t e iv ψ ( s(κ s) ABC 4 1,T + ie iv + κ3 + s 2κ + i sκ sv i sg( nv ) ( nv e iv I t nv e ψ) iv = 0 4 Direct method ABC 2 2,T nψ + e iπ/4 1/2 t ψ + κ 2 ψ ( κ ABC 3 2,T e iπ/ Σ 2 ( s(κ s) ABC 4 2,T + i 2 ) I 1/2 + κ3 + 2 s κ 8 iπ/4 sg(v ) t ψ e 2 ) I tψ i sg( nv ) 4 1/2 ( ) V I t V ψ ( ) I 1/2 t e iv ψ nv I t ( nv ψ) = 0
41 ABC: Padé approach Gauge change ABC 1 1,P nψ ie ) iv i t + Σ (e iv ψ ABC 2 1,P + κ ( ) 2 ψ + sveiv s (i t + Σ ) 1/2 e iv ψ κ ) 2 eiv (i t + Σ ) 1 Σ (e iv ψ = 0 Direct method ABC 1 2,P nψ i i t + Σ + V ψ ABC 2 2,P + κ 2 ψ κ 2 (i t + Σ + V ) 1 Σ ψ = 0
42 Taylor approach conditions Approximations of 1/2 t, I 1/2 t, I t by discrete convolutions, linked to the Crank-Nicolson scheme trapezoidal rule [Schmidt - Yevick (97), Antoine - Besse (03)] 1/2 2 t f(t n ) t I 1/2 t t f(t n ) 2 I tf(t n ) t 2 n β n k f k k=0 n α n k f k k=0 n γ n k f k k=0 (α 0, α 1, α 2,...) = (1, 1, 12, 12, 38, 38 ),... β k = ( 1) k α k, k 0 (γ 0, γ 1, γ 2,...) = (1, 2, 2, 2,...)
43 Padé approach conditions Approximation of i t + Σ + V Rational approximation of the square root by Padé approximants [Bruneau - Di Menza (95), Szeftel (04)] m z Rm(z) = a m k k=0 m a m k dm k z + d m k=1 k In the ABC M 2,P conditions: i t + Σ + V R m (i t + Σ + V ) i t + Σ + V ψ R m ( m ) m a m k ψ a m k dm k (i t + Σ + V + d m k ) 1 ψ }{{} k=0 k=1 ϕ k Auxiliary functions ϕ k, solutions of a Schrödinger equation on Σ T.
44 Outline of the talk 1 Motivation 2 Numerical methods for NLS 3 Absorbing boundary conditions for the Schrödinger equations 4 Some numerical experiments.
45 Numerical application ( i ix 2 1d linear case u(x, t) = 4t + i exp k 0 x + k0 2t ) k 0 = 8, Ω = [ 5, 5], 4t + i N = 1024, t = t=0 t=0.15 t=0.30 t=0.45 Evolution of u u x
46 2D linear case Explicit solution (2D) u(x 1, x 2, t) = i ( i 4t exp i x2 1 + x ix ) it. i 4t Finite Element Approximation (P 1 ): Ω i = D(0, 10), 3278 triangles, t = t=0.25 t=0.35 t=0.50
47
48 2D linear and nonlinear case with potentials Linear, Circular domain, V (x, y) = 5(x 2 + y 2 ) Linear, Mediator shaped domain, V (x, y) = 5 x 2 + y 2 Nonlinear, Circular domain, V (x, y) = 0 Nonlinear, Circular domain, V (x, y) = x 2 + y 2
Artificial boundary conditions for dispersive equations. Christophe Besse
Artificial boundary conditions for dispersive equations by Christophe Besse Institut Mathématique de Toulouse, Université Toulouse 3, CNRS Groupe de travail MathOcéan Bordeaux INSTITUT de MATHEMATIQUES
More informationAbsorbing boundary conditions for Schrödinger equations. Christophe Besse 1
Absorbing boundary conditions for Schrödinger equations by Christophe Besse 1 1 Laboratoire Paul Painlevé, Université Lille 1, Équipe INRIA Simpaf Workshop Stochastic Schrödinger Equation Paris Collaborators:
More informationOpen boundary conditions and computational schemes for Schrödinger equations with general potentials and nonlinearities
1 Open boundary conditions and computational schemes for Schrödinger equations with general potentials and nonlinearities Xavier Antoine, Pauline Klein Institut Elie Cartan Nancy, Nancy-Université, CNRS
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 informationComputer Physics Communications. Computational methods for the dynamics of the nonlinear Schrödinger/Gross-Pitaevskii equations
Computer Physics Communications (13) 1 3 Computer Physics Communications Computational methods for the dynamics of the nonlinear Schrödinger/Gross-Pitaevskii equations Xavier Antoine a,b, Weizhu Bao 1c,d,
More informationAbsorbing boundary conditions for the two-dimensional Schrödinger equation with an exterior potential. Part I: construction and a priori estimates
Absorbing boundary conditions for the two-dimensional Schrödinger equation with an exterior potential. Part I: construction and a priori estimates Xavier ANTOINE Institut Elie Cartan Nancy, Nancy-Université,
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 information1. Introduction. Let us consider the cubic nonlinear Schrödinger equation u t + i u + iε u 2 u =0, x R 2,t>0,
SIAM J. NUMER. ANAL. Vol. 4, No. 1, pp. 26 4 c 22 Society for Industrial and Applied Mathematics ORDER ESTIMATES IN TIME OF SPLITTING METHODS FOR THE NONLINEAR SCHRÖDINGER EQUATION CHRISTOPHE BESSE, BRIGITTE
More informationA Friendly Review of Absorbing Boundary Conditions and Perfectly Matched Layers for Classical and Relativistic Quantum Waves Equations
A Friendly Review of Absorbing Boundary Conditions and Perfectly Matched Layers for Classical and Relativistic Quantum Waves Equations X. Antoine a, E. Lorin b,c and Q. Tang a a Institut Elie Cartan de
More informationNumerical Solution of Time-Dependent Nonlinear Schrödinger Equations Using Domain Truncation Techniques Coupled With Relaxation Scheme
Numerical Solution of Time-Dependent Nonlinear Schrödinger Equations Using Domain Truncation Techniques Coupled With Relaxation Scheme Xavier Antoine, Christophe Besse, Pauline Klein To cite this version:
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J NUMER ANAL Vol, No 3, pp 93 95 c 00 Society for Industrial and Applied Mathematics A RELAXATION SCHEME FOR THE NONLINEAR SCHRÖDINGER EQUATION CHRISTOPHE BESSE Abstract In this paper, we present
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 informationNumerical solution of the Schrödinger equation on unbounded domains
Numerical solution of the Schrödinger equation on unbounded domains Maike Schulte Institute for Numerical and Applied Mathematics Westfälische Wilhelms-Universität Münster Cetraro, September 6 OUTLINE
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 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 informationTowards accurate artificial boundary conditions for nonlinear PDEs through examples
Towards accurate artificial boundary conditions for nonlinear PDEs through examples Xavier Antoine, Christophe Besse, Jérémie Szeftel Abstract. The aim of this paper is to give a comprehensive review of
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 informationDiscrete Transparent Boundary Conditions for the Schrödinger Equation. Matthias Ehrhardt, Anton Arnold
Discrete Transparent Boundary Conditions for the Schrödinger Equation Matthias Ehrhardt, Anton Arnold Fakultät für Mathematik und Informatik, Universität des Saarlandes, Postfach 15 11 5, D-6641 Saarbrücken,
More informationNumerical methods for a fractional diffusion/anti-diffusion equation
Numerical methods for a fractional diffusion/anti-diffusion equation Afaf Bouharguane Institut de Mathématiques de Bordeaux (IMB), Université Bordeaux 1, France Berlin, November 2012 Afaf Bouharguane Numerical
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationA quantum heat equation 5th Spring School on Evolution Equations, TU Berlin
A quantum heat equation 5th Spring School on Evolution Equations, TU Berlin Mario Bukal A. Jüngel and D. Matthes ACROSS - Centre for Advanced Cooperative Systems Faculty of Electrical Engineering and Computing
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationBoundary problems for fractional Laplacians
Boundary problems for fractional Laplacians Gerd Grubb Copenhagen University Spectral Theory Workshop University of Kent April 14 17, 2014 Introduction The fractional Laplacian ( ) a, 0 < a < 1, has attracted
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 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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
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 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 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 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 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 informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
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 informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More informationComputation of Nonlinear Schrödinger Equation on an Open Waveguide Terminated by a PML
Copyright 20 Tech Science Press CMES, vol.7, no.4, pp.347-362, 20 Computation of Nonlinear Schrödinger Equation on an Open Waveguide Terminated by a PML Jianxin Zhu, Zheqi Shen Abstract: It is known that
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 informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationFractional order operators on bounded domains
on bounded domains Gerd Grubb Copenhagen University Geometry seminar, Stanford University September 23, 2015 1. Fractional-order pseudodifferential operators A prominent example of a fractional-order pseudodifferential
More informationBose-Einstein condensates in optical lattices: mathematical analysis and analytical approximate formulas
0.5 setgray0 0.5 setgray1 Bose-Einstein condensates in optical lattices: mathematical analysis and analytical approximate formulas IV EBED João Pessoa - 2011 Rolci Cipolatti Instituto de Matemática - UFRJ
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 informationLow dimensional interacting bosons
Low dimensional interacting bosons Serena Cenatiempo PhD defence Supervisors: E. Marinari and A. Giuliani Physics Department Sapienza, Università di Roma February 7, 2013 It was only in 1995 Anderson et
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
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 informationChristophe Besse 1 and Brigitte Bidégaray 1
Mathematical Modelling and Numerical Analysis ESAIM: MAN Modélisation Mathématique et Analyse Numérique Vol. 35, N o 1, 1, pp. 35 55 NUMERICAL STUDY OF SELF-FOCUSING SOLUTIONS TO THE SCHRÖDINGER-DEBYE
More informationSuperconductivity in domains with corners
Superconductivity in domains with corners Virginie BONNAILLIE-NOËL IRMAR, Université Rennes 1 and ENS Cachan Bretagne Each China Normal University, Shanghai 2007.5.15 Outline 1. Introduction 2. Linear
More informationPhase-field systems with nonlinear coupling and dynamic boundary conditions
1 / 46 Phase-field systems with nonlinear coupling and dynamic boundary conditions Cecilia Cavaterra Dipartimento di Matematica F. Enriques Università degli Studi di Milano cecilia.cavaterra@unimi.it VIII
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationIntegral Representation Formula, Boundary Integral Operators and Calderón projection
Integral Representation Formula, Boundary Integral Operators and Calderón projection Seminar BEM on Wave Scattering Franziska Weber ETH Zürich October 22, 2010 Outline Integral Representation Formula Newton
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More information8.324 Relativistic Quantum Field Theory II MIT OpenCourseWare Lecture Notes Hong Liu, Fall 2010 Lecture 3
Lecture 3 8.324 Relativistic Quantum Field Theory II Fall 200 8.324 Relativistic Quantum Field Theory II MIT OpenCourseWare Lecture Notes Hong Liu, Fall 200 Lecture 3 We begin with some comments concerning
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationWeek 4 Lectures, Math 6451, Tanveer
1 Diffusion in n ecall that for scalar x, Week 4 Lectures, Math 6451, Tanveer S(x,t) = 1 exp [ x2 4πκt is a special solution to 1-D heat equation with properties S(x,t)dx = 1 for t >, and yet lim t +S(x,t)
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 informationThe Chern-Simons-Schrödinger equation
The Chern-Simons-Schrödinger equation Low regularity local wellposedness Baoping Liu, Paul Smith, Daniel Tataru University of California, Berkeley July 16, 2012 Paul Smith (UC Berkeley) Chern-Simons-Schrödinger
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 informationSingularity Formation in Nonlinear Schrödinger Equations with Fourth-Order Dispersion
Singularity Formation in Nonlinear Schrödinger Equations with Fourth-Order Dispersion Boaz Ilan, University of Colorado at Boulder Gadi Fibich (Tel Aviv) George Papanicolaou (Stanford) Steve Schochet (Tel
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 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 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 informationSymmetry Preserving Numerical Methods via Moving Frames
Symmetry Preserving Numerical Methods via Moving Frames Peter J. Olver University of Minnesota http://www.math.umn.edu/ olver = Pilwon Kim, Martin Welk Cambridge, June, 2007 Symmetry Preserving Numerical
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 informationOrdinary differential equations. Phys 420/580 Lecture 8
Ordinary differential equations Phys 420/580 Lecture 8 Most physical laws are expressed as differential equations These come in three flavours: initial-value problems boundary-value problems eigenvalue
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 information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationA Review of Transparent and Artificial Boundary Conditions Techniques for Linear and Nonlinear Schrödinger Equations
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 4, No. 4, pp. 79-796 Commun. Comput. Phys. October 8 REVIEW ARTICLE A Review of Transparent and Artificial Boundary Conditions Techniques for Linear and Nonlinear
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
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 informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationNTNU Trondheim, Institutt for fysikk
NTNU Trondheim, Institutt for fysikk Examination for FY3464 Quantum Field Theory I Contact: Michael Kachelrieß, tel. 998971 Allowed tools: mathematical tables 1. Spin zero. Consider a real, scalar field
More informationBoundary problems for fractional Laplacians and other mu-transmission operators
Boundary roblems for fractional Lalacians and other mu-transmission oerators Gerd Grubb Coenhagen University Geometry and Analysis Seminar June 20, 2014 Introduction Consider P a equal to ( ) a or to A
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 informationExistence and stability of solitary-wave solutions to nonlocal equations
Existence and stability of solitary-wave solutions to nonlocal equations Mathias Nikolai Arnesen Norwegian University of Science and Technology September 22nd, Trondheim The equations u t + f (u) x (Lu)
More informationLecture II Search Method for Lax Pairs of Nonlinear Partial Differential Equations
Lecture II Search Method for Lax Pairs of Nonlinear Partial Differential Equations Usama Al Khawaja, Department of Physics, UAE University, 24 Jan. 2012 First International Winter School on Quantum Gases
More informationNonlinear wave-wave interactions involving gravitational waves
Nonlinear wave-wave interactions involving gravitational waves ANDREAS KÄLLBERG Department of Physics, Umeå University, Umeå, Sweden Thessaloniki, 30/8-5/9 2004 p. 1/38 Outline Orthonormal frames. Thessaloniki,
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 informationDerivation of the Nonlinear Schrödinger Equation. from First Principles. Theodore Bodurov
Annales de la Fondation Louis de Broglie, Volume 30, no 3-4, 2005 343 Derivation of the Nonlinear Schrödinger Equation from First Principles Theodore Bodurov Eugene, Oregon, USA, email: bodt@efn.org ABSTRACT.
More informationNonlinear Gamow Vectors in nonlocal optical propagation
VANGUARD Grant 664782 Nonlinear Gamow Vectors in nonlocal optical propagation M.C. Braidotti 1,2, S. Gentilini 1,2, G. Marcucci 3, E. Del Re 1,3 and C. Conti 1,3 1 Institute for Complex Systems (ISC-CNR),
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 informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
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 informationOn semilinear elliptic equations with measure data
On semilinear elliptic equations with measure data Andrzej Rozkosz (joint work with T. Klimsiak) Nicolaus Copernicus University (Toruń, Poland) Controlled Deterministic and Stochastic Systems Iasi, July
More informationThermodynamic limit for a system of interacting fermions in a random medium. Pieces one-dimensional model
Thermodynamic limit for a system of interacting fermions in a random medium. Pieces one-dimensional model Nikolaj Veniaminov (in collaboration with Frédéric Klopp) CEREMADE, University of Paris IX Dauphine
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
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 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 informationON SOME RESULTS FOR QUANTUM HYDRODYNAMICAL MODELS
ON SOME RESULTS FOR QUANTUM HYDRODYNAMICAL MODELS PAOLO ANTONELLI, LARS ERIC HIENTZSCH, PIERANGELO MARCATI, AND HAO ZHENG Abstract. In this paper we review some recent results on the existence of finite
More informationHARMONIC ANALYSIS. Date:
HARMONIC ANALYSIS Contents. Introduction 2. Hardy-Littlewood maximal function 3. Approximation by convolution 4. Muckenhaupt weights 4.. Calderón-Zygmund decomposition 5. Fourier transform 6. BMO (bounded
More informationPDEs, Homework #3 Solutions. 1. Use Hölder s inequality to show that the solution of the heat equation
PDEs, Homework #3 Solutions. Use Hölder s inequality to show that the solution of the heat equation u t = ku xx, u(x, = φ(x (HE goes to zero as t, if φ is continuous and bounded with φ L p for some p.
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 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 informationIntroduction to Pseudodifferential Operators
Introduction to Pseudodifferential Operators Mengxuan Yang Directed by Prof. Dean Baskin May, 206. Introductions and Motivations Classical Mechanics & Quantum Mechanics In classical mechanics, the status
More informationExponentially Accurate Semiclassical Tunneling Wave Functions in One Dimension
Exponentially Accurate Semiclassical Tunneling Wave Functions in One Dimension Vasile Gradinaru Seminar for Applied Mathematics ETH Zürich CH 8092 Zürich, Switzerland, George A. Hagedorn Department of
More informationA posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem
A posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem Larisa Beilina Michael V. Klibanov December 18, 29 Abstract
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 informationOn some weighted fractional porous media equations
On some weighted fractional porous media equations Gabriele Grillo Politecnico di Milano September 16 th, 2015 Anacapri Joint works with M. Muratori and F. Punzo Gabriele Grillo Weighted Fractional PME
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More information