Translator Disclaimer
22 October 2020 Far-infrared photodetection in graphene nanoribbon heterostructures with black-phosphorus base layers
Author Affiliations +

We propose far-infrared photodetectors with the graphene nanoribbon (GNR) array as the photosensitive element and the black phosphorus (bP) base layer (BL). The operation of these GNR infrared photodetectors (GNR-IPs) is associated with the interband photogeneration of the electron–hole pairs in the GNR array followed by the tunneling injection of either electrons or holes into a wide gap bP BL. The GNR-IP operating principle is akin to that of the unitraveling-carrier photodiodes based on the standard semiconductors. Due to a narrow energy gap in the GNRs, the proposed GNR-IPs can operate in the far-, mid-, and near-infrared spectral ranges. The cut-off photon energy, which is specified by the GNR energy gap (i.e., is dictated by the GNR width), can be in the far-infrared range, being smaller that the energy gap of the bP BL of ΔG  ≃  300  meV. Using the developed device models of the GNR-IPs and the GNR-IP terahertz photomixers, we evaluate their characteristics and predict their potential performance. The speed of the GNR-IP response is determined by rather short times: the photocarrier try-to-escape time and the photocarrier transit time across the BL. Therefore, the GNR-IPs could operate as terahertz photomixers. The excitation of the plasma oscillations in the GNR array might result in a strong resonant photomixing.



The quantization of the energy spectrum in the graphene nanoribbon (GNR) arrays, obtained by partitioning of the graphene layer (GL), results in the transformation of the two-dimensional gapless spectrum of the electrons and holes to the spectrum consisting of the sets of one-dimensional subbands. Such a quantization also leads to the formation of the energy gap in the GNRs.18 The transformation of the energy spectrum brings about substantial modifications of the radiation absorption characteristics of the GNRs compared with the uniform GLs. These modifications include the appearance of the absorption resonances associated with the radiative transitions between the hole sub-bands tops and the electron sub-bands bottoms and the appearance of the absorption “red” edge.6,8 These features of the GNR radiation absorption can be used in different infrared optoelectronic devices. The incorporation of the black-phosphorus (bP) layers, exhibiting the interband absorption of the photons with the energies from Ω300 to 2000 meV depending on the bP layer thickness, opens new prospects.919 The heterostructures, in which GLs or GNR arrays are integrated with the bP layers, can exhibit additional functionalities.2026 This is due to the special GL and GNR and the bP layers’ band alignment (the Dirac point in the GLs and GNRs corresponds to the energy gap in the bP) and relatively narrow energy gap (ΔG300  meV) in the bP layers with a sufficiently large number of atomic sheets.27 The GL- and GNR-heterostructures with the black-arsenic (bAs) layers having even smaller band gaps2832 demonstrate similar properties.

The GNR infrared photodetectors have been (GNR-IPs) proposed and analyzed previously6,8,33 using the interband transitions. In those IPs, the carrier photogeneration leads to an increase in the GNR longitudinal conductivity. In contrast, here, we consider the GNR-based vertical IPs using the radiation interband absorption in the GNR arrays. We focus our analysis on the GNR-IPs with the undoped (weakly doped) GNR array sandwiched between the top cap layer (CL) doped by donors (with the donor surface density Σd) and the bP base layer (BL) doped by acceptors (with the surface density Σa). The BL is lying on the contact layer (substrate) of p+-type (see Fig. 1). The CL can be either the n-bP BL or made of other materials, in particular, hexagonal boron nitride (hBN). An example is the GNR-IP based on n-bP/GNR/p-bP/p+-bP heterostructure or on an n-hBN/GNR/p-bP/p+-bP heterostructure. The GNR array is supplied by the side n-type contact to each GNR, so that the GNR array with the contact forms a fork-like structure as shown in Fig. 1. The donors in the CL and acceptors in the BL can be spatially distributed uniformly across these layers or be concentrated in the delta-doped sheets. The GNR array plays the role of the photoemitter replacing the GL-photoemitter.26 The GNR-IP operation is associated with the emission of the photoexcited holes from the GNR due to tunneling into p-BL contributing to the terminal current.

Fig. 1

Schematic view of the GNR-IP device structure.


For the detection of radiation with Ω>300  meV, one can use more or less standard photodiodes using the interband transitions in bulk bP layers.16 We limit our consideration to the photon energies Ω<ΔG, where ΔG300  meV is the energy gap in the bP BL. The band diagram of the GNR-IP with n-CL and p-BL is shown in Fig. 2. Such a choice of the device structure is associated with a relative smallness of the valence band off-set ΔV at the GNR/bP interface in comparison with the pertinent conduction band off-set ΔC. This implies that the IPs of the n-bP/GNR/p-bP/p+-bP type can have smaller cut-off photon energy (in the far-infrared or even in the terahertz ranges) compared to the IPs of the p-bP/GNR/n-bP/n+-bP type. The GNR-IPs based on heterostructures of the latter type can be considered analogously.

