SEMI-CLASSICAL MODELS FOR THE SCHRÖDINGER EQUATION WITH PERIODIC POTENTIALS AND BAND CROSSINGS. Lihui Chai. Shi Jin. and Qin Li
|
|
- Elinor Norman
- 5 years ago
- Views:
Transcription
1 Manuscript submitted to AIMS Journals Volume X, Number X, XX X doi:.3934/... pp. X XX SEMI-CLASSICAL MODELS FO THE SCHÖDINGE EQUATION WITH PEIODIC POTENTIALS AND BAND COSSINGS Lihui Chai Department of Mathematical Science Tsinghua University Bejing 84, China Shi Jin Department of Mathematics, Institute of Natural Sciences, and MOE Key Lab in Scientific and Engineering Computing Shanghai Jiao Tong University Shanghai 4, China and Department of Mathematics University of Wisconsin-Madison Madison, WI 5376, USA and Qin Li Department of Mathematics University of Wisconsin-Madison Madison, WI 5376, USA (Communicated by the associate editor name Abstract. The Bloch decomposition plays a fundamental role in the study of quantum mechanics and wave propagation in periodic media. Most of the homogenization theory developed for the study of high frequency or semi-classical limit for these problems assumes no crossing of the Bloch bands, resulting in classical Liouville equations in the limit along each Bloch band. In this article, we derive semi-classical models for the equation in periodic media that take into account band crossing, which is important to describe quantum transitions between Bloch bands. Our idea is still based on the Wigner transform (on the Bloch eigenfunctions, but in taking the semi-classical approimation, we retain the off-diagonal entries of the Wigner matri, which cannot be ignored near the point of band crossing. This results in coupled inhomogeneous Liouville systems that can suitably describe quantum tunneling between bands that are not well-separated. We also develop a domain decomposition method that couples these semi-classical models with the classical Liouville equations (valid away from zones of band crossing for a multiscale computation. Solutions of these models are numerically compared with those of the Schröding equation to justify the validity of these new models for band-crossings. Mathematics Subject Classification. Primary: 58F5, 58F7; Secondary: 53C35. Key words and phrases. equation, Periodic, Semi-classical models, Inter-band transition. This work was partially supported by NSF grants DMS-4546 and NSF DMS NMS grant DMS-79.
2 L. CHAI, S. JIN AND Q. LI. Introduction. The linear equation with a periodic potential is an important model in solid state physics. It describes the motion of electrons in a crystal with a lattice structure. We consider the following one dimensional equation ( ( iεφ ε t = ε φε + V + U( φ ε, t >,. ( ε Here, φ ε is the comple-valued wave function, ε is the dimensionless rescaled Planck constant, V ( is a periodic potential with the lattice L = πz, so V ( + ν = V ( for any and ν L. U( is a smooth eternal potential function. We denote the dual lattice of L by L and L = Z. The fundamental domain of L is (, π and the first Brillouin zone is B = (,. In the semi-classical regime ε, the equation without a periodic potential has a semi-classical limit governed by the Liouville equation [35]. For the equation with a periodic potential, the Liouville equation can also be obtained along each Bloch band [], and it was justified rigorously in [3, 4, 5] for the case U. The result can be generalized to the case when a weak random potential [] and nonlinear self-consistent potential [3, 4] are presented. It has been shown that in the crystal, the electrons remain in a certain quantum subsystem, move along the m-th band and the dynamics is given by ẋ = k E m (k, where E m is the energy corresponding to the m-th Bloch band [5]. Higher order correction relevant to Berry phase can be included, see [8, 9, ]. All of these results use the adiabatic assumption, namely, different Bloch bands are well-separated and there is no band-crossing. The non-adiabatic or diabatic effect should be considered whenever the transitions between energy levels of the quantum system play an important role. This may happen when the gap between the energy levels becomes small enough in comparison to the scaled Planck constant ε. The well-known Landau-Zener formula [, 37] describes the asymptotic effect of avoided crossings in various specific situations. The study of such quantum tunnelings is important in many applications, from quantum dynamics in chemical reaction [34], semiconductors to Bose-Einstein condensation [6]. While in the case of band separation there have been significant mathematical progress in understanding the semi-classical limit [, 3, 5, 8, 9], as well as numerical methods that utilizes the Bloch decomposition [7], there has been little mathematical and computational works for the band-miing case. In the contet of surface hopping method, associated with the Born-Oppenheimer approimation, the Landau-Zener phenomenon has been studied computationally by Tully etc. [34, 33, 3, ] and mathematically [5,, 3, ]. For a Liouville equation based computational approach, see [9], and a quantum-classical model for surface hopping, see [6, 6]. In this paper, we use the Wigner-Bloch theory to derive semi-classical models for the equation with periodic potentials ( that account for bandcrossing. Without band-crossing, the classical Liouville equation can be obtained along each Bloch band, as ε. This only includes the diagonal entries of the Wigner-Bloch matri (the Wigner transform of the Bloch functions which is valid away from the crossing zone. In our semi-classical models we include the leading order of the off-diagonal entries, resulting in system of comple, inhomogeneous Liouville equations. These systems contain terms that describe for transitions between bands, as well as Berry phase information which is related to the quantum Hall effect [3]. These systems are still hyperbolic, with oscillatory forcing terms
3 SEMI-CLASSICAL MODELS FO BAND COSSINGS 3 that converge (in the weak sense to zero in the semi-classical limit away from band-crossing zones so the classical Liouville equation can be recovered. Several numerical eperiments using these models produce numerical solutions that adequately describe the quantum transitions between bands, when compared with the solutions of the original equation (. The computational cost of these semi-classical models, while considerably lower than that of the original equation, is higher than the adiabatic Liouville equations without band-crossing. In order to further reduce the computational cost, a hybrid method that couples the classical Liouville equation away from the crossing zone to the new semi-classical models to be used in the crossing zones is introduced. A couple of other more efficient computational approaches for linear periodic potentials, are also introduced in section 5. This paper is organized as follows. In Section, we introduce the Wigner transform and the Bloch eigenvalue problem that are the two tools to study the semiclassical limit of the equation with periodic potentials. We also review the classical limit without band-crossing which gives rise to the classical Liouville equation for each Bloch band. In Section 3 we derive the new semi-classical Liouville systems that account for quantum transitions between different bands. In Section 4 we introduce a domain-decomposition based hybrid model that couples the classical Liouville system away from the band crossing zone to the new semiclassical models used in the crossing zones. Some numerical eamples are presented in Section 5 to validate these semi-classical models for quantum transitions between Bloch bands. More efficient numerical methods were introduced for simpler linear potentials. Section 6 concludes this paper.. The classical limit without band-crossing... The Wigner transform. Define the asymmetric Wigner transformation as in [], dy W ε (t,, k = π eiky φ ε (t, εy φ ε (t,. ( where φ ε is the comple conjugate of φ ε, and φ ε is the solution to the equation (, then one obtains the Wigner equation: W ε t + k W ε + iε W ε = e iµ/ε ˆV (µ [Wε (, k µ W ε (, k] iε µ L + dω iε π eiω Û(ω [W ε (, k εω W ε (, k], where Û(ω is the Fourier transform of U(y: Û(ω = dy e iωy U(y, ω, and ˆV (µ is the discrete Fourier coefficients of V (y: ˆV (µ = π π dy e iµy V (y, µ L. Here (, π is the fundamental domain of the lattice L. (3
4 4 L. CHAI, S. JIN AND Q. LI Denote z = ε as the fast variable. To separate the dependence on both the slow and the fast variables, one can write W ε (t,, k as W ε (t,, z, k, and replace the spatial derivative by + ε z. Then (3 becomes: ( t W ε + k + W ε + iε ( ε z + W ε ε z = e iµ/ε ˆV (µ [Wε (, k µ W ε (, k] (4 iε µ L + dω iε π eiω Û(ω [W ε (, k εω W ε (, k]. Assume U( is smooth enough, one throws away high order terms, and (4 becomes: W ε t + k W ε U W ε k + W ε i z = ε LW ε, (5 where the skew symmetric operator L is given by Lf(z, k = k f z + i f z i µ L e iµ/ε ˆV (µ [f(, k µ f(, k]. Asymptotically one epands W ε by: W ε (t,, k = W ( t,, ε, k + εw ( t,, ε, k + and plugs this ansatz into (4. In the order of O( ε and O( respectively, one gets: W t + k W LW =, U W k + W i z = LW. (6a (6b This implies that W is in the kernel of L. We seek a good basis of kerl in section., and leave the eact formulation of W to section.3... The Bloch eigenvalue problem. The eigenfunctions of L are constructed by studying the following eigenvalue problem: Ψ(z, p + V (zψ(z, p z = E(pΨ(z, p, (7a Ψ(z + ν, p = e ipν Ψ(z, p, ν L, (7b Ψ (z + ν, p z = Ψ eipν (z, p, ν L. z (7c With each specific p, one constructs a boundary condition (7b and (7c and solves the eigenvalue problem. Denote E m (p as the m-th eigenvalue with multiplicity r m, and Ψ α m, with α =,..., r m, as the associated α-th eigenfunction. Ψ m is usually called the m-th Bloch eigenfunction [5]. For this problem, one can easily check the following properties: (a The eigenvalues E m (p are L periodic in p and have constant finite multiplicity outside a closed zero-measure subset F of p. Outside F, one orders the eigenvalues as E (p < E (p <... < E m (p <... with E m (p as m, uniformly in p.
5 SEMI-CLASSICAL MODELS FO BAND COSSINGS 5 (b For any p (,, {Ψ α m (, p} forms a complete orthonormal basis in L (, π, i.e. (Ψ α m, Ψ β n := π dz π Ψα m(z, p Ψ β n(z, p = δ mn δ αβ. (8 (c For all φ L (, one has the following Bloch decomposition: r m φ( = c α m(pψ α m(, p dp (9 Define m α= where c α m is the Bloch coefficient: c α m(p = φ( Ψ α m(, p d. B Φ α m(z, p = e ipz Ψ α m(z, p, then Φ α m is a z periodic function with period π. For any p (,, {Φ α m (, p} forms a complete orthonormal basis in L (, π. For simplicity, throughout this paper we make the following assumption: Assumption. The multiplicity for each eigenvalue E m is outside F. We will drop the superscript α when there is no ambiguity..3. The classical limit without band-crossing. In this section, we review the classical limit without band-crossing []. For that, we need the following assumption: Assumption. The set F is a empty set, i.e. eigenvalues are strictly apart from each other everywhere in p, namely E (p < E (p <... < E j (p <.... By taking the Wigner transformation on the Bloch eigenfunctions, one obtains a basis on the phase space. Define the z periodic functions Q mn (z, k by: Q mn (z, k = Q mn (z, µ k, p k = π where k is an arbitrary real number and is decomposed as: k = p k + µ k, p k B, µ k L. Lemma.. Define the inner product, : f, g := µ L π dy π eiky Ψ m (z y, p k Ψ n (z, p k, ( dz π f(z, µḡ(z, µ, f, g L( (, π, l (L, then for any p ( /, /, {Q mn (,, p} forms a complete orthonormal basis in L ( (, π, l (L. Proof. The orthonormal condition Q mn, Q jl = δ mj δ nl, ( can be proved by simply using (8. To prove the completeness, it is sufficient to show that: if there eists an f L ( (, π, l (L, such that f, Q mn = for all m, n N, then f(z, µ. Assume f, Q mn = for all m, n N. By the definition of Q mn, µ π π dydz (π f(z, µeiµy Φ m (z y, p Φ n (z, p =, m, n N.
6 6 L. CHAI, S. JIN AND Q. LI Since {Φ n (, p} forms a complete orthonormal basis in L (, π, the above equality implies that thus µ π which implies that f(z, µ. A straightforward computation gives: dy π f(z, µeiµy Φ m (z y, p =, m N, f(z, µe iµy, µ LQ mn (z, k = LQ mn (z, µ, p k = i(e m (p E n (pq mn (z, µ, p k. ( Apparently, under Assumptions and, ( gives: Therefore, from (6a W has the form kerl = span{q mm, m =,, }. (3 W (t,, z, k = m σ mm (t,, p k Q mm (z, µ k, p k, (4 with σ mm representing the epansion coefficients. To derive the equation for them, one plugs it back into (6b, and takes the inner product with Q mm on both sides, the right hand side vanishes due to the skew symmetry of L and (: The left hand side, on the other hand, gives: LW, Q mm = W, LQ mm =. (5 t W + k W U k W + i z W, Q mm = [ ] t σ nn Q nn, Q mm + σ nn kq nn, Q mm n [ U ( p σ nn Q nn, Q mm + σ m p Q nn, Q mm ] (6 n = t σ mm + p E m σ mm U p σ mm. Thus one can combine (5 and (6, and obtain: t σ mm + p E m σ mm U p σ mm =. (7 In the derivation of (6, the following equalities were used : kq mn, Q jl = i ( δnl ( z Ψ m, Ψ j + δ mj ( z Ψ l, Ψ n, (8a p Q mn, Q jl = δ nl ( p Φ m, Φ j + δ mj (Φ l, p Φ n, (8b z Q mn, Q jl = δ nl ( z Ψ m, Ψ j δ mj ( z Ψ l, Ψ n, (8c p E m δ mj + ( p Φ m, Φ j (E m E j = i ( z Ψ m, Ψ j, (8d ( p Φ m, Φ m + (Φ m, p Φ m =. (8e The details can also be found in []. Notice that (7 is the classical Liouville equation for each Bloch band. emark. Similar results were rigorously proved in [5, 3].
7 SEMI-CLASSICAL MODELS FO BAND COSSINGS 7 3. Asymptotic models for the band-to-band transition. This section is for the case when Assumption is not satisfied. Physically, ε is small but nonzero, and electrons can tunnel across bands for all p. But this tunneling is negligible when bands are far away from each other, the so-called adiabatic assumption. This is the basis for the asymptotically epansion (4, which throws away the dependence on the band-transition terms Q mn (m n of W. However, if there eists m n and a point p c, such that E m (p c E n (p c O(ε γ with γ being a real number larger than zero, then LQ mn O(ε γ, so the asymptotic epansion (6 does not hold any more. When this happens, physically one observes significant tunneling effect. We seek asymptotic models to handle the band-to-band transition in this section. 3.. The derivation of a two-band semi-classical Liouville system using the asymmetric Wigner transform. In order to handle the band-to-band transition phenomena, we come back to the asymptotic model (5 and use the following epression (compare with (4!: W ε = m σ mm Q mm + m n σ mn Q mn. (9 Without loss of generality, we tackle a two-band problem. Define the crossing point as p c = arg min p { E E } and assume and p c =. Plugging (9 into (5 and taking the inner product with Q mn as mentioned in Sec..3, one gets a system for σ mn s: σ t = U σ + E p [( Φ p, Φ + ( Ψ i z, Ψ ( σ + Φ, Φ p σ ] σ, U σ p σ t = U σ + E p [( Φ p, Φ + ( Ψ i z, Ψ ( σ + Φ, Φ p σ ] σ, U σ p σ t = U σ + E p [( Φ, Φ p + ( Ψ i z, Ψ σ + ( Φ p, Φ σ U σ + σ p + ie E ε σ + ( Φ, Φ p σ ( Φ p, Φ ] σ, σ t = U σ + E p [( Φ, Φ p + ( Ψ i z, Ψ σ + ( Φ p, Φ σ U σ + σ p + ie E ε σ + ( Φ, Φ p σ ( Φ p, Φ ] σ.
8 8 L. CHAI, S. JIN AND Q. LI This system can be written in vector form as: t σ + A σ + B p σ = BCσ + id ε σ (a where σ = ( σ σ σ σ T (b B = U I, D = diag (, E E, E E, p E ψ A = ψ p E p E ψ, ψ p E C = φ φ φ φ φ φ φ φ φ φ φ φ, (c (d (e ψ mn = i ( z Ψ m, Ψ n, and φ mn = ( p Φ m, Φ n, (f Where the superscript T denotes the matri transpose, I is the identity matri. Generally, ψ = ψ, φ = φ are comple-valued quantities, and φ, φ are purely imaginary, so A = A is Hermitian, and C = C is anti-hermitian (where the superscript denotes the matri conjugate transpose. Note that Uφ mm is the so-called Berry-connection [7, 9, 3]. emark. According to the adiabatic theory [36, 7], the state initially in one eigenstate Φ m (p( will remain in an instantaneous eigenstate Φ m (p(t when time t is sufficiently slow. One can write the state at time t as ( φ m (t = e iγm(t ep i ε t dt E m (p(t Φ m (p(t, m =,, ( where p(t is the trajectory satisfying ṗ(t = U. The second eponential in ( is known as the dynamical phase factor, and γ m in the first eponential is the path integral of the Berry connection, i.e. γ m = i φ mm (p dp, ( which is called the Berry phase. Obviously this term cancel out in the density function φ m φm, however if considering the inter-band density φ m φn with m n, then one obtains { ( φ m (p(t φ n (p(t = ep i γ m (t γ n (t ε Φ m (p(t Φ n (p(t, t dt ( E m (p(t E n (p(t } (3
9 SEMI-CLASSICAL MODELS FO BAND COSSINGS 9 and the phase in the above equation i d ( γ m (t γ n (t t dt ( E m (p(t E n (p(t dt ε = dp ( φmm (p(t φ nn (p(t i Em (p(t E n (p(t dt ε( = U (φ mm φ nn i ( ε (E m E n = BC + id. ε In comparison with (, one can see that it is just the coefficient of the σ mn term on the right hand side of (a for σ mn. This shows that the Berry phase plays an important role in the inter-band transitions. It can be shown the system ( is a hyperbolic system (see the details in Appendi A. Since the source matri C = C, all of its eigenvalues are purely imaginary. Thus this system has unique bounded solution if the initial data are bounded [9] uniformly in ε. For later convenience, we call the semi-classical Liouville system ( system, and this is the system that will be discussed and numerically solved in the paper. One also needs to equip it with appropriate initial condition. Choose the initial data of the equation as two wave packets along the two Bloch bands in the following form [7, 7]: mn (4 ( ( φ I = a ( Φ ε, S ( e is(/ε + a ( Φ ε, S ( e is(/ε. (5 Then, the initial data of the Wigner function, for ε, has the approimation: W I (, z, k a ( W (z, k + a ( W (z, k + a (a ( ( W (z, k + W (z, k, (6 with W mn (z, k = dy π eiky Φ m (z y, S ( εy Φ n (z, S (e i(s( εy S(/ε. Using Taylor epansion on S ( εy S ( and Φ m (z y, S ( εy, one gets dy W mn (z, k = π ei(k S(y Φ m (z y, S ( Φ n (z, S ( + O(ε, then by ignoring the O(ε term, and using the periodicity of Φ m (z, p on z, one can change the integral into a summation of integrals from to π: W mn (z, k = ν L π dy π ei(k S((y+ν Φ m (z y, S ( Φ n (z, S (. Applying the equality e ikν = ν L µ L δ(k + µ,
10 L. CHAI, S. JIN AND Q. LI one gets W mn (z, k = µ L δ(k + µ S = µ L δ(p k + µ S = µ L δ(p k + µ S π π π π = δ(p k + µ S µ L = ( pk + µ S Qmn (z, µ k, p k. µ L δ dy π ei(k Sy Φ m (z y, S Φ n (z, S dy π eiky Ψ m (z y, S Ψ n (z, S dy π eiky Ψ m (z y, p k + µ Ψ n (z, p k + µ dy π eiky Ψ m (z y, p k Ψ n (z, p k In the above derivation, from the second line to the third line, we use the fact that δ(p p f(p g(p dp = f(p g(p = δ(p p f(pg(p dp to replace the argument S to p k + µ. Without loss of generality, we assume that S ( /, /, then (7 becomes W mn (z, k = δ ( p k S Qmn (z, µ k, p k. (8 Compare (8 with (6, one has the initial data for σ: (7 σ(,, p = δ (p S ( ( a a a a a a T. (9 3.. The semi-classical Liouville system using the symmetric Wigner transform. All the analysis above is done in the framework of the asymmetric Wigner transform (. One could also use the symmetric Wigner transform: Wε s (t,, k = dy π eiky φ ε ( t, εy φε ( t, + εy The derivation is similar, thus we skip the details, and give a list of the results:. The Wigner equation corresponding (3 is Wε s + k W ε s t = [ e iµ/ε ˆV (µ Wε s iε µ L + dω iε π eiω Û(ω (, k µ. ( Wε s, k + µ ] [ ( Wε s, k εω ( Wε s, k + εω ].. Corresponding to the asymptotic Wigner equation for asymmetrical transformation (5, one has: Wε s + k W ε s t U Wε s k = ε Ls Wε s, (3 where the skew symmetric operator L s is given by L s f(z, k = k f z e iµ/ε ˆV (µ [f(, k µ/ f(, k + µ/]. i µ L
11 SEMI-CLASSICAL MODELS FO BAND COSSINGS 3. Corresponding to (6, one has the O ( ε and O( epansions: W s t + k W s L s W s =, (3a U W s k = Ls W s. (3b 4. Same as in (, one has the following symmetrical definition for Q mm Q s mn(z, k = Q s mn(z, µ, p = π They are eigenfunctions of L s dy ( π ei(p+µy Ψ m z y (, p Ψn z + y, p. L s Q s mn = i(e m E n Q s mn. 5. If the eigenvalues {E n } are well separated, i.e. E m E n, for m n, the solution to (3a is: W s = m σ s mmq s mm. By taking the inner product with Q s mm on both sides of (3b, one obtains the same classical Liouville equations for σ s mm as in (7 t σ s mm + p E m σ s mm U p σ s mm =. 6. If some bands touch at point p c, the solution to (3 is given by: W s ε = m σ mm Q s mm + m n σ mn Q s mn. In the two-band case, m, n =,, then σ s mn is governed by, t σ s + A s σ s + B s p σ s = B s C s σ s + ids ε σs, (3 where σ s = ( σ s σ s σ s σ s T, then B s = B, C s = C and D s = D are the same as the asymmetric case, while A s is given by p E ψ ψ A s = ψ E +E p ψ ψ E +E p ψ ψ ψ p E and ψ mn, φ mn are given by (f. Noted that A s = (A s. We call this new system obtained by the symmetric Wigner transform (3 the Liouville-S system. Apparently the only difference from lies in the transport matrices A A s. But they share the same weak limit as ε (see Appendi B for detail. This formally suggests the similar behavior of σ, σ and σ, s σ, s which is confirmed by numerical results in Section A multiscale domain decomposition method.
12 L. CHAI, S. JIN AND Q. LI 4.. Asymptotic behavior of the two-band semi-classical Liouville systems. Away from p c, E and E are well-separated, and the i ε terms for the transition coefficients σ and σ in /S lead to high oscillations, thus as ε, the system formally goes to its weak limit: the classical one (7. For ε >, around p c, however, both σ and σ become significant, and the band-to-band transition is no longer negligible. In fact, based on the distance to the crossing point p c, one could obtain some asymptotic properties of the transition coefficients σ mn (m n. Assume that the initial data for the transition coefficients are all zero, and that U does not change sign in time for all, take U > for eample, then: Case : If p C ε for C = O(, then σ and σ are of o( ε; Case : If C ε < p < C ε, then σ and σ are of O( ε, and σ and σ are slowly varying, i.e. t σ O ( ε, and t σ O ( ε ; Case 3: If p C ε, σ and σ are highly oscillatory with mean. We leave the justification for a simpler model problem to Appendi B. emark 3. The assumption σ mn (t = = (m n is a reasonable assumption. In fact, given arbitrary initial condition, one could check that away from p c, the weak limit of σ mn is always zero, as ε, as can be seen in Appendi B. So numerically we treat the initial data for both σ and σ zero, given the initial velocity S ( away from the crossing point, i.e. σ(,, p δ (p S ( ( a a T, if S ( p c. (33 This assumption is intuitive and empirical, but it does give us some convenience in solving the /S numerically. In fact, the numerical eamples provided in Section 5 indeed show that the band-to-band transition is captured very well with initial data ( A domain decomposition method. Clearly, one has several observations in hand: the classical Liouville is an approimation (in weak limit to away from the crossing point; σ and σ are slowly varying in a neighborhood of p c =, with a small amplitude before the characteristic hitting p c, and a rapid oscillation after that. These observations motivate a domain decomposition method in p-space. The idea is, away from p c, when the classical Liouville equation (7 is a good approimation, we solve this set of equations, and then switch back to when the approimation breaks down around p c (in the O(ε neighborhood of p c. The gain is obvious: numerically it is much easier and more efficient to solve the classical one, thus this approach saves a great amount of computational cost than solving everywhere. Based on the asymptotic properties of the transition coefficients, we propose the following: Given a fied spatial point, the sign of U determines the traveling direction of wave in p. Assume U >, i.e. the wave we study is right-going: Classical regions: p < C ε and p > C ε: In this region, coarse mesh independent on ε is used to solve the classical Liouville system (7, and σ and σ are set to be zero.
13 SEMI-CLASSICAL MODELS FO BAND COSSINGS 3 Semi-classical region: p B \ {Classical region}: Solve. The incoming boundary conditions for σ and σ are set to be zero, and the incoming boundary condition for σ and σ are the inflow boundary condition. A fine mesh is used with and t much less than ε. In case of U <, and the wave if left-going, boundary condition can be set up in the same way. emark 4. Our analysis is based on the regularity of the coefficient matri C in. But usually, the value of C s element can be of O(/δ where δ = min p E (p E (p is the minimal band gap. So C will be large if the minimal band gap δ is small, and the numerical discretization in the semi-classical region should resolve this small parameter δ. In the interested regime ε δ, where the transition cannot be ignored for avoided crossing when E (p c E (p c O( ε, o( ε mesh is enough. 5. Numerical eamples. In this section we solve the and Liouville-S numerically. Firstly in Section 5., we present a numerical method for the equation with a linear eternal potential U(. In this special case, we provide an efficient solver without using the domain decomposition method. In Section 5., the domain decomposition method is applied for general U(. For both eamples, we use the Mathieu model, i.e. the periodic potential is V (z = cos z. The first eight Bloch eigenvalues are shown in Figure. Apparently, th band 5th band E p (a E p (b Figure. The eigenvalues of the Mathieu Model, V ( = cos. some eigenvalues get very close to each other around p =, ±.5. We will focus on the 4th and 5th bands. Denote Ψ and Ψ as the Bloch functions corresponding to the 4th and 5th bands respectively. For comparison, we will compare the numerical results to the ones given by the original equation, computed through the methods given in [7] with mesh size and time step much smaller than ε. The minimum gap between the 4th and 5th bands is.47, located at p =. The gap is considered small enough so that the quantum effect can be seen for the ε being used.
14 4 L. CHAI, S. JIN AND Q. LI 5.. Linear U(. We deal with the linear eternal potential in this section: U( = U β. Then the equation is: [ ( iεφ ε t = ε ] φε + cos + (U β φ ε, (34 ε with the initial data given as a wave packet along the 4th band: ( φ I = a ( Φ ε, S ( e is(/ε, with S ( = p. (35 Correspondingly, the becomes: with given by: t σ + β p σ = σ, σ(,, p = σ I (, p (36a = a ( δ (p S ( ( T, (36b = βc + id ε A. (37 One encounters two computational challenges here. Firstly, one needs to numerically resolve the rapid oscillation. We will present an efficient way to overcome this difficulty by following the characteristics in Section 5... Secondly, the initial data contains a delta function. Usually one uses a Gaussian function with small variance to approimate it, and the error is related to this variance. As will be discussed in Section 5.., in some special cases, this can be avoided by using the singularity decomposition idea of [8]. emark 5. Here we only discuss. Liouville-S can be computed similarly A Fourier transform based integration method. Let the Fourier transform of f(t,, p with respect to be ˆf(t, η, p = e iη f(t,, pd. Taking this transform on (36a, one gets: t ˆσ + β p ˆσ = ( βc + idε iηa ˆσ =: ˆˆσ In this special case, β is a constant so the characteristic line can be obtained analytically. and consequently one can avoid using the mesh thus removes the difficulty due to the high oscillation introduced by the term id ε, as will be clear in the following. We take the first time step t [, t] for eample. Along the characteristic line p(t = p + βt, one evaluate ˆσ and ˆ at (t, p(t, η and has: dˆσ dt = ˆ ˆσ. Solution to this ODE system satisfies ˆσ (t =ˆσ ( + t t iη p E ˆσ (tdt ( β(φˆσ φ ˆσ iη ψ ˆσ dt, (38a
15 SEMI-CLASSICAL MODELS FO BAND COSSINGS 5 and where: and ˆσ (t =e t Kε (τdτ ˆσ ( + t e t s Kε (τdτ ( G(sˆσ (s H(sˆσ (s ds, K ε = i ε (E E iη p E β(φ φ, For t small, an approimation to(38 is G = β φ iηψ, H = β φ. ˆσ (t ˆσ ( iηt p E ( p(t ˆσ (t + ( βφ ( iηψ ( t ˆσ dt βφ ( ˆσ (t ˆσ (e t Kε (τdτ + ( G(ˆσ ( H(ˆσ ( t t e t s Kε (τdτ ds. ˆσ dt, Plug (39b into (39a, to evaluate ˆσ ( t and ˆσ ( t, one needs to compute: F := e t F := F := Kε (τdτ, t t t e t Kε (τdτ dt, e t s Kε (τdτ ds dt. (38b (39a (39b However, to compute F, F and F is not easy since their integrands are highly oscillatory. But if one chooses β t = p, then at each time step one follows eactly the characteristics, so p(t always lie on the grid points, thus F, F and F are time independent, and one only needs to compute them once (with a highly resolved calculation. Similar analysis can be carried out for ˆσ and ˆσ. We prove the stability of this method by using it on a simpler model problem, and it will be justified in Appendi C A singularity decomposition idea. To handle the delta function in the initial condition (36b, one usually approimates it with a Gaussian function with small variance, and numerical error was determined by the width of the Gaussian. As stated before, in some special cases, this error could be avoided, and the eample we are discussing here is for when S ( p = const, for which, we apply the singularity decomposition method introduced in [8] to reduce the error. Write the ansatz of ˆσ mn (t,, p as: σ(t,, p = ω(t,, p δ ( θ(t, p (4 in which: θ(t, p = p (p + βt, which solves the Liouville equation ω satisfies the same equation as σ: t θ + β p θ =, t ω + β p ω = ω. (4
16 6 L. CHAI, S. JIN AND Q. LI These can be proved by simple derivations. Formally, one has σ t = ( ω ω δ(θ = t t δ(θ + ω δ (θ θ t = ω δ(θ β ω p δ(θ βω δ (θ = ω δ(θ β ω p δ(θ βω δ (θ θ p = ω δ(θ β ( ω δ(θ p = σ β σ p. The equalities above should be understood in the distributional sense. The decomposition (4 enables one to solve for ω and θ separately with good (bounded initial data a ( and S ( respectively. The equation ω satisfies is the same as the one for σ, thus the numerical method introduced in Section 5.. can be used. In the final output, one needs to get back to σ using (4, so a discrete delta approimation is only needed at the output time, not during time evolution Numerical eperiments. We show the numerical results of the /S with the following data β =, p 5( π =.5, a ( = ep (. (4 We compute the density, the cumulative density function (c.d.f. and mass in the st band, defined respectively by: ρ ε = φ ε, γ ε = ρ ε (ydy, m (t = P φ ε (t, d, (43 where P n is the projection onto the nth band: P n φ( = dp dyφ(y Ψ n (y, pψ n (, p, φ L (, m N. B The two integrals in (43 are calculated by the midpoint quadrature rule numerically. Figure shows the density and c.d.f. computed for the equation, the and the Liouville-S respectively at t =.5. The results match quite well. Figure 3 shows the evolution of m as a function of time t. One can see the total mass on the first band jumps down at around t =.5, when the momentum p reaches p c =, reflecting the 4th-to-5th band transition. The eperiment also shows that smaller ε gives smaller transition rate. Note that some small oscillations occur around the crossing region. They are related to the interference phenomena, and are usually called the Stueckelberg oscillation [8, 3, 33]. Define L error in the cumulative distribution function (c.d.f. [4, ]: Err ε ( (t = ρ ε S (t, z ρ ε L(t, z dz d, (44 where ρ ε S and ρε L denote the density calculated by the equation and the Liouville system respectively. Numerically we compute (44 using the midpoint
17 SEMI-CLASSICAL MODELS FO BAND COSSINGS Liouville-S.8 Liouville-S ρ ε γ ε (a ε = Liouville-S.8 Liouville-S ρ ε.5.5 γ ε (b ε = Liouville-S.8 Liouville-S ρ ε.5.5 γ ε (c ε = Figure. The eample in Sec. 5..3, t =.5. The left and right column are for the position density ρ ε, and c.d.f. γ ε respectively. The solid line, the dash line and the dotted line are the numerical solutions to the, the and the Liouville-S respectively.
18 8 L. CHAI, S. JIN AND Q. LI.8 Liouville-S.8 Liouville-S.6.6 m.4 m t (a ε = t (b ε = 9.8 Liouville-S.9.8 Liouville-S m.6.4 m t (c ε = t (d ε = Figure 3. The eample in Sec : time evolution of m (t defined in (43. quadrature rule. Figure 4 shows this at time t =.5. As ε, the Liouville system gets more accurate, and the error decreases with the speed of O(ε. 5.. A domain decomposition computation. This section shows eamples with varying U. For this general case, the p-characteristic is no longer a straight line, and the fast solver in the previous section is no longer valid. To numerically solve, we use the domain decomposition method. The classical finite volume method is used for the convection terms. We compute the system with both a pure and a mied state initial data with: U( = sin A pure state initial data. In this eample, we use the same pure state initial data as in the previous eample (35,(4. Correspondingly, the initial data for the system is given by (36b. Numerically a Gaussian function centered at p with variance of ε/6 is used to approimate the δ function. Figures 5 and 6 show the density, the c.d.f at t =.5 and evolution of m (43 computed for both the equation and the system. The numerical results for the two systems agree well. Figure 7 gives decay of Err ε (44.
19 SEMI-CLASSICAL MODELS FO BAND COSSINGS 9.5. Liouville-S.5 Err ε ε Figure 4. The eample in Sec : Err ε as function of ε at t =.5. The numerical results show that the hybrid model can capture the band-to-band transition phenomena, and the error decays like O(ε. emark 6. As U varies with, the wave packet becomes decoherent. This weakens the interference phenomenon [8, 3, 33]. As one can see in Figure 6, the Stueckelberg oscillations around the crossing region is much weaker than those in the previous eample A mied state initial data. This eample is for the case when the initial data is a mied state: [ ( ( ] = a ( Φ ε, p e ip /ε + Φ ε, p e ip /ε, p =.5. φ I Correspondingly, the initial data for the semi-classical Liouville system should be: σ = a ( δ(p p [,,, ] T. Since p is away from the crossing point and σ and σ weakly converge to zero as ε, numerically, we regard them as zero and use σ a ( δ G (p p [,,, ] T as the initial condition, where δ G (p is a Gaussian function centered at zero. The density and the c.d.f. are computed for the and the, compared in Figure 8. Err ε as a function of ε is shown in Figure Conclusion. In this paper we derive semi-classical models for the linear equation with periodic potentials. These models take into account the band crossing, which is important to describe quantum transitions between different Bloch bands. Away from the band-crossing zones these models reduce (in the sense of weak limit to the classical Liouville system for each Bloch band. We also couple these semi-classical models (to be used near the crossing zones and the classical
20 L. CHAI, S. JIN AND Q. LI ρ ε γ ε (a ε = ρ ε γ ε (b ε = ρ ε.5 γ ε (c ε = Figure 5. The eample in Sec t =.5. The left and right columns show the position density ρ ε, and the c.d.f. γ ε respectively. Liouville equation (used away from the crossing zones for an efficient multiscale computation. Our numerical eperiments show that these semi-classical models provide correct quantum transitions near the crossing zones when compared with the direct simulation of the original equation.
21 m.6 m t (a ε = t (b ε = m.8.6 m t (c ε = t (d ε = Figure 6. The eample in Sec. defined in ( : time evolution of m (t Appendices A. Hyperbolicity of the two-band semi-classical Liouville systems. The semi-classical Liouville system ( is where t σ + A σ U p σ = Sσ (45 S = C + id ε, and σ, A, C, D are defined in (. Noted that A = A and S = S. To check the hyperbolicity, we separate the real and imaginary parts. σ = e σ + i Im σ, A = e A + i Im A, S = e S + i Im S, with e F and Im F denoting the real part and imaginary part of F respectively.
22 .6.4. Err ε ε Figure 7. The numerical eample in Sec. 5..: the L error Err ε between and the solution at t =.5. Define σ I = ( e σ, Im σ T, and [ ] [ e A Im A e S Im S A I =, and S Im A e A I = Im S e S ], then one can get t σ I + A I σ I U p σ I = S I σ I. (46 Since A = A, then e A = e A T, Im A = Im A T. Thus have that the matri A I = A T I is a symmetric matri, which implies the hyperbolicity of the system. In addition, since S = S, then e S = e S T, Im S = Im S T. Therefore the matri S I = SI T is a skew-symmetric matri. Similarly, one can obtain the hyperbolicity of system Liouvill-S (3. B. Some basic analysis of the semi-classical Liouville systems. To understand the asymptotic behavior of the solution to the system (, as mentioned in Section 4., we look at a simpler model system: t g + f + b( p g =, t f + a(p f + g + b( p f = i ε c(pf, (47 g(,, p = g I (, p, f(,, p = f I (, p. The initial conditions g I and f I are bounded smooth functions independent on ε, b > and the set of zeros for c(p: S c = {p : c(p = } is measured zero. It is easy to check that (47 is a linear hyperbolic system, and the solutions g and f are bounded uniformly in ε [9].
23 ρ ε γ ε (a ε = ρ ε γ ε (b ε = ρ ε γ ε (c ε = Figure 8. The eample in Sec. 5..: at time t =.5. The left and the right columns show the position density ρ ε, and the c.d.f. γ ε respectively. B.. Weak convergence. We consider the weak limit of the solution of (47 in this subsection. To do this, we introduce the inner product, as u, v = u(t,, p v(t,, p d dp dt.
24 Err ε ε Figure 9. The numerical eample in Sec. 5..: the L error Err ε as a function of ε at t =.5. Choose an arbitrary test function h C ( +, take the inner product on both side of (47 w.r.t h, one gets { t g, h f, h + b p g, h =, f, t h + af, h + g, h + bf, p h = i ε cf, h. (48 The derivatives in the equation of (48 are acted on the smooth function h, and the left side is bounded. One gets cf, h as ε for all h C ( +. Given that c is almost surely nonzero, and f is bounded, one gets f weakly. Combined with the first equation in (48, one gets t g + b( p g weakly. B.. Strong convergence: for constant b. In these two subsections, we formally prove that before getting close to the crossing region, c(p is assumed to be bigger than a constant c that is unrelated to ε. In this region, f is constantly small and controlled by O(ε. This subsection is for the case when the speed on p direction is a constant: b( = β. Along the p-characteristic line p(t = p + βt, one applies the Fourier transform to the -variable, and gets: where f(t, η = ( ĝ(t, η, p(t, ˆf(t, η, p(t T and ( η (t = η c(p(t/ε ηa(p(t d f = i(tf, (49 dt.
25 5 The two eigenvalues of (t are both real, and thus the system above has a bounded solution satisfying The equivalence between norms gives: ĝ(t, η + ˆf(t, η = ĝ I (η + ˆf I (η. ĝ(t, η + ˆf(t, η < C( ĝ I (η + ˆf I (η. Adopt it into the solution to (49, one gets t g(t, = π η ˆf(t, ηe iη dη η π ˆf(t, η dη ( C η ĝ I (η + ˆf I (η dη t g(t, = π η ˆf(t, ηe iη dη η π ˆf(t, η dη ( C η ĝ I (η + ˆf I (η dη If the initial conditions ĝ I and ˆf I are smooth enough, and decay fast as η, one could easily get that t g(t, and tg(t, are both bounded in time independent of ε, and thus g and g are slowly varying in time. emark 7. In the derivation above, we dropped the p(t-dependence in the functions for simplicity. The partial derivative in the t-variable t should be understood as taking along the p-characteristic line p(t, i.e. t g(t, = t g(t,, p + β p g(t,, p. Assume that f I (, one follows the characteristics of (t by solving ẋ(t = a(t, and gets: f ( t, (t, p(t t ( i t g ( = ep c (p(τ dτ t s, (t s, p(t s ds. ε t s By the assumption, before hitting the crossing region, c(p(τ > c >, then the stationary phase argument suggests that, given slowly varying g(t,, p(t, f = O(ε. The observations from the above two subsections suggest that f, and before getting close to the crossing region, f is as small as of O(ε. Based on these arguments, for the system (, and propose the following conjecture: if σ and σ are initially zero, then: Case : If p ε, then σ and σ are of o( ε; Case : If p [ ε, ε], then σ and σ are of O( ε, and slowly varying; Case 3: If p ε, σ and σ are highly oscillatory, and converge to weakly. C. The integration method of a simple model system. In this section, we apply the method in (39 onto a simple model to show stability. d dt f(t, p + t = (p + t f(t, p + t, where f(t, p = (g(t, p, f(t, p T, p = p + t, and ( r (p r (p = (p r (p r(p ε.
26 6 with r(p ε = r (p + i ε c(p, while r and r are purely imaginary, and c(p real and positive. r, r and r are independent on ε. Set up mesh as t j = j t, and p i = + (i p, with p and t being the mesh size. Denote g j i and f j i as the numerical result at (t j, p i, then (39 gives f(t, p(t = f j i e t rε (pi +τdτ g j i r (p i t g(t, p(t = g j i + g(t, p(t r (p(t + r (p i e t s rε (pi +τdτ ds, t (5a f(t j + s, p i + s ds (5b Plug (5a into (5b, and evaluate them at (t j+, p i+, one obtains: ( r (p i+ t ( t t g j+ i+ = g j i r (p i e t s rε (pi +τdτ ds dt Written in vector form gives with M i = +f j i r (p i t ( r(p i t t e ts r ε (pi +τdτ ds dt r (p i+ t e t rε (pi +τdτ dt. f j+ i+ = M i f j i, (5 r (p i t e t r ε (pi +τdτ dt r (p i+ t r (p i t e t s r ε (pi +τdτ ds e t r ε (pi +τdτ The following quantities in the matri M i should be evaluated very accurately: F = e t F = F = r ε (pi +τdτ, t t t e t rε (pi +τdτ dt, e t s rε (pi +τdτ ds dt. emark 8. These three quantities only depend on the mesh grid point inde i but not the time steps inde j, thus they only need to be computed once at the beginning of the computation. Given ε t, the integrands of F and F are highly oscillatory, and one can see that F O(ε and F O(ε. Simple calculation shows that M i can be written as M i = Ω M i, with ( ( Ω = diag r (p j+ t, r (p, Mi = i F r F. r F F With purely imaginary r and r ε, it is easy to prove that Ω, and M i ( + O(ε t, and thus M i ( + O(ε t. This implies asymptotic stability of the scheme (5 independent of ε..
27 7 EFEENCES [] N.W. Ashcroft and N.D. Mermin. Solid State Physics, Saunders College, 976. [] G. Bal, A. Fannjiang, G. Papanicolaou and L. yzhik, adiative transport in a periodic structure, Journal of Statistical Physics, 95 (999, [3] P. Bechouche. Semi-classical limits in a crystal with a coulombian self-consistent potential: Effective mass theorems, Asymptotic Analysis, 9 (999, [4] P. Bechouche, N.J. Mauser and F. Poupaud. Semiclassical limit for the -poisson equation in a crystal, Communications on Pure and Applied Mathematics, 54 (, [5] F. Bloch, Über die quantenmechanik der elektronen in kristallgittern, Zeitschrift für Physik A Hadrons and Nuclei, 5 (99, [6] A. Bohm, P.A. Mostafazadeh, H. Koizumi, Q. Niu and J. Zwanziger, The Geometric Phase in Quantum Systems: Foundations, Mathematical Concepts and Applications in Molecular and Condensed Matter Physics, Springer,. [7]. Carles, P.A. Markowich and C. Sparber, Semiclassical asymptotics for weakly nonlinear Bloch waves, Journal of Statistical Physics, 7 (4, [8] M.S. Child, Atom-Molecule Collision Theory, Plenum, New York, 979. [9] ichard Courant and David Hilbert, Methods of Mathematical Physics, Differential Equations, Wiley-VCH, 8. [] K. Drukker, Basics of surface hopping in mied quantum/classical simulations, Journal of Computational Physics, 53 (999, 5 7. [] W. E, J. Lu and X. Yang, Asymptotic analysis of quantum dynamics in crystals: The Bloch- Wigner transform, Bloch dynamics and berry phase, Acta Mathematicae Applicatae Sinica, English Series, (,. [] Clotilde Fermanian Kammerer and Caroline Lasser, Wigner measures and codimension two crossings, J. Math. Phys., 44 (3, 57 57, 3. [3] P. Gérard, P.A. Markowich, N.J. Mauser and F. Poupaud, Homogenization limits and Wigner transforms, Communications on Pure and Applied Mathematics, 5 (997, [4] L. Gosse and N.J. Mauser. Multiphase semiclassical approimation of an electron in a onedimensional crystalline lattice-iii. From ab initio models to WKB for -Poisson, Journal of Computational Physics, (6, [5] George A. Hagedorn, Molecular Propagation through Electron Energy Level Crossings, Mem. Amer. Math. Soc., 994. [6] I. Horenko, C. Salzmann, B. Schmidt and C. Schütte, Quantum-classical Liouville approach to molecular dynamics: Surface hopping Gaussian phase-space packets, The Journal of Chemical Physics, 7 (, [7] Z. Huang, S. Jin, P.A. Markowich and C. Sparber, A Bloch decomposition-based split-step pseudospectral method for quantum dynamics with periodic potentials, SIAM Journal on Scientific Computing, 9 (8, [8] S. Jin, H. Liu, S. Osher and Y.H.. Tsai, Computing multivalued physical observables for the semiclassical limit of the equation, Journal of Computational Physics, 5 (5, 4. [9] S. Jin, P. Qi and Z. Zhang, An Eulerian surface hopping method for the equation with conical crossings, SIAM Multiscale Modeling and Simulation, 9 (, [] S. Jin and D. Yin, Computational high frequency waves through curved interfaces via the Liouville equation and geometric theory of diffraction, Journal of Computational Physics, 7 (8, , 8. [] L. Landau. Zur Theorie der Energieubertragung II. Physik, Physics of the Soviet Union, (93, [] C. Lasser, T. Swart and S. Teufel, Construction and validation of a rigorous surface hopping algorithm for conical crossings. Communications in Mathematical Sciences, 5 (7, [3] C. Lasser and S. Teufel, Propagation through conical crossings: an asymptotic semigroup, Comm. Pure Appl. Math., 58 (5, [4] P.L. Lions and T. Paul, Sur les mesures de Wigner, ev. Mat. Iberoamericana, 9 (993,
28 8 [5] P.A. Markowich, N.J. Mauser and F. Poupaud, A Wigner-function approach to (semi classical limits: Electrons in a periodic potential, Journal of Mathematical Physics, 35 (994, [6] C.C. Martens and J.Y. Fang, Semiclassical-limit molecular dynamics on multiple electronic surfaces, Journal of Chemical Physics, 6 (997, [7] O. Morandi and F. Schürrer. Wigner model for quantum transport in graphene, Journal of Physics A: Mathematical and Theoretical, (44, [8] G. Panati, H. Spohn and S. Teufel, Effective dynamics for Bloch electrons: Peierls substitution and beyond, Communications in Mathematical Physics, 4 (3, [9] G. Panati, H. Spohn and S. Teufel, Motion of electrons in adiabatically perturbed periodic structures, in Analysis, Modeling and Simulation of Multiscale Problems (eds. A. Mielke, Springer Berlin Heidelberg, (6, [3] D.S. Sholl and J.C. Tully. A generalized surface hopping method, The Journal of Chemical Physics, 9 (998, :77 77, 998. [3] M. Sillanpää, T. Lehtinen, A. Paila, Y. Makhlin and P.J. Hakonen, Landau Zener interferometry in a Cooper-pair bo, Journal of Low Temperature Physics, 46 (7, [3] G. Sundaram and Q. Niu, Wave-packet dynamics in slowly perturbed crystals: Gradient corrections and Berry-phase effects, Phys. ev. B, ( [33] J.C. Tully, Molecular dynamics with electronic transitions, The Journal of Chemical Physics, 93 (99, 6 7. [34] J.C. Tully and.k. Preston, Trajectory surface hopping approach to nonadiabatic molecular collisions: the reaction of H with D, The Journal of Chemical Physics, 55 (97, [35] E. Wigner, On the quantum correction for thermodynamic equilibrium, Physical eview, 4 (93, [36] D. Xiao, M.C. Chang, and Q. Niu. Berry phase effects on electronic properties, eviews of Modern Physics, 8 (, [37] C. Zener, Non-adiabatic crossing of energy levels, Proceedings of the oyal Society of London. Series A, Containing Papers of a Mathematical and Physical Character, 37(833:696 7, 93. eceived ; revised. address: chailh8@mails.tsinghua.edu.cn address: jin@math.wisc.edu address: qinli@math.wisc.edu
Semiclassical Models for the Schrödinger Equation with Periodic Potentials and Band Crossings
Semiclassical Models for the Equation with Periodic Potentials and Band Crossings Lihui Chai, Shi Jin, Qin Li Contents Introduction 3 2 The classical limit without band-crossing 5 2. The Wigner transform...............................
More informationSemiclassical computational methods for quantum dynamics with bandcrossings. Shi Jin University of Wisconsin-Madison
Semiclassical computational methods for quantum dynamics with bandcrossings and uncertainty Shi Jin University of Wisconsin-Madison collaborators Nicolas Crouseilles, Rennes Mohammed Lemou, Rennes Liu
More informationModeling and Computational Methods for Multi-scale Quantum Dynamics and Kinetic Equations
Modeling and Computational Methods for Multi-scale Quantum Dynamics and Kinetic Equations By Qin Li A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy
More informationASYMPTOTIC ANALYSIS OF THE QUANTUM DYNAMICS IN CRYSTALS: THE BLOCH-WIGNER TRANSFORM AND BLOCH DYNAMICS. 1. Introduction
ASYMPTOTIC ANALYSIS OF THE QUANTUM DYNAMICS IN CRYSTALS: THE BLOCH-WIGNER TRANSFORM AND BLOCH DYNAMICS WEINAN E, JIANFENG LU, AND XU YANG Abstract. We study the semi-classical limit of the Schrödinger
More informationEulerian semiclassical computational methods in quantum dynamics. Shi Jin Department of Mathematics University of Wisconsin-Madison
Eulerian semiclassical computational methods in quantum dynamics Shi Jin Department of Mathematics University of Wisconsin-Madison outline Classical mechanics: Hamiltonian system, discontinuous Hamiltonian,
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 informationAsymptotic Analysis of the Quantum Dynamics in Crystals: The Bloch-Wigner Transform, Bloch Dynamics and Berry Phase
Asymptotic Analysis of the Quantum Dynamics in Crystals: The Bloch-Wigner Transform, Bloch Dynamics and Berry Phase Weinan E, Jianfeng Lu and Xu Yang Department of Mathematics, Princeton University, Princeton,
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationA Hybrid Schrödinger/Gaussian Beam Solver for Quantum Barriers and Surface Hopping
A Hybrid Schrödinger/Gaussian Beam Solver for Quantum Barriers and Surface Hopping Shi Jin Peng Qi May 13, 2011 Abstract In this paper, we propose a hybrid method coupling a Schrödinger solver and a Gaussian
More informationAn Eulerian surface hopping method for the Schrödinger equation with conical crossings
An Eulerian surface hopping method for the Schrödinger equation with conical crossings Shi Jin Peng Qi Zhiwen Zhang October 6, 9 Abstract In a nucleonic propagation through conical crossings of electronic
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 informationGaussian beams for high frequency waves: some recent developments
Gaussian beams for high frequency waves: some recent developments Jianliang Qian Department of Mathematics, Michigan State University, Michigan International Symposium on Geophysical Imaging with Localized
More informationThe Vlasov-Poisson Equations as the Semiclassical Limit of the Schrödinger-Poisson Equations: A Numerical Study
The Vlasov-Poisson quations as the Semiclassical Limit of the Schrödinger-Poisson quations: A Numerical Study Shi Jin, Xiaomei Liao and Xu Yang August 3, 7 Abstract In this paper, we numerically study
More informationGAUSSIAN BEAM METHODS FOR THE SCHRÖDINGER EQUATION WITH DISCONTINUOUS POTENTIALS
GAUSSIAN BEAM METHODS FOR THE SCHRÖDINGER EQUATION WITH DISCONTINUOUS POTENTIALS SHI JIN, DONGMING WEI, AND DONGSHENG YIN Abstract. We propose Eulerian and Lagrangian Gaussian beam methods for the Schrödinger
More informationGaussian beam methods for the Schrödinger equation with discontinuous potentials
Noname manuscript No. (will be inserted by the editor) Gaussian beam methods for the Schrödinger equation with discontinuous potentials Shi Jin Dongming Wei Dongsheng Yin Received: date / Accepted: date
More informationOn compact operators
On compact operators Alen Aleanderian * Abstract In this basic note, we consider some basic properties of compact operators. We also consider the spectrum of compact operators on Hilbert spaces. A basic
More informationA Padé approximation to the scalar wavefield extrapolator for inhomogeneous media
A Padé approimation A Padé approimation to the scalar wavefield etrapolator for inhomogeneous media Yanpeng Mi, Zhengsheng Yao, and Gary F. Margrave ABSTRACT A seismic wavefield at depth z can be obtained
More informationON THE BLOCH DECOMPOSITION BASED SPECTRAL METHOD FOR WAVE PROPAGATION IN PERIODIC MEDIA. version: February 15,
ON THE BLOCH DECOMPOSITION BASED SPECTRAL METHOD FOR WAVE PROPAGATION IN PERIODIC MEDIA ZHONGYI HUANG, SHI JIN, PETER A. MARKOWICH, AND CHRISTOF SPARBER Abstract. We extend the Bloch-decomposition based
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 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 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 informationM445: Heat equation with sources
M5: Heat equation with sources David Gurarie I. On Fourier and Newton s cooling laws The Newton s law claims the temperature rate to be proportional to the di erence: d dt T = (T T ) () The Fourier law
More informationFROZEN GAUSSIAN APPROXIMATION FOR HIGH FREQUENCY WAVE PROPAGATION
FROZEN GAUSSIAN APPROXIMATION FOR HIGH FREQUENCY WAVE PROPAGATION JIANFENG LU AND XU YANG Abstract. We propose the frozen Gaussian approimation for computation of high frequency wave propagation. This
More informationPROBLEMS In each of Problems 1 through 12:
6.5 Impulse Functions 33 which is the formal solution of the given problem. It is also possible to write y in the form 0, t < 5, y = 5 e (t 5/ sin 5 (t 5, t 5. ( The graph of Eq. ( is shown in Figure 6.5.3.
More informationThe effect of Group Velocity in the numerical analysis of control problems for the wave equation
The effect of Group Velocity in the numerical analysis of control problems for the wave equation Fabricio Macià École Normale Supérieure, D.M.A., 45 rue d Ulm, 753 Paris cedex 5, France. Abstract. In this
More informationProperties of the Scattering Transform on the Real Line
Journal of Mathematical Analysis and Applications 58, 3 43 (001 doi:10.1006/jmaa.000.7375, available online at http://www.idealibrary.com on Properties of the Scattering Transform on the Real Line Michael
More informationA Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
More informationA. Time dependence from separation of timescales
Lecture 4 Adiabatic Theorem So far we have considered time independent semiclassical problems. What makes these semiclassical is that the gradient term (which is multiplied by 2 ) was small. In other words,
More informationStatistical Geometry Processing Winter Semester 2011/2012
Statistical Geometry Processing Winter Semester 2011/2012 Linear Algebra, Function Spaces & Inverse Problems Vector and Function Spaces 3 Vectors vectors are arrows in space classically: 2 or 3 dim. Euclidian
More informationChapter 29. Quantum Chaos
Chapter 29 Quantum Chaos What happens to a Hamiltonian system that for classical mechanics is chaotic when we include a nonzero h? There is no problem in principle to answering this question: given a classical
More informationWave propagation in discrete heterogeneous media
www.bcamath.org Wave propagation in discrete heterogeneous media Aurora MARICA marica@bcamath.org BCAM - Basque Center for Applied Mathematics Derio, Basque Country, Spain Summer school & workshop: PDEs,
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 informationTaylor Series and Asymptotic Expansions
Taylor Series and Asymptotic Epansions The importance of power series as a convenient representation, as an approimation tool, as a tool for solving differential equations and so on, is pretty obvious.
More informationComputational High Frequency Wave Propagation
Computational High Frequency Wave Propagation Olof Runborg CSC, KTH Isaac Newton Institute Cambridge, February 2007 Olof Runborg (KTH) High-Frequency Waves INI, 2007 1 / 52 Outline 1 Introduction, background
More informationFourier transform of tempered distributions
Fourier transform of tempered distributions 1 Test functions and distributions As we have seen before, many functions are not classical in the sense that they cannot be evaluated at any point. For eample,
More informationSemiclassical limit of the Schrödinger-Poisson-Landau-Lifshitz-Gilbert system
Semiclassical limit of the Schrödinger-Poisson-Landau-Lifshitz-Gilbert system Lihui Chai Department of Mathematics University of California, Santa Barbara Joint work with Carlos J. García-Cervera, and
More informationSemi-Eulerian and High Order Gaussian Beam Methods for the Schrödinger Equation in the Semiclassical Regime
Commun. Comput. Phys. doi:.48/cicp.99.63s Vol. 9, No. 3, pp. 668-687 March Semi-Eulerian and High Order Gaussian Beam Methods for the Schrödinger Equation in the Semiclassical Regime Shi Jin,, Hao Wu and
More informationON EFFECTIVE HAMILTONIANS FOR ADIABATIC PERTURBATIONS. Mouez Dimassi Jean-Claude Guillot James Ralston
ON EFFECTIVE HAMILTONIANS FOR ADIABATIC PERTURBATIONS Mouez Dimassi Jean-Claude Guillot James Ralston Abstract We construct almost invariant subspaces and the corresponding effective Hamiltonian for magnetic
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 informationYURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL
Journal of Comput. & Applied Mathematics 139(2001), 197 213 DIRECT APPROACH TO CALCULUS OF VARIATIONS VIA NEWTON-RAPHSON METHOD YURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL Abstract. Consider m functions
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 informationHarmonic Analysis Homework 5
Harmonic Analysis Homework 5 Bruno Poggi Department of Mathematics, University of Minnesota November 4, 6 Notation Throughout, B, r is the ball of radius r with center in the understood metric space usually
More informationversion: January 21, Physical motivation
NUMERICAL SIMULATION OF THE NONLINEAR SCHRÖDINGER EQUATION WITH MULTI-DIMENSIONAL PERIODIC POTENTIALS ZHONGYI HUANG, SHI JIN, PETER A. MARKOWICH, AND CHRISTOF SPARBER Abstract. By extending the Bloch-decomposition
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationDedicated to Professor Georege Papanicolaou on the occasion of his 60th birthday
METHODS AND APPLICATIONS OF ANALYSIS. c 4 International Press Vol., No. 3, pp., September 4 6 NUMERICAL METHODS FOR THE NONLINEAR SCHRÖDINGER EQUATION WITH NONZERO FAR-FIELD CONDITIONS WEIZHU BAO Dedicated
More informationFourier Analysis Fourier Series C H A P T E R 1 1
C H A P T E R Fourier Analysis 474 This chapter on Fourier analysis covers three broad areas: Fourier series in Secs...4, more general orthonormal series called Sturm iouville epansions in Secs..5 and.6
More informationarxiv:math/ v1 [math.ap] 1 Jan 1992
APPEARED IN BULLETIN OF THE AMERICAN MATHEMATICAL SOCIETY Volume 26, Number 1, Jan 1992, Pages 119-124 A STEEPEST DESCENT METHOD FOR OSCILLATORY RIEMANN-HILBERT PROBLEMS arxiv:math/9201261v1 [math.ap]
More informationRe-Thinking Schroedinger Boundary Conditions: the Simple Rotator is Not Simple
Re-Thinking Schroedinger Boundary Conditions: the Simple Rotator is Not Simple Arthur Davidson ECE Dept., Carnegie Mellon University Pittsburgh, PA 1513 artdav@ece.cmu.edu ABSTRACT A contradiction arises
More information7.4. Why we have two different types of materials: conductors and insulators?
Phys463.nb 55 7.3.5. Folding, Reduced Brillouin zone and extended Brillouin zone for free particles without lattices In the presence of a lattice, we can also unfold the extended Brillouin zone to get
More informationTopics for the Qualifying Examination
Topics for the Qualifying Examination Quantum Mechanics I and II 1. Quantum kinematics and dynamics 1.1 Postulates of Quantum Mechanics. 1.2 Configuration space vs. Hilbert space, wave function vs. state
More informationSmall BGK waves and nonlinear Landau damping (higher dimensions)
Small BGK waves and nonlinear Landau damping higher dimensions Zhiwu Lin and Chongchun Zeng School of Mathematics Georgia Institute of Technology Atlanta, GA, USA Abstract Consider Vlasov-Poisson system
More informationTWO-SCALE ASYMPTOTIC EXPANSION : BLOCH-FLOQUET THEORY IN PERIODIC STRUCTURES
TWO-SCALE ASYMPTOTIC EXPANSION : BLOCH-FLOQUET THEORY IN PERIODIC STRUCTURES JACK ARBUNICH 1. Setting of Periodic Structures Our aim is to study an application of Bloch-Floquet Theory in the multiscale
More informationMixing Quantum and Classical Mechanics: A Partially Miscible Solution
Mixing Quantum and Classical Mechanics: A Partially Miscible Solution R. Kapral S. Nielsen A. Sergi D. Mac Kernan G. Ciccotti quantum dynamics in a classical condensed phase environment how to simulate
More informationNonlinear Geometric Optics Based Multiscale Stochastic Galerkin Methods for Highly Oscillatory Transport Equations with Random Inputs
Nonlinear Geometric Optics Based Multiscale Stochastic Galerkin Methods for Highly Oscillatory Transport Equations with Random Inputs Nicolas Crouseilles, Shi Jin, Mohammed Lemou, Liu Liu February 25,
More informationGeometric Modeling Summer Semester 2010 Mathematical Tools (1)
Geometric Modeling Summer Semester 2010 Mathematical Tools (1) Recap: Linear Algebra Today... Topics: Mathematical Background Linear algebra Analysis & differential geometry Numerical techniques Geometric
More informationCLASSICAL LIMIT FOR SEMI-RELATIVISTIC HARTREE SYSTEMS
CLASSICAL LIMIT FOR SEMI-RELATIVISTIC HARTREE SYSTEMS GONCA L. AKI, PETER A. MARKOWICH, AND CHRISTOF SPARBER Abstract. We consider the three-dimensional semi-relativistic Hartree model for fast quantum
More informationBlack holes as open quantum systems
Black holes as open quantum systems Claus Kiefer Institut für Theoretische Physik Universität zu Köln Hawking radiation 1 1 singularity II γ H γ γ H collapsing 111 star 1 1 I - future event horizon + i
More informationQuantum Dynamics. March 10, 2017
Quantum Dynamics March 0, 07 As in classical mechanics, time is a parameter in quantum mechanics. It is distinct from space in the sense that, while we have Hermitian operators, X, for position and therefore
More informationSemiclassical Electron Transport
Semiclassical Electron Transport Branislav K. Niolić Department of Physics and Astronomy, University of Delaware, U.S.A. PHYS 64: Introduction to Solid State Physics http://www.physics.udel.edu/~bniolic/teaching/phys64/phys64.html
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationOPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS
Stochastics and Dynamics, Vol. 2, No. 3 (22 395 42 c World Scientific Publishing Company OPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS Stoch. Dyn. 22.2:395-42. Downloaded from www.worldscientific.com by HARVARD
More informationarxiv:quant-ph/ v1 29 Mar 2003
Finite-Dimensional PT -Symmetric Hamiltonians arxiv:quant-ph/0303174v1 29 Mar 2003 Carl M. Bender, Peter N. Meisinger, and Qinghai Wang Department of Physics, Washington University, St. Louis, MO 63130,
More information1. Introduction We are interested in developing the Gaussian beam method for the Dirac equation in the semi-classical regime
GAUSSIAN BEAM METHODS FOR THE DIRAC EQUATION IN THE SEMI-CLASSICAL REGIME HAO WU, ZHONGYI HUANG, SHI JIN, AND DONGSHENG YIN Abstract. The Dirac equation is an important model in relativistic quantum mechanics.
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 informationSUMMARY INTRODUCTION THEORY
Seismic modeling using the frozen Gaussian approximation Xu Yang, University of California Santa Barbara, Jianfeng Lu, Duke University, and Sergey Fomel, University of Texas at Austin SUMMARY We adopt
More informationOPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS
Stochastics and Dynamics c World Scientific Publishing Company OPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS BRIAN F. FARRELL Division of Engineering and Applied Sciences, Harvard University Pierce Hall, 29
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationSYLLABUS FOR ENTRANCE EXAMINATION NANYANG TECHNOLOGICAL UNIVERSITY FOR INTERNATIONAL STUDENTS A-LEVEL MATHEMATICS
SYLLABUS FOR ENTRANCE EXAMINATION NANYANG TECHNOLOGICAL UNIVERSITY FOR INTERNATIONAL STUDENTS A-LEVEL MATHEMATICS STRUCTURE OF EXAMINATION PAPER. There will be one -hour paper consisting of 4 questions..
More informationSemi-Eulerian and High Order Gaussian Beam Methods for the Schrödinger Equation in the Semiclassical Regime
Semi-Eulerian and High Order Gaussian Beam Methods for the Schrödinger Equation in the Semiclassical Regime Shi Jin, Hao Wu, and Xu Yang October 9, 29 Abstract A novel Eulerian Gaussian beam method was
More informationNumerical techniques. Chapter Difference equations
Chapter 6 Numerical techniques The differential equations (4.61), (4.62) and (4.64), in combination with boundary conditions such as equations (4.65) (4.68), constitute a two point boundary value problem.
More informationWaves in Honeycomb Structures
Waves in Honeycomb Structures Michael I. Weinstein Columbia University Nonlinear Schrödinger Equations: Theory and Applications Heraklion, Crete / Greece May 20-24, 2013 Joint work with C.L. Fefferman
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More informationGaussian Beam Approximations
Gaussian Beam Approximations Olof Runborg CSC, KTH Joint with Hailiang Liu, Iowa, and Nick Tanushev, Austin Technischen Universität München München, January 2011 Olof Runborg (KTH) Gaussian Beam Approximations
More informationQuantized vortex stability and interaction in the nonlinear wave equation
Physica D 237 (2008) 2391 2410 www.elsevier.com/locate/physd Quantized vortex stability and interaction in the nonlinear wave equation Weizhu Bao a,b,, Rong Zeng c, Yanzhi Zhang a,1 a Department of Mathematics,
More informationBose Gases, Bose Einstein Condensation, and the Bogoliubov Approximation
Bose Gases, Bose Einstein Condensation, and the Bogoliubov Approximation Robert Seiringer IST Austria Mathematical Horizons for Quantum Physics IMS Singapore, September 18, 2013 R. Seiringer Bose Gases,
More informationPhysics 127b: Statistical Mechanics. Landau Theory of Second Order Phase Transitions. Order Parameter
Physics 127b: Statistical Mechanics Landau Theory of Second Order Phase Transitions Order Parameter Second order phase transitions occur when a new state of reduced symmetry develops continuously from
More informationChapter 3 Burgers Equation
Chapter 3 Burgers Equation One of the major challenges in the field of comple systems is a thorough understanding of the phenomenon of turbulence. Direct numerical simulations (DNS) have substantially
More informationMINIMAL NORMAL AND COMMUTING COMPLETIONS
INTERNATIONAL JOURNAL OF INFORMATION AND SYSTEMS SCIENCES Volume 4, Number 1, Pages 5 59 c 8 Institute for Scientific Computing and Information MINIMAL NORMAL AND COMMUTING COMPLETIONS DAVID P KIMSEY AND
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 informationPHYS Concept Tests Fall 2009
PHYS 30 Concept Tests Fall 009 In classical mechanics, given the state (i.e. position and velocity) of a particle at a certain time instant, the state of the particle at a later time A) cannot be determined
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationOutline. Escaping the Double Cone: Describing the Seam Space with Gateway Modes. Columbus Workshop Argonne National Laboratory 2005
Outline Escaping the Double Cone: Describing the Seam Space with Gateway Modes Columbus Worshop Argonne National Laboratory 005 Outline I. Motivation II. Perturbation Theory and Intersection Adapted Coordinates
More informationNumerical methods for hyperbolic systems with singular coefficients: well-balanced scheme, Hamiltonian preservation, and beyond
Proceedings of Symposia in Applied Mathematics Numerical methods for hyperbolic systems with singular coefficients: well-balanced scheme, Hamiltonian preservation, and beyond Shi Jin Abstract. This paper
More informationVisualizing the kinematics of relativistic wave packets
Visualizing the kinematics of relativistic wave packets Bernd Thaller Institute for Mathematics and Scientific Computation, University of Graz Heinrichstrasse 36, A-81 Graz, Austria Abstract. This article
More informationZEROS OF THE JOST FUNCTION FOR A CLASS OF EXPONENTIALLY DECAYING POTENTIALS. 1. Introduction
Electronic Journal of Differential Equations, Vol. 20052005), No. 45, pp. 9. ISSN: 072-669. URL: http://ejde.math.tstate.edu or http://ejde.math.unt.edu ftp ejde.math.tstate.edu login: ftp) ZEROS OF THE
More informationLecture 10: Whitham Modulation Theory
Lecture 10: Whitham Modulation Theory Lecturer: Roger Grimshaw. Write-up: Andong He June 19, 2009 1 Introduction The Whitham modulation theory provides an asymptotic method for studying slowly varying
More informationQuantum Physics III (8.06) Spring 2016 Assignment 3
Quantum Physics III (8.6) Spring 6 Assignment 3 Readings Griffiths Chapter 9 on time-dependent perturbation theory Shankar Chapter 8 Cohen-Tannoudji, Chapter XIII. Problem Set 3. Semi-classical approximation
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationFOURIER INVERSION. an additive character in each of its arguments. The Fourier transform of f is
FOURIER INVERSION 1. The Fourier Transform and the Inverse Fourier Transform Consider functions f, g : R n C, and consider the bilinear, symmetric function ψ : R n R n C, ψ(, ) = ep(2πi ), an additive
More informationC2: Band structure. Carl-Olof Almbladh, Rikard Nelander, and Jonas Nyvold Pedersen Department of Physics, Lund University.
C2: Band structure Carl-Olof Almbladh, Rikard Nelander, and Jonas Nyvold Pedersen Department of Physics, Lund University December 2005 1 Introduction When you buy a diamond for your girl/boy friend, what
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationResearch Article The Approximate Solution of Fredholm Integral Equations with Oscillatory Trigonometric Kernels
Applied Mathematics, Article ID 72327, 7 pages http://ddoiorg/055/204/72327 Research Article The Approimate Solution of Fredholm Integral Equations with Oscillatory Trigonometric Kernels Qinghua Wu School
More information2 Formal derivation of the Shockley-Read-Hall model
We consider a semiconductor crystal represented by the bounded domain R 3 (all our results are easily extended to the one and two- dimensional situations) with a constant (in space) number density of traps
More informationSemi-Classical Dynamics Using Hagedorn Wavepackets
Semi-Classical Dynamics Using Hagedorn Wavepackets Leila Taghizadeh June 6, 2013 Leila Taghizadeh () Semi-Classical Dynamics Using Hagedorn Wavepackets June 6, 2013 1 / 56 Outline 1 The Schrödinger Equation
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationSECOND ORDER ALL SPEED METHOD FOR THE ISENTROPIC EULER EQUATIONS. Min Tang. (Communicated by the associate editor name)
Manuscript submitted to AIMS Journals Volume X, Number X, XX 2X Website: http://aimsciences.org pp. X XX SECOND ORDER ALL SPEED METHOD FOR THE ISENTROPIC EULER EQUATIONS Min Tang Department of Mathematics
More information1. Introduction We are interested in the Gaussian beam methods for the numerical approximation of the Schrödinger equation:
COMMUN. MATH. SCI. Vol. 6, No. 4, pp. 995 2 c 28 International Press GAUSSIAN BEAM METHODS FOR THE SCHRÖDINGER EQUATION IN THE SEMI-CLASSICAL REGIME: LAGRANGIAN AND EULERIAN FORMULATIONS SHI JIN, HAO WU,
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 information