Transparent Quantum Graphs
J.R. Yusupov, K.K. Sabirov, M. Ehrhardt, D.U. Matrasulov

TL;DR
This paper develops a method for designing quantum graphs with transparent branching points using boundary conditions that simplify the analysis of quantum networks, applicable to star and tree topologies.
Contribution
It introduces a novel approach to derive transparent boundary conditions for quantum graphs, simplifying their analysis and design.
Findings
Derived simple constraints for transparent boundary conditions.
Applied the method to star and tree quantum graphs.
Extension to complex topologies is straightforward.
Abstract
We consider quantum graphs with transparent branching points. To design such networks, the concept of transparent boundary conditions is applied to the derivation of the vertex boundary conditions for the linear Schrodinger equation on metric graphs. This allows to derive simple constraints, which use equivalent usual Kirchhoff-type boundary conditions at the vertex to the transparent ones. The approach is applied to quantum star and tree graphs. However, extension to more complicated graph topologies is rather straight forward.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Transparent Quantum Graphs
J.R. Yusupov1, K.K. Sabirov2, M. Ehrhardt3 and D.U. Matrasulov1
1Turin Polytechnic University in Tashkent, 17 Niyazov Str., 100095, Tashkent, Uzbekistan
2Tashkent University of Information Technologies, 108 Amir Temur Str., 100200, Tashkent Uzbekistan
3Bergische Universität Wuppertal, Gaußstrasse 20, D-42119 Wuppertal, Germany
Abstract
We consider quantum graphs with transparent branching points. To design such networks, the concept of transparent boundary conditions is applied to the derivation of the vertex boundary conditions for the linear Schrödinger equation on metric graphs. This allows to derive simple constraints, which use equivalent usual Kirchhoff-type boundary conditions at the vertex to the transparent ones. The approach is applied to quantum star and tree graphs. However, extension to more complicated graph topologies is rather straight forward.
I Introduction
The problem of wave transport in branched structures and networks is of importance for many areas of contemporary physics, such as optics, fluid dynamics, condensed matter and polymers. Optical and quantum mechanical waves propagation in such systems appears e.g., in different branched waveguides BWG1 ; BWG2 ; BWG3 ; BWG4 ; BWG5 ; BWG6 ; BWG7 . Effective transfer of light, charge, heat, spin and signal in networks requires solving the problem of tunable wave dynamics. This can be achieved using the models, which provide realistic description of particle and wave transport in branched structures.
An important feature of particle and wave dynamics in networks is the transmission through the branching points, which is usually accompanied by the reflection (backscattering) of a wave at these points. Dominating of reflection compared to transmission implies large “resistivity” of a network with respect to the wave propagation. Therefore, it is important from the viewpoint of practical applications, to reduce such resistivity by providing a minimum of reflection, or by its absence. This task leads to the problem of tunable transport in branched structures, whose ideal result should be reflectionless transmission of the waves through the branching points of the structure. From the viewpoints of practical applications in condensed matter, such transmission implies ballistic transport of charge, spin, heat and other carriers in low-dimensional branched materials. The latter is of importance for the effective utilization of different functional materials.
Reflectionless transport of waves in optical fiber networks is one of the central tasks for the fiber optics, as many information-communication devices (e.g., computers, computer networks, telephones, etc.) use optical fiber networks for information (signal) transfer. Such networks are also used in different optoelectronic devices. High speed and lossless transfer of information in such devices require minimum of backscattering or its absence. Important areas, where the reflectionless, or ballistic transport in branched structures is required, are molecular electronics and conducting polymers. Both deal with the modeling of reflectionless wave transport in networks and require revealing the conditions providing such regime. Different quantum wire networks appearing in solid state physics are also potential structures, where the wave transport can be tuned from diffusive to ballistic regime using the transparent boundary conditions.
In this paper we study the problem of reflectionless transport in quantum networks, which are modeled in terms of the linear Schrödinger equations on metric graphs. To describe the regime, when there is no backscattering at the network branching points, we use the concept of transparent boundary conditions for the time-dependent Schrödinger equation studied earlier in detail in the Refs. Arnold1998 ; Ehrhardt1999 ; Ehrhardt2001 ; Ehrhardt2002 ; Arnold2003 ; Jiang2004 ; Antoine2008 ; Ehrhardt2008 ; Sumichrast2009 ; Antoine2009 ; Ehrhardt2010 ; Klein2011 ; Arnold2012 ; Feshchenko2013 ; Antoine2014 ; Petrov2016 . Combining transparent boundary conditions concept with that of vertex boundary conditions for the linear Schrödinger equation on metric graphs, we derive the constraints, providing the regime of reflectionless transmission of the waves through the graph branching points (vertices).
More precisely, we reveal the condition for the regime, at which Kirchhoff-type boundary conditions at the branching point become equivalent to the transparent boundary conditions. The condition is obtained in a simple form of the sum rule for the coefficients, which can characterize physical properties of the network branches. We note that the linear and nonlinear wave equations on metric graphs have been studied earlier in the context of quantum graphs Uzy1 ; Hul ; Kuchment04 ; Uzy2 ; Exner15 and soliton dynamics in networks Zarif ; zar2011 ; caputo14 ; Our2015 ; Adami16 ; Our1 ; Adami17 ; dimarecent ; Karim2018 ; KarimNLDE ; Our2018 . It was found in the Refs. Zarif ; Our1 ; KarimNLDE that transmission of the waves through the network branching point can be reflectionless, when provided certain constraints are fulfilled. However, no strict explanation for such effect have been presented in those studies.
This paper can be considered as the first step in the way for the formulation of clear, physically reasonable, strict mathematical conditions for the reflectionless wave transmission in network branching points in quantum regime. The paper is organized as follows. In the next section we briefly recall the concept of transparent boundary conditions for the linear Schrödinger equation on a line. In Section III we recall the concept of quantum graphs. Section IV provides extension of the concept of transparent boundary conditions to quantum networks described by the linear Schrödinger equation on metric graphs and presents some numerical results. Finally, Section V presents some concluding remarks.
II Transparent boundary conditions in quantum mechanics
The problem of transparent boundary conditions (TBC) for the wave equations has attracted much attention in different practically important contexts (see, e.g., papers Arnold1998 ; Ehrhardt1999 ; Ehrhardt2001 ; Ehrhardt2002 ; Arnold2003 ; Jiang2004 ; Antoine2008 ; Ehrhardt2008 ; Sumichrast2009 ; Antoine2009 ; Ehrhardt2010 ; Klein2011 ; Arnold2012 ; Feshchenko2013 ; Antoine2014 ; Petrov2016 for review). Such boundary conditions describe the absorption or reflectionless transmission of particles (waves) at the boundary of two domains. Approximate TBCs for the linear Schrödinger equation were formulated by Shibata Shibata and Kuska Kuska , where dispersion relations for the plane waves is approximated to derive TBCs. A strict mathematical analysis of TBC for different wave equations, including quantum mechanical Schrödinger equation can be found in the Refs. Ehrhardt2002 ; Arnold2003 ; Ehrhardt2008 ; Antoine2014 ; Petrov2016 .
The explicit form of such boundary conditions are much more complicated than those of Dirichlet, Neumann and Robin conditions. Therefore, their discretization causes a serious accuracy loss in the numerical computations and also often reduce the stability range of the overall scheme. Here we briefly recall the formulation of TBCs on a line following the Refs. Ehrhardt1999 ; Ehrhardt2001 ; Ehrhardt2002 ; Arnold2003 .
Let us consider the wave (particle) motion in a 1D domain and described by the following time-dependent Schrödinger equation (with ):
[TABLE]
with the initial condition given as
[TABLE]
where , and is an external potential. Our purpose is to formulate the boundary conditions for Eq. (1), which provide reflectionless transmission of the wave via the given points, [math] and . One of the prescriptions to design such boundary conditions was developed in the Refs. Ehrhardt1999 ; Ehrhardt2001 ; Ehrhardt2002 and uses the following two basic assumptions: the initial data is compactly supported in the computational domain , and the given external potential is constant outside this finite domain, i.e. for for .
Then we can separate the whole problem into the so-called “interior” and “exterior” problems. The “interior” problem is given by
[TABLE]
where denotes the Dirichlet-to-Neumann (DtN) maps at the boundaries, which can be obtained by solving the two “exterior” problems given as
[TABLE]
and it can analogously be done for .
Using the Laplace transformation
[TABLE]
the right “exterior” problem is transformed to
[TABLE]
The general solution for Eq. (5) can be written as
[TABLE]
Since its solutions have to decrease as , we obtain
[TABLE]
where denotes the square root with nonnegative real part.
Hence the Laplace-transformed Dirichlet-to-Neumann operator reads
[TABLE]
and is calculated analogously.
An inverse Laplace transformation yields the right TBC at :
[TABLE]
Similarly, the left TBC at is obtained as
[TABLE]
In this paper we apply this procedure to quantum graphs. It should be noted that no practical applications of the TBCs in physical systems have been considered so far. Here we will do this for quantum networks, which appear in different branches of optics, condensed matter and polymer physics.
III Quantum graphs
The concept of quantum graphs has been introduced first by Exner, Seba and Stovicek to describe free quantum motion on branched wires Exner1 . However, the pioneering treatment of the quantum mechanical motion in branched structures dates back to the Refs. Pauling ; Rud ; Alex .
Later Kostrykin and Schrader derived the general boundary conditions providing self-adjointness of the Schrödinger operator on graphs Kost . Relativistic quantum mechanics described by Dirac Bolte and Bogoliubov-de Gennes operators KarimBdG on graphs have been studied recently. Hul et al considered experimental realization of quantum graphs in optical microwave networks Hul . An important topic related to quantum graphs was studied in the context of quantum chaos theory and spectral statistics Uzy1 ; Bolte ; Uzy2 . Spectral properties and band structure of periodic quantum graphs also attracted much interest Exner15 ; Grisha . Different aspects of the Schrödinger operators on graphs have been studied in the Refs. Exner15 ; Grisha ; Mugnolo . Tunable directed transport in quantum driven graphs has been also studied in yusupov2015 , where the reflection at the branching point has been discussed.
Nonlinear extensions of the quantum graph concept have been studied in the context of soliton transport in networks in the Refs. Zarif ; zar2011 ; caputo14 ; Our2015 ; Our1 ; KarimNLDE ; Our2018 . In Zarif the nonlinear Schrödinger equation on metric graphs is studied and condition for integrability is derived in the form of a sum rule for nonlinearity coefficients. In zar2011 such study is extended to Ablowitz-Laddik equation. The stationary Schrödinger equation on metric graphs and standing wave solitons in networks are studied in Karim2013 ; Adami16 ; Adami17 . In Ourhe the linear and nonlinear head equations on metric star graphs is studied in the context of thermal diffusion. The nonlinear Schrödinger equation with the subcritical power nonlinearity on a metric star graph generalized Kirchhoff boundary conditions in dimarecent , where the stability of half-soliton solutions is also analyzed. Integrable sine-Gordon equation on metric graphs is studied in caputo14 ; Our1 ; Karim2018 . Linear and nonlinear systems of PDE on metric graphs are considered in Bolte ; KarimBdG ; KarimNLDE .
Quantum graphs, which are the one- or quasi-one dimensional branched quantum wires, can be modeled in terms of quantum mechanical wave equations on metric graphs by imposing the boundary conditions at the branching points (vertices) and bond ends. The metric graphs are the set of bonds with assigned length and which are connected to each other at the vertices. The connection rule is called topology of a graph and given in terms of the adjacency matrix Uzy1 ; Uzy2 :
[TABLE]
for .
The graph with the simplest topology is called star graph. It has one branching point, which connects three or more bonds (see, Fig. 1). The particle (wave) dynamics in a quantum graph is described in terms of the linear Schrödinger equation (with ), which can be written for each bond of a star graph with bonds as
[TABLE]
where is the wave function of th bond, is the potential acting on th bond and is the bond number. In the following, for simplicity we assume that . To solve Eq. (11) one needs to impose initial condition and boundary conditions at the branching point and bond ends. Most general boundary conditions providing self-adjointness of the Schrödinger operator on graphs have been derived in Kost . Physically reasonable vertex boundary conditions, which are consistent with these general ones can be chosen as the continuity of the wave function at the vertex and Kirchhoff’s rule Uzy1 ; Uzy2 :
[TABLE]
[TABLE]
According to the scattering approach Uzy1 this “natural” boundary conditions, given at a vertex , cause reflection of a wave packet with back-scattering amplitude and the reflection probability approaches 1 as the number of bonds increases. Other details of the solution of the stationary version of Eq. (11) can be found elsewhere (see, e.g. Uzy1 ; Uzy2 ). Here we will focus on the dynamical problem, which is described in terms of the time-dependent Schrödinger equation, by considering the wave transport in quantum graphs.
IV Transparent quantum graphs
Here we combine the above concepts of transparent boundary conditions and quantum graphs to design the vertex boundary conditions providing reflectionless transmission of waves through the graph branching points. In the following such graphs will be called “transparent quantum graphs”. Using such approach, below we show that Kirchhoff-type boundary conditions (weight continuity and current conservation) can become equivalent to transparent vertex boundary conditions under constraints in the form of a simple sum rule that is fulfilled for the weight coefficients. For simplicity (without losing generality) we consider first a star graph with three bonds. We assign to each bond a coordinate , which indicates the position along the bond: for bond it is and for they are . In the following we use the shorthand notation for and it is understood that is the coordinate on the bond to which the component refers. Then for the wave dynamics in such graph one can write the time-dependent Schrödinger equation given by (in the units )
[TABLE]
At the branching point (vertex) we impose the boundary conditions in the form of the wave function weight continuity given by
[TABLE]
and Kirchhoff type conditions, which is given as
[TABLE]
In the following we consider the wave going from the first to second and third bonds, i.e., the initial condition is compactly supported in the first bond. Then the “interior” problem for the first bond, can be written as
[TABLE]
Corresponding “exterior” problems for are given by
[TABLE]
Finally, using the Laplace transformation (4), the two “exterior” problems are transformed to
[TABLE]
The general solution of the exterior problems can be written as
[TABLE]
Since we have , we obtain
[TABLE]
Then
[TABLE]
At the vertex () for bonds (“input” from interior) using the continuity boundary conditions (15) we get
[TABLE]
The Laplace transformed current conservation (at ) takes the form
[TABLE]
Using the inverse transform we have
[TABLE]
with
The boundary condition given by (26) coincides with that in Eq. (9) and thereby providing reflectionless transmission for the bond , when , i.e. the following sum rule is fulfilled:
[TABLE]
Hence, the vertex boundary conditions given by Eqs. (15) and (16) become equivalent to the transparent vertex boundary conditions, provided the sum rule in Eq. (27) is fulfilled.
The numerical example in Fig. 2 shows a simulation of a right traveling Gaussian wave packet
[TABLE]
at four consecutive time steps. The Crank-Nicolson finite difference scheme with the space discretization and the time step has been used. It is clear from this figure that the wave completely transmits to the second and third bonds when time elapses. We note that the above boundary conditions provide conservation of the total norm, which is defined as the sum of partial norms for each bond. Detailed derivation of the norm conservation is provided in Appendix.
Fig. 3 presents plots of the partial and total norms for this case. In Fig. 4 the reflection coefficient determined as the ratio of the partial norm for the first bond to the total norm
[TABLE]
is plotted as a function of for the fixed values of and . As this plot shows, at the value of that provides fulfilling of the sum rule (27) the reflection coefficient becomes zero. This confirms once more that the sum rule in Eq. (27) makes the vertex boundary conditions in Eqs. (15) and (16) equivalent to the transparent ones.
We note that the above approach can be extended for other network branching topologies, e.g., for tree graph presented in Fig. 5. The graph consists of three “layers” , , , where run over the given bonds. On each bond , , we have the time-dependent Schrödinger equation given by (14), for which Kirchhoff-type boundary conditions are imposed at each vertex. Transparent boundary conditions can be imposed similarly to those for the star graph. Then, for all , , the and have to be determined from the sum rule like (27) at each vertex. For instance, at the three vertices in Fig. 5 we have
[TABLE]
Fulfilling of these sum rules turns the Kirchhoff-type boundary conditions equivalent to the transparent ones at the each vertex of a tree graph. The constants in Eq. (27) characterize the physical properties of the material, which is used to construct the network and can have different meaning for each specific system under consideration.
V Conclusions
In this paper we studied the problem of transparent quantum graphs by determining them as branched quantum wires providing reflectionless transmission of waves at the branching points. The boundary conditions for the time-dependent Schrödinger equation on graphs, providing such transmission, are formulated explicitly. A simple constraint that makes the usual Kirchhoff-type boundary conditions at the vertex equivalent to those of transparent ones is derived.
We have shown numerically for the star graph a reflectionless transmission of the Gaussian wave packet through the vertices, if these constraints are fulfilled. The approach can be extended straight forward for arbitrary graph topologies, which contain any subgraph connected to two or more outgoing, semi-infinite bonds. Stability analysis for the boundary conditions discretization scheme is not included into the above treatment and will be subject of forthcoming papers.
The motivation for the study of transparent quantum graphs comes from different practically important problems of optics, condensed matter and polymers. One of the problems allowing direct application the above model is charge carrier transport in conducting polymers, which are highly branched nanoscale structures. Reflectionless transmission of charge carriers in these structures causes high conductivity and optimal functioning of different organic electronic devices. One of such structures was discussed recently in the context of nonlinear charge carriers in branched conducting polymers Chsol . More attractive applications are possible also in linear fiber networks, where the wave dynamics is described in terms of the Helmholtz equation.
Finally, we note that extension of the above approach of other linear wave equations, such as heat, diffusion, Dirac and Klein-Gordon equations should be rather straightforward. Moreover, the method can be modified and applied for the nonlinear partial differential equations widely used in physics.
VI Acknowledgements
This work is supported by the grant of the Ministry for Innovation Development of Uzbekistan (Ref. Nr. BF2-022). The work of JRY is supported by Deutscher Akademischer Austauschdienst (DAAD).
Appendix A Norm conservation
For quantum star graph with three bonds, the total norm is given as
[TABLE]
where is the domain of the bond .
For Eq. (14) an easy calculation shows that the derivative of the norm with respect to reads
[TABLE]
Taking into account that vanishes as and vanish as we obtain
[TABLE]
and using boundary conditions (15) and (16) we get
[TABLE]
which means the preservation of the norm over time.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) R. Wilson, T.J. Karle, I. Moerman, T.F. Krauss, J. Opt. A: Pure Appl. Opt., 5 S 76 (2003).
- 2(2) P.I. Borel, , L.H. Frandsen, A. Harpøth, et. al. Electron. Lett., 41 69 (2005).
- 3(3) S.V. Boriskina, Opt. Express, 15 17371 (2007).
- 4(4) W. Yang, X. Chen, X. Shi, W. Lu, Physica B, 405 1832 (2010).
- 5(5) J.B. Driscoll, R.R. Grote, B. Souhan, Opt. Lett., 38 1854 (2013).
- 6(6) G. Berkolaiko, R. Carlson, S.A. Fulling and P. Kuchment (eds.), Quantum Graphs and Their Applications , Contemporary Mathematics 415, AMS, 2006.
- 7(7) P. Kuchment, ar Xiv preprint ar Xiv:0802.3442, 2008.
- 8(8) A. Arnold and M. Ehrhardt, J. Comput. Phys., 145(2) , 611-638 (1998).
