Schrödinger operators exhibiting parameter-dependent spectral transitions

Size: px
Start display at page:

Download "Schrödinger operators exhibiting parameter-dependent spectral transitions"

Transcription

1 Schrödinger operators exhibiting parameter-dependent spectral transitions Pavel Exner Doppler Institute for Mathematical Physics and Applied Mathematics Prague in collaboration with Diana Barseghyan, Andrii Khrabustovskyi, Vladimir Lotoreichik, and Miloš Tater A talk at the The fourth Najman Conference on Spectral Problems for Operators and Matrices Opatija, September 22, 2015 P. Exner: Spectral transitions... 4th Najman Conference September 22,

2 Talk outline The simplest example: Smilansky model P. Exner: Spectral transitions... 4th Najman Conference September 22,

3 Talk outline The simplest example: Smilansky model Spectral properties, the subcritical and supercritical case Numerical solution P. Exner: Spectral transitions... 4th Najman Conference September 22,

4 Talk outline The simplest example: Smilansky model Spectral properties, the subcritical and supercritical case Numerical solution A regular version of Smilansky model Choice of the potential Spectral properties, the subcritical and supercritical case Intervals and multiple channels P. Exner: Spectral transitions... 4th Najman Conference September 22,

5 Talk outline The simplest example: Smilansky model Spectral properties, the subcritical and supercritical case Numerical solution A regular version of Smilansky model Choice of the potential Spectral properties, the subcritical and supercritical case Intervals and multiple channels Another model: transition from a purely discrete to the real line Potential for which Weyl quantization fails Spectral properties, the subcritical and supercritical case The critical case P. Exner: Spectral transitions... 4th Najman Conference September 22,

6 Talk outline The simplest example: Smilansky model Spectral properties, the subcritical and supercritical case Numerical solution A regular version of Smilansky model Choice of the potential Spectral properties, the subcritical and supercritical case Intervals and multiple channels Another model: transition from a purely discrete to the real line Potential for which Weyl quantization fails Spectral properties, the subcritical and supercritical case The critical case Summary & open questions P. Exner: Spectral transitions... 4th Najman Conference September 22,

7 Smilansky model The model was originally proposed in [Smilansky 04] to describe a one-dimensional system interacting with a caricature heat bath represented by a harmonic oscillator. P. Exner: Spectral transitions... 4th Najman Conference September 22,

8 Smilansky model The model was originally proposed in [Smilansky 04] to describe a one-dimensional system interacting with a caricature heat bath represented by a harmonic oscillator. Mathematical properties of the model were analyzed in [Solomyak 04], [Evans-Solomyak 05], [Naboko-Solomyak 06]. More recently, time evolution in such a (slightly modified) model was analyzed [Guarneri 11] P. Exner: Spectral transitions... 4th Najman Conference September 22,

9 Smilansky model The model was originally proposed in [Smilansky 04] to describe a one-dimensional system interacting with a caricature heat bath represented by a harmonic oscillator. Mathematical properties of the model were analyzed in [Solomyak 04], [Evans-Solomyak 05], [Naboko-Solomyak 06]. More recently, time evolution in such a (slightly modified) model was analyzed [Guarneri 11] In PDE terms, the model is described through a 2D Schrödinger operator H Sm = 2 x ) ( 2 2 y 2 + y 2 + λyδ(x) on L 2 (R) with various modifications to be mentioned later. P. Exner: Spectral transitions... 4th Najman Conference September 22,

10 Smilansky model The model was originally proposed in [Smilansky 04] to describe a one-dimensional system interacting with a caricature heat bath represented by a harmonic oscillator. Mathematical properties of the model were analyzed in [Solomyak 04], [Evans-Solomyak 05], [Naboko-Solomyak 06]. More recently, time evolution in such a (slightly modified) model was analyzed [Guarneri 11] In PDE terms, the model is described through a 2D Schrödinger operator H Sm = 2 x ) ( 2 2 y 2 + y 2 + λyδ(x) on L 2 (R) with various modifications to be mentioned later. Due to a particular choice of the coupling the model exhibited a spectral transition with respect to the coupling parameter λ. P. Exner: Spectral transitions... 4th Najman Conference September 22,

11 A summary of results about the model Spectral transition: if λ > 2 the particle can escape to infinity along the singular channel in the y direction. In spectral terms, it corresponds to switch from a positive to a below unbounded spectrum at λ = 2. P. Exner: Spectral transitions... 4th Najman Conference September 22,

12 A summary of results about the model Spectral transition: if λ > 2 the particle can escape to infinity along the singular channel in the y direction. In spectral terms, it corresponds to switch from a positive to a below unbounded spectrum at λ = 2. At the heuristic level, the mechanism is easy to understand: we have an effective variable decoupling far from the x-axis and the oscillator potential competes there with the δ interaction eigenvalue 1 4 λ2 y 2. P. Exner: Spectral transitions... 4th Najman Conference September 22,

13 A summary of results about the model Spectral transition: if λ > 2 the particle can escape to infinity along the singular channel in the y direction. In spectral terms, it corresponds to switch from a positive to a below unbounded spectrum at λ = 2. At the heuristic level, the mechanism is easy to understand: we have an effective variable decoupling far from the x-axis and the oscillator potential competes there with the δ interaction eigenvalue 1 4 λ2 y 2. Eigenvalue absence: for any λ 0 there are no eigenvalues 1 2. If λ > 2, the point spectrum of H Sm is empty. P. Exner: Spectral transitions... 4th Najman Conference September 22,

14 A summary of results about the model Spectral transition: if λ > 2 the particle can escape to infinity along the singular channel in the y direction. In spectral terms, it corresponds to switch from a positive to a below unbounded spectrum at λ = 2. At the heuristic level, the mechanism is easy to understand: we have an effective variable decoupling far from the x-axis and the oscillator potential competes there with the δ interaction eigenvalue 1 4 λ2 y 2. Eigenvalue absence: for any λ 0 there are no eigenvalues 1 2. If λ > 2, the point spectrum of H Sm is empty. Existence of eigenvalues: for 0 < λ < 2 we have H Sm 0. The point spectrum is nonempty and finite, and N( 1 2, H Sm) 1 4 2(µ(λ) 1) holds as λ 2, where µ(λ) := 2/λ. P. Exner: Spectral transitions... 4th Najman Conference September 22,

15 Further results Absolute continuity: in the supercritical case λ > 2 we have σ ac (H Sm ) = R P. Exner: Spectral transitions... 4th Najman Conference September 22,

16 Further results Absolute continuity: in the supercritical case λ > 2 we have σ ac (H Sm ) = R Extension of the result to a two channel case with different oscillator frequencies [Evans-Solomyak 05] P. Exner: Spectral transitions... 4th Najman Conference September 22,

17 Further results Absolute continuity: in the supercritical case λ > 2 we have σ ac (H Sm ) = R Extension of the result to a two channel case with different oscillator frequencies [Evans-Solomyak 05] Extension to multiple channels on a system periodic in x [Guarneri 11]. In this paper the time evolution generated by H Sm is investigated and proposed as a model of wavepacket collapse. P. Exner: Spectral transitions... 4th Najman Conference September 22,

18 Further results Absolute continuity: in the supercritical case λ > 2 we have σ ac (H Sm ) = R Extension of the result to a two channel case with different oscillator frequencies [Evans-Solomyak 05] Extension to multiple channels on a system periodic in x [Guarneri 11]. In this paper the time evolution generated by H Sm is investigated and proposed as a model of wavepacket collapse. The above results have been obtained by a combination of different methods: a reduction to an infinite system of ODE s, facts from Jacobi matrices theory, variational estimates, etc. P. Exner: Spectral transitions... 4th Najman Conference September 22,

19 Further results Absolute continuity: in the supercritical case λ > 2 we have σ ac (H Sm ) = R Extension of the result to a two channel case with different oscillator frequencies [Evans-Solomyak 05] Extension to multiple channels on a system periodic in x [Guarneri 11]. In this paper the time evolution generated by H Sm is investigated and proposed as a model of wavepacket collapse. The above results have been obtained by a combination of different methods: a reduction to an infinite system of ODE s, facts from Jacobi matrices theory, variational estimates, etc. Before proceeding further, let show how the spectrum can be treated numerically in the subcritical case. P. Exner: Spectral transitions... 4th Najman Conference September 22,

20 Numerical search for eigenvalues In the halfplanes ±x > 0 the wave functions can be expanded using the transverse base spanned by the functions ψ n (y) = 1 2 n n! π e y 2 /2 H n (y) corresponding to the oscillator eigenvalues n + 1 2, n = 0, 1, 2,.... P. Exner: Spectral transitions... 4th Najman Conference September 22,

21 Numerical search for eigenvalues In the halfplanes ±x > 0 the wave functions can be expanded using the transverse base spanned by the functions ψ n (y) = 1 2 n n! π e y 2 /2 H n (y) corresponding to the oscillator eigenvalues n + 1 2, n = 0, 1, 2,.... Furthermore, one can make use of the mirror symmetry w.r.t. x = 0 and divide H λ into the trivial odd part H ( ) λ and the even part H (+) λ which is equivalent to the operator on L 2 (R (0, )) with the same symbol determined by the boundary condition f x (0+, y) = 1 αyf (0+, y). 2 P. Exner: Spectral transitions... 4th Najman Conference September 22,

22 Numerical solution, continued We substitute the Ansatz with κ n := n ɛ. f (x, y) = c n e κnx ψ n (y) n=0 P. Exner: Spectral transitions... 4th Najman Conference September 22,

23 Numerical solution, continued We substitute the Ansatz with κ n := n ɛ. f (x, y) = c n e κnx ψ n (y) n=0 This yields for solution with the energy ɛ the equation B λ c = 0, where c is the coefficient vector and B λ is the operator in l 2 with (B λ ) m,n = κ n δ m,n λ(ψ m, yψ n ). P. Exner: Spectral transitions... 4th Najman Conference September 22,

24 Numerical solution, continued We substitute the Ansatz with κ n := n ɛ. f (x, y) = c n e κnx ψ n (y) n=0 This yields for solution with the energy ɛ the equation B λ c = 0, where c is the coefficient vector and B λ is the operator in l 2 with (B λ ) m,n = κ n δ m,n λ(ψ m, yψ n ). Note that the matrix is in fact tridiagonal because (ψ m, yψ n ) = 1 2 ( n + 1 δm,n+1 + n δ m,n 1 ). P. Exner: Spectral transitions... 4th Najman Conference September 22,

25 Smilansky model eigenvalues 0.4 E n λ In most part of the subcritical region there is a single eigenvalue, the second one appears only at λ P. Exner: Spectral transitions... 4th Najman Conference September 22,

26 Smilansky model eigenvalues 0.4 E n λ In most part of the subcritical region there is a single eigenvalue, the second one appears only at λ The next thresholds are , , , ,... P. Exner: Spectral transitions... 4th Najman Conference September 22,

27 Smilansky model eigenvalues Close to the critical value, however, many eigenvalues appear which gradually fill the interval (0, 1 2 ) as the critical value is approached P. Exner: Spectral transitions... 4th Najman Conference September 22,

28 Their number is as predicted The dots mean the eigenvalue numbers, the red curve is the above mentioned asymptotics due to Solomyak P. Exner: Spectral transitions... 4th Najman Conference September 22,

29 Smilansky model ground state The numerical solution also indicates other properties, for instance, that the first eigenvalue behaves as ɛ 1 = 1 2 cλ4 + o(λ 4 ) as λ 0, with c P. Exner: Spectral transitions... 4th Najman Conference September 22,

30 In fact, we have c = Indeed, the relation B λ c = 0 can be written explicitly as µλ c0 λ + λ 2 2 cλ 1 = 0, kλ 2 2 cλ k 1 + k + 1λ k + µ λ ck λ + 2 ck+1 λ 2 = 0, k 1, where µ λ := 1 2 E 1(λ) and c λ = {c0 λ, cλ 1,... } is the corresponding normalized eigenvector of B λ. P. Exner: Spectral transitions... 4th Najman Conference September 22,

31 In fact, we have c = Indeed, the relation B λ c = 0 can be written explicitly as µλ c0 λ + λ 2 2 cλ 1 = 0, kλ 2 2 cλ k 1 + k + 1λ k + µ λ ck λ + 2 ck+1 λ 2 = 0, k 1, where µ λ := 1 2 E 1(λ) and c λ = {c0 λ, cλ 1,... } is the corresponding normalized eigenvector of B λ. Using the above relations and simple estimates, we get ck λ λ2 and c0 λ = 1 + O(λ 2 ) k=1 as λ 0+; hence we have in particular c λ 1 = λ O(λ2 ). P. Exner: Spectral transitions... 4th Najman Conference September 22,

32 In fact, we have c = The first of the above relation then gives as λ 0+, in other words µ λ = λ O(λ5 ) E 1 (λ) = 1 2 λ O(λ5 ). P. Exner: Spectral transitions... 4th Najman Conference September 22,

33 In fact, we have c = The first of the above relation then gives as λ 0+, in other words µ λ = λ O(λ5 ) E 1 (λ) = 1 2 λ O(λ5 ). And the mentioned coefficient is nothing else than P. Exner: Spectral transitions... 4th Najman Conference September 22,

34 Smilansky model eigenfunctions 0 1 = = y x = x = y x 10 5 x = = x x The first six eigenfunctions of HSm for λ = , in other words, λ = P. Exner: Spectral transitions... 4th Najman Conference September 22,

35 A regular version of Smilansky model Our next aim is to show that one can observe a similar effect for Schrödinger operators with regular potentials. Now, however, the coupling cannot be linear in y and the profile of the channel has to change with y. P. Exner: Spectral transitions... 4th Najman Conference September 22,

36 A regular version of Smilansky model Our next aim is to show that one can observe a similar effect for Schrödinger operators with regular potentials. Now, however, the coupling cannot be linear in y and the profile of the channel has to change with y. Recall that the effect comes from competition between the oscillator potential with the principal eigenvalue of the transverse part of the operator equal to 1 4 λ2 y 2. P. Exner: Spectral transitions... 4th Najman Conference September 22,

37 A regular version of Smilansky model Our next aim is to show that one can observe a similar effect for Schrödinger operators with regular potentials. Now, however, the coupling cannot be linear in y and the profile of the channel has to change with y. Recall that the effect comes from competition between the oscillator potential with the principal eigenvalue of the transverse part of the operator equal to 1 4 λ2 y 2. We replace the δ by a family of shrinking potentials whose mean matches the δ coupling constant, U(x, y) dx y. This can be achieved, e.g., by choosing U(x, y) = λy 2 V (xy) for a fixed function V. P. Exner: Spectral transitions... 4th Najman Conference September 22,

38 A regular version of Smilansky model Our next aim is to show that one can observe a similar effect for Schrödinger operators with regular potentials. Now, however, the coupling cannot be linear in y and the profile of the channel has to change with y. Recall that the effect comes from competition between the oscillator potential with the principal eigenvalue of the transverse part of the operator equal to 1 4 λ2 y 2. We replace the δ by a family of shrinking potentials whose mean matches the δ coupling constant, U(x, y) dx y. This can be achieved, e.g., by choosing U(x, y) = λy 2 V (xy) for a fixed function V. This motivates us to investigate the following operator on L 2 (R 2 ), H = 2 x 2 2 y 2 + ω2 y 2 λy 2 V (xy)χ { x a} (x), where ω, a are positive constants, χ { y a} is the indicator function of the interval ( a, a), and the potential V with supp V [ a, a] is a nonnegative function with bounded first derivative. P. Exner: Spectral transitions... 4th Najman Conference September 22,

39 A regular version of Smilansky model, continued By Faris-Lavine theorem the operator is e.s.a. on C0 same is true for its generalization, (R2 ) and the H = 2 x 2 2 y 2 + ω2 y 2 N λ j y 2 V j (xy)χ { x bj a j }(x) with a finite number of channels, where functions V j are positive with bounded first derivative, with the supports contained in (b j a j, b j + a j ) and such that supp V j supp V k = holds for j k. j=1 P. Exner: Spectral transitions... 4th Najman Conference September 22,

40 A regular version of Smilansky model, continued By Faris-Lavine theorem the operator is e.s.a. on C0 same is true for its generalization, (R2 ) and the H = 2 x 2 2 y 2 + ω2 y 2 N λ j y 2 V j (xy)χ { x bj a j }(x) with a finite number of channels, where functions V j are positive with bounded first derivative, with the supports contained in (b j a j, b j + a j ) and such that supp V j supp V k = holds for j k. Remark We note that the properties discussed below depend on the asymptotic behavior of the potential channels and would not change if the potential is modified in the vicinity of the x-axis, for instance, by replacing the above cut-off functions with χ y a (y) and χ y aj (y), respectively. j=1 P. Exner: Spectral transitions... 4th Najman Conference September 22,

41 Subcritical case To state the result we employ a 1D comparison operator L = L V, L = d2 dx 2 + ω2 λv (x) on L 2 (R) with the domain H 2 (R). What matters is the sign of its spectral threshold; since V is supposed to be nonnegative, the latter is a monotonous function of λ and there is a λ crit > 0 at which the sign changes. P. Exner: Spectral transitions... 4th Najman Conference September 22,

42 Subcritical case To state the result we employ a 1D comparison operator L = L V, L = d2 dx 2 + ω2 λv (x) on L 2 (R) with the domain H 2 (R). What matters is the sign of its spectral threshold; since V is supposed to be nonnegative, the latter is a monotonous function of λ and there is a λ crit > 0 at which the sign changes. Theorem (Barseghyan-E 14) Under the stated assumption, the spectrum of the operator H is bounded from below provided the operator L is positive. P. Exner: Spectral transitions... 4th Najman Conference September 22,

43 Proof outline It is sufficient to prove the claim for λ = 1. We employ Neumann bracketing, similarly as for the previous model. P. Exner: Spectral transitions... 4th Najman Conference September 22,

44 Proof outline It is sufficient to prove the claim for λ = 1. We employ Neumann bracketing, similarly as for the previous model. Let h n and h n be respectively the restrictions of operator H to the strips G n = R {y : ln n < y ln(n + 1)}, n = 1, 2,..., and G n, their mirror images w.r.t. y, with Neumann boundary conditions; then H (h n h n ) ; n=1 We find a uniform lower bound σ(h n ) and σ( h n ) as n. P. Exner: Spectral transitions... 4th Najman Conference September 22,

45 Proof outline It is sufficient to prove the claim for λ = 1. We employ Neumann bracketing, similarly as for the previous model. Let h n and h n be respectively the restrictions of operator H to the strips G n = R {y : ln n < y ln(n + 1)}, n = 1, 2,..., and G n, their mirror images w.r.t. y, with Neumann boundary conditions; then H (h n h n ) ; n=1 We find a uniform lower bound σ(h n ) and σ( h n ) as n. Using the assumptions about V we find ( ) ( ) 1 ln n V (xy) V (x ln n) = O, y 2 ln 2 n = O n ln n n for any (x, y) G n, and analogous relations for G n. P. Exner: Spectral transitions... 4th Najman Conference September 22,

46 Proof outline continued This yields for for any (x, y) G n ( ) ln n y 2 V (xy) ln 2 n V (x ln n) = O n P. Exner: Spectral transitions... 4th Najman Conference September 22,

47 Proof outline continued This yields ( ) ln n y 2 V (xy) ln 2 n V (x ln n) = O n for for any (x, y) G n which allows us to check that ( ) ln n inf σ(h n ) inf σ(l n ) + O, n where l n := 2 x 2 2 y 2 + ω 2 ln 2 n ln 2 n V (x ln n) on L 2 (G n ). P. Exner: Spectral transitions... 4th Najman Conference September 22,

48 Proof outline continued This yields ( ) ln n y 2 V (xy) ln 2 n V (x ln n) = O n for for any (x, y) G n which allows us to check that ( ) ln n inf σ(h n ) inf σ(l n ) + O, n where l n := 2 x 2 2 y 2 + ω 2 ln 2 n ln 2 n V (x ln n) on L 2 (G n ). The analogous relation holds for l n on L 2 ( G n ). It is important that all these operators have separated variables. P. Exner: Spectral transitions... 4th Najman Conference September 22,

49 Proof outline continued This yields ( ) ln n y 2 V (xy) ln 2 n V (x ln n) = O n for for any (x, y) G n which allows us to check that ( ) ln n inf σ(h n ) inf σ(l n ) + O, n where l n := 2 x 2 2 y 2 + ω 2 ln 2 n ln 2 n V (x ln n) on L 2 (G n ). The analogous relation holds for l n on L 2 ( G n ). It is important that all these operators have separated variables. Since the minimal eigenvalue of Neumann Laplacian d2 on the strips dy 2 ln n < y ln(n + 1), n = 1, 2,..., is zero, we have inf σ(l n ) = inf σ ( l n (1) ), where the last operator on L 2 (R) acts as l (1) n = d2 dx 2 + ω2 ln 2 n ln 2 n V (x ln n) P. Exner: Spectral transitions... 4th Najman Conference September 22,

50 Proof outline concluded Note that the cut-off function χ { x a} plays no role in the asymptotic estimate as it affects a finite number of terms only. P. Exner: Spectral transitions... 4th Najman Conference September 22,

51 Proof outline concluded Note that the cut-off function χ { x a} plays no role in the asymptotic estimate as it affects a finite number of terms only. t By the change of variable x = ln n the last operator is unitarily equivalent to ln 2 n L which is non-negative as long as L is non-negative. In the same way one proves that l n is non-negative ; this concludes the proof. P. Exner: Spectral transitions... 4th Najman Conference September 22,

52 Proof outline concluded Note that the cut-off function χ { x a} plays no role in the asymptotic estimate as it affects a finite number of terms only. t By the change of variable x = ln n the last operator is unitarily equivalent to ln 2 n L which is non-negative as long as L is non-negative. In the same way one proves that l n is non-negative ; this concludes the proof. In the same way one can treat systems restricted in the x direction: Corollary Let H be our operator on ( c, c) R for some c a with Dirichlet (Neumann, periodic) boundary conditions in the variable x. The spectrum of H is bounded from below if L 0 holds, where L is the comparison operator on L 2 ( c, c) with Dirichlet (respectively, Neumann or periodic) boundary conditions. P. Exner: Spectral transitions... 4th Najman Conference September 22,

53 Supercritical case Once the transverse channel principal eigenvalue dominates over the harmonic oscillator contribution, the spectral behavior changes: P. Exner: Spectral transitions... 4th Najman Conference September 22,

54 Supercritical case Once the transverse channel principal eigenvalue dominates over the harmonic oscillator contribution, the spectral behavior changes: Theorem (Barseghyan-E 14) Under our hypotheses, σ(h) = R holds if inf σ(l) < 0. P. Exner: Spectral transitions... 4th Najman Conference September 22,

55 Supercritical case Once the transverse channel principal eigenvalue dominates over the harmonic oscillator contribution, the spectral behavior changes: Theorem (Barseghyan-E 14) Under our hypotheses, σ(h) = R holds if inf σ(l) < 0. Proof relies on construction of an appropriate Weyl sequence: we have to find {ψ k } k=1 D(H) such that ψ k = 1 which contains no convergent subsequence, and at same time Hψ k µψ k 0 as k. P. Exner: Spectral transitions... 4th Najman Conference September 22,

56 Supercritical case Once the transverse channel principal eigenvalue dominates over the harmonic oscillator contribution, the spectral behavior changes: Theorem (Barseghyan-E 14) Under our hypotheses, σ(h) = R holds if inf σ(l) < 0. Proof relies on construction of an appropriate Weyl sequence: we have to find {ψ k } k=1 D(H) such that ψ k = 1 which contains no convergent subsequence, and at same time Hψ k µψ k 0 as k. The construction is rather technical and we sketch just the main steps. P. Exner: Spectral transitions... 4th Najman Conference September 22,

57 Proof outline The claim is invariant under scaling transformations, hence we may suppose inf σ(l) = 1. The spectral threshold is a simple isolated eigenvalue; we denote the corresponding normalized eigenfunction by h. P. Exner: Spectral transitions... 4th Najman Conference September 22,

58 Proof outline The claim is invariant under scaling transformations, hence we may suppose inf σ(l) = 1. The spectral threshold is a simple isolated eigenvalue; we denote the corresponding normalized eigenfunction by h. We want to show first that 0 σ ess (H). In fact, it would be enough for the proof to show that zero belongs to σ(h) but we get the stronger claim at no extra expense. P. Exner: Spectral transitions... 4th Najman Conference September 22,

59 Proof outline The claim is invariant under scaling transformations, hence we may suppose inf σ(l) = 1. The spectral threshold is a simple isolated eigenvalue; we denote the corresponding normalized eigenfunction by h. We want to show first that 0 σ ess (H). In fact, it would be enough for the proof to show that zero belongs to σ(h) but we get the stronger claim at no extra expense. We fix an ε > 0 and choose a natural k = k(ε) with which we associate a function χ k C0 2 (1, k) satisfying the following conditions k 1 1 z χ2 k (z) dz = 1 and k 1 z(χ k (z))2 dz < ε. P. Exner: Spectral transitions... 4th Najman Conference September 22,

60 Proof outline continued Such functions exist: as an example consider χ k (z) = 8 ln3 z ln 3 k χ {1 z 2 ln k 2 ln z k} (z) + ln k χ { k+1 z k 1} (z) +g k (z)χ { k<z< k+1} (z) + q k (z)χ {k 1<z k} (z), where g k and q k are interpolating functions chosen in such a way that χ k C0 2 (1, k), and define ( k ) 1/2 1 χ k (z) = z χ2 k (z) dz χ k (z). 1 P. Exner: Spectral transitions... 4th Najman Conference September 22,

61 Proof outline continued Such functions exist: as an example consider χ k (z) = 8 ln3 z ln 3 k χ {1 z 2 ln k 2 ln z k} (z) + ln k χ { k+1 z k 1} (z) +g k (z)χ { k<z< k+1} (z) + q k (z)χ {k 1<z k} (z), where g k and q k are interpolating functions chosen in such a way that χ k C0 2 (1, k), and define ( k ) 1/2 1 χ k (z) = z χ2 k (z) dz χ k (z). Given such functions χ k, put 1 ( ) ψ k (x, y) := h(xy) e iy 2 y /2 χ k + f (xy) n k y 2 e iy 2 /2 χ k ( y where f (t) := i 2 t2 h(t), t R, and n k N is a positive integer, which we choose using the following auxiliary result. P. Exner: Spectral transitions... 4th Najman Conference September 22, n k ),

62 Proof outline continued Lemma Let ψ k, k = 1, 2,..., as defined above; then for any given k one can achieve that ψ k L 2 (R 2 ) 1 2 holds by choosing n k large enough. P. Exner: Spectral transitions... 4th Najman Conference September 22,

63 Proof outline continued Lemma Let ψ k, k = 1, 2,..., as defined above; then for any given k one can achieve that ψ k L 2 (R 2 ) 1 2 holds by choosing n k large enough. We need one more auxiliary result: Lemma Let ψ k, k = 1, 2,..., be again functions defined above; then the inequality Hψ k 2 L 2 (R 2 ) < cε with a fixed constant c holds for k = k(ε). P. Exner: Spectral transitions... 4th Najman Conference September 22,

64 Proof outline continued Lemma Let ψ k, k = 1, 2,..., as defined above; then for any given k one can achieve that ψ k L 2 (R 2 ) 1 2 holds by choosing n k large enough. We need one more auxiliary result: Lemma Let ψ k, k = 1, 2,..., be again functions defined above; then the inequality Hψ k 2 L 2 (R 2 ) < cε with a fixed constant c holds for k = k(ε). Proofs are in both cases straightforward but rather tedious. P. Exner: Spectral transitions... 4th Najman Conference September 22,

65 Proof outline concluded Using the lemmata, we are able to complete the proof. We fix a sequence {ε j } j=1 such that ε j 0 holds as j and to any j we construct a function ψ k(εj ) in such a way that n k(εj ) > k(ε j 1 )n k(εj 1 ). P. Exner: Spectral transitions... 4th Najman Conference September 22,

66 Proof outline concluded Using the lemmata, we are able to complete the proof. We fix a sequence {ε j } j=1 such that ε j 0 holds as j and to any j we construct a function ψ k(εj ) in such a way that n k(εj ) > k(ε j 1 )n k(εj 1 ). The norms of Hψ k(εj ) are bounded from above with 9ε j on the right-hand side, and since the supports of ψ k(εj ), j = 1, 2,..., do not intersect each other by construction, their sequence converges weakly to zero. P. Exner: Spectral transitions... 4th Najman Conference September 22,

67 Proof outline concluded Using the lemmata, we are able to complete the proof. We fix a sequence {ε j } j=1 such that ε j 0 holds as j and to any j we construct a function ψ k(εj ) in such a way that n k(εj ) > k(ε j 1 )n k(εj 1 ). The norms of Hψ k(εj ) are bounded from above with 9ε j on the right-hand side, and since the supports of ψ k(εj ), j = 1, 2,..., do not intersect each other by construction, their sequence converges weakly to zero. This yields the sought Weyl sequence for zero energy; for any nonzero real number µ we use the same procedure replacing the above ψ k with ( ) y ψ k (x, y) = h(xy) e iɛµ(y) χ k + f (xy) ( ) y n k y 2 e iɛµ(y) χ k, n k where ɛ µ (y) := y µ t 2 + µ dt, and furthermore, the functions f, χ k are defined in the same way as above. P. Exner: Spectral transitions... 4th Najman Conference September 22,

68 Restricted motion In the supercritical case, too, the result extends to systems restricted in the x direction: Theorem Let H be the our operator on L 2 ( c, c) L 2 (R) for some c > 0 with Dirichlet condition at x = ±c and denote by L the corresponding Dirichlet operator on L 2 ( c, c). If the spectral threshold of L is negative, the spectrum of H covers the whole real axis. P. Exner: Spectral transitions... 4th Najman Conference September 22,

69 Restricted motion In the supercritical case, too, the result extends to systems restricted in the x direction: Theorem Let H be the our operator on L 2 ( c, c) L 2 (R) for some c > 0 with Dirichlet condition at x = ±c and denote by L the corresponding Dirichlet operator on L 2 ( c, c). If the spectral threshold of L is negative, the spectrum of H covers the whole real axis. Observing the domains of the quadratic form associated with such operators we get Corollary The claim of the above theorem remains valid if the Dirichlet boundary conditions at x = ±c are replaced by any other self-adjoint boundary conditions. P. Exner: Spectral transitions... 4th Najman Conference September 22,

70 The multichannel case The above results are interesting not only per se or to deal with the Guarneri-type periodic modification of the model. P. Exner: Spectral transitions... 4th Najman Conference September 22,

71 The multichannel case The above results are interesting not only per se or to deal with the Guarneri-type periodic modification of the model. Using a simple bracketing argument we can show how the spectral-regime transition looks like in the multichannel case: Theorem (Barseghyan-E 14) Let H be our operator with the potentials satisfying the stated assumptions, namely the functions V j are positive with bounded first derivative and supp V j supp V k = holds for j k. Denote by L j the comparison operator on L 2 (R) with the potential V j and set t V := min j inf σ(l j ). Then H is bounded from below if and only if t V 0 and in the opposite case its spectrum covers the whole real axis. P. Exner: Spectral transitions... 4th Najman Conference September 22,

72 Another model class Consider next a related family of systems in which the transition is even more dramatic passing from purely discrete spectrum in the subcritical case to the whole real line in the supercritical one. P. Exner: Spectral transitions... 4th Najman Conference September 22,

73 Another model class Consider next a related family of systems in which the transition is even more dramatic passing from purely discrete spectrum in the subcritical case to the whole real line in the supercritical one. Recall that there are situations where Weyl s law fails and the spectrum is discrete even if the classically allowed phase-space volume is infinite. A classical example due to [Simon 83] is a 2D Schrödinger operator with the potential V (x, y) = x 2 y 2 or more generally, V (x, y) = xy p with p 1. P. Exner: Spectral transitions... 4th Najman Conference September 22,

74 Another model class Consider next a related family of systems in which the transition is even more dramatic passing from purely discrete spectrum in the subcritical case to the whole real line in the supercritical one. Recall that there are situations where Weyl s law fails and the spectrum is discrete even if the classically allowed phase-space volume is infinite. A classical example due to [Simon 83] is a 2D Schrödinger operator with the potential V (x, y) = x 2 y 2 or more generally, V (x, y) = xy p with p 1. Similar behavior one can observe for Dirichlet Laplacians in regions with hyperbolic cusps see [Geisinger-Weidl 11] for recent results and a survey. Moreover, using the dimensional-reduction technique of Laptev and Weidl one can prove spectral estimates for such operators. P. Exner: Spectral transitions... 4th Najman Conference September 22,

75 Another model class Consider next a related family of systems in which the transition is even more dramatic passing from purely discrete spectrum in the subcritical case to the whole real line in the supercritical one. Recall that there are situations where Weyl s law fails and the spectrum is discrete even if the classically allowed phase-space volume is infinite. A classical example due to [Simon 83] is a 2D Schrödinger operator with the potential V (x, y) = x 2 y 2 or more generally, V (x, y) = xy p with p 1. Similar behavior one can observe for Dirichlet Laplacians in regions with hyperbolic cusps see [Geisinger-Weidl 11] for recent results and a survey. Moreover, using the dimensional-reduction technique of Laptev and Weidl one can prove spectral estimates for such operators. A common feature of these models is that the particle motion is confined into channels narrowing towards infinity. P. Exner: Spectral transitions... 4th Najman Conference September 22,

76 Adding potentials unbounded from below This may remain true even for Schrödinger operators with unbounded from below in which a classical particle can escape to infinity with an increasing velocity. P. Exner: Spectral transitions... 4th Najman Conference September 22,

77 Adding potentials unbounded from below This may remain true even for Schrödinger operators with unbounded from below in which a classical particle can escape to infinity with an increasing velocity. The situation changes, however, if the attraction is strong enough P. Exner: Spectral transitions... 4th Najman Conference September 22,

78 Adding potentials unbounded from below This may remain true even for Schrödinger operators with unbounded from below in which a classical particle can escape to infinity with an increasing velocity. The situation changes, however, if the attraction is strong enough As an illustration, let us analyze the following class of operators: ( L p (λ) : L p (λ)ψ = ψ + xy p λ(x 2 + y 2 ) p/(p+2)) ψ, p 1 on L 2 (R 2 ), where (x, y) are the standard Cartesian coordinates in R 2 and the parameter λ in the second term of the potential is non-negative; unless the value of λ is important we write it simply as L p. P. Exner: Spectral transitions... 4th Najman Conference September 22,

79 Adding potentials unbounded from below This may remain true even for Schrödinger operators with unbounded from below in which a classical particle can escape to infinity with an increasing velocity. The situation changes, however, if the attraction is strong enough As an illustration, let us analyze the following class of operators: ( L p (λ) : L p (λ)ψ = ψ + xy p λ(x 2 + y 2 ) p/(p+2)) ψ, p 1 on L 2 (R 2 ), where (x, y) are the standard Cartesian coordinates in R 2 and the parameter λ in the second term of the potential is non-negative; unless the value of λ is important we write it simply as L p. 2p Note that p+2 < 2 so the operator is e.s.a. on C 0 (R2 ) by Faris-Lavine theorem again; the symbol L p or L p (λ) will always mean its closure. P. Exner: Spectral transitions... 4th Najman Conference September 22,

80 The subcritical case The spectral properties of L p (λ) depend crucially on the value of λ and there is a transition between different regimes as λ changes. P. Exner: Spectral transitions... 4th Najman Conference September 22,

81 The subcritical case The spectral properties of L p (λ) depend crucially on the value of λ and there is a transition between different regimes as λ changes. Let us start with the subcritical case which occurs for small values of λ. To characterize the smallness quantitatively we need an auxiliary operator which will be an (an)harmonic oscillator Hamiltonian on line, H p : H p u = u + t p u on L 2 (R) with the standard domain. Let γ p be the minimal eigenvalue of this operator; in view of the potential symmetry we have γ p = inf σ(h p ), where H p : H p u = u + t p u on L 2 (R + ) with Neumann condition at t = 0. P. Exner: Spectral transitions... 4th Najman Conference September 22,

82 The subcritical case continued The eigenvalue γ p = inf σ(h p ) equals one for p = 2; for p it becomes γ = 1 4 π2 ; it smoothly interpolates between the two values. P. Exner: Spectral transitions... 4th Najman Conference September 22,

83 The subcritical case continued The eigenvalue γ p = inf σ(h p ) equals one for p = 2; for p it becomes γ = 1 4 π2 ; it smoothly interpolates between the two values. Since x p 1 χ [0,1] (x) we have γ p ɛ , where ɛ 0 is the ground-state energy of the rectangular potential well of depth one. P. Exner: Spectral transitions... 4th Najman Conference September 22,

84 The subcritical case continued The eigenvalue γ p = inf σ(h p ) equals one for p = 2; for p it becomes γ = 1 4 π2 ; it smoothly interpolates between the two values. Since x p 1 χ [0,1] (x) we have γ p ɛ , where ɛ 0 is the ground-state energy of the rectangular potential well of depth one. In fact, a numerical solution gives true minimum γ p attained at p 1.788; in the semilogarithmic scale the plot is as follows: γ p p P. Exner: Spectral transitions... 4th Najman Conference September 22,

85 The subcritical case continued The spectrum is naturally bounded from below and discrete if λ = 0; our aim is to show that this remains to be the case provided λ is small enough. P. Exner: Spectral transitions... 4th Najman Conference September 22,

86 The subcritical case continued The spectrum is naturally bounded from below and discrete if λ = 0; our aim is to show that this remains to be the case provided λ is small enough. Theorem (E-Barseghyan 12) For any λ [0, λ crit ], where λ crit := γ p, the operator L p (λ) is bounded from below for p 1; if λ < γ p its spectrum is purely discrete. P. Exner: Spectral transitions... 4th Najman Conference September 22,

87 The subcritical case continued The spectrum is naturally bounded from below and discrete if λ = 0; our aim is to show that this remains to be the case provided λ is small enough. Theorem (E-Barseghyan 12) For any λ [0, λ crit ], where λ crit := γ p, the operator L p (λ) is bounded from below for p 1; if λ < γ p its spectrum is purely discrete. Idea of the proof: Let λ < γ p. By minimax we need to estimate L p from below by a s-a operator with a purely discrete spectrum. To construct it we employ bracketing imposing additional Neumann conditions at concentric circles of radii n = 1, 2,.... P. Exner: Spectral transitions... 4th Najman Conference September 22,

88 The subcritical case continued The spectrum is naturally bounded from below and discrete if λ = 0; our aim is to show that this remains to be the case provided λ is small enough. Theorem (E-Barseghyan 12) For any λ [0, λ crit ], where λ crit := γ p, the operator L p (λ) is bounded from below for p 1; if λ < γ p its spectrum is purely discrete. Idea of the proof: Let λ < γ p. By minimax we need to estimate L p from below by a s-a operator with a purely discrete spectrum. To construct it we employ bracketing imposing additional Neumann conditions at concentric circles of radii n = 1, 2,.... In the estimating operators the variables decouple asymptotically and the spectral behavior is determined by the angular part of the operators. P. Exner: Spectral transitions... 4th Najman Conference September 22,

89 Subcritical behavior the proof Specifically, in polar coordinates we get direct sum of operators acting as L (1) n,pψ = 1 ( r ψ ) 1 2 ( ) ψ r 2p r r r n 2 ϕ p sin 2ϕ p λr 2p/(p+2) ψ on the annuli G n := {(r, ϕ) : n 1 r < n, 0 ϕ < 2π}, n = 1, 2,... with Neumann conditions imposed on G n. P. Exner: Spectral transitions... 4th Najman Conference September 22,

90 Subcritical behavior the proof Specifically, in polar coordinates we get direct sum of operators acting as L (1) n,pψ = 1 ( r ψ ) 1 2 ( ) ψ r 2p r r r n 2 ϕ p sin 2ϕ p λr 2p/(p+2) ψ on the annuli G n := {(r, ϕ) : n 1 r < n, 0 ϕ < 2π}, n = 1, 2,... with Neumann conditions imposed on G n. Obviously σ(l (1) n,p) is purely discrete for each n = 1, 2,..., hence it is sufficient to check that inf σ(l (1) n,p) holds as n. P. Exner: Spectral transitions... 4th Najman Conference September 22,

91 Subcritical behavior the proof Specifically, in polar coordinates we get direct sum of operators acting as L (1) n,pψ = 1 ( r ψ ) 1 2 ( ) ψ r 2p r r r n 2 ϕ p sin 2ϕ p λr 2p/(p+2) ψ on the annuli G n := {(r, ϕ) : n 1 r < n, 0 ϕ < 2π}, n = 1, 2,... with Neumann conditions imposed on G n. Obviously σ(l (1) n,p) is purely discrete for each n = 1, 2,..., hence it is sufficient to check that inf σ(l (1) n,p) holds as n. We estimate L (1) n,p from below by an operator with separating variables, note that the radial part does not contribute and use the symmetry of the problem; for ε (0, 1) the question is then to analyze ( L (2) n,p : L (2) n n,pu = u 2p p sin p 2x λ ) 1 ε n(4p+4)/(p+2) u on L 2 (0, π/4) with Neumann conditions, u (0) = u (π/4) = 0. P. Exner: Spectral transitions... 4th Najman Conference September 22,

92 Subcritical behavior proof continued We have n 2 inf σ(l (1) n,p) inf σ(l (2) n 1,p ) if n is large enough, specifically for n > ( 1 (1 ε) (p+2)/(4p+4)) 1, hence it is sufficient to investigate the spectral threshold µ n,p of L (2) n,p as n. P. Exner: Spectral transitions... 4th Najman Conference September 22,

93 Subcritical behavior proof continued We have n 2 inf σ(l (1) n,p) inf σ(l (2) n 1,p ) if n is large enough, specifically for n > ( 1 (1 ε) (p+2)/(4p+4)) 1, hence it is sufficient to investigate the spectral threshold µ n,p of L (2) n,p as n. The trigonometric potential can be estimated by a powerlike one with the similar behavior around the minimum introducing, e.g. ( ( ) p L (3) n,p := d2 2 dx 2 + n2p+2 x p χ (0,δ(ε)] (x) + χ [δ(ε),π/4)(x)) λ ε n (4p+4)/(p+2) π for small enough δ(ε) with Neumann boundary conditions at x = 0, 1 4 π, where we have denoted λ ε := λ(1 ε) p 1. P. Exner: Spectral transitions... 4th Najman Conference September 22,

94 Subcritical behavior proof continued We have n 2 inf σ(l (1) n,p) inf σ(l (2) n 1,p ) if n is large enough, specifically for n > ( 1 (1 ε) (p+2)/(4p+4)) 1, hence it is sufficient to investigate the spectral threshold µ n,p of L (2) n,p as n. The trigonometric potential can be estimated by a powerlike one with the similar behavior around the minimum introducing, e.g. ( ( ) p L (3) n,p := d2 2 dx 2 + n2p+2 x p χ (0,δ(ε)] (x) + χ [δ(ε),π/4)(x)) λ ε n (4p+4)/(p+2) π for small enough δ(ε) with Neumann boundary conditions at x = 0, 1 4 π, where we have denoted λ ε := λ(1 ε) p 1. We have L (2) n,p (1 ε) p L (3) n,p. To estimate the rhs by comparing the indicated potential contributions it is useful to pass to the rescaled variable x = t n (2p+2)/(p+2). P. Exner: Spectral transitions... 4th Najman Conference September 22,

95 Subcritical behavior proof concluded In this way we find that µ n,p := inf σ(l (3) n,p) satisfies µ n,p n 2 as n P. Exner: Spectral transitions... 4th Najman Conference September 22,

96 Subcritical behavior proof concluded In this way we find that µ n,p := inf σ(l (3) n,p) satisfies µ n,p n 2 as n Through the chain of inequalities we come to conclusion that inf σ(l (1) n,p) holds as n which proves discreteness of the spectrum for λ < γ p. P. Exner: Spectral transitions... 4th Najman Conference September 22,

97 Subcritical behavior proof concluded In this way we find that µ n,p := inf σ(l (3) n,p) satisfies µ n,p n 2 as n Through the chain of inequalities we come to conclusion that inf σ(l (1) n,p) holds as n which proves discreteness of the spectrum for λ < γ p. If λ = γ p the sequence of spectral thresholds no longer diverges but it remains bounded from below and the same is by minimax principle true for the operator L p (λ). P. Exner: Spectral transitions... 4th Najman Conference September 22,

98 Subcritical behavior proof concluded In this way we find that µ n,p := inf σ(l (3) n,p) satisfies µ n,p n 2 as n Through the chain of inequalities we come to conclusion that inf σ(l (1) n,p) holds as n which proves discreteness of the spectrum for λ < γ p. If λ = γ p the sequence of spectral thresholds no longer diverges but it remains bounded from below and the same is by minimax principle true for the operator L p (λ). Remark It is natural to conjecture that σ(l p (γ p )) R +. There may be a negative discrete spectrum in the critical case; we return to this question a little later. P. Exner: Spectral transitions... 4th Najman Conference September 22,

99 The supercritical case Theorem (E-Barseghyan 12) The spectrum of L p (λ), p 1, is unbounded below from if λ > λ crit. P. Exner: Spectral transitions... 4th Najman Conference September 22,

100 The supercritical case Theorem (E-Barseghyan 12) The spectrum of L p (λ), p 1, is unbounded below from if λ > λ crit. Idea of the proof: Similar as above with a few differences: P. Exner: Spectral transitions... 4th Najman Conference September 22,

101 The supercritical case Theorem (E-Barseghyan 12) The spectrum of L p (λ), p 1, is unbounded below from if λ > λ crit. Idea of the proof: Similar as above with a few differences: now we seek an upper bound to L p (λ) by a below unbounded operator, hence we impose Dirichlet conditions on concentric circles P. Exner: Spectral transitions... 4th Najman Conference September 22,

102 The supercritical case Theorem (E-Barseghyan 12) The spectrum of L p (λ), p 1, is unbounded below from if λ > λ crit. Idea of the proof: Similar as above with a few differences: now we seek an upper bound to L p (λ) by a below unbounded operator, hence we impose Dirichlet conditions on concentric circles the estimating operators have now a nonzero contribution from the radial part, however, it is bounded by π 2 independently of n P. Exner: Spectral transitions... 4th Najman Conference September 22,

103 The supercritical case Theorem (E-Barseghyan 12) The spectrum of L p (λ), p 1, is unbounded below from if λ > λ crit. Idea of the proof: Similar as above with a few differences: now we seek an upper bound to L p (λ) by a below unbounded operator, hence we impose Dirichlet conditions on concentric circles the estimating operators have now a nonzero contribution from the radial part, however, it is bounded by π 2 independently of n the negative λ-dependent term now outweights the anharmonic oscillator part so that inf σ(l (1,D) n,p ) holds as n P. Exner: Spectral transitions... 4th Najman Conference September 22,

104 The supercritical case Theorem (E-Barseghyan 12) The spectrum of L p (λ), p 1, is unbounded below from if λ > λ crit. Idea of the proof: Similar as above with a few differences: now we seek an upper bound to L p (λ) by a below unbounded operator, hence we impose Dirichlet conditions on concentric circles the estimating operators have now a nonzero contribution from the radial part, however, it is bounded by π 2 independently of n the negative λ-dependent term now outweights the anharmonic oscillator part so that inf σ(l (1,D) n,p ) holds as n Using suitable Weyl sequences similar to those the previous model, however, we are able to get a stronger result: Theorem (Barseghyan-E-Khrabustovskyi-Tater 15) σ(l p (λ)) = R holds for any λ > γ p and p > 1. P. Exner: Spectral transitions... 4th Najman Conference September 22,

105 Spectral estimates: bounds to eigenvalue sums Let us return to the subcritical case and define the following quantity: α := 1 ( 1 + ) > γ 1 p 2 P. Exner: Spectral transitions... 4th Najman Conference September 22,

Schrödinger operators exhibiting a sudden change of the spectral character

Schrödinger operators exhibiting a sudden change of the spectral character Schrödinger operators exhibiting a sudden change of the spectral character Pavel Exner Doppler Institute for Mathematical Physics and Applied Mathematics Prague in collaboration with Diana Barseghyan,

More information

Quantum systems exhibiting parameter-dependent spectral transitions

Quantum systems exhibiting parameter-dependent spectral transitions Quantum systems exhibiting parameter-dependent spectral transitions Pavel Exner Doppler Institute for Mathematical Physics and Applied Mathematics Prague in collaboration with Diana Barseghyan, Andrii

More information

Schrödinger operators exhibiting a sudden change of the spectral character

Schrödinger operators exhibiting a sudden change of the spectral character Schrödinger operators exhibiting a sudden change of the spectral character Pavel Exner Doppler Institute for Mathematical Physics and Applied Mathematics Prague in collaboration with Diana Barseghyan,

More information

Smilansky-Solomyak model with a δ -interaction

Smilansky-Solomyak model with a δ -interaction Smilansky-Solomyak model with a δ -interaction Jiří Lipovský University of Hradec Králové, Faculty of Science jiri.lipovsky@uhk.cz joint work with P. Exner Ostrava, April 6, 2018 Jiří Lipovský Smilansky-Solomyak

More information

On the ground state of quantum graphs with attractive δ-coupling

On the ground state of quantum graphs with attractive δ-coupling On the ground state of quantum graphs with attractive δ-coupling Pavel Exner a,b,, Michal ex a,c a Doppler Institute for Mathematical Physics and Applied Mathematics, Czech Technical University, Břehová

More information

ON THE REMOVAL OF FINITE DISCRETE SPECTRUM BY COEFFICIENT STRIPPING

ON THE REMOVAL OF FINITE DISCRETE SPECTRUM BY COEFFICIENT STRIPPING ON THE REMOVAL OF FINITE DISCRETE SPECTRUM BY COEFFICIENT STRIPPING BARRY SIMON Abstract. We prove for a large class of operators, J, including block Jacobi matrices, if σ(j) \ [α, β] is a finite set,

More information

arxiv: v1 [math-ph] 23 Mar 2014

arxiv: v1 [math-ph] 23 Mar 2014 March 25, 2014 2:57 WSPC Proceedings - 9.75in x 6.5in delta line arxiv page 1 1 Spectral asymptotics for δ interaction supported by a infinite curve arxiv:1403.5798v1 [math-ph] 23 Mar 2014 Michal Jex Doppler

More information

ON THE SPECTRUM OF NARROW NEUMANN WAVEGUIDE WITH PERIODICALLY DISTRIBUTED δ TRAPS

ON THE SPECTRUM OF NARROW NEUMANN WAVEGUIDE WITH PERIODICALLY DISTRIBUTED δ TRAPS ON THE SPECTRUM OF NARROW NEUMANN WAVEGUIDE WITH PERIODICALLY DISTRIBUTED δ TRAPS PAVEL EXNER 1 AND ANDRII KHRABUSTOVSKYI 2 Abstract. We analyze a family of singular Schrödinger operators describing a

More information

Asymptotic behaviour of the heat equation in twisted waveguides

Asymptotic behaviour of the heat equation in twisted waveguides Asymptotic behaviour of the heat equation in twisted waveguides Gabriela Malenová Faculty of Nuclear Sciences and Physical Engineering, CTU, Prague Nuclear Physics Institute, AS ČR, Řež Graphs and Spectra,

More information

SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS

SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties

More information

Straight quantum waveguide with Robin boundary conditions

Straight quantum waveguide with Robin boundary conditions Straight quantum waveguide with Robin boundary conditions by Martin JÍLEK Supervisor: David KREJČIŘÍK 2 Consultant: Miloš Tater 2 Abstract We investigate spectral properties of a quantum particle confined

More information

BOUNDARY VALUE PROBLEMS ON A HALF SIERPINSKI GASKET

BOUNDARY VALUE PROBLEMS ON A HALF SIERPINSKI GASKET BOUNDARY VALUE PROBLEMS ON A HALF SIERPINSKI GASKET WEILIN LI AND ROBERT S. STRICHARTZ Abstract. We study boundary value problems for the Laplacian on a domain Ω consisting of the left half of the Sierpinski

More information

Stability of the Magnetic Schrödinger Operator in a Waveguide

Stability of the Magnetic Schrödinger Operator in a Waveguide Communications in Partial Differential Equations, 3: 539 565, 5 Copyright Taylor & Francis, Inc. ISSN 36-53 print/153-4133 online DOI: 1.181/PDE-5113 Stability of the Magnetic Schrödinger Operator in a

More information

Page 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19

Page 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 information

Essential Spectra of complete manifolds

Essential Spectra of complete manifolds Essential Spectra of complete manifolds Zhiqin Lu Analysis, Complex Geometry, and Mathematical Physics: A Conference in Honor of Duong H. Phong May 7, 2013 Zhiqin Lu, Dept. Math, UCI Essential Spectra

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,

More information

Schrödinger operators with singular interactions on conical surfaces

Schrödinger operators with singular interactions on conical surfaces Schrödinger operators with singular interactions on conical surfaces Vladimir Lotoreichik Nuclear Physics Institute CAS, Řež, Czech Repbulic Bilbao, 14.05.2015 Vladimir Lotoreichik (NPI CAS) 14.05.2015

More information

PHY 407 QUANTUM MECHANICS Fall 05 Problem set 1 Due Sep

PHY 407 QUANTUM MECHANICS Fall 05 Problem set 1 Due Sep Problem set 1 Due Sep 15 2005 1. Let V be the set of all complex valued functions of a real variable θ, that are periodic with period 2π. That is u(θ + 2π) = u(θ), for all u V. (1) (i) Show that this V

More information

Asymptotic distribution of eigenvalues of Laplace operator

Asymptotic distribution of eigenvalues of Laplace operator Asymptotic distribution of eigenvalues of Laplace operator 23.8.2013 Topics We will talk about: the number of eigenvalues of Laplace operator smaller than some λ as a function of λ asymptotic behaviour

More information

UNIQUENESS OF POSITIVE SOLUTION TO SOME COUPLED COOPERATIVE VARIATIONAL ELLIPTIC SYSTEMS

UNIQUENESS OF POSITIVE SOLUTION TO SOME COUPLED COOPERATIVE VARIATIONAL ELLIPTIC SYSTEMS TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9947(XX)0000-0 UNIQUENESS OF POSITIVE SOLUTION TO SOME COUPLED COOPERATIVE VARIATIONAL ELLIPTIC SYSTEMS YULIAN

More information

Waves on 2 and 3 dimensional domains

Waves on 2 and 3 dimensional domains Chapter 14 Waves on 2 and 3 dimensional domains We now turn to the studying the initial boundary value problem for the wave equation in two and three dimensions. In this chapter we focus on the situation

More information

Singular Schrödinger operators with interactions supported by sets of codimension one

Singular Schrödinger operators with interactions supported by sets of codimension one Singular Schrödinger operators with interactions supported by sets of codimension one Pavel Exner Doppler Institute for Mathematical Physics and Applied Mathematics Prague in collaboration with Jussi Behrndt,

More information

Eigenvalues of Robin Laplacians on infinite sectors and application to polygons

Eigenvalues of Robin Laplacians on infinite sectors and application to polygons Eigenvalues of Robin Laplacians on infinite sectors and application to polygons Magda Khalile (joint work with Konstantin Pankrashkin) Université Paris Sud /25 Robin Laplacians on infinite sectors 1 /

More information

Partial Differential Equations

Partial Differential Equations Partial Differential Equations Spring Exam 3 Review Solutions Exercise. We utilize the general solution to the Dirichlet problem in rectangle given in the textbook on page 68. In the notation used there

More information

The Sommerfeld Polynomial Method: Harmonic Oscillator Example

The Sommerfeld Polynomial Method: Harmonic Oscillator Example Chemistry 460 Fall 2017 Dr. Jean M. Standard October 2, 2017 The Sommerfeld Polynomial Method: Harmonic Oscillator Example Scaling the Harmonic Oscillator Equation Recall the basic definitions of the harmonic

More information

1 Lyapunov theory of stability

1 Lyapunov theory of stability M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability

More information

REVIEW OF DIFFERENTIAL CALCULUS

REVIEW OF DIFFERENTIAL CALCULUS REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be

More information

Periodic Schrödinger operators with δ -potentials

Periodic Schrödinger operators with δ -potentials Networ meeting April 23-28, TU Graz Periodic Schrödinger operators with δ -potentials Andrii Khrabustovsyi Institute for Analysis, Karlsruhe Institute of Technology, Germany CRC 1173 Wave phenomena: analysis

More information

Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.

Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip

More information

GEOMETRIC VERSUS SPECTRAL CONVERGENCE FOR THE NEUMANN LAPLACIAN UNDER EXTERIOR PERTURBATIONS OF THE DOMAIN

GEOMETRIC VERSUS SPECTRAL CONVERGENCE FOR THE NEUMANN LAPLACIAN UNDER EXTERIOR PERTURBATIONS OF THE DOMAIN GEOMETRIC VERSUS SPECTRAL CONVERGENCE FOR THE NEUMANN LAPLACIAN UNDER EXTERIOR PERTURBATIONS OF THE DOMAIN JOSÉ M. ARRIETA AND DAVID KREJČIŘÍK Abstract. We analyze the behavior of the eigenvalues and eigenfunctions

More information

C.7. Numerical series. Pag. 147 Proof of the converging criteria for series. Theorem 5.29 (Comparison test) Let a k and b k be positive-term series

C.7. Numerical series. Pag. 147 Proof of the converging criteria for series. Theorem 5.29 (Comparison test) Let a k and b k be positive-term series C.7 Numerical series Pag. 147 Proof of the converging criteria for series Theorem 5.29 (Comparison test) Let and be positive-term series such that 0, for any k 0. i) If the series converges, then also

