Splittings, commutators and the linear Schrödinger equation
|
|
- Sophia Park
- 5 years ago
- Views:
Transcription
1 UNIVERSITY OF CAMBRIDGE CENTRE FOR MATHEMATICAL SCIENCES DEPARTMENT OF APPLIED MATHEMATICS & THEORETICAL PHYSICS Splittings, commutators and the linear Schrödinger equation Philipp Bader, Arieh Iserles, Karolina Kropielnicka & Pranav Singh Minneapolis May 202
2 THE PROBLEM We consider the linear Schrödinger equation i h u t = h2 2 u V (x)u, t 0, x, 2m x2 where 0 < h, and the interaction potential V is a smooth periodic function, given with smooth periodic initial and boundary conditions. We can contemplate this equation on a multivariate torus, i h u t = h2 u V (x)u. 2m However, our purpose here is to explore ideas, rather than presenting them in the greatest possible generality. It is convenient to set ω = h, whereby the equation becomes u t = u icω 2 + iωv (x)u, c = x2 2m. 2
3 EXPONENTIALS AND SPLITTINGS... The conventional approach: replace 2 u x2 by a linear combination of function values (or Fourier modes). This results in the linear system whose exact solution is y = (ω K + ωd)y, t 0, y(0) = y 0, y(t) = e t(ω K+ωD) y 0. It is usual to split the exponential, e.g. with the Strang splitting e 2 tω K e tωd e 2 tω K = e t(ω K+ωD) + O ( t 3) or higher-order splittings. Typically, such high-order splittings are obtained either via the Yošida device or similar techniques, and result in palindromic expressions of the form e α tω K e β tωd e α 2tω K e β rtωdb e α 2tω K e β tωd e α tω K. {}}{{}}{{}}{ (Palindromy implies even order and has other welcome features.) 3
4 In one sense the above expression is good: it separates scales. Thus, exponentials with ω and ω are kept apart. In another sense it is bad: we need plenty of exponentials to attain reasonable order. Typically, the order of semidiscretization is very high: it is wasteful to accompany it with very low order of the time-stepping method! And in yet another sense it is ugly: ideally, we would have liked to combine numerics with asymptotics, namely for the arguments of exponentials to become progressively smaller while in the present case they are all of roughly the same order of magnitude. The challenge is to develop a methodology which attains all the advantages without any disadvantages: A high-order splitting that separates scales and produces an asymptotic expansion in increasing powers of h. 4
5 ALGEBRA OF OPERATORS The vector field is a linear combination of two linear operators: 2 x and (multiplication by) V. Let us consider the free Lie algebra F generated by 2 x and V. Note that [V, x] 2 = ( xv 2 ) 2( x V ) x, [[V, x], 2 x] 2 = ( xv 4 ) + 4( xv 3 ) x + 4( xv 2 ) x, 2 [[V, x], 2 V ] = 2( x V ) 2, [[[V, x], 2 x], 2 x] 2 = ( xv 6 ) 6( xv 5 ) x 2( xv 4 ) x 2 8( xv 3 ) x 3 and so on. In general, F g, where g is the Lie algebra g = n k=0 y k (x) k x : n Z +, y 0, y,..., y n smooth & periodic ht n k=0 is the height of an element in g. y k (x) x k = n. 5
6 PROPOSITION ht([x, Y ]) = (ht(x) + ht(y ) ) + for all X, Y g. COROLLARY Any nested commutator in F which has at least two more Vs than 2 xs is necessarily zero. For example, out of the 4 terms in the Hall basis of all the elements of F with at most 7 letters, are all zero. [[[V, x 2 ], V ], V ], [[[[V, 2 x ], V ], V ], V ], [[[[[V, 2 x ], 2 x ], V ], V ], V ], [[[[[V, x 2 ], V ], V ], V ], V ], [[[[V, 2 x ], V ], V ], [V, 2 x ]], [[[[[[V, x 2 ], 2 x ], V ], V ], V ], V ], [[[[[[V, 2 x ], V ], V ], V ], V ], V ], [[[[[V, x 2 ], V ], V ], V ], [V, 2 x ]], [[[[V, 2 x ], V ], V ], [[V, 2 x ], 2 x ]], [[[[V, x 2 ], V ], V ], [[V, 2 x ], V ]] But this is not the real reason why all this is important! The real reason is that V never walks on its own, it is always multiplied by a large parameter ω. The elimination of ω-rich terms allows us to derive asymptotic splittings! 6
7 THE SYMMETRIC ZASSENHAUS SPLITTING We commence from the symmetric Baker Campbell Hausdorff formula where e tx e ty e tx = e sbch(t;x,y ), sbch(t; X, Y ) = t(2x + Y ) 6 t3 [Y, [X, Y ]] + t 5 ( [X, [X, [X, [X, Y ]]]] [Y, [Y, [Y, [X, Y ]]]] 90 [X, [Y, [Y, [X, Y ]]]] + 45 [Y, [X, [X, [X, Y ]]]] + 60 [X, [X, [Y, [X, Y ]]]] + 30 [Y, [Y, [X, [X, Y ]]]]) + O ( t 7). It is possible to expand as far as we wish, noting that all the expansion terms live in the free Lie algebra generated by the letters X and Y. Such algebras possess convenient bases that can be constructed algorithmically, e.g. the Hall basis, the Lyndon basis and the Dynkin basis. 7
8 The classical Zassenhaus splitting is where e t(x+y ) = e tx e ty e t2 U 2 (X,Y ) e t3 U 3 (X,Y ) e t4 U 4 (X,Y ), U 2 (X, Y ) = 2 [X, Y ], U 3 (X, Y ) = 3 [Y, [X, Y ]] + 6 [X, [X, Y ]], U 4 (X, Y ) = 24 [[[X,Y ],X],X] 8 [[[X,Y ],X],Y ] 8 [[[X,Y ],Y ],Y ]. We are, however, interested in symmetric (i.e., palindromic) splittings. Moreover, powers of t are misleading: we need to reckon with three parameters: ω (large!), the number N of degrees of freedom once we semi-discretize (large!) and the time step t (small!). It makes sense to reduce this to a single currency by requiring N O(ω ρ ), t O ( ω σ) for some ρ, σ > 0. Note that x scales like O(N) = O(ω ρ ). In the sequel we assume that ρ = σ = 2. 8
9 SYMMETRIC ZASSENHAUS Setting for simplicity c =, we seek a splitting of the form e i( t)(ω 2 x +ωv ) e R 0e R e Rs e T s+e Rs e R e R 0, {}}{{}}{{}}{ where (recalling, in line with our currency conversion, that x O ( ω /2) ) R k = R k ( t, ω, V ) = O ( ω k+/2), k = 0,,..., s, T s+ = T s+ ( t, ω, V ) = O ( ω s /2). There are two options: we can start with the derivative or the potential. In this talk we ll start with the potential. To this end we let X = 2 τωv, Y = τω 2 x + τωv, where τ = i t = O ( ω /2) and rewrite the sbch formula as e Y = e X e sbch(x,y ) e X. 9
10 In what follows we restrict ourselves to an O ( ω 7/2) expansion, i.e. s = 2. After long calculation: Expanding sbch( 2 τωv, τω 2 x + τω) with the sbch formula and replacing elements from F by terms from g; 2 Throwing away all zero terms, but also all terms which are O ( ω α) for α 7 2 ; 3 Expressing everything in the Hall basis we obtain the asymptotic expansion sbch(x, Y ) = τω C 2 τ 3 ω C 4 24 τ 3 ωc τ 5 ω C τ 5 ωc 80 τ 5 ω C τ 5 ωc τ 7 ωc τ 7 ωc τ 7 ωc τ 7 ωc 40 + O ) (ω 7/2, 0
11 where all that survives of the Hall basis is C = 2 x, C 4 = ( 4 x V ) + 4( 3 x V ) x + 4( 2 x V ) 2 x, C 5 = 2( x V ) 2, C 0 = {6( 5 x V )( xv ) + 2( 4 x V )( 2 x V ) + 8( 3 x V )2 } 24{( 4 x )( xv ) + ( 3 x V )( 2 x V )} x 24( 3 x )( xv ) 2 x, C = 8( 2 x V )( xv ) 2, C 3 = {2( 5 x V )( xv ) 4( 4 x )( 2 x V ) 4( 3 x V )2 } + 8{( 4 x V )( xv ) ( 3 x V )( 2 x V )} x + 4( 3 x V )( xv ) 2 x, C 4 = 4( 2 x V )( xv ) 2, C 27 = 48( 3 x V )( xv ) 3, C 32 = 6( 3 x V )( xv ) ( 2 x V )2 ( x V ) 2, C 35 = 8( 3 x V )( xv ) 3 6( 2 x V )2 ( x V ) 2, C 40 = 32( 2 x V )2 ( x V ) 2. 4 We next express the remaining commutators as terms in the larger Lie algebra g. But here we have a problem!!! Ultimately, we will replace derivatives by (finite-dimensional) differentiation matrices. Even derivatives yield symmetric matrices and multiplication by i results in skew-hermitian matrices good! But, by the same token, odd derivatives skew-symmetric matrices (after multiplication by i) Hermitian matrices bad! We must somehow get rid of odd derivatives!!!
12 5 Replacement of odd derivatives: The main idea is to replace odd derivatives by linear combinations of even derivatives: x [ x ] y(x) x = 2 y(x) 3 x = y (x) 2 x x and so on. The outcome is 0 y(ξ) dξ 2 x 2 y (x) x [ x 0 x 0 y(ξ) dξ 0 y(ξ) dξ 4 x + 4 y (x) 2 2 x[y (x) ] y(ξ) dξ ] sbch(x, Y ) = τω 2 x 2 τ 3 ω { ( 4 x V ) x [( 2 x V ) ] + 2( 2 x V ) 2 x ] + 2 τ 3 ω( x V ) τ 5 ω {6( 5 x V )( xv ) + 2( 4 x V )( 2 x V ) + 4( 3 x V )2 2 2 x [( 3 x V )( xv ) ] 2( 3 x V )( xv ) 2 x } 90 τ 5 ω( 2 x V )( xv ) 2 90 τ 5 ω {3( 5 x V )( xv ) 2( 4 x V )( 2 x V ) 4( 3 x V )2 + 2 x 2 [( 3 x V )( xv ) ] 4 x 2 [( 2 x V )2 ] + 4( x 3 )( xv ) x 2 8( 2 x )2 x 2 } + 72 τ 5 ω( x 2 V )( xv ) τ 7 ω( x 3 V )( xv ) τ 7 ω{2( x 3 V )( xv ) 3 + 4( x 2 V )2 ( x V ) 2 } τ 7 ω{( x 3 V )( xv ) 3 + 2( x 2 V )2 ( x V ) 2 } ) τ 7 ω( x 2 V )2 ( x V ) 2 + O (ω 7/2., 2
13 6 Arranging in decreasing powers of ω: = sbch(x, Y ) ω /2 {}}{ τω 2 x + 2 τ 2 ω( x V ) 2 } ω 3/2 {}}{ 360 τ 5 ω( xv 2 )( x V ) 2 6 τ 3 ω { x[( 2 xv 2 ) ] + ( xv 2 ) x} 2 ω { 5/2 }}{ 2 τ 3 ω ( xv 4 ) ω 5/2 {}}{ 80 τ 5 ω { 5 x[( 2 xv 3 )( x V ) ] + 4 x[( 2 xv 2 ) 2 ] 5( x)( 3 x V ) x 2 + 4( xv 2 ) 2 x} 2 ω { 5/2 }}{ 520 τ 7 ω{09( xv 3 )( x V ) ( xv 2 ) 2 ( x V ) 2 } +O ( ω 7/2). 3
14 We now set R 0 = X old and X new = 2 τω 2 x 24 τ 3 ω( x V ) 2, Y new = sbch(x old, Y old ). We repeat steps 6 and this results in {}}{ sbch(x, Y ) = 360 τ 5 ω( x 2 V )( xv ) 2 6 τ 3 ω { x 2 [( 2 x V ) ] + ( 2 x V ) 2 x } ω 5/2 {}}{ 2 τ 3 ω ( 4 x V ) ω 3/2 ω 5/2 {}}{ 80 τ 5 ω { 5 2 x [( 3 x V )( xv ) ] x [( 2 x V )2 ] 5( 3 x )( xv ) 2 x + 4( 2 x V )2 2 x } ω 5/2 {}}{ ) 520 τ 7 ω{09( x 3 V )( xv ) ( x 2 V )2 ( x V ) 2 } +O (ω 7/2. 4
15 Thus, R = X old and we again set X new as minus half the leading term of Y new = sbch(x old, Y old ). This is the last such step, leading to R 2 and the remainder, T 3. To sum up, R 0 = 2 τωv = O( ω /2), R = 2 τω 2 x + 24 τ 3 ω( x V ) 2 = O ( ω /2), R 2 = 2 τ 3 ω { 2 x[( 2 xv ) ] + ( 2 xv ) 2 x} τ 5 ω( 2 xv )( x V ) 2 = O ( ω 3/2), T 3 = 2 τ 3 ω ( 4 xv ) + 80 τ 5 ω { 5 2 x[( 3 xv )( x V ) ] x[( 2 xv ) 2 ] 5( 3 x)( x V ) 2 x + 4( 2 xv ) 2 2 x} τ 7 ω{09( 3 xv )( x V ) ( 2 xv ) 2 ( x V ) 2 } = O ( ω 5/2). 5
16 Note that, once we replace derivatives by (even-degree!) differentiation matrices, e R 0 is a diagonal matrix, while matrix/vector products with R, R 2 and T 3 can be computed with FFTs as part of a Krylov-subspace-based computation of their exponentials. Stability Suppose that each (2r)th derivative is replaced by the Hermitian differentiation matrix K 2r. Since τ = i t always features with an odd power, K 2r, as well as the diagonal matrix that we obtain from discretizing multiplication operators, are multiplied by ±i and become skew-hermitian. For example R 0 2 ( t)ωid V, R 2 ( t)ω ik 2 24 ( t)3 ωid V 2 x. However... Because of replacement of odd derivatives, we have expressions of the form idh and ihd where D is diagonal and H Hermitian and they are not skew-hermitian! Is this a problem? 6
17 Not at all! Recalling that each term in F is multiplied by ±i, we note that FLA(i 2 x, iv ) su. All our operations, using sbch but also the replacement of odd derivatives, are Lie-algebra compliant (linear combinations and commutators), therefore everything we have obtained in the course of our algorithm stays within su. Thus, as a sanity check, K symmetric, D diagonal i(kd + DK) skew-hermitian, hence R 2 2 ( t)3 ω i(k 2 D V 2 x + D V 2 x K 2 ) ( t)5 ωid Vxx V 2 x T 3 2 τ 3 ω D 4 x V + 80 τ 5 ω {4(K 2 D ( 2 x V ) 2 + D ( 2 x V ) 2 K 2 ) 5(K 2 D ( 3 x V )( x V ) + D ( 3 x V )( xv ) K 2)} τ 7 ω{09d ( 3 x V )( x V ) D ( 2 x V ) 2 ( x V ) 2 } and it is easy to verify that bothr 2 and T 3 are skew-hermitian. THEOREM It is true that R k su, k = 0,,..., s, and T s+ su. 7
18 LEADING WITH THE DERIVATIVE We have commenced the splitting with X = τωv an alternative is to start from X = τω x 2. On the face of it, this makes little sense, because we want to start with the larger O ( ω /2) term. However, let s try... R 0 = 2 τω 2 x = O ( ω /2), R = 2 τωv = O( ω /2), R 2 = 2 τ 3 ω( x V ) 2 = O ( ω /2), R 3 = 24 τ 3 ω { 2 x[( 2 xv ) ] + ( 2 xv ) 2 x} 480 τ 5 ω( 2 xv )( x V ) 2 = O ( ω 3/2), T 4 = 24 τ 3 ω ( xv 4 ) τ 5 ω { 30 2 x[( xv 3 )( x V ) ] + 30 ( 3 xv )( x V ) x x[( xv 2 ) 2 ] 60 ( 2 xv ) 2 x} 2 τ 7 ω{ ( 3 xv )( x V ) ( 2 xv ) 2 ( x V ) 2 } = O ( ω 5/2). On the face of it, more terms. However... R 0 is a circulant and e R 0 can be computed with a single FFT, R and R 2 are diagonal only R 3 = O ( ω 3/2) and T 4 = O ( ω 5/2) require Krylov subspace methods to compute their exponential! 8
19 SPACE DISCRETIZATION In principle, we have three options to discretise 2s x to spectral accuracy: Spectral methods: We project in L 2 the solution on Nth-degree trigonometric polynomials: the unknowns are Fourier coefficients, K is diagonal and D a circulant. Spectral collocation: We interpolate at equally-spaced points by N th-degree trigonometric polynomials: the unknowns are nodal values, K is a circulant and D is diagonal. Pseudospectral method: We wrap around a finite-difference method an infinite number of times: again, the unknowns are nodal values, K is a circulant and D is diagonal. Except that, having committed already an O ( ω 7/2) error, we don t need spectral accuracy! 9
20 FINITE DIFFERENCES We discretise in space to the same O ( ω 7/2) = O ( ( x) 7) (actually, O ( ( x) 8) ) accuracy using finite differences: u m u(m x, ), m = N,..., N, x = N+, and 2 u m ( x) 2( 560 u m u m 3 5 u m u m u m Therefore: u m+ 5 u m u m u m+4).! K is a 9-diagonal circulant, and! D is a diagonal matrix. In that case all matrix-vector products are O(N) operations! 20
21 u(x) = 2 + sin πx 2
22 u(x) = e cos πx 22
23 So, does the error converge to a multiple of some mystery function? Yes: in general, for the rth finite difference of order 2r we have r m= r u (x) ( x) 2 r m= r a m u(x + m x) a m cos mθ = θ 2 + c r θ 2r+2 + O ( θ 2r+4), whereby the error is c r ( x) 2r u (2r+2) (x) + O ( ( x) 2r+2). In our case c 4 = 350, therefore the error is 350 ( x)8 u (0) (x) + O ( ( x) 0), x = N
24 The error (in blue) and the error estimate (in red) for u(x) = (2 + sin πx) and ω =
25 COMPUTING THE EXPONENTIALS Our method produces a dichotomy: Arguments of the exponentials are either trivial or small! They are trivial if the matrix is diagonal or is just a scaled differentiation matrix in both cases computing the exponential is straightforward. They are small if the argument is A = O ( ω α) for some α > 0. In that case the mth term in the Krylov basis {v, Av, A 2 v,, A 3 v...} is O ( ω mα). Therefore, using estimates of Hochbruck, Lubich & Selhofer, it is possible to prove that we need ridiculously small dimension of a Krylov subspace to compute e A v to O ( ω 7/2), say. Leading by derivative, R 3 = O ( ω 3/2) and T 4 = O ( ω 5/2) require dimensions 3 and 2 respectively! 25
26 CHALLENGES Other pairs (ρ, σ). While the generalisation is straightforward (as a concept: the computation is fiendish), the interesting bit is to explore what are good choices of ρ, σ (0, ). Time-dependent potentials. Here we need to marry Magnus expansions with symmetric Zassenhaus. Preliminary results indicate that, although the algebra is formidable, this can be done: the underlying Free Lie algebra is much more complicated but still embedded in g and height restrictions apply. Nonlinear Schrödinger. Can we extend our method to the nonlinear Schrödinger equation i hu t = h2 2m u V (x)u + λ u 2 u? Straightforward generalisation does not work but is there a clever way this set of ideas can be extended? Other equations. Can this work with the wave equation? With nonlinear wave equations, e.g. Klein Gordon? With Hamiltonian ODEs of the form H(p, q) = H (p, q) + H 2 (p, q), where H 2 H? 26
27 Happy birthday Peter!!! 27
Efficient methods for time-dependence in semiclassical Schrödinger equations
Efficient methods for time-dependence in semiclassical Schrödinger equations Philipp Bader, Arieh Iserles, Karolina Kropielnicka & Pranav Singh November 25, 25 Abstract We build an efficient and unitary
More informationOn skew-symmetric differentiation matrices
On skew-symmetric differentiation matrices Arieh Iserles January 9, 20 Abstract The theme of this paper is the construction of finite-difference approximations to the first derivative in the presence of
More information7. Baker-Campbell-Hausdorff formula
7. Baker-Campbell-Hausdorff formula 7.1. Formulation. Let G GL(n,R) be a matrix Lie group and let g = Lie(G). The exponential map is an analytic diffeomorphim of a neighborhood of 0 in g with a neighborhood
More informationRepresentation theory and quantum mechanics tutorial Spin and the hydrogen atom
Representation theory and quantum mechanics tutorial Spin and the hydrogen atom Justin Campbell August 3, 2017 1 Representations of SU 2 and SO 3 (R) 1.1 The following observation is long overdue. Proposition
More informationCommutator-free Magnus Lanczos methods for the linear Schrödinger equation
Commutator-free Magnus Lanczos methods for the linear Schrödinger equation Arieh Iserles, Karolina Kropielnicka & Pranav Singh October 3, 27 AMS Mathematics Subject Classification: Primary 65M7, Secondary
More informationQUADRATIC ELEMENTS IN A CENTRAL SIMPLE ALGEBRA OF DEGREE FOUR
QUADRATIC ELEMENTS IN A CENTRAL SIMPLE ALGEBRA OF DEGREE FOUR MARKUS ROST Contents Introduction 1 1. Preliminaries 2 2. Construction of quadratic elements 2 3. Existence of biquadratic subextensions 4
More informationIntroduction - Motivation. Many phenomena (physical, chemical, biological, etc.) are model by differential equations. f f(x + h) f(x) (x) = lim
Introduction - Motivation Many phenomena (physical, chemical, biological, etc.) are model by differential equations. Recall the definition of the derivative of f(x) f f(x + h) f(x) (x) = lim. h 0 h Its
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
More informationQuantum Theory and Group Representations
Quantum Theory and Group Representations Peter Woit Columbia University LaGuardia Community College, November 1, 2017 Queensborough Community College, November 15, 2017 Peter Woit (Columbia University)
More informationMATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES
MATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES J. WONG (FALL 2017) What did we cover this week? Basic definitions: DEs, linear operators, homogeneous (linear) ODEs. Solution techniques for some classes
More information1 Smooth manifolds and Lie groups
An undergraduate approach to Lie theory Slide 1 Andrew Baker, Glasgow Glasgow, 12/11/1999 1 Smooth manifolds and Lie groups A continuous g : V 1 V 2 with V k R m k open is called smooth if it is infinitely
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
More informationSupplementary Notes March 23, The subgroup Ω for orthogonal groups
The subgroup Ω for orthogonal groups 18.704 Supplementary Notes March 23, 2005 In the case of the linear group, it is shown in the text that P SL(n, F ) (that is, the group SL(n) of determinant one matrices,
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationLinear Algebra, Summer 2011, pt. 2
Linear Algebra, Summer 2, pt. 2 June 8, 2 Contents Inverses. 2 Vector Spaces. 3 2. Examples of vector spaces..................... 3 2.2 The column space......................... 6 2.3 The null space...........................
More informationNumerical Simulation of Spin Dynamics
Numerical Simulation of Spin Dynamics Marie Kubinova MATH 789R: Advanced Numerical Linear Algebra Methods with Applications November 18, 2014 Introduction Discretization in time Computing the subpropagators
More informationAccuracy, Precision and Efficiency in Sparse Grids
John, Information Technology Department, Virginia Tech.... http://people.sc.fsu.edu/ jburkardt/presentations/ sandia 2009.pdf... Computer Science Research Institute, Sandia National Laboratory, 23 July
More informationDefinition 5.1. A vector field v on a manifold M is map M T M such that for all x M, v(x) T x M.
5 Vector fields Last updated: March 12, 2012. 5.1 Definition and general properties We first need to define what a vector field is. Definition 5.1. A vector field v on a manifold M is map M T M such that
More informationQM and Angular Momentum
Chapter 5 QM and Angular Momentum 5. Angular Momentum Operators In your Introductory Quantum Mechanics (QM) course you learned about the basic properties of low spin systems. Here we want to review that
More informationPart 3.3 Differentiation Taylor Polynomials
Part 3.3 Differentiation 3..3.1 Taylor Polynomials Definition 3.3.1 Taylor 1715 and Maclaurin 1742) If a is a fixed number, and f is a function whose first n derivatives exist at a then the Taylor polynomial
More informationGenerators for Continuous Coordinate Transformations
Page 636 Lecture 37: Coordinate Transformations: Continuous Passive Coordinate Transformations Active Coordinate Transformations Date Revised: 2009/01/28 Date Given: 2009/01/26 Generators for Continuous
More informationNumber Systems III MA1S1. Tristan McLoughlin. December 4, 2013
Number Systems III MA1S1 Tristan McLoughlin December 4, 2013 http://en.wikipedia.org/wiki/binary numeral system http://accu.org/index.php/articles/1558 http://www.binaryconvert.com http://en.wikipedia.org/wiki/ascii
More informationAn Introduction to Partial Differential Equations
An Introduction to Partial Differential Equations Ryan C. Trinity University Partial Differential Equations Lecture 1 Ordinary differential equations (ODEs) These are equations of the form where: F(x,y,y,y,y,...)
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND First Printing, 99 Chapter LINEAR EQUATIONS Introduction to linear equations A linear equation in n unknowns x,
More informationLinear maps. Matthew Macauley. Department of Mathematical Sciences Clemson University Math 8530, Spring 2017
Linear maps Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 8530, Spring 2017 M. Macauley (Clemson) Linear maps Math 8530, Spring 2017
More informationEigenvalues and Eigenvectors
/88 Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Eigenvalue Problem /88 Eigenvalue Equation By definition, the eigenvalue equation for matrix
More information2014 Summer Review for Students Entering Algebra 2. TI-84 Plus Graphing Calculator is required for this course.
1. Solving Linear Equations 2. Solving Linear Systems of Equations 3. Multiplying Polynomials and Solving Quadratics 4. Writing the Equation of a Line 5. Laws of Exponents and Scientific Notation 6. Solving
More informationCompact schemes for laser-matter interaction in Schrödinger equation
Compact schemes for laser-matter interaction in Schrödinger equation Arieh Iserles, Karolina Kropielnicka & Pranav Singh April 21, 218 AMS Mathematics Subject Classification: Primary 65M7, Secondary 35Q41,
More informationPoisson Brackets and Lie Operators
Poisson Brackets and Lie Operators T. Satogata January 22, 2008 1 Symplecticity and Poisson Brackets 1.1 Symplecticity Consider an n-dimensional 2n-dimensional phase space) linear system. Let the canonical
More informationLecture for Week 2 (Secs. 1.3 and ) Functions and Limits
Lecture for Week 2 (Secs. 1.3 and 2.2 2.3) Functions and Limits 1 First let s review what a function is. (See Sec. 1 of Review and Preview.) The best way to think of a function is as an imaginary machine,
More informationDesigning Information Devices and Systems I Fall 2018 Lecture Notes Note Positioning Sytems: Trilateration and Correlation
EECS 6A Designing Information Devices and Systems I Fall 08 Lecture Notes Note. Positioning Sytems: Trilateration and Correlation In this note, we ll introduce two concepts that are critical in our positioning
More informationDesigning Information Devices and Systems I Fall 2018 Lecture Notes Note Positioning Sytems: Trilateration and Correlation
EECS 6A Designing Information Devices and Systems I Fall 08 Lecture Notes Note. Positioning Sytems: Trilateration and Correlation In this note, we ll introduce two concepts that are critical in our positioning
More informationThe Spinor Representation
The Spinor Representation Math G4344, Spring 2012 As we have seen, the groups Spin(n) have a representation on R n given by identifying v R n as an element of the Clifford algebra C(n) and having g Spin(n)
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 9, February 8, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 9, February 8, 2006 The Harmonic Oscillator Consider a diatomic molecule. Such a molecule
More informationHints and Selected Answers to Homework Sets
Hints and Selected Answers to Homework Sets Phil R. Smith, Ph.D. March 8, 2 Test I: Systems of Linear Equations; Matrices Lesson. Here s a linear equation in terms of a, b, and c: a 5b + 7c 2 ( and here
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 informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More information[Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty.]
Math 43 Review Notes [Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty Dot Product If v (v, v, v 3 and w (w, w, w 3, then the
More informationAlgebra. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.
This document was written and copyrighted by Paul Dawkins. Use of this document and its online version is governed by the Terms and Conditions of Use located at. The online version of this document is
More informationNEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING Fernando Casas Universitat Jaume I, Castellón, Spain Fernando.Casas@uji.es (on sabbatical leave at DAMTP, University of Cambridge) Edinburgh,
More informationStatistics 3657 : Moment Generating Functions
Statistics 3657 : Moment Generating Functions A useful tool for studying sums of independent random variables is generating functions. course we consider moment generating functions. In this Definition
More informationInverse Variation. y varies inversely as x. REMEMBER: Direct variation y = kx where k is not equal to 0.
Inverse Variation y varies inversely as x. REMEMBER: Direct variation y = kx where k is not equal to 0. Inverse variation xy = k or y = k where k is not equal to 0. x Identify whether the following functions
More informationFourier series on triangles and tetrahedra
Fourier series on triangles and tetrahedra Daan Huybrechs K.U.Leuven Cambridge, Sep 13, 2010 Outline Introduction The unit interval Triangles, tetrahedra and beyond Outline Introduction The topic of this
More informationz x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.
Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These
More informationCS205B / CME306 Homework 3. R n+1 = R n + tω R n. (b) Show that the updated rotation matrix computed from this update is not orthogonal.
CS205B / CME306 Homework 3 Rotation Matrices 1. The ODE that describes rigid body evolution is given by R = ω R. (a) Write down the forward Euler update for this equation. R n+1 = R n + tω R n (b) Show
More informationSCHUR-WEYL DUALITY FOR U(n)
SCHUR-WEYL DUALITY FOR U(n) EVAN JENKINS Abstract. These are notes from a lecture given in MATH 26700, Introduction to Representation Theory of Finite Groups, at the University of Chicago in December 2009.
More informationPolynomials and Polynomial Equations
Polynomials and Polynomial Equations A Polynomial is any expression that has constants, variables and exponents, and can be combined using addition, subtraction, multiplication and division, but: no division
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationRelevant sections from AMATH 351 Course Notes (Wainwright): Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls):
Lecture 5 Series solutions to DEs Relevant sections from AMATH 35 Course Notes (Wainwright):.4. Relevant sections from AMATH 35 Course Notes (Poulin and Ingalls): 2.-2.3 As mentioned earlier in this course,
More informationPRELIMINARY THEORY LINEAR EQUATIONS
4.1 PRELIMINARY THEORY LINEAR EQUATIONS 117 4.1 PRELIMINARY THEORY LINEAR EQUATIONS REVIEW MATERIAL Reread the Remarks at the end of Section 1.1 Section 2.3 (especially page 57) INTRODUCTION In Chapter
More informationINSTANTON MODULI AND COMPACTIFICATION MATTHEW MAHOWALD
INSTANTON MODULI AND COMPACTIFICATION MATTHEW MAHOWALD () Instanton (definition) (2) ADHM construction (3) Compactification. Instantons.. Notation. Throughout this talk, we will use the following notation:
More informationCOLLEGE ALGEBRA. Paul Dawkins
COLLEGE ALGEBRA Paul Dawkins Table of Contents Preface... iii Outline... iv Preliminaries... 7 Introduction... 7 Integer Exponents... 8 Rational Exponents...5 Radicals... Polynomials...30 Factoring Polynomials...36
More informationπ(x Y ) and π(k); Semi-Direct Products
π(x Y ) and π(k); Semi-Direct Products Richard Koch January 29, 2006 1 Group Products Recall that if H and K are groups, the new group H K is the set of all ordered pairs (h, k) with h H and k K. The group
More informationPartial Differential Equations Separation of Variables. 1 Partial Differential Equations and Operators
PDE-SEP-HEAT-1 Partial Differential Equations Separation of Variables 1 Partial Differential Equations and Operators et C = C(R 2 ) be the collection of infinitely differentiable functions from the plane
More informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationLinear Algebra Notes. Lecture Notes, University of Toronto, Fall 2016
Linear Algebra Notes Lecture Notes, University of Toronto, Fall 2016 (Ctd ) 11 Isomorphisms 1 Linear maps Definition 11 An invertible linear map T : V W is called a linear isomorphism from V to W Etymology:
More informationCombinatorial Structures
Combinatorial Structures Contents 1 Permutations 1 Partitions.1 Ferrers diagrams....................................... Skew diagrams........................................ Dominance order......................................
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More information1 Matrices and vector spaces
Matrices and vector spaces. Which of the following statements about linear vector spaces are true? Where a statement is false, give a counter-example to demonstrate this. (a) Non-singular N N matrices
More informationConsistent Histories. Chapter Chain Operators and Weights
Chapter 10 Consistent Histories 10.1 Chain Operators and Weights The previous chapter showed how the Born rule can be used to assign probabilities to a sample space of histories based upon an initial state
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More informationLast Update: April 7, 201 0
M ath E S W inter Last Update: April 7, Introduction to Partial Differential Equations Disclaimer: his lecture note tries to provide an alternative approach to the material in Sections.. 5 in the textbook.
More informationHomework Solutions: , plus Substitutions
Homework Solutions: 2.-2.2, plus Substitutions Section 2. I have not included any drawings/direction fields. We can see them using Maple or by hand, so we ll be focusing on getting the analytic solutions
More informationPartial Differential Equations
Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives with respect to those variables. Most (but
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More information1 Fields and vector spaces
1 Fields and vector spaces In this section we revise some algebraic preliminaries and establish notation. 1.1 Division rings and fields A division ring, or skew field, is a structure F with two binary
More informationBASICS OF ALGEBRA M.K. HOME TUITION. Mathematics Revision Guides. Level: GCSE Foundation Tier
Mathematics Revision Guides Basics of Algebra Page 1 of 7 M.K. HOME TUITION Mathematics Revision Guides Level: GCSE Foundation Tier BASICS OF ALGEBRA Version: 3.1 Date: 01-12-2013 Mathematics Revision
More informationA polynomial expression is the addition or subtraction of many algebraic terms with positive integer powers.
LEAVING CERT Honours Maths notes on Algebra. A polynomial expression is the addition or subtraction of many algebraic terms with positive integer powers. The degree is the highest power of x. 3x 2 + 2x
More informationa 11 x 1 + a 12 x a 1n x n = b 1 a 21 x 1 + a 22 x a 2n x n = b 2.
Chapter 1 LINEAR EQUATIONS 11 Introduction to linear equations A linear equation in n unknowns x 1, x,, x n is an equation of the form a 1 x 1 + a x + + a n x n = b, where a 1, a,, a n, b are given real
More informationMATH 583A REVIEW SESSION #1
MATH 583A REVIEW SESSION #1 BOJAN DURICKOVIC 1. Vector Spaces Very quick review of the basic linear algebra concepts (see any linear algebra textbook): (finite dimensional) vector space (or linear space),
More informationB 1 = {B(x, r) x = (x 1, x 2 ) H, 0 < r < x 2 }. (a) Show that B = B 1 B 2 is a basis for a topology on X.
Math 6342/7350: Topology and Geometry Sample Preliminary Exam Questions 1. For each of the following topological spaces X i, determine whether X i and X i X i are homeomorphic. (a) X 1 = [0, 1] (b) X 2
More informationUNIFORM BOUNDS FOR BESSEL FUNCTIONS
Journal of Applied Analysis Vol. 1, No. 1 (006), pp. 83 91 UNIFORM BOUNDS FOR BESSEL FUNCTIONS I. KRASIKOV Received October 8, 001 and, in revised form, July 6, 004 Abstract. For ν > 1/ and x real we shall
More informationQuantum Mechanics - I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras. Lecture - 16 The Quantum Beam Splitter
Quantum Mechanics - I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 16 The Quantum Beam Splitter (Refer Slide Time: 00:07) In an earlier lecture, I had
More informationIn this chapter we study several functions that are useful in calculus and other areas of mathematics.
Calculus 5 7 Special functions In this chapter we study several functions that are useful in calculus and other areas of mathematics. 7. Hyperbolic trigonometric functions The functions we study in this
More information6.5 Trigonometric Equations
6. Trigonometric Equations In this section, we discuss conditional trigonometric equations, that is, equations involving trigonometric functions that are satisfied only by some values of the variable (or
More informationHandout 2: Invariant Sets and Stability
Engineering Tripos Part IIB Nonlinear Systems and Control Module 4F2 1 Invariant Sets Handout 2: Invariant Sets and Stability Consider again the autonomous dynamical system ẋ = f(x), x() = x (1) with state
More informationComputational Linear Algebra
Computational Linear Algebra PD Dr. rer. nat. habil. Ralf-Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2018/19 Part 4: Iterative Methods PD
More informationSurprising Sinc Sums and Integrals
David Borwein University of Western Ontario Peter Borwein Conference May 12-16, 28 1 Motivation and preliminaries. This talk is based on material in a paper to appear shortly in MAA MONTHLY with the above
More informationINTRODUCTION TO ALGEBRAIC GEOMETRY
INTRODUCTION TO ALGEBRAIC GEOMETRY WEI-PING LI 1 Preliminary of Calculus on Manifolds 11 Tangent Vectors What are tangent vectors we encounter in Calculus? (1) Given a parametrised curve α(t) = ( x(t),
More informationES.1803 Topic 7 Notes Jeremy Orloff. 7 Solving linear DEs; Method of undetermined coefficients
ES.1803 Topic 7 Notes Jeremy Orloff 7 Solving linear DEs; Method of undetermined coefficients 7.1 Goals 1. Be able to solve a linear differential equation by superpositioning a particular solution with
More informationRings If R is a commutative ring, a zero divisor is a nonzero element x such that xy = 0 for some nonzero element y R.
Rings 10-26-2008 A ring is an abelian group R with binary operation + ( addition ), together with a second binary operation ( multiplication ). Multiplication must be associative, and must distribute over
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 informationChapter 2 The Group U(1) and its Representations
Chapter 2 The Group U(1) and its Representations The simplest example of a Lie group is the group of rotations of the plane, with elements parametrized by a single number, the angle of rotation θ. It is
More informationDifferential of the Exponential Map
Differential of the Exponential Map Ethan Eade May 20, 207 Introduction This document computes ɛ0 log x + ɛ x ɛ where and log are the onential mapping and its inverse in a Lie group, and x and ɛ are elements
More information2 Systems of Linear Equations
2 Systems of Linear Equations A system of equations of the form or is called a system of linear equations. x + 2y = 7 2x y = 4 5p 6q + r = 4 2p + 3q 5r = 7 6p q + 4r = 2 Definition. An equation involving
More informationPascal s Triangle on a Budget. Accuracy, Precision and Efficiency in Sparse Grids
Covering : Accuracy, Precision and Efficiency in Sparse Grids https://people.sc.fsu.edu/ jburkardt/presentations/auburn 2009.pdf... John Interdisciplinary Center for Applied Mathematics & Information Technology
More informationQuadratic reciprocity (after Weil) 1. Standard set-up and Poisson summation
(September 17, 010) Quadratic reciprocity (after Weil) Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ I show that over global fields (characteristic not ) the quadratic norm residue
More informationLinear Algebra, Summer 2011, pt. 3
Linear Algebra, Summer 011, pt. 3 September 0, 011 Contents 1 Orthogonality. 1 1.1 The length of a vector....................... 1. Orthogonal vectors......................... 3 1.3 Orthogonal Subspaces.......................
More informationQuantum Field Theory III
Quantum Field Theory III Prof. Erick Weinberg March 9, 0 Lecture 5 Let s say something about SO(0. We know that in SU(5 the standard model fits into 5 0(. In SO(0 we know that it contains SU(5, in two
More informationLie Groups for 2D and 3D Transformations
Lie Groups for 2D and 3D Transformations Ethan Eade Updated May 20, 2017 * 1 Introduction This document derives useful formulae for working with the Lie groups that represent transformations in 2D and
More informationMath Exam 2, October 14, 2008
Math 96 - Exam 2, October 4, 28 Name: Problem (5 points Find all solutions to the following system of linear equations, check your work: x + x 2 x 3 2x 2 2x 3 2 x x 2 + x 3 2 Solution Let s perform Gaussian
More information1 GSW Sets of Systems
1 Often, we have to solve a whole series of sets of simultaneous equations of the form y Ax, all of which have the same matrix A, but each of which has a different known vector y, and a different unknown
More informationNotes on SU(3) and the Quark Model
Notes on SU() and the Quark Model Contents. SU() and the Quark Model. Raising and Lowering Operators: The Weight Diagram 4.. Triangular Weight Diagrams (I) 6.. Triangular Weight Diagrams (II) 8.. Hexagonal
More informationQUIVERS, REPRESENTATIONS, ROOTS AND LIE ALGEBRAS. Volodymyr Mazorchuk. (Uppsala University)
QUIVERS, REPRESENTATIONS, ROOTS AND LIE ALGEBRAS Volodymyr Mazorchuk (Uppsala University) . QUIVERS Definition: A quiver is a quadruple Q = (V, A, t, h), where V is a non-empty set; A is a set; t and h
More informationThe double layer potential
The double layer potential In this project, our goal is to explain how the Dirichlet problem for a linear elliptic partial differential equation can be converted into an integral equation by representing
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 information(x 1 +x 2 )(x 1 x 2 )+(x 2 +x 3 )(x 2 x 3 )+(x 3 +x 1 )(x 3 x 1 ).
CMPSCI611: Verifying Polynomial Identities Lecture 13 Here is a problem that has a polynomial-time randomized solution, but so far no poly-time deterministic solution. Let F be any field and let Q(x 1,...,
More informationAssignment 10. Arfken Show that Stirling s formula is an asymptotic expansion. The remainder term is. B 2n 2n(2n 1) x1 2n.
Assignment Arfken 5.. Show that Stirling s formula is an asymptotic expansion. The remainder term is R N (x nn+ for some N. The condition for an asymptotic series, lim x xn R N lim x nn+ B n n(n x n B
More information