Convergence analysis of high-order commutator-free quasi-magnus exponential integrators for non-autonomous linear Schrödinger equations
|
|
- Richard Sherman
- 5 years ago
- Views:
Transcription
1 IMA Journal of Numerical Analysis 8) Page of doi:.9/imanum/drn Convergence analysis of high-order commutator-free quasi-magnus exponential integrators for non-autonomous linear Schrödinger equations SERGIO BLANES Universitat Politècnica de València, Instituto de Matemática Multidisciplinar, 46 Valencia, Spain. FERNANDO CASAS Universitat Jaume I, IMAC and Departament de Matemàtiques, 7 Castellón, Spain. fernando.casas@uji.es CESÁREO GONZÁLEZ Universidad de Valladolid, Departamento de Matemática Aplicada, 47 Valladolid, Spain. ome@am.uva.es MECHTHILD THALHAMMER Leopold Franzens Universität Innsbruck, Institut für Mathematik, 6 Innsbruck, Austria. mechthild.thalhammer@uibk.ac.at [Received on 7 October 8] This work is devoted to the derivation of a convergence result for high-order commutator-free quasi- Magnus CFQM) exponential integrators applied to non-autonomous linear Schrödinger equations; a detailed stability and local error analysis is provided for the relevant special case, where the Hamilton operator comprises the Laplacian and a regular space-time-dependent potential. In the context of nonautonomous linear ordinary differential equations, CFQM exponential integrators are composed of exponentials involving linear combinations of certain values of the associated time-dependent matrix; this approach extends to non-autonomous linear evolution equations given by unbounded operators. An inherent advantage of CFQM exponential integrators over other time integration methods such as Runge Kutta methods or Magnus integrators is that structural properties of the underlying operator family affecting stability are well-preserved; this characteristic is confirmed by a theoretical analysis ensuring unconditional stability in the underlying Hilbert space and the full order of convergence under low regularity requirements on the initial state. Due to the fact that convenient tools for products of matrix exponentials such as the Baker Campbell Hausdorff formula involve infinite series and thus cannot be applied in connection with unbounded operators, a certain complexity in the investigation of higher-order CFQM exponential integrators for Schrödinger equations is related to an appropriate treatment of compositions of evolution operators; an effective concept for the derivation of a local error expansion relies on suitable linearisations of the evolution equations for the exact and numerical solutions, representations by the variation-of-constants formula, and Taylor series expansions of parts of the integrands, where the arising iterated commutators determine the regularity requirements on the problem data. Keywords: Non-autonomous linear evolution equations, Schrödinger equations, Quantum systems, Time integration methods, Exponential integrators, Magnus integrators, Commutator-free quasi-magnus exponential integrators, Stability, Local error, Convergence AMS: 65L5, 65M, 65J c The author 8. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.
2 of S. BLANES, F. CASAS, C. GONZÁLEZ, M. THALHAMMER. Introduction Commutator-free quasi-magnus exponential integrators. Numerical experiments for nonautonomous linear differential equations confirm the favourable stability and error behaviour of exponential methods based on the Magnus expansion in comparison with classical time integration methods, see for instance BLANES, CASAS, OTEO, ROS 9); BLANES, MOAN 6); HOCHBRUCK, LU- BICH ); ISERLES, KROPIELNICKA, SINGH 8) and references given therein; however, in the context of large-scale applications, it is advantageous to avoid the actual computation and storage of commutators. As illustrated in ALVERMANN, FEHSKE ); ALVERMANN, FEHSKE, LITTLEWOOD ); BADER, BLANES, KOPYLOV 8); BLANES, CASAS, THALHAMMER 7), the replacement of a single matrix exponential involving iterated commutators by compositions of matrix exponentials leads to numerical approximations that are potentially superior to those obtained by standard Magnus integrators; we refer to this type of time integration methods as commutator-free quasi-magnus CFQM) exponential integrators. In our former works BLANES, CASAS, THALHAMMER 7, 8), the focus is on the study of CFQM exponential integrators applied to dissipative quantum systems or partial differential equations of parabolic type, respectively; during the development of these contributions, we also performed numerical tests for Hamiltonian systems and Schrödinger equations. Our comparisons confirmed the excellent behaviour of CFQM exponential integrators and motivated further theoretical investigations. Stability and error analysis. In the present work, our objective is the derivation of a convergence result for high-order CFQM exponential integrators in the context of non-autonomous linear Schrödinger equations; we thereby proceed and complete the investigations in BLANES, CASAS, THALHAMMER 7, 8) on the design, numerical comparison, and convergence analysis of CFQM exponential integrators for ordinary and partial differential equations of parabolic type. For the sake of concretion, we provide a detailed stability and local error analysis for the relevant special case, where the Hamilton operator comprises the Laplacian and a regular space-time-dependent potential. For non-autonomous linear ordinary differential equations, CFQM exponential integrators are composed of exponentials involving linear combinations of certain values of the associated time-dependent matrix; it is straightforward to extend this approach to non-autonomous linear evolution equations that are defined by unbounded operators. Concerning stability, an inherent advantage of CFQM exponential integrators over other classes of time integration methods such as Runge Kutta methods or standard Magnus integrators is that structural properties of the underlying operator family are well-preserved; this characteristic is confirmed by our theoretical analysis, which implies unconditional stability in the underlying Hilbert space and the full order of convergence under low regularity requirements on the initial state. Convenient tools for analysing products of matrix exponentials such as the Baker Campbell Hausdorff formula involve infinite series and thus cannot be applied in connection with unbounded operators; consequently, a certain complexity in the local error analysis of higher-order CFQM exponential integrators for Schrödinger equations is related to an appropriate treatment of compositions of evolution operators. We explore an effective approach that relies on suitable linearisations of the evolution equations satisfied by the exact and numerical solutions, representations by the variation-ofconstants formula, and Taylor series expansions of parts of the integrands; under sufficient regularity requirements on the initial state and the space-time-dependent potential, boundedness of the arising iterated commutators and well-definedness of the resulting local error expansion is ensured. Outline. The present manuscript has the following structure. In Section, we formulate non-
3 CONVERGENCE ANALYSIS OF CFQM EXPONENTIAL INTEGRATORS of autonomous linear Schrödinger equations as abstract evolution equations on Hilbert spaces and specify fundamental assumptions on the defining Hamilton operators. In Section, we introduce CFQM exponential integrators and recapitulate second-order and fourth-order example methods. In Section 4, we provide a detailed stability and local error analysis for the special case, where the Schrödinger equation is given by the Laplacian and a regular space-time-dependent potential. In Section 5, we indicate the extension to the general case.. Class of evolution equations Analytical framework. Our analytical framework for a convergence analysis of CFQM exponential integrators applied to Schrödinger equations relies on the theory of self-adjoint operators on Hilbert spaces; for a comprehensive treatment of basic mathematical means and the principles of quantum theory, we refer to BONGAARTS 5); ENGEL, NAGEL ). Evolution equation. We study the initial value problem for a non-autonomous linear evolution equation of Schrödinger type, cast into the abstract form { u t) = At)ut) = i Ht)ut), t t,t ),.a) ut ) given; throughout, we employ the fundamental assumption that the arising time-dependent Hamilton operator is self-adjoint on the underlying Hilbert space Ht) : D H H, t [t,t ],.b) and that the associated evolution operator is unitary ut) = E t t,t )ut ), E t t,t ) H H =, t [t,t ]..c) Unitarity. For a regular solution to.a), differentiation with respect to time confirms.c); indeed, with the help of the chain rule and due to the required self-adjointness of the Hamilton operator, the following relations implying the unitarity of the evolution operator are obtained d dt ut) H = R u t) ) ut) H = R i Ht)ut) ut) )H =, t t,t ), ut) H = ut ).) H, t [t,t ]. Applications. In view of concrete applications, we primarily study the relevant special case, where the arising Hamilton operator comprises the Laplacian and a space-time-dependent potential, see Section 4; in this situation, it is natural to consider square-integrable complex-valued functions that are defined on the Euclidean space or on a cartesian product of bounded intervals, respectively. That is, the underlying Hilbert space is given by a Lebesgue space H = L Ω,C), Ω R d, and the common domain of the defining operator family is related to a subspace of a Sobolev space D H Ω,C) = W Ω,C). We point out that in the context of Schrödinger equations generally no additional consistency conditions have to be taken into account.
4 4 of S. BLANES, F. CASAS, C. GONZÁLEZ, M. THALHAMMER. Class of time integration methods CFQM exponential integrators. Our main objective is to provide a rigorous convergence analysis for commutator-free quasi-magnus CFQM) exponential integrators applied to evolution equations of Schrödinger type.). As standard in a time-stepping approach, we choose time grid points with corresponding time stepsizes and determine approximations to the exact solution values by recurrence t < < t n < < t N = T, τ n = t n+ t n, u n+ = S τ n,t n )u n ut n+ ) = E τ n,t n )ut n ), n {,,...,N };.a) for a CFQM exponential integrator, the numerical evolution operator is given by the composition of several exponentials that involve linear combinations of certain values of the defining operator S τ n,t n ) = J j= B j τ n,t n ) = K k= a jk At n + c k τ n ), j {,...,J}, e τ n B j τ n,t n ) = e τ n B J τ n,t n) e τ n B τ n,t n ), n {,,...,N }..b) A proper choice of the underlying quadrature nodes and of the associated real method coefficients a jk R, b j = c k [,], k {,...,K},.c) K k= a jk, j {,...,J}, k {,...,K},.d) permits to attain the desired order of consistency p N ; that is, in the context of nonstiff nonautonomous linear ordinary differential equations with associated families of regular time-dependent matrices, an expansion of the local error in the first step with respect to the time stepsize leads to a relation of the form L τ,t ) = S τ,t ) E τ,t ) = O Examples. The simplest instance of a second-order CFQM exponential integrator.) is the exponential midpoint rule p =, J = K =, c =, a = b =,.) S τ n,t n ) = e τ n At n + τn), n {,,...,N }. The two-stage Gaussian quadrature rule leads to a canonical fourth-order CFQM exponential integrator τ p+ ). p = 4, J = K =, α = 6, c = α, c = + α, a = a = 4 + α, a = a = 4 α, b = a + a = b = a + a =, B j τ n,t n ) = a j At n + c τ n ) + a j At n + c τ n ), j {,}, S τ n,t n ) = e τ n B τ n,t n ) e τ n B τ n,t n ), n {,,...,N }..) For additional examples and detailed information on order conditions, we refer to BLANES, CASAS, THALHAMMER 7, 8).
5 CONVERGENCE ANALYSIS OF CFQM EXPONENTIAL INTEGRATORS 5 of Stability. Within the analytical framework introduced in Section, Stone s Theorem implies that any one-parameter group generated by the defining time-dependent linear operator is unitary on the underlying Hilbert space, i.e. e σ At) H H =, σ R, t [t,t ],.4) see ENGEL, NAGEL ). By means of this fundamental result, we prove that time integration methods of the format.) and in particular the schemes.) and.) are unconditionally stable, in consistence with the characteristics of the exact evolution operator.), see Section 4.. Remark. It is natural to compare CFQM exponential integrators of the form.) to Magnus integrators, which are given by a single exponential involving iterated comutators of certain values of the defining operator, see BLANES, CASAS, OTEO, ROS 9); BLANES, MOAN 6); HOCHBRUCK, LUBICH ); ISERLES, KROPIELNICKA, SINGH 8) and references given therein. For both classes of methods, it is ensured that the numerical evolution operator is unitary and hence preserves fundamental properties of the exact evolution operator associated with a Schrödinger equation. As indicated before, for CFQM exponential integrators, this result relies on an application of Stone s Theorem; for a Magnus integrator, it is used that the commutator is hermitian, which in essence follows from the required self-adjointness of the Hamilton operator [As),At) ] w As)At) v )H = ) ) w At)As) v = Hs)Ht)v w ) H H + Ht)Hs)v w ) H = Ht)v ) Hs)w H + Hs)v ) Ht)w H = v ) Ht)Hs)w H + v Hs)Ht)w )H = v At)As) As)At) ) w )H = v [ As),At) ] ) w However, a basic difference between CFQM exponential integrators and Magnus integrators concerns well-definedness and stability. An inherent advantage of CFQM exponential integrators is that the linear combinations arising in.) preserve the structure of the underlying operator family, and this is reflected in a favourable stability behaviour; contrary, the stability of higher-order Magnus integrators is only ensured in a spatially semi-discrete setting and under major stepsize restrictions, see HOCHBRUCK, LUBICH ). Schrödinger versus parabolic equations. In order to ensure the stability of CFQM exponential integrators for non-autonomous linear evolution equations of Schrödinger type, we have to restrict ourselves to schemes involving real coefficients. In the context of dissipative quantum systems and parabolic evolution equations, however, it is of interest to study and design CFQM exponential integrators with complex method coefficients, see BLANES, CASAS, THALHAMMER 7, 8). 4. Special case Situation. In this section, we focus on non-autonomous linear Schrödinger equations, where the Hamilton operator is given by the Laplacian with respect to the spatial variables and a regular realvalued space-time-dependent potential acting as multiplication operator Ht) = +V t), V t) = V,t) : Ω R, t [t,t ]; 4.) in view of the abstract formulation.), we again suppress the space dependence. The consideration of this special case is of relevance in applications, see BADER, BLANES, KOPYLOV 8); HOCHBRUCK, LUBICH ); ISERLES, KROPIELNICKA, SINGH 8) and references given therein. H.
6 6 of S. BLANES, F. CASAS, C. GONZÁLEZ, M. THALHAMMER Notation. With regard to the formulation and derivation of our main result, we henceforth employ convenient abbreviations for Sobolev spaces and related norms H m = H m Ω,C) = W m Ω,C), m N, GW l,m [t,t ],Ω,K) = max G j) t) L [t,t ],W m Ω,K)), l,m N, j {,,...,l} where K = R or K = C, respectively; evidently, H coincides with the underlying Lebesgue space H = L Ω,C). 4. Convergence result Theoretical result. In the following, we state a theoretical result, which confirms the favourable convergence behaviour of CFQM exponential integrators for evolution equations of Schrödinger type; compared to other classes of time integration methods, stability is ensured without any restriction on the time stepsizes, and the full order of convergence is obtained under low regularity requirements on the initial state. Theorem Assume that the considered non-autonomous linear Schrödinger equation.) is defined by a Hamilton operator of the special structure 4.). Provided that the arising real-valued space-timedependent potential and the prescribed initial state fulfill suitable regularity requirements such that C V = max V W l,p l) [t,t ],Ω,R) <, ut ) H p <, l {,...,p} any pth-order CFQM exponential integrator.) satisfies a global error estimate of the form un ut ) H u ut ) H +C τmax p ut ) H p, τ max = max { τ n : n {,,...,N } } ; the arising constant depends on an upper bound for the time stepsizes, the length of the considered time interval [t,t ], the bound C V, and the method coefficients, but it is independent of the number of time steps n and the actual time stepsizes. Proof. Our proof relies on stability bounds and local error estimates, which we deduce in the subsequent Sections 4. and 4.. The constant arising in 4.9) in particular depends on bounds for certain space and time derivatives of the potential on the interval [t n,t n+ ], where n {,,...,N }; estimation by C V leads to a constant that is independent of n {,,...,N }. Similarly, we estimate the constant in 4.) by a constant that depends on the length of the time interval [t,t ] and the uniform bound C V. In order to relate the global error to local errors, we employ the telescopic identity N u N ut ) = n= S τ n,t n ) u ut ) ) N N + n= ν=n+ S τ ν,t ν ) S τ n,t n ) E τ n,t n ) ) E t n t,t )ut ); due to the fact that the numerical evolution operator associated with a CFQM exponential integrator involving real coefficients is unitary S τn,t n ) H H =, n {,,...,N },
7 CONVERGENCE ANALYSIS OF CFQM EXPONENTIAL INTEGRATORS 7 of see 4.), straightforward estimation in the underlying Hilbert space yields un ut ) H u ut ) H + N n= S τ n,t n ) E τ n,t n ) ) E t n t,t )ut ) H. In our expansion of the local error operator with respect to the time stepsize, iterated commutators of the Laplacian and the potential up to order p arise naturally and lead to an estimate of the form S τ n,t n ) E τ n,t n ) H H p C τn p+, n {,,...,N }, see 4.9); for this reason, we have to embed the decisive contributions S τn,t n ) E τ n,t n ) ) E t n t,t )ut ), n {,,...,N }, in the Sobolev space H p. Using the stability bound for the exact evolution operator E tn t,t ) H p H p C, n {,,...,N }, see 4.) and 4.6), we obtain the relation un ut ) H u ut ) H +C N n= τn p+ ut ) H p, which implies the stated global error estimate. 4. Stability Stability in Sobolev spaces Exact evolution operator). The exact solution to the non-autonomous linear Schrödinger equation.) with Hamilton operator given by 4.) inherits the spatial regularity of the potential and the initial state. As already discussed in Section, the evolution operator associated with a self-adjoint Hamilton operator is unitary on the underlying Hilbert space E t t,t ) H H =, t [t,t ], 4.a) see.c) and.). Moreover, an estimate of the form E t t,t ) H m H m +C m τ, ) V W,m t t [t,t],ω,r) ) C, t t τ, t [t,t ], m N, 4.b) holds with a constant that depends on an upper bound for the length of the considered time interval and uniform bounds for certain space derivatives of the potential; the form of the bound and in particular the factor t t are essential in view of the derivation of an analogous stability bound for compositions of the numerical evolution operator. In order to exemplify this result, we consider a non-autonomous linear Schrödinger equation in a single space dimension on the unbounded domain t ψx,t) = i Hx,t)ψx,t) = i xx ψx,t) +V x,t)ψx,t) ), x,t) R t,t );
8 8 of S. BLANES, F. CASAS, C. GONZÁLEZ, M. THALHAMMER the first and second spatial derivatives of the solution satisfy the related evolution equations χ = x ψ, χ = xx ψ, t χ m x,t) = i Hx,t) χ m x,t) + R m x,t) ), m {,}, { x V x,t)ψx,t) if m =, R m x,t) = x,t) R t,t ). x V x,t) χ x,t) + xx V x,t)ψx,t) if m =, A representation by means of the variation-of-constants formula together with the application of the unitarity of the evolution operator.c) leads to the estimate χ,t) L Ω,C) χ,t ) t L Ω,C) + R,s) L Ω,C) ds t χ,t ) L Ω,C) + t t ) x V L ψ,t Ω [t,t],r) ) L Ω,C), t [t,t ]; as a consequence, we obtain the relation χ,t) L Ω,C) χ,t ) t L Ω,C) + R,s) L Ω,C) ds t χ,t ) L Ω,C) + t t ) x V L sup χ Ω [t,t],r),s) L Ω,C) s [t,t] + t t ) xx V ψ,t L ) Ω [t,t],r) L Ω,C) χ,t ) L Ω,C) + t t ) x V χ L,t Ω [t,t],r) ) L Ω,C) + t t ) t t ) x V L Ω [t,t],r) + xx V L Ω [t,t],r)) ψ,t ) L Ω,C), t [t,t ]. Altogether, this yields the following bounds in the associated Sobolev spaces ψ,t) H Ω,C) + x V L Ω [t,t],r) )) ψ,t t t ) H Ω,C), ψ,t) H Ω,C) +C τ, x V L Ω [t,t],r), xx V ) L t t Ω [t,t],r) )) ψ,t ) H Ω,C), t t τ, t [t,t ], which are in accordance with 4.). The extension to higher-order derivatives and higher space dimensions is straightforward. Stability in Sobolev spaces Numerical evolution operator). Similar considerations permit to deduce a bound for the numerical evolution operator defining a CFQM exponential integrator.); more precisely, provided that the method coefficients are real, the analogous relations S τn,t n ) H H =, S τ n,t n ) H m H m +C m τ n, ) V W,m [t n τ,t n+ ],Ω,R) n, τ n τ n, m N, 4.a) n {,,...,N },
9 CONVERGENCE ANALYSIS OF CFQM EXPONENTIAL INTEGRATORS 9 of hold with a constant that depends on an upper bound for the time stepsize, uniform bounds for certain space derivatives of the potential, and the method coefficients. By means of the relation + x e x for x, we obtain a stability estimate for compositions N n= S τ n,t n ) H m H m N e C mτ n, V W,m [tn,tn+ ],Ω,R) )τn e C mτ, V W,m [t,t ],Ω,R) )T t ) n= C, max { τ n : n {,,...,N } } τ, m N, 4.b) with constant depending in particular on an upper bound for the time stepsizes, the length of the time interval [t,t ], and uniform bounds for certain space derivatives of the potential. For a space-time-dependent Hamilton operator of the special structure 4.), the associated timeindependent operators rewrite as B j τ n,t n ) = K k= a jk At n + c k τ n ) = i b j + K k= j {,...,J}, n {,...,N }. ) a jk V t n + c k τ n ), If a jk R for any j,k) {,...,J} {,...,K}, Stone s Theorem ensures unitarity of the arising evolution operators with respect to the norm of the underlying Lebesgue space e σb jτ n,t n ) H H =, σ R, j {,...,J}, n {,...,N }, 4.4a) and thus unconditional stability of the numerical method follows S τn,t n ) H H =, n {,...,N }, see also.4) and compare with 4.). Slight modifications of the arguments used in the exemplification of 4.) imply the validity of the following estimate on bounded time intervals e σb j τ n,t n ) H m H m +C m τ n, ) V W,m [t n τ,t n+ ],Ω,R) n, σ τ n τ n, 4.4b) j {,...,J}, n {,,...,N }, m N, and justify the stated result 4.). 4. Local error Approach. In the following, we deduce suitable local error expansions for CFQM exponential integrators.) applied to.), making use of the special structure of the Hamilton operator 4.). A certain complexity in the treatment of higher-order schemes concerns the correct handling of the arising product of exponentials, since tools that are helpful for matrix exponentials such as the Baker Campbell Hausdorff BCH) formula involve infinite series and thus cannot be applied in the context of unbounded operators. Our approach in the spirit of LUNARDI 995) relies on linearisations of the evolution equations for the exact and numerical solutions as well as representations by the variationof-constants formula; together with Taylor series expansions of parts of the integrands and bounds for naturally arising iterated commutators, the stated result is obtained.
10 of S. BLANES, F. CASAS, C. GONZÁLEZ, M. THALHAMMER Fourth-order scheme. Due to the fact that the treatment of a general high-order CFQM exponential integrator would involve a considerable amount of technicalities, we provide a rigorous derivation of the local error expansion for the fourth-order scheme based on two Gaussian nodes, see.) and.); it suffices to study the error of.) in the first time step τ = τ >, L τ,t ) = S τ,t ) E τ,t ), S τ,t ) = e τ B τ,t ) e τ B τ,t ), B j τ,t ) = i b j + a j V t + c τ) + a j V t + c τ) ), j {,}, α = 6, c = α, c = + α, a = a = 4 + α, a = a = 4 α, b = a + a = b = a + a =. 4.5) Simplified notation. Henceforth, we specify the arising bounds for space and time derivatives of the potential, but we suppress the dependence of constants on upper bounds for the time stepsize; this slightly simplified notation is also justified by practical implementations, where commonly τ, ). For instance, we use that the relation +C m τ, V ) W,m τ C [t,t +τ],ω,r) m V ) W,m 4.6) [t,t +τ],ω,r) is valid with generic constant, see 4.) and 4.). Auxiliary abbreviations and estimates. In view of linearisations of the underlying evolution equations for the exact and numerical solutions about the midpoint of the first subinterval and representations by the variation-of-constants formula, we denote A = A t + τ ) = i +V t + τ ) ), Rt) = At) A, t [t,t + τ], S j τ,t ) = b j B j τ,t ) A = a j Rt + c τ) + a j Rt + c τ) ) 4.7a), j {,}. We make use of the fact that the remainder is determined by values of the potential such that the following auxiliary estimate holds with generic constant Rt) = i V t) V t + τ ) ), Rt) H m H m + S j τ,t ) H m H m C m V W,m [t,t +τ],ω,r), 4.7b) t [t,t + τ], j {,}, m N ; in certain places, to prove that a term has the correct order with respect to the time stepsize, we instead use an integral representation involving the first time derivative t Rt) = i V ζ ) dζ, t + τ Rt) H m H m + S j τ,t ) H m H m C m V W,m [t,t +τ],ω,r) τ, Evolution equations and linearisations. composed by two exponentials t [t,t + τ], j {,}, m N. 4.7c) The numerical evolution operator associated with 4.5) is S τ,t ) = S τ,τ,t )S τ,τ,t ), S j σ,τ,t ) = e σ B jτ,t ), σ R, j {,}, 4.8a)
11 CONVERGENCE ANALYSIS OF CFQM EXPONENTIAL INTEGRATORS of which satisfy autonomous linear evolution equations d dσ S jσ,τ,t ) = B j τ,t )S j σ,τ,t ) = A S j σ,τ,t ) + S jτ,t )S j σ,τ,t ), σ,τ), j {,}. For the exact evolution operator, we employ the analogous identity E τ,t ) = E τ,t + τ ) E τ ),t 4.8b) 4.9a) and make use of the fact that the associated evolution equation rewrites as d dσ E σ,ϑ ) = Aϑ + σ)e σ,ϑ ) = A E σ,ϑ ) + Rϑ + σ)e σ,ϑ ), σ, τ ), ϑ { t,t + τ } 4.9b) ; evidently, a linear transformation yields the equivalent formulation d dσ E σ,ϑ ) = A E σ,ϑ ) + Rϑ + σ)e σ,ϑ ), σ,τ), ϑ { t,t + τ } 4.9c). Auxiliary initial value problem. value problem The evolution equations in 4.8) and 4.9) suggest to study the initial { d dσ Y σ,z,z) = Z Y σ,z,z) + zσ)y σ,z,z), σ,σ ), Y,Z,z) = I, 4.a) and to distinguish the two cases { Z = A, z = S jτ,t ), j {,}, σ = τ, Z = A, zσ) = Rϑ + σ), ϑ { t,t + τ }, σ [,σ ], σ = τ ; 4.b) for the sake of clarity, we indicate the dependence of the solution on the defining linear operator and the defining function. Our previous considerations imply relations of the form { Cm V zσ) H m H m W,m [t,t +τ],ω,r), C m V W,m [t,t +τ],ω,r) τ, m N, e σ Z H H =, Y σ,z,z) H H =, e σ Z V ) H m H m C W m,m, [t,t +τ],ω,r) Y σ,z,z) V ) H m H m C W m,m, [t,t +τ],ω,r) σ, σ [,σ ], m N, with generic constants, see Section 4., and recall 4.6) as well as 4.7). 4.) Solution representation and first expansion. In view of 4.5), we deduce a stepwise expansion of the solution to 4.). Our starting point is a representation based on the variation-of-constants formula Y σ,z,z) = e σ Z + R ) σ,z,z), R ) σ,z,z) = σ e σ σ )Z zσ )Y σ,z,z) dσ, 4.)
12 of S. BLANES, F. CASAS, C. GONZÁLEZ, M. THALHAMMER involving the exponential and a remainder in integral form. The relations in 4.) imply that the following bound holds with p = R p+) σ,z,z) H H C V ) W, τ p+, 4.) [t,t +τ],ω,r) which explains the role of the upper index. By inserting the analogous identity σ Y σ,z,z) = e σ Z + e σ σ )Z zσ )Y σ,z,z) dσ, σ [,σ ], into 4.), we obtain the alternative expansion Y σ,z,z) = e σ Z + I ) σ,z,z) + R 4) σ,z,z); 4.4a) the principal term and the remainder are given by I ) R 4) σ,z,z) = e σ Z σ and fulfill the estimate σ σ,z,z) = e σ Z e σ Z zσ )e σ Z dσ, σ I l) σ,z,z) H H e σ Z zσ )e σ Z e σ Z zσ )Y σ,z,z) dσ dσ, V ) C W, τ l, [t,t +τ],ω,r) C V ) W, τ l, [t,t +τ],ω,r) 4.4b) 4.5) with l = as well as 4.) with p =. In the context of the fourth-order scheme 4.5), we additionally employ the expansion Y σ,z,z) = e σ Z + I ) σ,z,z) + I 4) σ,z,z) + R 6) σ,z,z); 4.6a) the principal term and the remainder I 4) σ,z,z) = e σ Z σ R 6) σ,z,z) = e σ Z σ σ σ σ e σ Z zσ )e σ Z e σ Z zσ )e σ Z dσ dσ, e σ Z zσ )e σ Z e σ Z zσ )e σ Z e σ Z zσ )Y σ,z,z) dσ dσ dσ, 4.6b) satisfy 4.5) with l = 4 and 4.) with p = 4, respectively. Accordingly, for the local error analysis of a general pth-order CFQM exponential integrator, an auxiliary expansion with remainder fulfilling 4.) is used. First expansion for product. twofold product In view of 4.8) and 4.9), we specify the expansion obtained for a Y σ,z,z )Y σ,z,z ) = I ) σ,z,z,z ) + R 6) σ,z,z,z ) 4.7a)
13 CONVERGENCE ANALYSIS OF CFQM EXPONENTIAL INTEGRATORS of that involves a linear operator Z and two functions z,z satisfying 4.). A straightforward calculation shows that the principal term is given by I ) σ,z,z,z ) = e σ Z + e σ Z I ) σ,z,z ) + I ) σ,z,z )e σ Z the remainder comprises several terms + e σ Z I 4) σ,z,z ) + I 4) σ,z,z )e σ Z + I ) σ,z,z )I ) σ,z,z ); R 6) σ,z,z,z ) = R 6) σ,z,z )Y σ,z,z ) + I 4) σ,z,z )R ) σ,z,z ) + I ) σ,z,z )R 4) σ,z,z ) + e σ Z R 6) σ,z,z ), which have the desired order with respect to the time step size R 6) σ,z,z,z ) H H C V ) W, τ 6, [t,t +τ],ω,r) see 4.) and 4.5). 4.7b) 4.7c) 4.7d) Taylor series expansion. With regard to the integrands in 4.4) and 4.6), we perform a Taylor series expansion of a function of the form gσ) = e σ Z Ge σ Z, Z = α +V t + τ ) ), G = α V t) V t + τ ) ), 4.8a) σ R, α,α C, t [t,t + τ]; by means of the following compact notation for iterated commutators ad Z G) = G, ad Z G) = ad Z G) = [ G,Z ] = GZ Z G, the values of the function and its derivatives read as ad l Z G) = [ ad l Z G),Z ], l N, g l) σ) = e σ Z ad l Z G)e σ Z, σ R, l N. 4.8b) For the convenience of the reader, we recall the Taylor series expansion of order m N with remainder in integral form gσ) = m l= l! σ σ ) l g l) σ ) + m! σ σ ) m+ θ) m g m+) σ + θ σ σ ) ) dθ, σ,σ R; the specification to 4.8) leads to the relation e σ Z Ge σ Z = e σ Z m l= l! σ σ ) l ad l Z G) ) + m! σ σ ) m+ θ) m e θ σ σ )Z ad m+ G)e θ σ σ )Z dθ e σ Z, σ R. Z
14 4 of S. BLANES, F. CASAS, C. GONZÁLEZ, M. THALHAMMER Appropriate choices of the center σ R are given below. Commutators. In order to characterise the regularity requirements that apply accordingly to the above expansion, it remains to study the principal terms in the iterated commutators; as the extension to higher space dimensions is straightforward, it suffices to consider a single space dimension. A simple calculation shows that the first commutator of the Laplacian and a multiplication operator is given by a first-order differential operator involving coefficients that depend on the first and second derivatives of the multiplication operator; more precisely, we have [ xx,v t) ] v L Ω,C) x V t) x L v Ω,R) L Ω,C) + xx V t) L Ω,R) vl Ω,C), v H Ω,C), t [t,t + τ]; the second commutator satisfies the relation [ xx, [ xx,v t) ]] v L Ω,C) xx V t) L Ω,R) xx v L Ω,C) + 4 xxx V t) L Ω,R) x v L Ω,C) + xxxx V t) L v Ω,R) L Ω,C), v H Ω,C), t [t,t + τ]. By induction, we obtain a bound of the form involving a generic constant. ad m Z V t) V t + τ )) H H m t [t,t + τ], m N, { Cm V W,m [t,t +τ],ω,r), C m V W,m [t,t +τ],ω,r) τ, Expansion of principal term. In order to deduce a suitable representation of the principal term in 4.7), we apply the above stated relations. More precisely, we make use of the following Taylor series expansions e σ Z zσ)e σ Z = e σ Z zσ)e σ Z + r ) σ,σ,z,z), with remainders satisfying the bounds = e σ Z zσ) + σ σ )ad Z zσ) ) + σ σ ) ad Z zσ) ) ) e σ Z + r 4) σ,σ,z,z), σ [,σ ], r ) σ,σ,z,z) = σ σ ) e σ ) +θ σ σ ))Z ad Z zσ) e σ +θ σ σ ))Z dθ, r ) σ,σ,z,z) H H C V W, [t,t +τ],ω,r) τ, r 4) σ,σ,z,z) = σ σ ) θ) e σ +θ σ σ ))Z ad ) Z zσ) e σ +θ σ σ ))Z dθ, r 4) σ,σ,z,z) H H C V W,6 [t,t +τ],ω,r) τ4, σ [,σ ]. On the one hand, integration implies I ) σ,σ,z,z) = I ) σ,σ,z,z) + R ) σ,σ,z,z) = I ) σ,σ,z,z) + I ) σ,σ,z,z) + I 4) σ,σ,z,z) + R 5) σ,σ,z,z),
15 CONVERGENCE ANALYSIS OF CFQM EXPONENTIAL INTEGRATORS 5 of where the principal terms are given by I ) I ) σ σ,σ,z,z) = e σ σ )Z zσ ) dσ e σ Z, σ σ,σ,z,z) = e σ σ )Z σ σ )ad Z zσ ) ) dσ e σ Z, I 4) σ,σ,z,z) = e σ σ )Z σ and the remainders have the desired orders σ σ ) ad Z zσ ) ) dσ e σ Z, R ) σ σ,σ,z,z) = e σ Z r ) σ,σ,z,z) dσ, ) R σ,σ,z,z) H H C V W, [t,t +τ],ω,r) τ, R 5) σ σ,σ,z,z) = e σ Z r 4) σ,σ,z,z) dσ, 5) R σ,σ,z,z) H H C V W,6 [t,t +τ],ω,r) τ5, see 4.4); on the other hand, we obtain the following representation I 4) σ,σ,z,z) = I 4) 4 σ,σ,z,z) + R 5) 4 σ,σ,z,z), I 4) 4 σ,σ,z,z) = e σ σ )Z σ σ zσ ) zσ ) dσ dσ e σ Z, involving a remainder R 5) σ σ 4 σ,σ,z,z) = e σ Z r ) σ,σ,z,z)e σ Z zσ )e σ Z ) + e σ Z zσ )e σ Z r ) σ,σ,z,z) dσ dσ that fulfills an estimate of the form R 5) 4 σ,σ,z,z) H H C V ) W, τ 5 [t,t +τ],ω,r) with generic constant, see 4.) and 4.6). Choices of centers. As described below, suitable choices of the centers in the Taylor series expansions permit to extract a common exponential factor; accordingly, we rewrite the expansion 4.7) as Y σ,z,z )Y σ,z,z ) = I ) 5 σ,ζ,z,z,z ) + R 5) 5 σ,ζ,z,z,z ) with principal term given by I ) 5 σ,ζ,z,z,z ) = e σ Z + e σ Z I ) σ,ζ,z,z ) + I ) σ,ζ,z,z ) + I 4) σ,ζ,z,z ) ) + I ) σ,ζ,z,z ) + I ) σ,ζ,z,z ) + I 4) σ,ζ,z,z ) ) e σ Z + e σ Z I 4) 4 σ,ζ,z,z ) + I 4) 4 σ,ζ 4,Z,z )e σ Z + I ) σ,ζ 5,Z,z )I ) σ,ζ 6,Z,z );
16 6 of S. BLANES, F. CASAS, C. GONZÁLEZ, M. THALHAMMER the remainder comprises the terms R 5) 5 σ,ζ,z,z,z ) = R 6) σ,z,z,z ) + e σ Z R 5) σ,ζ,z,z ) + R 5) σ,ζ,z,z )e σ Z + e σ Z R 5) 4 σ,ζ,z,z ) + R 5) 4 σ,ζ 4,Z,z )e σ Z + R ) σ,ζ 5,Z,z )I ) σ,z,z ) + I ) σ,ζ 5,Z,z )R ) σ,ζ 6,Z,z ) and thus satisfies an estimate of the form 5) R 5 σ,ζ,z,z,z ) V ) H H C W,6 τ 5, [t,t +τ],ω,r) see again 4.). Local error expansion. By means of our previous considerations, we obtain a suitable local error expansion for the fourth-order scheme 4.5); more precisely, based on the identity L τ,t ) = S τ,t ) E τ,t ) = S τ,τ,t )S τ,τ,t ) E τ,t + τ ) E τ ),t = Y τ, A, S τ,t ) ) Y τ, A, S τ,t ) ) Y τ,a,r t + τ + ))) Y τ,a,r t + ) )) = Y τ, A, S τ,t ) ) Y τ, A, S τ,t ) ) Y τ, A, R t + τ + ))) Y τ, A, R t + ))), see 4.8), 4.9), and 4.), we arrive at the relation The first summand comprises several terms I ) L τ,t ) = I ) 6 τ,t,ζ,a,v ) + R 5) 6 τ,t,ζ,a,v ). 6 τ,t,ζ,a,v ) = I ) 5 τ,ζ, A, S τ,t ), S τ,t ) ) I ) 5 τ,ζ, A, R t + τ + )), R t + ))) = I ) 6 τ,t,ζ,ζ,a,v ) + I ) 6 τ,t,ζ,ζ,a,v ) + I ) 6 τ,t,ζ,ζ,a,v ) + I 4) 64 τ,t,ζ,ζ 4,ζ 5,ζ 6,A,V ) I 4) 65 τ,t,ζ,ζ 4,ζ 5,ζ 6,A,V ) which have the desired order four, as we justify in detail below; the remainder R 5) 6 τ,t,ζ,a,v ) = R 5) 5 τ,ζ, A, S τ,t ), S τ,t ) ) R 5) 5 τ,ζ, A, R t + τ + )), R t + ))) is well-defined when considered as operator from H to H and satisfies the bound R 5) 6 τ,t,ζ,a,v ) H H C V ) W,6 τ 5. [t,t +τ],ω,r)
17 CONVERGENCE ANALYSIS OF CFQM EXPONENTIAL INTEGRATORS 7 of As indicated before, we make use of the fact that the particular choices ζ =, ζ = τ, ζ = τ, ζ 4 =, ζ 5 =, ζ 6 = τ, of the centers in the Taylor series expansions employed for the numerical and exact evolution operators permit to extract common exponential factors. On the one hand, we obtain the relations I ) I ) 6 τ,t,ζ,ζ,a,v ) = e τ A I ) τ,ζ, A, S τ,t ) ) I ) τ,ζ, A, R t + )))) + τ,ζ, A, S τ,t ) ) I ) τ,ζ, A, R t + τ + )))) e τ A I ) = eτ A I τ,t,v ), τ I τ,t,v ) = S τ,t ) R t + σ ) ) dσ + I ) τ S τ,t ) R t + τ + σ ) ) dσ, 6 τ,t,ζ,ζ,a,v ) = e τ A I ) τ,ζ, A, S τ,t ) ) I ) τ,ζ, A, R t + )))) + τ,ζ, A, S τ,t ) ) I ) τ,ζ, A, R t + τ + )))) e τ A = 4 eτ A ad A I τ,t,v ) ), I τ,t,v ) = I ) τ I 4) σ S τ,t ) R t + σ ) ) dσ + τ σ + τ) S τ,t ) R t + τ + σ ) ) dσ, 6 τ,t,ζ,ζ,a,v ) = e τ A I 4) τ,ζ, A, S τ,t ) ) I 4) τ,ζ, A, R t + )))) + τ,ζ, A, S τ,t ) ) I 4) τ,ζ, A, R t + τ + )))) e τ A = 6 eτ A ad A I τ,t,v ) ), τ I τ,t,v ) = S τ,t ) R t + σ ) ) dσ + σ on the other hand, we have I 4) 64 τ,t,ζ,ζ 4,ζ 5,ζ 6,A,V ) = e τ A I 4) 4 τ,ζ, A, S τ,t ) ) + I 4) 4 τ,ζ4, A, S τ,t ) ) e τ A + I ) τ,ζ5, A, S τ,t ) ) I ) τ,ζ6, A, S τ,t ) ) = 4 e τ A I 4 τ,t,v )e τ A, I 4 τ,t,v ) = τ σ τ σ + τ) S τ,t ) R t + τ + σ ) ) dσ ; S τ,t )S τ,t ) + S τ,t )S τ,t ) ) τ τ dσ dσ + S τ,t ) dσ S τ,t ) dσ,
18 8 of S. BLANES, F. CASAS, C. GONZÁLEZ, M. THALHAMMER I 4) 65 τ,t,ζ,ζ 4,ζ 5,ζ 6,A,V ) = e τ A I 4) 4 τ,ζ, A, R t + ))) + I 4) 4 τ,ζ 4, A, R t + τ + ))) e τ A + I ) τ,ζ 5, A, R t + τ + ))) I ) τ,ζ 6, A, R t + ))) = 4 e τ A I 5 τ,t,v )e τ A, I 5 τ,t,v ) = τ σ + τ R t + σ ) R t + σ ) + R t + τ + σ ) R t + τ + σ ) ) dσ dσ R t + τ + σ ) τ dσ R t + σ ) dσ. Straightforward Taylor series expansions of the arising integrands lead to the auxiliary estimates I τ,t,v ) H H C V W 4, [t,t +τ],ω,r) τ5, ad A I τ,t,v ) ) H H C V W, [t,t +τ],ω,r) τ5, ad A I τ,t,v ) ) H H C V W,4 [t,t +τ],ω,r) τ5, I4 τ,t,v ) H H + I5 τ,t,v ) H H C V W, [t,t +τ],ω,r) τ6, see also 4.7), and at once imply ) I 6 τ,t,ζ,a,v ) H H C V W 4, [t,t +τ],ω,r), V W, [t,t +τ],ω,r), V ) W,4 τ 5. [t,t +τ],ω,r) As a consequence, we arrive at the following result L τ,t ) H H C max l {,,,4} ) V W l,4 l) τ 5. [t,t +τ],ω,r) Local error estimate. Altogether, this proves that the fourth-order CFQM exponential integrator.) satisfies the local error bound S τn,t n ) E τ n,t n ) ) H H p C p max V W l {,...,p} l,p l) [t n τ p+,t n+ ],Ω,R) n C τn p+, 4.9) n {,,...,N }, with p = 4 and in particular explains the dependence of the constant arising in Theorem on certain space and time derivatives of the potential. Exponential midpoint rule. It is straightforward to apply the above approach in the context of the exponential midpoint rule; as the scheme involves a single exponential, the calculations are considerably simplified. 5. General case For the sake of completeness, we include the extension of Theorem for Hamilton operators without special structure. The derivation of the local error expansion based on linearisations of the underlying
19 REFERENCES 9 of evolution equations, representations by the variation-of-constants formula, and Taylor series expansions of parts of the integrands in principle applies to the general case as well, see Section 4. ; however, for the specification of regularity assumptions on the problem data, more detailed information on the form of the considered Hamiltonians would be required. Theorem Assume that the operator familiy defining the non-autonomous linear Schrödinger equation.) and its exact solution are sufficiently regular such that compositions of the form A l) s)u m) t) with l {,,..., p} and m {,,..., p } remain bounded. Then, any pth-order CFQM exponential integrator.) satisfies a global error estimate of the form un ut ) H u ut ) H +C τ p max ; the arising constant depends on the data of the problem, the length of the considered time interval, and the method coefficients, but it is independent of the number of time steps n and the maximal time stepsize. Illustrations. Numerical experiments that confirm the global error bound of Theorem and Theorem, respectively, are found in BLANES, CASAS, THALHAMMER 7, 8), see also BADER, BLANES, KOPYLOV 8); due to the fact that the exponential midpoint rule is an instance of a Magnus integrator and a CFQM exponential integrator, we also refer to the numerical example given in HOCHBRUCK, LUBICH ). For vanishing dissipation parameter δ =, the evolution equation considered in BLANES, CASAS, THALHAMMER 7) is of Schrödinger type; the presented results in particular verify the order of convergence and provide numerical evidence that CFQM exponential integrators are superior in efficiency compared to explicit Runge Kutta methods and Magnus integrators for differential equations involving a higher degree of freedom. Acknowledgements Part of this work was developed during a research stay at the Wolfgang Pauli Institute Vienna; the authors are grateful to the director Norbert Mauser and the staff members for their support and hospitality. The first three authors acknowledge funding by the Ministerio de Economía y Competitividad Spain) through project MTM P AEI/FEDER, UE). References A. ALVERMANN, H. FEHSKE, High-order commutator-free exponential time-propagation of driven quantum systems. J. Comp. Phys. /5 ) A. ALVERMANN, H. FEHSKE, P. B. LITTLEWOOD, Numerical time propagation of quantum systems in radiation fields. New. J. Phys. 4 ) 58 pp). PH. BADER, S. BLANES, N. KOPYLOV, Exponential propagators for the Schrödinger equation with a timedependent potential. J. Chem. Phys. 48 8) 449. S. BLANES, F. CASAS, J. A. OTEO, J. ROS, The Magnus expansion and some of its applications. Physics Reports 47/5-6 9) 5 8. S. BLANES, F. CASAS, M. THALHAMMER, High-order commutator-free quasi-magnus exponential integrators for non-autonomous linear evolution equations. Comp. Phys. Commun. 7) 4 6. S. BLANES, F. CASAS, M. THALHAMMER, Convergence analysis of high-order commutator-free quasi-magnus exponential integrators for nonautonomous linear evolution equations of parabolic type. IMA J. Numer. Anal. 8/ 8)
20 of REFERENCES S. BLANES, P. C. MOAN, Fourth- and sixth-order commutator-free Magnus integrators for linear and non-linear dynamical systems. App. Num. Math. 56/ 6) P. BONGAARTS, Quantum Theory. A Mathematical Approach. Springer, Cham, 5. K. J. ENGEL, R. NAGEL, One-Parameter Semigroups for Linear Evolution Equations. Springer, New York,. M. HOCHBRUCK, CH. LUBICH, On Magnus integrators for time-dependent Schrödinger equations. SIAM J. Numer. Anal. 4/ ) A. ISERLES, K. KROPIELNICKA, P. SINGH, Compact schemes for laser-matter interaction in Schrödinger equation based on effective splittings of Magnus expansion. Comp. Phys. Commun. 4 9) 95. E. KYOSEVA, N. VITANOVA, B. SHORE, Physical realization of coupled Hilbert-space mirrors for quantum-state engineering. J. Modern Optics 54 7) A. LUNARDI, Analytic Semigroups and Optimal Regularity in Parabolic Problems. Birkhäuser, Basel, 995.
Favourable time integration methods for. Non-autonomous evolution equations
Favourable time integration methods for non-autonomous evolution equations International Conference on Scientific Computation and Differential Equations (SciCADE) Session on Numerical integration of non-autonomous
More 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 information1. Introduction. In the present paper, we consider a nonautonomous differential equation involving a time-dependent linear operator A
SIAM J. NUMER. ANAL. Vol. 44, No. 2, pp. 851 864 c 2006 Society for Industrial and Applied Mathematics A FOURTH-ORDER COMMUTATOR-FREE EXPONENTIAL INTEGRATOR FOR NONAUTONOMOUS DIFFERENTIAL EQUATIONS MECHTHILD
More informationSymplectic propagators for the Kepler problem with time-dependent mass
Celestial Mechanics and Dynamical Astronomy manuscript No. (will be inserted by the editor) Symplectic propagators for the Kepler problem with time-dependent mass Philipp Bader Sergio Blanes Fernando Casas
More informationSymplectic time-average propagators for the Schödinger equation with a time-dependent Hamiltonian
Symplectic time-average propagators for the Schödinger equation Symplectic time-average propagators for the Schödinger equation with a time-dependent Hamiltonian Sergio Blanes,, a Fernando Casas, 2, b
More informationDIMENSION SPLITTING FOR TIME DEPENDENT OPERATORS. Eskil Hansen and Alexander Ostermann
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume X, Number 0X, XX 00X pp. 1 11 DIMENSION SPLITTING FOR TIME DEPENDENT OPERATORS Eskil Hansen and Alexander Ostermann Institut
More informationSANGRADO PAGINA 17CMX24CM. PhD Thesis. Splitting methods for autonomous and non-autonomous perturbed equations LOMO A AJUSTAR (AHORA 4CM)
SANGRADO A3 PAGINA 17CMX24CM LOMO A AJUSTAR (AHORA 4CM) We have considered the numerical integration of non-autonomous separable parabolic equations using high order splitting methods with complex coefficients.
More informationA posteriori error estimation for Magnus-type integrators
ASC Report No. /28 A posteriori error estimation for Magnus-type integrators W. Auzinger, H. Hofsta tter, O. Koch, M. Quell and M. Thalhammer Institute for Analysis and Scientific Computing Vienna University
More informationA pedagogical approach to Magnus expansion
A pedagogical approach to Magnus expansion S Blanes F Casas J A Oteo and J. Ros Instituto de Matemática Multidisciplinar, Universidad Politécnica de Valencia, E-46022 Valencia, Spain Departament de Matemàtiques,
More informationCOMPUTING THE MATRIX EXPONENTIAL WITH AN OPTIMIZED TAYLOR POLYNOMIAL APPROXIMATION
COMPUTING THE MATRIX EXPONENTIAL WITH AN OPTIMIZED TAYLOR POLYNOMIAL APPROXIMATION PHILIPP BADER, SERGIO BLANES, FERNANDO CASAS Abstract. We present a new way to compute the Taylor polynomial of the matrix
More informationBIHARMONIC WAVE MAPS INTO SPHERES
BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.
More informationS. Sánchez, F. Casas & A. Fernández
New analytic approximations based on the Magnus expansion S. Sánchez, F. Casas & A. Fernández Journal of Mathematical Chemistry ISSN 259-9791 Volume 49 Number 8 J Math Chem (211) 49:1741-1758 DOI 1.17/s191-11-9855-
More informationSplitting and composition methods for the time dependent Schrödinger equation
Splitting and composition methods for the time dependent Schrödinger equation S. Blanes, joint work with F. Casas and A. Murua Instituto de Matemática Multidisciplinar Universidad Politécnica de Valencia,
More informationSplittings, commutators and the linear Schrödinger equation
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,
More informationPhilippe Chartier 1, Florian Méhats 2, Mechthild Thalhammer 3 and Yong Zhang Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher CONVERGENCE OF MULTI-REVOLUTION COMPOSITION TIME-SPLITTING METHODS FOR HIGHLY OSCILLATORY
More informationHigh-order actions and their applications to honor our friend and
High-order actions and their applications to honor our friend and collaborator Siu A. Chin Universitat Politecnica de Catalunya, Barcelona, Spain Lecture: Magnus Expansion and Suzuki s Method Outline of
More informationLie algebraic quantum control mechanism analysis
Lie algebraic quantum control mechanism analysis June 4, 2 Existing methods for quantum control mechanism identification do not exploit analytic features of control systems to delineate mechanism (or,
More informationExponential integrators
Exponential integrators Marlis Hochbruck Heinrich-Heine University Düsseldorf, Germany Alexander Ostermann University Innsbruck, Austria Helsinki, May 25 p.1 Outline Time dependent Schrödinger equations
More informationSchrödinger equation
Splitting methods for the time dependent Schrödinger equation S. Blanes, Joint work with F. Casas and A. Murua Instituto de Matemática Multidisciplinar Universidad Politécnica de Valencia, SPAIN Workshop:
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 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 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 informationBLOCH FUNCTIONS ON THE UNIT BALL OF AN INFINITE DIMENSIONAL HILBERT SPACE
BLOCH FUNCTIONS ON THE UNIT BALL OF AN INFINITE DIMENSIONAL HILBERT SPACE OSCAR BLASCO, PABLO GALINDO, AND ALEJANDRO MIRALLES Abstract. The Bloch space has been studied on the open unit disk of C and some
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 informationCOMPOSITION MAGNUS INTEGRATORS. Sergio Blanes. Lie Group Methods and Control Theory Edinburgh 28 June 1 July Joint work with P.C.
COMPOSITION MAGNUS INTEGRATORS Sergio Blanes Departament de Matemàtiques, Universitat Jaume I Castellón) sblanes@mat.uji.es Lie Group Methods and Control Theory Edinburgh 28 June 1 July 24 Joint work with
More informationEfficient 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 informationQuantum Computing Lecture 3. Principles of Quantum Mechanics. Anuj Dawar
Quantum Computing Lecture 3 Principles of Quantum Mechanics Anuj Dawar What is Quantum Mechanics? Quantum Mechanics is a framework for the development of physical theories. It is not itself a physical
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 informationarxiv: v1 [math.na] 5 Apr 2016
THE BCH-FORMULA AND ORDER CONDITIONS FOR SPLITTING METHODS arxiv:1604.01190v1 [math.na] 5 Apr 2016 WINFRIED AUZINGER, WOLFGANG HERFORT, OTHMAR KOCH, AND MECHTHILD THALHAMMER Abstract. As an application
More informationSMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract
SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm Abstract Dedicated to Prof. Michele Caputo on the occasion of his 8th birthday We consider ordinary fractional
More informationIMPROVED HIGH ORDER INTEGRATORS BASED ON MAGNUS EPANSION S. Blanes ;, F. Casas, and J. Ros 3 Abstract We build high order ecient numerical integration
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Improved High Order Integrators Based on Magnus Expansion S. Blanes, F. Casas and J. Ros DAMTP 999/NA8 August, 999 Department of Applied Mathematics and
More informationTHE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS
THE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS MARIANNA KHANAMIRYAN Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road,
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationHigh order splitting methods for analytic semigroups exist
BIT manuscript No. (will be inserted by the editor High order splitting methods for analytic semigroups exist Eskil Hansen Alexander Ostermann Received: date / Accepted: date Abstract In this paper, we
More informationMinimal periods of semilinear evolution equations with Lipschitz nonlinearity
Minimal periods of semilinear evolution equations with Lipschitz nonlinearity James C. Robinson a Alejandro Vidal-López b a Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. b Departamento
More informationarxiv: v2 [math-ph] 15 Jun 2012
Efficient computation of the Zassenhaus formula Fernando Casas 1 Ander Murua 2 Mladen Nadinic 3 ariv:1204.0389v2 [math-ph] 15 Jun 2012 June 18, 2012 Abstract A new recursive procedure to compute the Zassenhaus
More informationNOTES ON CALCULUS OF VARIATIONS. September 13, 2012
NOTES ON CALCULUS OF VARIATIONS JON JOHNSEN September 13, 212 1. The basic problem In Calculus of Variations one is given a fixed C 2 -function F (t, x, u), where F is defined for t [, t 1 ] and x, u R,
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 informationLeast Squares Based Self-Tuning Control Systems: Supplementary Notes
Least Squares Based Self-Tuning Control Systems: Supplementary Notes S. Garatti Dip. di Elettronica ed Informazione Politecnico di Milano, piazza L. da Vinci 32, 2133, Milan, Italy. Email: simone.garatti@polimi.it
More informationSwitching, sparse and averaged control
Switching, sparse and averaged control Enrique Zuazua Ikerbasque & BCAM Basque Center for Applied Mathematics Bilbao - Basque Country- Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ WG-BCAM, February
More informationUnitary transformations depending on a small parameter
Unitary transformations depending on a small parameter Fernando Casas Fernando.Casas@mat.uji.es www.gicas.uji.es Departament de Matemàtiques Institut de Matemàtiques i Aplicacions de Castelló (IMAC) Universitat
More informationAdaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear Degenerate Reaction-Diffusion Problems
EPJ Web of Conferences 108, 0019 (016) DOI: 10.1051/ epjconf/ 0161080019 C Owned by the authors, published by EDP Sciences, 016 Adaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear
More 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 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 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 information#A69 INTEGERS 13 (2013) OPTIMAL PRIMITIVE SETS WITH RESTRICTED PRIMES
#A69 INTEGERS 3 (203) OPTIMAL PRIMITIVE SETS WITH RESTRICTED PRIMES William D. Banks Department of Mathematics, University of Missouri, Columbia, Missouri bankswd@missouri.edu Greg Martin Department of
More informationarxiv: v1 [math.na] 12 Oct 2018
arxiv:181.5412v1 [math.na] 12 Oct 218 Sixth-order schemes for laser matter interaction in Schrödinger equation Pranav Singh October 15, 218 AMS Mathematics Subject Classification: Primary 65M7, Secondary
More informationExponential Runge-Kutta methods for parabolic problems
Exponential Runge-Kutta methods for parabolic problems Marlis Hochbruck a,, Alexander Ostermann b a Mathematisches Institut, Heinrich-Heine Universität Düsseldorf, Universitätsstraße 1, D-4225 Düsseldorf,
More informationSome asymptotic properties of solutions for Burgers equation in L p (R)
ARMA manuscript No. (will be inserted by the editor) Some asymptotic properties of solutions for Burgers equation in L p (R) PAULO R. ZINGANO Abstract We discuss time asymptotic properties of solutions
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationGALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS
GALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND CHARALAMBOS MAKRIDAKIS Abstract. We consider discontinuous as well as continuous Galerkin methods for the time discretization
More informationOptimal primitive sets with restricted primes
Optimal primitive sets with restricted primes arxiv:30.0948v [math.nt] 5 Jan 203 William D. Banks Department of Mathematics University of Missouri Columbia, MO 652 USA bankswd@missouri.edu Greg Martin
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 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 informationResolvent Algebras. An alternative approach to canonical quantum systems. Detlev Buchholz
Resolvent Algebras An alternative approach to canonical quantum systems Detlev Buchholz Analytical Aspects of Mathematical Physics ETH Zürich May 29, 2013 1 / 19 2 / 19 Motivation Kinematics of quantum
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationA POSTERIORI ERROR ESTIMATES FOR THE BDF2 METHOD FOR PARABOLIC EQUATIONS
A POSTERIORI ERROR ESTIMATES FOR THE BDF METHOD FOR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND PANAGIOTIS CHATZIPANTELIDIS Abstract. We derive optimal order, residual-based a posteriori error estimates for
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationSYNCHRONIZATION OF NONAUTONOMOUS DYNAMICAL SYSTEMS
Electronic Journal of Differential Equations, Vol. 003003, No. 39, pp. 1 10. ISSN: 107-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu login: ftp SYNCHRONIZATION OF
More informationA max-plus based fundamental solution for a class of infinite dimensional Riccati equations
A max-plus based fundamental solution for a class of infinite dimensional Riccati equations Peter M Dower and William M McEneaney Abstract A new fundamental solution for a specific class of infinite dimensional
More informationOn the efficiency of the Peaceman-Rachford ADI-dG method for wave-type methods
On the efficiency of the Peaceman-Rachford ADI-dG method for wave-type methods Marlis Hochbruck, Jonas Köhler CRC Preprint 2017/34 (revised), March 2018 KARLSRUHE INSTITUTE OF TECHNOLOGY KIT The Research
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationSharp Sobolev Strichartz estimates for the free Schrödinger propagator
Sharp Sobolev Strichartz estimates for the free Schrödinger propagator Neal Bez, Chris Jeavons and Nikolaos Pattakos Abstract. We consider gaussian extremisability of sharp linear Sobolev Strichartz estimates
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 informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
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 informationOptimal shape and position of the support for the internal exact control of a string
Optimal shape and position of the support for the internal exact control of a string Francisco Periago Abstract In this paper, we consider the problem of optimizing the shape and position of the support
More informationReproducing formulas associated with symbols
Reproducing formulas associated with symbols Filippo De Mari Ernesto De Vito Università di Genova, Italy Modern Methods of Time-Frequency Analysis II Workshop on Applied Coorbit space theory September
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
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 informationGlobal Solutions for a Nonlinear Wave Equation with the p-laplacian Operator
Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator Hongjun Gao Institute of Applied Physics and Computational Mathematics 188 Beijing, China To Fu Ma Departamento de Matemática
More informationHOPF S DECOMPOSITION AND RECURRENT SEMIGROUPS. Josef Teichmann
HOPF S DECOMPOSITION AND RECURRENT SEMIGROUPS Josef Teichmann Abstract. Some results of ergodic theory are generalized in the setting of Banach lattices, namely Hopf s maximal ergodic inequality and the
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationFactorization of p-dominated Polynomials through L p -Spaces
Michigan Math. J. 63 (2014), 345 353 Factorization of p-dominated Polynomials through L p -Spaces P. Rueda & E. A. Sánchez Pérez Introduction Since Pietsch s seminal paper [32], the study of the nonlinear
More informationarxiv:physics/ v1 [math-ph] 17 May 1997
arxiv:physics/975v1 [math-ph] 17 May 1997 Quasi-Exactly Solvable Time-Dependent Potentials Federico Finkel ) Departamento de Física Teórica II Universidad Complutense Madrid 84 SPAIN Abstract Niky Kamran
More informationTHE Hamilton equations of motion constitute a system
Proceedings of the World Congress on Engineering 0 Vol I WCE 0, July 4-6, 0, London, U.K. Systematic Improvement of Splitting Methods for the Hamilton Equations Asif Mushtaq, Anne Kværnø, and Kåre Olaussen
More informationThe Redner Ben-Avraham Kahng coagulation system with constant coefficients: the finite dimensional case
The Redner Ben-Avraham Kahng coagulation system with constant coefficients: the finite dimensional case F.P. da Costa, J.T. Pinto and R. Sasportes Abstract. We study the behaviour as t of solutions (c
More informationOn the bang-bang property of time optimal controls for infinite dimensional linear systems
On the bang-bang property of time optimal controls for infinite dimensional linear systems Marius Tucsnak Université de Lorraine Paris, 6 janvier 2012 Notation and problem statement (I) Notation: X (the
More information2.1 Green Functions in Quantum Mechanics
Chapter 2 Green Functions and Observables 2.1 Green Functions in Quantum Mechanics We will be interested in studying the properties of the ground state of a quantum mechanical many particle system. We
More informationMulti-product expansion, Suzuki s method and the Magnus integrator for solving time-dependent problems
Multi-product expansion, Suzuki s method and the Magnus integrator for solving time-dependent problems Jürgen Geiser and Siu A. Chin geiser@mathematik.hu-berlin.de Abstract. In this paper we discuss the
More informationLinearly dependent vectorial decomposition of clutters
Linearly dependent vectorial decomposition of clutters Jaume Martí-Farré 1,2 Departament de Matemàtica Aplicada IV Universitat Politècnica de Catalunya, BarcelonaTech Barcelona, Spain Abstract This paper
More informationFeynman s path integral approach to quantum physics and its relativistic generalization
Feynman s path integral approach to quantum physics and its relativistic generalization Jürgen Struckmeier j.struckmeier@gsi.de, www.gsi.de/ struck Vortrag im Rahmen des Winterseminars Aktuelle Probleme
More informationPropagators for TDDFT
Propagators for TDDFT A. Castro 1 M. A. L. Marques 2 A. Rubio 3 1 Institut für Theoretische Physik, Fachbereich Physik der Freie Universität Berlin 2 Departamento de Física, Universidade de Coimbra, Portugal.
More informationAn inverse scattering problem for short-range systems in a time-periodic electric field. François Nicoleau
An inverse scattering problem for short-range systems in a time-periodic electric field. François Nicoleau Laboratoire Jean Leray UMR CNRS-UN 6629 Département de Mathématiques 2, rue de la Houssinière
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationNonarchimedean Cantor set and string
J fixed point theory appl Online First c 2008 Birkhäuser Verlag Basel/Switzerland DOI 101007/s11784-008-0062-9 Journal of Fixed Point Theory and Applications Nonarchimedean Cantor set and string Michel
More informationPERIODIC SOLUTIONS WITH NONCONSTANT SIGN IN ABEL EQUATIONS OF THE SECOND KIND
PERIODIC SOLUTIONS WITH NONCONSTANT SIGN IN ABEL EQUATIONS OF THE SECOND KIND JOSEP M. OLM, XAVIER ROS-OTON, AND TERE M. SEARA Abstract. The study of periodic solutions with constant sign in the Abel equation
More informationFRACTIONAL ERROR ESTIMATES OF SPLITTING SCHEMES FOR THE NONLINEAR SCHRÖDINGER EQUATION
FRACTIONAL ERROR ESTIMATES OF SPLITTING SCHEMES FOR THE NONLINEAR SCHRÖDINGER EQUATION JOHANNES EILINGHOFF, ROLAND SCHNAUBELT, AND KATHARINA SCHRATZ Abstract. We investigate the Lie and the Strang splitting
More 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 informationThe heat equation in time dependent domains with Neumann boundary conditions
The heat equation in time dependent domains with Neumann boundary conditions Chris Burdzy Zhen-Qing Chen John Sylvester Abstract We study the heat equation in domains in R n with insulated fast moving
More informationGreedy control. Martin Lazar University of Dubrovnik. Opatija, th Najman Conference. Joint work with E: Zuazua, UAM, Madrid
Greedy control Martin Lazar University of Dubrovnik Opatija, 2015 4th Najman Conference Joint work with E: Zuazua, UAM, Madrid Outline Parametric dependent systems Reduced basis methods Greedy control
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 informationGeneralized Schrödinger semigroups on infinite weighted graphs
Generalized Schrödinger semigroups on infinite weighted graphs Institut für Mathematik Humboldt-Universität zu Berlin QMATH 12 Berlin, September 11, 2013 B.G, Ognjen Milatovic, Francoise Truc: Generalized
More informationTrotter s product formula for projections
Trotter s product formula for projections Máté Matolcsi, Roman Shvydkoy February, 2002 Abstract The aim of this paper is to examine the convergence of Trotter s product formula when one of the C 0-semigroups
More informationAbstract In this paper, we consider bang-bang property for a kind of timevarying. time optimal control problem of null controllable heat equation.
JOTA manuscript No. (will be inserted by the editor) The Bang-Bang Property of Time-Varying Optimal Time Control for Null Controllable Heat Equation Dong-Hui Yang Bao-Zhu Guo Weihua Gui Chunhua Yang Received:
More informationON WEIGHTED INEQUALITIES FOR FRACTIONAL INTEGRALS OF RADIAL FUNCTIONS. dy, 0 < γ < n. x y
ON WEIGHTED INEQUALITIES FOR FRACTIONAL INTEGRALS OF RADIAL FUNCTIONS PABLO L. DE NÁPOLI, IRENE DRELICHMAN, AND RICARDO G. DURÁN Abstract. We prove a weighted version of the Hardy-Littlewood-Sobolev inequality
More informationUnitary Dynamics and Quantum Circuits
qitd323 Unitary Dynamics and Quantum Circuits Robert B. Griffiths Version of 20 January 2014 Contents 1 Unitary Dynamics 1 1.1 Time development operator T.................................... 1 1.2 Particular
More informationSemigroups of Operators
Lecture 11 Semigroups of Operators In tis Lecture we gater a few notions on one-parameter semigroups of linear operators, confining to te essential tools tat are needed in te sequel. As usual, X is a real
More informationComplete orbits for twist maps on the plane: the case of small twist
Complete orbits for twist maps on the plane: the case of small twist Markus Kunze 1 & Rafael Ortega 1 Universität Duisburg-Essen, Fakultät für Mathematik, D - 45117 Essen, Germany Departamento de Matemática
More information