More information

here, this space is in fact infinite-dimensional, so t σ ess. Exercise Let T B(H) be a self-adjoint operator on an infinitedimensional

here, this space is in fact infinite-dimensional, so t σ ess. Exercise Let T B(H) be a self-adjoint operator on an infinitedimensional 15. Perturbations by compact operators In this chapter, we study the stability (or lack thereof) of various spectral properties under small perturbations. Here s the type of situation we have in mind:

More information

THE ESSENTIAL SELF-ADJOINTNESS OF SCHRÖDINGER OPERATORS WITH OSCILLATING POTENTIALS. Adam D. Ward. In Memory of Boris Pavlov ( )

THE ESSENTIAL SELF-ADJOINTNESS OF SCHRÖDINGER OPERATORS WITH OSCILLATING POTENTIALS. Adam D. Ward. In Memory of Boris Pavlov ( ) NEW ZEALAND JOURNAL OF MATHEMATICS Volume 46 (016), 65-7 THE ESSENTIAL SELF-ADJOINTNESS OF SCHRÖDINGER OPERATORS WITH OSCILLATING POTENTIALS Adam D. Ward (Received 4 May, 016) In Memory of Boris Pavlov

More information

The Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:

The 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 information

MATH 31BH Homework 1 Solutions

MATH 31BH Homework 1 Solutions MATH 3BH Homework Solutions January 0, 04 Problem.5. (a) (x, y)-plane in R 3 is closed and not open. To see that this plane is not open, notice that any ball around the origin (0, 0, 0) will contain points

