Graphene nanoribbons (GNRs) are promising components in future nanoelectronics due to the large mobility of graphene electrons and their tunable electronic band gap in combination with recent experimental developments of on-surface chemistry strategies for their growth. Here, we explore a prototype 4-terminal semiconducting device formed by two crossed armchair GNRs (AGNRs) using state-of-the-art first-principles transport methods. We analyze in detail the roles of intersection angle, stacking order, inter-GNR separation, GNR width, and finite voltages on the transport characteristics. Interestingly, when the AGNRs intersect at , electrons injected from one terminal can be split into two outgoing waves with a tunable ratio around 50% and with almost negligible back-reflection. The split electron wave is found to propagate partly straight across the intersection region in one ribbon and partly in one direction of the other ribbon, i.e., in analogy with an optical beam splitter. Our simulations further identify realistic conditions for which this semiconducting device can act as a mechanically controllable electronic beam splitter with possible applications in carbon-based quantum electronic circuits and electron optics. We rationalize our findings with a simple model suggesting that electronic beam splitters can generally be realized with crossed GNRs.
I. INTRODUCTION
The wave nature of electrons that propagate coherently in ballistic, one-dimensional conductors has certain qualities in common with photons propagating in vacuum.1 This analogy has spawned the field of electron quantum optics, in which a number of optical setups have been realized in the form of their electronic counterparts, such as the Hanbury Brown and Twiss geometry for studies of Fermion anti-bunching2,3 and the two-particle Aharanov-Bohm effect4–6 as well as Mach–Zehnder interferometry with charged quasiparticles.7,8 The advent of coherent single-particle sources9–13 and entangled electron pair generators14,15 has further provided exciting possibilities for novel quantum technologies and information processing.
A fundamental component for such electron quantum optics is the need for semi-transparent “mirrors,” i.e., electronic beam splitters. Currently, most experiments1 rely on mesoscopic devices based on high-mobility two-dimensional electron gases in the quantum Hall effect regime, in which the electron transport occurs by chiral edge channels that are generally protected against backscattering. In this case, a beam splitter is realized with a quantum point contact that is tuned via electrostatic gates such that only one quantum transport channel transmits with probability T = 0.5. However, a drawback of the technology in the quantum Hall regime is the need for low temperatures and high magnetic fields which severely limits possible applications outside of the laboratory.
Graphene nanoribbons (GNRs)16–18 have some highly desirable properties for their use in molecular-scale electronics devices—they can be designed with specific band gaps19–21 and long defect-free samples can now be fabricated with both armchair (AGNR)22 and zigzag (ZGNR) edge topology23 via on-surface synthesis. However, in the standard bottom-up approach it is difficult to fully explore the GNR electronic properties due to interactions with the metallic substrates used for the synthesis. Very recently, this drawback has been bypassed using synthesis on semi-conducting substrates24,25 and by post-synthesis transfer to an insulating substrate.23 Manipulation of single GNRs has also been demonstrated with scanning probe microscopy,26,27 which opens the possibility to build novel electronic networks with GNRs. Simple 4-terminal tunneling junctions can be fabricated by crossing 1D-structures such as carbon nanotubes28,29 or GNRs.30 Indeed, in the context of electron quantum optics, it was very recently theoretically proposed that two crossed ZGNRs could act as an electronic beam splitter.31
The quantum transport properties of GNR-based devices have been extensively studied with first-principles methods, for instance in the contexts of chemical functionalization,32 optical excitations,33,34 thermoelectrics,35–37 local current-density patterns,38 vibrational excitations,39 and spin-scattering in ZGNRs40–42 and hydrogenated AGNRs.43,44 Various multi-terminal GNR geometries have also been addressed, both in-plane GNR devices45–48 and tunneling junctions formed between GNRs.47,49–52 Finite-bias calculations in a multi-terminal context were pioneered by Saha et al.53 and are becoming increasingly accessible in first-principles transport codes, such as the post-processing tool Gollum54 and the open-source, self-consistent methods of TranSiesta.55,56
In this manuscript, we employ state-of-the-art first-principles methods to study the transport properties of tunneling junctions formed by two crossed AGNRs. Earlier studies have explored similar systems,49,52 but these did not account for the charge redistribution in the junction at finite bias. We analyze in detail the roles of intersection angle, stacking order, inter-GNR separation, GNR width, and finite voltages in this effective 4-terminal device. Interestingly, we discovered that when the two AGNRs cross at an intersection angle , a substantial current can be passed from one ribbon to the other and, more specifically, that electrons injected from one terminal can be split into two outgoing waves with a tunable ratio around 50% and with almost negligible back-reflection. We quantify how this inter-GNR tunneling mechanism depends on the precise atomic arrangement and demonstrate how this enables our device to be tuned and controlled to act as an electronic beam splitter. We further propose a simple model to understand qualitatively the critical role of the intersection angle, which points toward the possibility that electronic beam splitters can be realized with GNRs of different chiralities and widths. We therefore speculate that such GNR-based beam splitters could find applications in electron quantum optics at the nanoscale.
II. METHODOLOGY
A. Multiterminal DFT-NEGF
The calculations presented here were performed using the Siesta/TranSiesta packages55–57 that are based on density functional theory (DFT) and nonequilibrium Green’s functions (NEGF), a combination that is referred to as DFT-NEGF. The TranSiesta code was recently generalized to deal with multi-terminal devices in complex geometries, i.e., to allow any number of electrodes pointing in arbitrary directions.56 Following Saha et al.,53 our multi-terminal system is defined by an expanded scattering region that includes the connections to the electrodes and a central region which is chosen such that any two terminals only interact through it. Each semi-infinite terminal j is assumed to be in thermal equilibrium characterized by a chemical potential μj. The transport properties at the steady state are obtained within the NEGF approach58–60 by the propagator through the scattering region Gr which, at energy E, is given by
with . Here S and H are the scattering region overlap and Hamiltonian matrices, respectively, and the jth lead retarded self-energy that introduces the effect of connecting the jth electrode to the central region. On the one hand, when a bias voltage is applied to an electrode it is assumed that its energy levels are rigidly shifted. Therefore, each electrode j has a chemical potential defined by eV, where EF is the Fermi energy of the combined system in equilibrium, V is the applied bias window (the maximum absolute potential difference between any two terminals), and is a proportionality factor that defines the chemical potential of the jth electrode in terms of V. The central region, on the other hand, will have the charge distribution modified due to the connection to biased electrodes, which is then determined self-consistently within the DFT-NEGF procedure.56
In a multi-terminal setup, it is a non-trivial task to determine the electrostatic potential which solves the Poisson equation and fulfills the boundary conditions imposed by all electrodes. In our calculations, we use the box approximation,56 which consists of reinforcing the potential difference between the electrodes at each self-consistent step. This is done by adding the chemical potential to the periodic solution of the Poisson equation at the region belonging to the jth electrode, and with a redefinition of the common energy reference at each iteration step. The box approximation, particularly when combined with semi-conducting low-dimensional electrodes as in the present case, can potentially create an abrupt behavior of the potential at the boundaries between the electrodes and the central region. However, in the calculations presented here, only a modest charge accumulation occurs at the central-region/lead boundary. Even at the largest applied bias less than , atoms accumulate at each side of the boundary, producing a negligible scattering, as demonstrated by the fact that varying the locations of the central-region/lead boundaries had only a negligible effect on the results.
Once the DFT-NEGF self-consistency is achieved, one can compute the transport properties. The current flowing out of a given electrode j depends on the transmission probabilities of electrons being scattered to any of the other electrodes j′. This is expressed in terms of the multi-terminal Landauer-Büttiker formula61
where
is the jth electrode spectral function,
the level-width function, and the Fermi-Dirac distribution. The factor 2 is due to the spin degeneracy in a spin-less treatment.
Finally, to analyze the electron transport properties of multi-terminal devices in real space, we calculate the so-called bond currents,62 i.e., the amount of current flowing from atom α to β. For scattering states originating from the jth electrode, the bond current is defined as
where [Ea, Eb] characterizes the energy window of interest. A summation over orbitals belonging to atoms α and β is implicit in Eq. (5).
B. Details of calculations
Our device, shown in Fig. 1, is comprised of two infinite H-passivated AGNRs rotated by an angle θ with respect to each other. In the majority of the paper, we considered a crossed structure made by 14-AGNRs (armchair graphene nanoribbons with a width of W = 14 carbon atoms), except for the section where we study the role of the ribbon width. Each GNR in the scattering region thus bridges two semi-infinite electrodes, one on each side of the intersection, i.e., the system has a total of four terminals. All calculations were therefore performed using the vdW density functional of Dion et al.63 with the modified exchange by Klimeš, Bowler, and Michaelides64 since the description of dispersive interactions is crucial to describe interlayer binding and density rearrangement.65 The core electrons were described by non-conserving Troullier-Martins pseudopotentials66 and a double-ζ basis set was used to expand the valence-electron wavefunctions.57 The fineness of the real space grid and the orbital radii were defined using, respectively, a 350 Ry energy cutoff and a 30 meV energy shift.57
First we allowed one ribbon, with axis along the direction, to fully relax using the conjugate gradient method with a force tolerance of 5 meV/Å. The relaxed structure was then duplicated and translated along the direction by 3.34 Å (lowest energy distance67,76,77 for ) to explore the dependence of the transport properties on the other geometrical parameters (e.g., angle, stacking) defining our device. We additionally considered the dependence of the transport properties on small variations of the distance between the ribbons.
Each 14-AGNR consists of 640 atoms in the scattering region and, altogether, the system is described by a total of 9280 orbitals with the chosen basis set. The electrode region j, i.e., where the jth semi-infinite lead is coupled to the system (boxes in Fig. 1(a)), is defined by 64 atoms and is described by a chemical potential . The system configuration (relative position and rotation) is defined by the angle θ between the edge vectors (black arrows at Fig. 1(b)) and the relative position between one reference atom and its replica. In order to uniquely define the different structures, we choose the reference atom as the fifth carbon atom along one edge (red dots at Fig. 1(b)). The system is thus geometrically defined by . In what follows, if only θ is explicitly specified then it is understood that the duplicated ribbon was rotated with respect to the center of mass of the portion of the ribbon in the central region, i.e., that shown in Fig. 1.
In our simulations, we considered an electronic temperature of T = 300 K. For the electrode calculations we used 60 k-points along the periodic direction. A level broadening of eV was considered in the electrodes, while eV was used for the contour integrations over the complex plane.55,56 The self-consistency cycle was stopped when the difference between each element of the density matrix changed by less than 10−6.
III. RESULTS
Throughout the paper we will use intra-GNR to refer to events on the same ribbon (such as the transmission between electrodes belonging to the same ribbon) and inter-GNR for events involving the two different ribbons (such as the transmission from one ribbon to the other). Also, we will refer to the AGNR attached to the electrodes 1 and 2 as GNR12 and, analogously, the ribbon attached to the electrodes 3 and 4 as GNR34.
A. Band structure and zero-bias transmission of the isolated 14-AGNR
A natural starting point for investigating the crossbar system lies in understanding the properties of a single GNR. As seen in Fig. 2(a), our periodic calculations predict the 14-AGNR to be a semiconductor with a small band gap of Eg = 132 meV, consistent with the result expected for a width of type W = 3p + 2.20 This opening of the gap as compared to bulk graphene occurs due to the one-dimensional confinement. A direct correspondence between the band structure and the zero-bias transmission of a pristine 14-AGNR (Fig. 2(b), dashed green line) can be made. As for any one-dimensional pristine structure of atomic-scale cross-section, the GNRs have a step-like transmission T(E) = NT(E), where NT(E) is the number of conductance channels (or, equivalently, number of bands) available at a given energy E. Thus, for the single 14-AGNR around the Fermi level (in a range of approximately eV) we first find a small region of zero transmission, due to the gap, and a plateau of transmission 1 associated with the highest valence (VB) and lowest conduction band (CB), respectively, at larger negative and positive energies. The behavior of the transmission is very symmetric with respect to the middle of the band gap, reflecting the approximate electron-hole symmetry of the band structure in the system.
B. Effect of inter-ribbon voltage
We now start analyzing the effect of the scattering due to the interaction between the ribbons both in the intra- and inter-GNR transport characteristics. In Fig. 2(b) we present the intra-GNR transmissions T12 (blue) and T34 (dashed red line) obtained for the structure as a function of an inter-GNR voltage V14, i.e., the bias is applied to create a potential difference between the two GNRs (, with and ). As mentioned above, the zero-bias transmission of a pristine 14-AGNR (dashed green) serves as a reference.
The main effect observed in the intra-GNR transmissions is a rigid shift of the ribbons’ electronic levels by . At higher energies, far from EF, dips on the intra-GNR transmission are observed, which are related to the increase of the backscattering probability, or reflection function, (for electrode 1 see the orange dashed curves in Fig. 2(b)).
For the device, the intra-GNR transmission is considerably larger than the inter-GNR, as can be seen in Fig. 2(c) (notice the logarithmic scale in this figure, see also the top panels in Fig. 3). With increasing bias, the major effect in the inter-GNR is a widening of the transmission gap around EF (Fig. 2(c)), which is proportional to the energy difference between the position of the CB of GNR12 (whose levels were shifted up in energy by the applied bias) and that of the VB from GNR34 (whose levels were shifted down in energy). When the applied bias achieves the same order of magnitude as the energy gap (Eg = 132 meV), the VB from GNR12 reaches the CB of GNR34, which gives rise to an inter-GNR transmission at EF, as shown in Fig. 2(c) for V14 = 0.15 V. For higher bias, e.g., V14 = 0.5 V, the overlap between the GNR12 VB and the GNR34 CB increases and, as a result, an inter-GNR transmission plateau is formed around EF that widens with the applied voltage.
The inter-GNR transmissions T14 and T23 (as well as T13 and T24) exhibit a very similar behavior, which is due to the high degree of symmetry of the system. This symmetry becomes even more evident for devices with and, therefore, we show only the inter-GNR transmissions T13 and T14 from here on. Note, however, that they are not exactly equivalent because the 14-AGNR (see inset to Fig. 2(a)) does not possess mirror symmetry along the axis defining its extended direction.
C. Role of intersection angle
An interesting phenomenon is discovered when one varies the intersection angle θ between the GNRs. In Fig. 3(a), we show the zero bias intra-GNR transmissions T12 and T34 calculated for different angles , , , , and . Again the pristine 14-AGNR transmission (dashed green) is included for reference. Essentially, one observes an overall reduction of the intra-GNR transmission with decrease θ. The lowest transmission values close to EF were obtained with the structure, exactly where one finds a closer matching between the honeycomb lattice of both ribbons in the crossing region. This decrease of the intra-GNR transmission with θ also translates into the opposite behavior of the inter-GNR transmission T14, which tends to increase (Fig. 3(b)). The effect is particularly dramatic for the case, where one finds that of one inter-GNR transmission channel is open in the energy window eV. Surprisingly, the devices with angles only slightly away from θ = 60° ( and ) exhibit a T14 that is at least one order of magnitude smaller in the mentioned energy range. An additional interesting observation is that the device with exhibits a larger inter-GNR transmission T14 for the VB than for the CB, while for the situation is reversed.
One important property observed for all considered rotation angles is the low reflection probability around the Fermi energy (see, for instance, the electrode 1 reflection function R1(E) in Fig. 3(a)), indicating that in the absence of external potential low energy electrons can propagate with negligible backscattering.
In Fig. 3(c) we present the calculated current I14 flowing from electrode 1 to 4 as a function of an inter-GNR applied bias V14. The structure stands out when compared to all other cases, showing an inter-GNR current I14 higher by one order of magnitude, in accordance with what one would expect from the zero bias transmission analysis. The red arrow in Fig. 3(c) indicates the onset of the inter-GNR current at meV (the non-zero values of I14 below the onset observed for 60° are attributed to the small broadening used in the calculations). We note that our results for are at variance with the negative differential resistance reported in Ref. 49 for the same system beyond V14 = 0.2 V. Within our treatment (going beyond some of the approximations applied in that reference) no such feature was found, even at elevated voltages. This suggests that an accurate treatment of the electrostatics and finite-bias self-consistency is an important ingredient in the description of the device characteristics.
At , the inter-GNR transport proves to be more than just a secondary effect. Rather it is as significant as the direct intra-GNR transport. Moreover, comparing the inter-GNR transmissions T13 and T14 (Fig. 3(b)), one can predict that the scattering states from electrode 1 that are transmitted to the crossing GNR will propagate most likely towards the electrode 4 rather than 3 for all , an effect that is most evident for .
This prognosis is confirmed with the bond currents from electrode 1 calculated with an inter-GNR voltage of V14 = 0.5 V (, with and ) and integrated over the energy window eV. On the one hand, with a setup (Fig. 4(a)), all scattering states from electrode 1 almost fully propagate towards terminal 2 and essentially no current flows to the crossing ribbon.68 On the other hand, for (Fig. 4(b)) only about half of the states propagate towards terminal 2, while the other half is transmitted through the crossing to electrode 4, and no current flows from 1 to 3.
D. Operation of one GNR as a gate electrode
In this section, we study the transport characteristics in the crossing as a function of inter- and intra-GNR voltages. The intra-GNR voltage V12 was applied only among the electrodes 1 and 2, i.e., , while the electrodes 3 and 4 were maintained at the same chemical potential, . The inter-GNR voltage V14 was defined by the difference between the chemical potentials from electrode 1 and 4, i.e., . Therefore, within this setup one can investigate how the GNR34 can act as a gate to the current flowing through GNR12 in crossing systems presenting low inter-GNR transmission, such as . Moreover, this allows one to tune the current splitting on devices with higher inter-GNR transmission, which is the case for .
In Fig. 5 we present the different components for the current with the variation of V12 and V14, for both and devices. The intra-GNR current I12 presents an onset at meV and is more sensitive for (top left panel in Fig. 5) where it clearly increases fast with V12 but slower with V14, indicating that the GNR34 produces only a weak gating effect on the current flowing through GNR12. We note that this weak gating effect is also in contrast with the calculations reported in Ref. 49, showing a current variation of several orders of magnitude with the inter-GNR bias.
The inter-GNR current components I13, I14, I23, and I24 for 90° (left panels in Fig. 5) are all negligible compared to the intra-GNR I12, and essentially no change is observed within the applied bias range. For the device (right panels in Fig. 5) the inter-GNR currents I14 and I23 present the same order of magnitude for V12 = 0. When a finite intra-GNR voltage is applied, the main effect observed is that the current flowing to GNR34 arises more from electrode 1 and less from 2, meaning that the electron splitting can be tuned combining V12 and V14.
E. Analysis of the scattering potential at the crossing
In order to characterize the change of the scattering properties in the crossing region, we present here results for the distribution of the electrostatic potential in the central region. This is defined as the Hartree potential plus the local pseudopotential describing the electron-ion interaction within the Siesta/TranSiesta packages.55–57 Figs. 6(a)–6(e) show the electrostatic potential at the middle plane between the two ribbons resulting from our calculations at different angles θ and without applied voltage (). For all angles the potential is higher (more repulsive for electrons) in the crossing region, and with the highest values in regions where the lattices of the two ribbons match. Accordingly, for the case, where the lattices happen to match within the entire intersection region, the potential reveals “bumps” over the entire crossing, which might be interpreted as a source of larger scattering (and, thus, a harder barrier) for propagating electrons. Thus, one might be tempted to assign to this larger corrugation of the effective electron potential the simultaneous decrease of the intra-GNR and increase of the inter-GNR scattering at .
Fig. 6(f) explores whether this scattering effect can strongly modify the current at finite bias. In this figure we show the difference of the electrostatic potential for an inter-GNR voltage V14 = 0.5 V (, with and ) and zero bias () calculations for the device. This plot only reveals smooth changes in the self-consistent electrostatic potential due to the applied bias. In particular, we do not find noticeable changes in the crossing region. This indicates that the electron scattering at the crossing will not be drastically modified by the inter-GNR bias, in agreement with the general trends observed for the transmission functions presented so far.
F. Role of inter-ribbon distance
The above analysis of the electrostatic potential in the crossing suggests that the overlap of carbon π-orbitals may produce a strong effect on the potential distribution and, thus, on the scattering at the crossing, increasing the inter-GNR transmission. In order to test this hypothesis and to have a better understanding of the observed transport properties, we now consider the role of the inter-ribbon distance on the ratio between the intra- and inter-GNR transport.
In Fig. 7 is presented the electrode 1 zero bias reflection and transmissions as a function of the inter-ribbon distance around our reference value c = 3.34 Å.67,76,77 We analyze two extreme cases, namely (the case with higher intra-GNR and lower inter-GNR transmissions) and (with lower intra-GNR and higher inter-GNR transmissions). Inside the varying interval of Å, almost no change is observed in the transmission for the 90° device close to EF. For higher energies, eV, the decreasing distance between the GNRs infers a stronger scattering effect, which is expressed in terms of the reflection probability R1. The dependence with the distance is significantly different for the 60° case (Fig. 7, on the right). As the distance between GNRs decreases, we observe a clear increase of the inter-GNR transmission. This takes place at the expense of the intra-GNR transport, which gets drastically reduced. We note that this result indicates that the transmission in a device could be tuned by applying an external force to the junction. The feasibility of this kind of electromechanical switching has been also suggested for crossed carbon nanotubes.29
G. Role of lattice matching
To investigate the role of lattice matching in the crossing region, another set of calculations was performed by translating one GNR with respect to the other while keeping the inter-GNR distance fixed at c = 3.34 Å. In Fig. 8 electrode 1 zero bias reflection and transmissions for structures with and are presented for six different stacking configurations each. Very little difference is observed in the transmissions among the translated structures with 90° (left panels in Fig. 8). This is consistent with the idea of the overlap between the π orbitals being the key parameter, since the average overlap does not vary much when translating the GNRs with . In other words, an average of the different stackings between carbon atoms in the two ribbons is always sampled when the ribbons cross at and, therefore, a small shift of the ribbons’ positions does not qualitatively change the situation.
In contrast, for the 60° case the results show a strong change already close to EF (right panels in Fig. 8). Two particular stackings can be highlighted: AA (where the carbons from the different ribbons lie on top of each other in the crossing region) and AB (when half of the carbons lie on top of other carbons, while half of them reside on the center of the other GNR hexagons and, therefore, do not overlap). These two stackings correspond exactly to the maximum (AA) and minimum (AB) possible overlap between the carbon atoms in the crossing. Accordingly, the AA case presents the highest/lowest inter-/intra-GNR transmission, while for the AB stacking one finds one of the lowest/highest inter-/intra-GNR transmissions.
H. Role of GNR width
AGNRs of different width can be classified into three distinct family structures obeying the following hierarchy for the gap size with characterizing the number of carbon atoms across the ribbon.20,21 Within a given family the gap size is inversely proportional to the width. Given such hierarchy, one natural question is whether these three categories share similar transport characteristics for the crossed four-terminal structures.
The 14-AGNR considered up to this point belongs to the category 3p + 2 with smallest energy gaps. Therefore, in order to explore the role of the ribbon width we selected the comparable 13-AGNR (3p + 1) and 12-AGNR (3p) as representatives of the other families. The calculated electrode 1 reflection and transmission functions at zero bias of these two new devices are shown in Figs. 9(b) and 9(c) for rotation angles of 90° and 60°. To facilitate a comparison, the corresponding results for 14-AGNR are shown again in Fig. 9(a).
Comparing first the single pristine ribbon transmission curves (dashed green lines) for the different considered widths, one notes that indeed the energy gap follows the hierarchy mentioned above. In all devices with the intra-GNR transmission T12 shows little scattering effects at lower energies and, therefore, almost no reflection at electrode 1 as well as no inter-GNR transmissions T13 and T14 are observed within an energy window eV. On the other hand, for a significant decrease in T12 is observed in all cases.
The transmission properties of the crossed 12-AGNR are rather similar to the 14-AGNR device previously discussed. For the intra-GNR transmission T12 shows a reduction of about 50% when compared to the same structure at 90° for energies between the second sub-bands, followed by an increase of the same order in the inter-GNR transmission T14, while T13 presents a negligible value. However, for the 12-AGNR case a higher backscattering is noticed on R1 at energies right below EF when compared to the 14-AGNR device.
Some qualitative changes are found for the crossed 13-AGNR compared to the other two cases. The overall effect of significant inter-GNR transmission at 60° is preserved, but compared to the other widths, the crossed 13-AGNR presents a larger degree of backscattering already at lower energies. In the narrow energy interval of eV ( eV) corresponding to the onset of the VB (CB) until the onset of VB − 1 (CB + 1) almost no inter-GNR transmission is observed and most electrons are scattered back to electrode 1. The inter-GNR transmission becomes relevant for energies lower (higher) than the VB − 1 (CB + 1) onset.
Finally, besides the transport characteristics of devices with ribbons belonging to different families, another question is whether the inter-GNR transmission at 60° survives also in the limit of very narrow ribbons. One case that can be fabricated with high control nowadays is the 7-AGNR,22 with a width of just nm. The results obtained from a zero bias calculation of such a crossed 7-AGNR structure are presented in Fig. 9(d). Similar to the 13-AGNR case (both belong to the 3p + 1 family), a large backscattering is observed in the energy interval between the onsets of VB and VB − 1 as well as between CB and CB + 1. However, for energies below (above) the onset of VB − 1 (CB + 1) little backscattering is observed and the inter-GNR transmission T14 becomes close to the desired 50%.
IV. DISCUSSION
A. Quantifying lattice matching with registry index
The results presented in Sec. III support the hypothesis that the overlap between π-orbitals in the intersection determines the scattering properties in our system. However, if this simple picture would be enough to explain all the observed phenomena, one could in principle quantify the amount of scattering using some measure that characterizes the overlap for each structure. The registry index (RI) can be used to provide such a measure (see Ref. 69 for a detailed review). The idea is to consider a circle around each carbon atom belonging to the crossing region and compute the overlapping area SCC between the circles from different GNRs. For graphene-like materials, it has been shown that the ideal circle radius to be considered corresponds to half of a C–C covalent bond length in graphene (i.e., 0.71 Å). The RI is then defined as , where SAA and SAB are, respectively, the maximum and minimum possible overlaps between the GNR orbitals.
Table I shows the RI calculated for structures with different intersection angles together with the special cases of with AA and AB stacking. Considering all the atoms in the crossing region (RItot), the highest value is obtained for with AA stacking (RItot = 1.0) and the minimum one for with AB stacking (RItot = 0.0), as one would expect from the definition above. These values qualitatively describe the changes of the inter-GNR currents among the 60° cases [A versus A]. Moreover, the calculated RItot exhibits a qualitative agreement with the transport properties for , , , and devices. However, a discrepancy occurs when all cases are contemplated, since all devices with , including the AB stacking with RItot = 0.0, present higher inter-GNR current than all the other devices with different θ, for which RItot > 0.0.
Device . | 90° . | 80° . | 70° . | 60° . | 50° . | 60° AA . | 60° AB . |
---|---|---|---|---|---|---|---|
RItot | 0.08 | 0.11 | 0.17 | 0.45 | 0.44 | 1.00 | 0.00 |
RIedge | 0.00 | 0.24 | 0.35 | 0.12 | 0.55 | 1.00 | 0.15 |
Device . | 90° . | 80° . | 70° . | 60° . | 50° . | 60° AA . | 60° AB . |
---|---|---|---|---|---|---|---|
RItot | 0.08 | 0.11 | 0.17 | 0.45 | 0.44 | 1.00 | 0.00 |
RIedge | 0.00 | 0.24 | 0.35 | 0.12 | 0.55 | 1.00 | 0.15 |
One could consider that, for example, only the GNR edges in the crossing are relevant for describing the scattering properties. Hence, the registry index can be calculated only considering the overlaps from carbons belonging to the GNR edges (RIedge). This will change the RI to quantitative different values, see Table I. However, none of the registry index values do qualitatively describe the angle dependency on the current when all cases are taken into account.
B. Simple model for inter-ribbon tunneling
The analysis in Sec. IV A shows that the overlap of π orbitals in the crossing region cannot alone account for the observed inter-GNR transmission and, in particular, explain what makes a device with a very interesting and effective candidate as a beam splitter.
The key to understand the physical origin of this important effect is to consider the tunneling probability between the relevant electronic states in each of the two crossing GNRs. To do this properly, it is necessary to take into account not only the overlaps between atomic orbitals in neighboring structures but also the relative phases and amplitudes with which these orbitals participate in those scattering states. The VB and CB states of the 14-AGNR calculated at the Γ point with DFT are shown in Figs. 10(a) and 10(b). These two states are representatives for the available electron bands in the energy range of eV around EF, and share a common structure with four nodal planes along the ribbon. Reminiscent of the electron states at the Dirac point in graphene,70 these states reflect the characteristic ratio
between the electron momentum along the ribbon axis () and perpendicular to it (). The AGNR states can be interpreted as quantized states of a graphene layer and, thus, they can be qualitatively represented by a propagating wave along the ribbon axis with a Bloch wavevector together with a particle-in-a-box state corresponding to a wavevector in the confined (perpendicular) direction.71 Ignoring the details of the wavefunction inside the graphene unit cell, and just focusing on the envelope wavefunction, we just approximate this situation using plane-waves, in which case we are left with the traditional bands of a free-electron wire, i.e., we have
where depends on the GNR width W and the quantum number n (positive integer). As can be seen in Fig. 10(c), the description of the scattering states of the AGNRs using this simple model is qualitatively correct.
The second AGNR is described similarly, but with rotated wavefunctions where the wavevectors in the two ribbons are related via the rotation matrix defined for a clockwise rotation angle θ,
In the spirit of perturbation theory, the inter-GNR tunneling probability49,52,72 is assumed to be proportional to the modulus square of the overlap between the two wavefunctions,
These overlaps can readily be evaluated numerically as shown in Fig. 10(d) as a function of , θ, and the ratio . For the fundamental mode n = 1, the two wavefunctions have a significant overlap in a large part of the parameter space. In the limit of (), where the two GNRs are aligned in parallel (antiparallel), the overlap goes to infinity because of the diverging integration area. As n increases, the region with a significant overlap shrinks towards one universal curve. This situation corresponds to the wavevector matching condition
which in turn yields the relationship
According to Eq. (6) this simply corresponds to , i.e., the exact condition for a maximal inter-GNR tunneling as found in our simulations for AGNR.73
The meaning of condition Eq. (11) can be easily understood. In an elastic scattering process between two extended layers interacting weakly, we expect both the energy and the momentum of electrons to be conserved. In the present case, due to the finite cross-section of the ribbons, the momentum conservation condition will be relaxed for the scattering between ribbons crossing with arbitrary angles. Thus, it is only an approximate condition that becomes more relevant the wider the ribbons under consideration, i.e., the larger the crossing area. Thus, Eq. (11) is just reminiscent of the momentum conservation condition in the case of inter-GNR elastic scattering events.
Notice that this simplified model cannot account for the dependence of the current on the stacking of the GNRs. In order to do so, in addition to the phases carried by the envelope wavefunctions it is necessary to account for the structure of the wavefunctions inside the graphene unit cell (and take into account for the overlaps between π orbitals and their relative phases within the unit cell). This explains the partial success of the RI that allows rationalizing the changes of the inter-GNR transport for a fixed angle θ. However, the main effect of the rotation angle can be accounted for by our simplified model based on a description of the electronic states as plane waves.
Although Eq. (7) can be a good approximation for the CB and VB of AGNRs, the situation is more complicated for nanoribbons of different chiralities and, in particular, for ZGNRs.31,70,71,74 While qualitatively correct at intermediate energies, the quantized graphene bands fail to describe important features of the low energy spectrum of ZGNRs, as the appearance of the edge states at the Fermi level. Therefore, we would need to consider states sufficiently far from the Fermi energy in order to apply a simple quantization model to ZGNRs. In doing so, the relation between parallel and perpendicular momentum must be reversed for small energy states with respect to that presented in Eq. (6) and we obtain
Thus, for ZGNRs Eq. (12) gives rise to a maximum tunneling probability for , in agreement with the results reported in Ref. 31 using a π-orbital tight-biding model of the system.
Our results present clear connections with previous work investigating the modifications of bilayer graphene band structure as a function of the rotation angle of the two layers.75 The above argumentation explains why the alignment of the honeycomb lattices of the AGNR ribbons at intersection radically increases the inter-ribbon interaction regardless of stacking, why this also happens for ZGNRs, although there the electron scattering takes place preferentially at 120°. It also explains the higher inter-GNR conductance at 30° and 90° reported in Ref. 47 for crossed AGNR/ZGNR devices.
V. CONCLUSIONS
In this paper we studied the electronic and transport properties of a 4-terminal junction defined by two crossed AGNRs from first-principles with the Siesta/TranSiesta codes.55–57 Our research comprises a detailed investigation of the system behavior under the variation of structural parameters, such as intersection angle, inter-GNR distance, stacking order, and GNR width, as well as its response under nonequilibrium conditions by considering different setups with a finite voltage applied between the electrodes.
First we focused on the analysis of a crossed device constituted by two 14-AGNRs. Varying the intersection angle between the GNRs we found two extreme cases, namely and with low and high inter-GNR transmission, respectively. Remarkably, for the 60° case we found that the inter-GNR transmission channel is close to 50% and the reflection negligible over a relatively large energy window of eV around the Fermi energy without an applied voltage. Moreover, for all considered cases with the majority of inter-transmitted electrons propagate only in one direction in the other ribbon. Those findings indicate that semiconducting crossed AGNR structures are interesting candidates to be incorporated in quantum electronics devices. In particular, we showed that a system with can operate as an electronic beam splitter where the ratio of intra/inter transmission can be tuned by changing the inter-GNR distance, i.e., it can be mechanically controlled by applying an external force to the junction.
Interestingly, we have shown that this result is not unique for the 14-AGNR, but the overall effect of a significant inter-GNR transmission at 60° is preserved for crossed structures formed by AGNRs of other widths, namely 13, 12, and 7-AGNRs.
We also explored how the crossed structures behave with biased electrodes. Applying an inter-GNR bias voltage, the 60° configuration is again distinguished with an inter-GNR current one order of magnitude higher than all the other considered intersection angles. When one AGNR is subjected to an intra-GNR bias voltage, changing the inter-GNR voltage produces only a weak gating effect on the 90° devices, but reveals the possibility of tuning the current splitting on the 60° case.
Analyzing the electrostatic potential, we found that the lattice matching on the crossing region plays an important role on the scattering properties. Indeed, a significant change of the transmission probabilities is observed by varying the stacking order on the 60° device. These results suggest that the overlap of carbon π orbitals is another essential parameter to the scattering process. The structures’ registry indices indicate that the amount of π-orbital overlap in the crossing can qualitatively describe the changes in the inter-GNR currents among 60° cases as well as explain the trend in the transport properties for structures with . However, the registry index does not describe the inter-GNR transmission in a general fashion. To this extent, we presented a simple model based on a description of the electronic states as plane waves that captures the effect of the angle, which in turn can be related to the approximate conservation of electron momentum during inter-GNR elastic scattering events. We show how our model explains the role of the intersection angle in crossed GNRs of different chiralities.
The emerging picture from the combination of Ref. 31 and the work presented here is that GNRs with different chiralities and widths may be combined in nanoscale crossbar junctions which should allow, under suitable control of the intersection angle, to construct effective and tunable electronic beam splitters.
ACKNOWLEDGMENTS
The authors acknowledge financial support from FP7 FET-ICT “Planar Atomic and Molecular Scale devices” (PAMS) project (funded by the European Commission under Contract No. 610446), the Spanish Ministerio de Economia y Competitividad (MINECO) (Grant No. MAT2013-46593-C6-2-P), the Basque Departamento de Educación, and the UPV/EHU (Grant No. IT-756-13). N.P. acknowledges financial support from No. EU H2020 Project No. 676598, “MaX: Materials at the eXascale” Center of Excellence in Supercomputing Applications.
REFERENCES
Using the same functional and basis set, we found that for bilayer graphene the equilibrium distance c between the two layers is 3.486 Å for AA stacking and 3.294 Å for AB stacking. For a crossed GNR system with one GNR rotated by 90° with respect to the other, we found 3.339 Å as the lowest energy distance, which lies in between these two extreme values for a bilayer graphene. Moreover, the lattice parameter calculated for graphite with different stackings ( Å, Å and Å), which are in close agreement with experimental measurements of c = 3.354 Å,76,77 also presents an interlayer distance between the bilayer graphene values for AA and AB stacking. Therefore, we found it reasonable to use the distance c = 3.339 Å as the reference value for all our crossed structures.
For a better visualization of the bond currents, a cutoff was applied so that all bond currents below a given threshold are not printed.
Note that another matching condition is given by the relation , from where one obtains a maximum for , which taking into account the minus sign from corresponds to the same preferential scattering at 60°.