Fig. 2

(a) The GNR-IP band diagram with tunneling of the photoexcited holes under bias voltage and (b) interband electron transitions between the sub-bands (n=1n=1 and Δ<Ω<2Δ) in the GNR valence and conduction bands. The open and opaque circles correspond to the holes and electrons, respectively.


At Σ=ΣdΣa+κBLV/4πeW>0, where κBL and W are the BL dielectric constant and thickness, respectively, and V is the bias voltage, the Fermi level in the GNRs is in the conduction band, i.e., εGNRF>0 [see Fig. 2(a)].

The operation principle of the GNR-IPs resembles that of the unitraveling-carrier photodiodes3436 with a bulk absorbing layer serving as the emitter with the diffusive supply of the carriers injected into the BL.

Since the electric field in the p-BL is high, all the holes injected from the GNR can reach the p+ substrate. As a result, the current is determined solely by the injection from the GNR [as schematically shown in Fig. 1(a)], so that the specific of the hole transport in the p-BL (ballistic or drift) does not affect the steady-state characteristics. However, the hole transit time across the p-BL can be fairly different in the cases of ballistic and drift transport (see Sec. 4). In this regard, the situation is analogous to that in the unitraveling-carrier photodiodes.


Main Equations

The electron and hole, ε(p) and ε+(p), energy spectra in GNRs can be presented as the sets of the one-dimensional sub-bands68 [see Fig. 2(b)]

Eq. (1)

where p=py is the carrier momentum along the GNR (in the y direction), vW108  cm/s is the characteristic carrier velocity in GLs (and GNRs), and n=1,2,3, is the index of the sub-band in the valence band (lower sign) and the conduction band (upper sign). The energy spectrum given by Eq. (1) corresponds to the GNR energy gap Δ=2πvW/d and the carrier effective mass (near the sub-band’s top and bottom) m=Δ/2vW2. At d=15 to 90 nm, in agreement with the experimental data,3 one obtains Δ(30030)  meV and m(0.010.002)m0,7 where m0 is the electron mass in vacuum.

Under the steady-state irradiation by the photon flux IΩ with the photon energy Ω, the average photocurrent density jΩ and the GNR-IP responsivity RΩ=jΩ/ΩIΩ are given by (see also the Appendix)

Eq. (2)

Here, e is the electron charge, βΩ is the ratio of the flux of the incident photons absorbed in the GNR array due to the interband transitions to the flux of the incident photons IΩ (the average photon interband absorption coefficient of the GNR), and ΘΩ is the probability of the photoexcited holes escape from the GNR. As seen from Eq. (2), the photocurrent in the GNR-IPs with the vertical injection is determined not only by the absorption coefficient, but also by the escape probability of the photoexcited carriers (holes) ΘΩ (as in the interband GL-IDs,26 see also references therein). For the normal radiation incidence, the GNR array interband absorption coefficient is given by6,8

Eq. (3)

Here, β=πe2/cκ)π/137κ, c is the speed of light, κ is the refractive index determined by the dielectric constants of the CL and BL, κCL and κBL, respectively, and the factor d/D is the fraction of the GNR-IP area covered by the GNR array (the GNR aspect ratio). The function

Eq. (4)

accounts for the Pauli principle for the interband transitions in the GNR (with the Fermi energy εGNRF counted from the Dirac point at the temperature T), whereas the function

Eq. (5)

describes the probability of the interband transitions in the GNR with the electron and holes states given by Eq. (1) and the Gaussian broadening (or smearing) of these states due the collisions, characterized by γ=Γ/Δ, where Γ/τ with τ being the electron and hole momentum relaxation time.

The quantity ΘΩ in Eq. (2) is the tunneling escape probability of the photoexcited holes having the energy εΩ=(ΩΔ)/2 to escape from the GNR.

We estimate ΘΩ as

Eq. (6)

Here, ηΩ=1Ω/2ΔV, Etunn=42mzΔV3/2/3e is the characteristic tunneling field, mz is the hole effective mass in the BL in the direction perpendicular to the atomic sheets, τrelax is the relaxation time of the photoexcited holes in the GNR due to the recombination and the energy relaxation on acoustic phonons, and τesc is the photoexcited holes’ try-to-escape time. The electric field at the GNR array plane EGNR determining the shape of the tunneling barriers depends on the potential drop V between the collector and the GNR. Considering the BL doping and neglecting the quantum capacitance,37 we arrive at the following equation:

Eq. (7)

Here, Ea=4πeΣa/κBL and Ed=4πeΣd/κCL.

Setting ΔV=100  meV and mz0.28m0,10 we obtain Etunn1.19×106  V/cm. For Σa=2.2×1012  cm2 and κBL=4, we obtain Ea1×106  V/cm.


Spectral Characteristics of the GNR-IP Responsivity and Detectivity

According to Eq. (2), the GNR-IP responsivity RΩ—one of the main IP characteristics—can be presented as

Eq. (8)


Eq. (9)


Figure 3 shows the spectral characteristics of the GNR-IP responsivity calculated for parameters of the GNR-IP structure (Δ, Γ, εGNRF, and T). The following parameters are assumed: β=1.15×102, (κ=4), ΔV=100  meV, τesc/τrelax0.1, Etunn1.19×106  V/cm (mz2.5×1028  g), and EGNR=1102  V/μm. The temperature is assumed to be T=26  meV (T300  K) except for Fig. 3(d). The parameter Γ was assumed to be Γ=1 and 5 meV, which corresponds to the realistic momentum relaxation times τ6.25×1013  s and 1.25×1013  s, respectively. At Δ=100  meV, the range of the Fermi energies εGNRF=25 to 50 meV corresponds to the ΣdΣa(1.8to4.4)×1012  cm2. For the degenerate electron gas in the GNRs (sufficiently large Σ), the Fermi energy in the GNRs εGNRFΔ(Σd2/2)2Σ/Δ (see the Appendix).

Fig. 3

The responsivity RΩ versus photon energy Ω of GNR-IPs with (a) Γ=5  meV and (b) Γ=1  meV at different energy gaps Δ (different GNR widths), (c) Δ=75  meV, Γ=5  meV, and εGNRF=25  meV at different Fermi energies in the GNR εGNRF, and (d) Δ=75  meV, Γ=5  meV, and εGNRF=25  meV at different temperatures T. Dotted lines correspond to the GL-IP responsivity.


At Δ=75  meV, the range of the Fermi energies εGNRF=25 to 50 meV corresponds to ΣdΣa(2.4to5.9)×1012  cm2.

As seen from the comparison of the dependences for GNR-IPs with those for the GL-IP (without the partitioning of the GL into the GNRs, i.e., with Δ=0, and presented by the dotted lines) shown in Figs. 3(a) and 3(b), the former can exhibit much higher RΩ when ΩΔ.

As follows from Figs. 3(c) and 3(d), an increase in the Fermi energy εGNRF and in temperature T leads to a smaller responsivity. This is attributed to stronger restrictions for the radiative interband transitions imposed by the Fermi principle [described by the factor SΩ given by Eq. (4)].

The GNR-IP detectivity is given by

Eq. (10)

where j0 is the dark-current density. As follows from Eq. (10), the shape of the spectral characteristic of the detectivity repeats that for the responsivity (DΩ*RΩ). Since DΩ* depends also on the dark-current density DΩ*1/j0, it is sensitive to the parameters determining the tunneling of equilibrium carriers from the GNR into the BL, first of all, to the energy states near the Fermi energy εGNRF. Considering that the probability of the hole tunneling escape from their Fermi energy can be presented as [compare with Eq. (6)]

Eq. (11)

where η0=1±εGNRF/ΔV (the sign “+” corresponds the situation when the Fermi level is in the conduction band, so that the GNR array is charged predominately by electrons, and the sign “−” corresponds the situation when the Fermi level is in the GNR valence band), in the most interesting case (Σ>0, i.e., the GNR Fermi level is in the conduction band) for the dark-current density we obtain

Eq. (12)


Eq. (13)


Eq. (14)

At εGNRF>0, the exponential factor in Eq. (14), and hence the GNR-IP detectivity can be fairly large. Although, when εGNRF becomes larger than Ω/2, SΩ and D¯Ω* can exhibit a steep drop in line with the pertinent responsivity roll-off seen in Fig. 3(b). This indicated that the selection of the CL and BL parameters, where first of all their doping opens the possibilities of the GNR-IP detectivity enhancement. One needs to stress that the proper choosing of the bias voltage might also be important for the GNR-ID responsivity and detectivity.


Detection of Modulated Radiation

If the incident radiation comprises the steady-state and modulation components, IΩ and δIΩωexp(iωt), respectively, where δIΩω and ω are the amplitude and frequency of the modulation components, considering the Shockley–Ramo theorem,38,39 the ac terminal current density is given by [see Eqs. (29) and (30) in the Appendix]

Eq. (15)


Eq. (16)

which replaces ΘΩ in Eq. (26) and τtr is the hole transit time across the BL. The factor [exp(iωτtr)1]/iωτtr in the right-hand side of Eq. (8) describes the hole transit-time effect (transit delay) in the BL.

As follows from Eqs. (8), (15), and (16), for the responsivity to the modulation signals, RΩω=δjΩω/ΩδIΩω one can arrive at

Eq. (17)


Eq. (18)

is the transit-time factor. Hence,

Eq. (19)

According to Eq. (19), the roll-off of the relative responsivity |RΩω/RΩ| with increasing modulation frequency ω is primarily determined by the transit time τtr if the try-to-escape time τesc is shorter than the transit time of the holes τr

Eq. (20)



Plasmonic Response

Depending on the sign of Σ=ΣdΣaκBLV/4πeW, the one-dimensional electron or hole gases are formed in each GNR. The collective response of these electron or hole systems can lead to pronounced plasmonic phenomena.

The ac current of the photoexcited holes emitted from the GLs causes ac hole injection from the side contacts to the GNRs. The ac current along the GNRs (in the y direction, see Fig. 1) could lead to the excitation of the plasmonic oscillations, i.e., the self-consistent spatial–temporal variation of the electron or hole density and the GNR array potential in the GNR array δφΩω(y)exp(iωt).33,4042 This plasmonic oscillation is the standing wave formed by the plasmonic waves propagated along the GNR in the opposite directions (i.e., with the wave numbers ±ky). As a result, accounting for the plasmonic oscillations in the GNR, the averaged current density can be presented as [see the Appendix, Eq. (37)]

Eq. (21)

Here, 1/τdβΩIΩ is given by Eq. (34) in the Appendix. Equation (21) yields the following equation for the responsivity:

Eq. (22)


Eq. (23)

is the factor describing the contribution of the plasmonic resonances [compare Eqs. (17) and (22)] and ωp=ππe2WΣ/κBLmL2 is the plasmonic frequency of the carriers (holes) induced in the GNRs, where ΣΣaΣd+κBLV/4πeW is the hole density in the GNR array induced by the doping of the CL and BL and the bias voltage, and mΔ/2vW2 is the hole effective mass in GNRs. Considering the relation between m and Δ, we find

Eq. (24)

Setting EGNR=106  V/cm, Δ=100  meV, W=0.1  μm, and L=2.5  μm, we obtain ωp/2π1.4  THz.

If ωτ1, τd1ωp, one obtains Pω1, and Eq. (22) turns to Eq. (17). At the combined resonance Ω/Δ1 and ω/ωp1, the responsivity |RΩω| exhibits a particularly high maximum.

Figure 4 shows examples of the modulation characteristics |RΩω/RΩ| as functions of the modulation frequency ω/2π. These characteristics account for the delay in the photoexcited holes escape from the GNRs, the hole transit-time delay, and the plasmonic effects.

Fig. 4

The responsivity |RΩω| normalized by |RΩ| versus the modulation frequency ω/2π for (a) ωp/2π=1.5  THz and τ=0.625  ps, (b) ωp/2π=1.5  THz and τ=1.875  ps, and (c) ωp/2π=2.5  THz and τ=1.875  ps. Dashed lines correspond to the modulation without the excitation of plasmonic oscillation.


The following parameters are assumed: ωp/2π=1.5 and ωp/2π=2.5  THz, τtr=0.5  ps, τd=10  ps, τesc=τ=0.625  ps and τesc=τ=1.875  ps (corresponding to Γ=/τ=1  meV and Γ=1/3  meV, respectively). The frequency dependences shown in Fig. 4 reflect the interplay of the transit-time and plasmonic effects.

As seen from Fig. 4, |RΩω/RΩ| and hence |RΩω| exhibit the maxima corresponding to the plasmonic resonance ω/ωp. The responsivity |RΩω| becomes particularly large at the combined resonances ΩnΔ/ and ω/ωp1. The resonance effects under consideration are pronounced when the carrier momentum relaxation time τ is sufficiently long.

The inclusion of the plasmonic effects enables the possibility of the GNR-IPs operation in the terahertz range of the modulation frequencies, in particular, as the terahertz photomixers.



Both the realization of the resonant interband intersubband radiative transitions and the possibility of the pronounced plasmonic oscillations require relatively perfect GNRs with high carrier mobilities along the GNRs or their ballistic transport, which are achievable in real GNR-based heterostructures. The electron and hole mobilities on the order of (2000to2500)  cm2/Vs measured in GNRs43,44 are sufficiently large to support resonant interband transitions and plasmonic oscillations.

As seen from Figs. 3(a) and 3(b), the comparison of the vertical GNR-IP and GL-IP responsivities’ dotted lines in [Figs. 3(a) and 3(b)], RΩ and RΩ, at Ω=Δ yields RΩ/RΩ1. The same is valid for RΩω/RΩω. Thus, the GNR-IPs can markedly surpass the GL-IPs in the responsivity if the GNRs are sufficiently perfect with long electron and hole momentum relaxation times, in particular, τ=(1.25to6.25)×1013  s [corresponding to Γ=(1to5)  meV], assumed in Figs. 3(a) and 3(b).

Despite a similarity of the GNR-IP and the standard unitraveling-carrier photodiodes, there are the following distinctions enabling the advantages of the former due to the following: the operation in the far-infrared range because of a narrow energy gap in the GNRs and the BLs, an increase photoexcitation efficiency associated with the resonant intersubband transitions, and a shortened escape time of the photoexcited holes.

