🔬 mesoscale physics

Enabling the bulk photovoltaic effect in centrosymmetric materials through an external electric field

This paper presents a practical approach using a static out-of-plane electric field to dress Wannier-interpolated Hamiltonians, demonstrating how this method enables and tunes bulk photovoltaic effects, including the emergence of shift currents in centrosymmetric 2D materials like MoS $_2$ , by explicitly modifying orbital hybridization and band dispersion beyond perturbative treatments.

Original authors: Guilherme J. Inacio, Juan José Esteve-Paredes, Maurício F. C. Martins Quintela, Wendel S. Paz, Juan José Palacios

Published 2026-02-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Guilherme J. Inacio, Juan José Esteve-Paredes, Maurício F. C. Martins Quintela, Wendel S. Paz, Juan José Palacios

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you have a sandwich made of two slices of bread (the layers of a material) with a filling in between. In the world of physics, this is a centrosymmetric material—it's perfectly balanced. If you look at it from the top, it looks exactly the same as if you flipped it upside down. Because of this perfect balance, if you shine a bright light on it, it usually just absorbs the energy or reflects it. It doesn't generate a steady flow of electricity (a "photocurrent") on its own.

This paper is about a clever trick to break that balance and make the sandwich generate electricity using light, simply by applying a gentle "push" from the outside.

Here is the breakdown of their discovery using simple analogies:

1. The Problem: The Perfectly Balanced Seesaw

Think of the material's electrons as kids sitting on a seesaw.

In a normal, balanced material (like a double-layer sandwich): The seesaw is perfectly level. If you push down on one side with light, the other side pushes back equally. Nothing moves forward; there is no net flow of electricity.
In a naturally unbalanced material: The seesaw is already tilted. Light can easily make the kids slide to one side, creating a current.

For a long time, scientists thought you couldn't get electricity from the balanced seesaw without physically breaking the sandwich (changing its structure) or twisting it.

2. The Solution: The "Electric Gate"

The authors found a way to tilt the seesaw without breaking it. They applied a static electric field (a constant, invisible push) from the top, like a gate bias.

The Analogy: Imagine holding a ruler vertically and pressing down on the top of the balanced seesaw. Even though the seesaw is still "structurally" the same, the pressure forces the kids to shift their weight. The balance is broken.
The Result: Now, when you shine light on it, the electrons don't just vibrate in place; they slide in a specific direction, creating a usable electric current. This is called the Bulk Photovoltaic Effect.

3. The Secret Sauce: "Dressing" the Electrons

The paper introduces a new way of doing the math to predict this.

Old Way (The "Weak Push"): Previous scientists treated the electric field as a tiny, temporary nudge. They assumed the material barely changed. This worked for weak effects but failed when the push got strong.
New Way (The "Full Outfit"): The authors treat the electric field as if it is dressing the electrons in a new outfit before the light even hits them. They change the "clothing" (the electronic state) of the material first.
- Why this matters: Just like wearing a heavy coat changes how you move, this "dressing" changes how the electrons interact with each other and how they jump between layers. It allows the material to react strongly to the light, even if the electric push isn't huge.

4. The Three Scenarios Tested

The team tested this on three different types of "sandwiches" (Molybdenum Disulfide, or MoS2):

The Single Slice (Monolayer): It's already a bit unbalanced. The electric push adds a tiny bit more tilt. The effect is small, like adding a grain of sand to a scale.
The Balanced Double Slice (2H Bilayer): This is the star of the show. It was perfectly balanced (no electricity). The electric push broke the symmetry, and suddenly, a strong current appeared. It's like taking a flat, still pond and tilting it just enough to make a river flow.
The Tilted Double Slice (3R Bilayer): This one was already tilted (naturally unbalanced). The electric push could either make the tilt steeper (more current) or flatten it out (less current). It's like a volume knob: you can turn the sound up or down, or even mute it completely if you push the knob the other way.

5. The "Saturation" Limit