More information

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential

More information

An introduction to Birkhoff normal form

An 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 information

ON THE SPECTRUM OF THE DIRICHLET LAPLACIAN IN A NARROW STRIP

ON THE SPECTRUM OF THE DIRICHLET LAPLACIAN IN A NARROW STRIP ON THE SPECTRUM OF THE DIRICHLET LAPLACIAN IN A NARROW STRIP LEONID FRIEDLANDER AND MICHAEL SOLOMYAK Abstract. We consider the Dirichlet Laplacian in a family of bounded domains { a < x < b, 0 < y < h(x)}.

More information

Spotlight on Laplace s Equation

Spotlight on Laplace s Equation 16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian

More information

Uniqueness of ground state solutions of non-local equations in R N

Uniqueness of ground state solutions of non-local equations in R N Uniqueness of ground state solutions of non-local equations in R N Rupert L. Frank Department of Mathematics Princeton University Joint work with Enno Lenzmann and Luis Silvestre Uniqueness and non-degeneracy

More information

Quantum waveguides: mathematical problems

Quantum waveguides: mathematical problems Quantum waveguides: mathematical problems Pavel Exner in part in collaboration with Denis Borisov and Olaf Post exner@ujf.cas.cz Department of Theoretical Physics, NPI, Czech Academy of Sciences and Doppler

