🔬 mesoscale physics

Charge state and hydration of vacancies in silicon and silicon carbide in contact with water: A grand-canonical first-principles stability map

This study employs grand-canonical first-principles calculations to map the stability of silicon and silicon carbide vacancies in contact with water, revealing that their charge states and hydration levels are jointly determined by pH and electrode potential, with hydrogen binding cooperatively to silicon vacancies but failing to passivate carbon vacancies without cathodic bias.

Original authors: Abhirup Patra

Published 2026-10-27

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Abhirup Patra

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

Inside the solid materials that power our modern electronics, tiny missing atoms play a massive role. These missing spots, known as vacancies, act like empty seats in a crowded theater, disrupting the flow of electricity and light. For decades, scientists have understood how these vacancies behave in a dry, controlled environment, describing their electrical state based on a single factor: the energy level of the electrons surrounding them. However, the real world is rarely dry. Many of these materials, including silicon and silicon carbide, are frequently exposed to water, whether during the manufacturing process, in humid air, or when used in chemical sensors. When water touches these crystals, the missing atoms do not just sit empty; they begin to interact with the water molecules, grabbing hydrogen and oxygen atoms to fill the void. This interaction changes the vacancy's electrical personality entirely, turning a simple empty spot into a complex chemical structure. Understanding this transformation is critical because it dictates how long electronic devices last, how efficiently they convert light into electricity, and whether they can survive in wet environments without corroding.

A researcher has now mapped out exactly how these vacancies change when they meet water, creating a detailed guide that predicts their behavior under different conditions. Instead of looking at the vacancy in isolation, the scientist treated it as a traveler moving through a landscape defined by two main forces: the electrical potential of the material and the chemical environment of the water. They used powerful computer simulations to calculate the energy of every possible version of a vacancy, from a bare, empty spot to one filled with various combinations of hydrogen and oxygen atoms. By finding the version with the lowest energy for every specific condition, they drew a stability map. This map acts like a weather chart for the atomic world, showing which form of the vacancy will exist at any given pH level (how acidic or basic the water is) and any given electrical voltage.

The results reveal a surprising story about how these materials react to water. In silicon, a vacancy never remains bare when it is in contact with water; it is energetically impossible for the empty spot to stay empty. Instead, it immediately begins to collect atoms from the water. The researcher found that hydrogen atoms bind to the vacancy in a cooperative manner, meaning that once the first hydrogen attaches, it becomes easier for the next ones to follow. This leads to a rapid transformation where the vacancy goes straight from being empty to being fully covered by four hydrogen atoms, effectively becoming a quiet, inert structure. However, this full coverage only happens if the water is chemically reducing, which usually requires a specific electrical push. In more neutral water, the vacancy prefers to grab an oxygen atom instead, forming a structure known as the A center. This oxygen-filled vacancy is stable and electrically active, sitting in the middle of the material's energy gap. Only when the water is made more acidic or the electrical voltage is shifted in a specific direction does the vacancy switch to grabbing a hydroxyl group, or finally, to accepting enough hydrogen to become fully passivated.

The story changes slightly when looking at silicon carbide, a material prized for its durability in high-power electronics. This material has two different types of missing atoms: one where a silicon atom is gone and another where a carbon atom is missing. The researcher discovered that hydrogen behaves very differently at these two sites. At the silicon vacancy, hydrogen binds tightly, forming a strong bond similar to what happens in pure silicon. But at the carbon vacancy, which is responsible for a notorious defect that limits the performance of silicon carbide devices, hydrogen refuses to bind strongly. Instead of forming a tight bond, the hydrogen atom sits loosely between two silicon neighbors, barely holding on. In fact, under normal conditions, the hydrogen is not bound at all and would drift away. This suggests that simply exposing silicon carbide to water or hydrogen gas will not automatically fix the carbon vacancy defect. The defect remains active and capable of disrupting the material's performance unless other factors, such as multiple hydrogen atoms or oxygen, intervene in ways that were not fully captured in this specific simulation.

The study also highlights the importance of the electrical environment. In the simulations, the researcher found that the electrical state of the material, often described by the Fermi level, has very little effect on which vacancy form is stable in silicon. This is because the energy levels that would normally change with voltage are pushed so high in the material's energy gap that they do not influence the chemical reaction with water. The chemical environment—specifically the pH and the presence of water—dominates the outcome. This finding simplifies the picture for engineers: in wet conditions, the chemical nature of the water is the primary driver of how these defects behave, rather than the electrical tuning of the device.

These findings provide a clear, thermodynamic picture of what happens when solid materials meet water. The simulations show that the bare, electrically active vacancy is a transient state that disappears almost instantly in an aqueous environment. It is replaced by decorated complexes that are chemically stable. For silicon, this means the material naturally tends toward an oxygen-filled state in neutral water, while full protection against electrical activity requires specific, more extreme conditions. For silicon carbide, the picture is more complex, with the carbon vacancy resisting simple hydrogen passivation. The researcher notes that while their calculations provide a solid thermodynamic baseline, the actual speed at which these changes happen in the real world depends on how fast atoms can move, a factor that requires further study. Nevertheless, this work offers a fundamental map for predicting how these critical materials will age and perform when exposed to the wet, electrically active environments they often encounter.

Technical Summary: Charge State and Hydration of Vacancies in Silicon and Silicon Carbide in Contact with Water

Problem Statement
The thermodynamic state of point defects in semiconductors is traditionally described solely by the Fermi level ( $E_F$ ). However, in wet-etched, electrochemically processed, or humid environments, vacancies in materials like silicon (Si) and silicon carbide (SiC) exchange electrons with the host while simultaneously adsorbing hydrogen and oxygen from water. The supply of these species is governed by pH and electrode potential. Existing treatments often separate electronic and chemical reservoirs, failing to capture the coupled evolution of a vacancy from a bare, electrically active charge state to a decorated, passivated complex as the aqueous environment changes. This work addresses the need for a unified framework to determine the stable charge state and hydrogen/oxygen content of vacancies under these conditions.

