STRANG SPLITTING METHODS FOR A QUASILINEAR SCHRÖDINGER EQUATION - CONVERGENCE, INSTABILITY AND DYNAMICS 1. INTRODUCTION.
|
|
- Judith Dixon
- 5 years ago
- Views:
Transcription
1 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 Schrödinger equations. We establish the convergence of the scheme for solutions with small initial data. We analyze the linear instability of the numerical scheme, which explains the numerical blow-up of large data solutions and connects to analytical breakdown of regularity of solutions to quasilinear Schrödinger equations. Numerical tests are performed for a modified version of the superfluid thin film equation. 1. INTRODUCTION Consider a general quasi-linear Schrödinger equation (1) i u t = u + u f ( u ) + ug ( u ) g ( u ), for f, g : R R. Such equations can be written as i u t + a j k (u) j k u = F (u, u), () u(0, x) = u 0 (x) u : R R d C m with small initial data in a space with relatively low Sobolev regularity but with some extra decay assumptions. Here a : C m (C m ) d R d d, are smooth functions which we will assume satisfy F : C m (C m ) d C m a(y, z) = I d +O( y + z ), F (y, z) = O( y 3 + z 3 ) near (y, z) = (0,0). Quasilinear equations of this form have arisen in several models. See [40 for a thorough list, but we mention here works related to the superfluid thin-film equation [31, modeling ultrashort pulse lasers [14, 15, and time dependent density functional theory [11. The model we will consider here numerically equates to setting g (s) = f (s) = s, and hence (3) i u t = u + u u + u ( u ). This is a pseudo-attractive version of the superfluid thin-film equation, which is given by (4) i u t = u + u u u ( u ) and can be seen as a leading order contribution to the ultra-short pulse laser models from [14, 15. Existence of solutions to quasilinear equations have been studied analytically in several case, Date: March 31,
2 JIANFENG LU AND J.L. MARZUOLA see [14, 15, 7, 8, 9, 35, 36, 40 and many others. The reason we choose to study (3) is that while similar to (4) in that it is guaranteed to have small data local well-posedness from [36 and hence can be used to verify our numerical convergence results for general quasilinear models, the dynamics of (3) can lead to a breakdown of regularity due to a non-positive definite conserved energy. The model (4) on the other hand has a positive energy quantity and as a result much more stable dynamics. The nonlinear flow will allow interesting singularities to form in the evolution for large enough initial data. In particular, we observe blow-up at a particular amplitude threshold, but these singularities are representative of a breakdown of regularity in the higher derivatives and hence not the standard self-similar style blow-up from the semilinear Schrödinger equation. Such a threshold was observed as an obstruction to local well-posedness using Nash-Moser type arguments in [3. We show analytically that this mechanism for instability is inherited by the Strang-Splitting scheme through a rigorous convergence result and analysis of a finite frequency approximation. Moreover, we observe numerically that this threshold for ill-posedness arises in several different types of initial configuration and is rather robust. However, we note that this threshold is not the numerically observed sharp threshold for long-time well-posedness, as indeed the dynamics are able to drive nearby solutions to this critical amplitude. These features of (3) will be explored in Section. Let us consider the nonlinear part of the equation (5) i v t = v f ( v ) + v g ( v ) g ( v ). Taking the complex conjugate, we have i sv t = sv f ( v ) + sv g ( v ) g ( v ). We calculate i t v = i sv t v + i v t sv (6) = v f ( v ) + v g ( v ) g ( v ) v f ( v ) v g ( v ) g ( v ) = 0, and hence under the evolution (5) the amplitude is conserved. This will be a key property used to develop the numerical scheme. To construct a stable numerical scheme, we consider a Strang splitting method for the quasilinear Schrödinger equation, which is a composition of the exact flows of the differential equations (7) i t u = u and (8) i t u = u f ( u ) + ug ( u ) g ( u ).
3 STRANG SPLITTING FOR QLS 3 More concretely, we approximate u(t n ) with t n = nτ for a step size τ > 0 by u n via u n+1/ = e i τ u n ; (9) u + n+1/ = u n+1/ exp( iτ ( f ( u n+1/ ) + g ( u n+1/ ) g ( u n+1/ ) )) ; u n+1 = e i τ u + n+1/. We note that the scheme is explicit and symmetric, thanks to the amplitude preserving property (6) of (5). One can use a Fourier pseudo-spectral method for the spatial discretization, and hence the flow exp( i τ ) can be efficiently calculated using fast Fourier transform (FFT), and the flow (8) amounts to changing the phase of the solution on each mesh point. Due to the advantage of being structure-preserving, the Strang splitting [43 and higher order splitting schemes (e.g. [44, 50) have been widely applied to a variety of nonlinear Schrödinger equations, mainly semilinear Schrödinger equations, modeling monochromatic light in nonlinear optics, Bose-Einstein condensates, as well as envelope solutions for surface wave trains in fluids. See for example [, 3, 5, 6, 7, 8, 9, 1, 17, 19,, 3, 37, 38, 39, 4, 49. While we focus on the Strang splitting scheme for quasilinear Schrödinger equations, let us also mention that many other time discretization approaches to solve non-linear evolution equations have been developed, including Crank-Nicholson type schemes (see e.g. [41 and also [1 and [0 for applications in studying numerical blow-ups and nonlinear scattering), Magnus expansion approaches ([34 and also the recent review article [10), exponential time-differencing schemes (see e.g. [13, 4), implicit-explicit methods (see e.g. [4), the comparison study in [45, and many others. The convergence of splitting schemes for semilinear Schrödinger equation was analyzed in [16, 18, 1, 33, 4, 46. In the present work, we extend the previous works to the quasi-linear Schrödinger equation. We will prove the convergence of the time-splitting method to the original evolution for the superfluid thin-film equation. The analysis follows the ideas in the seminal contribution by Lubich in [33, where the main tools are the calculus of Lie derivatives. We will emphasize on the regularity of the time flow, for which the behavior of the quasilinear Schrödinger equation is different from the semilinear ones. This is a Lie theoretic approach to the continuous time approximation and Sobolev-based well-posedness results of the second author with J. Metcalfe and D. Tataru in order to model small initial data solutions of finite time intervals [35, 36. The scheme is symplectic and is stable within a range of parameters, motivated by the analysis in [36, where the analysis is done purely in Sobolev spaces H s for s sufficiently large. In addition the Strang splitting method converges in the order τ for time-step τ. Moreover, we are able to extend a linear instability observed in a quasilinear Schrödinger equation in [3 to the numerical scheme used to approximate it, justifying the accuracy of a numerically observed blow-up. We study the dynamics of this blow-up solution using both Gaussian and plane-wave configurations of initial data to observe that the threshold for instability is not the sharp global well-posedness threshold for the equation and can indeed be reached through frequency dynamics on lower-amplitude solutions. This is not the standard blow-up
4 4 JIANFENG LU AND J.L. MARZUOLA through re-scaling of a nonlinear state, but is a frequency instability of sorts that causes high frequencies to grow exponentially in a method akin to a backwards heat map. The result is laid out as follows. We begin with a numerical study of the equation with of the modified superfluid film equation in 1d (3) using the Strang splitting scheme (9). To understand the convergence of the scheme, we prove the convergence of the numerical scheme for small data in Section 3, by using Lie theoretic results and necessary multilinear estimates. To understand the blow-up behavior observed for large enough data, we analyze the stability and instability of the scheme in Section 4. In Section 5, we discuss the regularity of the time flow of the quasilinear Schrödinger equation and the time-splitting scheme. In [3, an L threshold for local wellposedness was observed through use of Frechet derivatives in a Nash-Moser scheme. We will show this similarly arises in analysis about exact plane-wave solutions on the torus using analysis similar to that of [49. In order to establish the differentiability of the numerical solution with respect to time to sufficiently high accuracy, we rely on bounds in a much stronger topology in space.. NUMERICAL RESULTS To test the Strang splitting scheme (9) for quasilinear Schrödinger equations and numerically study the regularity breakdown, we consider the modified superfluid thin film equation (see (3) and also [31, 40) given by (10) i u t + u xx = u u + ( u ) xx u on the domain ( π, π with periodic boundary condition. We have done similar computations for the ultra-short pulse laser equation as described in [14, 15, but no further interesting features of the numerical analysis arose, so we do not present them here for clarity of exposition..1. Symmetric Gaussian initial condition. Let us consider initial conditions given by (11) u(0, x) = ae x /(σ ), on the domain ( π,π with periodic boundary condition. Here σ is the width and a is the amplitude of the Gaussian profile. We calculate the solution up to time T = π/4 with parameters a = 1/5 and σ = 1/5 for the initial condition. A Fourier pseudo-spectral method with N = 56 spatial grid points is used. To verify the second order accuracy of the time splitting scheme, we choose different number of time steps and estimate the error by comparing the numerical solutions to a solution with N t = The results in Table 1 and Figure 1 confirm the second order convergence. Provided the solution remains in H 1, the PDE (10) conserves mass and energy given by (1) M(u) = u dx, E(u) = 1 u x dx + 1 u 4 dx 1 (13) ( u ) x dx. 4 4
5 STRANG SPLITTING FOR QLS 5 N t L -norm H 1 -seminorm e e e e e e e e e e 06 TABLE 1. Numerical error of the time splitting scheme for initial data (11) with a = 1/5 and σ = 1/5 at T = π/ error L error H 1 error N t FIGURE 1. Log-log plot of the numerical error in Table 1 measured in L norm and H 1 seminorm. The numerical scheme conserves the mass conservation law. While there is no energy conservation law [, Table 1, numerically the energy is observed to remain conserved with very high accuracy, which is confirmed in Figure. Due to the nonlinearity of the equation (10), the problem becomes more stiff for initial condition with larger amplitude or derivatives. For the family of Gaussian initial data (11), this means to increase a or reduce σ. We next consider the example with a = 0.65, σ = 1/10 and T = π/4. The problem is considerably more difficult than the previous choice of parameters. We refine the spatial discretization to N = 51 to resolve the oscillatory profile of the solution. The numerical error can be found in Table. We still observe second order accuracy, though in this case, the time step size cannot be too large, otherwise the numerical scheme becomes unstable. We remark that to make the scheme more stable, it is possible to apply Fourier spectrum truncation to eliminate spurious Fourier components of the numerical solution, as introduced in [30. At each time step, we set to zero all Fourier coefficients with amplitude below a certain threshold δ times the maximum amplitude of Fourier coefficients. In practice, for this example, we find the threshold δ = 1e 3 makes the scheme stable with N t = 000 (recall that the solution is not stable for N t = without Fourier truncation). On the other hand, the filtering might introduce inconsistency to the numerical results.
6 6 JIANFENG LU AND J.L. MARZUOLA conserved quantity Mass Energy t FIGURE. Mass and energy conservation for the numerical solution. The initial condition is given by (11) with a = 1/5 and σ = 1.5. The numerical solution is calculated up to time T = 4π and with N t = 3000 time steps. Energy and mass are recorded every 100 time steps. N t L -norm H 1 -seminorm unstable e e e e e e e e e e 04 TABLE. Numerical error of the time splitting scheme for initial data (11) with a = 0.65 and σ = 1/10 at T = π/4. The error is estimated by comparing with the numerical solution with N t = If we further increase the amplitude of the initial condition, the numerical results indicate a blow-up behavior for the PDE. We increase the amplitude to a = 0.65 while keeping σ = 1/10. The numerical solution is calculated up to T = Figure 3 shows max u(, t) as a function of t for different choices of time step sizes. The sudden jump and exponential increase of the magnitude of maximum around t = indicates a numerical blow-up of the solution. Note that the onset point of the behavior does not depend on the choice of time step size, indicating that this is not due to numerical instability of the time integration. Here we have chosen N = 4096 spatial grid points. The blow-up behavior persists for further refinement of the spatial discretization. To investigate more closely the above observed blow-up, we study the solution around x = 0 and the time when the blow-up occurs. Note that since the initial condition is symmetric, the
7 STRANG SPLITTING FOR QLS 7.5 max u t x max u t x 10 3 FIGURE 3. max u(, t) as a function of time for Gaussian initial condition with a = 0.65 and σ = 1/10. Five time step sizes are taken, correspond to N t = 10000,0000,40000,80000, and (blue, green, red, cyan, and purple curves) for total simulation time T = π. The bottom panel zooms in the region near the numerical blow-up. The dashed horizontal line indicates the level /. solution has the symmetry u(x, t) = u( x, t). Hence, we plot the absolute value of the solution in Figure 4. The numerical simulation indicates that the solution develops a focusing peak at x = 0 with amplitude close to /. The numerical results suggest that the solution to the PDE becomes unstable for this family of initial conditions when the amplitude reaches around /. To further confirm this, we compare the results for the initial condition with a = 0.65 and σ = 1/10, the solution stays below the amplitude of / as in Figure 5. Numerically, no blow-up is observed for a = The instability for large time step size is caused by pollution in the Fourier spectrum, but not the intrinsic instability of solutions to the PDE.
8 8 JIANFENG LU AND J.L. MARZUOLA 0.75 Abs u(x, t) x FIGURE 4. Snapshots of the absolute value of the numerical solutions around the numerical blow-up. The solution squeezes as time increases and leads to a blow-up. The reference horizontal line is plotted at the value / max u t x 10 3 FIGURE 5. The maximum magnitude of u as a function in time for initial condition with a = 0.65 and σ = 1/10. Compare with the quite different behavior in Figure 3... Plane wave initial conditions. The only exact solution we are aware of for the superfluid equation (10) is the family of wave trains: (14) u(x, t) = a expi (kx ωt). This is a solution to (10) provided that (15) ω = k + a. Since u(x, t) = a for the solution (14) at any x and t, the splitting error of the Strang splitting scheme vanishes, as the potential commutes with operator.
9 STRANG SPLITTING FOR QLS 9 We study the instability by adjusting the amplitude a of the initial data u(x,0) = a exp(i kx) of the solution (14). Figure 6 shows the simulation results for two solutions with initial conditions given by (14) with a = / 10 8 and a = /+10 8 respectively. Even though the amplitudes of the two solutions only differ by 10 8, the behavior of the numerical solutions are completely different. While the numerical solution for the former is stable and accurate, the local truncation error kicks off instability in the latter case. This indicates again that / is the threshold of instability. FIGURE 6. Numerical solution for (14) with a = / 10 8 and a = /+10 8 respectively. The solution to the PDE is unstable when the amplitude is larger than /. We also study multiple Fourier mode solutions to observe if non-local interactions can vary the blow-up profile. Hence, given a pseudospectral discretization scheme keeping the first N Fourier modes, we take initial data of the form (16) u(x,0) = a k j x, j =1expi
10 10 JIANFENG LU AND J.L. MARZUOLA for 0 k 1 k j N. The blow-up behavior of these solutions become more complicated, in particular, oscillations begin to factor around the blow-up after an initial exponential growth of the maximum amplitude, see Figure 7. However, it seems that generically / is still a threshold for blow-up. max u max u t t 1 max u t max u t x 10 3 FIGURE 7. max u(, t) as a function of time for multiple Fourier Mode initial condition with a = 0.65 and frequencies,8 in the -mode setting (top), and,8,14,0 in the 4-mode setting (bottom). Four time step sizes are taken, correspond to N t = 10000,0000,40000, and (blue, green, red, and cyan curves) for total simulation time T = The right panel zooms in the region near the numerical blow-up. The dashed horizontal line indicates the level /. 3. CONVERGENCE OF TIME-SPLITTING SCHEME In this section, we prove the convergence of the Strang splitting scheme for small data. We suppose that the solution u(t) to the modified superfluid thin film equation (10) in 1d is in H 7 for 0 t T, and set Theorem 1. The numerical solution u n given by the Strang splitting scheme (9) with time step size τ > 0 has a second-order error bound in H 1 (17) u n u(t n ) H 1 C (m 7,T )τ, where t n = nτ T and m 7 = max 0 t T u(t) H 7.
11 STRANG SPLITTING FOR QLS 11 Remark. The small data local existence in H s for the cubic quasilinear nonlinear terms is established in Marzuola-Metcalfe-Tataru [35, 36 (See also the works of Poppenberg [40, Kenig-Ponce- Vega and Kenig-Ponce-Rolvung-Vega [5, 6, 7, 8, 9). In particular, in the case of equation (3) in dimension d, there exists a local in time solution in H σ as long as σ > d+5 provided u 0 is sufficiently small. Therefore, the regularity assumption m 7 < holds for sufficiently small data in H σ for σ sufficiently large (and hence in L ). This is far from sharp however and much work must be done to explore the threshold between well-posedness and blow-up. Proof. We will actually prove the result for arbitrary spatial dimension d. The convergence proof follows a Lie theoretic idea of Lubich [33 for semilinear nonlinear Schrödinger equations. Before we begin let us take (18) m k = max u(t) H 0 t T k, k max( 7, d ɛ ), for and ɛ > 0 such that u the solution to (1) with small initial data can be defined in H k using [35, 36. We can approximate the solution through the continuous time generators of the split step equations: (19) (0) i ψ = ψ i ψ = V [ψψ where (1) V [ψ = ψ + ( ψ ). The generators of the split step method can thus be described as exponential maps of the vector fields given by () (3) ˆT (ψ) = i ψ, ˆV (ψ) = iv [ψψ = i [ ψ + ( ψ )ψ. The key estimate we will require is of the type (4) (uv)w H s C u H s+ 4+d +ɛ v s+ 4+d H +ɛ w H s for s > 0 using the L L L L Hölder s inequality and the Sobolev embedding for L. Note, we could likely sharpen our results with multilinear estimates such as (5) (uv)w H s C u H s+ 6+d 3 v H s+ 6+d 3 w H s+ d 3 for s 0 using L 6 L 6 L 6 L Hölder s inequality and the Sobolev embedding for L 6 as well as moving to L p based spaces and applying the techniques of Strichartz estimates, etc. However, for simplicity of exposition we will simply use the Sobolev embedding estimate (4). Before computing Lie derivatives, we want to understand the stability of the evolution generated by ˆV. To do this, we study (6) i ν = V [ψν, ν(0) = ψ.
12 1 JIANFENG LU AND J.L. MARZUOLA For ψ sufficiently regular, it is possible to show that the evolution varies continuously with the choice of initial data in a weak topology. In particular, we can show that given i ν = V [ψν, ν(0) = ψ, i µ = V [φµ, µ(0) = φ, then by looking at the difference of these two evolutions, expanding V [φµ V [ψν = (V [φ V [ψ)µ V [ψ(ν µ) and applying (4) we have (7) µ(t) ν(t) H s ψ φ H s +C 1 t ψ φ H s+ d+4 +ɛ +C t 0 µ(s) ν(s) H s, for any ɛ > 0 and s > d+4 chosen sufficiently large to control the evolution, where C 1, C both depend upon M s = max 0 t T { φ H s, ψ H s }. As a result, a Gronwall type argument shows (8) µ(t) ν(t) H s e C 0τ ψ φ H s+ d+5 where C 0 depends upon m k. Now, to compare the full evolution to the split-step method, we must compute the Lie commutators between generating vector fields: (9) Hence, we observe [ ˆT, ˆV ψ = ( ψ ψ ( ψ )ψ ) [ ψ( ψψ ψ ψ) [ ( ψ ψ)ψ (ψ ψ)ψ + ( ψ ) ψ. (30) [ ˆT, ˆV (ψ) H 1 C ψ 3 max(5, 4+d H +). In addition, we then can easily compute (31) [ ˆT,[ ˆT, ˆV (ψ) H 1 C ψ 3 max(7, 6+d H +). Setting the vector field Ĥ = ˆT + ˆV, the underlying idea is that the evolution of the full quasilinear Schrödinger equation given by the exact evolution (3) ψ(τ) = exp(τd H )Id(ψ 0 ) when well defined can be compared through a double Duhamel expansion to the split-step generator (33) ψ SS (τ) = exp( 1 τd T )exp(τd V )exp( 1 τd T )Id(ψ 0 ), the error terms of which can be written using the Lie commutators. Indeed, the error estimates come from successive application of the quadrature first order error formula (34) τf ( 1 τ 1 τ) f (s)ds = τ κ 1 (θ)f (θτ)dθ and the second-order error formula (35) τf ( 1 τ 1 τ) f (s)ds = τ 3 κ (θ)f (θτ)dθ
13 STRANG SPLITTING FOR QLS 13 where κ 1 (θ) and κ (θ) are the Peano kernels for the midpoint rule and f (s) = exp((τ s)d T )D V exp(sd T )Id(ψ 0 ) and hence f (s) = e i s [ ˆT, ˆV e i (τ s) ψ 0, f (s) = e i s [ ˆT,[ ˆT, ˆV e i (τ s) ψ 0. Note, the Peano kernels are defined as the integral kernels of the linear transformation L : C k+1 [0,T R such that k f (j ) (0) L(f ) = f T j = 1 j! k! j =0 T 0 κ k (s)f (k+1) (s)d s, hence we observe that using the mid-point rule the f (t/) term vanishes explaining why there is not a quadratic term in (35), though the expressions (34) and (35) can still vary due to the nature of the error term expansions in both cases. Hence, it is essential that for the below we can prove that for our approximation we have f (s) C 3, which very much relates to the analyticity of the linear Schrödinger evolution kernel in the Strang splitting scheme as in particular a generic quasilinear Schrödinger flow cannot be shown to be more than C 0 by the purely dispersive techniques in [35, 36. We will come back to this in more detail in Section 5. Applying (3), (33), (30), (31) and (7) in succession as in [33 gives (36) u n u(t n ) H 1 C (m K0,T )τ for t n = nτ T and K 0 = max(7, d+6 + ɛ). In order to obtain the τ convergence here, it is important to compute the double commutator bound leading to (33) in order to expand out to 3rd order in the Lie derivatives. However, we note the same quadratic convergence would hold in L with only K 0 = max(5, d+5 + ɛ) as then the double Duhamel commutator would not be required and we would be mostly restricted by the well-posedness threshold for (1). Remark. So far we have considered the convergence of the time-splitting flow to the flow of the original PDE. We further discretize the spatial degree of freedom using a Fourier pseudo-spectral method. The convergence of the fully discretized scheme follows if we can show that the fully discretized scheme converges to the time-splitting flow. Though this is beyond the scope of our aim in this work. See for instance [18, 4, 47, 48 for analysis of fully discretized scheme for semilinear Schrödinger equations. Remark. Note, we have chosen a model example that demonstrates the techniques quite nicely in our choice of (19) and (0), but the flow of the metric term and the nonlinear term in () are both well defined and have well defined Lie derivatives provided one can prove the existence of
14 14 JIANFENG LU AND J.L. MARZUOLA solutions. Hence, the splitting process can be done in complete generality using the method (37) (38) i ψ = i g i j [ψ j ψ i ψ = V [ψ, ψψ. Each evolution would require time-dependent kinetic energy terms and multilinear estimates akin to those in [35, 36. However, this gives a nonlinear dependence in the Hamiltonian related to the evolution of (37). The nonlinearity causes troubles in designing structure preserving schemes and also difficulty associated with the regularity of the solution map with respect to time in order to approximate using Lie derivatives. 4. STABILITY AND INSTABILITY OF THE NUMERICAL SCHEME As discussed in Section, for large data, we observe blow-up behavior in the numerical study. In this section, we will investigate the numerical instability of the scheme, which shed some light on the blow-up behavior Linear stability analysis for wave train. For the uniform wave trains solution (14), we study the stability for perturbations around the solution. Consider a perturbed solution of the form (39) u(x, t) = u 0 (x, t)(1 + ε(x, t)), where u 0 is the plane-wave solution ae i (kx ωt) and ε 1. To the leading order, we get (40) iε t + i kε x + ε xx = a (ε +sε) + a (ε +sε) xx. Let us expand ε in Fourier series (with ξ n = n for ε periodic on [0,π) (41) ε(x, t) = ε n (t)exp(iξ n x). n= The equation of ε can then be written as a system of ODEs ( ) ( ) d εn εn (4) = G n, dt ŝε n ŝε n where (43) G n = i ( kξn ξ n a + a ξ n a + a ξ n a a ξ n kξ n + ξ n + a a ξ n The eigenvalues of G n, λ n is given by (44) λ n = i kξ n ± ξ n ξ n a + a ξ n The solution becomes unstable if one of the eigenvalues has a positive real eigenvalue or equivalently (45) a ξ n a ξ n > 0 A sufficient condition for stability is (46) a /. )
15 STRANG SPLITTING FOR QLS 15 Note that the stability threshold / agrees with the numerical observations in Section. 4.. Linear Stability Analysis for the Strang splitting algorithm. To study the stability of general initial data, where explicit solution is not available, we linearize around a solution to (1). Locally, the linear instability is essentially equivalent to the plane wave instability observed in Section., and analyzed in the previous subsection. The key in the observation is that the instability occurs for all k. This has been done in [3, Equation (8), where one observes that perturbation around a solution u = w + z leads to (47) z t = a [M w z +G w z + H w z + f (t), z(0) = g, where for w = w 1 + i w we observe M w = [ w1 w w 1 1 w 1 w 1 w, which has determinant 1 w. The matrix functions G w and H w come from the linearization and will be expressed in full in (48) below. Using this linearization and a Frechet based iteration argument in the space H, the authors then show local well-posedness for small data solutions to equation of the form (1). To understand the instability of the numerical scheme, we linearize the discretized Strang- Splitting algorithm and show that the linear instability threshold in the continuous problem exists in the discretized version as well. Letting u = w + z for some solution w of (4) and generally write h(x, t) = h 1 (x, t) + i h (x, t) for any complex function h, we have that the linearized continuous PDE (47) takes the form (see also [3) [ [ z1 w1 w (w = 1) [ [ z1 w z (1 w t 1 ) w + + w 1 w w1 + w (48) 1w z [3w1 + w w 1w [ [ 4w w 1 4w w z1 + [4w 1 w 1 [4w 1 w z [ w w 1 j =1 + (w j w j )) + w w [ j =1 (w j w j ) + w 1 w 1 w 1 w [ z1 z [ z1 z. We then have [ w1 w w M w = 1 1 w1 w 1 w [ 4w w1 T 4w w T, G w = 4w 1 w1 T 4w 1 w T, H w = [ w + w 1 w + w w 1 H 1 H 1 [w 1 w + w 1 w where H 1 = w 1 + w + w 1 w 1 + w 1 w 1 + 4w w + w w, H 1 = 3w 1 + w + 4w 1 w 1 + w 1 w 1 + w w + w w + 4w 1 w 1.
16 16 JIANFENG LU AND J.L. MARZUOLA form We wish to compare the linearization of the full PDE to the discretized linearization of the (49) z n+1/ = e i τ z n ; z + n+1/ = exp( iτ ( w n+1/ + ( w n+1/ ) ))( ) I d iτ H w z n+1/ n+1/ ; z n+1 = e i τ z + n+1/, where taking w n+1/ = w 1 + i w, we have [ w1 H w = w w [ w n+1/ w1 w + + w 1 w w1 + w 1w [3w1 + w w 1w [ 4w w 1 4w n w + [4w 1 w 1 [4w 1 w [ w w 1 j =1 + (w j w j ) + w w [ j =1 (w j w j ) + w 1 w 1 [w 1 w. To address the linear stability of the Strang-Splitting scheme, consider a linearly unstable mode corresponding to (48) such that z is regular. The linearized splitting scheme to the leading order works as (50) ( ( )) z n+1/ = I d + 1 τ z n ; τ 0 ( z + w1 n+1/ (I = w n w d + τ )) n w1 z n+1/ n w 1 w ; n ( ( )) z n+1 = I d + 1 τ z + n+1/ τ 0, and hence ( w1 w (w (51) z n+1 = z n + τ 1) ) (1 w1 ) w z n + O (τ ). 1w As a result, if w 1, w were constant (as we can assume only locally with any accuracy) and w 1 > 0, we observe that each Fourier mode z n+1,k can be approximated by the linearized dynamical system (5) z n+1,k = z n,k + τ ( w1 w k (w ) 1)k (w 1 1)k w 1 w k which has eigenvalues 1 ± τk w 1 and hence would clearly grow exponentially for w > /. For sufficiently large k, we observe that all the Fourier modes of w1, w are small perturbations and hence exponential growth will occur just as in the backwards heat flow generated from the continuous approximation. Indeed, for τ sufficiently small, we have that the linearized Strang-Splitting flow well approximates the unstable backwards heat flow and hence displays linear instability. However, of course the nonlinear effects are ignored in this computation. z n,k
17 STRANG SPLITTING FOR QLS 17 To make this more precise, we observe from [36 for the full PDE model (3) that we can construct initial data for (4) having has L norm larger than / but sufficiently localized in frequency such that the solution exists locally in time. A simple example of a solution with a global existence time and initial L norm larger than / is an exact plane wave solution for the periodic problem. However, perturbations of such exact solutions are still linearly unstable as calculated in Section 4.1. The backward heat equation that represents the linearization of the continuous model exists locally in time when considering frequency localized data, however, this time scale will depend on the frequency cut-off and potentially be quite short given the nonlinear interactions. Using Theorem 1, we observe that in the semi-discretized equation choosing τ 1 sufficiently small compared to the scale of local existence for (4), the numerical solution computed using the Strang-Splitting scheme is O (τ ). Of course, from (51), on a single time-step, the Strang-Splitting solution has a polynomial instability. However, the linearized equations in (48) give exponential growth dynamics for the full solution on the time scale of local existence. Since on the scale of existence, the numerical solution remains O (τ ), the linear instability is inherited by the numerical solution over repeated iterations of the time step. 5. REGULARITY OF THE TIME EVOLUTION In the analysis of the convergence of the Strang splitting scheme, we have used the analyticity of the linear Schrödinger evolution. This is in general however not true for the quasilinear Schrödinger evolution. The Strang splitting scheme actually regularizes the time flow of the original PDE. In this section, we give some further discussion for the regularity of the time evolution. Let us show that the solution map of the quasilinear Schrödinger evolution is continuous in time. The continuity partially hinges upon the proof of uniqueness for the evolution of (1). In particular, take two solutions to (1), say u 1 and u. Setting v = u 1 u and linearizing (1), we have an equation of the form { i vt + g j k (u) j k v +V v +W v = 0, v(0, x) = u 1 (0) u (0) with V = V (u 1, u 1,u, u ), W = h(u 1, u 1,u, u ) + g (u 1,u ) u 1 for functions V, h and g related to Taylor expanding the metric and the nonlinearity. Then, to solve this linear equation, we use [36, Proposition 5.1 (see also [35, Proposition 5.) to show that the weak Lipschitz bound (53) v L H σ v(0) H σ holds for any 0 σ s 1 via energy estimates on the linearized equation, where we recall that the initial condition lies in H s for s > d+5. In the well-posedness result for the linearized version of (1), see Proposition 5.1 of [36 for instance, we have at most continuity of the solution map with respect to time in the H s norm however. The key ideas to the proof follow from the theory of frequency envelopes as discussed in both [35, 36, Sections and 5, where it is proven that the size of a dyadic frequency component of
18 18 JIANFENG LU AND J.L. MARZUOLA the solution to (1) in a natural energy space is bounded by a uniform constant times the corresponding dyadic frequency component of the initial data in H s. To be more precise, we shall use a Littlewood-Paley decomposition of the spatial frequencies, S i (D) = 1, i=0 where S i localizes to frequency ξ [ i 1, i+1 for i > 0 and to frequencies ξ for i = 0. By a frequency envelope, we recall from [35, Section.4 we mean that given a translation invariant space U such that u U S k u U, a frequency envelope for u in U is a positive sequence a j so that k=0 (54) S j u U a j u U, a j 1. We say that a frequency envelope is admissible if a 0 1 and it is slowly varying, a j δ j k a k, j,k 0, 0 < δ 1. An admissible frequency envelope always exists, say by (55) a j = δj + u 1 U max δ j k S k u U. k Abusing notation and avoiding for simplicity the atomic space formulations in [35, 36, we rely upon a uniform bound over the evolution such that effectively (56) u L H s u 0 H s. We note that the L H s norm appears the estimate here is due to the cubic interactions in the nonlinearity and the compactness of our domain, otherwise one must enforce further summability as in [35. A key estimate is the following proposition. Proposition (Proposition 5.3, [35; Proposition 5.4, [36). Let u be a small data solution to (1), which satisfies (56). Let {a j } be an admissible frequency envelope for the initial data u 0 in H s. Then {a j } is also a frequency envelope for u in L H s. Once we have Proposition, the continuity of the solution map can be established as is Section 5.7 of [35. Namely, we consider a sequence of initial data {u n 0 } u 0 in H s. Frequency envelope bounds can then be chosen such that there exists a uniform N ɛ for which a (n) N ɛ ɛ for all n, which gives a uniform upper bound by Proposition on the high frequencies of each corresponding solution u (n) to (1) with initial data u (n) 0 in the L H s norm. Separating into low and high frequencies, using the smallness of the high frequencies and the uniform convergence in weaker Sobolev norms provided by (53), the result follows. Let us emphasize however that we generally gain no more than continuity of the solution map from such arguments. Hence, in order to accurately compare the flow of the full solution map
19 STRANG SPLITTING FOR QLS 19 defined by (1) and that of the Strang splitting method, we rely on differentiating the equation and the balancing of spatial and time regularity, as in Section CONCLUSION In this work, we have studied the Strang splitting scheme for quasilinear nonlinear Schrödinger equations. The splitting scheme is proved to have second order convergence for small data. We further investigate the regularity of the time flow and the instability of the numerical scheme which leads to blow-ups observed numerically. Our work is motivated by numerical approaches towards time-dependent density functional theory computations as discussed in the documentation of the software package Octopus 1 and also [11 and references therein. The mathematical analysis and numerical schemes for timedependent density functional theory will be an exciting direction to explore in the future. Acknowledgements. The first author was supported in part by Alfred P. Sloan Foundation and by National Science Foundation under grant number DMS He would like to thank Weizhu Bao for helpful discussion. The second author was supported as a guest lecturer at Karlsruhe Institute of Technology in Summer 013 as well as by NSF grant DMS and wishes to especially thank his collaborators Jason Metcalfe and Daniel Tataru for introducing him to quasilinear Schrödinger theory. We also wish the thank the anonymous reviewers for many helpful comments and generally improving the exposition in the manuscript. REFERENCES [1 G.D. Akrivis, V. A. Dougalis, O. A. Karakashian, and W. R. McKinney: Numerical approximation of blow-up of radially symmetric solutions of the nonlinear Schrödinger equation, SIAM J. Sci. Comput. 5, No. 1 (003), [ X. Antoine, W. Bao and C. Besse: Computational methods for the dynamics of the nonlinear Schrödinger/Gross- Pitaevskii equations (013). [3 H. Appel and E.K.U. Gross: Static and time-dependent many-body effects via density-functional theory, Quantum simulations of complex many-body systems: From theory to algorithms, NIC Series 10, John von Neumann Institute for Computing Press (00), [4 U.M. Ascher, S.J. Ruuth, and B.T.R. Wetton: Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations, SIAM. J. Numer. Anal. 3 (1995), [5 W. Bao and Y. Cai: Mathematical theory and numerical methods for Bose-Einstein condensates, Kin. Rel. Mod. 6, No. 1 (013), [6 W. Bao, D. Jaksch, and P.A. Markowich: Numerical solution of the Gross-Pitaevskii equation for Bose-Einstein condensation, J. Comput. Phys. 187 (003), [7 W. Bao, S. Jin, and P.A. Markowich: Numerical study of time-splitting spectral discretizations of nonlinear Schrödinger equations in the semiclassical regimes, SIAM J. Sci. Comput. 5 (003), [8 W. Bao, N. Mauser, and H.P. Stimming: Effective one particle quantum dynamics of electrons: A numerical study of the Schrödinger-Poisson-Xα model, Comm. Math. Sci. 1 (003), [9 W. Bao and J. Shen: A fourth-order time-splitting Laguerre-Hermite pseudospectral method for Bose-Einstein condensates, SIAM J. Sci. Comput. 6 (005), [10 S. Blanes, F. Casas, J.A. Oteo and J. Ros: The Magnus expansion and some of its applications, Phys. Rep. 470 (009),
20 0 JIANFENG LU AND J.L. MARZUOLA [11 A. Castro, M.A.L. Marques and A. Rubio: Propagators for the time-dependent Kohn-Sham equations, J. Chem. Phys. 11, No. 8 (004), 345. [1 S.A. Chin: Higher-order splitting algorithms for solving the nonlinear Schrödinger equation and their instabilities, Phys. Rev. E 76 (007) [13 S.M. Cox and P.C. Matthews: Exponential time differencing for stiff systems. J. Comput. Phys. 176 (00), [14 A. De Bouard, N. Hayashi and J.-C. Saut: Global existence of small solutions to a relativistic nonlinear Schrödinger equation. Comm. Math. Phys. 189, (1997), [15 A. De Bouard, N. Hayashi, P.I. Naumkin and J.-C. Saut: Scattering problem and asymptotics for a relativistic nonlinear Schrödinger equation. Nonlinearity 1, (1999), [16 S. Descombes and M. Thalhammer: The Lie-Trotter splitting for nonlinear evolutionary problems with critical parameters: a compact local error representation and application to nonlinear Schrödinger equations in the semiclassical regime, IMA J. Numer. Anal. 33 (013) [17 E. Faou, V. Gradinaru, and Ch. Lubich: Computing semi-classical quantum dynamics with Hagedorn wavepackets, SIAM J. Sci. Comput [18 L. Gauckler: Convergence of a split-step Hermite method for the Gross-Pitaevskii equation, IMA J. Numer. Anal. 31 (011) [19 R.H. Hardin, F.D. Tappert: Applications of the split-step Fourier method to the numerical solution of nonlinear and variable coefficient wave equations. SIAM Rev. 15 (1973), 43. [0 J. Holmer, J. Marzuola and M. Zworski: Soliton splitting by external delta potentials, J. Nonlin. Sci. 17, No. 4 (007), [1 T. Jahnke and Ch. Lubich: Error bounds for exponential operator splittings, BIT 40 (000), [ S. Jin: Schrödinger equation: Computation, Encyclopedia of Applied and Computational Mathematics (ed. by B. Engquist), to appear. [3 S. Jin, P.A. Markowich, and C. Sparber: Mathematical and computational methods for semiclassical Schrödinger equations, Acta Numer. 0 (011), [4 A.-K. Kassam and L.N. Trefethen: Fourth-order time-stepping for stiff PDEs, SIAM J. Sci. Comput. 6 (005), [5 C. E. Kenig, G. Ponce, and L. Vega: Small solutions to nonlinear Schrödinger equations. Ann. Inst. H. Poincaré Anal. Non Linéaire 10 (1993), [6 C. E. Kenig, G. Ponce, and L. Vega: Smoothing effects and local existence theory for the generalized nonlinear Schrödinger equations. Invent. Math. 134 (1998), [7 C. E. Kenig, G. Ponce, and L. Vega: The Cauchy problem for quasi-linear Schrödinger equations. Invent. Math. 158 (004), [8 C. E. Kenig, G. Ponce, C. Rolvung, and L. Vega: The general quasilinear ultrahyperbolic Schrödinger equation. Adv. Math. 196, No. (005), [9 C. E. Kenig, G. Ponce, C. Rolvung, and L. Vega: Variable coefficient Schrödinger flows for ultrahyperbolic operators. Adv. Math. 06, No. (006), [30 R. Krasny: A study of singularity formation in a vortex sheet by the point-vortex approximation. J. Fluid Mech. 167 (1986), [31 S. Kurihara: Large-amplitude quasi-solitons in superfluid films. J. Phys. Soc. Japan 50 (1981), [3 H. Lange, M. Poppenberg and H. Teismann: Nash Moser methods for the solution of quasilinear Schrödinger equations. Comm. PDE, 4, No. 7-8 (007), [33 Ch. Lubich: On splitting methods for Schrödinger-Poisson and cubic nonlinear Schrödinger equations. Math. Comp. 77, No. 64 (008), [34 W. Magnus: On the exponential solution of differential equations for a linear operator. Comm. Pure Appl. Math. 7 (1954), [35 J. Marzuola, J. Metcalfe, and D. Tataru: Quasilinear Schrödinger equations I: Small data and quadratic interactions. Adv. Math., 31, No. (01),
21 STRANG SPLITTING FOR QLS 1 [36 J. Marzuola, J. Metcalfe, and D. Tataru: Quasilinear Schrödinger equations II: Small data and cubic nonlinearities. to appear in Kyoto J. Math. (013), preprint available at arxiv: [37 R.I. MacLachlan and R. Quispel: Splitting methods, Acta Numer., 11 (00), [38 D. Pathria and J. L. Morris: Pseudo-spectral solution of nonlinear Schrödinger equations, J. Comput. Phys., 87 (1990), [39 V.M. Pérez-García and X. Liu: Numerical methods for the simulation of trapped nonlinear Schrödinger systems, Appl. Math. Comp. 144 (003), [40 M. Poppenberg: On the local well-posedness of Quasilinear Schrödinger equations in arbitrary space dimension. J. Diff. Eqns., 17 (001), [41 J.M. Sanz-Serna: Methods for the numerical solution of the nonlinear Schrödinger equation, Math. Comp., 43 (1984), 1-7. [4 J. Shen and Z.-Q. Wang: Error analysis of the Strang time-splitting Laguerre Hermite/Hermite collocation methods for the Gross-Pitaevskii equation, to appear in Found. Comp. Math. (013). [43 G. Strang: On the construction and comparison of difference schemes. SIAM J. Numer. Anal., 5 (1968), [44 M. Suzuki: Improved Trotter-like formula. Phys. Lett. A, 180 (1993), [45 T.R. Taha and M.J. Ablowitz: Analytical and numerical aspects of certain nonlinear evolution equations. II Numerical, nonlinear Schrödinger equation, J. Comput. Phys., 55 (1984), [46 M. Thalhammer: High-order exponential operator splitting methods for time-dependent Schrödinger equations, SIAM J. Numer. Anal., 46 (008), [47 M. Thalhammer: Convergence Analysis of High-Order Time-Splitting Pseudospectral Methods for Nonlinear Schrödinger Equations, SIAM J. Numer. Anal., 50, No. 6 (01), [48 H. Hofstätter, O. Koch and M. Thalhammer: High-order exponential operator splitting methods for time-dependent Schrödinger equations, Numer. Math., 46 (013), [49 J.A.C. Weideman and B.M. Herbst: Split-step methods for the solution of the nonlinear Schrödinger equation, SIAM J. Numer. Anal., 3 (1986), [50 H. Yoshida: Construction of higher order symplectic integrators, Phys. Lett. A, 150 (1990), DEPARTMENT OF MATHEMATICS AND DEPARTMENT OF PHYSICS, DUKE UNIVERSITY, BOX 9030, DURHAM, NC address: jianfeng@math.duke.edu DEPARTMENT OF MATHEMATICS, UNC-CHAPEL HILL, CB#350 PHILLIPS HALL, CHAPEL HILL, NC address: marzuola@math.unc.edu
High-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 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 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 informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationCUTOFF RESOLVENT ESTIMATES AND THE SEMILINEAR SCHRÖDINGER EQUATION
CUTOFF RESOLVENT ESTIMATES AND THE SEMILINEAR SCHRÖDINGER EQUATION HANS CHRISTIANSON Abstract. This paper shows how abstract resolvent estimates imply local smoothing for solutions to the Schrödinger equation.
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 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 informationSCATTERING FOR THE TWO-DIMENSIONAL NLS WITH EXPONENTIAL NONLINEARITY
SCATTERING FOR THE TWO-DIMENSIONAL NLS WITH EXPONENTIAL NONLINEARITY S. IBRAHIM, M. MAJDOUB, N. MASMOUDI, AND K. NAKANISHI Abstract. We investigate existence and asymptotic completeness of the wave operators
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 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 informationTHE RADIUS OF ANALYTICITY FOR SOLUTIONS TO A PROBLEM IN EPITAXIAL GROWTH ON THE TORUS
THE RADIUS OF ANALYTICITY FOR SOLUTIONS TO A PROBLEM IN EPITAXIAL GROWTH ON THE TORUS DAVID M AMBROSE Abstract A certain model for epitaxial film growth has recently attracted attention, with the existence
More informationGlobal well-posedness for semi-linear Wave and Schrödinger equations. Slim Ibrahim
Global well-posedness for semi-linear Wave and Schrödinger equations Slim Ibrahim McMaster University, Hamilton ON University of Calgary, April 27th, 2006 1 1 Introduction Nonlinear Wave equation: ( 2
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 informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
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 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 informationDispersive numerical schemes for Schrödinger equations
Dispersive numerical schemes for Schrödinger equations Enrique Zuazua joint work with L. Ignat zuazua@bcamath.org Basque Center for Applied Mathematics (BCAM), Bilbao, Basque Country, Spain IMA Workshop:
More informationA REMARK ON AN EQUATION OF WAVE MAPS TYPE WITH VARIABLE COEFFICIENTS
A REMARK ON AN EQUATION OF WAVE MAPS TYPE WITH VARIABLE COEFFICIENTS DAN-ANDREI GEBA Abstract. We obtain a sharp local well-posedness result for an equation of wave maps type with variable coefficients.
More informationEnergy transfer model and large periodic boundary value problem for the quintic NLS
Energy transfer model and large periodic boundary value problem for the quintic NS Hideo Takaoka Department of Mathematics, Kobe University 1 ntroduction This note is based on a talk given at the conference
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 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 informationDISPERSIVE EQUATIONS: A SURVEY
DISPERSIVE EQUATIONS: A SURVEY GIGLIOLA STAFFILANI 1. Introduction These notes were written as a guideline for a short talk; hence, the references and the statements of the theorems are often not given
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 informationGlobal well-posedness for KdV in Sobolev spaces of negative index
Electronic Journal of Differential Equations, Vol. (), No. 6, pp. 7. ISSN: 7-669. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Global well-posedness for
More informationComputer simulation of multiscale problems
Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20
More informationA PHYSICAL SPACE PROOF OF THE BILINEAR STRICHARTZ AND LOCAL SMOOTHING ESTIMATES FOR THE SCHRÖDINGER EQUATION
A PHYSICAL SPACE PROOF OF THE BILINEAR STRICHARTZ AND LOCAL SMOOTHING ESTIMATES FOR THE SCHRÖDINGER EQUATION TERENCE TAO Abstract. Let d 1, and let u, v : R R d C be Schwartz space solutions to the Schrödinger
More informationSANGRADO PAGINA 17CMX24CM. PhD Thesis. Splitting methods for autonomous and non-autonomous perturbed equations LOMO A AJUSTAR (AHORA 4CM)
SANGRADO A3 PAGINA 17CMX24CM LOMO A AJUSTAR (AHORA 4CM) We have considered the numerical integration of non-autonomous separable parabolic equations using high order splitting methods with complex coefficients.
More informationNUMERICAL SIMULATIONS OF THE ENERGY-SUPERCRITICAL NONLINEAR SCHRÖDINGER EQUATION
Journal of Hyperbolic Differential Equations Vol. 7, No. 2 (2010) 279 296 c World Scientific Publishing Company DOI: 10.1142/S0219891610002104 NUMERICAL SIMULATIONS OF THE ENERGY-SUPERCRITICAL NONLINEAR
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 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 informationA Numerical Study of the Focusing Davey-Stewartson II Equation
A Numerical Study of the Focusing Davey-Stewartson II Equation Christian Klein Benson Muite Kristelle Roidot University of Michigan muite@umich.edu www.math.lsa.umich.edu/ muite 20 May 2012 Outline The
More informationDispersive Equations and Hyperbolic Orbits
Dispersive Equations and Hyperbolic Orbits H. Christianson Department of Mathematics University of California, Berkeley 4/16/07 The Johns Hopkins University Outline 1 Introduction 3 Applications 2 Main
More informationStrichartz Estimates for the Schrödinger Equation in Exterior Domains
Strichartz Estimates for the Schrödinger Equation in University of New Mexico May 14, 2010 Joint work with: Hart Smith (University of Washington) Christopher Sogge (Johns Hopkins University) The Schrödinger
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 informationExponential integrators and functions of the matrix exponential
Exponential integrators and functions of the matrix exponential Paul Matthews, Stephen Cox, Hala Ashi and Linda Cummings School of Mathematical Sciences, University of Nottingham, UK Introduction to exponential
More informationSemi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations
Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang
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 informationBielefeld Course on Nonlinear Waves - June 29, Department of Mathematics University of North Carolina, Chapel Hill. Solitons on Manifolds
Joint work (on various projects) with Pierre Albin (UIUC), Hans Christianson (UNC), Jason Metcalfe (UNC), Michael Taylor (UNC), Laurent Thomann (Nantes) Department of Mathematics University of North Carolina,
More informationSpectral theory for magnetic Schrödinger operators and applicatio. (after Bauman-Calderer-Liu-Phillips, Pan, Helffer-Pan)
Spectral theory for magnetic Schrödinger operators and applications to liquid crystals (after Bauman-Calderer-Liu-Phillips, Pan, Helffer-Pan) Ryukoku (June 2008) In [P2], based on the de Gennes analogy
More informationSharp blow-up criteria for the Davey-Stewartson system in R 3
Dynamics of PDE, Vol.8, No., 9-60, 011 Sharp blow-up criteria for the Davey-Stewartson system in R Jian Zhang Shihui Zhu Communicated by Y. Charles Li, received October 7, 010. Abstract. In this paper,
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 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 informationBound-state solutions and well-posedness of the dispersion-managed nonlinear Schrödinger and related equations
Bound-state solutions and well-posedness of the dispersion-managed nonlinear Schrödinger and related equations J. Albert and E. Kahlil University of Oklahoma, Langston University 10th IMACS Conference,
More informationCurriculum Vitae for Hart Smith - Updated April, 2016
Curriculum Vitae for Hart Smith - Updated April, 2016 Education Princeton University, 1984 1988. Ph.D. in Mathematics, January, 1989. Thesis supervisor: E. M. Stein. Thesis title: The subelliptic oblique
More informationGevrey regularity in time for generalized KdV type equations
Gevrey regularity in time for generalized KdV type equations Heather Hannah, A. Alexandrou Himonas and Gerson Petronilho Abstract Given s 1 we present initial data that belong to the Gevrey space G s for
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 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 informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
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 informationIntegration of Vlasov-type equations
Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Plasma the fourth state of matter 99% of the visible matter in the universe is made of plasma
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 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 informationStrauss conjecture for nontrapping obstacles
Chengbo Wang Joint work with: Hart Smith, Christopher Sogge Department of Mathematics Johns Hopkins University Baltimore, Maryland 21218 wangcbo@jhu.edu November 3, 2010 1 Problem and Background Problem
More informationBLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED
BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED TAOUFIK HMIDI AND SAHBI KERAANI Abstract. In this note we prove a refined version of compactness lemma adapted to the blowup analysis
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Modified Equation of a Difference Scheme What is a Modified Equation of a Difference
More informationRecent developments on the global behavior to critical nonlinear dispersive equations. Carlos E. Kenig
Recent developments on the global behavior to critical nonlinear dispersive equations Carlos E. Kenig In the last 25 years or so, there has been considerable interest in the study of non-linear partial
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 informationANALYTIC SMOOTHING EFFECT FOR NONLI TitleSCHRÖDINGER EQUATION IN TWO SPACE DIMENSIONS. Citation Osaka Journal of Mathematics.
ANALYTIC SMOOTHING EFFECT FOR NONLI TitleSCHRÖDINGER EQUATION IN TWO SPACE DIMENSIONS Author(s) Hoshino, Gaku; Ozawa, Tohru Citation Osaka Journal of Mathematics. 51(3) Issue 014-07 Date Text Version publisher
More informationOn the stability of filament flows and Schrödinger maps
On the stability of filament flows and Schrödinger maps Robert L. Jerrard 1 Didier Smets 2 1 Department of Mathematics University of Toronto 2 Laboratoire Jacques-Louis Lions Université Pierre et Marie
More informationNONLINEAR PROPAGATION OF WAVE PACKETS. Ritsumeikan University, and 22
NONLINEAR PROPAGATION OF WAVE PACKETS CLOTILDE FERMANIAN KAMMERER Ritsumeikan University, 21-1 - 21 and 22 Our aim in this lecture is to explain the proof of a recent Theorem obtained in collaboration
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 informationON THE REGULARITY OF WEAK SOLUTIONS OF THE 3D NAVIER-STOKES EQUATIONS IN B 1
ON THE REGULARITY OF WEAK SOLUTIONS OF THE 3D NAVIER-STOKES EQUATIONS IN B 1, A. CHESKIDOV AND R. SHVYDKOY ABSTRACT. We show that if a Leray-Hopf solution u to the 3D Navier- Stokes equation belongs to
More informationSome Collision solutions of the rectilinear periodically forced Kepler problem
Advanced Nonlinear Studies 1 (2001), xxx xxx Some Collision solutions of the rectilinear periodically forced Kepler problem Lei Zhao Johann Bernoulli Institute for Mathematics and Computer Science University
More informationSolutions to the Nonlinear Schrödinger Equation in Hyperbolic Space
Solutions to the Nonlinear Schrödinger Equation in Hyperbolic Space SPUR Final Paper, Summer 2014 Peter Kleinhenz Mentor: Chenjie Fan Project suggested by Gigliola Staffilani July 30th, 2014 Abstract In
More informationA Time-Splitting Spectral Method for the Generalized Zakharov System in Multi-dimensions
A Time-Splitting Spectral Method for the Generalized Zakharov System in Multi-dimensions Shi Jin and Chunxiong Zheng Abstract The generalized Zakharov system couples a dispersive field E (scalar or vectorial)
More informationRANDOM PROPERTIES BENOIT PAUSADER
RANDOM PROPERTIES BENOIT PAUSADER. Quasilinear problems In general, one consider the following trichotomy for nonlinear PDEs: A semilinear problem is a problem where the highest-order terms appears linearly
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 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 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 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 informationGrowth of Sobolev norms for the cubic NLS
Growth of Sobolev norms for the cubic NLS Benoit Pausader N. Tzvetkov. Brown U. Spectral theory and mathematical physics, Cergy, June 2016 Introduction We consider the cubic nonlinear Schrödinger equation
More informationNonlinear Schrödinger Equation BAOXIANG WANG. Talk at Tsinghua University 2012,3,16. School of Mathematical Sciences, Peking University.
Talk at Tsinghua University 2012,3,16 Nonlinear Schrödinger Equation BAOXIANG WANG School of Mathematical Sciences, Peking University 1 1 33 1. Schrödinger E. Schrödinger (1887-1961) E. Schrödinger, (1887,
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 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 informationMIT (Spring 2014)
18.311 MIT (Spring 014) Rodolfo R. Rosales May 6, 014. Problem Set # 08. Due: Last day of lectures. IMPORTANT: Turn in the regular and the special problems stapled in two SEPARATE packages. Print your
More informationSMOOTHING ESTIMATES OF THE RADIAL SCHRÖDINGER PROPAGATOR IN DIMENSIONS n 2
Acta Mathematica Scientia 1,3B(6):13 19 http://actams.wipm.ac.cn SMOOTHING ESTIMATES OF THE RADIAL SCHRÖDINGER PROPAGATOR IN DIMENSIONS n Li Dong ( ) Department of Mathematics, University of Iowa, 14 MacLean
More informationIntegro-differential equations: Regularity theory and Pohozaev identities
Integro-differential equations: Regularity theory and Pohozaev identities Xavier Ros Oton Departament Matemàtica Aplicada I, Universitat Politècnica de Catalunya PhD Thesis Advisor: Xavier Cabré Xavier
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 informationDunkl operators and Clifford algebras II
Translation operator for the Clifford Research Group Department of Mathematical Analysis Ghent University Hong Kong, March, 2011 Translation operator for the Hermite polynomials Translation operator for
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationComplex Behavior in Coupled Nonlinear Waveguides. Roy Goodman, New Jersey Institute of Technology
Complex Behavior in Coupled Nonlinear Waveguides Roy Goodman, New Jersey Institute of Technology Nonlinear Schrödinger/Gross-Pitaevskii Equation i t = r + V (r) ± Two contexts for today: Propagation of
More informationUNIQUENESS OF POSITIVE SOLUTION TO SOME COUPLED COOPERATIVE VARIATIONAL ELLIPTIC SYSTEMS
TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9947(XX)0000-0 UNIQUENESS OF POSITIVE SOLUTION TO SOME COUPLED COOPERATIVE VARIATIONAL ELLIPTIC SYSTEMS YULIAN
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 informationOn the local existence for an active scalar equation in critical regularity setting
On the local existence for an active scalar equation in critical regularity setting Walter Rusin Department of Mathematics, Oklahoma State University, Stillwater, OK 7478 Fei Wang Department of Mathematics,
More informationGlobal unbounded solutions of the Fujita equation in the intermediate range
Global unbounded solutions of the Fujita equation in the intermediate range Peter Poláčik School of Mathematics, University of Minnesota, Minneapolis, MN 55455, USA Eiji Yanagida Department of Mathematics,
More informationWell-posedness for the Fourth-order Schrödinger Equations with Quadratic Nonlinearity
Well-posedness for the Fourth-order Schrödinger Equations with Quadratic Nonlinearity Jiqiang Zheng The Graduate School of China Academy of Engineering Physics P. O. Box 20, Beijing, China, 00088 (zhengjiqiang@gmail.com)
More informationUNIFORM DECAY OF SOLUTIONS FOR COUPLED VISCOELASTIC WAVE EQUATIONS
Electronic Journal of Differential Equations, Vol. 16 16, No. 7, pp. 1 11. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu UNIFORM DECAY OF SOLUTIONS
More informationAdaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear Degenerate Reaction-Diffusion Problems
EPJ Web of Conferences 108, 0019 (016) DOI: 10.1051/ epjconf/ 0161080019 C Owned by the authors, published by EDP Sciences, 016 Adaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear
More informationarxiv: v1 [physics.flu-dyn] 14 Jun 2014
Observation of the Inverse Energy Cascade in the modified Korteweg de Vries Equation D. Dutykh and E. Tobisch LAMA, UMR 5127 CNRS, Université de Savoie, Campus Scientifique, 73376 Le Bourget-du-Lac Cedex,
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
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 informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationOn the Linearization of Second-Order Dif ferential and Dif ference Equations
Symmetry, Integrability and Geometry: Methods and Applications Vol. (006), Paper 065, 15 pages On the Linearization of Second-Order Dif ferential and Dif ference Equations Vladimir DORODNITSYN Keldysh
More informationHeat equations with singular potentials: Hardy & Carleman inequalities, well-posedness & control
Outline Heat equations with singular potentials: Hardy & Carleman inequalities, well-posedness & control IMDEA-Matemáticas & Universidad Autónoma de Madrid Spain enrique.zuazua@uam.es Analysis and control
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 information2 Formal derivation of the Shockley-Read-Hall model
We consider a semiconductor crystal represented by the bounded domain R 3 (all our results are easily extended to the one and two- dimensional situations) with a constant (in space) number density of traps
More informationA Composite Runge Kutta Method for the Spectral Solution of Semilinear PDEs
Journal of Computational Physics 8, 57 67 (00) doi:0.006/jcph.00.77 A Composite Runge Kutta Method for the Spectral Solution of Semilinear PDEs Tobin A. Driscoll Department of Mathematical Sciences, University
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 information