Ministry of Education and Science of the Russian Federation Saint Petersburg National Research University of Information
Technologies, Mechanics, and Optics
NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS
2015, volume 6 (1)
Special Issue
INTERNATIONAL CONFERENCE "MATHEMATICAL CHALLENGE OF QUANTUM TRANSPORT IN NANOSYSTEMS - 2014. PIERRE DUCLOS WORKSHOP"
Наносистемы: физика, химия, математика 2015, том 6, № 1
N A N O NANOSYSTEMS: Ф &Х &М PHYSICS, CHEMISTRY, MATHEMATICS
ADVISORY BOARD MEMBERS Chairman: V.N. Vasiliev (St. Petersburg, Russia), V.M. Ievlev (Voronezh), P.S. Kop’ev(St. Petersburg, Russia), V.N. Parmon (Novosibirsk), A.I. Rusanov (St. Petersburg, Russia).
EDITORIAL BOARD Editor-in-Chief: N.F. Morozov (St. Petersburg, Russia) Deputy Editor-in-Chief: I.Yu. Popov (St. Petersburg, Russia)
Section Co-Editors: Physics – V.M. Uzdin (St. Petersburg, Russia), Chemistry, material science –V.V. Gusarov (St. Petersburg, Russia), Mathematics – I.Yu. Popov (St. Petersburg, Russia).
Editorial Board Members: V.M. Adamyan (Odessa, Ukraine); O.V. Al’myasheva (St. Petersburg, Russia); A.P. Alodjants (Vladimir, Russia); S. Bechta (Stockholm, Sweden); J. Behrndt (Graz, Austria); M.B. Belonenko (Volgograd, Russia); V.G. Bespalov (St. Petersburg, Russia); J. Brasche (Clausthal, Germany); A. Chatterjee (Hyderabad, India); S.A. Chivilikhin (St. Petersburg, Russia); A.V. Chizhov (Dubna, Russia); P.P. Fedorov (Moscow, Russia); E.A. Gudilin (Moscow, Russia); H. Jónsson (Reykjavik, Iceland); A.A. Kiselev (Madison, USA); Yu.S. Kivshar (Canberra, Australia); S.A. Kozlov (St. Petersburg, Russia); P.A. Kurasov (Stockholm, Sweden); A.V. Lukashin (Moscow, Russia); V.A. Margulis (Saransk, Russia); I.V. Melikhov (Moscow, Russia); G.P. Miroshnichenko (St. Petersburg, Russia); I.Ya. Mittova (Voronezh, Russia); H. Neidhardt (Berlin, Germany); K. Pankrashkin (Orsay, France); B.S. Pavlov (Auckland, New Zealand); A.V. Ragulya (Kiev, Ukraine); V. Rajendran (Tamil Nadu, India); A.A. Rempel (Ekaterinburg, Russia); V.P. Romanov (St. Petersburg, Russia); V.Ya. Rudyak (Novosibirsk, Russia); D Shoikhet (Karmiel, Israel); P Stovicek (Prague, Czech Republic); V.M. Talanov (Novocherkassk, Russia); A.Ya. Vul’ (St. Petersburg, Russia); V.A. Zagrebnov (Marseille, France).
Editors: I.V. Blinova; A.I. Popov; A.I. Trifanov; E.S. Trifanova (St. Petersburg, Russia), R. Simoneaux (Philadelphia, Pennsylvania, USA).
Address: University ITMO, Kronverkskiy pr., 49, St. Petersburg 197101, Russia. Phone: +7(812)232-67-65, Journal site: http://nanojournal.ifmo.ru/, E-mail: [email protected]
AIM AND SCOPE The scope of the journal includes all areas of nano-sciences. Papers devoted to basic problems of physics, chemistry, material science and mathematics inspired by nanosystems investigations are welcomed. Both theoretical and experimental works concerning the properties and behavior of nanosystems, problems of its creation and application, mathematical methods of nanosystem studies are considered. The journal publishes scientific reviews (up to 30 journal pages), research papers (up to 15 pages) and letters (up to 5 pages). All manuscripts are peer-reviewed. Authors are informed about the referee opinion and the Editorial decision. Content
From the Editorial Board 5
INTERNATIONAL CONFERENCE "MATHEMATICAL CHALLENGE OF QUANTUM TRANSPORT IN NANOSYSTEMS - 2014. PIERRE DUCLOS WORKSHOP"
INVITED SPEAKERS
H. Neidhardt, L. Wilhelm, V.A. Zagrebnov A new model for quantum dot light emitting-absorbing devices: proofs and supplements 6
K. Pankrashkin On the Robin eigenvalues of the Laplacian in the exterior of a convex polygon 46
CONTRIBUTED TALKS
E.V. Alexandrov Thermally induced transitions and minimum energy paths for magnetic systems 57
S.A. Avdonin, A.S. Mikhaylov, V.S. Mikhaylov On some applications of the boundary control method to spectral estimation and inverse problems 63
C. Cacciapuoti On the derivation of the Schrödinger equation with point-like nonlinearity 79
A.E. Ivanova, S.A. Chivilikhin, I.Yu. Popov, A.V. Gleim On the possibility of using optical y-splitter in quantum random number generation systems based on fluctuations of vacuum 95
A. Mantile A linearized model of quantum transport in the asymptotic regime of quantum wells 100
T.F. Pankratova Tunneling in multidimensional wells 113
A.G. Petrashen, N.V. Sytenko Polarization characteristics of radiation of atomic ensemble under coherent excitation in the presence of a strong magnetic field 122
A.B. Plachenov, B.T. Plachenov Border effect, fractal structures and brightness oscillations in non-silver photographic materials 133
A.I. Popov, I.S. Lobanov, I.Yu. Popov, T.V. Gerya On the Stokes flow computation algorithm based on woodbury formula 140
M.A. Skryabin Strong solutions and the initial data space for some non-uniformly parabolic equations 146
FROM THE EDITORIAL BOARD
In this issue we publish the Proceedings of the International Conference Mathematical Challenge of Quantum Transport in Nanosystems. “Pierre Duclos Workshop” organized by Saint Petersburg National Research University of Information Technologies, Mechanics, and Optics in September 2014. This workshop was the eighth of a series which has the aim to bring together specialists in nanosystem modeling, mathematicians and condensed matter physicists.
PROGRAM COMMITTEE Vladimir N. Vasilyev (ITMO University, Russia) – Chairman Igor Yu. Popov (ITMO University, Russia) – Vice-Chairman Sergey N. Naboko (St. Petersburg State University, Russia) Boris S. Pavlov (Massey University, New Zealand) George P. Miroshnichenko (ITMO University, Russia)
LOCAL ORGANIZING COMMITTEE Leonid M. Studenikin – Chairman Igor Yu. Popov – Vice-Chairman Igor S. Lobanov Ekaterina S. Trifanova Irina V. Blinova Alexander I. Trifanov Maxim A. Skryabin Anton I. Popov
INVITED SPEAKERS Pavel Exner (Prague, Czech Republic) Hagen Neidhardt (Berlin, Germany) Valentin Zagrebnov (Marseille, France) Konstantin Pankrashkin (Paris, France)
THE MAIN TOPICS OF THE CONFERENCE: Spectral theory Scattering Quantum transport Quantum communications and computations
NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS, 2015, 6 (1), P. 6–45
A NEW MODEL FOR QUANTUM DOT LIGHT EMITTING-ABSORBING DEVICES: PROOFS AND SUPPLEMENTS
1H. Neidhardt, 2L. Wilhelm, 3V.A. Zagrebnov 1,2WIAS Berlin, Mohrenstr. 39, 10117 Berlin, Germany 3Institut de Mathématiques de Marseille - UMR 7373 CMI - Technopôle Château-Gombert 39, rue F. Joliot Curie 13453 Marseille Cedex 13, France [email protected] [email protected]
PACS 85.60.Jb, 88.40.H, 85.35.Be, 73.63.Kv, 42.50.Pq DOI 10.17586/2220-8054-2015-6-1-6-45
Motivated by the Jaynes-Cummings (JC) model, we consider here a quantum dot coupled simultaneously to a reservoir of photons and to two electric leads (free-fermion reservoirs). This new Jaynes-Cummings-Leads (JCL) model permits a fermion current through the dot to create a photon flux, which describes a light-emitting device. The same model is also used to describe the transformation of a photon flux into a current of fermions, i.e. a quantum dot light-absorbing device. The key tool to obtain these results is the abstract Landauer-Büttiker formula. Keywords: Landauer-Büttiker formula, Jaynes-Cummings model, coupling to leads, light emission, solar cells. Received: 22 December 2014
1. Introduction In the following the fermion current going through a quantum dot is analyzed as a function of (1) the electro-chemical potentials on leads and of (2) the contact with the external photon reservoir. Although the latter is the canonical JC-interaction, the coupling of the JC model with leads needs certain precautions, if we wish to remain in the one-particle quantum mechanical Hamiltonian approach and scattering theory framework. To this end, we proposed a new Jaynes-Cummings- Leads (JCL-) model [19]. This model makes possible a photon flux into the resonator, created by a fermion current through the dot; i.e. it describes a light-emitting device, as well as to transform the external photon flux into a current of fermions, which corresponds to a quantum dot light- absorbing device. The paper is an extended version of [19], which means that the JCL-model, as well as all theorems, corollaries of this article one can already find in [19], however, without any proofs. In the following article, we are going to close this gap and give full proofs of all statements. In doing so, we have added some statements absent in [19]. The paper is organized as follows. The JCL-model is introduced and discussed in Sections 2.1-2.7. For simplicity, we choose for the lead Hamiltonians the one-particle discrete Schrödinger operators, with constant one-site (electric) potentials on each of leads. In Section 2.5, we show that the our model fits into the framework of trace-class scattering. In Section 2.7, we verify Quantum dot light emitting-absorbing devices 7 the important property that the coupled Hamiltonian has no singular continuous spectrum. Our main tool for analysis of different currents is an abstract Landauer-Büttiker-type formula applied in Sections 3.1 and 3.2 to the case of the JCL-model. This allows us to calculate the outgoing flux of photons induced by electric current via leads. This corresponds to a light-emitting device. We also found that pumping the JCL quantum dot with a photon flux from a resonator may induce a current of fermions into the leads. This reversing imitates a quantum light-absorbing cell device. These are the main properties of our model and the main application of the Landauer-Büttiker- type formula of Sections 3.1 and 3.2. They are presented in Sections 4 and 5, where we distinguish contact-induced and photon-induced fermion currents. To describe the results of Sections 4 and 5, note that in our setup, the sample Hamiltonian is a two-level quantum dot decoupled from the one-mode photon resonator. Then, the unperturbed Hamiltonian H0 describes a collection of four totally decoupled sub-systems: a sample, a resonator and two leads. The perturbed Hamiltonian H is a fully coupled system, and the feature of our model is that it is totally (i.e. including the leads) embedded in the external electromagnetic field of the resonator. Hence, each electron can be interpreted as a fermion with internal harmonic degrees of freedom, or a Fermi-particle caring its individual photon cloud. Similar to the “Black Box” system-leads (SL-) model {HSL,H0}, [1], [4], it turns out that the JCL-model also fits into the framework of the abstract Landauer-Büttiker formula, and in particular, is a trace-class scattering system {HJCL = H,HSL}. The current in the SL-model c is called the contact-induced current Jel. This current was the subject of numerous papers, see e.g. [1,5], or [4]. Note that the current Jel is due to the difference of the electro-chemical potentials between two leads, but it may be zero even if this difference is not null [12, 13]. The fermion current in the JCL-model, takes into account the effect of the electron-photon interaction under the assumption that the leads are already coupled. This is called the photon- ph induced component Jel of the total current. To the best of our knowledge, the present paper is the first which rigorously studies this phenomenon. We show that the total free-fermion current Jel in the JCL-model decomposes into a sum of the contact- and the photon-induced currents: c ph Jel := Jel + Jel . An extreme case is when the contact-induced current is zero, but the photon- induced component is not, c.f. Section 5.1. In this case, the flux of photons Jph out of the quantum dot (sample) is also non-zero, i.e. the dot serves as a light emitting device, c.f. Section 5.2. In ph general the Jph 6= 0 only when the photon-induced component Jel 6= 0. By choosing the parameters of the model in an suitable manner, one can get either a photon emitting or a photon absorbing system. Hence, the JCL-model can be used either as a light ph emission device or as a solar cell. Proofs of explicit formulas for fermion and photon currents Jl,el, Jph are contained in Sections 4 and 5. Note that the JCL-model is called mirror symmetric if (roughly speaking) one can inter- change left and right leads and the JCL-model remains unchanged. In Section 5 we discuss a surprising example of a mirror symmetric JCL-model such that the free-fermion current is zero but the model is photon emitting. This peculiarity is due to a specific choice of the photon-electron interaction, which produces fermions with internal harmonic degrees of freedom.
2. Jaynes-Cummings quantum dot coupled to leads 2.1. Jaynes-Cummings model The starting point for the construction of our JCL-model is the quantum optics Jaynes-Cummings Hamiltonian HJC . Its simplest version is a two-level system (quantum dot) with the energy spacing 2 ε, defined by Hamiltonian hS on the Hilbert space hS = C , see e.g. [16]. It is assumed that this system is “open” and interacts with the one-mode ω photon resonator with Hamiltonian hph. 8 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
Since mathematically hph coincides with a quantum harmonic oscillator, the Hilbert space ph of the resonator is the boson Fock space h = F+(C) over C and hph = ω b∗b . (2.1) Here, b∗ and b are verifying the Canonical Commutation Relations (CCR) creation and annihila- 2 tion operators with domains in F+(C) ' ` (N0). Operator (2.1) is self-adjoint on its domain: ph 2 X 2 2 dom(h ) = (k0, k1, k2,...) ∈ ` (N0): n |kn| < ∞ . n∈N0 2 Note that the canonical basis {φn := (0, 0, . . . , kn = 1, 0,...)}n∈N0 in ` (N0) consists of eigen- ph vectors of operator (2.1): h φn = nω φn. To model the two-level system with the energy spacing ε, one fixes in C2 two ortho-normal 0! 1! vectors {eS, eS}, for example eS := and eS := , which are eigenvectors of Hamiltonian 0 1 0 1 1 0 S S hS with eigenvalues {λ0 = 0, λ1 = ε}. To this end, we set: 1 0! h := ε , (2.2) S 0 0 and we introduce two ladder operators: 0 1! 0 0! σ+ := and σ− := . (2.3) 0 0 1 0
+ − S + S S − S − S S Then, one obtains hS = ε σ σ as well as e1 = σ e0 , e0 = σ e1 and σ e0 = 0. So, e0 is the ground state of Hamiltonian hS. Note that the non-interacting Jaynes-Cummings Hamiltonian JC JC ph 2 H0 resides in the space H = hS ⊗ h = C ⊗ F+(C) and it is defined as the matrix operator: JC ph ph H0 := hS ⊗ Ih + IhS ⊗ h . (2.4) ph ph Here, Ih denotes the identity operator in the Fock space h , whereas IhS stays for the identity matrix in the space hS. With operators (2.3), the interaction VSb between quantum dot and photons (bosons) in the resonator is defined (in the rotating-wave approximation [16]) by the operator: + − ∗ VSb := gSb (σ ⊗ b + σ ⊗ b ) . (2.5) Operators (2.4) and (2.5) define the Jaynes-Cummings model Hamiltonian: JC HJC := H0 + VSb , (2.6) JC which is a self-adjoint operator on the common domain dom(H0 ) ∩ dom(VSb). The standard interpretation of HJC is that (2.6) describes an “open” two-level system interacting with an external one-mode electromagnetic field [16]. Since the one-mode resonator is able to absorb infinitely many bosons, this interpretation sounds reasonable, but one can see that the spectrum σ(HJC ) of the Jaynes-Cummings model is + − ∗ ph discrete. To this end, note that the so-called number operator NJC := σ σ ⊗ Ih + IhS ⊗ b b commutes with HJC . Then, since for any n > 0: JC S S JC S Hn>0 := {ζ0e0 ⊗ φn + ζ1e1 ⊗ φn−1}ζ0,1∈C , Hn=0 := {ζ0e0 ⊗ φ0}ζ0∈C , JC JC are eigenspaces of operator NJC , they reduce HJC , i.e. HJC : Hn → Hn . Note that JC L JC JC H = n>0 Hn , where each Hn is invariant subspace of operator (2.6). Therefore, it has Quantum dot light emitting-absorbing devices 9 the representation: M (n) (0) HJC = HJC , n > 1 ,HJC = 0 . (2.7) n∈N0 (n) JC Here operators HJC are the restrictions of HJC , which act in each Hn as follows:
(n) S S H (ζ0 e ⊗ φn + ζ1 e ⊗ φn−1) JC 0 1 (2.8) √ S √ S = [ζ0nω + ζ1gSb n] e0 ⊗ φn + [ζ1(ε + (n − 1)ω) + ζ0gSb n] e1 ⊗ φn−1 . S (n) (n) Hence, the spectrum σ(HJC ) = n>0 σ(HJC ). By virtue of (2.8), the spectrum σ(HJC ) is defined (n) for n > 1 by eigenvalues E(n) of two-by-two matrix HcJC acting on the coefficient space {ζ0, ζ1}: ! √ ! ! ! (n) ζ1 ε + (n − 1) ω gSb n ζ1 ζ1 HcJC = √ = E(n) . (2.9) ζ0 gSb n nω ζ0 ζ0
Then, (2.7) and (2.9) imply that the spectrum of the Jaynes-Cummings model Hamiltonian HJC is pure point: 1 q σ(H ) = σ (H ) = {0} ∪ [ nω + (ε − ω) ± (ε − ω)2/4 + g2 n . JC p.p. JC 2 Sb n∈N
This property evidently persists for any system Hamiltonian hS with discrete spectrum and linear interaction (2.5) with a finite mode photon resonator [16]. We resume the above observations concerning the Jaynes-Cummings model, which is our starting point, by following remarks: (a) The standard Hamiltonian (2.6) describes instead of flux only oscillations of photons be- tween resonator and quantum dot, i.e. the system hS is not “open” enough. (b) Since one our aim is to model a light-emitting device, the system hS needs an external source of energy to pump it into the dot, which will be transformed by interaction (2.5) into the outgoing photon current by pumping the resonator. (c) To reach this aim we extend the standard Jaynes-Cummings model to our JCL-model by attaching to the quantum dot hS (2.2) two leads, which are (infinite) reservoirs of free fermions. Manipulating with electro-chemical potentials of fermions in these reservoirs we can force one of them to inject fermions in the quantum dot, whereas another one to absorb the fermions out the quantum dot with the same rate. This current of fermions throughout the dot will pump the dot and induce a photon current according scenario (b). (d) The most subtle point is to invent a leads-dot interaction VlS, which ensures the above mechanism and which is simple enough that one would still be able to treat this JCL- model using our extension of the Landauer-Büttiker formalism.
2.2. The JCL-model First, let us make some general remarks and formulate certain indispensable conditions when one follows the modeling (d). (1) Note that since the Landauer-Büttiker formalism [13] is essentially a scattering theory on a contact between two subsystems, it is developed only on a “one-particle” level. This allows one to study with this formalism only ideal (non-interacting) many-body systems. We impose this condition on many-body fermion systems (electrons) in two leads. Thus, only direct interaction between different components of the system: dot-photons VSb and electron-dot VlS are allowed. 10 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
(2) It is well-known that fermion reservoirs are technically simpler to treat than boson ones [13]. Moreover, in the framework of our model, it is also very natural since we study electric current, although produced by “non-interacting electrons”. So, below we use fermions/electrons synonymously. (3) In spite of the precautions formulated above, the first difficulty to consider in an ideal many- body system interacting with quantized electromagnetic field (photons) is induced indirect interaction. If electrons can emit and absorb photons, it is possible for one electron to emit a photon that another electron absorbs, thus creating an indirect photon-mediated electron- electron interaction. This interaction makes it impossible to develop the Landauer-Büttiker formula, which requires a non-interacting framework. Assumption 2.1. To solve this difficulty, we forbid in our model the photon-mediated interaction. To this end, we assume that every electron (in leads and in dot) interacts with its own distinct copy of the electromagnetic field. So, considering electrons together with their photon fields as non- interacting “composed particles”, allows us to apply the Landauer-Büttiker approach. Formally, it corresponds to the “one-electron” Hilbert space hel ⊗ hph, where hph is the Hilbert space of the individual photon field. The fermion description of composed-particles hel ⊗ hph corresponds to el ph the antisymmetric Fock space F−(h ⊗ h ). The composed-particle assumption 2.1 allows us to use the Landauer-Büttiker formalism developed for ideal many-body fermion systems. Now, we come closer to the formal description of our JCL-model with two (infinite) leads and a one-mode quantum resonator. Recall that the Hilbert space of the Jaynes-Cummings Hamiltonian with two energy levels JC 2 is H = C ⊗F+(C). The boson Fock space is constructed from a one-dimensional Hilbert space, since we consider only photons of a single fixed frequency. We model the electrons in the leads as free fermions residing on discrete semi-infinite lattices. Thus: el 2 2 2 el el h = ` (N) ⊕ C ⊕ ` (N) = hl ⊕ hS ⊕ hr , (2.10) el is the one-particle Hilbert space for the electrons and for the dot. Here, hα , α ∈ {l, r}, are the 2 respective Hilbert spaces of the left and right lead, while hS = C is the Hilbert space of the quantum dot. We denote by: α S 1 {δn }n∈N, {δn }j=0, el the canonical basis consisting of individual lattice sites of hα , α ∈ {l, r}, and of hS, respectively. ph 2 With the Hilbert space for photons, h = F+(C) ' ` (N0), we define the Hilbert space of the full system, i.e. quantum dot with leads and with the photon field, as follows: el ph 2 2 2 2 H = h ⊗ h = ` (N) ⊕ C ⊕ ` (N) ⊗ ` (N0). (2.11) Remark 2.2. Note that the structure of full space (2.11) takes into account the condition 2.1 and produces composed fermions via the last tensor product. It also manifests that electrons in the dot as well as those in the leads are composed with photons. This is different than the picture imposed by the the Jaynes-Cummings model, when only the dot is composed with photons: 2 2 2 2 JC 2 2 H = ` (N) ⊕ C ⊗ ` (N0) ⊕ ` (N) , H = C ⊗ ` (N0) , (2.12) JC ph see (2.4), (2.5) and (2.6), where H = hS ⊗ h . The next step is a choice of interactions between subsystems: dot-resonator-leads. According to (2.10), the decoupled leads-dot Hamiltonian is the matrix operator: el hl 0 0 ul el 2 2 h0 = 0 hS 0 on u = uS , {uα ∈ ` (N)}α∈{l,r} , uS ∈ C , el 0 0 hr ur Quantum dot light emitting-absorbing devices 11
el D where hα = −∆ + vα with a constant potential bias vα ∈ R, α ∈ {l, r}, and hS can be any S S S self-adjoint two-by-two matrix with eigenvalues {λ0 , λ1 := λ0 + ε}, ε > 0, and eigenvectors S S D 2 {e0 , e1 }, cf (2.2). Here, ∆ denotes the discrete Laplacian on ` (N) with homogeneous Dirichlet boundary conditions given by: D (∆ f)(x) := f(x + 1) − 2f(x) + f(x − 1), x ∈ N, D 2 dom(∆ ) := {f ∈ ` (N0): f(0) := 0}, which is obviously a bounded self-adjoint operator. Notice that σ(∆D) = [0, 4]. We define the lead-dot interaction for coupling gel ∈ R by the matrix operator acting in (2.10) as follows: S l 0 h·, δ0 iδ1 0 l S r S vel = gel h·, δ1iδ0 0 h·, δ1iδ1 , (2.13) S r 0 h·, δ1 iδ1 0 where non-trivial off-diagonal entries are projection operators in the Hilbert space (2.10) with the el S S el scalar product u, v 7→ hu, vi for u, v ∈ h . Here, {δ0 , δ1 } is ortho-normal basis in hS , which S S in general may be different from {e0 , e1 }. Hence, interaction (2.13) describes quantum tunneling l r between leads and the dot via contact sites of the leads, which are supports of δ1 and δ1. Then we define the Hamiltonian for the system of interacting leads and dot as hel := el el el el h0 + vel. Here, both h0 and h are bounded self-adjoint operators on h . Recall that photon Hamiltonian in the one-mode resonator is defined by operator hph = ∗ 2 ωb b with domain in the Fock space F+(C) ' ` (N0), (2.1). We denote the canonical basis in 2 ph ` (N0) by {Υn}n∈N0 . Then for the spectrum of h one obviously gets: ph ph [ σ(h ) = σpp(h ) = {nω}. (2.14) n∈N0
We introduce the following decoupled Hamiltonian H0, which describes the system when the leads are decoupled from the quantum dot and the electron does not interact with the photon field: el ph H0 := H0 + H , (2.15) where el el ph ph H0 := h0 ⊗ Ihph and H := Ihel ⊗ h . ph el ph The operator H0 is self-adjoint on dom(H0) = dom(Ihel ⊗h ). Recall that h0 and h are bounded el el self-adjoint operators. Hence, H0 and H are semi-bounded from below, which yields that H0 is semi-bounded from below. The interaction of the photons and the electrons in the quantum dot is given by the coupling of the dipole moment of the electrons to the electromagnetic field in the rotating wave approxima- tion. Namely: S S S S ∗ Vph = gph (·, e0 )e1 ⊗ b + (·, e1 )e0 ⊗ b , (2.16) for some coupling constant gph ∈ R. The total Hamiltonian is given by: el ph H := H + H + Vph = H0 + Vel + Vph, (2.17) el el where H := h ⊗ Ihph and Vel := vel ⊗ Ihph . In the following, we call S = {H,H0} the Jaynes-Cummings-leads system, in short JCL- model, which we are going to analyze. In particular, we are interested in the electron and photon currents for that system. The analysis will be based on the abstract Landauer-Büttiker formula, cf. [1, 13].
Lemma 2.3. H is bounded from below self-adjoint such that dom(H) = dom(H0). 12 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
Proof. Let c > 2. Then, 2 ∗ 2 −1 2 kbΥnk 6 kb Υnk = n + 1 6 c n + c, n ∈ N0. ph ph Consider elements f ∈ hS ⊗ h ∩ dom(Ihel ⊗ h ) with the following: X f = βjlej ⊗ Υl, j ∈ {0, 1}, l ∈ N0, j,l JC el ph 2 P 2 ∗ 2 P 2 2 which are dense in H := hS ⊗h . Then, kfk = j,l|βjl| and k(Ihel ⊗b b)fk = j,l=1|βjl| l . We obtain the following: S S 2 X 2 2 k((·, e1 )e0 ⊗ b)fk 6 |βjl| kbΥlk 6 j,l X 2 −1 2 −1 ∗ 2 2 |βjl| (c l + c) = c k(Ihel ⊗ b b)fk + ckfk j,l Similarly, S S ∗ 2 −1 ∗ 2 2 k((·, e1 )e0 ⊗ b )fk 6 c k(Ihel ⊗ b b)fk + c kfk . ph If c > 2 is large enough, then we obtain that Vph is dominated by H with relative bound less than el ph one. Hence, H is self-adjoint and dom(H0) = dom(H). Since H0 and Vel are bounded and H el ph is self-adjoint and bounded from below, it follows that H = H0 + H + Vel + Vph is bounded from below [17, Thm. V.4.1]. 2.3. Time reversible symmetric systems A system described by the Hamiltonian H is called time reversible symmetric if there is a conju- gation Γ defined on H such that ΓH = HΓ. Recall that Γ is a conjugation if the conditions Γ2 = I and (Γf, Γg) = (f, g), f, g ∈ H. ph ph Let hn , n ∈ N0, the subspace spanned by the eigenvector Υn in h . We set: el ph Hnα := hα ⊗ hn , n ∈ N0, α ∈ {l, r}. (2.18) Notice that M H = Hnα n∈N0,α∈{l,r} Definition 2.4. The JCL-model is called time reversible symmetric if there is a conjugation Γ acting on H such that H and H0 are time reversible symmetric and the subspaces Hnα , n ∈ N0, α ∈ {l, r}, reduces Γ. el el Example 2.5. Let γα and γS be conjugations defined by the following: el el γα fα := fα := {fα(k)}k∈N, fα ∈ hα , α ∈ {l, r}, and ! ! el el fS(0) fS(0) γS fS = γS := fS(1) fS(1) el el el el We set γ := γl ⊕ γS ⊕ γr . Further, we set: ph ph γ ψ := ψ = {ψ(n)}n∈N0 , ψ ∈ h . Let Γ := γel ⊗ γph. One easily checks that Γ is a conjugation on H = hel ⊗ hph. el ph Lemma 2.6. Let γα , α ∈ {S, l, r}, and γ be given by Example 2.5. el S S el S S (i) If the conditions γS e0 = e0 and γS e1 = e1 are satisfied, then H0 is time reversible
symmetric with respect to Γ and, moreover, the subspaces Hnα , n ∈ N0, α ∈ {l, r}, reduces Γ. Quantum dot light emitting-absorbing devices 13
el S S el S S (ii) If in addition the conditions γS δ0 = δ0 and γS δ1 = δ1 are satisfied, then JCL-model is time reversible symmetric. Proof. (i) Obviously we have el el el el ph ph ph ph γα hα = hα γα , α ∈ {l, r}, and γ h = h γ . el S S el S S el el el el el el el el If γS e0 = e0 and γS e1 = e1 are satisfied, then γS hS = hS γS which yields γ h0 = h0 γ and, el el el ph ph ph hence, ΓH0 = ΓH0. Since γ hα = hα and γ h = h one gets ΓHnα = Hnα which shows that
Hnα reduces Γ. el α α el S S (ii) Notice that γα δ1 = δ1 , α ∈ {l, r}. If in addition the conditions γS δ0 = δ0 and el S S el el el el el el γS δ1 = δ1 are satisfied, then γ vel = velγ is valid, which yields γ h = h γ . Hence, ΓH = HΓ. Together with (i), this proves that the JCL-model is time reversible symmetric. Choosing the following: 1! 0! 1 1! 1 1 ! eS := , eS := , δS := √ , δS := √ (2.19) 0 0 1 1 0 2 1 1 2 −1 el S S el S S el S S el S S one satisfies the condition γS e0 = e0 and γS e1 = e1 as well as γS e0 = e0 and γS e1 = e1 . 2.4. Mirror symmetric systems A unitary operator U acting on H is called a mirror symmetry if the following conditions are met: 0 0 UHnα = Hnα0 , α, α ∈ {l, r}, α 6= α JC JC JC el ph are satisfied. In particular, this yields UH = H , H := hS ⊗ h . Definition 2.7. The JCL-model is called mirror symmetric if there is a mirror symmetry commut- ing with H0 and H.
One can easily verify that if H0 is mirror symmetric, then 0 0 Hnα0 U = UHnα , n ∈ N0, α, α ∈ {l, r}, α 6= α , where el ph el 0 0 Hn := h ⊗ I ph + I el ⊗ h = h + nω, n ∈ 0, α, α ∈ {l, r}, α 6= α . α α hn hα n α N el In particular, this yields that vα = vα0 . Moreover, one gets UHS = HSU where HS := hS ⊗Ihph + ph Ihel ⊗ h . Notice that if H and H0 commute with the same mirror symmetry U, then also the operator el ph Hc := h ⊗ Ihph + Ihel ⊗ h also commutes with U, i.e, is mirror symmetric. S S S Example 2.8. Let S = {H,H0} be the JCL-model. Let vl = vr and let e0 and e1 as well as δ0 S and δ1 be given by (2.19). We set: el S S el S S uS e0 := e0 and uS e1 = −e1 , (2.20) as well as ph −inπ u Υn = e Υn, n ∈ N0. (2.21) el ph JC Obviously, US := uS ⊗u defines a unitary operator on H . Straightforward computation shows that: USHS = HSUS and USVph = VphUS. (2.22) Furthermore, we set: el l r el r l urlδn := δn, and ulrδn = δn, n ∈ N, (2.23) 14 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov and el 0 0 ulr el el u := 0 uS 0 . el ulr 0 0 We thus have the following:
el ∗ S l fl < fS, (uS ) δ0 > δ1 el el ∗ l S el ∗ r S vel u fS = < fr, (ulr) δ1 > δ0 + < fl, (url) δ1 > δ1 (2.24) el ∗ S r fr < fS, (uS ) δ1 > δ1
δS := √1 (eS + eS) δS := √1 (eS − eS) Since 0 2 0 1 and 1 2 0 1 we get from (2.20)
el ∗ S S el ∗ S S (uS ) δ0 = δ1 and (uS ) δ1 = δ0 . (2.25) Obviously, we then have el ∗ l r el ∗ r l (ulr) δ1 = δ1 (url) δ1 = δ1. (2.26) Inserting (2.25) and (2.26) into (2.24), we find
S l fl < fS, δ1 > δ1 el r S l S vel u fS = < fr, δ1 > δ0 + < fl, δ1 > δ1 (2.27) S r fr < fS, δ0 > δ1 Further, we have: S l fl < fS, δ1 > δ1 el l S r S u vel fS = < fl, δ1 > δ1 + < fr, δ1 > δ0 . (2.28) S r fr < fs, δ0 > δ1 el el el ph Comparing (2.27) and (2.28), we get u vel = velu . Setting U := u ⊗ u one immediately proves that UH0 = H0U and UH = HU. Since UHnα = Hnα0 , it is satisfied that S is mirror symmetric.
We note that Example 2.8 S is also time-reversible symmetric.
2.5. Spectral properties of H: first part
In the following, our goal is to apply the Landauer-Büttiker formula to the JCL-model. By Lp(H), 1 6 p 6 ∞, we denote in the following the Schatten-v.Neumann ideals. −1 −1 Proposition 2.9. If S = {H,H0} is the JCL-model, then (H + i) − (H0 + i) ∈ L1(H). In ac ac particular, the absolutely continuous parts H and H0 are unitarily equivalent. Proof. We have
−1 −1 −1 −1 (H + i) − (H0 + i) = (H0 + i) V (H + i) = −1 −1 −1 −1 −1 (H0 + i) V (H0 + i) − (H0 + i) V (H0 + i) V (H + i) , where V = H − H0 = Vel + Vph. Taking into account Lemma 2.3, it suffices to prove that −1 −1 ph ph (H0 + i) V (H0 + i) ∈ L1(H). Using the spectral decomposition of h with respect to h = L hph hph n∈N0 n , where n are the subspaces spanned by Υn, we obtain the following: −1 M el −1 (H0 + i) = (h + nω + i) ⊗ I ph . (2.29) 0 hn n∈N0 Quantum dot light emitting-absorbing devices 15
−1 −1 −1 −1 We have (H0 + i) V (H0 + i) = (H0 + i) (Vel + Vph)(H0 + i) . Since vel is a finite rank ph operator, we have kvelkL1 < ∞. Furthermore, hn is obviously one-dimensional for any n ∈ N0. Hence, kI ph kL = 1. From (2.29) and Vel = vel ⊗ I ph , we obtain the following: hn 1 h −1 −1 X el −1 el −1 k(H0 + i) Vel(H0 + i) kL1 = k(h0 + nω + i) vel(h0 + nω + i) kL1 n∈N0 X el −2 6 k(h0 + nω + i) k kvelkL1 n∈N0 el Since h0 is bounded, we get: q −1 el −1 2 −1 k(h0 + nω + i) k = sup (λ + nω) + 1 6 c(n + 1) , (2.30) el λ∈σ(h0 ) −1 −1 for some c > 0. This immediately implies that k(H0 + i) Vel(H0 + i) kL1 < ∞. −1 −1 ph ph We are going to handle (H0 + i) Vph(H0 + i) . Let pn be the projection from h onto ph hn . We have the following: −1 S S −1 (H0 + i) (·, e0 )e1 ⊗ b (H0 + i) X el −1 S S el −1 ph ph = (h0 + mω + i) (·, e0 )e1 (h0 + nω + i) ⊗ pm bpn m,n∈N0 X el −1 S S el −1 √ = (h0 + (n − 1)ω + i) (·, e0 )e1 (h0 + nω + i) ⊗ nΥn−1h·, Υni n∈N
From (2.30), we get the following: √ √ n el −1 S S el −1 2 (h0 + (n − 1)ω + i) (·, e0 )e1 (h0 + nω + i) ⊗ nΥnh·, Υni 6 c , L1 n(n + 1) n ∈ N, which yields: √ ∞ n k(H + i)−1 (·, eS)eS ⊗ b (H + i)−1k c2 X < ∞. 0 0 1 0 L1 6 n(n + 1) n∈N Since −1 S S ∗ −1 −1 S S −1 k(H0 + i) (·, e1 )e0 ⊗ b (H0 + i) kL1 = k(H0 + i) (·, e0 )e1 ⊗ b (H0 + i) kL1 , −1 −1 one gets (H0 + i) Vph(H0 + i) ∈ L1(H), which completes the proof. el Thus, the JCL-model S = {H,H0} is a L1-scattering system. Let us recall that hα = D el el 2 −∆ + vα, α ∈ {l, r}, on hl = hr = ` (N). Lemma 2.10. Let α ∈ {l, r}. We have the following: el el σ(hα ) = σac(hα ) = [vα, 4 + vα]. el The normalized generalized eigenfunctions of hα are given by: − 1 2 − 1 gα(x, λ) = π 2 (1 − (−λ + 2 + vα) /4) 4 sin arccos((−λ + 2 + vα)/2)x for x ∈ N, λ ∈ (vα, 4 + vα). Proof. We prove the absolute continuity of the spectrum by showing that:
{gα(x, λ) | λ ∈ (−2, 2)} is a complete set of generalized eigenfunctions. Note that it suffices to prove the lemma for ((∆D + 2)f)(x) = f(x + 1) + f(x − 1), f(0) = 0. 16 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
The lemma then follows by replacing λ with −λ + 2 + vα. Let λ ∈ (−2, 2) and − 1 2 − 1 g∆D (x, λ) = π 2 (1 − λ /4) 4 sin arccos(λ/2)x .
Note that g∆D (0, λ) = 0, when the boundary condition is satisfied. We substitute µ = arccos(λ/2) ∈ (0, π), i.e. λ = 2 cos(µ) and obtain the following: sin(µ(x + 1)) + sin(µ(x − 1)) = 2 sin(µx) cos(µ), 2 when g∆D (x, λ) satisfies the eigenvalue equation. It is obvious that g∆D (·, λ) ∈/ ` (N0) for λ ∈ (−2, 2). To complete the proof of the lemma, it remains to show the ortho-normality and the completeness. For the ortho-normality, we have to show that X g∆D (x, λ)g∆D (x, ν) = δ(λ − ν). x∈N ∞ Let ψ ∈ C0 (−2, 2) . We use the substitution µ = arccos(ν/2) and the relation
2 − 1 sin(arccos(y)) = (1 − y ) 2 to obtain the following: Z 2 X dν g∆D (x, λ)g∆D (x, ν)ψ(ν) −2 x∈N Z π sin(µ) sin arccos(λ/2)x sin(µx) −1 X = 2π dµ 1 1 ψ(2 cos(µ)) 0 2 2 x∈N (sin(µ)) (sin(arccos(λ/2))) 1 Z π 2 −1 X (sin(µ)) i(arccos(λ/2)−µ)x = (2π) dµ 1 e + 0 2 x∈N (sin(arccos(λ/2))) e−i(arccos(λ/2)−µ)x − ei(arccos(λ/2)+µ)x − e−i(arccos(λ/2)+µ)x ψ(2 cos(µ)).
Observe that for the Dirichlet kernel: X (eixy + e−ixy) − 1 = 2π δ(y), x∈N0 when Z 2 X dν g∆D (x, λ)g∆D (x, ν)ψ(ν) = −2 x∈N 1 Z π (sin(µ)) 2 dµ 1 (δ(arccos(λ/2) − µ) + δ(arccos(λ/2) + µ)) ψ(2 cos(µ)) = ψ(λ). 0 (sin(arccos(λ/2))) 2 In the second equality we use that the summand containing δ(arccos(λ/2) + µ) is zero since both arccos(λ/2) > 0 and µ > 0. Thus, the generalized eigenfunctions are orthonormal. Finally, using once more the substitution µ = arccos(ν/2), we obtain the following: Z 2 dν g∆D (x, ν)g∆D (y, ν) = −2 Z 2 − 1 dν 1 − (ν/2)2 2 sin arccos(ν/2)x sin arccos(ν/2)y = −2 Z π −1 −1 2π dµ (sin(µ)) sin(µ)sin(µx) sin(µy) = δxy, 0 for x, y ∈ N, when the family of generalized eigenfunctions is also complete. Quantum dot light emitting-absorbing devices 17
From these two lemmas, we obtain the following corollary that gives us the spectral prop- erties of H0.
Proposition 2.11. Let S = {H,H0} be the JCL-model. Then, σ(H0) = σac(H0)∪σpp(H0), where [ σac(H0) = [vl + nω, vl + 4 + nω] ∪ [vr + nω, vr + 4 + nω] n∈N0 and [ S σpp(H0) = {λj + nω : j = 0, 1}. n∈N0 S The eigenvectors are given by ge(m, n) = em ⊗ Υn, m = 0, 1, n ∈ N0. The generalized eigenfunc- tions are given by geα(·, λ, n) = gα(·, λ − nω) ⊗ Υn for λ ∈ σac(H0), n ∈ N0, α ∈ {l, r}.
Proof. It is well known (see e.g. [15]) that for two self-adjoint operators A and B with σsc(A) = σsc(B) = ∅, we have σsc(A ⊗ 1 + 1 ⊗ B) = ∅, σac(A ⊗ 1 + 1 ⊗ B) = σac(A) + σ(B) ∪ σ(A) + σac(B) and σpp(A ⊗ 1 + 1 ⊗ B) = σpp(A) + σpp(B).
Furthermore, if ψA(λA) and ψB(λB) are (generalized) eigenfunctions of A and B, respectively, then ψA(λA) ⊗ ψB(λB) is a (generalized) eigenfunction of A ⊗ I + I ⊗ B for the (generalized) eigenvalue λA + λB. el ph The lemma follows now with A = h0 and B = h using Lemmata 2.10 and (2.14) and S S S S S the fact that hS has eigenvectors {e0 , e1 } with eigenvalues {λ0 , λ1 = λ0 + ε}. 2.6. Spectral representation For the convenience of the reader, we define here what we mean under a spectral representation ac of the absolutely continuous part K0 of a self-adjoint operator K0 on a separable Hilbert space K. Let k be an auxiliary separable Hilbert space. We consider the Hilbert space L2(R, dλ, k). By M, we define the multiplication operator induced by the independent variable λ in L2(R, dλ, k). ac 2 ac 2 Let Φ: K (K0) −→ L (R, dλ, k) be an isometry acting from K (K0) into L (R, dλ, k) such that: ac Φdom(K0 ) ⊆ dom(M) and ac ac MΦf = ΦK0 f, f ∈ dom(K0 ). Obviously, the orthogonal projection P := ΦΦ∗ commutes with M which yields the existence of a measurable family, {P (λ)}λ∈R, such that: 2 (P fb)(λ) = P (λ) fb(λ), fb ∈ L (R, λ, k). We set L2(R, dλ, k(λ)) := PL2(R, λ, k), k(λ) := P (λ)k, and call the triplet ac 2 Π(K0 ) := {L (R, dλ, k(λ)), M, Φ} ac 2 ac a spectral representation of K0 . If {L (R, dλ, k(λ)), M, Φ} is a spectral representation of K , ac 2 ac ∗ then K is unitarily equivalent M0 := M L (R, dλ, k(λ)). Indeed, one has ΦK0 Φ = M0. ac ac The function ξK0 (λ) := dom(k(λ)), λ ∈ R, is called the spectral multiplicity function of K0 . 0 ξac (λ) ∞ λ ∈ Notice that 6 K0 6 for R. el For α ∈ {l, r}, the generalized eigenfunctions of hα define generalized Fourier transforms el el el,ac el 2 by φα : hα = hα (hα ) → L ([vα, vα + 4]) and el X el (φα fα)(λ) = gα(x, λ)fα(x), fα ∈ hα . (2.31) x∈N0 18 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
We then set: el C λ ∈ [vα, vα + 4] hα (λ) := (2.32) 0 λ ∈ R \ [vα, vα + 4]. el 2 el el One can easily verify that Π(hα ) = {L (R, dλ, hα (λ)), M, φα } is a spectral representation of el el,ac el hα = hα , α = l, r, where we always assumed implicitly that (φα fα)(λ) = 0 for λ ∈ R \ [vα, vα + 4]. Setting: el hl (λ) el 2 h (λ) := ⊕ ⊆ C , λ ∈ R, (2.33) el hr (λ) and introducing the map: el hl el el,ac el 2 el φ : h (h0 ) = ⊕ −→ L (R, dλ, h (λ)), (2.34) el hr defined by: el ! ! el φl fl fl φ f := el , where f := (2.35) φr fr fr el,ac 2 el el we obtain a spectral representation Π(h0 ) = {L (R, dλ, h (λ)), M, φ } of the absolutely con- el,ac el el el ac tinuous part h0 = hl ⊕ hr of h0 . One easily verifies that 0 ξ el (λ) 2 for λ ∈ R. 6 h0 6 Introducing: el el λmin := min{vl, vr} and λmax := max{vl + 4, vr + 4}, (2.36) ac el el one easily verifies that ξ el (λ) = 0 for λ ∈ R \ [λmin, λmax]. h0 Notice, if vr + 4 6 vl, then , λ ∈ [v , v + 4] ∪ [v , v + 4], hel(λ) = C r r l l {0}, otherwise
el ac which shows that h0 has a simple spectrum. In particular, it holds ξ el (λ) = 1 for λ ∈ [vr, vr + h0 ac 4] ∪ [vl, vl + 4] and otherwise ξ el (λ) = 0. h0 h 2 2 L h h 2 Let us introduce the Hilbert space := l (N0, C ) = n∈N0 n, n := C , n ∈ N0. el Regarding h (λ − nω) as a subspace of hn, one regards: M el h(λ) := hn(λ), hn(λ) := h (λ − nω), λ ∈ R, (2.37) n∈N0 as a measurable family of subspaces in h. Notice that 0 6 dim(h(λ)) < ∞, λ ∈ R. We consider the Hilbert space L2(R, dλ, h(λ)). ac 2 Furthermore, we introduce the isometric map Φ: H(H0 ) −→ L (R, dλ, h(λ)) defined by el ! M (φl fl(n))(λ − nω) (Φf)(λ) = el , λ ∈ R, (2.38) (φr fr(n))(λ − nω) n∈N0 where hel ⊗ hph f (n)! l n M l ∈ M hel,ac(hel) ⊗ hph = M ⊕ f (n) 0 n n∈ r n∈ n∈ el ph N0 N0 N hr ⊗ hn h L hph hph ph where ph = n∈N0 n and n is the subspace spanned by the eigenvectors Υn of h . One easily ac ac 2 verifies that Φ is an isometry acting from H (H0 ) onto L (R, dλ, h(λ)). Quantum dot light emitting-absorbing devices 19
2 ac Lemma 2.12. The triplet {L (R, dλ, h(λ)), M, Φ} forms a spectral representation of H0 , that is, Π(Hac) = {L2( , dλ, h(λ)), M, Φ} where there is a constant d ∈ such that 0 ξac (λ) 0 R N0 6 H0 6 el el λmax−λmin el el 2dmax for λ ∈ R where dmax := ω and λmax and λmin are given by (2.36). ac Proof. It remains to be shown that Φ transforms H0 into the multiplication operator M. We have el ! ac M (hl fl)(n) + nωfl(n) H0 f = el (hr fr)(n) + nωfr(n) n∈N0 which yields the following:
el el el ! ac M (φl (hl fl)(n))(λ − nω) + nω(φl fl(n))(λ − nω) (ΦH0 f)(λ) = el el el (φr (hr fr)(n))(λ − nω) + nω(φr fr(n))(λ − nω) n∈N0 el ! M λ(φl fl(n))(λ − nω) = el = (MΦf)(λ), λ ∈ R. λ(φr fr(n))(λ − nω) n∈N0 which proves the desired property. el el One easily checks that h(λ) might only be only non-trivial if λ−nω ∈ [λmin, λmax]. Hence, we obtain that h(λ) is non-trivial if the condition: λ − λel λ − λel max n min ω 6 6 ω is satisfied. Hence, ( λ − λel λ − λel ) 0 ξac (λ) 2 card n ∈ : max n min , λ ∈ . 6 H0 6 N0 ω 6 6 ω R or ( λel − λel ) 0 ξac (λ) 2card n ∈ : 0 n max max , λ ∈ . 6 H0 6 N0 6 6 ω R 0 ξac (λ) d λ ∈ Hence 6 H0 6 max for R.
In the following we denote the orthogonal projection from h(λ) onto hn(λ) by Pn(λ), h L h P λ ∈ R, cf (2.37). Since (λ) = n∈N0 n(λ) we have Ih(λ) = n∈N0 Pn(λ), λ ∈ R. Further, we introduce the following subspaces:
el hnα (λ) := hα (λ − nω), λ ∈ R, n ∈ N0. Notice that: M hn(λ) = hnα (λ), λ ∈ R, n ∈ N0. α∈{l,r}
By Pnα (λ) we denote the orthogonal projection from h(λ) onto hnα (λ), λ ∈ R. Clearly, we have P Pn(λ) = α∈{l,r} Pnα (λ), λ ∈ R.
ac Example 2.13. In general, the direct integral Π(H0 ) can be very complicated, in particular, the structure of h(λ) given by (2.37) is difficult to analyze. However, there are interesting simple cases: el 2 (i) Let v = vl = vr and 4 6 ω. In this case we have h (λ) = C for [v, v + 4] and 2, λ ∈ [v + nω, v + nω + 4], n ∈ , h(λ) = C N0 {0}, otherwise. 20 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
(ii) Let vr = 0, vl = 4, ω0 = 4. Then el h (λ) = , λ ∈ [0, 4), r C el 2 hlr(λ) = C , λ ∈ [4, 8), h(λ) = el 2 hrl(λ) = C , λ ∈ [8, 12), ··· where el hα (λ) el 0 0 hαα0 (λ) = ⊕ , α, α ∈ {l, r}, α 6= α . el hα0 (λ) Hence, dim(h(λ)) = 2 for λ > 4. ac ac Let Z be a bounded operator acting on H (H0) and commuting with H0 . Since Z com- ac mutes with H0 there is a measurable family {Z(λ)}λ∈R of bounded operators acting on h(λ) such ac that Z is unitarily equivalent to the multiplication operator induced by {Z(λ)}λ∈R in Π(H0 ). We then set: Z (λ) := P (λ)Z(λ) h (λ), λ ∈ , m, n ∈ , α, ∈ {l, r}. mαnκ mα nκ R N0 κ Let Z := P ZP where P is the orthogonal projection from H onto H ⊆ Hac(H ), mαnκ mα nκ mα mα 0 cf. (2.18). Clearly, the multiplication operator induced {Z (λ)} in Π(Hac) is unitarily mαnκ λ∈R 0 equivalent to Z . mαnκ Since, by Lemma 2.12, h(λ) is a finite dimensional space, the operators Z(λ) are finite dimensional ones and we can introduce the following quantity: σ (λ) := tr(Z (λ)∗Z (λ)), λ ∈ , m, n ∈ , α, ∈ {l, r}. mαnκ mαnκ mαnκ R N0 κ
Lemma 2.14. Let H0 be the self-adjoint operator defined by (2.15) on H. Furthermore, let Z be a ac ac bounded operator on H (H0) commuting with H0
(i) Let Γ be a conjugation on H, cf. Section 2.3. If Γ commutes with H0 and Pnα , n ∈ N0, α ∈ {l, r} and ΓZΓ = Z∗ holds, then σ (λ) = σ (λ), λ ∈ . mαnκ nκ mα R (ii) Let U be a mirror symmetry on H. If U commutes with H and Z, then σ (λ) = 0 mαnκ 0 0 0 0 σm 0 n 0 (λ), λ ∈ , m, n ∈ 0, α, α , , ∈ {l, r}, α 6= α , 6= . α κ R N κ κ κ κ ac Proof. (i) Since Γ commutes with H0 the conjugation Γ is reduce by H (H0). So without loss of ac generality, we assume that Γ acts on H (H0). We set Γnα := Γ Hnα . Notice that: M Γ = Γnα . n∈N0,α∈{l,r}
There is a measurable family {Γ(λ)}λ∈R of conjugations such that the multiplication operator ac induced by {Γ(λ)}λ∈R in Π(H0 ) is unitarily equivalent to Γ. Moreover, since Γ commutes with
Pnα , we see that the multiplication operator induced by the measurable family:
Γnα (λ) := Γ(λ) hnα (λ), λ ∈ R, m ∈ N0, α ∈ {l, r}, ∗ ∗ is unitarily equivalent to Γnα . Using ΓZΓ = Z we get Γmα Zmαn Γn = Z . Hence, κ κ nκ mα Γ (λ)Z (λ)Γ (λ) = Z (λ)∗, λ ∈ . (2.39) mα mαnκ nκ nκ mα R If X is a trace class operator, then tr(ΓXΓ) = tr(X). Using that we find σ (λ) = tr(Γ (λ)Z (λ)∗Z (λ)Γ (λ)) = mαnκ nκ mαnκ mαnκ nκ tr(Γ (λ)Z (λ)∗Γ Γ Z (λ)Γ (λ)) nκ mαnκ mα mα mαnκ nκ Quantum dot light emitting-absorbing devices 21
From (2.39), we obtain the following: σ (λ) = tr(Z (λ)Z (λ)∗) = σ (λ), λ ∈ , mαnκ nκ mα nκ mα nκ mα R which proves (i). ac (ii) Again, without loss of generality we can assume that U acts only H (H0). Since U commutes with H0, there is a measurable family {U(λ)}λ∈R of unitary operators acting on h(λ) such that the multiplication operator induced by {U(λ)}λ∈R is unitarily equivalent to U. Since
UHnα = Hnα0 we have U(λ)hnα (λ) = hnα0 (λ), λ ∈ R. Hence, σ (λ) = tr(U(λ)Z (λ)∗Z (λ)U(λ)∗) = mαnκ mαnκ mαnκ tr(U(λ)Z (λ)∗U(λ)∗U(λ)Z (λ)U(λ)∗). mα,nκ mα,nκ Hence, ∗ ∗ ∗ σmαn (λ) = tr(Pn 0 U(λ)Z(λ) U(λ) Pm 0 (λ)U(λ)Z(λ)U(λ) Pn 0 (λ)). κ κ α κ Since U commutes with Z, we find that: ∗ σmαn (λ) = tr(Pn 0 Z(λ) Pm 0 (λ)Z(λ)Pn 0 (λ)) = σm 0 n 0 (λ), λ ∈ , κ κ α κ α κ R which proves (ii). 2.7. Spectral properties of H: second part
Since we have full information for the spectral properties of H0, we can use this to show that H has no singular continuous spectrum. Crucial for that is the following lemma: with the help of [6, Cor. IV.15.19], which establishes existence and completeness of wave operators and absence of singular continuous spectrum through a time-falloff method. We cite it as a Lemma for convenience, with slight simplifications that suffice for our purposes.
Lemma 2.15 ( [6, Corollary IV.15.19]). Let {H0,H} be a scattering system and let Λ be a closed ac countable set. Let F+ and F− be two self-adjoint operators such that F+ + F− = PH0 and
∓itH0 ±itH0 s − lim e F±e = 0. t→∞ −1 −1 ac If (H − i) − (H0 − i) ∈ L∞(H), (1 − PH0 )γ(H0) ∈ L∞(H), and Z ±∞ −1 −1 −itH0 dt (H0 − i) − (H − i) e γ(H0)F± < ∞ 0 ∞ for all γ ∈ C0 (R \ Λ), then W±(H,H0) exist and are complete and σsc(H) = σsc(H0) = ∅. Furthermore, each eigenvalue of H and H0 in R \ Λ is of finite multiplicity and these eigenvalues accumulate the most at points of Λ or at ±∞. We already know that the wave operators exist and are complete since the resolvent differ- ence is trace class. Hence, we need Lemma 2.15 only to prove the following proposition. Proposition 2.16. The Hamiltonian H defined by (2.17) has no singular continuous spectrum, that is, σsc(H) = ∅. 2 2 Proof. At first we have to construct the operators F±. To this end, let F : L (R) → L (R) be the usual Fourier transform, i.e Z 1 −iµx 2 (Ff)(µ) := fb(µ) := √ e f(x)dx, f ∈ L (R, dx), µ ∈ R. 2π R 2 2 Further, let Π± be the orthogonal projection onto L (R±) in L (R). We set: ∗ ∗ F± = Φ FΠ±F Φ, 22 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov where Φ is given by (2.38). We immediately obtain F− + F+ = Pac(H0). We still have to show that: ∓itH0 ∗ ∗ ±itH0 s − lim ke Φ FΠ±F Φe fk = 0 t→∞ ac for f ∈ H (H0). We prove the relation only for F+ since the proof for F− is essentially identical. We have the following:
1 Z ∗ itH0 − 2 i(x+t)µ Π+F Φe f (x) = (2π) χR+ (x) dµ e fb(µ) = χR+ (x)ψ(x + t) R with ψ = Ffb. Now,
−itH0 ∗ ∗ itH0 2 ke Φ FΠ+F Φe fk = ∞ Z 2 Z 2 t→∞ ∗ itH0 2 kΠ+F Φe fk = dx ψ(x + t) = dx ψ(x) −→ 0. R+ t −1 −1 Concerning the compactness condition, we already know that (H − i) − (H0 − i) ∈ L1(H) ⊂ L∞(H) from Proposition 2.9. Let [ Λ := {vl + nω, vr + nω, vl + 4 + nω, vr + 4 + nω}, n∈N0 which is closed and countable. We know from Corollary 2.11 that H0 has no singular continuous spectrum and the eigenvalues are of finite multiplicity. It follows then that (1 − Pac(H0))γ(H0) is ∞ compact for every γ ∈ C0 (R \ Λ). The remaining assumption of Lemma 2.15 is as follows: Z ±∞ −1 −1 −itH0 dt (H − i) − (H0 − i) γ(H0)e F± < ∞. 0 If we can prove this, then we immediately obtain that H has no singular continuous spectrum. −1 −1 −1 −1 −1 Now (H − i) − (H0 − i) = (H − i) (Vel + Vph)(H0 − i) . But (H − i) is bounded, ac el el ph ran(F±) ⊂ H (H0) = (hl ⊕ hr ) ⊗ h , ac and VphP (H0) = 0. Also, Vel = vel ⊗ Ihph and ⊥ l r ker(vel) ⊂ Cδ1 ⊕ hS ⊕ Cδ1. Hence, it suffices to prove: Z ±∞ α −1 −itH0 dt P1 (H0 − i) γ(H0)e F± < ∞, 0 α α α el α ∈ {l, r}, where P1 = p1 ⊗ Ihph and p1 is the orthogonal projection onto hα . In the following we treat only the case F+. The calculations for F− are completely analogous. We use that Φ maps ac H0 into the multiplication operator M induced by λ. Hence, we get the following:
α −itH0 ∗ α ∗ −itH0 ∗ P1 γe(H0)e Φ Ff = P1 Φ Φγe(H0)e Φ Ff = Z Z 2 1 , − 1 X −iλ(x+t) 2 = (2π) 2 dλ gα(1, λ − nω)γe(λ) dx e f(x) δα,n R+ n∈N0 −1 where supp (f) ⊆ R+, γe(λ) := (λ − i) γ(λ), λ ∈ R, and δα,n := [vα + nω0, vα + nω + 4]. Notice ∞ that γe(λ) ∈ C0 (R \ Λ). We find: Z Z −iλ(x+t) dλ gα(1, λ − nω)γe(λ) dx e f(x) = δj,n R+ Z vα+4 Z −i(λ+nω)(x+t) dλ gα(1, λ)γe(λ + nω) dx e f(x) vα R+ Quantum dot light emitting-absorbing devices 23 which yields:
α ∗ −itH0 ∗ P1 Φ Φγe(H0)e Φ Ff = 2 1 1 Z vα+4 Z 2 − X −i(λ+nω0)(x+t) (2π) 2 dλ gα(1, λ)γe(λ + nω0) dx e f(x) . vα R+ n∈N0 P Since the support of γ(λ) is compact, we see that the sum n∈N0 is finite. Changing the integrals, we get: Z Z −iλ(x+t) dλ gα(1, λ − nω)γe(λ) dx e f(x) = δα,n R+ Z Z vα+4 −inω0(x+t) −iλ(x+t) dx f(x)e dλ gα(1, λ)γe(λ + nω)e R+ vα Integrating by parts m-times, we obtain: Z Z −iλ(x+t) dλ gα(1, λ − nω)γe(λ) dx e f(x) = δα,n R+
Z −inω(x+t) Z vα+4 m m e −iλ(x+t) d (−i) dx f(x) m dλ e m (gα(1, λ)γe(λ + nω)) R+ (x + t) vα dλ Hence, !2 Z Z 2 Z 1 −iλ(x+t) 2 dλ gα(1, λ − nω)γe(λ) dx e f(x) 6 Cn dx |f(x)| m δα,n R+ R+ (x + t) which yields: Z Z 2 1 −iλ(x+t) 2 2 dλ gα(1, λ − nω)γe(λ) dx e f(x) 6 Cn (2m−1) kfk δα,n R+ t for m ∈ N where: Z vα+4 m d Cn := dλ m gα(1, λ)γe(λ + nω) . vα dλ Notice that Cn = 0 for sufficiently large n ∈ N. Therefore, 1/2 1 α −itH0 ∗ X 2 2 P γ(H0)e Φ Ff C kfk, f ∈ L ( +, dx), 1 e 6 n tm−1/2 R n∈N0
α −itH0 1 which shows that P1 γe(H0)e F+ ∈ L (R+, dt) for m > 2. 3. Landauer-Büttiker formula and applications 3.1. Landauer-Büttiker formula The abstract Landauer-Büttiker formula can be used to calculate currents through devices. Usually one considers a pair S = {K,K0} be of self-adjoint operators where the unperturbed Hamiltonian K0 describes a totally decoupled system, that means, the inner system is closed and the leads are decoupled from it, while the perturbed Hamiltonian K describes the system where the leads are coupled to the inner system. An important component is system S = {K,K0}, which represents a complete scattering or even a trace class scattering system. In [1], an abstract Landauer-Büttiker formula was derived in the framework of a trace class scattering theory for semi-bounded self-adjoint operators which allows one to reproduce the results of [18] and [7] rigorously. In [13], the results of [1] were generalized to non-semi- bounded operators. Following [1], we consider a trace class scattering system S = {K,K0}. We 24 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov recall that S = {K,K0} is called a trace class scattering system if the resolvent difference of K and K0 belongs to the trace class. If S = {K,K0} is a trace class scattering system, then the wave operators W±(K,K0) exist and are complete. The scattering operator is defined by ∗ S(K,K0) := W+(K,K0) W−(K,K0). The main components, besides the trace class scattering system S = {K,K0}, are the density and the charge operators ρ and Q, respectively. The density operator ρ is a non-negative bounded self-adjoint operator commuting with K0. The charge Q is a bounded self-adjoint operator commuting also with K0. If K has no singular continuous spectrum, then the current related to the density operator ρ and the charge Q is defined as follows: S ∗ Jρ,Q = −i tr (W−(K,K0)ρW−(K,K0) [K,Q]) , (3.1) where [K,Q] is the commutator of K and Q. In fact, the commutator [K,Q] might be not defined. In this case, the regularized definition: 1 1 J S = −i tr W (K,K )(I + K2)ρW (K,K )∗ [K,Q] , (3.2) ρ,Q − 0 0 − 0 K − i K + i 2 is used, where it is assumed that (I + K0 )ρ is a bounded operator. Since the condition −1 −1 (H − i) [H,Q](H + i) ∈ L1(H) is satisfied, definition (3.2) makes sense. By L1(H) is the ideal of trace class operators is denoted. Let K0 be self-adjoint operator on the separable Hilbert space K. We call ρ be a density operator for K0 if ρ is a bounded non-negative self-adjoint operator commuting with K0. Since ρ ac commutes with K0, one sees that ρ leaves invariant the subspace K (K0). We then set ac ρac := ρ K (K0). call ρac the ac-density part of ρ. A bounded self-adjoint operator, Q, commuting with K0, is called a charge. If Q is the charge, then: ac Qac := Q K (K0), is called its ac-charge component. ac 2 ac Let Π(K0 ) = {L (R, dλ, k(λ)), M, Φ} be a spectral representation of K0 . If ρ is a density operator, then there is a measurable family {ρac(λ)}λ∈R of bounded self-adjoint operators such that the multiplication operator: 2 (Mρac fb)(λ) := ρac(λ) fb(λ), fb ∈ dom(Mρac ) := L (R, dλ, k(λ)), ∗ is unitarily equivalent to the ac-part ρac, that is, Mρac = ΦρacΦ . In particular, this yields that: kρ (λ)k = kρ k ac {ρ (λ)} ess-sup λ∈R ac B(k(λ) ac B(K (K0)). In the following, we call ac λ∈R the density matrix of ρac. Similarly, one obtains that if Q is a charge, then there is a measurable family {Qac(λ)}λ∈R of bounded self-adjoint operators, such that the multiplication operator:
(MQac fb)(λ) := Qac(λ) fb(λ), 2 2 fb ∈ dom(Qac) := {f ∈ L (R, dλ, k(λ)) : Qac(λ) fb(λ) ∈ L (R, dλ, k(λ))}, ∗ is unitarily equivalent to Qac, i.e. MQac = ΦQacΦ . In particular, one has:
kQ (λ)k = kQ k ac . ess-sup λ∈R ac B(k(λ)) ac B(K (K0)) (3.3)
If Q is a charge, then the family {Qac(λ)}λ∈R is called the charge matrix of the ac-component of Q.
Let S = {K,K0} be a trace scattering system. By {S(λ)}λ∈R, we denote the scattering matrix, which corresponds to the scattering operator S(K,K0) with respect to the spectral repre- ac ac sentation Π(K0 ). The operator T := S(K,K0) − P (K0) is called the transmission operator. Quantum dot light emitting-absorbing devices 25
By {T (λ)}λ∈R, we denote the transmission matrix which is related to the transmission operator. Scattering and transmission matrices are related by S(λ) = Tk(λ) + T (λ) for a.e. λ ∈ R. Notice that T (λ) belongs for to the trace class a.e. λ ∈ R.
Theorem 3.1 ( [13, Corollary 2.14]). Let S := {K,K0} be a trace class scattering system and let ac {S(λ)}λ∈R be the scattering matrix of S with respect to the spectral representation Π(K0 ). Fur- thermore, let ρ and Q be density and charge operators and let {ρac(λ)}λ∈R and {Qac(λ)}λ∈R be ac the density and charge matrices of the ac-components ρac and charge Qac with respect to Π(K0 ), 2 S respectively. If (I + K0 )ρ is bounded, then the current Jρ,Q defined by (3.2) admits the represen- tation: Z S 1 ∗ Jρ,Q = tr ρac(λ)(Qac(λ) − S (λ)Qac(λ)S(λ)) dλ, (3.4) 2π R S where the integrand on the right side and the current Jρ,Q satisfy the estimate: ∗ |tr (ρac(λ)(Qac(λ) − S (λ)Qac(λ)S(λ)))| 6 (3.5)
4kρ(λ)kL(k(λ))kT (λ)kL1(k(λ))kQ(λ)kL(k(λ)), for a.e. λ ∈ R and S −1 −1 |Jρ,Q| 6 C0k(H + i) − (H0 + i) kL1(K), (3.6) 2 2 where C0 := π k(1 + H0 )ρkL(K). In applications, not every charge Q is a bounded operator. We say the self-adjoint operator −p Q commuting with K0 is a p-tempered charge if Q(H0 − i) is a bounded operator for p ∈ N0. ac As above, we can introduce Qac := Q dom(Q) ∩ K (K0). It follows that QEK0 (∆) is a bounded operator for any bounded Borel set ∆. This yields that the corresponding charge matrix,
{Qac(λ)}λ∈R, is a measurable family of bounded self-adjoint operators such that: (1 + λ2)p/2kQ (λ)k < ∞. ess-sup λ∈R ac L(k(λ)) S To generalize the current Jρ,Q to tempered charges Q, one uses the fact that Q(∆) := QEK0 (∆) is S a charge for any bounded Borel set ∆. Hence, the current Jρ,Q(∆) is well-defined by (3.2) for any bounded Borel set ∆. Using Theorem 3.1 one gets that for p-tempered charges, the limit S S Jρ,Q := lim Jρ,Q(∆) (3.7) ∆→R p+2 exists, provided (H0 − i) ρ is a bounded operator. This gives rise to the following corollary:
Corollary 3.2. Let the assumptions of the Theorem 3.1 be satisfied. If for some p ∈ N0 the p+2 operator (H0 − i) ρ is bounded and Q is a p-tempered charge for K0, then the current defined by (3.7) admits the representation (3.4), where the right hand side of (3.4) satisfies the estimate S (3.5). Moreover, the current Jρ,Q can be estimated in the following manner: S −1 −1 |Jρ,Q| 6 Cpk(H + i) − (H0 + i) kL1(K), (3.8) 2 2 p+2/2 2 −p/2 where Cp := π k(1 + H0 ) ρkL(K)kQ(I + H0 ) kL(K). At first glance, the formula (3.4) is not very similar to the original Landauer-Büttiker for- mula of [7,18]. To make the formula more convenient, we recall that a standard application exam- ple for the Landauer-Büttiker formula is the so-called black-box model, cf. [1]. In this case, the Hilbert space K is given by: N M K = KS ⊕ Kj, 2 6 N < ∞. (3.9) j=1 26 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov and K0 by: N M K0 = KS ⊕ Kj, 2 6 N < ∞. (3.10) j=1
The Hilbert space KS is called the sample or dot and KS is the sample or dot Hamiltonian. The Hilbert spaces Kj are called reservoirs or leads and Kj are the reservoir or lead Hamiltonians. For simplicity, we assume that the reservoir Hamiltonians Kj are absolutely continuous and the sample Hamiltonian KS has a point spectrum. A typical choice for the density operator is: N M ρ = fS(KS) ⊕ fj(Kj), (3.11) j=1 where fS(·) and fj(·) are non-negative bounded Borel functions, and for the charge: N M Q = gS(Hs) ⊕ gj(Hj), (3.12) j=1 where gS(·) and gj(·) a bounded Borel functions. Making this choice the Landauer-Büttiker for- mula (3.4) takes the form: N Z S 1 X Jρ,Q = (fj(λ) − fk(λ))gj(λ)σjk(λ)dλ, (3.13) 2π j,k=1 R where ∗ σjk(λ) := tr(Tjk(λ) Tjk(λ)), j, k = 1, . . . , N, λ ∈ R, (3.14) are called the total transmission probability from reservoir k to reservoir j, cf. [1]. We call it the cross-section of the scattering process going from channel k to channel j at energy λ ∈ R. {Tjk(λ)}λ∈R is called the transmission matrix from channel k to channel j at energy λ ∈ R with ac respect to the spectral representation Π(K0 ). We note that {Tjk(λ)}λ∈R corresponds to the trans- mission operator: ac Tjk := PjT (K,K0)Pk,T (K,K0) := S(K,K0) − P (K0), (3.15) acting from the reservoir k to reservoir j where T (K,K0) is called the transmission operator. Let S {T (λ)}λ∈R be the transmission matrix. Following [1], the current Jρ,Q given by (3.13) is directed from the reservoirs into the sample. ∗ The quantity kT (λ)kL2 = tr(T (λ) T (λ)) is well-defined and is called the cross-section of the scattering system S at energy λ ∈ R. Notice that: N ∗ X σ(λ) = kT (λ)kL2 = tr(T (λ) T (λ)) = σjk(λ). λ ∈ R, j,k=1
We point out that the channel cross-sections σjk(λ) admit the property: N N X X σjk(λ) = σkj(λ), λ ∈ R, (3.16) j=1 j=1 which is a consequence of the unitarity of the scattering matrix. Moreover, if there is a conjugation J, such that KJ = JK and K0J = JK0 holds, that is, if the scattering system S is time reversible symmetric, then we have even more, namely, it holds that:
σjk(λ) = σkj(λ), λ ∈ R. (3.17) Usually, the Landauer-Büttiker formula (3.13) is used to calculated the electron current el entering the reservoir j from the sample. In this case one has to choose Q := Qj := −ePj where Quantum dot light emitting-absorbing devices 27
Pj is the orthogonal projection form K onto Kj and e > 0 is the magnitude of the elementary charge. This is equivalent to choosing gj(λ) = −e and gk(λ) = 0 for k 6= j, λ ∈ R. In doing so, we get the Landauer-Büttiker formula simplifying to: N Z S e X J el = − (fj(λ) − fk(λ))σjk(λ)dλ. (3.18) ρ,Qj 2π k=1 R To restore the original Landauer-Büttiker formula, one sets:
fj(λ) = f(λ − µj), λ ∈ R, (3.19) where µj is the chemical potential of the reservoir Kj and f(·) is a bounded non-negative Borel function called the distribution function. This gives rise to the following formula: N Z S e X J el = − (f(λ − µj) − f(λ − µk))σjk(λ)dλ. (3.20) ρ,Qj 2π k=1 R In particular, if we choose one: 1 f(λ) := f (λ) := , β > 0, λ ∈ , (3.21) FD 1 + eβλ R where fFD(·) is the Fermi-Dirac distribution function, and inserting (3.21) into (3.20) we obtain: N Z S e X J el = − (fFD(λ − µj) − fFD(λ − µk))σjk(λ)dλ. (3.22) ρ,Qj 2π k=1 R If we have only two reservoirs, then they are usually denoted by l (left) and r (right). Let j = l and k = r. Then, Z S e Jρ,Qel = − (fFD(λ − µl) − fFD(λ − µr))σlr(λ)dλ. (3.23) l 2π R J S 0 µ µ One easily checks that ρ,Ql 6 if l > r. That means, the current is leaving the left reservoir and is entering the right one which is in accordance with physical expectations.
el el el Example 3.3. Notice that sc := {h , h0 } is a L1 scattering system. The Hamiltonian h takes 2 2 into account the effect of coupling of reservoirs or leads hl := l (N) and hr := l (N) to the sample 2 hS = C which is also called the quantum dot. The Hamiltonians for the leads are given by: el D el hα = −∆ + vα, α = l, r. The sample or quantum dot Hamiltonian is given by hS . The wave operators are given by: el el ithel −ithel ac el w±(h , h ) := s- lim e e 0 P (h ). (3.24) 0 t→∞ 0 el el ∗ el el el,ac The scattering operator is given by sc := w+(h , h0 ) w−(h , h0 ). Let Π(h0 ) be the spectral el,ac el el representation of h0 introduced in Section 2.6. If ρ and q are density and charge operators for el h0 , then the Landauer-Büttiker formula takes the following form: Z sc 1 el el ∗ el Jρel,qel = tr ρac(λ) qac − sc(λ) qac(λ)sc(λ) , (3.25) 2π R el el where {sc(λ)}λ∈R, {q (λ)}λ∈R and {ρ (λ)}λ∈R are the scattering, charge and density matrices el,ac el 2 el with respect to Π(h0 ), respectively. The condition that ((h0 ) + Ihel )ρ is a bounded operator is el superfluous because h0 is a bounded operator. For the same reason we have that every p-tempered charge qel is in fact a charge, that means, qel is a bounded self-adjoint operator. el el The scattering system sc is a black-box model with reservoirs hl and hr . Choosing el el el el ρ = fl(h ) ⊕ fS(hS ) ⊕ fr(hr ), 28 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov where fα(·), α = l, r, are bounded Borel functions, and el el el el q = gl(hl ) ⊕ gS(hS ) ⊕ gr(hr ), where gα(·), α ∈ {l, r}, are locally bounded Borel functions, then from (3.13) it follows that: Z sc 1 X Jρel,qel = (fα(λ) − fκ(λ))gα(λ)σc(λ)dλ, 2π R α,κ∈{l,r} α6=κ where {σc(λ)}λ∈R is the channel cross-section from left to right and vice versa. Indeed, let {tc(λ)}λ∈R be the transition matrix which corresponds to the transition operator tc := sc − Ihel . el Clearly, one has tc(λ) = Ih(λ) − sc(λ), λ ∈ R. Let {pα (λ)}λ∈R be the matrix which corre- el el el c el el sponds to the orthogonal projection pα from h onto hα . Further, let trl(λ) := pr (λ)tc(λ)pl c el el and tlr := pl (λ)tc(λ)pr . Notice that both quantities are in fact scalar functions. Accordingly, c c c the channel cross-sections σlr(λ) and σrl(λ) at energy λ ∈ R are given by σc(λ) := σlr(λ) = c 2 c 2 c |tlr(λ)| = |trl(λ)| = σrl(λ), λ ∈ R. In particular, if gl(λ) = 1 and gr = 0, then: Z sc 1 Jρel,qel = (fl(λ) − fr(λ))σc(λ)dλ, (3.26) l 2π R el el sc and ql := pl . Following [1], J el el denotes the current entering the quantum dot from the left ρ ,ql lead. 3.2. Application to the JCL-model
Let S = {H,H0} now be the JCL-model. Furthermore, let ρ and Q be the density operator and a S charge for H0, respectively. Under these assumptions, the current Jρ,Q is defined by: 1 1 J S := −itr W (H,H )(I + H2)ρW (H,H )∗ [H,Q] , (3.27) ρ,Q − 0 0 − 0 H − i H + i p+2 and admits representation (3.4). If Q is a p-tempered charge and (H0−i) ρ is a bounded operator, S then the current Jρ,Q is defined in accordance with (3.7) and the Landauer-Büttiker formula (3.4) is also valid. We introduce the intermediate scattering system Sc := {H,Hc}, where: el ph Hc := h ⊗ Ihph + Ihel ⊗ h = H0 + Vel.
The Hamiltonian Hc describes the coupling of the leads to the quantum dot, but under the assump- tion that the photon interaction is not switched on. Accordingly, Sph := {H,Hc} and Sc := {Hc,H0} are L1-scattering systems. The ac corresponding scattering operators are denoted by Sph and Sc, respectively. Let Π(Hc ) = 2 ac {L (R, dλ, hc(λ)), M, Φc} of Hc be a spectral representation of Hc. The scattering matrix of ac the scattering system {H,Hc} with respect to Π(Hc ) is denoted by {Sph(λ)}λ∈R. The scattering ac 2 matrix of the scattering system {Hc,H0} with respect to Π(H0 ) = {L (R, dλ, h0(λ)), M, Φ0} is denoted by {Sc(λ)}λ∈R. Since Sc is a L1-scattering system, the wave operators W±(Hc,H0) exist and are complete ∗ and since ΦcW±(Hc,H0)Φ0 commutes with M, there are measurable families {W±(λ)}λ∈R of isometries acting from h0(λ) onto hc(λ) for a.e. λ ∈ R such that: ∗ 2 (ΦcW±(Hc,H0)Φ0 fb)(λ) = W±(λ) fb(λ), λ ∈ R, fb ∈ L (R, dλ, h0(λ)).
The families {W±(λ)}λ∈R are called wave matrices. ∗ Straightforward computation shows that Sbph := W+(Hc,H0) SphW+(Hc,H0) commutes ac with H0. Hence, with respect to the spectral representation Π(H0 ), the operator Sbph is unitarily Quantum dot light emitting-absorbing devices 29 equivalent to a multiplication induced by a measurable family { Sbph (λ)}λ∈R of unitary operators in h0(λ). Straightforward computation shows that: ∗ Sbph(λ) = W+(λ) Sph(λ)W+(λ), (3.28) for a.e. λ ∈ R. Roughly speaking, { Sbph (λ)}λ∈R is the scattering matrix of Sph with respect to the ac spectral representation Π(H0 ). Furthermore, let c ∗ ρ := W−(Hc,H0)ρW−(Hc,H0) (3.29) and c ∗ Q := W+(Hc,H0)QW+(Hc,H0) . (3.30) The operators ρc and Qc are the density and tempered charge operators for the scattering system c c c Sph. Indeed, one easily verifies that ρ and Q are commute with Hc. Moreover, ρ is non-negative. Furthermore, if Q is a charge, then Qc is also a charge. This gives rise to the introduction of c Sc currents Jρ,Q := Jρ,Q, c ∗ 1 1 Jρ,Q := −itr W−(Hc,H0)ρW−(Hc,H0) [Hc,Q] , (3.31) Hc − i Hc + i
ph Sph and Jρ,Q := Jρc,Qc , 1 1 J ph := −itr W (H,H )ρcW (H,H )∗ [H,Qc] , (3.32) ρ,Q − c − c H − i H + i p+2 which are well defined. If Q is p-tempered charge and (H0 − i) ρ is a bounded operator, then c p+2 c one easily checks that Q is a p-tempered charge and (Hc − i) ρ is a bounded operator. Hence Sc the definition of the currents Jρc,Qc can be extended to this case and the Landauer-Büttiker formula (3.4) holds. c c Finally, we note that the corresponding matrices {ρac(λ)}λ∈R and {Qac(λ)}λ∈R are related to the matrices {ρac(λ)}λ∈R and {Qac(λ)}λ∈R by c ∗ c ∗ ρac(λ) = W−(λ)ρac(λ)W−(λ) and Qac(λ) = W+(λ)Qac(λ)W+(λ) (3.33) for a.e. λ ∈ R.
Proposition 3.4 (Current decomposition). Let S = {H,H0} be the JCL-model. Furthermore, let ρ and Q be the density operator and a p-tempered charge, p ∈ N0, for H0, respectively. If p+2 (H0 − i) ρ is a bounded operator, then the decomposition, S c ph Jρ,Q = Jρ,Q + Jρ,Q, (3.34) c ph holds where Jρ,Q and Jρ,Q are given by (3.31) and (3.32). In particular, let {Sc(λ)}λ∈R, {ρac(λ)}λ∈R and {Qac(λ)}λ∈R be scattering, density and ac c charge matrices of Sc, ρ and Q with respect to Π(H0 ) and let {Sph(λ)}λ∈R, {ρac(λ)}λ∈R and c {Qac(λ)}λ∈R be the scattering, density and charge matrices of the scattering operator Sph, density operator ρc, cf. (3.29), and charge operator Qc, cf. (3.30), with respect to the spectral representa- ac tion Π(Hc }. Then, the following representations: Z c 1 ∗ Jρ,Q := tr(ρac(λ)(Qac(λ) − Sc(λ) Qac(λ)Sc(λ))dλ, (3.35) 2π R Z ph 1 c c ∗ c Jρ,Q := tr(ρac(λ)(Qac(λ) − Sph(λ) Qac(λ)Sph(λ)))dλ, (3.36) 2π R take place. 30 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
Proof. Since Sc and Sph are L1-scattering systems from Theorem 3.1 the representations (3.35) and (3.36) are easily follow. Taking into account (3.33), we get the following: c c ∗ c tr(ρac(λ)(Qac(λ) − Sph(λ) Qac(λ)Sph(λ))) = ∗ ∗ c tr(W−(λ)ρacW−(λ) (W+(λ)Qac(λ)W+(λ) − Sph(λ) Qac(λ)Sph(λ))). ∗ Using Sc(λ) = W+(λ) W−(λ) we find that: c c ∗ c tr(ρac(λ)(Qac(λ) − Sph(λ) Qac(λ)Sph(λ))) = tr (ρac(λ)× (3.37) ∗ ∗ ∗ ∗ (Sc(λ) Qac(λ)Sc(λ) − W−(λ) Sph(λ) W+(λ)Qac(λ)W+(λ) Sph(λ)W−(λ))) .
Since {Hc,H0} and {H,Hc} are L1-scattering systems, the existence of the wave opera- tors W±(H,Hc) and W±(Hc,H0) follows. Using the chain rule, we find W±(H,H0) = W±(H,Hc)W±(Hc,H0) which yields: ∗ S = W+(H,H0) W+(H,H0) ∗ ∗ = W+(Hc,H0) W+(H,Hc)W−(H,Hc)W−(Hc,H0) = W+(Hc,H0) SphW−(Hc,H0).
Hence, the scattering matrix {S(λ)}λ∈R of {H,H0} admits the representation ∗ S(λ) = W+(λ) Sph(λ)W−(λ), λ ∈ R. (3.38) Inserting (3.38) into (3.37), we get the following: Z ph 1 ∗ ∗ Jρ,Q = tr(ρac(λ)(Sc(λ) Qac(λ)Sc(λ) − S(λ) Qac(λ)S(λ)))dλ (3.39) 2π R Using (3.39), we obtain the following: Z c ph 1 ∗ Jρ,Q + Jρ,Q = tr(ρac(λ)(Qac(λ) − S(λ) Qac(λ)S(λ)))dλ. 2π R Finally, taking into account (3.4), we obtain (3.34). Remark 3.5. c (i) The current Jρ,Q is due to the coupling of the leads to the quantum dot and it is therefore called the contact induced current. ph (ii) The current Jρ,Q is due to the interaction of photons with electrons and it is called the photon induced current. Notice the this current is calculated under the assumption that the leads are already in contact with the dot. Corollary 3.6. Let the assumptions of Proposition 3.4 be satisfied. With respect to the spectral ac ac ph representation Π(H0 ) of H0 the photon induced current Jρ,Q can be represented by: Z ph 1 ∗ Jρ,Q := tr( ρbac(λ)(Qac(λ) − Sbph(λ) Qac(λ) Sbph (λ)))dλ, (3.40) 2π R where the measurable families { Sbph(λ) }λ∈R and { ρbac(λ) }λ∈R are given by (3.28) and ∗ ρbac(λ) := Sc(λ)ρac(λ)Sc(λ) λ ∈ R, (3.41) respectively.
∗ Proof. Using (3.33) and Sc(λ) = W+(λ) W−(λ), we find: c c ∗ c tr(ρac(λ)(Qac(λ) − Sph(λ) Qac(λ)Sph(λ))) = ∗ ∗ ∗ ∗ tr (Sc(λ)ρac(λ)Sc(λ) (Qac(λ) − W+(λ) Sph(λ) W+(λ)Qac(λ)W+(λ) Sph(λ)W+(λ))) . Quantum dot light emitting-absorbing devices 31
Taking into account the representations (3.28) and (3.41), we get the following: c c ∗ c tr(ρac(λ)(Qac(λ) − Sph(λ) Qac(λ)Sph(λ))) = ∗ ∗ tr(Sc(λ)ρac(λ)Sc(λ) (Qac(λ) − Sbph(λ) Qac(λ) Sbph(λ) )), which immediately yields (3.40).
Remark 3.7. In the following, we call { ρbac(λ) }λ∈R, cf. (3.41), the photon modified electron density matrix. Notice that { ρbac(λ) }λ∈R might be non-diagonal, even if the electron density matrix {ρac(λ)}λ∈R is diagonal. 4. Analysis of currents
c ph In the following, we analyze currents Jρ,Q and Jρ,Q under the assumption that ρ and Q have the tensor product structure: ρ = ρel ⊗ ρph and Q = qel ⊗ qph, (4.1) el ph el ph el where ρ and ρ as well as q and q are density operators and (tempered) charges for h0 and hph, respectively. Since ρph commutes with hph, which is discrete, the operator ρphhas the form: ph X ph ρ = ρ (n)(·, Υn)Υn, (4.2) n∈N0 where ρph(n) are non-negative numbers. Similarly, qph can be represented as: ph X ph q = q (n)(·, Υn)Υn, (4.3) n∈N0 where qph(n) are real numbers.
Lemma 4.1. Let S = {H,H0} be the JCL-model. Assume that ρ 6= 0 and Q have the structure el el el (4.1) where ρ is a density operator and q is a charge for h0 . p+2 (i) The operator (H0 − i) ρ, p ∈ N0, is bounded if and only if the condition: sup ρph(n)np+2 < ∞, (4.4) n∈N0 is satisfied. (ii) The charge Q is p-tempered if and only if: sup |qph(n)|n−p < ∞, (4.5) n∈N is valid
p+2 Proof. (i) The operator (H0 − i) ρ admits the representation: p+2 M ph el p+2 el (H0 − i) ρ = ρ (n)(h0 + nω − i) ρ . p∈N0 We have: p+2 ph el p+2 el k(H0 − i) ρkL(H) = sup ρ (n)k(h0 + nω − i) ρ kL(hel) (4.6) p∈N0
= sup ρph(n)np+2n−(p+2) (hel + nω − i)p+2ρel . 0 L(hel) p∈N0 −(p+2) el p+2 el p+2 el Since limn→∞ n (h + nω − i) ρ = ω kρ kL(hel), we obtain for sufficiently 0 L(hel) large n ∈ N0 that: p+2 ω el −(p+2) el p+2 el kρ k el n k(h + nω − i) ρ k el . 2 L(h ) 6 0 L(h ) 32 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
Using that and (4.6), we immediately obtain (4.4). Conversely, from (4.6) and (4.4), we obtain that p+2 (H0 − i) ρ is a bounded operator. (ii) As above, we have: −p M ph el Q(H0 − i) = q (n)q . n∈N0 Hence: −p ph el el −p kQ(H0 − i) kL(H) = sup |q (n)|kq (h0 + nω − i) kL(hel). n∈N0 p el −p −p el Since limn→∞ n k(h0 +nω −i) kL(hel) = ω kq kL(hel), we similarly obtain, as above, that (4.5) holds. The converse is obvious.
4.1. Contact induced current
Let us recall that Sc = {Hc,H0} is a L1-scattering system. Straightforward computation shows that: el el W±(Hc,H0) = w±(h , h0 ) ⊗ Ihph , el el where w±(h , h0 ) is given by (3.24). Hence: el el ∗ el el Sc = sc ⊗ Ihph , where sc := w+(hc , h0 ) w−(hc , h0 ).
Proposition 4.2. Let S = {H,H0} be the JCL-model. Assume that ρ and Q are given by (4.1) el el el ph ph ph where ρ and q are density and charge operators for h0 and ρ and q for h , respectively. If c for some p ∈ N0 the conditions (4.4) and (4.5) are satisfied, then the current Jρ,Q is well defined and admits the representation:
c sc X ph ph Jρ,Q = γJρel,qel , γ := q (n)ρ (n), (4.7) n∈N0
sc ph ph c sc where Jρel,qel is defined by (3.2). In particular, if tr(ρ ) = 1 and q = Ihph , then Jρ,Q = Jρel,qel .
p+2 Proof. First, we note that by lemma 4.1 the operator (H0 −i) ρ is bounded and Q is p-tempered. Sc Hence, the current Jρ,Q is correctly defined and the Landauer-Büttiker formula (3.4) is valid. ac With respect to the spectral representation Π(H0 ) of Lemma 2.12, the charge matrix el ph {Qac(λ)}λ∈R of Qac = qac ⊗ q admits the representation: M el ph Qac(λ) = qac(λ − nω)q (n), λ ∈ R. (4.8) n∈N0
Since Sc = sc ⊗ Ihph , the scattering matrix {Sc(λ)}λ∈R admits the representation: M Sc(λ) = sc(λ − nω), λ ∈ R. n∈N0 Hence: ∗ Qac(λ) − Sc(λ) Qac(λ)Sc(λ) = (4.9) M ph el ∗ el q (n) qac(λ − nω) − sc(λ − nω) qac(λ − ωn)sc(λ − nω) . n∈N0
Moreover, the density matrix {ρac(λ)}λ∈R admits the representation: M ph el ρac(λ) = ρ (n)ρac(λ − nω) (4.10) n∈N0 Quantum dot light emitting-absorbing devices 33
Inserting (4.10) into (4.9) we find the following: ac ∗ ρ (λ)(Qac(λ) − Sc(λ) Qac(λ)Sc(λ)) = M ph ph el el ∗ el q (n)ρ (n)ρac(λ − nω) qac(λ − ωn) − sc(λ − nω) qac(λ − ωn)sc(λ − nω) n∈N0 P ph ph Since γ = n∈N0 q (n)ρ (n) is absolutely convergent by (4.4) and (4.5), we obtain that: ac ∗ tr (ρ (λ)(Qac(λ) − Sc(λ) Qac(λ)Sc(λ))) = (4.11) X ph ph el el ∗ el q (n)ρ (n)tr ρac(λ − nω) qac(λ − ωn) − sc(λ − nω) qac(λ − ωn)sc(λ − nω) n∈N0 Clearly, we have: el el ∗ el tr ρac(λ − nω) qac(λ − ωn) − sc(λ − nω) qac(λ − ωn)sc(λ − nω) 6 el el 4kρac(λ − nω)kL(hn(λ))kqac(λ − nω)kL(hn(λ)), λ ∈ R. We insert (4.11) into the Landauer-Büttiker formula (3.35). Using (4.4) and (4.5) as well as Z el el kρac(λ)kL(hn(λ))kqac(λ)kL(hn(λ))dλ < ∞, R we see that we can interchange the integral and the sum. By doing so, we get: Z c X ph ph 1 el Jρ,Q = q (n)ρ (n) tr ρac(λ − nω)× 2π R n∈N0 el ∗ el qac(λ − ωn) − sc(λ − nω) qac(λ − ωn)sc(λ − nω) dλ. Using (3.25) we prove (4.7). ph P ph ph ph ph If tr(ρ ) = 1, then N0 ρ (n) = 1. Furthermore, if ρ = Ih , then q (n) = 1. Hence, γ = 1. 4.2. Photon induced current ph To calculate the current Jρ,Q, we use representation (3.40). We then set:
ph Sbmn (λ) := Pm(λ) Sbph (λ) hn(λ), λ ∈ R, where { Sbph (λ)}λ∈R is defined by (3.28) and Pm(λ) is the orthogonal projection from h(λ), cf. el (2.37), onto hm(λ) := h (λ − mω), λ ∈ R.
Proposition 4.3. Let S = {H,H0} be the JCL-model. Assume that ρ and Q are given by (4.1) el el el ph ph ph where ρ and q are density and charge operators for h0 and ρ and q for h , respectively. If ph for some p ∈ N0 the conditions (4.4) and (4.5) are satisfied, then the current Jρ,Q is well-defined and it admits the following representation: Z ph X ph X ph 1 el Jρ,Q = ρ (m) q (n) dλ tr ρbac(λ − mω) × (4.12) 2π R m∈N0 n∈N0 el ph ∗ el ph qac(λ − nω)δmn − Sbnm (λ) qac(λ − nω) Sbnm (λ) ,
el where { ρbac(λ) }λ∈R is the photon modified electron density defined, cf. (3.41), which takes the following form: el el ∗ ρbac(λ) = sc(λ)ρ (λ)sc(λ) , λ ∈ R. (4.13) 34 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
p Proof. By Lemma 4.1 we get that that the charge Q is p-tempered and (H0 − i) ρ is a bounded ph Sph operator. By Corollary 3.2, the current Jρ,Q := Jρc,Qc is well-defined. ∗ Since Qac(λ) − Sbph(λ) Qac(λ) Sbph (λ) is a trace class operator for λ ∈ R, we get from (3.40) and (4.10) that:
∗ tr ρbac(λ) Qac(λ) − Sbph(λ) Qac(λ) Sbph (λ) = X ph el ∗ ρ (m)tr ρb (λ − mω) Pm(λ) Qac(λ) − Sbph(λ) Qac(λ) Sbph (λ) Pm(λ) m∈N0 Furthermore, we have:
∗ Pm(λ) Qac(λ) − Sbph(λ) Qac(λ) Sbph (λ) Pm(λ) ph el ∗ = q (m) q (λ − mω) − Pm(λ) Sbph(λ) Qac(λ) Sbph (λ) Pm(λ) ph el X ph ph ∗ el ph = q (m)q (λ − mω) − q (n) Sbnm(λ) q (λ − nω) Sbnm(λ) , n∈N0
bph ∗ b P for λ ∈ R where Snm(λ) := Pn(λ) Sph (λ)Pm(λ), λ ∈ R. Notice that n∈N0 is a sum with a finite number of summands. Hence: ∗ X ph X ph tr ρbac(λ) Qac(λ) − Sbph(λ) Qac(λ) Sbph (λ) = ρ (m) q (n)× m∈N0 n∈N0 el el ph ∗ el ph tr ρb (λ − mω) q (λ − mω)δmn − Sbnm(λ) q (λ − nω) Sbnm(λ) We are going to show that Z X ph X ph el ρ (m) |q (n)| tr ρb (λ − mω) × R m∈N0 n∈N0 el ph ∗ el ph q (λ − mω)δmn − Sbnm(λ) q (λ − nω) Sbnm(λ) dλ < ∞. Clearly, one has the following estimate: el el ph ∗ el ph |tr ρb (λ − mω) q (λ − mω)δmn − Sbnm(λ) q (λ − nω) Sbnm(λ) 6 el el el 2k ρb (λ − mω) kL(hm(λ)) kq (λ − mω)kL(hm(λ))δnm + kq (λ − nω)kL(hn(λ)) . Furthermore, we get: Z el el k ρb (λ − mω) kL(hm(λ))kq (λ − mω)kL(hm(λ))δnm 6 λ∈R Z el el k ρb (λ) kL(hm(λ))kq (λ)kL(hm(λ))dλ λ∈R and Z el el k ρb (λ − mω) kL(hm(λ))kq (λ − nω)kL(hn(λ))dλ 6 R Z el el kqackL(hel) k ρb (λ − (m − n)ω) kL(hm−n(λ))dλ λ∈R If the conditions (4.4) and (4.5) are satisfied, then Z X ph ph el el ρ (m)|q (m)| k ρb (λ) kL(hm(λ))kq (λ)kL(hm(λ))dλ < ∞. R m∈N0 Quantum dot light emitting-absorbing devices 35
Furthermore, we obtain: Z X ph X ph el ρ (m) |q (n)| k ρb (λ − (m − n)ω) kL(hm−n(λ))dλ 6 λ∈R m∈N0 n∈N0 el X ph X ph (vmax − vmin + 4)kρackL(hel) ρ (m) |q (n)| < ∞, m∈N0 |m−n|6dmax where dmax is introduced by Lemma 2.12. To prove the following: X ρph(m) X |qph(n)| < ∞, m∈N0 |m−n|6dmax we again use (4.4) and (4.5). The last step allows us to interchange the integral and the sums, which immediately proves (4.12).
Corollary 4.4. Let S = {H,H0} be the JCL-model. We assume that ρ and Q are given by (4.1), el el el ph ph ph where ρ and q are density and charge operators for h0 and ρ and q for h , respectively. If el el el el ρ is an equilibrium state, i.e. ρ = f (h0 ), then: Z ph X ph 1 ph el ph el Jρ,Q = q (n) ρ (n)f (λ − nω) − ρ (m)f (λ − mω) × 2π R m,n∈N0 ph ∗ el ph tr Sbnm (λ) qac(λ − nω) Sbnm (λ) dλ. (4.14) Proof. From (4.12), we obtain the following: Z ph X ph X ph 1 el Jρ,Q = q (n) ρ (m) dλ f (λ − mω)× 2π R n∈N0 m∈N0 el ph ∗ el ph tr qac(λ − nω)δmn − Sbnm (λ) qac(λ − nω) Sbnm (λ) . Hence, Z ph X ph 1 X ph el Jρ,Q = q (n) dλ ρ (m)f (λ − mω)× 2π R n∈N0 m∈N0 el ph ∗ el ph tr qac(λ − nω)δmn − Sbnm (λ) qac(λ − nω) Sbnm (λ) . This gives the following: Z ph X ph 1 ph el el Jρ,Q = q (n) dλ ρ (n)f (λ − nω)tr qac(λ − nω) − (4.15) 2π R n∈N0 X ph el ph ∗ el ph ρ (m)f (λ − mω)tr Sbnm (λ) qac(λ − nω) Sbnm (λ) . m∈N0 Since X ph el ph ∗ el ph ρ (m)f (λ − mω)tr Sbnm (λ) qac(λ − nω) Sbnm (λ) = m∈N0 X ph el ph el ph ∗ el ph ρ (m)f (λ − mω) − ρ (n)f (λ − nω) tr Sbnm (λ) qac(λ − nω) Sbnm (λ) + m∈N0 ph el X ph ∗ el ph ρ (n)f (λ − nω) tr Sbnm (λ) qac(λ − nω) Sbnm (λ) , m∈N0 then inserting this into (4.15) we obtain (4.14). 36 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
5. Electron and photon currents 5.1. Electron current To calculate the electron current induced by contacts and photon contact, we make the following choice throughout this section. We set el el ph el el ph Qα := qα ⊗ q , qα := −epα and q := Ihph , α ∈ {l, r}, (5.1) el el el where pα denotes the orthogonal projection from h onto hα . By e > 0, we denote the magnitude el el el of the elementary charge. Since pα commutes with hα , one can easily verify that Qα commutes el el with H0, which shows that Qα is a charge. Following [1], the flux related to Qα gives us the S el electron current J el entering the lead α from the sample. Notice Q = −ePα where Pα is ρ,Qα α el ph ph the orthogonal projection from H onto Hα := hα ⊗ h . Since q = Ihph , condition (4.5) is immediately satisfied for any p > 0. Let f(·): R −→ R be a non-negative bounded measurable function. We set: el el el el el el ρ = ρl ⊕ ρS ⊕ ρr , ρα := f(hα − µα), α ∈ {l, r}, (5.2) el ph and ρ = ρ ⊗ ρ . By µα the chemical potential of the lead α is denoted. In applications, one sets f(λ) := fFD(λ), λ ∈ R, where fFD(λ) is the so-called Fermi-Dirac distribution given by (3.21). el el ph If β = ∞, then fFD(λ) := χR− (λ), λ ∈ R. Notice that [ρ , p ] = 0. For ρ , we choose the Gibbs state: 1 ph ph 1 ρph := e−βh ,Z = tr(e−βh ) = . (5.3) Z 1 − e−βω ph −βω −βhph ph ph Hence, ρ = (1 − e )e . If β = ∞, then ρ := (·, Υ0)Υ0. Clearly, tr(ρ ) = 1. We note ph −βω −nβω that ρ (n) = (1 − e )e , n ∈ N0, satisfies the condition (4.4) for any p > 0. Accordingly, el ph ρ0 = ρ ⊗ ρ is the density operator for H0. el el Definition 5.1. Let S = {H,H0} be the JCL-model. If Q := Qα , where Qα is given by (5.1), el ph el ph el S and ρ := ρ0 := ρ ⊗ ρ , where ρ and ρ are given by (5.2) and (5.3), then J el := J el ρ0,Qα ρ0,Qα c ph is called the electron current entering the lead α. The currents J el and J el are called the ρ0,Qα ρ0,Qα contact-induced and photon-induced electron currents. 5.1.1. Contact induced electron current. The following proposition immediately follows from Proposition 4.2.
Proposition 5.2. Let S = {H,H0} be the JCL-model. Then the contact induced electron current c c sc J el , α ∈ {l, r}, is given by J el = J el el . In particular, one has: ρ0,Qα ρ0,Qα ρ ,qα Z c e J el = − (f(λ − µ ) − f(λ − µ )σ (λ)dλ, α, ∈ {l, r}, α 6= , (5.4) ρ0,Qα α κ c κ κ 2π R where {σc(λ)}λ∈R is the channel cross-section from left to the right of the scattering system sc = el el {h , h0 }, cf. Example 3.3.
ph c sc Proof. Since tr(ρ ) = 1 it follows from Proposition 4.2 that J el = J el el . From (3.26), cf. ρ0,Qα ρ ,qα Example 3.3, we find (5.4). c If µl > µr and f(·) is decreasing, then J el < 0. Hence, the electron contact current ρ0,Ql is going from the left lead to the right which is in accordance with the physical expectations. In particular, this is valid for the Fermi-Dirac distribution.
el ph Proposition 5.3. Let S = {H,H0} be the JCL-model. Further, let ρ and ρ be given by (5.2) el and (5.3), respectively. If the charge Qα is given by (5.1), then the following holds: Quantum dot light emitting-absorbing devices 37
c (E) If µl = µr, then J el = 0, α ∈ {l, r}. ρ0,Qα c (S) If vl vr + 4, then J el = 0, α ∈ {l, r}, even if µl 6= µl. > ρ0,Qα S S S S c (C) If e = δ and e = δ , then J el = 0, α ∈ {l, r}, even if µl 6= µl. 0 0 1 1 ρ0,Qα c Proof. (E) If µl = µr, then f(λ − µl) = f(λ − µr). Applying formula (5.4) we obtain J el = 0. ρ0,Qα el,ac (S) If vl > vr + 4, then h0 has simple spectrum. Hence the scattering matrix {sc(λ)}λ∈R el el of the scattering system sc = {h , h0 } is a scalar function which immediately yields σc(λ) = 0, c λ ∈ , which yields J el = 0. R ρ0,Qα (C) In this case, the Hamiltonian hel decomposes into the direct sum of two non-interacting el el Hamiltonians. Hence, the scattering matrix of {sc(λ)}λ∈R of the scattering system sc = {h , h0 } c is diagonal, which immediately yields J el = 0. ρ0,Qα
5.1.2. Photon induced electron current. It is hopeless to analyze the properties of (4.12) if we el make no assumptions concerning ρ and the scattering operator sc. The simplest assumptions is el el el that ρ and sc commute. In this case, we get ρb (λ) = ρ (λ), λ ∈ R. el Lemma 5.4. Let S = {H,H0} be the JCL-model. Furthermore, let ρ be given by (5.2). If one el of the cases (E), (S) or (C) of Proposition 5.3 is realized, then the ρ and sc commute.
el el el Proof. If (E) holds, then ρ = f(h0 ) which yields [ρ , sc] = 0. If (S) is valid, then the scattering el matrix {sc(λ)}λ∈R is a scalar function which shows [ρ , sc] = 0. Finally, if (C) is realized, then el el the scattering matrix {sc(λ)}λ∈R diagonal. Since the ρ is given by (5.2) we get [ρ , sc] = 0.
ph We will now calculate the current J el , see (4.12). Clearly, we have Pα(λ) = ρ0,Qα P el n∈N0 pα (λ − nω) and Ih(λ) = Pl(λ) + Pr(λ), λ ∈ R. We then set: el Pnα (λ) := Pα(λ)Pn(λ) = Pn(λ)Pα(λ) = pα (λ − nω), α ∈ {l, r}, n ∈ N0, λ ∈ R. In the following we use the notation Tbph (λ) = Sbph (λ) − Ih(λ), λ ∈ R, where
{ Tbph(λ) }λ∈R is called the transition matrix and { Sbph(λ) }λ∈R is given by (3.28). We set: ph Tb (λ) := Pkα (λ) Tbph (λ)Pm (λ), λ ∈ , α, ∈ {l, r}, k, m ∈ 0, kαmκ κ R κ N and: ph ph ph σ (λ) = tr( Tb (λ)∗ Tb (λ)), λ ∈ , (5.5) bkαmκ kαmκ kαmκ R which is the cross-section between the channels kα and mκ.
Proposition 5.5. Let S = {H,H0} be the JCL-model. el el (i) If ρ commutes with the scattering operator sc and q , then: Z ph X e ph ph ph J el = − ρ (n)f(λ − µ − nω) − ρ (m)f(λ − µ − mω) σ (λ) dλ. ρ0,Q α κ bnαm α 2π R κ m,n∈N0 κ∈{l,r} (5.6) (ii) If in addition S = {H,H0} is time reversible symmetric, then Z ph X e ph ph ph J el = − ρ (n)f(λ − µα − nω) − ρ (m)f(λ − µα0 − mω) σ (λ) dλ, ρ0,Q bnαmα0 α 2π R m,n∈N0 (5.7) α, α0 ∈ {l, r}, α 6= α0. 38 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
Proof. (i) Let us assume that: el X el q = gκ(hκ ), κ∈{l,r} Notice that el X el qac(λ) = gκ(λ)pκ (λ), λ ∈ R. (5.8) κ∈{l,r} ph Inserting (5.8) into (4.12) and using q = Ihph , we obtain the following: Z ph X ph X 1 Jρ0,Q = ρ (m) dλ φα(λ − mω)gκ(λ − nω)× 2π R m∈N0 n∈N0 α∈{l,r} κ∈{l,r} el el bph ∗ el bph tr pα (λ − mω) pκ (λ − nω)δmn − Snm (λ) pκ (λ − nω) Snm (λ) , where, for simplicity, we have set:
φα(λ) := f(λ − µα), λ ∈ R, n ∈ N0, α ∈ {l, r}. (5.9) Accordingly, we get: 1 Z J ph = X ρph(n) dλ φ (λ − nω)g (λ − nω)tr pel(λ − nω) − ρ0,Q κ κ κ 2π R n∈N0 κ∈{l,r} Z X X ph 1 ρ (m) dλ φα(λ − mω)gκ(λ − nω)× (5.10) 2π R n∈N0 m∈N0 κ∈{l,r} α∈{l,r} el bph ∗ el bph el tr pα (λ − mω) Snm (λ) pκ (λ − nω) Snm (λ)pα (λ − mω) . ph Since the scattering matrix { Sb (λ)}λ∈R is unitary, we have: X el el bph ∗ el bph el pκ (λ − nω) = pκ (λ − nω) Smn (λ) pα (λ − mω) Smn (λ)pκ (λ − nω), (5.11) m∈N0 α∈{l,r} for n ∈ N0 and κ ∈ {l, r}. Inserting (5.11) into (5.10), we find that: Z ph X X ph 1 Jρ0,Q = ρ (n) dλ φκ(λ − nω)gκ(λ − nω)× 2π R n∈N0 m∈N0 κ∈{l,r} α∈{l,r} el bph ∗ el bph el tr pκ (λ − nω) Snm (λ) pα (λ − mω) Smn (λ)pκ (λ − nω) − Z X X ph 1 ρ (m) dλ φα(λ − mω)gκ(λ − nω)× 2π R n∈N0 m∈N0 κ∈{l,r} α∈{l,r} el bph ∗ el bph el tr pα (λ − mω) Snm (λ) pκ (λ − nω) Snm (λ)pα (λ − mω) . Using the notation (5.5), we find the following: Z ph X X ph 1 ph Jρ0,Q = ρ (n) dλ φκ(λ − nω)gκ(λ − nω) σbmαn (λ) − 2π R κ n∈N0 m∈N0 κ∈{l,r} α∈{l,r} Z X X ph 1 ph ρ (m) dλ φα(λ − mω)gκ(λ − nω) σbn mα (λ) . 2π R κ n∈N0 m∈N0 κ∈{l,r} α∈{l,r} Quantum dot light emitting-absorbing devices 39
By (3.16), we find that: X σph (λ) = X σph (λ) λ ∈ . bmαnκ bnκ mα R m∈N0 m∈N0 α∈{l,r} α∈{l,r} Using that, we obtain the following: Z ph X 1 Jρ0,Q = × (5.12) 2π R m,n∈N0 α,κ∈{l,r} ph ph ph ρ (n)φ (λ − nω) − ρ (m)φα(λ − mω) g (λ − nω) σ (λ) dλ. κ κ bnκ mα
Setting gα(λ) = −e and gκ(λ) ≡ 0, κ 6= α, we obtain (5.6). (ii) Straightforward computation shows that: Z X ph ph ph ρ (n)f(λ − µa − nω) − ρ (m)f(λ − µa − mω) σbnαmα (λ) dλ = R n,m∈N0 Z X ph ph ph ρ (m)f(λ − µa − mω) − ρ (n)f(λ − µa − nω) σbmαnα (λ) dλ. R n,m∈N0 ph ph Since σmαnα (λ) = σnαmα (λ), λ ∈ R, we obtain the following: Z X ph ph ph ρ (n)f(λ − µa − nω) − ρ (m)f(λ − µa − mω) σbnαmα (λ) dλ = R n,m∈N0 Z X ph ph ph − ρ (n)f(λ − µa − nω) − ρ (m)f(λ − µa − mω) σbnαmα (λ) dλ R n,m∈N0 which yields: Z X ph ph ph ρ (n)f(λ − µa − nω) − ρ (m)f(λ − µa − mω) σbnαmα (λ) dλ = 0. R n,m∈N0 Using that, we immediately obtain the representation (5.7) from (5.6).
Corollary 5.6. Let S = {H,H0} be theJCL-model. (i) If the cases cases (E), (S) or (C) of Proposition 5.3 are realized, then the representation (5.6) holds. (ii) If the case (E) of Proposition 5.3 is realized and the system S = {H,H0} is time reversible symmetric, then Z ph X e ph ph ph J el = − (ρ (n)f(λ − µ − nω) − ρ (m)f(λ − µ − mω)) σ (λ)dλ (5.13) ρ0,Q bnαmα0 α 2π R m,n∈N0 0 n ∈ N0, α ∈ {l, r} where µ := µl = µr and α 6= α . (iii) If the case (E) of Proposition 5.3 is realized and the system S = {H,H0} is time reversible ph and mirror symmetric, then J el = 0. ρ0,Qα Proof. (i) The statement follows from Proposition 5.5(i) and Lemma 5.4. (ii) Setting µα = µα0 formula (5.13) follows (5.7). (iii) If S = {H,H0} is time reversible and mirror symmetric, we obtain from Lemma ph ph 0 0 2.14 (ii) that σ (λ) = σ (λ), λ ∈ , n, m ∈ 0, α, α ∈ {l, r}, α 6= α . Using that, we bnαmα0 bnα0 mα R N obtain from (5.13) the following: Z ph X e ph ph ph J el = − (ρ (n)f(λ − µ − nω) − ρ (m)f(λ − µ − mω)) σ (λ)dλ. ρ0,Q bnα0 mα α 2π R m,n∈N0 40 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
Interchanging m and n, we obtain: Z ph X e ph ph ph J el = − (ρ (m)f(λ − µ − mω) − ρ (n)f(λ − µ − nω)) σ (λ)dλ. ρ0,Q bmα0 nα α 2π R m,n∈N0 Using that S is time reversible symmetric we get from Lemma 2.14 (i) that: Z ph X e ph ph ph J el = − (ρ (m)f(λ − µ − mω) − ρ (n)f(λ − µ − nω)) σ (λ)dλ. ρ0,Q bnαmα0 α 2π R m,n∈N0 ph ph ph which shows that J el = −J el . Hence J el = 0. ρ0,Qα ρ0,Qα ρ0,Qα c We note that by Proposition 5.3, the contact induced current is zero, i.e. J el = 0. Hence, ρ0,Qα S if the S is time reversible and mirror symmetric, then the total current is zero, i.e. J el = 0. ρ0,Qα
Remark 5.7. Let the case (E) of Proposition 5.3 be realized, that is, µl = µr. Moreover, we assume for simplicity that 0 =: vr 6 v := vl. ph ph (i) If β = ∞, then ρ (n) = δ0n, n ∈ 0. From (5.6), we immediately obtain that J el el = 0. N ρ ,Qα That means, if the temperature is zero, then the photon-induced electron current is zero. (ii) The photon-induced electron current might be zero even if β < ∞. Indeed, let S = {H,H0} be time reversible symmetric and let the case (E) be realized. If ω > v + 4 and el el el h (λ) := hn (λ) = h (λ − nω), n ∈ N0. Hence, one always has n = m in formula (5.13), ph which immediately yields J el = 0. ρ0,Qα (iii) The photon-induced electron current might be different than zero. Indeed, let S = {H,H0} be time reversible symmetric and let v = 2 and ω = 4, then one sees that to calculate the ph J el , one has to consider m = n + 1 in formula (5.13). Therefore, we find that: ρ0,Ql
ph X e Jρ ,Qel = − × 0 l 2π n∈N0 Z dλ ρph(n)f(λ − µ − nω) − ρph(n + 1)f(λ − µ − (n + 1)ω) σph (λ). bnl (n+1)r R ph If ρ is given by (5.3) and f(λ) = fFD(λ), cf. (3.21), then one easily verifies that ∂ ρph(x)f (λ − µ − xω) < 0, x, µ, λ ∈ . ∂x FD R ph Hence, ρ (n)fFD(λ − µ − nω) is decreasing in n ∈ N0 for λ, µ ∈ R which yields ph ph ph ρ (n)f(λ − µ − nω) − ρ (n + 1)f(λ − µ − (n + 1)ω) > 0. Therefore, J el 6 0 ρ0,Ql which means that the photon-induced current leaves the left-hand side and enters the right- ph ph hand side. In fact, J el = 0 implies that σbn (n+1) (λ) = 0 for n ∈ N0 and λ ∈ R, which ρ0,Ql l r means that there is no scattering from the left-hand side to the right one and vice versa which can be excluded generically. 5.2. Photon current The photon current is related to the charge by equation: ph Q := Q = −Ihel ⊗ n, ∗ ph where n = dΓ(1) = b b is the photon number operator on h = F+(C), which is self-adjoint ph ph and commutes with h . It follows that Q is also self-adjoint and commutes with H0. It is not ph ph −1 bounded, but since dom(n) = dom(h ), it is immediately obvious that Q (H0 +θ) is bounded, Quantum dot light emitting-absorbing devices 41
ac when N is a tempered charge. Its charge matrix with respect to the spectral representation Π(H0 ) of Lemma 2.12 is given by: ph M Qac (λ) = − nPn(λ). n∈N0 el We recall that Pn(λ) is the orthogonal projection from h(λ) onto hn(λ) = h (λ − nω), λ ∈ R. We will now calculate the photon current or, as it is also known, the photon production rate.
5.2.1. Contact induced photon current. The following proposition is, in fact, in accordance with the physical intuition.
c Proposition 5.8. Let = {H,H } be the JCL-model. Then J ph = 0. S 0 ρ0,Q
el sc Proof. We note that qac(λ) = Ihel(λ), λ ∈ R. Inserting this into (3.25), we obtain Jρel,qel = 0. c Applying Proposition 4.2 we prove J ph = 0. ρ0,Q The result reflects the fact that the lead contact does not contribute to the photon current, which is plausible from the physical point of view.
5.2.2. Photon current. From Proposition 5.8, we see that only the photon-induced photon current ph S S ph ph J ph contributes to the photon current J ph . Since J ph = J ph , we call J ph simply ρ0,Q ρ0,Q ρ0,Q ρ0,Q ρ0,Q the photon current. ph el Using the notation Tbnm (λ) := Pn(λ) Tbph (λ) h (λ − mω), λ ∈ R, m, n ∈ N0. We set: ph ph Tenm(λ) = Tbnm (λ)sc(λ − mω), λ ∈ R, m, n ∈ N0, (5.14) and ph ph el Te (λ) := Pn (λ)Te (λ) h (λ − mω), λ ∈ , (5.15) nκ mα κ nm α R ph ph ∗ ph m, n ∈ 0, α, ∈ {l, r}, as well as σ (λ) := tr(Te (λ) Te (λ)), λ ∈ . N κ enκ mα nκ mα nκ mα R
Proposition 5.9. Let S = {H,H0} be the JCL-model. (i) Then: Z ph X ph 1 ph J ph = (n − m)ρ (m) f(λ − µ − mω)σ (λ)dλ. (5.16) ρ0,Q α en mα 2π R κ m,n∈N0 α,κ∈{l,r}
el (ii) If ρ commutes with sc, then: Z ph X ph 1 ph J ph = (n − m)ρ (m) f(λ − µ − mω) σ (λ)dλ. (5.17) ρ0,Q α bn mα 2π R κ m,n∈N0 α,κ∈{l,r}
el (iii) If ρ commutes with sc and S = {H,H0} is time reversible symmetric, then: Z ph X 1 J ph = dλ× (5.18) ρ0,Q 2π R m,n∈N0,n>m κ,α∈{l,r} ph ph ph (n − m) ρ (m)f(λ − µα − mω) − ρ (n)f(λ − µ − nω) σ (λ), κ bnκ mα where α0 ∈ {l, r} and α0 6= α. 42 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov
Proof. (i) From (4.12) we get
ph X ph J ph = − nρ (m)× ρ0,Q m,n∈N0 Z 1 el ph ∗ el ph dλ tr ρbac(λ − mω) Pn(λ)δmn − Sbnm (λ) qac(λ − nω) Sbnm (λ) . 2π R Hence: Z ph X ph 1 el ph ∗ ph J ph = − mρ (m) tr ρ (λ − mω) P (λ) − Sb (λ) P (λ) Sb (λ) dλ+ ρ0,Q bac m mm m mm 2π R m∈N0 Z X ph 1 el ph ∗ ph nρ (m) tr ρbac(λ − mω) Sbnm(λ) Pn(λ) Sbnm(λ) dλ. 2π R m,n∈N0 m6=n P Using the relation Pm(λ) = Ih(λ) − n∈N0,m6=n Pn(λ), λ ∈ R, we obtain the following: Z ph X ph 1 el ph ∗ ph J ph = − mρ (m) tr ρ (λ − mω) Sb (λ) P (λ) Sb (λ) dλ+ ρ0,Q bac nm n nm 2π R m,n∈N0 m6=n Z X ph 1 el ph ∗ ph nρ (m) tr ρbac(λ − mω) Sbnm(λ) Pn(λ) Sbnm(λ) dλ. 2π R m,n∈N0 m6=n
Since Tbph (λ) = Sbph (λ) − Ih(λ), λ ∈ R, we find Z ph X ph 1 el ph ∗ ph J ph = − (m − n)ρ (m) tr ρ (λ − mω) Tb (λ) Tb (λ) dλ. ρ0,Q bac nm nm 2π R m,n∈N0 Using (4.13) and definition (5.14) one readily sees that: Z ph X ph 1 el ph ∗ ph J ph = − (m − n)ρ (m) tr ρ (λ − mω)Te (λ) Te (λ) dλ . ρ0,Q ac nm nm 2π R m,n∈N0
el el el el Since ρac = ρl ⊕ ρr where ρα is given by (5.2), we find the following: Z ph X ph 1 ph ∗ ph J ph = − (m − n)ρ (m) f(λ − µ − mω)tr Te (λ) Te (λ) dλ, ρ0,Q α n mα n mα 2π R κ κ m,n∈N0 α,κ∈{l,r} where we have used (5.15). Using σph (λ) = tr(Teph (λ)∗Teph (λ)), we prove (5.16). enκ mα nκ mα nκ mα el el el (ii) If ρac commutes with sc, then ρbac (λ) = ρac(λ), λ ∈ R, which yields that one can replace σph (λ) by σph (λ), λ ∈ . Therefore, (5.17) holds. enκ mα bnκ mα R (iii) Clearly, we have:
ph J ph = (5.19) ρ0,Q Z X ph 1 ph (n − m)ρ (m) f(λ − µα − mω) σbn mα (λ)dλ + 2π R κ m,n∈N0,n>m α,κ∈{l,r} Z X ph 1 ph (n − m)ρ (m) f(λ − µα − mω) σbn mα (λ)dλ . 2π R κ m,n∈N0,n Moreover, a straightforward computation shows the following: Z X ph 1 ph (n − m)ρ (m) f(λ − µα − mω) σbn mα (λ)dλ = 2π R κ m,n∈N0,n Since S = {H,H0} is time reversible symmetric, we find the following: Z X ph 1 ph (n − m)ρ (m) f(λ − µα − mω) σbmαn (λ)dλ = (5.20) 2π R κ m,n∈N0,n Corollary 5.10. Let S = {H,H0} be the JCL-model and let f = fFD. If case (E) of Proposition ph 5.3 is realized and = {H,H } is time reversible symmetric, then J ph 0. S 0 ρ0,Q > Proof. We set µ := µl = µr. One has: ρph(m)f(λ − µ − mω) − ρph(n)f(λ − µ − nω) = −mβω −(n−m)βω e (1 − e )fFD(λ − µ − mω)fFD(λ − µ − nω) > 0, ph for n > m. From (5.18), we see that J ph 0. ρ0,Q > ph Remark 5.11. We will now comment the results. If J ph 0, then system is called light ρ0,Q > S ph emitting. Similarly, if J ph 0, then we call it light absorbing. Of course if is light emitting ρ0,Q 6 S ph and absorbing, then J ph = 0. ρ0,Q ph (i) If β = ∞, then ρ (m) = δ0m, m ∈ N0. Inserting this into (5.16), we get: Z ph X 1 ph J ph = n f(λ − µα)σ (λ)dλ 0 ρ0,Q enκ 0α > 2π R n∈N0 α,κ∈{l,r} Hence, the system S is light emitting. (ii) Let us show S might be light emitting even if β < ∞. We consider the case (E) of Proposition 5.3. If S is time reversible symmetric, then it follows from Corollary 5.10 that the system is light emitting. ph If the system S is time reversible and mirror symmetric, then J el = 0, α ∈ {l, r}, ρ0,Qα c S by Corollary 5.6(iii) . Since J el = 0 by Proposition 5.3, we get that J el = 0 but ρ0,Q ρ0,Qα the photon current is larger than zero. So our JCL-model is light emitting by a zero total S electron current J el . ρ0,Qα Let vr = 0, vl = 2 and ω = 4. Hence S is not mirror symmetric. Then, we get from ph ph Remark 5.7(iii) that J el = −Jρ ,Qel 6 0. Hence, there is an electron current from the ρ0,Ql 0 r c S left to the right lead. Notice that by Proposition 5.3 J el = 0. Hence, J el 6 0. ρ0,Ql ρ0,Ql 44 H. Neidhardt, L. Wilhelm, V.A. Zagrebnov (iii) To realize a light absorbing situation, we consider the case (S) of Proposition 5.3 and assume that S is time reversible symmetric. Notice that by Lemma 5.4, sc commutes with ρel. We make the following choices: vr = 0, vl > 4, ω = vl, µl = 0, µr = ω = vl. It follows out that with respect to the representation (5.18) one has only to m = n − 1, κ = r and α = l. Hence, ph X 1 J ph = × ρ0,Q 2π n∈N Z dλ ρph(n − 1)f(λ − (n − 1)ω) − ρph(n)f(λ − (n + 1)ω) σph (λ) bnl(n−1)r R Since, f(λ) = fFD(λ), we find: ρph(n − 1)f(λ − (n − 1)ω) − ρph(n)f(λ − (n + 1)ω) = ρph(n − 1)f(λ − (n − 1)ω)f(λ − (n + 1)ω) × 1 + eβ(λ−(n+1)ω) − e−βω(1 + eβ(λ−ω(n−1))) , or ρph(n − 1)f(λ − (n − 1)ω) − ρph(n)f(λ − (n + 1)ω) = ρph(n − 1)f(λ − (n − 1)ω)f(λ − (n + 1)ω)(1 − e−βω)(1 − eβ(λ−ωn)). ph ph Since λ − nω > 0 we find ρ (n − 1)f(λ − (n − 1)ω) − ρ (n)f(λ − (n + 1)ω) 6 0 which ph yields J ph 0. ρ0,Q 6 ph 0 To calculate J el , we use formula (5.7). Setting α = l, we obtain α = r, which yields ρ0,Ql ph X e Jρ ,Qel = − × 0 l 2π m,n∈N0 Z dλ ρph(n)f(λ − µ − nω) − ρph(m)f(λ − µ − mω) σph (λ) . r l bnlmr R σph (λ) = 0 σph (λ) = 0 m 6= n + 1 n ∈ One verifies that b0l0r and bnlmr for , N. Hence, ph X e Jρ ,Qel = − × 0 l 2π n∈N Z dλ ρph(n)f(λ − µ − nω) − ρph(n − 1)f(λ − µ − (n + 1)ω) σph (λ) , r l bnl(n+1)r R Since µr = ω and µl = 0, we find: ph X e Jρ ,Qel = − × 0 l 2π n∈N Z f(λ − (n + 1)ω)ρph(n − 1)(1 − e−βω) σph (λ) dλ 0. bnl(n+1)r 6 R Hence, there is electron current flowing from the left to right induced by photons. We recall c that J el = 0. ρ0,Ql Quantum dot light emitting-absorbing devices 45 Acknowledgments The first two authors would like to thank the University of Aalborg and the Centre de Physique Théorique - Luminy for hospitality and financial support. We thank Horia D. Cornean for discus- sions concerning the JCL-model and Igor Yu.Popov for the opportunity to present this complete version of our results. The work on this paper has been partially supported by the ERC grant no.267802 "AnaMultiScale". References [1] W. Aschbacher, V. Jakšic,´ Y. Pautrat, and C.-A. Pillet. Transport properties of quasi-free fermions. J. Math. Phys., 48(3), P. 032101, 28 (2007). [2] S. Attal, A. Joye, and C.-A. Pillet, editors. Open quantum systems. I, volume 1880 of Lecture Notes in Mathe- matics. Springer-Verlag, Berlin (2006). [3] S. Attal, A. Joye, and C.-A. Pillet, editors. Open quantum systems. II, volume 1881 of Lecture Notes in Mathe- matics. Springer-Verlag, Berlin, 2006. The Markovian approach, Lecture notes from the Summer School held in Grenoble, June 16–July 4, 2003. [4] S. Attal, A. Joye, and C.-A. Pillet, editors. Open quantum systems. III, volume 1882 of Lecture Notes in Math- ematics. Springer-Verlag, Berlin, 2006. Recent developments, Lecture notes from the Summer School held in Grenoble, June 16–July 4, 2003. [5] M. Baro, H.-Chr. Kaiser, H. Neidhardt, and J. Rehberg. A quantum transmitting Schrödinger-Poisson system. Rev. Math. Phys., 16(3), P. 281–330 (2004). [6] H. Baumgärtel and M. Wollenberg. Mathematical scattering theory. Akademie-Verlag, Berlin (1983). [7] M. Büttiker, Y. Imry, R. Landauer, and S. Pinhas. Generalized many-channel conductance formula with applica- tion to small rings. Phys. Rev. B, 31(10), P. 6207–6215 (May 1985). [8] H. D. Cornean, P. Duclos, G. Nenciu, and R. Purice. Adiabatically switched-on electrical bias and the Landauer- Büttiker formula. J. Math. Phys., 49(10), P. 102106, 20 (2008). [9] H. D. Cornean, P. Duclos, and R. Purice. Adiabatic non-equilibrium steady states in the partition free approach. Ann. Henri Poincaré, 13(4), P. 827–856 (2012). [10] H. D. Cornean, C. Gianesello, and V. A. Zagrebnov. A partition-free approach to transient and steady-state charge currents. J. Phys. A, 43(47), P. 474011, 15 (2010). [11] H. D. Cornean, A. Jensen, and V. Moldoveanu. A rigorous proof of the Landauer-Büttiker formula. J. Math. Phys., 46(4), P. 042106, 28 (2005). [12] H. D. Cornean, A. Jensen, and V. Moldoveanu. The Landauer-Büttiker formula and resonant quantum transport. Lecture Notes in Phys., 690, pages 45–53. Springer, Berlin (2006). [13] H. D. Cornean, H. Neidhardt, L. Wilhelm, and V. A. Zagrebnov. The Cayley transform applied to non-interacting quantum transport. J. Funct. Anal., 266(3), P. 1421–1475 (2014). [14] H. D. Cornean, H. Neidhardt, and V. A. Zagrebnov. The effect of time-dependent coupling on non-equilibrium steady states. Ann. Henri Poincaré, 10(1), P. 61–93 (2009). [15] M. Damak. On the spectral theory of tensor product Hamiltonians. J. Operator Theory, 55(2), P. 253–268 (2006). [16] C. Gerry and P. Knight. Introductory Quantum Optics. Cambridge University Press, Cambridge (2005). [17] T. Kato. Perturbation theory for linear operators. Classics in Mathematics. Springer-Verlag, Berlin (1995). Reprint of the 1980 edition. [18] R. Landauer. Spatial Variation of Currents and Fields Due to Localized Scatterers in Metallic Conduction. IBM J. Res. Develop., 1(3), P. 223–231 (1957). [19] H. Neidhardt, L. Wilhelm, and V. A. Zagrebnov. A new model for quantum dot light emitting-absorbing devices. J. Math. Phys., Anal., Geom., 10(3), P. 350–385 (2014). [20] G. Nenciu. Independent electron model for open quantum systems: Landauer-Büttiker formula and strict posi- tivity of the entropy production. J. Math. Phys., 48(3), P. 033302, 8 (2007). NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS, 2015, 6 (1), P. 46–56 ON THE ROBIN EIGENVALUES OF THE LAPLACIAN IN THE EXTERIOR OF A CONVEX POLYGON Konstantin Pankrashkin Laboratoire de mathematiques,´ Universite´ Paris-Sud Batimentˆ 425, 91405 Orsay Cedex, France [email protected] PACS 41.20.Cv, 02.30.Jr, 02.30.Tb DOI 10.17586/2220-8054-2015-6-1-46-56 2 Let Ω ⊂ R be the exterior of a convex polygon whose side lengths are `1, . . . , `M . For a real constant α, let Ω Hα denote the Laplacian in Ω, u 7→ −∆u, with the Robin boundary conditions ∂u/∂ν = αu at ∂Ω, where Ω Ω ν is the outer unit normal. We show that, for any fixed m ∈ N, the mth eigenvalue Em(α) of Hα behaves Ω 2 D −1/2 D as Em(α) = −α + µm + O(α ) as α → +∞, where µm stands for the mth eigenvalue of the operator 00 D1 ⊕· · ·⊕DM and Dn denotes the one-dimensional Laplacian f 7→ −f on (0, `n) with the Dirichlet boundary conditions. Keywords: eigenvalue asymptotics, Laplacian, Robin boundary condition, Dirichlet boundary condition. Received: 5 November 2014 1. Introduction 1.1. Laplacian with Robin boundary conditions Let Ω ⊂ Rd, d ≥ 2, be a connected domain with a compact Lipschitz boundary ∂Ω. Ω For α > 0, let Hα denote the Laplacian u 7→ −∆u in Ω with the Robin boundary conditions Ω ∂u/∂ν = αu at ∂Ω, where ν stands for the outer unit normal. More precisely, Hα is the self-adjoint operator in L2(Ω) generated by the sesquilinear form: ZZ Z Ω 2 2 Ω 1,2 hα (u, u) = |∇u| dx − α |u| dσ, D(hα ) = W (Ω). Ω ∂Ω Here and below, σ denotes the (d − 1)-dimensional Hausdorff measure. Ω One checks in the standard way that the operator Hα is semibounded from below. If Ω is bounded (i.e. Ω is an interior domain), then it has a compact resolvent, and we denote Ω by Em(β), m ∈ N, its eigenvalues taken according to their multiplicities and enumerated in the non-decreasing order. If Ω is unbounded (i.e. Ω is an exterior domain), then the Ω essential spectrum of Hα coincides with [0, +∞), and the discrete spectrum consists of Ω finitely many eigenvalues which will be denoted again by Em(α), m ∈ {1,...,Kα}, and enumerated in the non-decreasing order taking into account the multiplicities. Ω We are interested in the behavior of the eigenvalues Em(α) for large α. It seems that the problem was introduced by Lacey, Ockedon, Sabina [11] when studying a reaction- diffusion system. Giorgi and Smits [6] studied a link to the theory of enhanced surface su- perconductivity. Recently, Freitas and Krejciˇ rˇ´ık [10] and then Pankrashkin and Popoff [15] studied the eigenvalue asymptotics in the context of the spectral optimization. Robin-Laplacian eigenvalues 47 Let us list some available results. Under various assumptions one showed the asymptotics of the form: Ω 2 2 Em(α) = −CΩα + o(α ) as α tends to +∞, (1) where CΩ ≥ 1 is a constant depending on the geometric properties of Ω. Lacey, Ockedon, 4 Sabina [11] showed (1) with m = 1 for C compact domains, for which CΩ = 1, and for triangles, for which CΩ = 2/(1 − cos θ), where θ is the smallest corner. Lu and Zhu [13] 1 showed (1) with m = 1 and CΩ = 1 for compact C smooth domains, and Daners and Kennedy [2] extended the result to any fixed m ∈ N. Levitin and Parnovski [12] showed (1) with m = 1 for domains with piecewise smooth compact Lipschitz boundaries. They proved, in particular, that if Ω is a curvilinear polygon whose smallest corner is θ, then for θ < π there holds CΩ = 2/(1 − cos θ), otherwise CΩ = 1. Pankrashkin [14] considered two- dimensional domains with a piecewise C4 smooth compact boundary and without convex Ω 2 2/3 corners, and it was shown that E1 (α) = −α − γα + O(α ), where γ is the maximum of the signed curvature at the boundary. Exner, Minakov, Parnovski [4] showed that for 4 Ω 2 2/3 compact C smooth domains the same asymptotics Em(α) = −α − γα + O(α ) holds for any fixed m ∈ N. Similar results were obtained by Exner and Minakov [3] for a class of two-dimensional domains with non-compact boundaries and by Pankrashkin and Popoff [15] for C3 compact domains in arbitrary dimensions. Cakoni, Chaulet, Haddar [1] studied the asymptotic behavior of higher eigenvalues. 1.2. Problem setting and the main result The computation of further terms in the eigenvalue asymptotics needs more pre- cise geometric assumptions. To our knowledge, such results are available for the two- dimensional case only. Helffer and Pankrashkin [9] studied the tunneling effect for the eigenvalues of a specific domain with two equal corners, and Helffer and Kachmar [8] considered the domains whose boundary curvature has a unique non-degenerate maximum. The machinery of both papers is based on the asymptotic properties of the eigenfunctions: it was shown that the eigenfunctions corresponding to the lowest eigenvalues concentrate near the smallest convex corner at the boundary or, if no convex corners are present, near the point of the maximum curvature, and this is used to obtain the corresponding eigenvalue asymptotics. The aim of the present note is to consider a new class of two-dimensional domains Ω. Namely, our assumption is as follows: 2 The domain R \ Ω is a convex polygon (with straight edges). Such domains are not covered by the above cited works: all the corners are non-convex, and the curvature is constant on the smooth part of the boundary, and it is not clear how the eigenfunctions are concentrated along the boundary. We hope that our result will be of use for the understanding of the role of non-convex corners. In order to formulate the main result, we first need some notation. We denote the 2 2 vertices of the polygon R \ Ω by A1,...,AM ∈ R , M ≥ 3, and assume that they are enumerated is such a way that the boundary ∂Ω is the union of the M line segments Ln := [An,An+1], n ∈ {1,...,M}, where we denote AM+1 := A1, A0 := AM . It is also assumed that there are no artificial vertices, i.e. that An ∈/ [An−1,An+1] for any n ∈ {1,...,M}. Furthermore, we denote by `n the length of the side Ln, and by Dn the Dirichlet 00 2 Laplacian f 7→ −f on (0, `n) viewed as a self-adjoint operator in L (0, `n). The main result of the present note is as follows: 48 K. Pankrashkin Theorem 1. For any fixed m ∈ N there holds: 1 EΩ (α) = −α2 + µD + O √ as α tends to +∞, m m α D where µm is the mth eigenvalue of the operator D1 ⊕ · · · ⊕ DM . The proof is based on the machinery proposed by Exner and Post [5] to study the convergence on graph-like manifolds. In reality, our construction appears to be quite similar to that of Post [16], which was used to study decoupled waveguides. We remark that due to the presence of non-convex corners the domain of the oper- Ω 2,2 ator Hα contains singular functions and is not included in W (Ω), see e.g. Grisvard [7]. This does not produce any difficulties, as our approach is purely variational and is entirely based on analysis of the sesqulinear form. 2. Preliminaries 2.1. Auxiliary operators 2 For α > 0, we denote by Tα the following self-adjoint operator in L (R+): 00 2,2 0 Tαv = −v , D(Tα) = v ∈ W (R+): v (0) + αv(0) = 0 . It is well known that: −αs 2 2 e spec Tα = {−α } ∪ [0, +∞), ker(T + α ) = Cϕα, ϕα(s) := √ . (2) 2α The sesqulinear form tα for the operator Tα looks as follows: ∞ Z 0 2 2 1,2 tα(v, v) = v (s) ds − α v(0) , D(tα) = W (R+). 0 1,2 Lemma 2. For any v ∈ W (R+), there holds: ∞ ∞ ∞ ∞ Z Z 2 Z Z 2 1 0 2 2 2 2 v(s) ds − ϕα(s)v(s)ds ≤ v (s) ds − α v(0) + α v(s) ds . α2 0 0 0 0 2 2 Proof. We denote by P the orthogonal projector on ker(Tα + α ) in L (R+), then by the spectral theorem, we have: 2 2 tα(v, v) + α kP vk = tα(v − P v, v − P v) ≥ 0, for any v ∈ D(tα). As ϕα is normalized, there holds ∞ Z ϕα(x)v(x)dx = kP vk, 0 and we arrive at the conclusion. Another important estimate is as follows, see Lemmas 2.6 and 2.8 in [12]: Lemma 3. Let Λ ⊂ R2 be an infinite sector of opening θ ∈ (0, 2π), then for any ε > 0 and any function v ∈ W 1,2(Λ) there holds: Z ZZ ZZ 2 2 2 Cθ 2 , θ ∈ (0, π), |v| ds ≤ ε |∇v| dx + |v| dx with Cθ = 1 − cos θ (3) ε ∂Λ Λ Λ 1, θ ∈ [π, 2π). Robin-Laplacian eigenvalues 49 2.2. Decomposition of Ω Let us proceed with a decomposition of the domain Ω which will be used through the 1 2 proof. Let n ∈ {1,...,M}. We denote by Sn and Sn the half-lines originating respectively at An and An+1, orthogonal to Ln and contained in Ω. By Πn, we denote the half-strip 1 2 bounded by the half-lines Sn and Sn and the line segment Ln, and by Λn we denote the 2 1 infinite sector bounded by the lines Sn−1 and Sn and contained in Ω. The constructions are illustrated in Fig. 1. We note that the 2M sets Λn and Πn, n ∈ {1,...,M}, are SM SM non-intersecting and that Ω = n=1 Λn ∪ n=1 Πn. We deduce from Lemma 3: Lemma 4. There exists a constant C > 0 such that for any ε > 0, any n ∈ {1,...,M} 1,2 and any v ∈ W (Λn) there holds Z ZZ 1 ZZ |v|2dσ ≤ Cε |∇v|2dx + |v|2dx . ε2 ∂Λn Λn Λn Furthermore, for each n ∈ {1,...,M} we denote by Θn the uniquely defined isome- try R2 → R2 such that: An = Θn(0, 0) and Πn = Θn (0, `n) × R+ : We remark that due to the spectral properties of the above operator Tα, see (2), we have, 1,2 for any u ∈ W (Πn), `n ∞ `n `n ∞ Z Z ∂ 2 Z 2 Z Z 2 u Θ (t, s) ds dt − α u Θ (t, s) dt + α2 u Θ (t, s) ds dt ∂s n n n 0 0 0 0 0 `n ∞ ∞ Z Z ∂ 2 2 Z 2 = u Θ (t, s) ds − α u Θ (t, 0) + α u Θ (t, s) ds dt ≥ 0, ∂s n n n 0 0 0 which implies, in particular, the following: `n ∞ Z Z ∂ 2 0 ≤ u Θ (t, s) ds dt ∂t n 0 0 `n ∞ `n ∞ Z Z ∂ 2 Z Z ∂ 2 ≤ u Θ (t, s) ds dt + u Θ (t, s) ds dt ∂t n ∂s n 0 0 0 0 `n `n ∞ Z 2 Z Z 2 2 − α u Θn(t, 0) dt + α u Θn(t, s) ds dt 0 0 0 ZZ Z ZZ = |∇u|2dx − α |u|2dσ + α2 |u|2dx. (4) Πn Ln Πn 2.3. Eigenvalues and identification maps We will use an eigenvalue estimate which is based on the max-min principle and is just a suitable reformulation of Lemma 2.1 in [5] or of Lemma 2.2 in [16]: 50 K. Pankrashkin FIG. 1. Decomposition of the domain Proposition 5. Let B and B0 be non-negative self-adjoint operators acting respectively in Hilbert space H and H0 and generated by sesqulinear form b and b0. Choose m ∈ N and assume that the operator B has at least m eigenvalues λ1 ≤ · · · ≤ λm < inf specess B and that the operator B0 has a compact resolvent. If there exists a linear map J : D(b) → D(b0) −1 (identification map) and two constants δ1, δ2 > 0 such that δ1 ≤ (1 + λm) , and that for any u ∈ D(b) there holds: 2 2 2 kuk − kJuk ≤ δ1 b(u, u) + kuk , 0 2 b (Ju, Ju) − b(u, u) ≤ δ2 b(u, u) + kuk , then 0 (λmδ1 + δ2)(1 + λm) λm ≤ λm + , 1 − (1 + λm)δ1 0 0 where λm is the mth eigenvalue of the operator B . 3. Proof of Theorem 1 3.1. Dirichlet-Neumann bracketing Consider the following sesqulinear form: M ZZ M ZZ Z Ω,D X 2 X 2 2 hα (u, u) = |∇u| dx + |∇u| dx − α |u| dσ , n=1 n=1 Λn Πn Ln M M Ω,D M 1,2 M 1,2 D(hα ) = W0 (Λn) ⊕ Wf0 (Πn), n=1 n=1 1,2 1,2 1 2 Wf0 (Πn) := f ∈ W (Πn): f = 0 at Sn ∪ Sn , Ω,D 2 Ω,D and denote by Hα the associated self-adjoint operator in L (Ω). Clearly, the form hα is Ω a restriction of the initial form hα , and due to the max-min principle we have: Ω Ω,D Em(α) ≤ Em (α), Robin-Laplacian eigenvalues 51 Ω,D Ω,D where Em (α) is the mth eigenvalue of Hα (as soon at it exists). On the other hand, we have the decomposition: M M Ω,D M D M D Hα = − ∆n ⊕ Gn,α, n=1 n=1 D 2 D where (−∆n ) is the Dirichlet Laplacian in L (Λn) and Gn,α is the self-adjoint operator in 2 L (Πn) generated by the sesquilinear form: ZZ Z D 2 2 D 1,2 gn,α(u, u) = |∇u| dx − α |u| dσ, D(gn,α) = Wf0 (Πn). Πn Ln Consider the following unitary maps: 2 2 Un : L (Πn) → L (0, `n) × R+ ,Unf := f ◦ Θn, n ∈ {1,...,M}, D ∗ D then it is straightforward to check that UnGn,αUn = Dn⊗1+Tα⊗1. As the operators (−∆n ) Ω,D Ω,D 2 D are non-negative, it follows that specess Hα = [0, +∞) and that Em (α) = −α + µm, which gives the majoration: Ω 2 D Em(α) ≤ −α + µm, (5) D 2 for all m with µm < α . In particular, the inequality (5) holds for any fixed m as α tends to +∞. Similarly, we introduce the following sesquilinear form: M ZZ M ZZ Z Ω,N X 2 X 2 2 hα (u, u) = |∇u| dx + |∇u| dx − α |u| dσ , n=1 n=1 Λn Πn Ln M M Ω,N M 1,2 M 1,2 D(hα ) = W (Λn) ⊕ W (Πn), n=1 n=1 Ω,N 2 and denote by Hα the associated self-adjoint operator in L (Ω). Clearly, the initial form Ω N,Ω hα is a restriction of the form hα , and due to the max-min principle we have: Ω,N Ω Em (α) ≤ Em(α), Ω,N Ω,N where Em (α) is the mth eigenvalue of Hα , and the inequality holds for those m for Ω which Em(α) exists. On the other hand, we have the decomposition: M M Ω,N M N M N Hα = − ∆n ⊕ Gn,α, n=1 n=1 N 2 N where (−∆n ) denotes the Neumann Laplacian in L (Λn) and Gn,α is the self-adjoint oper- 2 ator in L (Πn) generated by the sesquilinear form ZZ Z N 2 2 N 1,2 gn,α(u, u) = |∇u| dx − α |u| dσ, D(gn,α) = W (Πn). Πn Ln N ∗ 00 There holds UnGn,αUn = Nn ⊗ 1 + Tα ⊗ 1, where Nn is the operator f 7→ −f on (0, `n) with the Neumann boundary condition viewed as a self-adjoint operator in the Hilbert 2 N space L (0, `n), n ∈ {1,...,M}. The operators (−∆n ) are non-negative, and we have Ω,N Ω,N 2 N N specess Hα = [0, +∞) and Em (α) = −α + µm, where µm is the mth eigenvalue of the operator N1 ⊕ · · · ⊕ NM . Thus, we obtain the minorations: Ω 2 Ω 2 N Hα ≥ −α and Em(α) ≥ −α + µm, (6) 52 K. Pankrashkin which holds for any fixed m as α tends to +∞. By combining the inequalities (5) and (6) we also obtain the rough estimate: Ω 2 Em(α) = −α + O(1) for any fixed m and for α tending to +∞. (7) 3.2. Construction of an identification map In order to conclude the proof of Theorem 1, we will apply Proposition 5 to the operators: Ω 2 0 B = Hα + α ,B = D1 ⊕ · · · ⊕ Dn, which will allow us to obtain another inequality between the quantities: Ω 2 0 D λm = Em(α) + α , λm = µm. Note that for any fixed m ∈ N, one has λm = O(1) for large α, see (7). Therefore, it is sufficient to construct an identification map J = Jα as in Proposition 5 with δ1 + δ2 = O(α−1/2). Recall that the respective forms b and b 0 in our case are given by: Ω 2 2 Ω 1,2 b(u, u) = hα (u, u) + α kuk , D(b) = D(hα ) = W (Ω), ` M Z n 0 X 0 2 0 1,2 b (f, f) = fn(t) dt, D(b ) = f = (f1, . . . , fM ): fn ∈ W0 (0, `n) . n=1 0 Here and below, by kuk we mean the usual norm in L2(Ω). The positivity of b 0 is obvious, and the positivity of b follows from (6). Consider the maps: ∞ Z 1,2 2 Pn,α : W (Πn) → L (0, `n), (Pn,αu)(t) = ϕα(s)u Θn(t, s) ds, n ∈ {1,...,M}. 0 1,2 1,2 If u ∈ W (Ω), then u ∈ W (Πn) for any n ∈ {1,...,M}, and one can estimate, using the Cauchy-Schwarz inequality: ∞ ∞ Z Z 2 2 2 2 (Pn,αu)(0) + (Pn,αu)(`n) ≤ u Θn(0, s) ds + u Θn(`n, s) ds 0 0 Z Z = |u|2dσ + |u|2dσ. 1 2 Sn Sn 2 1 −1 As Sn−1 ∪ Sn = ∂Λn, we can use Lemma 4 with ε = α , which gives: M M Z Z X 2 2 X 2 2 (Pn,αu)(0) + (Pn,αu)(`n) ≤ |u| dσ + |u| dσ n=1 n=1 1 2 Sn Sn M M X Z C X ZZ ZZ = |u|2dσ ≤ |∇u|2dx + α2 |u|2dx . (8) α n=1 n=1 ∂Λn Λn Λn For each n ∈ {1,...,M}, we introduce a map: ` π : (0, ` ) → {0, ` }, π (t) = 0 for t < n , π (t) = ` otherwise, n n n n 2 n n Robin-Laplacian eigenvalues 53 ` and choose a function: ρ ∈ C∞ [0, ` ] with ρ (0) = ρ (` ) = 1 and ρ n = 0. n n n n n n 2 Finally, we define: M 1,2 M 2 Jα : W (Ω) → L (0, `n), (Jαu)n(t) = (Pn,αu)(t) − (Pn,αu) πn(t) ρn(t). n=1 1,2 1,2 We remark that (Jαu)n ∈ W0 (0, `n) for any u ∈ W (Ω) and n ∈ {1,...,M}, i.e. Jα maps D(b) into D(b 0) and will be used as an identification map. 3.3. Estimates for the identification map Take any δ > 0. Using the following inequality: 1 (a + a )2 ≥ (1 − δ)a2 − a2, a , a ≥ 0, 1 2 1 δ 2 1 2 we estimate: M M `n ZZ ZZ Z 2 2 2 X 2 X 2 kuk − kJαuk = |u| dx + |u| dx − Pn,αu (t) − Pn,αu π(t) ρ(t) dt n=1 n=1 Λn Πn 0 ` ` M ZZ M ZZ Z n Z n X 2 X 2 2 1 2 ≤ |u| dx + |u| dx − (1 − δ) Pn,αu (t) dt + Pn,αu π(t) ρ(t) dt δ n=1 n=1 Λn Πn 0 0 ` M ZZ M ZZ Z n X 2 X 2 2 = |u| dx + |u| dx − Pn,αu (t) dt n=1 n=1 Λn Πn 0 ` ` M Z n M Z n X 2 1 X 2 +δ Pn,αu (t) dt + Pn,αu π(t) ρn(t) dt δ n=1 0 n=1 0 =: I1 + I2 + I3 + I4. We have the trivial inequality: M 1 X ZZ ZZ I ≤ |∇u|2dx + α2 |u|2dx . 1 α2 n=1 Λn Λn To estimate the term I2, we use Lemma 2 and then (4): M `n ∞ ∞ Z Z Z 2 X 2 I2 = u Θn(t, s) ds − ϕα(s)u Θn(t, s) ds dt n=1 0 0 0 M `n ∞ ∞ Z Z 2 Z 1 X ∂ 2 2 2 ≤ u Θn(t, s) ds − α u Θn(t, 0) + α u Θn(t, s) ds dt α2 ∂s n=1 0 0 0 M 1 X ZZ Z ZZ ≤ |∇u|2dx − |u|2dσ + α2 |u|2dx , α2 n=1 Πn Ln Πn 54 K. Pankrashkin which gives: 1 I + I ≤ hΩ(u, u) + α2kuk2 . 1 2 α2 α Furthermore, with the help of the Cauchy-Schwarz inequality, we have: M `n ∞ M Z Z 2 ZZ X X 2 2 I3 ≤ δ u Θn(t, s) ds dt = δ |u| dx ≤ δkuk , n=1 n=1 0 0 Πn To estimate the last term, I4, we introduce the following constant: ` Z n n 2 o R := max ρn(t) dt : n ∈ {1,...,M} , 0 then, using first the estimate (8), and then the inequality (4), M R X 2 I ≤ sup P u π (t) 4 δ n,α n n=1 t∈(0,`n) M R X 2 2 ≤ (Pn,αu)(0) + (Pn,αu)(`n) δ n=1 M RC X ZZ ZZ ≤ |∇u|2dx + α2 |u|2dx δα n=1 Λn Λn M RC X ZZ ZZ ≤ |∇u|2dx + α2 |u|2dx δα n=1 Λn Λn M X ZZ Z ZZ + |∇u|2dx − |u|2dσ + α2 |u|2dx n=1 Πn Ln Πn RC = hΩ(u, u) + α2kuk . δα α Choosing δ = α−1/2 and summing up the four terms, we see that: c c kuk2 − kJ uk2 ≤ √1 hΩ(u, u) + α2kuk2 + kuk2 ≡ √1 b(u, u) + kuk2 , α α α α with a suitable constant c1 > 0. 0 Now, we need to compare b (Jαu, Jαu) and b(u, u). Take δ ∈ (0, 1) and use the inequality: 2 (a + a )2 ≤ (1 + δ)a2 + a2, a , a ≥ 0, 1 2 1 δ 2 1 2 Robin-Laplacian eigenvalues 55 Then: M `n Z 2 0 X 0 0 Ω 2 2 b (Jαu, Jαu) − b(u, u) = (Pn,αu) − ρ n (Pn,αu) ◦ πn dt − hα (u, u) + α kuk n=1 0 M `n M `n X Z 2 2 X Z 2 ≤(1 + δ) (P u)0 dt + ρ 0 (P u) ◦ π dt n,α δ n n,α n n=1 0 n=1 0 M X ZZ ZZ − |∇u|2dx + α2 |u|2dx n=1 Λn Λn M ZZ Z ZZ X 2 2 2 2 − |∇u| dx − |u| dσ + α |u| dx (9) n=1 Πn Ln Πn M `n M `n X Z 2 2 X Z 2 ≤ (1 + δ) (P u)0 dt + ρ 0 (P u) ◦ π dt n,α δ n n,α n n=1 0 n=1 0 M X ZZ Z ZZ − |∇u|2dx − |u|2dσ + α2 |u|2dx . n=1 Πn Ln Πn Using first the Cauchy-Schwarz inequality and then inequality (4), we have: `n `n ∞ Z Z Z 2 ZZ Z ZZ 0 2 ∂ 2 2 2 2 (Pn,αu) dt ≤ u Θn(t, s) ds dt ≤ |∇u| dx − |u| dσ + α |u| dx. ∂t 0 0 0 Πn Ln Πn Substituting the last inequality into (9), we arrive at: M ZZ Z ZZ 0 X 2 2 2 2 b (Jαu, Jαu) − b(u, u) ≤ δ |∇u| dx − |u| dσ + α |u| dx n=1 Πn Ln Πn (10) M `n 2 X Z 2 + ρ0 (P u) ◦ π dt. δ n n,α n n=1 0 Furthermore, using the constant: ` Z n 0 n 0 2 o R := max ρn(t) dt : n ∈ {1,...,M} , 0 and the inequality (8), we have: M `n M Z 2 X 0 0 X 2 ρn (Pn,αu) ◦ πn dt ≤ R sup (Pn,αu) πn(t) t∈(0,`n) n=1 0 n=1 M 0 M ZZ ZZ 0 X 2 2 R C X 2 2 2 ≤ R (Pn,αu)(0) + (Pn,αu)(`n) ≤ |∇u| dx + α |u| dx . α n=1 n=1 Λn Λn 56 K. Pankrashkin The substitution of this inequality into (10) and the choice δ = α−1/2 then leads to: c c b 0(Ju, Ju) − b(u, u) ≤ √2 hΩ(u, u) + α2kuk2 ≤ √2 b(u, u) + kuk2 , α α α with a suitable constant c2 > 0. By Proposition 5, for any fixed m ∈ N and for large α, we D Ω 2 −1/2 have the estimate µm ≤ Em(α) + α + O(α ). The combination with (5) gives the result. Acknowledgments The work was partially supported by ANR NOSEVOL (ANR 2011 BS01019 01) and GDR Dynamique quantique (GDR CNRS 2279 DYNQUA). References [1] F. Cakoni, N. Chaulet, H. Haddar. On the asymptotics of a Robin eigenvalue problem. Comptes Rendus Math., 351 (13–14), P. 517–521 (2013). [2] D. Daners, J.B. Kennedy. On the asymptotic behaviour of the eigenvalues of a Robin problem. Differen- tial Integr. Equ., 23 (7–8), P. 659–669 (2010). [3] P. Exner, A. Minakov. Curvature-induced bound states in Robin waveguides and their asymptotical properties. J. Math. Phys., 55 (12), P. 122101 (2014). [4] P. Exner, A. Minakov, L. Parnovski. Asymptotic eigenvalue estimates for a Robin problem with a large parameter. Portugal. Math., 71 (2), P. 141–156 (2014). [5] P. Exner, O. Post. Convergence of spectra of graph-like thin manifolds. J. Geom. Phys., 54 (1), P. 77– 115 (2005). [6] T. Giorgi, R. Smits. Eigenvalue estimates and critical temperature in zero fields for enhanced surface superconductivity. Z. Angew. Math. Phys., 58 (2), P. 224–245 (2007). [7] P. Grisvard. Singularities in boundary value problems. Recherches en Mathematiques´ Appliques, No. 22 (1992). [8] B. Helffer, A. Kachmar. Eigenvalues for the Robin Laplacian in domains with variable curvature. arXiv:1411.2700 (2014). [9] B. Helffer, K. Pankrashkin. Tunneling between corners for Robin Laplacians. J. London Math. Soc., 91 (1), P. 225248 (2015). [10] P. Freitas, D. Krejciˇ rˇ´ık. The first Robin eigenvalue with negative boundary parameter. arXiv:1403.6666, (2014). [11] A.A. Lacey, J.R. Ockendon, J. Sabina. Multidimensional reaction diffusion equations with nonlinear boundary conditions. SIAM J. Appl. Math., 58 (5), P. 1622–1647 (1998). [12] M. Levitin, L. Parnovski. On the principal eigenvalue of a Robin problem with a large parameter. Math. Nachr., 281 (2), P. 272–281 (2008). [13] Y. Lou, M. Zhu. A singularly perturbed linear eigenvalue problem in C1 domains. Pacific J. Math., 214 (2), P. 323–334 (2004). [14] K. Pankrashkin. On the asymptotics of the principal eigenvalue for a Robin problem with a large parameter in planar domains. Nanosystems: Phys. Chem. Math., 4 (4), P. 474–483 (2013). [15] K. Pankrashkin, N. Popoff. Mean curvature bounds and eigenvalues of Robin Laplacians. arXiv:1407.3087, (2014). [16] O. Post. Branched quantum wave guides with Dirichlet boundary conditions: the decoupling case. J. Phys. A, 38 (22), P. 4917–4932 (2005). NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS, 2015, 6 (1), P. 57–62 THERMALLY INDUCED TRANSITIONS AND MINIMUM ENERGY PATHS FOR MAGNETIC SYSTEMS E. V. Alexandrov ITMO University, St. Petersburg, Russia [email protected] PACS 75.10.Hk, 82.20.Kh, 02.60.-x DOI 10.17586/2220-8054-2015-6-1-57-62 Thermally induced magnetic transitions are rare events as compared with vibrations of individual magnetic mo- ments. Timescales for these processes differ by 10 orders of magnitude or more. Therefore, the standard Monte- Carlo simulation is not suitable for the theoretical description of such phenomena. However, a statistical approach based on transition state theory is applicable for calculations of the transition rates. It presupposes finding the minimum energy path (MEP) between stable magnetic states on the multidimensional energy surface of the sys- tem. A modification of the Nudged Elastic Band (NEB) method for finding the energy barriers between states is suggested. A barrier on the energy surface corresponds to the difference between maximum energy along the MEP (highest saddle point) and the initial state minimum. The NEB procedure is implemented for spin rotations in Cartesian representation with geometric constraint on the magnitude of the magnetic moment. In this case, the effective magnetic forces are restricted to the tangent plane of the magnetic momentum vector. Keywords: Potential energy surface, minimum energy path, nudged elastic band method, numerical optimization, quick min-mode. Received: 31 December 2014 1. Introduction In recent years, there has been large interest in the study of magnetization reversal processes under the action of thermal fluctuations and external perturbations. The stability of magnetic states is of great importance for research and development of miniaturized magnetic storage devices [1]. For estimation of the rates of activated transitions it is essential to know the response of the system to random perturbations. Direct Monte-Carlo simulation of temper- ature induced magnetization reversal is effectively impossible due to the time-scale difference between high-frequency vibrations of individual moments and low-frequency magnetic transi- tions between magnetic states of the whole system. However, a statistical approach known as ‘transition state theory’ [2, 3] can be used in this case. This approach describes transitions through description of behavior of a system in the least observable ‘transition state’. Stable states of the system correspond to the minima on the multidimensional energy surface of the system considered as a function of all parameters determining the magnetic configuration. For magnetic systems, these parameters have to determine all spin vectors. These could be Cartesian coordinates of the magnetic moments. However, due to the large difference in the time scale of longitudinal and transverse relaxation of the magnetic moment, the length can be taken to be fixed and thereby decrease the number of independent variables. Transition rates depend on the height of energy barriers, which correspond to saddle points on minimum energy paths (MEP), as well as on the shape of the energy surface at the minimum and saddle points [2]. A popular and efficient approach used for calculation of MEP for atomic rearrangements is the nudged elastic band (NEB) method [3]. A transition path is represented there by a chain of images (replicas) of the system between two local minima on the potential energy surface; each image 58 E. V. Alexandrov corresponds to specific configuration of the system with known potential energy. The images are connected with artificial springs to keep them well separated and prevent them from falling into the minima. The images give a discrete representation of the path on the energy surface. An iterative algorithm moves the images under the action of a force that contains the transverse component of the real force and parallel component of the spring force. Detailed description of this procedure, taking into account the constraint on the modulus of magnetic moments, is given below. 2. Method 2.1. Minimization technique The system is defined by the Hamiltonian which gives the energy of the system for any configuration of the spin vectors. It can be either an analytical function of spin vectors, or numerical procedure, such as DFT calculations [4]. A constraint on the length of the spin vec- tors confining them to the surface of a sphere or a self-consistent determination of the modulus of magnetic moments require specific measures to be incorporated into the implementation of the NEB method. Common constrained-optimization techniques, such as Lagrange multipliers method, are not suitable in this case due to the computational cost which tends to be high for non-analytical Hamiltonians [3,5]. In addition, in most cases pure analytical methods for finding extrema are unsuitable due to the computational difficulty for evaluating the second derivatives of the Hamiltonian. Thereby, the implementation method itself should preferably incorporate these constraints. One way to do this is to use a coordinate system that is bounded to sur- face of the sphere, such as spherical or stereographical coordinates [5]. This approach is quite successful; however, it suffers from topological artifacts of a sphere, such as coordinate peri- odicity and problems with pole regions. As an alternative to using sphere-specific coordinates and unconstrained transformations, we propose to use Cartesian coordinates and sphere-specific constrained transformations, e.g. rotations. Using this approach in the optimization procedure allows one to overcome most of the topological artifacts of a sphere and prevent violation of native constraints. We use a minimization technique known as ‘quick-min’ or ‘velocity projection opti- mization’ [4, 5], which is essentially a gradient descent method, optimized with accumulation of velocity in direction of minimization movement. An iteration step can be formulated in four stages as follows: 1 ~a = − ∇E S~ , i m i ~a + ~a ~v = ~v + i i−1 ∆t, i i−1 2 h~ai|~vii χ (h~ai0|~vi0i) ~vi = k~vik h~ai0|~vi0i χ (h~ai0|~vi0i)~ai0 = 2 ~ai, k~aik ~a S~ = S~ + ~v ∆t + i ∆t2, i+1 i i 2 where m and ∆t are optimization parameters and χ is Heaviside function. The method is inspired by the velocity Verlet integration of motion equation. The native constraint requires the force to be restricted to the tangent plane of the current spin vector. If so, movement on tangent plane can be calculated. An infinitesimal movement in the tangent plane is equal to movement on the surface of the sphere; thereby we can transform that shift into rotation in Thermally induced transitions and minimum energy paths for magnetic systems 59 movement plane. We can construct rotation matrix in axis-angle representation from ~n – normal to movement plane and θ – rotation angle: R (~n,θ) = 2 2 nx + cos θ (1 − nx) nxny (1 − cos θ) − nz sin θ nxnz (1 − cos θ) + ny sin θ 2 2 . (1) nxny (1 − cos θ) + nz sin θ ny + cos θ 1 − ny nynz (1 − cos θ) − nx sin θ 2 2 nxnz (1 − cos θ) − ny sin θ nynz (1 − cos θ) + nx sin θ nz + cos θ (1 − nz) The minimization step is performed by using matrix (1) in following form: ~ ~ Si+1 = R (~ni, θi) · Si. In principle, instead of rotation matrices, another type of object can be used – unit quaternions, also known as versors. From ~n, θ, quaternion describing the same rotation is defined as follows: θ θ q (s,~v) = cos , ~n sin . (2) 2 2 If quaternion is defined as in (2), then minimization step can be written in following form: ~ ~ −1 Si+1 = qi · Si · qi , ~ where Si is converted to quaternion by adding unit scalar part to it. Accurately implemented, both rotation matrix approach and quaternion approach have similar computational cost. Quaternions have major computational benefit: the composition of the rotations can be calculated by simple multiplication of corresponding quaternions, which is much faster than multiplication of rotation matrices, but the difference was not significant in the present applications. However, it seems worth looking into yet another spin representation for this task: each spin vector as unit vector of a surface of a sphere can be constructed by rotation from a selected start vector. Therefore, it is possible to use rotations as coordinates. It may be of large interest to implement the procedure described above solely using quaternion terms. 2.2. Calculation of minimum energy path The NEB method is implemented as minimization of all intermediate images of the path, using the specially modified force. In addition to the true force, each image is affected by forces of artificial springs that connect it to its neighbors. This force for each spring is determined from length of arc segment between two spin vectors in different images. To prevent the ‘corner-cutting’ effect [3], nudging direction is introduced using ‘upwind tangent” approach [6]. In this approach, nudging direction is chosen as direction to highest energy neighbor. After that, forces are projected dependent on that direction: ~ ~ ~ FΣ = FT ⊥ + FE||. With this modified force, the iterative minimization procedure moves the images to the MEP. The saddle point, which corresponds to maximum energy on MEP, can be estimated from two neighboring images with highest energy. In addition, ‘climbing-image’ (CI) modification [7] can be used: The highest energy image marked as ‘climbing’, is free from spring forces and parallel component of the true force is inverted to move that image higher along MEP to the saddle point. 60 E. V. Alexandrov 3. Test cases Two simple systems were chosen as test cases. Both represent a single atom with two degrees of freedom to make it easier to visualize. In that case potential energy surface can be shown as a two-dimensional surface in three-dimensional space. In these particular cases, spherical coordinates were used for plotting: φ – azimuth angle and θ – elevation angle. The potential energy surfaces, plotted as function of these angles, are shown in Fig. 1 and Fig. 2. In both cases there are two minima corresponding to φ = θ = 0 and φ = π, θ = 0. The initial path corresponds to the chain shown with empty circles and is chosen to be the same in both cases. Images of a system in that state were intentionally distributed non-equidistantly through predetermined randomized procedure to show that NEB procedure is not constricted to a limited set of initial path configurations. The final MEP, obtained using the algorithm described above, corresponds to the chain shown with solid circles. The saddle point is denoted by a cross mark. FIG. 1. Potential energy surface for single spin with spin Hamiltonian (3). Calcu- lation of MEP without Climbing Image. Saddle point is directly between minima. Maximums are in pole regions (θ = ±π/2) The first system is a single spin with Hamiltonian containing only internal anisotropy part: 4 0 0 ~T ~ H = S 0 1 0 S. (3) 0 0 −9 For this system, the MEP is expected to be a straight line between the local minima. In Fig. 1, this case is displayed, with even number of intermediary images and without CI. The saddle point is estimated to be between the two images with highest energy on opposite slopes. Due to symmetry, the point in the center is the saddle point. Also, if an odd number of intermediate images were used, then the middle image would be expected to get to saddle point even without using CI. Note that in Fig. 1, the distances between adjacent images in the final configuration of the path are equal. The second system is single spin with Gaussian Hamiltonian: Thermally induced transitions and minimum energy paths for magnetic systems 61 FIG. 2. Potential energy surface for single spin with Gaussian Hamiltonian (4). Calculation of MEP with Climbing Image. Saddle point is exactly in climbed image. There are maximums in pole regions (θ = ±π/2), and there is maximum directly between minima (ϕ = π/2, θ = 0), which results in curved form of MEP X αi H = w exp − , i 0.35 i n o w = −5, −5, 1, 3.5, 4, 4 , (4) −1 1 0 0 0 0 P = 0 , 0 , −1 , 1 , 0 , 0 , 0 0 0 0 −1 1 ~ where αi is the angle between Pi and spin vector. For this system, the MEP is non-trivial; it is a curved line since now there is a maximum in the energy at the center. In Fig. 2, the results of a calculation which was carried out with CI turned on. The saddle point is estimated to be exactly at the climbing image. The usage of CI results in splitting the band in two parts, separated by saddle point, with a different distance between adjacent images on the two sides. Due to symmetry, if either number of images were odd or CI were turned off, then the images would have even distribution of distances between them. Note that in Fig. 1, there are two possible MEPs, and in Fig. 2, there are four. This illustrates that the NEB method finds the local MEP closest to the initial path, not necessarily the global MEP. The final configuration depends in that way on the initial path. These cases show that the NEB implementation described above is in accordance with theoretical description [3] and behaves in the expected way. There are no problems specific to pole regions or periodicity of coordinates, which is good for application to magnetic systems. In the future, this implementation should be compared with other implementations in terms of speed of convergence. Also, it may be interesting to try a pure quaternion implementation. 62 E. V. Alexandrov Acknowledgements This work was partially supported by RFBR grant No. 14-02-00102 and Nordic-Russian Training Network for Magnetic Nanotechnology (NCM-RU10121). Helpful discussions with Pavel Bessarab, Valery Uzdin and Hannes Jonsson´ are gratelfully acknowledged. References [1] Loth S., Baumann S., et al. Bistability in atomic-scale antiferromagnets. Science, 335, P. 196–198 (2012). [2] Bessarab P.F., Uzdin V.M., Jonsson´ H. Harmonic transition-state theory of thermal spin transitions. Phys. Rev. B, 85, P. 184409 (2012). [3] Jonsson´ H., Mills G., Jacobsen K.W. Nudged elastic band method for finding minimum energy paths of transitions. Proceedings of the Conference ‘International School of Physics. Classical and quantum dynamics in condensed phase simulations’, Lerici, Italy, July 1997, P. 385–404 (1997). [4] Henkelman G., Johannesson´ G., Jonsson´ H. Methods for finding saddle points and minimum energy paths, in Progress on Theoretical Chemistry and Physics. Kluwer Academic Publishers, Dordrecht, the Netherlands, 2000, P. 269–300 (2000). [5] Bessarab P.F., Uzdin V.M., Jonsson´ H. Potential energy surfaces and rates of spin transitions. Zeitschrift fur Physikalische Chemie, 227, P. 1543 (2013). [6] Henkelman G., Jonsson´ H. Improved tangent estimate in the NEB method for finding minimum energy paths and saddle points. J. Chem. Phys., 113, P. 9978 (2000). [7] Henkelman G., Uberuaga B., Jonsson´ H. A Climbing-Image NEB Method for finding saddle points and minimum energy paths. J. Chem. Phys., 113, P. 9901 (2000). NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS, 2015, 6 (1), P. 63–78 ON SOME APPLICATIONS OF THE BOUNDARY CONTROL METHOD TO SPECTRAL ESTIMATION AND INVERSE PROBLEMS S. A. Avdonin1, A. S. Mikhaylov2,3, V.S. Mikhaylov2,4 1Department of Mathematics and Statistics, University of Alaska, Fairbanks, USA 2St. Petersburg Department of V.A. Steklov Institute of Mathematics of the Russian Academy of Sciences, St. Petersburg, Russia 3St. Petersburg State University, Faculty of Mathematics and Mechanics, St. Petersburg, Russia 4St. Petersburg State University, Faculty of Physics, St. Petersburg, Russia [email protected], [email protected], [email protected] PACS 02.30 Zz, 02.30 Yy, 02.30 Nw DOI 10.17586/2220-8054-2015-6-1-63-78 We consider applications of the Boundary Control (BC) method to generalized spectral estimation problems and to inverse source problems. We derive the equations of the BC method for these problems and show that the solvability of these equations crucially depends on the controllability properties of the corresponding dynamical system and properties of the corresponding families of exponentials. Keywords: spectral estimation problem, boundary control method, identification problem, inverse problem, Schrodinger¨ equation, hyperbolic system. Received: 5 December 2014 1. Introduction The classical spectral estimation problem consists of recovering the coefficients an, λk, k = 1,...,N,N ∈ N, of a signal N X λkt s(t) = ake , t > 0 n=1 from the given observations s(j), j = 0,..., 2N − 1, where the coefficients ak, λk may be arbitrary complex numbers. The literature describing various methods for solving the spectral estimation problem is very extensive: see for example the list of references in [1, 2]. In these papers a new approach to this problem was proposed: a signal s(t) was treated as a kernel of a certain convolution operator corresponding to an input-output map for some linear discrete- time dynamical system. While the system realized from the input-output map is not unique, the coefficients an and λn can be determined uniquely using the non-selfadjoint version of the boundary control method [3]. In [4, 8], this approach has been generalized to the infinite-dimensional case: more precisely, the problem of the recovering the coefficients ak, λk ∈ C, k ∈ N, of the given signal: ∞ X λkt S(t) = ak(t)e , t ∈ (0, 2T ), (1.1) k=1 from the given data S ∈ L2(0, 2T ) was considered. In [4], the case ak ∈ C has been treated, PLk−1 i i in [8] the case when for each k, ak(t) = i=0 akt are polynomials of the order Lk − 1 with i complex valued coefficients ak was studied. 64 S. A. Avdonin, A. S. Mikhaylov, V.S. Mikhaylov Recently, it was observed [9, 15] that the results of [4, 8] are closely related to the dynamical inverse source problem: let H be a Hilbert space, A be an operator in H with the domain D(A), Y be another Hilbert space, O : H ⊃ D(O) 7→ Y be an observation operator (see [18]). Given the dynamical system in H: ( u − Au = 0, t > 0, t (1.2) u(0) = a, we denote by ua its solution, and by y(t) := (Oua)(t) the observation (output of this system). The operator that realizes the correspondence a 7→ (Oua)(t) is called the observation operator T O : H 7→ L2(0,T ; Y ). We fix some T > 0 and assume that y(t) ∈ L2(0,T ; Y ). One can pose the following questions: what information on the operator A could be recovered from the observation y(t)? We mention works on the multidimensional inverse problems for the Schrodinger,¨ heat and wave equations by one measurement, concerning this subject. Some of the results (for the Schrodinger¨ equation) are given in [9, 10, 16]. To answer this question in the abstract setting, in [15] the authors derived the version of the BC-method equations under the condition that A is self-adjoint and Y = R. In the present paper, we address the same question without the assumption about selfajointness of A . The possibility of recovering the spectral data from the dynamical one is well-known for the dynamical system with a boundary control [11, 12]. We extend these ideas to the case of the dual (observation) system. The solvability of the BC-method equations for the spectral estimation problem critically depends on the properties of corresponding exponential family. The solvability of the BC- method equations for system (1.2) depends on the controllability properties of the dual system. We point out the close relation between these two problems: they both leads to essentially the same equations (see section 4 for applications), and conditions for the solvability of these equations are the same (on the connections between the controllability of a dynamical systems and properties of exponential families see [5]). In the second section, we outline the solution for the spectral estimation problem in infinite dimensional spaces (see [8] for details). In the third section, we derive the equations of the BC-method for problem (1.2), extending the results of [15] to the case of non self-adjoint operator. Also, we answer the question on the extension of the observation y(t) = (Oa)(t). The last section is devoted to the applications to inverse problem by one measurement of the Schrodinger¨ equation on the interval and to the problem of extending the inverse data for the first order hyperbolic system on the interval, see also [4, 7–9]. 2. The spectral estimation problem in infinite dimensional spaces The problem is set up in the following way: given the signal (1.1), S ∈ L2(0, 2T ), for T > 0, to recover the coefficients ak(t), λk, k ∈ N. Below, we outline the procedure of recovering unknown parameters, for the details see [8]. We consider the dynamical systems in a complex Hilbert space H: x˙(t) = Ax(t) + bf(t), t ∈ (0,T ), x(0) = 0. (2.1) y˙(t) = A∗y(t) + dg(t), t ∈ (0,T ), y(0) = 0, (2.2) ∞ here b, d ∈ H, f, g ∈ L2(0,T ), and we assume that the spectrum of the operator A, {λk}k=1 is not simple. We denote the algebraic multiplicity of λk by Lk, k ∈ N, and also assume that the i set of all root vectors {φk}, i = 1,...,Lk, k ∈ N, forms a Riesz basis in H. Here, the vectors i Lk from the chain {φk}i=1, k ∈ N, satisfy the equations: 1 i i−1 (A − λk) φk = 0, (A − λk) φk = φk , 2 6 i 6 Lk. On some applications of the Boundary Control method to spectral estimation and . . . 65 ∗ ∞ i The spectrum of A is {λk}k=1 and the root vectors {ψk}, i = 1,...,Lk, k ∈ N, also form a Riesz basis in H and satisfy the equations: