Fundamentals · Materials and Interfaces

Graphene-Silicon Heterostructures and Charge Separation at Interfaces

Every device that turns energy into electricity has the same critical point: the place where positive and negative charges are pulled apart before they can recombine. In a solar cell it is the junction between two differently doped regions of silicon. When graphene meets silicon, a new kind of junction appears, one whose barrier can be tuned and whose atomically thin top layer reacts to heat, light and electric fields. Such interfaces are attractive for converting weak ambient energy, but they raise a hard problem. A junction separates charge well when a strong, steady input drives it. Ambient inputs are small, varied and fluctuating, often close to the thermal noise of the material itself. Turning such inputs into a net current requires asymmetry, a departure from equilibrium and careful control of losses at every interface.

The Established Approaches

Semiconductor physics provides a mature description of junctions.

The p-n junction. Shockley’s theory explains how a junction between p-type and n-type material forms a depletion region with a built-in electric field that sweeps photogenerated carriers apart [1]. The current-voltage behaviour follows the diode equation, I = I_s · (exp(qV / n·k·T) − 1), with saturation current I_s and ideality factor n [2]. This is the basis of the silicon solar cell.

The Schottky junction. A metal on a semiconductor forms a barrier φ_B that also rectifies [2]. Graphene acts as a tunable semimetal in this role: graphene-on-silicon Schottky junctions work as solar cells [3], and gating graphene shifts its Fermi level and therefore the barrier height, as shown in the graphene barristor [4].

Van der Waals heterostructures. Two-dimensional materials can be stacked layer by layer with clean interfaces, which allows band alignments to be designed rather than accepted [5].

Rectifying small and fast signals. Geometric diodes made from asymmetric graphene structures rectify at very high frequencies without a conventional junction [6]. The theory of Brownian motors explains how spatial asymmetry combined with a drive away from equilibrium turns unbiased fluctuations into directed transport [7]. Measurements on freestanding graphene connected to diodes confirm that fluctuation-driven currents can be modelled precisely, with the laws of thermodynamics holding throughout [8].

The limits are equally clear. A junction at thermal equilibrium delivers no net power. Rectification efficiency drops sharply when the input signal is small compared with the thermal voltage k·T/q, about 26 mV at room temperature. Each interface adds resistance and recombination losses.

The Open Question

Junction models describe one driving input at a time, usually light, under conditions far from the noise floor. They do not describe an interface exposed to several weak inputs at once: thermal fluctuations, electromagnetic fields and mechanical strain acting together, each near the noise level, in a stack of many graphene-silicon interfaces. Which share of the combined input ends up as separated charge, how the asymmetry of the stack sets the direction of the current, and how losses add up across many interfaces are not captured by the single-junction picture. An accounting framework for the whole stack is needed.

A Generalized Approach

A generalized approach to this problem is represented by the Schubart Master Formula, which sums the coupled inputs over the active material and applies one total efficiency:

Schubart Master Formula: P(t) = η · ∫_V Φ_eff(r, t) · σ_eff(E) dV, with r as bold position vector

P(t) = η · ∫_V Φ_eff(r, t) · σ_eff(E) dV

The derivation, notation and open questions are published at pt-schubart.com. The model is introduced in Mathematical Models for Multi-Source Energy Harvesting.

Definitions

  • P(t): instantaneous electrical output, the usable electrical power delivered at time t. Unit: watt (W). Energy over an interval is the time integral of P(t), in joules or watt-hours.
  • η: total conversion efficiency, dimensionless, between 0 and 1. It combines coupling, resonance, rectification and electrical losses.
  • Φ_eff(r, t): effective ambient flux density at position r and time t, the local strength of all ambient inputs acting on the material.
  • σ_eff(E): energy-dependent effective interaction cross-section, how strongly the material couples to an input of energy E.
  • ∫_V: volume integral over the active material.

Together, Φ_eff · σ_eff is the coupled power per unit volume (W/m³), so the integral over V yields watts. How that product divides between the two factors depends on the input channel, so no separate units are assigned to them here.

Assumptions and Boundaries

The Schubart Master Formula is a proposed generalized framework, not an established equation of physics. Its assumptions and limits:

  • The material is an open, non-equilibrium system that exchanges energy with its surroundings.
  • Output is bounded by coupled input power: P_out ≤ ΣP_in.
  • Charge separation and rectification enter through η. At thermal equilibrium, no interface delivers net power.
  • The active volume is concentrated near the interfaces, where built-in fields act; bulk regions far from a junction contribute little.
  • Interface resistance and recombination reduce η and grow with the number of interfaces; they must be measured, not assumed.

Worked Example

All values below are illustrative values, chosen to show the method. They are not measured results.

A stack contains 50 graphene-silicon interfaces, each 10 cm × 10 cm. Assume that the charge-separating region around each interface is 100 nm deep. The active volume per interface is 0.01 m² × 10⁻⁷ m = 10⁻⁹ m³, and the total is V = 5 × 10⁻⁸ m³.

With a coupled power density Φ_eff · σ_eff of 20,000 W/m³ near the interfaces, ΣP_in = 1 mW. With η = 0.15, the output is P = 150 µW, within the boundary (150 µW ≤ 1 mW).

If interface losses grew with the number of interfaces so that η fell to 0.10 for a stack of 100, the output would be 100 × 10⁻⁹ m³ × 20,000 W/m³ × 0.10 = 200 µW: more interfaces, but less than double the output. The example shows why interface quality matters as much as interface count.

Where It Leads

The material behind these interfaces is described in Graphene in Energy Conversion: From Charge Mobility to Output. How stacks become modules, and modules become devices and systems, is set out in From Microscopic Conversion to Grid Scale and From Sensor to Grid.

References

  1. W. Shockley, “The theory of p-n junctions in semiconductors and p-n junction transistors,” Bell System Technical Journal 28, 435–489 (1949).
  2. S. M. Sze and K. K. Ng, Physics of Semiconductor Devices, 3rd edition, Wiley (2007).
  3. X. Li, H. Zhu, K. Wang et al., “Graphene-on-silicon Schottky junction solar cells,” Advanced Materials 22, 2743–2748 (2010).
  4. H. Yang, J. Heo, S. Park et al., “Graphene barristor, a triode device with a gate-controlled Schottky barrier,” Science 336, 1140–1143 (2012).
  5. A. K. Geim and I. V. Grigorieva, “Van der Waals heterostructures,” Nature 499, 419–425 (2013).
  6. Z. Zhu, S. Joshi, S. Grover and G. Moddel, “Graphene geometric diodes for terahertz rectennas,” Journal of Physics D: Applied Physics 46, 185101 (2013).
  7. P. Hänggi and F. Marchesoni, “Artificial Brownian motors: controlling transport on the nanoscale,” Reviews of Modern Physics 81, 387–442 (2009).
  8. P. M. Thibado, J. C. Neu, P. Kumar, S. Singh and L. L. Bonilla, “Charging capacitors from thermal fluctuations using diodes,” Physical Review E 108, 024130 (2023).
  9. Schubart Master Formula: derivation, notation and open questions, pt-schubart.com.

This page is published by the Neutrino Engineering District. Back to Fundamentals