Open Access


  • Marco Pala1,
  • Hermann Sellier2,
  • Benoit Hackens3,
  • Frederico Martins3,
  • Vincent Bayot2, 3 and
  • Serge Huant2Email author
Nanoscale Research Letters20127:472

DOI: 10.1186/1556-276X-7-472

Received: 16 July 2012

Accepted: 2 August 2012

Published: 22 August 2012




Braess paradox Mesoscopic physics Congested networks Scanning gate microscopy


Adding a new road to a congested road network can paradoxically lead to a deterioration of the overall traffic situation, i.e., longer trip times for individual road users, or, in reverse, blocking certain streets in a complex road network can surprisingly reduce congestion[1]. This counterintuitive behavior has been known as the Braess paradox[2, 3]. Later extended to networks in classical physics such as electrical or mechanical networks[4, 5], this paradox lies in the fact that adding extra capacity to a congested network can degrade counterintuitively its overall performance.

Known so far in classical networks only, we have recently extended the concept of the Braess paradox to the mesoscopic world[6]. By combining quantum simulations of a model system and scanning gate microscopy[711], we have discovered that an analog of the Braess paradox can occur in mesoscopic electron networks, where transport is governed by quantum mechanics. To explore the possibility of a mesoscopic Braess paradox, we had set up a simple two-path network in the form of a hollow rectangular corral connected to a source and a drain via two openings, with dimensions such that the embedded two-dimensional electron gas (2DEG) is in the ballistic and coherent regimes of electron transport at 4.2 K. The short wires in the initial corral, Figure1a, were narrower than the long wires in order to behave as congested constrictions for propagating electrons (see below). Branching out this basic network by adding a central wire as shown in Figure1a opens an additional path to the electrons. Also, we have used scanning gate microscopy[711] to partially block by local gate effects the electron transmission through this additional path. Doing so should intuitively result in a decreased total current transmitted through the device since one electron path partly loses efficiency, but we counterintuitively found, both numerically and experimentally, that it is exactly the opposite behavior that can actually take place[6].
Figure 1

Evidence for the key role of network congestion in occurrence of mesoscopic Braess paradox. (a, b, c) Network geometries: all networks have a central (additional) branch of 160-nm width and 320-nm wide openings. In (a), W = 140 nm, L = 180 nm. In (b), W = 560 nm, L = 180 nm. In (c), W = 560 nm, L = 500 nm. (d, e, f) The current transmitted through the networks as a function of the tip position scanned along the median lines of the networks (red lines in (a, b, c)). The source-drain voltage applied to the networks is Vds = 1 mV, and the potential applied to the tip is −1 V (see text). Fermi wavelength (λF) = 47 nm, T = 4.2 K.

A key ingredient in the occurrence of classical Braess paradoxes is network congestion. Our previous work was made on a congested mesoscopic network, and it indeed exhibited a marked paradoxical behavior. In this letter, we study numerically in more detail the effect of congestion by simulating three rectangular corrals of different dimensions, i.e., different degrees of congestion. We show that releasing congestion considerably relaxes the paradoxical behavior. Simulations of the spatial distribution of the current density inside the networks for different positions of the local gate help to interpret our predictions in terms of current redistribution inside the network.


Theoretical details

The three simulated networks are shown in Figure1a,b,c. The narrowest network in Figure1a is nearly identical to the network simulated in our previous work[6], apart from slightly larger openings (320 nm instead of 300 nm). Its dimensions are chosen such that the electron flow is congested. Indeed, in a system where electrons can be backscattered solely by the walls defining the structure geometry, a sufficient condition to reach congestion is obtained when the number of conducting modes allowed by internal constrictions is smaller than the number of conducting modes in the external openings, which implies 2 W < W0 , where W and W0 denote the widths of the lateral arms (both of the same width) and of the external openings (of equal widths too), respectively. In turn, increasing W such that 2 W > W0 , as shown in Figure1b, progressively relaxes congestion since all conducting modes injected by the openings can be admitted in the lateral arms. Starting from the network of Figure1b, we will further relax the congestion by increasing the widths L of the horizontal long arms, as shown in Figure1c.