Comparing the responsivity RΩω of the vertical GNR-IP under consideration with the responsivity R¯Ωω of the lateral GNR-IP with the p–i–n structure having the depleted GNR array (like those considered in Refs. 8 and 4547) and exploiting the variation of the GNR conductivity under irradiation (with the same lateral sizes), at low- and high-modulation frequencies ω one can find, respectively, RΩ/R¯Ω2 and RΩωR¯Ωω2τ¯tr/τtr. Here, τ¯tr/τtr(L/W)1 is the ratio of the lateral (along the GNRs) and vertical (across the BL) transit times. The responsivity of the lateral GNR-IPs with the n-contacts and the energy barrier formed by the GNR between these contacts6,8,33 can include the photoconductive gain increasing their responsivity, but at the expense of a lower speed operation.

The structure of the GNR-IP can be inverted placing the BL on the device top. Such GNR-IPs can be made on the hBN substrate.

Due to the similarity of the band alignment of the GNR with the bP and bAs layers, the obtained results can be applied to the GNR-IPs with the b-As layers considering the pertinent quantitative difference in ΔC and ΔV.

The resonant response of the GNR-IPs is revealed when the smearing of the carrier spectra is relatively small. This requires a smallness of γ (i.e., sufficiently long momentum relaxation time τ) and a small overlap of the wave functions in the neighboring GNRs. The latter implies that the spacing between these GNRs l=(Dd) should not be too small. This condition looks like exp[22mz(ΔVΔ/2)(Dd)/]1. For the GNR-IPs operating at the photon energies Ω100  meV, d4  nm, the above inequality is satisfied if Dd>0.4  nm. This means that the spacing between the GNRs can be fairly small, so that d/D1 and the potential sag in the inter-GNR region is small. All these indicate that the lateral nonuniformity of the potential in the GNR array plane can be very small.



We proposed and evaluated interband far-infrared photodetectors with the GNR array and the bP BL—GNR-IPs. Due to the specifics of the energy spectra of graphene and bP, the GNR-IPs can operate at photon energies much smaller than 300 meV. Due to a strong interband resonant absorption in the GNRs, the GNR-IPs can surpass the interband GL-IPs and different intersub-band IPs in responsivity and detectivity. The GNR-IPs are analogous to the unitraveling-carrier photodiodes and can operate at fairly high-modulation frequencies of the incident radiation. The possibility of the excitation of plasmonic oscillations in the GNR array of the device might lead to a pronounced response when the radiation modulation frequency is close to the frequency of the plasmonic oscillations, which is typically in the terahertz range. This enables the use of the GNR-IPs as terahertz photomixers. The GNR-IP response can be particularly strong at the combined resonance—the radiation frequency corresponds to the resonant interband transitions between the sub-bands, whereas the radiation modulation frequency corresponds to the plasmonic resonance.




Steady-State Photocurrent

Considering that the incident radiation comprises the steady-state and modulation components, IΩ and δIΩωexp(iωt), respectively, where δIΩω and ω are the amplitude and frequency of the modulation component, the dc and ac components of surface density, ΣΩ and δΣΩω, of the photoexcited carriers in the GNR array averaged in the lateral directions ΣΩ+δΣΩω is governed by the following equations:

Eq. (25)


Eq. (26)

The parameters in Eqs. (25) and (26) were given in the main text.

The dc components of the average current density are given by

Eq. (27)

Equations (25) and (27) yield Eq. (2) with the factors determined by Eqs. (3)–(6).


Photocurrent at Modulated Radiation

The ac component of the average terminal current density is equal to

Eq. (28)

Here, τtr is the transit time of the carriers crossing the BL. This time can be estimated as τtr=W/vS, where vS is the average (or saturation) velocity in the BL. The factor containing τtr describes the effect of transit-time finiteness on the induced terminal current in line with the Shockley–Ramo theorem.38,39 Equation (28) does not account for the possible excitation of plasma oscillations (spatiotemporal variations of the carrier density in the GNR-array and of the electric field in the BL), which could result in the temporal variation of EGNR and in the displacement terminal current. Such a simplification of Eq. (28) is justified when ω is far from the characteristic plasma frequencies.

Equations (26) and (28) yield

Eq. (29)


Eq. (30)



Plasmonic Oscillations in GNR Arrays

The amplitude of the plasmonic mode potential δφΩω is governed by the equation, which is the consequence of the electron transport equation along the GNR and the Poisson equation (similar to that used previously for different plasmonic structures)33,40,41

Eq. (31)

with the boundary conditions δφΩω|y=0=0 (at the contact) and dδφΩω/dy|y=L=0 (at free GNR ends). Here, s=4πe2WΣ/κBLm is the characteristic plasma velocity in the GNR array separated from the highly conducting p+ substrate by the BL of the thickness W,7 Σ is the hole density in the GNR array, m is the effective mass in the GNRs, and L is the length of the GNRs (spacing between the side contacts to the GNRs).

Since the ac current δjΩω apart from the current due to the variation of the radiation intensity includes the component arisen due to the variations of the tunneling probability associated with the change of the electric field at the GNR array δEGNR in the case of the plasmonic oscillations, Eq. (29) should be replaced by

Eq. (32)