More information

Existence and Uniqueness

Existence and Uniqueness Chapter 3 Existence and Uniqueness An intellect which at a certain moment would know all forces that set nature in motion, and all positions of all items of which nature is composed, if this intellect

More information

Asymptotics of generalized eigenfunctions on manifold with Euclidean and/or hyperbolic ends

Asymptotics of generalized eigenfunctions on manifold with Euclidean and/or hyperbolic ends Asymptotics of generalized eigenfunctions on manifold with Euclidean and/or hyperbolic ends Kenichi ITO (University of Tokyo) joint work with Erik SKIBSTED (Aarhus University) 3 July 2018 Example: Free

More information

Aharonov and Bohm vs. Welsh eigenvalues

Aharonov and Bohm vs. Welsh eigenvalues Aharonov and Bohm vs. Welsh eigenvalues P. EXNER 1, and S. KONDEJ 3 1 Doppler Institute for Mathematical Physics and Applied Mathematics, Czech Technical University in Prague, Břehová 7, 11519 Prague,

More information

Equivariant self-similar wave maps from Minkowski spacetime into 3-sphere

Equivariant self-similar wave maps from Minkowski spacetime into 3-sphere Equivariant self-similar wave maps from Minkowski spacetime into 3-sphere arxiv:math-ph/99126v1 17 Oct 1999 Piotr Bizoń Institute of Physics, Jagellonian University, Kraków, Poland March 26, 28 Abstract