The transport properties of these structures are simulated within an exact numerical approach based on the Keldysh Green’s function formalism. A thermal average is performed around the Fermi energy EF at the temperature T = 4.2 K. We adopt a mesh size of Δx = Δy = 2.5 nm. The Green’s function of the system is computed in the real space representation that allows us to take into account all possible conducting and evanescent modes. Moreover, in order to reduce the computational time and memory requirements, we exploit a recursive algorithm, which is based on the Dyson equation[6, 9].

In this framework, the current densities along the x-axis (transport direction) and the y-axis (transverse direction) between two adjacent nodes read as follows:
J i , i + 1 ; k , k = 4 e h d ω Re H i , i + 1 ; k , k G i + 1 , i ; k , k < ω ,
J i , i ; k , k + 1 = 4 e h d ω Re H i , i ; k + 1 , k G i , i ; k + 1 , k < ω ,

where H i,i';k,k' represents the Hamiltonian discretized on the local basis, and G< i,i';k,k' (ω) is the ‘lesser-than Green’s function’[9] in the real space representation and energy domain.

The tip-induced potential is simulated by considering a point-like gate voltage of −1 V placed at 100 nm above the 2DEG, which corresponds to a lateral extension of ≈ 400 nm for the tip-induced potential perturbation at the 2DEG level.

Results and discussion

The key role of congestion in the network

Figure1d,e,f shows the current flowing through the structures depicted in Figure1a,b,c, respectively, as a function of the tip position scanned along the median lines (red lines). Figure1d shows the occurrence of an analog of the classical Braess paradox in a congested mesoscopic network as a distinctive current peak centered at Ytip = 0 nm. When the tip-induced potential closes the central wire connecting the two openings in Figure1a, the current is counterintuitively increased. However, Figure1e,f shows that as soon as the condition for congestion is relaxed, allowing a larger number of conducting channels to propagate in the region inside the structure, the paradox disappears, and the total current exhibits a maximum when the tip is placed over the two antidots.

In order to microscopically study this behavior, we have simulated in Figure2 the spatial distribution of the absolute value of the current |J| inside the three structures for Ytip = 0 nm and Ytip = −400 nm. When comparing Figure2a and Figure2d for the congested structure, we can notice that the opening of a third central wire connecting the contacts has a twofold effect. The first consequence is to create a direct connection between the source and the drain, which should positively contribute to the total current flowing through the system. The second one is to generate alternative paths that trap electrons in the central region and should promote a longer stay inside the network. We believe that this second effect is the one responsible for the decrease of the total current as long as the third wire is opened. The comparison of Figure2a and Figure2d is indeed very instructive, and in particular, the behavior of the current through the right path paradoxically decreases while the depleting tip moves away. This behavior clearly indicates that the current contribution of trapped electrons around the right antidot compensates partially the initial current. This effect is only partly replicated in the networks of Figure1b,c, whose current redistributions are shown in Figure2b,e and Figure2c,f, respectively. In these cases, the reopening of the third wire, obtained by placing the tip over the antidot, induces a number of new internal paths, which are small compared to the large number of semiclassical trajectories already present in the lateral arms. Therefore, the closing of the central path implies only a small current increase in Figure1e,f around the position Ytip = 0 nm, which is not sufficient to overcome the current at Ytip = −400 nm.
Figure 2

Current redistribution in the mesoscopic networks. All figures depict contour plots of the spatial distribution of the current density. (a, b, c) The tip, marked by a red dot, is positioned above the middle of the networks, i.e., above the center of the additional arm. (d, e, f) The depleting tip is positioned above the center of the left-hand side antidot. Fermi wavelength (λF) = 47 nm, Vds = 1 mV, and T = 4.2 K.

The robustness of the paradox

Finally, in order to test the robustness of our results, we simulated the non-congested structure of Figure1c at different Fermi wavelengths (λF = 57, 47, and 38 nm). This is shown in Figure3. The behavior of the three curves is qualitatively very similar: they present two regions of maximum current when the gated tip is placed over the two antidots, allowing the passage of electrons through the central path, but they also show a local increase in current around Ytip = 0, when the tip closes the central path. This is a signature that the mechanism responsible for the occurrence of the paradox in the congested structure of Figure1a, even if still present, is not predominant with respect to the direct coupling between the two contacts provided by the third wire.
Figure 3