The paper also discovered a limit.

The Analogy: Imagine pushing a car. At first, the harder you push, the faster it goes. But eventually, you hit a wall of air resistance. Pushing harder doesn't make it go faster; it just gets exhausting.
The Physics: As the electric field gets very strong, the material's internal structure changes so much that the "current" stops growing and actually starts to drop. The electrons get "stuck" in a new configuration where they can't slide as easily.

Why This Matters

This is a big deal for the future of technology:

New Solar Cells: It means we can take materials that are naturally symmetrical (and usually useless for solar power) and turn them into efficient power generators just by applying a voltage.
Tunable Devices: We can build sensors or switches where we control the electricity flow not by changing the hardware, but by simply turning a "voltage dial."
Better Predictions: The new math method the authors developed allows scientists to accurately predict how these materials will behave in strong fields, which is crucial for designing the next generation of ultra-fast, light-powered computers.

In a nutshell: The authors found a way to use an invisible electric "hand" to tilt a perfectly balanced material, allowing light to push electrons into a steady flow of electricity. They also built a new mathematical "map" to predict exactly how much electricity will flow, helping us design better, smarter solar and light-based devices.

1. Problem Statement

The bulk photovoltaic effect (BPVE), specifically the shift current, is a second-order nonlinear optical phenomenon that typically requires the absence of inversion symmetry (non-centrosymmetric materials). While strategies exist to induce second-order responses in nominally centrosymmetric 2D materials (e.g., strain, stacking engineering, or electric fields), existing theoretical approaches have significant limitations:

Perturbative Limitations: Previous methods often treat the static electric field ( $E_{DC}$ ) only as a perturbation (e.g., via third-order tensors $\sigma^{(3)}$ ). This fails to capture non-perturbative effects such as band-gap renormalization, field-induced hybridization of Bloch states, and the eventual saturation of the response at experimentally realistic field strengths.
Model Limitations: Tight-binding models often neglect off-diagonal position matrix elements (intersite dipoles), treating the field merely as a layer-resolved on-site potential. This omission can artificially suppress the coupling between the static field and electronic states, particularly in planar systems or centrosymmetric multilayers where interlayer hybridization is crucial.

The authors aim to develop a practical, non-perturbative approach to model and predict how a static out-of-plane electric field (gate bias) can induce and tune nonlinear photocurrents in layered materials, including centrosymmetric bilayers like 2H-MoS $_2$ .

2. Methodology

The authors propose a non-perturbative "field-dressing" approach implemented within a Wannier-interpolated Hamiltonian framework.