so that Eq. (31) becomes

Eq. (33)


Eq. (34)

Solving Eq. (33) with the above boundary conditions, we obtain

Eq. (35)

so that

Eq. (36)

where ωp=πs/2L=ππe2WΣ/κBLmL2 is the plasmonic frequency.

Adding the density of the displacement current δiΩω=iω(κBL/4πW)δφΩω to δjΩω, we find

Eq. (37)



Fermi Energy in the GNR Array

If Σ=ΣdΣaκBLV/4πeW>0, i.e., at a stronger doping of the CL than the BL and at not too high bias voltages, the GNRs are occupied primarily by the induced electrons. At sufficiently large Σ, the electron gas in the GNRs can be degenerate. In this case, the linear electron density in GNRs is given by

Eq. (38)

Here, pF is the electron momentum at the Fermi energy. The average surface electron density is equal Σ=ρ/D, so that

Eq. (39)


At not too high electron densities when the electrons occupy only the lowest sub-band n=1, considering Eq. (1), we have

Eq. (40)


Eq. (41)



The work at RIEC, UoA, IUHFSE, and UB was supported by the Japan Society for Promotion of Science (KAKENHI Grant No. 16H06361), and the RIEC Nation-Wide Collaborative Research (Projects Nos. R03/A12 and H31/A01), Japan. The work at MIPT was supported by the Russian Foundation for Basic Research (Grant Nos. 18-07-01379 and 18-29-02089), Russia. The work at RPI was supported by the US Army Research Laboratory Cooperative Research Agreement (project monitor Dr. Meredith Reed) and the US Office of Naval Research (project monitor Dr. Paul Maki).



C. Berger et al., “Electronic confinement and coherence in patterned epitaxial graphene,” Science, 312 1191 (2006). SCIEAS 0036-8075 Google Scholar


Y.-W. Son, M. L. Cohen and S. G. Louie, “Half-metallic graphene nanoribbons,” Nature, 444 347 (2006). Google Scholar


M. Y. Han et al., “Energy band-gap engineering of graphene nanoribbons,” Phys. Rev. Lett., 98 206805 (2007). PRLTAO 0031-9007 Google Scholar


L. Yang et al., “Quasiparticle energies and band gaps in graphene nanoribbons,” Phys. Rev. Lett., 99 186801 (2007). PRLTAO 0031-9007 Google Scholar


L. Brey and H. A. Fertig, “Elementary electronic excitations in graphene nanoribbons,” Phys. Rev. B, 75 125434 (2007). Google Scholar


V. Ryzhii et al., “Device model for graphene nanoribbon phototransistor,” Appl. Phys. Express, 1 063002 (2008). APEPC4 1882-0778 Google Scholar


V. V. Popov et al., “Oblique terahertz plasmons in graphene nanoribbon arrays,” Phys. Rev. B, 81 073404 (2010). Google Scholar


V. Ryzhii et al., “Terahertz and infrared photodetectors based on multiple graphene layer and nanoribbon structures,” Opto-Electron. Rev., 20 15 –20 (2012). Google Scholar


R. W. Keyes, “The electrical properties of black phosphorous,” Phys. Rev., 92 580 (1953). PHRVAO 0031-899X Google Scholar


A. Morita, “Semiconducting black phosphorus,” Appl. Phys. A, 39 227 (1986). Google Scholar


H. Asahina and A. Morita, “Band structure and optical properties of black phosphorus,” J. Phys. C, 17 1839 (1984). Google Scholar


X. Ling et al., “The renaissance of black phosphorus,” Proc. Natl. Acad. Sci. U. S. A., 112 4523 (2015). Google Scholar


H. Liu et al., “Phosphorene: an unexplored 2D semiconductor with a high hole mobility,” ACS Nano, 8 4033 (2014). ANCAC3 1936-0851 Google Scholar


Y. Xu, J. Dai and X. C. Zeng, “Electron-transport properties of few-layer black phosphorous,” J. Phys. Chem. Lett., 6 1996 (2015). Google Scholar


F. Xia, H. Wang and Y. Jia, “Rediscovering black phosphorous as an anisotropic layered material for optoelectronics and electronics,” Nat. Commun., 5 4458 (2014). NCAOBW 2041-1723 Google Scholar


M. Engel, M. Steiner and P. Avouris, “A black phosphorus photo-detector for multispectral high-resolution imaging,” Nano Lett., 14 6414 (2014). Google Scholar


A. Castellanos-Gomez, “Black phosphorus: narrow gap, wide applications,” J. Phys. Chem. Lett., 6 4280 (2015). JPCLCD 1948-7185 Google Scholar


M. Batmunkh, M. Bat-Erdene and J. G. Shapter, “Phosphorene and phosphorene-based materials—prospects for future applications,” Adv. Mater., 28 8586 (2016). Google Scholar


J. Qiao et al., “High-mobility transport anisotropy and linear dichroism in few-layer black phosphorus,” Nat. Commun., 5 4475 (2014). NCAOBW 2041-1723 Google Scholar