More information

One-dimensional Schrödinger equation

One-dimensional Schrödinger equation Chapter 1 One-dimensional Schrödinger equation In this chapter we will start from the harmonic oscillator to introduce a general numerical methodology to solve the one-dimensional, time-independent Schrödinger

More information

CHAPTER 1: Review (See also the Precalculus notes at

CHAPTER 1: Review (See also the Precalculus notes at CHAPTER 1: Review (See also the Precalculus notes at http://www.kkuniyuk.com) TOPIC 1: FUNCTIONS (Chapter 1: Review) 1.01 PART A: AN EXAMPLE OF A FUNCTION Consider a function f whose rule is given by f

More information

LEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.

LEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates. LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is

More information

arxiv: v1 [physics.comp-ph] 28 Jan 2008

arxiv: v1 [physics.comp-ph] 28 Jan 2008 A new method for studying the vibration of non-homogeneous membranes arxiv:0801.4369v1 [physics.comp-ph] 28 Jan 2008 Abstract Paolo Amore Facultad de Ciencias, Universidad de Colima, Bernal Díaz del Castillo

More information

Title Project Summary

Title Project Summary Title Project Summary The aim of the project is an estimation theory of those special functions of analysis called zeta functions after the zeta function of Euler (1730). The desired estimates generalize

More information

The following definition is fundamental.

The following definition is fundamental. 1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic

More information

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points Nonlinear Control Lecture # 3 Stability of Equilibrium Points The Invariance Principle Definitions Let x(t) be a solution of ẋ = f(x) A point p is a positive limit point of x(t) if there is a sequence

More information

A Minimal Uncertainty Product for One Dimensional Semiclassical Wave Packets

A Minimal Uncertainty Product for One Dimensional Semiclassical Wave Packets A Minimal Uncertainty Product for One Dimensional Semiclassical Wave Packets George A. Hagedorn Happy 60 th birthday, Mr. Fritz! Abstract. Although real, normalized Gaussian wave packets minimize the product

More information

Ordinary Differential Equations II

Ordinary Differential Equations II Ordinary Differential Equations II February 23 2017 Separation of variables Wave eq. (PDE) 2 u t (t, x) = 2 u 2 c2 (t, x), x2 c > 0 constant. Describes small vibrations in a homogeneous string. u(t, x)

More information

Eigenvalues and eigenfunctions of the Laplacian. Andrew Hassell

Eigenvalues and eigenfunctions of the Laplacian. Andrew Hassell Eigenvalues and eigenfunctions of the Laplacian Andrew Hassell 1 2 The setting In this talk I will consider the Laplace operator,, on various geometric spaces M. Here, M will be either a bounded Euclidean

More information

08a. Operators on Hilbert spaces. 1. Boundedness, continuity, operator norms

08a. Operators on Hilbert spaces. 1. Boundedness, continuity, operator norms (February 24, 2017) 08a. Operators on Hilbert spaces Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/real/notes 2016-17/08a-ops

More information

Ground-State Energies of Coulomb Systems and Reduced One-Particle Density Matrices

Ground-State Energies of Coulomb Systems and Reduced One-Particle Density Matrices Ground-State Energies of Coulomb Systems and Reduced One-Particle Density Matrices Heinz Siedentop Chennai, August 16, 2010 I. Introduction Atomic Schrödinger operator H N,Z := N n=1 self-adjointly realized

More information

A Variational Analysis of a Gauged Nonlinear Schrödinger Equation

A Variational Analysis of a Gauged Nonlinear Schrödinger Equation A Variational Analysis of a Gauged Nonlinear Schrödinger Equation Alessio Pomponio, joint work with David Ruiz Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Variational and Topological

More information

Gradient Estimates and Sobolev Inequality

Gradient Estimates and Sobolev Inequality Gradient Estimates and Sobolev Inequality Jiaping Wang University of Minnesota ( Joint work with Linfeng Zhou) Conference on Geometric Analysis in honor of Peter Li University of California, Irvine January

More information

Elements of linear algebra

Elements of linear algebra Elements of linear algebra Elements of linear algebra A vector space S is a set (numbers, vectors, functions) which has addition and scalar multiplication defined, so that the linear combination c 1 v

More information

Estimates on Neumann eigenfunctions at the boundary, and the Method of Particular Solutions" for computing them

Estimates on Neumann eigenfunctions at the boundary, and the Method of Particular Solutions for computing them Estimates on Neumann eigenfunctions at the boundary, and the Method of Particular Solutions" for computing them Department of Mathematics Australian National University Dartmouth, July 2010 Outline Introduction

More information

Ari Laptev and Timo Weidl. Department of Mathematics Royal Institute of Technology SE Stockholm Sweden

Ari Laptev and Timo Weidl. Department of Mathematics Royal Institute of Technology SE Stockholm Sweden RECENT RESULTS ON LIEB-THIRRING INEQUALITIES Ari Laptev and Timo Weidl Department of Mathematics Royal Institute of Technology SE-10044 Stockholm Sweden laptev@math.kth.se weidl@math.kth.se June 3, 2000

More information

Simple Examples on Rectangular Domains

Simple Examples on Rectangular Domains 84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation

More information

LAPLACE EQUATION. = 2 is call the Laplacian or del-square operator. In two dimensions, the expression of in rectangular and polar coordinates are

LAPLACE EQUATION. = 2 is call the Laplacian or del-square operator. In two dimensions, the expression of in rectangular and polar coordinates are LAPLACE EQUATION If a diffusion or wave problem is stationary (time independent), the pde reduces to the Laplace equation u = u =, an archetype of second order elliptic pde. = 2 is call the Laplacian or

More information

Diffusion on the half-line. The Dirichlet problem

Diffusion on the half-line. The Dirichlet problem Diffusion on the half-line The Dirichlet problem Consider the initial boundary value problem (IBVP) on the half line (, ): v t kv xx = v(x, ) = φ(x) v(, t) =. The solution will be obtained by the reflection

More information

Linear Algebra. Session 12

Linear Algebra. Session 12 Linear Algebra. Session 12 Dr. Marco A Roque Sol 08/01/2017 Example 12.1 Find the constant function that is the least squares fit to the following data x 0 1 2 3 f(x) 1 0 1 2 Solution c = 1 c = 0 f (x)

More information

Laplace s Equation. Chapter Mean Value Formulas

Laplace s Equation. Chapter Mean Value Formulas Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic

More information

Mathematical Methods for Physics and Engineering

Mathematical Methods for Physics and Engineering Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory

More information

WEYL S LAW FOR HYPERBOLIC RIEMANN SURFACES

WEYL S LAW FOR HYPERBOLIC RIEMANN SURFACES WEYL S LAW OR HYPERBOLIC RIEMANN SURACES MATTHEW STEVENSON Abstract. These are notes for a talk given in Dima Jakobson s class on automorphic forms at McGill University. This is a brief survey of the results

More information

SOLUTIONS TO SELECTED PROBLEMS FROM ASSIGNMENTS 3, 4

SOLUTIONS TO SELECTED PROBLEMS FROM ASSIGNMENTS 3, 4 SOLUTIONS TO SELECTED POBLEMS FOM ASSIGNMENTS 3, 4 Problem 5 from Assignment 3 Statement. Let be an n-dimensional bounded domain with smooth boundary. Show that the eigenvalues of the Laplacian on with

More information

We will outline two proofs of the main theorem:

We will outline two proofs of the main theorem: Chapter 4 Higher Genus Surfaces 4.1 The Main Result We will outline two proofs of the main theorem: Theorem 4.1.1. Let Σ be a closed oriented surface of genus g > 1. Then every homotopy class of homeomorphisms

More information

Stationary radial spots in a planar threecomponent reaction-diffusion system

Stationary radial spots in a planar threecomponent reaction-diffusion system Stationary radial spots in a planar threecomponent reaction-diffusion system Peter van Heijster http://www.dam.brown.edu/people/heijster SIAM Conference on Nonlinear Waves and Coherent Structures MS: Recent

More information

Entire functions, PT-symmetry and Voros s quantization scheme

Entire functions, PT-symmetry and Voros s quantization scheme Entire functions, PT-symmetry and Voros s quantization scheme Alexandre Eremenko January 19, 2017 Abstract In this paper, A. Avila s theorem on convergence of the exact quantization scheme of A. Voros

More information

REVIEW: The Matching Method Algorithm

REVIEW: The Matching Method Algorithm Lecture 26: Numerov Algorithm for Solving the Time-Independent Schrödinger Equation 1 REVIEW: The Matching Method Algorithm Need for a more general method The shooting method for solving the time-independent

More information

21 Laplace s Equation and Harmonic Functions

21 Laplace s Equation and Harmonic Functions 2 Laplace s Equation and Harmonic Functions 2. Introductory Remarks on the Laplacian operator Given a domain Ω R d, then 2 u = div(grad u) = in Ω () is Laplace s equation defined in Ω. If d = 2, in cartesian

More information

Daniel M. Oberlin Department of Mathematics, Florida State University. January 2005

Daniel M. Oberlin Department of Mathematics, Florida State University. January 2005 PACKING SPHERES AND FRACTAL STRICHARTZ ESTIMATES IN R d FOR d 3 Daniel M. Oberlin Department of Mathematics, Florida State University January 005 Fix a dimension d and for x R d and r > 0, let Sx, r) stand

More information

Spectral stability of periodic waves in dispersive models

Spectral stability of periodic waves in dispersive models in dispersive models Collaborators: Th. Gallay, E. Lombardi T. Kapitula, A. Scheel in dispersive models One-dimensional nonlinear waves Standing and travelling waves u(x ct) with c = 0 or c 0 periodic

More information

u( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)

u( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3) M ath 5 2 7 Fall 2 0 0 9 L ecture 4 ( S ep. 6, 2 0 0 9 ) Properties and Estimates of Laplace s and Poisson s Equations In our last lecture we derived the formulas for the solutions of Poisson s equation

More information

which implies that we can take solutions which are simultaneous eigen functions of

which implies that we can take solutions which are simultaneous eigen functions of Module 1 : Quantum Mechanics Chapter 6 : Quantum mechanics in 3-D Quantum mechanics in 3-D For most physical systems, the dynamics is in 3-D. The solutions to the general 3-d problem are quite complicated,

More information

Stable solitons of the cubic-quintic NLS with a delta-function potential

Stable solitons of the cubic-quintic NLS with a delta-function potential Stable solitons of the cubic-quintic NLS with a delta-function potential François Genoud TU Delft Besançon, 7 January 015 The cubic-quintic NLS with a δ-potential We consider the nonlinear Schrödinger

More information

2.4 Eigenvalue problems

2.4 Eigenvalue problems 2.4 Eigenvalue problems Associated with the boundary problem (2.1) (Poisson eq.), we call λ an eigenvalue if Lu = λu (2.36) for a nonzero function u C 2 0 ((0, 1)). Recall Lu = u. Then u is called an eigenfunction.

More information

Functional Analysis Review

Functional Analysis Review Outline 9.520: Statistical Learning Theory and Applications February 8, 2010 Outline 1 2 3 4 Vector Space Outline A vector space is a set V with binary operations +: V V V and : R V V such that for all

More information

On Schrödinger equations with inverse-square singular potentials

On Schrödinger equations with inverse-square singular potentials On Schrödinger equations with inverse-square singular potentials Veronica Felli Dipartimento di Statistica University of Milano Bicocca veronica.felli@unimib.it joint work with Elsa M. Marchini and Susanna

More information

1 Distributions (due January 22, 2009)

1 Distributions (due January 22, 2009) Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =

More information

A Fourier-Bessel Expansion for Solving Radial Schrödinger Equation in Two Dimensions

A Fourier-Bessel Expansion for Solving Radial Schrödinger Equation in Two Dimensions A Fourier-Bessel Expansion for Solving Radial Schrödinger Equation in Two Dimensions H. TAŞELI and A. ZAFER Department of Mathematics, Middle East Technical University, 06531 Ankara, Turkey ABSTRACT The

More information

Landesman-Lazer type results for second order Hamilton-Jacobi-Bellman equations

Landesman-Lazer type results for second order Hamilton-Jacobi-Bellman equations Author manuscript, published in "Journal of Functional Analysis 258, 12 (2010) 4154-4182" Landesman-Lazer type results for second order Hamilton-Jacobi-Bellman equations Patricio FELMER, Alexander QUAAS,

More information

Physics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory

Physics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three

More information

Analysis-3 lecture schemes

Analysis-3 lecture schemes Analysis-3 lecture schemes (with Homeworks) 1 Csörgő István November, 2015 1 A jegyzet az ELTE Informatikai Kar 2015. évi Jegyzetpályázatának támogatásával készült Contents 1. Lesson 1 4 1.1. The Space

More information

Part 2 Continuous functions and their properties

Part 2 Continuous functions and their properties Part 2 Continuous functions and their properties 2.1 Definition Definition A function f is continuous at a R if, and only if, that is lim f (x) = f (a), x a ε > 0, δ > 0, x, x a < δ f (x) f (a) < ε. Notice

More information

Lecture 3 : Bifurcation Analysis

Lecture 3 : Bifurcation Analysis Lecture 3 : Bifurcation Analysis D. Sumpter & S.C. Nicolis October - December 2008 D. Sumpter & S.C. Nicolis General settings 4 basic bifurcations (as long as there is only one unstable mode!) steady state

More information

PCA with random noise. Van Ha Vu. Department of Mathematics Yale University

PCA with random noise. Van Ha Vu. Department of Mathematics Yale University PCA with random noise Van Ha Vu Department of Mathematics Yale University An important problem that appears in various areas of applied mathematics (in particular statistics, computer science and numerical

More information

On Chern-Simons-Schrödinger equations including a vortex point

On Chern-Simons-Schrödinger equations including a vortex point On Chern-Simons-Schrödinger equations including a vortex point Alessio Pomponio Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Workshop in Nonlinear PDEs Brussels, September 7

More information

Quantum Mechanics Solutions

Quantum Mechanics Solutions Quantum Mechanics Solutions (a (i f A and B are Hermitian, since (AB = B A = BA, operator AB is Hermitian if and only if A and B commute So, we know that [A,B] = 0, which means that the Hilbert space H

More information

Resolvent expansions for the Schrödinger operator on a graph with rays

Resolvent expansions for the Schrödinger operator on a graph with rays Resolvent expansions for the Schrödinger operator on a graph with rays Kenichi ITO (Kobe University) joint work with Arne JENSEN (Aalborg University) 6 September 2017 Purpose of this talk: Resolvent expansion

More information

SPECTRAL PROPERTIES AND NODAL SOLUTIONS FOR SECOND-ORDER, m-point, BOUNDARY VALUE PROBLEMS

SPECTRAL PROPERTIES AND NODAL SOLUTIONS FOR SECOND-ORDER, m-point, BOUNDARY VALUE PROBLEMS SPECTRAL PROPERTIES AND NODAL SOLUTIONS FOR SECOND-ORDER, m-point, BOUNDARY VALUE PROBLEMS BRYAN P. RYNNE Abstract. We consider the m-point boundary value problem consisting of the equation u = f(u), on

More information

MATH 220: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY

MATH 220: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY MATH 22: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY When discussing separation of variables, we noted that at the last step we need to express the inhomogeneous initial or boundary data as

More information

September Math Course: First Order Derivative

September Math Course: First Order Derivative September Math Course: First Order Derivative Arina Nikandrova Functions Function y = f (x), where x is either be a scalar or a vector of several variables (x,..., x n ), can be thought of as a rule which

More information

Algebras of singular integral operators with kernels controlled by multiple norms

Algebras of singular integral operators with kernels controlled by multiple norms Algebras of singular integral operators with kernels controlled by multiple norms Alexander Nagel Conference in Harmonic Analysis in Honor of Michael Christ This is a report on joint work with Fulvio Ricci,

More information

SCHRÖDINGER OPERATORS ON GRAPHS AND GEOMETRY III. GENERAL VERTEX CONDITIONS AND COUNTEREXAMPLES.

SCHRÖDINGER OPERATORS ON GRAPHS AND GEOMETRY III. GENERAL VERTEX CONDITIONS AND COUNTEREXAMPLES. ISSN: 1401-5617 SCHRÖDINGER OPERATORS ON GRAPHS AND GEOMETRY III. GENERAL VERTEX CONDITIONS AND COUNTEREXAMPLES. Pavel Kurasov, Rune Suhr Research Reports in Mathematics Number 1, 2018 Department of Mathematics

More information

Nodal lines of Laplace eigenfunctions

Nodal lines of Laplace eigenfunctions Nodal lines of Laplace eigenfunctions Spectral Analysis in Geometry and Number Theory on the occasion of Toshikazu Sunada s 60th birthday Friday, August 10, 2007 Steve Zelditch Department of Mathematics

More information