Reduced models and numerical analysis in molecular quantum dynamics I. Variational approximation
|
|
- Imogen Gallagher
- 5 years ago
- Views:
Transcription
1 Reduced models and numerical analysis in molecular quantum dynamics I. Variational approximation Christian Lubich Univ. Tübingen Bressanone/Brixen, 14 February 2011
2 Computational challenges and buzzwords multiscale high-dimensional highly oscillatory structure preservation interdisciplinary large-scale... all that and more in computational quantum dynamics
3 Programme Variational approximation: Monday Multiconfiguration methods: Tuesday Semiclassical wavepackets: Thursday
4 Reference for this short course EMS, Zurich, 2008
5 Setting Full quantum dynamics: Schrödinger equation i ψ = Hψ, ψ = ψ(x 1,..., x N, t) Intermediate models: via Dirac Frenkel variational principle Born Oppenheimer, MCTDH, Gaussian wavepackets,... Classical molecular dynamics: Newtonian eqs. of motion Mẍ = V (x), x = (x 1 (t),..., x N (t))
6 Outline The Dirac-Frenkel variational approximation principle Abstract formulation Examples: Born Oppenheimer, Hartree, Gaussians Structural properties Approximation properties Variational splitting integrators Preparation: Strang splitting for TDSE Variational splitting: abstract formulation Variational splitting for Gaussian wavepacket dynamics
7 Outline The Dirac-Frenkel variational approximation principle Abstract formulation Examples: Born Oppenheimer, Hartree, Gaussians Structural properties Approximation properties Variational splitting integrators Preparation: Strang splitting for TDSE Variational splitting: abstract formulation Variational splitting for Gaussian wavepacket dynamics
8 Dirac-Frenkel principle: abstract formulation i ψ = Hψ on H = L 2 (R 3N ) M H approximation manifold T u M tangent space at u M, required to be complex linear approximate wavefunction u = u(, t) M determined from i u Hu, v = 0 for all v T u M Galerkin method on a state-dependent approximation space Dirac 1930, Frenkel 1934 (TDHF)
9 Examples Born Oppenheimer Hartree Gaussian wavepackets
10 Molecular Schrödinger equation N nuclei, L electrons x 1,..., x N y 1,..., y L x n, y l R 3 wavefunction Ψ = Ψ(x 1,..., x N, y 1,..., y L, t) i Ψ = HΨ with Hamiltonian H = T + U: T = T N + T e = N n=1 2 2M n xn L l=1 2 2m y l kinetic energy U = U ee + U en + U NN, U ee = j<l e 2 y j y l,... potential CO 2 : N + L = = 25 (x, y) R 75
11 Adiabatic approximation (Born Oppenheimer) H = T N + H e with H e = T e + U, U = U(x, y) H e = H e (x) on L 2 y parametric dependence on nuclear coord. x electronic eigenvalue problem: H e (x)φ(x, ) = E(x)Φ(x, ) variational approximation on M = {u : u(x, y) = ψ(x)φ(x, y), ψ L 2 x} reduced Schrödinger equation for the nuclei: i ψ = (T N + E +... )ψ, ψ = ψ(x 1,..., x N, t) get rid of the electrons (x 1,..., x N, y 1,..., y L )
12 Hartree approximation nuclear SE i ψ = (T + V )ψ, ψ = ψ(x 1,..., x N, t) T = n T n = n V = V (x 1,..., x N ) E variational approximation on 2 2M n xn M = {u : u(x) = φ 1 (x 1 )... φ N (x N )} nonlinear system of low-dimensional SE for φ n : i φ n = (T n + V n )φ n, φ n = φ n (x n, t) (n = 1,..., N) V n (x n ) = k n φ k V k n φ k x n mean field potential separate the nuclei TDH, TDSCF
13 Gaussian wavepackets nuclear SE iɛ ψ = (T + V )ψ, ψ = ψ(x 1,..., x N, t) T = n ɛ2 2M n xn, V = V (x 1,..., x N ) variational approximation by complex Gaussians: M = { u : u(x) = N n=1 ( i exp an x n q n ɛ( 2 ) )} + p n (x n q n ) + c n parameters q n, p n R 3, a n = α n + iβ n C, c n = γ n + iδ n C note q n = u x n u position, p n = u iɛ n u momentum motivation: Gaussian wavepackets exact for quadratic potentials parametrize the wavefunction Heller 1976
14 GWP equations of motion ODEs for parameters q, p, a = α + iβ, c = γ + iδ: q = p/m ṗ = V ȧ = 2a2 m 1 6 V ċ = p 2 2m V + 3ɛia m + ɛ 8β V with V = u V u = V (x) exp ( 2β ɛ x q 2 2δ ) ɛ dx, etc. note: for ɛ 0 obtain classical dynamics q = p/m, ṗ = V (q) semiclassical dynamics Heller 1976
15 Structural properties of variational approximation i u Hu, v = 0 for all v T u M
16 Orthogonal projection Minimum defect condition: determine approximation t u(t) M such that equivalently: u = ϑ T u M with ϑ 1 Hu = min! i Re u 1 i Hu, v = 0 for all v T um or u = P(u) 1 i Hu : orthogonal projection Frenkel 1934 Dirac
17 Cover illustration: tangent space projection u = P(u) 1 i Hu
18 Symplectic projection Im u 1 i Hu, v = 0 for all v T um symplectic 2-form on L 2 (R 3N ): ω(ξ, η) = 2 Im ξ, η Hamilton function: H(u) = Hu, u Hamiltonian system on the manifold M: Consequences: ω( u, v) = dh(u)v for all v T u M energy conservation: H(u(t)) = Const. symplecticity of the flow: for ξ(t), η(t) tangent vectors propagated by the linearized flow along u(t), ω(ξ(t), η(t)) = Const.
19 Action principle action functional S(u) = t1 t 0 i u(t) Hu(t), u(t) dt made stationary with respect to variations of paths on the manifold M with fixed end-points: t1 0 = δs(u) = 2 Im u(t) 1 i Hu(t), δu(t) dt. t 0 Dirac s quantum-mechanical analogue of Hamilton s principle in classical mechanics.
20 Poisson structure coordinates y on M: u = χ(y) total energy: K(y) = u H u structure matrix: B(y) = ( 2 Im χ (y) χ (y) ) 1 equations of motion in coordinates y: ẏ = B(y) K(y) Poisson system: bracket {F, G} = F T B G is skew-symmetric and satisfies the Jacobi and Leibniz identities Poisson flow ϕ t : {F ϕ t, G ϕ t } = {F, G} ϕ t F, G non-canonical Hamiltonian system
21 Example: Gauss as Poisson Gaussian wavepacket: y = (p, q, α, β, γ, δ) Poisson system ẏ = B(y) K(y) with p B(y) = 1 4β ɛ 0 β N(y) 0 0 4β2 3ɛ 0 β 0 p T β 0 4 ɛ 0 0 β ɛ 0 K(y) = u H u, N(y) = u 2 Both K and N are conserved quantities
22 Approximation properties u(t) ψ(t)?
23 Approximation error A posteriori error bound u(t) ψ(t) Quasi-optimality t 0 dist ( 1 i Hu(τ), T u(τ)m ) dτ t u(t) ψ(t) d(t) + Ce ct d(τ) dτ with d(t) = dist(ψ(t), M) the best-approximation error, under assumptions that include Hartree, Gaussians,... 0 L. 2005
24 Example: Approximation by Gaussians u(t) ψ(t) t 0 dist ( 1 iε Hu(τ), T u(τ)m ) dτ local harmonic approximation for Gaussian u centred at q H = T + V with V = U q + W q note Tu T u M and U q u T u M, therefore dist ( 1 iε Hu, T um ) = dist ( 1 iε W qu, T u M ) = O(ε 1/2 ) and hence u(t) ψ(t) C t ε 1/2 cf. Hagedorn 1980
25 Outline The Dirac-Frenkel variational approximation principle Abstract formulation Examples: Born Oppenheimer, Hartree, Gaussians Structural properties Approximation properties Variational splitting integrators Preparation: Strang splitting for TDSE Variational splitting: abstract formulation Variational splitting for Gaussian wavepacket dynamics
26 Strang splitting for the Schrödinger equation i ψ = Hψ, H = T + V (low-dimensional) TDSE ψ(t + t) = exp( i th)ψ(t) ψ n+1 = exp( i 2 tv ) exp( i tt ) exp( i 2 tv ) ψ n ( ) exp( itv )ψ (x) = e itv (x) ψ(x) ( ) F exp( itt )ψ (ξ) = e it ξ 2 (Fψ)(ξ) space discretization FFT Splitting on sparse grid in higher dimension: Gradinaru 2007
27 Error bounds for the Strang splitting without space discretization, for bounded smooth potential V : ψ n ψ(t n ) C t 2 t n max ψ(t) H 2 cf. Crank-Nicolson scheme ψ n+1 = (1 + i 2 th) 1 (1 i 2 th)ψ n : ψ n ψ(t n ) C t 2 t n max d 3 ψ L dt 3 2 Jahnke & L. 2000
28 Variational splitting: the formal game variational approximation: H = T + V i u Hu, v = 0 for all v T u M time step u n u n+1 via 1. half-step with V : u n+1/2 solution at t/2 of i u V u, v = 0 for all v T u M 2. full step with T : u + n+1/2 solution at t of i u T u, v = 0 for all v T u M 3. half-step with V : u n+1 solution at t/2 of i u V u, v = 0 for all v T u M
29 Variational splitting: examples Gaussian wavepackets MCTDH: tomorrow
30 Variational splitting for Gaussian wavepacket dynamics by magic: Substeps for T and V can be solved exactly and explicitly! Faou & L. 2006
31 Gaussian wavepacket integrator (1) First half-step with V : p n+1/2 = p n t 2 V n α + n = α n t 12 V n γ + n = γ n + t ɛ 16β n V n
32 Gaussian wavepacket integrator (2) Step with T : From p n+1/2, a n + = α n + + iβ n, c n + = γ n + + iδ n get q n+1, an+1 = α n+1 + iβ n+1, cn+1 = γ n+1 + iδ n+1: q n+1 = q n + t m p n+1/2 a n+1 = a + n /( ) t m a+ n c n+1 = c + n iɛ log (1 + 2 t m a+ n )
33 Gaussian wavepacket integrator (3) Second half-step with V : p n+1 = p n+1/2 t 2 V n+1 α n+1 = α n+1 t 12 V n+1 γ n+1 = γ n+1 + t ɛ 16β n+1 V n+1 explicit step, because V n+1 and V n+1 depend only on q n+1, β n+1, δ n+1 determined before
34 Properties of Gaussian wavepacket integrator (1) explicit integrator of order 2 time-reversible integrator Poisson/symplectic integrator: y n+1 = Φ t (y n ) with {F Φ t, G Φ t } = {F, G} Φ t F, G (as a composition of exact Poisson flows) L 2 norm of wavepacket is conserved exactly linear and angular momentum are conserved exactly in systems with translational/rotational symmetry
35 Properties of Gaussian wavepacket integrator (2) long-time near-conservation of energy: K(y) = u H u is conserved up to O( t 2 ) over exponentially long times t e c/ t uniformly in ɛ (use backward error analysis: numerical map is Poisson and therefore the almost-exact flow of a modified Poisson system) classical limit ɛ 0: yields Verlet method p n+1/2 = p n t 2 V (q n) q n+1 = q n + t m p n+1/2 p n+1 = p n+1/2 t 2 V (q n+1) standard integrator of classical MD
36 Outline The Dirac-Frenkel variational approximation principle Abstract formulation Examples: Born Oppenheimer, Hartree, Gaussians Structural properties Approximation properties Variational splitting integrators Preparation: Strang splitting for TDSE Variational splitting: abstract formulation Variational splitting for Gaussian wavepacket dynamics
37
38 2 x 10 3 Error x 10 3 Error x 10 3 Error T x 10 4 Hénon-Heiles, t = 0.1, ε = 10 2, 10 3, 10 4
39 ε = 0.01 Verlet ε = q2 2 3 ε = q1 3D Kepler, t = 0.1, classical limit toward the Verlet integrator
40 Reduced models and numerical analysis in molecular quantum dynamics II. Multiconfiguration methods Christian Lubich Univ. Tübingen Bressanone/Brixen, 15 February 2011
41 Outline MCTDH MCTDHF Variational splitting integrator Modelling error of MCTDH
42 Outline MCTDH MCTDHF Variational splitting integrator Modelling error of MCTDH
43 Recap: Dirac-Frenkel principle, abstract formulation i ψ = Hψ on H = L 2 (R 3N ) M H approximation manifold T u M tangent space at u M, required to be complex linear approximate wavefunction u = u(, t) M determined from i u Hu, v = 0 for all v T u M Galerkin method on a state-dependent approximation space
44 Recap: Time-dependent Hartree method i ψ = (T + V )ψ, ψ = ψ(x 1,..., x N, t) T = n T n = n 2 2M n xn variational approximation on manifold of Hartree products M = {u : u(x) = φ 1 (x 1 )... φ N (x N )} nonlinear system of low-dimensional SE for φ n : i φ n = (T n + V n )φ n, φ n = φ n (x n, t) (n = 1,..., N) V n (x n ) = k n φ k V k n φ k x n mean field potential separation of variables TDH
45 MCTDH Multi-Configuration Time-Dependent Hartree method approximate by linear combination of Hartree products ψ(x, t) r 1 j 1 =1... r N j N =1 a j1...j N (t) φ (1) j 1 (x 1, t)... φ (N) j N (x N, t) mutually orthogonal single-particle functions φ (n) 1,..., φ(n) r n core tensor (a j1...j N ) of full multilinear rank (r 1,..., r N ) H.D. Meyer et al , MCTDH book 2009
46 MCTDH equations of motion from Dirac Frenkel variational principle on MCTDH manifold: coupled system of ODEs and low-dimensional nonlinear PDEs i da J dt = L Φ J H Φ L a L (n) φ(n) iρ t = (I P (n) ) H xn φ (n) with Φ J (x, t) = N n=1 φ(n) j n (x n, t) for multi-indices J = (j 1,..., j N ) Meyer, Manthe & Cederbaum 1990
47 Density matrices ρ (n) (t) nth unfolding of coefficient tensor (a j1...j N ) (for n = 1,..., N): A (n) C rn k n r k set ρ (n) = A (n) A (n) Crn rn hermitian, positive semi-definite ρ (n) positive definite A (n) has full rank (a j1...j N ) has full n-rank
48 MCTDH equations of motion from Dirac Frenkel variational principle on MCTDH manifold: coupled system of ODEs and low-dimensional nonlinear PDEs i da J dt = L Φ J H Φ L a L (n) φ(n) iρ t = (I P (n) ) H xn φ (n) with Φ J (x, t) = N n=1 φ(n) j n (x n, t) for multi-indices J = (j 1,..., j N ) Meyer, Manthe & Cederbaum 1990
49 Existence and regularity assume smooth and bounded potential: V BC 2 initial data in the Sobolev space H 2 There exists a unique strong solution of the MCTDH equations a J C 2 ([0, t ), C), φ (n) j n C 1 ([0, t ), L 2 ) C([0, t ), H 2 ), where either t = or ρ (n) (t) becomes singular as t t. The solution depends on the initial data Lipschitz continuously with respect to the H 2 norm on every compact interval. Moreover, u defined by a J, φ (n) j n indeed solves the Dirac Frenkel variational equation. Koch & L. 2007
50 Outline MCTDH MCTDHF Variational splitting integrator Modelling error of MCTDH
51 Time-dependent Hartree-Fock method fermions: antisymmetric wavefunction ψ(x 1,..., x N, t) consider Coulomb interaction V = j<k 1 x j x k variational approximation on manifold of Slater determinants M = {u : u = φ 1 φ N 1 N! det ( φ j (x k ) ) } nonlinear system of low-dimensional SE for orbitals φ j : i φ j = (T + V j )φ j, V j mean field potential TDHF, Dirac 1930
52 MCTDHF Multi-Configuration Time-Dependent Hartree-Fock method approximate by linear combination of Slater determinants ψ(, t) a j1...j N (t) φ j1 (, t)... φ jn (, t) 1 j 1 <...<j N K mutually orthogonal orbitals φ 1,..., φ K core tensor (a j1...j N ) of full multilinear rank equivalent to MCTDH with antisymmetric core tensor
53 MCTDHF equations of motion from Dirac Frenkel variational principle on MCTDHF manifold: coupled system of ODEs and low-dimensional nonlinear PDEs i da J dt iρ φ t = L Φ J H Φ L a L = (I P) H MF φ with Φ J = N n=1 φ j n for multi-indices J = (j 1 < < j N ) like MCTDH equations, simplified
54 Existence and regularity assume Coulomb interaction initial data in H 2 There exists a unique strong solution of the MCTDHF equations a J C 2 ([0, t ), C), φ j C 1 ([0, t ), L 2 ) C([0, t ), H 2 ), where either t = or ρ (n) (t) becomes singular as t t. Moreover, u defined by a J, φ j indeed solves the Dirac Frenkel variational equation. Bardos, Catto, Mauser & Trabelsi 2010 Koch & L. 2010
55 Outline MCTDH MCTDHF Variational splitting integrator Modelling error of MCTDH
56 Difficulties with standard time integration methods explicit methods: (e.g., explicit Runge-Kutta) step size restriction t c x 2 for stability implicit methods: (e.g., Crank-Nicolson) computationally too expensive, infeasible aim for unconditionally stable explicit method
57 Variational splitting: the formal game variational approximation: H = T + V i u Hu, v = 0 for all v T u M time step u n u n+1 via 1. half-step with V : u n+1/2 solution at t/2 of i u V u, v = 0 for all v T u M 2. full step with T : u + n+1/2 solution at t of i u T u, v = 0 for all v T u M 3. half-step with V : u n+1 solution at t/2 of i u V u, v = 0 for all v T u M
58 Variational splitting for MCTDH and MCTDHF Step with T : since Tu T u M for u M, integration step with T decouples into ȧ J = 0 and i φ(n) j (x n, t) = 2 xn φ (n) j (x n, t) t 2M n low-dimensional free Schrödinger equations FFT Step with V : MCTDH(F) for Hamiltonian V instead of H larger time steps independently of the space discretization explicit, unconditionally stable scheme L. 2004
59 Error of variational splitting for MCTDH(F) second-order error bound u n u(t n ) C t 2 max u(t) H 2 same as for Strang splitting for the linear Schrödinger equation basic tool: Lie-commutator bounds for vector fields corresponding to T and V L (MCTDH with bounded smooth potential) Koch & L (MCTDHF with Coulomb potential)
60 Outline MCTDH MCTDHF Variational splitting integrator Modelling error of MCTDH
61 Convergence of MCTDH to Schrödinger solution??? Naive expectation: Taking more and more terms in the linear combination of Hartree products yields an ever better accuracy Obstructions: approximation properties of basis of Hartree products Φ J (, t)? ill-conditioned coefficient tensor (a J ) ill-conditioned density matrices ρ (n) = A (n) A (n)
62 Quasi-optimality of variational approximation u(t) MCTDH approximation, ψ(t) exact wave function MCTDH error bounded by best-approximation error: t u(t) ψ(t) d(t) + C κ e cκt d(τ) dτ with d(t) = dist(ψ(t), M) and κ curvature of manifold M 0 but κ cond(ρ (n) ) 1/2 as number of configurations L. 2005
63 Objectives error bounds in the case of ill-conditioned density matrices convergence as number of configurations under appropriate assumptions Conte & L. 2010
64 Approximability assumption Exact wavefunction can be written as ψ(t) = v(t) + e(t) where v(t) M has small defect: i v (, t) Hv(, t) t ε Small error e(t) and small defect (ε) for some linear combination of Hartree products
65 Assumptions on the coefficient tensor B (n) (t) nth unfolding of coefficient tensor (b j1...j N (t)) of v(t) M, for n = 1,..., N Bound of the pseudo-inverse (large: small δ!): B (n) (t) 2 = 1 σ rn (B (n) (t)) δ 1 Small time derivatives of the coefficient tensor in components that correspond to small singular values: B (n) (t) Ḃ (n) (t) c, 2
66 Error bound u(t) ψ(t) e(t) + 2tε for t = O(δ/ε) Bound independent of δ (ill-conditioning) for fixed rank, can achieve δ ε by perturbing v(t): t = O(1) Convergence: if we assume δ ε, we can admit arbitrary number of configurations: ε 0
67 Proof of the error bound use estimates for the tangent space projection P(u) to show a quadratic differential inequality d dt u v C δ u v 2 + ε uses bound of the potential V at present no error bound for MCTDHF with Coulomb potential
68 Outline MCTDH MCTDHF Variational splitting integrator Modelling error of MCTDH
69 Computing semiclassical quantum dynamics using Hagedorn wavepackets Christian Lubich, Univ. Tübingen joint work with Erwan Faou and Vasile Gradinaru Bressanone/Brixen, 17 February 2011
70 Outline The Schrödinger equation in the semi-classical regime Hagedorn wavepackets A splitting method for time integration
71 Outline The Schrödinger equation in the semi-classical regime Hagedorn wavepackets A splitting method for time integration
72 Schrödinger equation in semi-classical scaling iε ψ ε2 (x, t) = t 2m xψ(x, t) + V (x)ψ(x, t) for the wavefunction ψ = ψ(x, t), x = (x 1,..., x N ) R N, t 0 initial value problem: ψ specified at time t = 0 SE for the nuclei in a molecule 0 < ε 1
73 Computational challenges solutions are highly oscillatory with wavelengths ε localized with width ε, with velocity 1 wavepacket near classical motion m q = V (q) semi-classical quantum dynamics high dimension: N = 3 n particles no grids! (neither full nor sparse)
74 Rescue? wavefunction is well approximated by complex Gaussian polynomial Hagedorn wavepackets
75 Outline The Schrödinger equation in the semi-classical regime Hagedorn wavepackets A splitting method for time integration
76 Complex Gaussians in Hagedorn s parametrization ϕ 0 [q, p, Q, P](x) = (πε) N/4 (det Q) 1/2 ( i exp 2ε (x q)t PQ 1 (x q) + i ) ε pt (x q), q R N position, p R N momentum Q, P complex N N matrices such that ( ) Re Q Im Q Y = is symplectic: Y T JY = J for J = Re P Im P ( ) 0 I I 0 Consequence: PQ 1 is complex symmetric with positive definite imaginary part Hagedorn 1980
77 Hagedorn wavepackets L 2 -orthonormal set of functions ϕ k (x) = ϕ k [q, p, Q, P](x) for multi-indices k = (k 1,..., k N ), constructed recursively: define the raising operator R = (R j ) = 1 ) (P (x q) + Q ( iε x p) 2ε With j = ( ) the jth unit vector, set ϕ k+ j = 1 kj + 1 R jϕ k. ϕ k are polynomials of degree k k N multiplied with the Gaussian ϕ 0 (N = 1: Hermite functions). Hagedorn 1998
78 Recursive evaluation Q( kj + 1 ϕ k+ j (x) ) N j=1 = 2 ε (x q)ϕ k(x) Q ( kj ϕ k j (x)) N j=1
79 Approximate wavefunction by Hagedorn wavepacket ψ(x, t) e is(t)/ε k K c k (t) ϕ k [q(t), p(t), Q(t), P(t)](x) over multi-index set K in low dimensions, full cube: k j K (j = 1,..., N) in moderate dimensions, hyperbolic cross: (1 + k 1 )... (1 + k N ) K in high dimensions, axes: k j > 0 only for a single component j in each k problem-adapted moving basis functions
80
81 Outline The Schrödinger equation in the semi-classical regime Hagedorn wavepackets A splitting method for time integration
82 Recap: Schrödinger equation with the Hamiltonian iε ψ t = Hψ H = T + V composed of the kinetic energy operator and a smooth potential T = ε2 2m x V = V (x).
83 Bits and pieces H = T + U q(t) + W q(t) We can solve exactly the free Schrödinger equation, with the wavefunction remaining in the Hagedorn wavepacket form with unaltered coefficients c k. For a quadratic potential, we can solve exactly the potential equation with the wavefunction remaining in the Hagedorn wavepacket form with the same coefficients c k. For the non-quadratic remainder, we compute the variational approximation of the potential equation on the linear space spanned by the functions ϕ k with fixed parameters q, p, Q, P, letting the coefficients c k vary.
84 Free Schrödinger equation iε ψ t = ε2 2m ψ A time-dependent Hagedorn wavepacket solves the free Schrödinger equation with modified positions q(t) = q(0) + t m p(0) Q(t) = Q(0) + t m P(0) and unchanged momenta p, P and unchanged coefficients c k. change only position q and Q and phase S
85 Quadratic potential iε ψ t = Uψ For a quadratic potential U(x), a time-dependent Hagedorn wavepacket solves the equation with modified momenta p(t) = p(0) t U(q(0)) P(t) = P(0) t 2 U(q(0))Q(0) and unchanged positions q and Q and unchanged coefficients c k. change only momentum p and P and phase S
86 Galerkin approximation for the remainder iε ψ t = W ψ, W = W (x) fix Gauss parameters q, p, Q, P in ϕ k (x) = ϕ k [q, p, Q, P](x) Galerkin condition: determine u(x, t) = k K c k(t)ϕ k (x) from ϕ k, iε t u Wu = 0 k K
87 Galerkin approximation for the remainder (ctd.) Galerkin condition determines the coefficient vector c = (c k ) as c(t) = exp ( it ) ε F c(0) with the Hermitian matrix F = (f kl ), f kl = W (x) ϕ k (x) ϕ l (x) dx R N The integrals are non-oscillatory, approximated by sparse Gauss Hermite quadrature. F = O(ε 3/2 ) if the quadratic Taylor polynomial of W at q vanishes. Therefore, exp ( it ε F ) c(0) is computed efficiently using just a few Lanczos iterations with F. change only coefficients c k
88 Time-stepping algorithm start from position q 0, momentum p 0, phase S 0, width matrices Q 0, P 0 satisfying the symplecticity condition, and coefficients c 0 k ψ(x, t 0 ) u 0 (x) = e is0 /ε k K c 0 k ϕ k[q 0, p 0, Q 0, P 0 ](x) determine approximation u 1 (x) of the same form after time step t using a splitting algorithm
89 Splitting algorithm 1. Half-step of kinetic part: updates q 1/2, Q 1/2, S 1/2,. 2. Full step of potential part: split the potential V (x) = U 1/2 (x) + W 1/2 (x) into its quadratic Taylor polynomial U 1/2 (x) at q 1/2 and the remainder solve with quadratic potential U 1/2 : updates p 1, P 1, S 1/2,+ Galerkin approximation for the non-quadratic remainder W 1/2 : update coefficients c 1 k 3. Half-step of kinetic part: updates q 1, Q 1, S 1.
90 Properties preserves the symplecticity relation of the matrices Q and P preserves the L 2 norm of the wavepacket time-reversible method for position q and momentum p: Störmer-Verlet method for the corresponding classical Hamiltonian system limit of taking the full basis set ϕ k with all k N N : Strang splitting of the Schrödinger equation robust in the semi-classical limit ε 0: approximation in the potential part becomes exact for ε 0, while the kinetic part is solved exactly for all ε.
91 Comparison with the Fourier method Error vs. time in a 2D example (Henon-Heiles) Hagedorn with 20 basis functions, Fourier with 2 20 basis functions
92 Error behaviour Maximum error vs. number of basis functions at t = 1 and t = 5.
93 Computing time vs. dimension CPU time with a sparse Hagedorn wavepacket at fixed K = 8 in dimensions 2 to 12.
94 Flying carpet Squared absolute values of the approximate wave function evaluated on the flying carpet of quadrature points.
95 References G.A. Hagedorn, Comm. Math. Phys. 71, (1980) Ann. Phys. 269: (1998) C.L., blue book, Eur. Math. Soc., Zurich, 2008 E. Faou, V. Gradinaru, C.L., SIAM J. Sci. Comp. 31, (2009) V. Gradinaru, G.A. Hagedorn, A. Joye, JCP 132, (2010)
Semi-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 informationDynamical low-rank approximation
Dynamical low-rank approximation Christian Lubich Univ. Tübingen Genève, Swiss Numerical Analysis Day, 17 April 2015 Coauthors Othmar Koch 2007, 2010 Achim Nonnenmacher 2008 Dajana Conte 2010 Thorsten
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 informationNon Adiabatic Transitions in a Simple Born Oppenheimer Scattering System
1 Non Adiabatic Transitions in a Simple Born Oppenheimer Scattering System George A. Hagedorn Department of Mathematics and Center for Statistical Mechanics, Mathematical Physics, and Theoretical Chemistry
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 informationNon Adiabatic Transitions near Avoided Crossings: Theory and Numerics
Non Adiabatic Transitions near Avoided Crossings: Theory and Numerics Raoul Bourquin a, Vasile Gradinaru a, George A. Hagedorn b April 8, 2011 Abstract We present a review of rigorous mathematical results
More informationTheoretical Photochemistry SoSe 2014
Theoretical Photochemistry SoSe 2014 Lecture 9 Irene Burghardt (burghardt@chemie.uni-frankfurt.de) http://www.theochem.uni-frankfurt.de/teaching/ Theoretical Photochemistry 1 Topics 1. Photophysical Processes
More informationZ. Zhou On the classical limit of a time-dependent self-consistent field system: analysis. computation
On the classical limit of a time-dependent self-consistent field system: analysis and computation Zhennan Zhou 1 joint work with Prof. Shi Jin and Prof. Christof Sparber. 1 Department of Mathematics Duke
More informationNumerical tensor methods and their applications
Numerical tensor methods and their applications 14 May 2013 All lectures 4 lectures, 2 May, 08:00-10:00: Introduction: ideas, matrix results, history. 7 May, 08:00-10:00: Novel tensor formats (TT, HT,
More informationMolecular propagation through crossings and avoided crossings of electron energy levels
Molecular propagation through crossings and avoided crossings of electron energy levels George A. Hagedorn Department of Mathematics and Center for Statistical Mechanics and Mathematical Physics Virginia
More informationProject: Vibration of Diatomic Molecules
Project: Vibration of Diatomic Molecules Objective: Find the vibration energies and the corresponding vibrational wavefunctions of diatomic molecules H 2 and I 2 using the Morse potential. equired Numerical
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
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 informationIntroduction to Electronic Structure Theory
Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these
More informationHierarchical Matrices in Computations of Electron Dynamics
Hierarchical Matrices in Computations of Electron Dynamics Othmar Koch Christopher Ede Gerald Jordan Armin Scrinzi July 20, 2009 Abstract We discuss the approximation of the meanfield terms appearing in
More informationCh 125a Problem Set 1
Ch 5a Problem Set Due Monday, Oct 5, 05, am Problem : Bra-ket notation (Dirac notation) Bra-ket notation is a standard and convenient way to describe quantum state vectors For example, φ is an abstract
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationAttempts at relativistic QM
Attempts at relativistic QM based on S-1 A proper description of particle physics should incorporate both quantum mechanics and special relativity. However historically combining quantum mechanics and
More informationDerivation of Mean-Field Dynamics for Fermions
IST Austria Mathematical Many Body Theory and its Applications Bilbao, June 17, 2016 Derivation = show that solution of microscopic eq. is close to solution of effective eq. (with fewer degrees of freedom)
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 informationExponentially Accurate Semiclassical Tunneling Wave Functions in One Dimension
Exponentially Accurate Semiclassical Tunneling Wave Functions in One Dimension Vasile Gradinaru Seminar for Applied Mathematics ETH Zürich CH 8092 Zürich, Switzerland, George A. Hagedorn Department of
More informationSuperintegrability in a non-conformally-at space
(Joint work with Ernie Kalnins and Willard Miller) School of Mathematics and Statistics University of New South Wales ANU, September 2011 Outline Background What is a superintegrable system Extending the
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 informationPHY 407 QUANTUM MECHANICS Fall 05 Problem set 1 Due Sep
Problem set 1 Due Sep 15 2005 1. Let V be the set of all complex valued functions of a real variable θ, that are periodic with period 2π. That is u(θ + 2π) = u(θ), for all u V. (1) (i) Show that this V
More information1 Time-Dependent Two-State Systems: Rabi Oscillations
Advanced kinetics Solution 7 April, 16 1 Time-Dependent Two-State Systems: Rabi Oscillations a In order to show how Ĥintt affects a bound state system in first-order time-dependent perturbation theory
More informationSection 9 Variational Method. Page 492
Section 9 Variational Method Page 492 Page 493 Lecture 27: The Variational Method Date Given: 2008/12/03 Date Revised: 2008/12/03 Derivation Section 9.1 Variational Method: Derivation Page 494 Motivation
More informationSETTING AND ANALYSIS OF THE MULTICONFIGURATION TIME-DEPENDENT HARTREE-FOCK EQUATIONS
SETTING AND ANALYSIS OF THE MULTICONFIGURATION TIME-DEPENDENT HARTREE-FOCK EQUATIONS CLAUDE BARDOS, ISABELLE CATTO, NORBERT J. MAUSER, AND SABER TRABELSI Abstract. In this paper we motivate, formulate
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More information5.4 Given the basis e 1, e 2 write the matrices that represent the unitary transformations corresponding to the following changes of basis:
5 Representations 5.3 Given a three-dimensional Hilbert space, consider the two observables ξ and η that, with respect to the basis 1, 2, 3, arerepresentedby the matrices: ξ ξ 1 0 0 0 ξ 1 0 0 0 ξ 3, ξ
More informationQuestioning Quantum Mechanics? Kurt Barry SASS Talk January 25 th, 2012
Questioning Quantum Mechanics? Kurt Barry SASS Talk January 25 th, 2012 2 Model of the Universe Fundamental Theory Low-Energy Limit Effective Field Theory Quantum Mechanics Quantum Mechanics is presently
More informationWave function methods for the electronic Schrödinger equation
Wave function methods for the electronic Schrödinger equation Zürich 2008 DFG Reseach Center Matheon: Mathematics in Key Technologies A7: Numerical Discretization Methods in Quantum Chemistry DFG Priority
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 informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationLecture I: Constrained Hamiltonian systems
Lecture I: Constrained Hamiltonian systems (Courses in canonical gravity) Yaser Tavakoli December 15, 2014 1 Introduction In canonical formulation of general relativity, geometry of space-time is given
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 informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationMassachusetts Institute of Technology Physics Department
Massachusetts Institute of Technology Physics Department Physics 8.32 Fall 2006 Quantum Theory I October 9, 2006 Assignment 6 Due October 20, 2006 Announcements There will be a makeup lecture on Friday,
More information4 Quantum Mechanical Description of NMR
4 Quantum Mechanical Description of NMR Up to this point, we have used a semi-classical description of NMR (semi-classical, because we obtained M 0 using the fact that the magnetic dipole moment is quantized).
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 informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
More informationIntroduction to Integrability
Introduction to Integrability Problem Sets ETH Zurich, HS16 Prof. N. Beisert, A. Garus c 2016 Niklas Beisert, ETH Zurich This document as well as its parts is protected by copyright. This work is licensed
More informationLecture 3 Dynamics 29
Lecture 3 Dynamics 29 30 LECTURE 3. DYNAMICS 3.1 Introduction Having described the states and the observables of a quantum system, we shall now introduce the rules that determine their time evolution.
More informationDynamical low rank approximation in hierarchical tensor format
Dynamical low rank approximation in hierarchical tensor format R. Schneider (TUB Matheon) John von Neumann Lecture TU Munich, 2012 Motivation Equations describing complex systems with multi-variate solution
More informationMean-Field Limits for Large Particle Systems Lecture 2: From Schrödinger to Hartree
for Large Particle Systems Lecture 2: From Schrödinger to Hartree CMLS, École polytechnique & CNRS, Université Paris-Saclay FRUMAM, Marseilles, March 13-15th 2017 A CRASH COURSE ON QUANTUM N-PARTICLE DYNAMICS
More informationL 2 Geometry of the Symplectomorphism Group
University of Notre Dame Workshop on Innite Dimensional Geometry, Vienna 2015 Outline 1 The Exponential Map on D s ω(m) 2 Existence of Multiplicity of Outline 1 The Exponential Map on D s ω(m) 2 Existence
More informationAn Exactly Solvable 3 Body Problem
An Exactly Solvable 3 Body Problem The most famous n-body problem is one where particles interact by an inverse square-law force. However, there is a class of exactly solvable n-body problems in which
More informationGEOMETRIC QUANTIZATION
GEOMETRIC QUANTIZATION 1. The basic idea The setting of the Hamiltonian version of classical (Newtonian) mechanics is the phase space (position and momentum), which is a symplectic manifold. The typical
More informationPHY 396 K. Problem set #5. Due October 9, 2008.
PHY 396 K. Problem set #5. Due October 9, 2008.. First, an exercise in bosonic commutation relations [â α, â β = 0, [â α, â β = 0, [â α, â β = δ αβ. ( (a Calculate the commutators [â αâ β, â γ, [â αâ β,
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
More informationThe Time{Dependent Born{Oppenheimer Approximation and Non{Adiabatic Transitions
The Time{Dependent Born{Oppenheimer Approximation and Non{Adiabatic Transitions George A. Hagedorn Department of Mathematics, and Center for Statistical Mechanics, Mathematical Physics, and Theoretical
More informationKet space as a vector space over the complex numbers
Ket space as a vector space over the complex numbers kets ϕ> and complex numbers α with two operations Addition of two kets ϕ 1 >+ ϕ 2 > is also a ket ϕ 3 > Multiplication with complex numbers α ϕ 1 >
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 informationMatrix stabilization using differential equations.
Matrix stabilization using differential equations. Nicola Guglielmi Universitá dell Aquila and Gran Sasso Science Institute, Italia NUMOC-2017 Roma, 19 23 June, 2017 Inspired by a joint work with Christian
More informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
More informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
More information3: Many electrons. Orbital symmetries. l =2 1. m l
3: Many electrons Orbital symmetries Atomic orbitals are labelled according to the principal quantum number, n, and the orbital angular momentum quantum number, l. Electrons in a diatomic molecule experience
More informationDynamical Low-Rank Approximation to the Solution of Wave Equations
Dynamical Low-Rank Approximation to the Solution of Wave Equations Julia Schweitzer joint work with Marlis Hochbruck INSTITUT FÜR ANGEWANDTE UND NUMERISCHE MATHEMATIK 1 KIT Universität des Landes Baden-Württemberg
More informationStructure relations for the symmetry algebras of classical and quantum superintegrable systems
UNAM talk p. 1/4 Structure relations for the symmetry algebras of classical and quantum superintegrable systems Willard Miller miller@ima.umn.edu University of Minnesota UNAM talk p. 2/4 Abstract 1 A quantum
More informationwhere P a is a projector to the eigenspace of A corresponding to a. 4. Time evolution of states is governed by the Schrödinger equation
1 Content of the course Quantum Field Theory by M. Srednicki, Part 1. Combining QM and relativity We are going to keep all axioms of QM: 1. states are vectors (or rather rays) in Hilbert space.. observables
More information1 The postulates of quantum mechanics
1 The postulates of quantum mechanics The postulates of quantum mechanics were derived after a long process of trial and error. These postulates provide a connection between the physical world and the
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationChaos, Quantum Mechanics and Number Theory
Chaos, Quantum Mechanics and Number Theory Peter Sarnak Mahler Lectures 2011 Hamiltonian Mechanics (x, ξ) generalized coordinates: x space coordinate, ξ phase coordinate. H(x, ξ), Hamiltonian Hamilton
More information2 Electronic structure theory
Electronic structure theory. Generalities.. Born-Oppenheimer approximation revisited In Sec..3 (lecture 3) the Born-Oppenheimer approximation was introduced (see also, for instance, [Tannor.]). We are
More informationLecture #8: Quantum Mechanical Harmonic Oscillator
5.61 Fall, 013 Lecture #8 Page 1 Last time Lecture #8: Quantum Mechanical Harmonic Oscillator Classical Mechanical Harmonic Oscillator * V(x) = 1 kx (leading term in power series expansion of most V(x)
More informationIntroduction to Modern Quantum Field Theory
Department of Mathematics University of Texas at Arlington Arlington, TX USA Febuary, 2016 Recall Einstein s famous equation, E 2 = (Mc 2 ) 2 + (c p) 2, where c is the speed of light, M is the classical
More informationHamiltonian Systems of Negative Curvature are Hyperbolic
Hamiltonian Systems of Negative Curvature are Hyperbolic A. A. Agrachev N. N. Chtcherbakova Abstract The curvature and the reduced curvature are basic differential invariants of the pair: Hamiltonian system,
More informationComputational Methods of Physical Problems with Mixed Scales
Computational Methods of Physical Problems with Mixed Scales Shi Jin University of Wisconsin-Madison IPAM Tutorial, March 10-12, 2008 Scales and Physical Laws Figure from E & Engquist, AMS Notice Connections
More informationSecond Quantization Method for Bosons
Second Quantization Method for Bosons Hartree-Fock-based methods cannot describe the effects of the classical image potential (cf. fig. 1) because HF is a mean-field theory. DFF-LDA is not able either
More informationCOULOMB SYSTEMS WITH CALOGERO INTERACTION
PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Mathematical Sciences 016, 3, p. 15 19 COULOMB SYSTEMS WITH CALOGERO INTERACTION P h y s i c s T. S. HAKOBYAN, A. P. NERSESSIAN Academician G. Sahakyan
More informationComputing High Frequency Waves By the Level Set Method
Computing High Frequency Waves By the Level Set Method Hailiang Liu Department of Mathematics Iowa State University Collaborators: Li-Tien Cheng (UCSD), Stanley Osher (UCLA) Shi Jin (UW-Madison), Richard
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Using Quadrature-Based Projection Methods Julian Köllermeier, Manuel Torrilhon 29th International Symposium on Rarefied Gas Dynamics, Xi an July 14th, 2014 Julian Köllermeier, Manuel Torrilhon 1 / 12 Outline
More informationHigh-order time-splitting methods for Schrödinger equations
High-order time-splitting methods for Schrödinger equations and M. Thalhammer Department of Mathematics, University of Innsbruck Conference on Scientific Computing 2009 Geneva, Switzerland Nonlinear Schrödinger
More 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 informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
More informationUniformly accurate averaging numerical schemes for oscillatory evolution equations
Uniformly accurate averaging numerical schemes for oscillatory evolution equations Philippe Chartier University of Rennes, INRIA Joint work with M. Lemou (University of Rennes-CNRS), F. Méhats (University
More informationSuperintegrable 3D systems in a magnetic field and Cartesian separation of variables
Superintegrable 3D systems in a magnetic field and Cartesian separation of variables in collaboration with L. Šnobl Czech Technical University in Prague GSD 2017, June 5-10, S. Marinella (Roma), Italy
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationG : Quantum Mechanics II
G5.666: Quantum Mechanics II Notes for Lecture 5 I. REPRESENTING STATES IN THE FULL HILBERT SPACE Given a representation of the states that span the spin Hilbert space, we now need to consider the problem
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationQuantum Field Theory III
Quantum Field Theory III Prof. Erick Weinberg January 19, 2011 1 Lecture 1 1.1 Structure We will start with a bit of group theory, and we will talk about spontaneous symmetry broken. Then we will talk
More informationBRST and Dirac Cohomology
BRST and Dirac Cohomology Peter Woit Columbia University Dartmouth Math Dept., October 23, 2008 Peter Woit (Columbia University) BRST and Dirac Cohomology October 2008 1 / 23 Outline 1 Introduction 2 Representation
More informationTHE SIEGEL UPPER HALF SPACE IS A MARSDEN WEINSTEIN QUOTIENT: SYMPLECTIC REDUCTION AND GAUSSIAN WAVE PACKETS
THE SIEGEL UPPER HALF SPACE IS A MARSDEN WEINSTEIN QUOTIENT: SYMPLECTIC REDUCTION AND GAUSSIAN WAVE PACKETS TOMOKI OHSAWA Abstract. We show that the Siegel upper half space Σ d is identified with the Marsden
More informationTowards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity
Towards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity Melvin Leok Mathematics, University of Michigan, Ann Arbor. Joint work with Mathieu Desbrun, Anil Hirani, and Jerrold
More informationThe Density Matrix for the Ground State of 1-d Impenetrable Bosons in a Harmonic Trap
The Density Matrix for the Ground State of 1-d Impenetrable Bosons in a Harmonic Trap Institute of Fundamental Sciences Massey University New Zealand 29 August 2017 A. A. Kapaev Memorial Workshop Michigan
More informationMean-Field Theory: HF and BCS
Mean-Field Theory: HF and BCS Erik Koch Institute for Advanced Simulation, Jülich N (x) = p N! ' (x ) ' 2 (x ) ' N (x ) ' (x 2 ) ' 2 (x 2 ) ' N (x 2 )........ ' (x N ) ' 2 (x N ) ' N (x N ) Slater determinant
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
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 informationQUALIFYING EXAMINATION, Part 2. Solutions. Problem 1: Quantum Mechanics I
QUALIFYING EXAMINATION, Part Solutions Problem 1: Quantum Mechanics I (a) We may decompose the Hamiltonian into two parts: H = H 1 + H, ( ) where H j = 1 m p j + 1 mω x j = ω a j a j + 1/ with eigenenergies
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 informationIntroduction to Quantum Mechanics PVK - Solutions. Nicolas Lanzetti
Introduction to Quantum Mechanics PVK - Solutions Nicolas Lanzetti lnicolas@student.ethz.ch 1 Contents 1 The Wave Function and the Schrödinger Equation 3 1.1 Quick Checks......................................
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 informationAN ACCELERATED SURFACE-HOPPING METHOD FOR COMPUTATIONAL SEMICLASSICAL MOLECULAR DYNAMICS. Laren K. Mortensen
AN ACCELERATED SURFACE-HOPPING METHOD FOR COMPUTATIONAL SEMICLASSICAL MOLECULAR DYNAMICS by Laren K. Mortensen A senior thesis submitted to the faculty of Brigham Young University in partial fulfillment
More informationPhase-space approach to bosonic mean field dynamics
field dynamics sept. 5th, 2013 Outline 1 2 3 4 N-body bosonic problem Let Z = L 2 (X, m; C) be the one particle space. e.g. (X, m) measured add. group. Z separable. X R d, X Z d, X = Z d /(KZ) d. The bosonic
More informationSymmetric multistep methods for charged-particle dynamics
Version of 12 September 2017 Symmetric multistep methods for charged-particle dynamics Ernst Hairer 1, Christian Lubich 2 Abstract A class of explicit symmetric multistep methods is proposed for integrating
More informationIntroduction to Hamiltonian systems
Introduction to Hamiltonian systems Marlis Hochbruck Heinrich-Heine Universität Düsseldorf Oberwolfach Seminar, November 2008 Examples Mathematical biology: Lotka-Volterra model First numerical methods
More informationarxiv:hep-th/ v1 8 Mar 1995
GALILEAN INVARIANCE IN 2+1 DIMENSIONS arxiv:hep-th/9503046v1 8 Mar 1995 Yves Brihaye Dept. of Math. Phys. University of Mons Av. Maistriau, 7000 Mons, Belgium Cezary Gonera Dept. of Physics U.J.A Antwerpen,
More information