V. Ryzhii et al., “Electrical modulation of terahertz radiation using graphene-phosphorene heterostructures,” Semicond. Sci. Technol., 33 124010 (2018). SSTEET 0268-1242 Google Scholar


V. Ryzhii et al., “Negative terahertz conductivity at vertical carrier injection in a black-arsenic-phosphorus-graphene heterostructure integrated with a light-emitting diode,” IEEE J. Sel. Top. Quantum Electron., 25 2000209 (2019). IJSQEN 1077-260X Google Scholar


V. Ryzhii et al., “Negative terahertz conductivity and amplification of surface plasmons in graphene-black phosphorus injection laser heterostructures,” Phys. Rev. B, 100 115436 (2019). Google Scholar


V. Ryzhii et al., “Real-space-transfer mechanism of negative differential conductivity in gated graphene-phosphorene hybrid structures: phenomenological heating model,” J. Appl. Phys., 124 114501 (2018). JAPIAU 0021-8979 Google Scholar


V. Ryzhii et al., “Negative photoconductivity and hot-carrier bolometric detection of terahertz radiation in graphene-phosphorene hybrid structures,” J. Appl. Phys., 125 151608 (2019). JAPIAU 0021-8979 Google Scholar


M. Y. Morozov et al., “Optical pumping in graphene-based terahertz/far-infrared superluminescent and laser heterostructures with graded-gap black-PxAs1x absorbing-cooling layers,” Opt. Eng., 59 061606 (2019). Google Scholar


V. Ryzhii et al., “Far-infrared photodetectors based on graphene/black-AsP heterostructures,” Opt. Express, 28 2480 –2498 (2020). OPEXFF 1094-4087 Google Scholar


Y. Cai, C. Zhang and Y.-W. Zhang, “Layer-dependent band alignment and work function of few-layer phosphorene,” Sci. Rep., 4 6677 (2014). SRCEC3 2045-2322 Google Scholar


B. Liu et al., “Black arsenic–phosphorus: layered anisotropic infrared semiconductors with highly tunable compositions and properties,” Adv. Mater., 27 4423 –4429 (2015). ADVMEW 0935-9648 Google Scholar


M. Long et al., “Room temperature high-detectivity mid-infrared photodetectors based on black arsenic phosphorus,” Sci. Adv., 3 (6), e1700589 (2017). STAMCV 1468-6996 Google Scholar


S. Yuan et al., “Air-stable room-temperature mid-infrared photodetectors based on hBN/black arsenic phosphorus/hBN heterostructures,” Nano Lett., 18 3172 –3179 (2018). NALEFD 1530-6984 Google Scholar


L. Yu et al., “Electrically tunable optical properties of few-layer black arsenic phosphorus,” Nanotechnology, 29 484001 (2018). NNOTER 0957-4484 Google Scholar


Y. Chen et al., “Black arsenic: a layered semiconductor with extreme in-plane anisotropy,” Adv. Mater., 30 1800754 (2018). Google Scholar


V. Ryzhii et al., “Graphene nanoribbon phototransistor: proposal and analysis,” Jpn. J. Appl. Phys., 48 04C144 (2009). Google Scholar


T. Ishibashi and H. Ito, “Uni-traveling-carrier photodiodes,” 83 –87 Lake Tahoe, California (1997). Google Scholar


H. Ito et al., “Photonic generation of continuous THz wave using uni-traveling-carrier photodiode,” J. Lightwave Technol., 23 4016 –4021 (2005). JLTEDG 0733-8724 Google Scholar


A. Wakatsuki, Y. Muramoto and T. Ishibashi, “Development of terahertz-wave photomixer module using a uni-traveling-carrier photodiode,” (2012). Google Scholar


S. Luryi, “Quantum capacitance devices,” Appl. Phys. Lett., 52 501 –503 (1988). APPLAB 0003-6951 Google Scholar


W. Shockley, “Currents to conductors induced by a moving point charge,” J. Appl. Phys., 9 635 (1938). JAPIAU 0021-8979 Google Scholar


S. Ramo, “Currents induced by electron motion,” Proc. IRE, 27 584 –585 (1939). Google Scholar


V. Ryzhii and M. Shur, “Tunneling- and barrier-injection transit-time mechanisms of terahertz plasma instability in high-electron mobility transistors,” Semicond. Sci. Technol., 17 1168 –1171 (2002). SSTEET 0268-1242 Google Scholar


V. Ryzhii et al., “Plasma mechanism of terahertz photomixing in high-electron mobility transistor under interband photoexcitation,” J. Appl. Phys., 92 5756 –5760 (2002). JAPIAU 0021-8979 Google Scholar


M. Ryzhii, V. Ryzhii and M. S. Shur, “Effect of near-ballistic photoelectron transport on resonant plasma-assisted photomixing in high-electron mobility btransistors,” Semicond. Sci. Technol., 19 S74 –S76 (2004). SSTEET 0268-1242 Google Scholar