Methodology
The author develops a grand-canonical first-principles model to minimize the grand potential ( $\Omega$ ) of a vacancy complex over its charge state ( $q$ ) and hydrogen/oxygen decoration ( $m_H, m_O$ ).

Grand Potential Formulation: The formation energy is expressed as a grand potential relative to electron, host, hydrogen, and oxygen reservoirs. The chemical potentials of hydrogen ( $\mu_H$ ) and oxygen ( $\mu_O$ ) are tied to a water reservoir via the Computational Hydrogen Electrode (CHE) formalism, linking them to pH and electrode potential ( $U$ ).
Computational Setup: Calculations utilize spin-polarized Density Functional Theory (DFT) with the PBE exchange-correlation functional in 64-atom supercells (2 $\times$ 2 $\times$ 2 repetitions of the conventional cell).
Corrections and Refinements:
- Finite-Size Corrections: Makov–Payne monopole corrections and potential alignment terms are applied to charged supercells.
- Vibrational Free Energy: Harmonic vibrational frequencies are calculated for bound H and O atoms to include zero-point energy and thermal free energy corrections.
- Reference States: Hydrogen and oxygen chemical potentials are referenced to $H_2$ and $H_2O$ molecules, avoiding the overbinding issues of $O_2$ in semilocal functionals.
Materials Studied: The model is parameterized for the silicon vacancy ( $V_{Si}$ ) and its H, O, and OH complexes, as well as the silicon ( $V_{Si}$ ) and carbon ( $V_C$ ) vacancies in 3C-SiC.

Key Results

Silicon Vacancy ( $V_{Si}$ ) in Water:
- Cooperative Hydrogenation: Hydrogen binds cooperatively to the silicon vacancy. The formation energy decreases monotonically with each added hydrogen, leading directly from the bare vacancy to the fully passivated $VH_4$ complex. Intermediate hydrogenated species ( $VH$ , $VH_2$ , $VH_3$ ) are thermodynamically unstable against disproportionation.
- Stability Map: In contact with water, the bare vacancy is never the lowest-energy state.
  - Near neutral pH, the vacancy binds oxygen to form the closed-shell $V-O$ (A center).
  - The paramagnetic $V-OH$ complex forms in acidic conditions or under mild cathodic bias.
  - Full passivation to $VH_4$ requires a cathodic bias (near the hydrogen evolution potential).
- Fermi Level Dependence: Within the PBE functional, the stability boundaries are largely independent of the Fermi level because the acceptor levels of the vacancy and its hydrides lie at or near the conduction-band edge (which is underestimated by PBE).
Silicon Carbide (3C-SiC) Vacancies:
- Asymmetric Hydrogen Binding: The two sublattice vacancies exhibit fundamentally different hydrogen binding behaviors.
  - Silicon Vacancy ( $V_{Si}$ ): The first hydrogen binds strongly (1.97 eV) via a standard C–H bond.
  - Carbon Vacancy ( $V_C$ , the Z1/2 center): The first hydrogen does not form a terminal Si–H bond but instead forms a weak three-center Si–H–Si bridge. This complex is unbound with respect to $H_2$ gas ( $E_{bind} \approx -0.13$ eV) and requires a significant cathodic bias to form.
- Implications for Z1/2: Since the Z1/2 center is the carbon vacancy, a single hydrogen atom does not effectively deactivate it. The stability of multiply hydrogenated complexes or oxygen/hydroxyl species at $V_C$ remains an open question within this specific study's scope.
Validation and Limitations:
- The DFT results reproduce experimental lattice constants, the bridging geometry of the A center, and Si–H stretch frequencies (within ~6.5% without scaling) for $VH_1$ – $VH_4$ .
- The study acknowledges limitations: the PBE functional underestimates band gaps, leading to semiquantitative absolute transition levels and a failure to recover the negative- $U$ ordering of the silicon vacancy levels observed experimentally. The 64-atom cell size also limits the convergence of formation energies compared to larger cells.

Significance and Claims
The paper claims to provide the first stability map that simultaneously tracks a single bulk point defect's charge state and chemical decoration (H/O content) as a function of both Fermi level and aqueous environment (pH/potential).

Unified Framework: It successfully integrates charged-defect thermodynamics, the computational hydrogen electrode, and aqueous stability analysis into a single grand-canonical framework.
Thermodynamic Ground States: It establishes that in aqueous environments, the "bare" electrically active vacancy is not the thermodynamic ground state; instead, oxygen or hydrogen passivation dominates depending on the electrochemical conditions.
Material Specificity: The work highlights a critical asymmetry in SiC, where the carbon vacancy (Z1/2) resists simple hydrogen passivation compared to the silicon vacancy, suggesting that the deactivation of deep levels in SiC power devices via hydrogenation may require more complex mechanisms (e.g., multiple hydrogens or oxygen incorporation) than previously assumed for single defects.
Applicability: The author suggests this approach is applicable to other materials in contact with water (oxides, nitrides, perovskites, 2D materials) where charge state and hydration jointly control corrosion, stability, or dopant passivation.

The study concludes that while the current PBE-based model provides a robust thermodynamic baseline, resolving the precise coupling between the Fermi level and hydration (particularly for negative- $U$ centers) will require hybrid functionals. Experimental validation, such as comparing DLTS signals before and after controlled aqueous treatments, is proposed as a direct method to settle questions regarding the deactivation of the Z1/2 center.

More like this