Robustness of the results. The current transmitted through the network of Figure1c as a function of the tip position scanned along the median line (red lines in Figure1c) for three different Fermi wavelengths (λF = 57, 47, and 38 nm). The potential applied to the tip is −1 V, Vds = 1 mV, and T = 4.2 K.


In this letter, we have studied the geometric conditions of mesoscopic networks for the occurrence of a quantum analog of the Braess paradox, known previously for classical systems only. By analyzing the spatial distribution of current density in different structures, we have shown that congested structures are the most suitable geometries to the occurrence of such a counterintuitive phenomenon. This is reminiscent to what is known for the classical paradoxes, in particular, for the historic road-network Braess paradox.



This work has been supported by the French Agence Nationale de la Recherche (MICATEC project), the FRFC (grant no. 2.4.546.08.F) and FNRS (grant no., and by the Belgian Science Policy (Program IAP-6/42). Vincent Bayot acknowledges support from the Grenoble Nanosciences Foundation (Scanning-Gate Nanoelectronics project).

Authors’ Affiliations

Institut Néel, CNRS and Université Joseph Fourier
IMCN/NAPS, UCLouvain, 2 Chemin du Cyclotron


  1. Youn H, Gastner MT, Jeong H: Price of anarchy in transportation networks: efficiency and optimality transport. Phys Rev Lett 2008, 101: 128701.View ArticleGoogle Scholar
  2. Braess D: Uber ein paradoxon der verkehrsplanung. Unternehmensforschung 1968, 12: 258–268.Google Scholar
  3. Braess D, Nagurney A, Wakolbinger T: On a paradox of traffic planning. Transp Sci 2005, 39: 446–450. 10.1287/trsc.1050.0127View ArticleGoogle Scholar
  4. Cohen JE, Horowitz P: Paradoxical behavior of mechanical and electrical networks. Nature 1991, 352: 699–701. 10.1038/352699a0View ArticleGoogle Scholar
  5. Penchina CM, Penchina LJ: The Braess paradox in mechanical, traffic, and other networks. Am J Phys 2003, 71: 479–482. 10.1119/1.1538553View ArticleGoogle Scholar
  6. Pala MG, Baltazar S, Liu P, Sellier H, Hackens B, Martins F, Bayot V, Wallart X, Desplanque L, Huant S: Transport inefficiency in branched-out mesoscopic networks: An analog of the Braess paradox. Phys Rev Lett 2012, 108: 076802.View ArticleGoogle Scholar
  7. Hackens B, Martins F, Ouisse T, Sellier H, Bollaert S, Wallart X, Cappy A, Chevrier J, Bayot V, Huant S: Imaging and controlling electron transport inside a quantum ring. Nature Phys 2006, 2: 826–830. 10.1038/nphys459View ArticleGoogle Scholar
  8. Martins F, Hackens B, Pala MG, Ouisse T, Sellier H, Wallart X, Bollaert S, Cappy A, Chevrier J, Bayot V, Huant S: Imaging electron wave functions inside open quantum rings. Phys Rev Lett 2007, 99: 136807.View ArticleGoogle Scholar
  9. Pala MG, Hackens B, Martins F, Sellier H, Bayot V, Huant S, Ouisse T: Local density of states in mesoscopic samples from scanning gate microscopy. Phys Rev B 2008, 77: 125310.View ArticleGoogle Scholar
  10. Hackens B, Martins F, Faniel S, Dutu CA, Sellier H, Huant S, Pala MG, Desplanque L, Wallart X, Bayot V: Imaging Coulomb islands in a quantum Hall interferometer. Nat Commun 2010, 1: 39.View ArticleGoogle Scholar
  11. Sellier H, Hackens B, Pala MG, Martins F, Baltazar S, Wallart X, Desplanque L, Bayot V, Huant S: On the imaging of electron transport in semiconductor quantum structures by scanning-gate microscopy: successes and limitations. Semicond Sci Technol 2011, 26: 064008. 10.1088/0268-1242/26/6/064008View ArticleGoogle Scholar


© Pala et al.; licensee Springer. 2012
