Two-Dimensional Wigner Crystal in Anisotropic Semiconductor

Xin Wan,R. N. Bhatt
DOI: https://doi.org/10.1103/PhysRevB.65.233209
1999-12-22
Abstract:We investigate the effect of mass anisotropy on the Wigner crystallization transition in a two-dimensional (2D) electron gas. The static and dynamical properties of a 2D Wigner crystal have been calculated for arbitrary 2D Bravais lattices in the presence of anisotropic mass, as may be obtainable in Si MOSFETs with (110) surface. By studying the stability of all possible lattices, we find significant change in the crystal structure and melting density of the electron lattice with the lowest ground state energy.
Mesoscale and Nanoscale Physics,Strongly Correlated Electrons
What problem does this paper attempt to address?
This paper aims to explore the influence of mass anisotropy on the Wigner crystallization transition in two - dimensional (2D) electron gases. Specifically, the authors studied the static and dynamic properties of different lattice structures in a two - dimensional electron system with anisotropic mass, and evaluated how these properties affect the stability of the crystal structure and its melting density. ### Research Background 1. **Metal - Insulator Transition (MIT)**: In the two - dimensional electron systems in semiconductor structures, the possible existence of metal - insulator transition (MIT) has been experimentally observed. However, according to the localization theory, MIT in two - dimensional systems under pure scattering potential does not exist. Therefore, strong Coulomb interaction may be the cause of this phenomenon. 2. **Wigner Crystallization**: In the low - density limit, two - dimensional electrons will form a triangular lattice (Wigner crystallization) due to electron - electron interactions. Ideally, the critical parameter \( r_s \) is predicted to be 37 ± 5. However, the presence of impurities may reduce this critical value. ### Research Questions The key questions that the paper attempts to solve are: - In a two - dimensional electron system with mass anisotropy, how will the critical density and crystal structure of Wigner crystallization change? - What is the influence of mass anisotropy on the ground - state energy of two - dimensional Wigner crystals? ### Main Findings 1. **Ground - State Energy Calculation**: By calculating the ground - state energy of any two - dimensional Bravais lattice, the authors found that mass anisotropy significantly changes the crystal structure and melting density. In particular, for the (110) surface in Si MOSFETs, mass anisotropy may lead to a significant reduction in the critical density. 2. **Stability and Lattice Selection**: In the presence of mass anisotropy, the most stable lattice structure is no longer the triangular lattice, but a quasi - one - dimensional lattice stretched along certain specific directions. For example, at \( r_s = 60 \), no stable lattice can be formed. 3. **Critical Density**: For an ideal clean two - dimensional anisotropic electron system (such as the silicon (110) surface), the critical density \( r_s \) of Wigner crystallization may be much lower than the value predicted for a homogeneous system, and is estimated to be 80. ### Conclusion The paper shows that mass anisotropy has an important influence on the Wigner crystallization of two - dimensional electron systems. Especially in the low - density limit, it may lead to a significant reduction in the critical density. This provides a new perspective for further understanding the metal - insulator transition in two - dimensional electron systems. ### Formula Summary - **Ground - State Energy Formula**: \[ E_s = -\frac{1}{N} \sum_{i < j} \frac{e^2}{|r_i - r_j|} = -\frac{2e^2}{\sqrt{v}} \left[ 2 - \sum_{r \neq 0} \phi^{-1/2}\left( \frac{\pi |r|^2}{v} \right) \right] \] where \( \phi_n(x) \) is the Misra function, defined as: \[ \phi_n(x) = \int_0^\infty dt \, t^n e^{-xt} \] - **Hamiltonian**: \[ H = \sum_i \sum_\alpha \frac{p_\alpha(r_i)^2}{2m_\alpha} + \frac{1}{2} \sum_{i,j} \sum_{\alpha,\beta} u_\alpha(r_i) \Phi_{\alpha\beta}(r_i, r_j) u_\beta(r_j) \] - **Dynamics Matrix**: \[ D(q) = \sum_{R_0} \Phi(R_0) e^{-iq \cdot R_0} \] - **Ground - State Energy Expression**: \[ E_0 = -\frac{2.21}{r_s} + \frac{1.63}{r_s^{3/2}} \] These