Hamiltonian Formulation:
- They utilize the length-gauge formulation where the total Hamiltonian includes a static out-of-plane field ( $E_{DC}$ ) and an optical probe field.
- Crucially, the static field is incorporated directly into the electronic ground state before optical excitation: $\hat{H}^{(DC)}_0 = \hat{H}_0 - e E_{DC} \hat{r}_z$ .
- Unlike minimal tight-binding models, they use Wannier functions derived from ab initio (DFT) calculations. This ensures the inclusion of off-diagonal (intersite) position matrix elements ( $\langle \alpha | \hat{r}_z | \beta \rangle$ ). These elements allow the field to modify orbital hybridization and band dispersion, not just the local potential.
Calculation of Shift Conductivity:
- The system's eigenstates are solved for the field-dressed Hamiltonian.
- Using the independent-particle approximation and the Sipe-Shkrebtii formalism, the second-order shift conductivity ( $\sigma^{(2)}$ ) is calculated using the field-dressed eigenvalues and eigenvectors.
- This allows for the evaluation of the full evolution of the response from weak fields (perturbative regime) to strong fields (saturation regime).
Connection to Perturbation Theory:
- The authors perform a Taylor expansion of the field-dressed $\sigma^{(2)}(E_{DC})$ around $E_{DC}=0$ .
- They analytically demonstrate that the first-order term in this expansion recovers the mixed third-order tensor $\sigma^{(3)}_{abcd}(0; \omega, -\omega, 0)$ , establishing a direct link between their non-perturbative method and existing perturbative theories.
Systems Studied:
- MoS $_2$ Monolayer: $D_{3h}$ symmetry (becomes $C_{3v}$ under field).
- 2H-MoS $_2$ Bilayer (AA' stacking): Centrosymmetric $D_{3d}$ symmetry (becomes $C_{3v}$ under field).
- 3R-MoS $_2$ Bilayer (AB stacking): Non-centrosymmetric $C_{3v}$ symmetry (remains $C_{3v}$ , but polarity is tuned).

3. Key Contributions

Non-Perturbative Framework: The development of a method to treat static electric fields non-perturbatively in Wannier Hamiltonians, capturing effects like band-gap renormalization and state hybridization that perturbative methods miss.
Role of Intersite Dipoles: Highlighting that off-diagonal position matrix elements are essential for correctly modeling field-induced symmetry breaking in layered materials, correcting deficiencies in minimal tight-binding models.
Unified Theory: Providing a unified description that connects the linear-in-field response (third-order optics) with the strong-field regime, explaining the physical origin of response saturation.
Activation of Centrosymmetric Materials: Demonstrating a practical route to generate large, tunable shift currents in centrosymmetric bilayers (2H-MoS $_2$ ) via a simple gate bias, effectively turning them into functional photovoltaic materials.

4. Key Results

2H-MoS $_2$ Bilayer (Centrosymmetric):
- At $E_{DC}=0$ , the shift current is zero due to inversion symmetry.
- Applying a perpendicular field breaks inversion symmetry, inducing a finite shift current that grows linearly at low fields ( $|E_{DC}| \lesssim 1$ mV/Å).
- Saturation: At higher fields, the conductivity saturates and eventually decreases. This is attributed to a competition between the growing shift vector and the weakening interband coupling (oscillator strength) as resonant transitions move away from high-symmetry valleys.
- The induced response is significantly stronger than in the monolayer, with in-plane and mixed components dominating.
MoS $_2$ Monolayer:
- The field activates new out-of-plane tensor components (e.g., $\sigma_{zyy}$ , $\sigma_{zzz}$ ) that are forbidden at zero field.
- However, the magnitude of the change is small because the monolayer is thin, limiting out-of-plane charge displacement. In-plane components remain largely insensitive to the field.
3R-MoS $_2$ Bilayer (Non-centrosymmetric):
- Possesses an intrinsic shift current at zero field.
- The external field acts as a tuning parameter rather than an enabler. Depending on the field polarity, it can either reinforce the built-in layer asymmetry (enhancing current) or counteract it (suppressing current).
- A "compensation effect" is observed where a specific field strength cancels the internal potential difference, temporarily suppressing the shift current (mimicking a restoration of mirror symmetry).
Field Dependence:
- The shift conductivity rises linearly at small fields and saturates at higher intensities.
- The band gap decreases linearly with the field (Stark effect), causing a redshift in optical peaks.

5. Significance

Device Applications: This work establishes a theoretical foundation for designing gate-tunable photodetectors and photovoltaic devices using centrosymmetric 2D materials (like 2H-MoS $_2$ ), which are often more stable and easier to fabricate than non-centrosymmetric counterparts.
Theoretical Advancement: It resolves the discrepancy between perturbative models and experimental observations of saturation at high fields, offering a more accurate tool for predicting nonlinear optical responses in strong-field regimes.
Unified Picture: By linking the non-perturbative field dressing to the perturbative third-order response, the paper clarifies the microscopic origins of field-induced nonlinearities, bridging the gap between weak-field linear optics and strong-field physics.
Experimental Relevance: The predicted field strengths for significant effects are within experimentally accessible ranges (below 0.1 mV/Å for relevant trends), suggesting immediate feasibility for experimental verification in van der Waals heterostructures.