Y. Li et al., “Graphene plasmonic devices for terahertz optoelectronics,” Nanophotonics, 9 1901 –1920 (2020). Google Scholar


K. Bang et al., “Effect of ribbon width on electrical transport properties of graphene nanoribbons,” Nano Convergence, 5 7 (2018). Google Scholar


V. Ryzhii et al., “Terahertz and infrared photodetection using p-i-n multiple-graphene-layer structures,” J. Appl. Phys., 107 054512 (2010). JAPIAU 0021-8979 Google Scholar


E. Ahmadi, A. Asgari and K. Ahmadiniar, “The optical responsivity in IR-photodetector based on armchair graphene nanoribbons with p–i–n structure,” Superlattices Microstruct., 52 605 –611 (2012). SUMIEK 0749-6036 Google Scholar


E. Ahmadi and A. Asgari, “Modeling of the infrared photodetector based on multi layer armchair graphene nanoribbons,” J. Appl. Phys., 113 093106 (2013). JAPIAU 0021-8979 Google Scholar


Maxim Ryzhii is currently a senior associate professor at the University of Aizu, Aizu-Wakamatsu, Japan. His research activity includes physics and computer modeling of optoelectronic and terahertz nanostructure devices, and computational modeling of biophysical systems. He is the author or coauthor of more than 180 journal publications and numerous conference papers. He is a senior member of the IEEE and a member of the APS.

Victor Ryzhii is currently a principle researcher at the Institute of Ultra-High Frequency Semiconductor Electronics and a visiting professor at the Research Institute of Electrical Communication, Tohoku University, Japan. He is the author and coauthor of more than 400 research publications in the theory and computer modeling of micro- and nanostructures and electronic and optoelectronic devices and holds 11 patents. He is a corresponding member of the RAS and a fellow of the IEEE and the APS.

Petr P. Maltsev is currently a scientific director at the Institute of Ultra-High Frequency Semiconductor Electronics of RAS, Moscow, Russia. He graduated from Tashkent University of Information Technologies in 1971. Since 1998, he has been a professor at Moscow State Institute of Radio Engineering, Electronics, and Automation. He has more than 250 papers in refereed journals, 40 patents, and 7 monographs.

Dmitry S. Ponomarev works as a deputy director at the Institute of Ultra-High Frequency Semiconductor Electronics of RAS, Moscow, Russia. He is also a leading researcher at A. M. Prokhorov General Physics Institute of RAS and Moscow Institute of Physics and Technology, Moscow, Russia. His research interests are focused on photonic devices, terahertz spectroscopy, and plasmonics.

Vladimir G. Leiman is currently a professor in the Department of General Physics and a principle researcher at the 2-D Material’s Optoelectronics Laboratory, Moscow Institute of Physics and Technology, Russia. His research interests include the theory of intensive electron beams in vacuum and plasma devices and plasmonic phenomena in electronic and optoelectronic devices. He is the author and coauthor of more than 100 journal publications and numerous conference papers. He is a member of the OSA.

Valery E. Karasik is a professor in the Department of Radioelectronics and a director at the Center of Photonics and Infrared Engineering of Bauman Moscow State Technical University. His interests include the research and design of laser devices with the focus on the laser target detection and location. He is a member of the Russian Academy of EII.

Vladimir Mitin is currently SUNY distinguished professor in the Department of Electrical Engineering at the University at Buffalo, the State University of New York. His fields of specialization are microelectronic, nanoelectronic, and optoelectronic devices and materials with special emphasis on sensors and detectors based on nanostructures. He has more than 300 publications including 10 patents, 4 monographs, and 6 textbooks. He is a fellow of the IEEE, the APS, and other professional societies.

Michael S. Shur is the Roberts professor in the Department of Electrical, Computer, and Systems Engineering, Rensselaer Polytechnic Institute, Troy, New York, USA. He is a cofounder of Electronics of the Future, Inc. He is a fellow the IEEE, APS, Electrochemical Society, and a foreign member of the Lithuanian Academy of Sciences and several other professional societies.

Taiichi Otsuji is a professor at the Research Institute of Electrical Communication, Tohoku University, Sendai, Japan. His current research interests include terahertz electronic and photonic materials/devices and their applications. He has authored and coauthored more than 400 research publications. He is a senior member of the Optical Society of America, a fellow of the IEEE and the APS, and several other professional societies.

© The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Maxim Ryzhii, Victor Ryzhii, Petr P. Maltsev, Dmitry S. Ponomarev, Vladimir G. Leiman, Valery E. Karasik, Vladimir Mitin, Michael S. Shur, and Taiichi Otsuji "Far-infrared photodetection in graphene nanoribbon heterostructures with black-phosphorus base layers," Optical Engineering 60(8), 082002 (22 October 2020).
Received: 26 May 2020; Accepted: 19 August 2020; Published: 22 October 2020

Back to Top