How does quenched disorder destroy crystalline coherence in two-dimensional electron solids?

Disorder shortens crystalline coherence rather than switching it off

The paper asks how far crystalline organization survives once randomly positioned impurities are present, and how that distance shrinks as disorder grows. Coulomb repulsion favors a triangular arrangement, while impurities favor an irregular density pattern, so the two compete over the same electrons. The answer reported here is a crossover, not a boundary: increasing disorder weakens crystalline correlations, broadens the structure factor, and reduces the measured coherence length toward the pristine lattice spacing.

That framing matters because it rules out the intuitive picture of a finite-disorder line separating a perfect crystal from an amorphous solid. The authors connect their numerical evolution to the Imry-Ma and Larkin picture, in which arbitrarily weak random pinning destabilizes pristine translational order on a sufficiently long scale, and stronger disorder shortens the distortion scale further. Weakly pinned domains can therefore still look crystalline over a finite observation window, which is why the paper treats the diagrams as degrees of noncrystallinity rather than as a phase boundary.

Two microscopic models test the same structural question

The paper runs two parallel calculations with separate Hamiltonians, clean reference states, disorder control parameters, and numerical checks. The quantum study uses continuum Hartree-Fock theory to follow the self-consistent electron density in a short-range random potential, with a supercell of nine electrons spanning three primitive lattice vectors in each direction of a triangular crystal. A direct-only Hartree calculation with exactly the same disorder and discretization acts as a control that isolates the role of exchange, since the clean direct-only Hartree density is uniform while the selected Hartree-Fock states at r_s = 3, 4, 8, and 16 are crystalline.

The classical study uses molecular dynamics and zero-temperature relaxation to follow point electrons with Coulomb interactions and open electrostatic boundaries in randomly positioned attractive wells. It resolves particle configurations, coordination defects, positional and orientational correlations, and cage escape. The two disorder protocols are deliberately complementary: Hartree-Fock increases the amplitude of one fixed Gaussian random potential, while molecular dynamics increases the concentration of attractive wells in one nested Poisson impurity field while holding their physical depth and radius fixed. Neither calculation averages over disorder realizations, so the two disorder axes are defined separately and the comparison concerns a shared physical mechanism rather than a numerical conversion between them.

Quantum density weight leaves the pristine Bragg vectors

In the quantum calculation, the elastic density intensity at a nonzero supercell reciprocal vector measures the static spatial organization of the self-consistent density. It is the coherent contribution built from the expectation value of the density operator with a per-cell normalization, not the full equal-time quantum structure factor, which would also contain a connected exchange and fluctuation term that this work does not evaluate. With that definition in hand, increasing disorder transfers weight from the pristine sixfold shell to other supercell reciprocal vectors, and periodic boundary conditions keep these spectral lines discrete with no artificial peak broadening applied.

The paper also reports an operational intra-cell coherence length. At r_s = 4 it changes from an unresolved value above 1.5a at weak disorder to approximately 1.1a and 0.9a at W/E_F = 2 and 4. The supercell spans only three pristine lattice spacings along each primitive direction, so the density-coherence diagnostic is evaluated within half that range and is explicitly distinguished from an asymptotic correlation length. The quantum map describes the computed density within one repeated cell and does not establish the thermodynamic decay of correlations in independently disordered cells.

The classical structure factor broadens as pin density rises

In the classical calculation at electron density n_e = n_0, the structure factor evolves from six sharp peaks into diffuse scattering as impurity density increases, and both positional and orientational correlations weaken. The measured translational coherence decreases from more than 5a at weak pinning to approximately 2.7a, 1.2a, and 0.8a at impurity concentrations c = 1/8, 1/4, and 4. The strong-disorder endpoint is a mechanically stable amorphous configuration with only local packing order, and the bandwidth checks establish an order-of-scale result rather than a precise distance below one lattice constant.

A two-density map at fixed well radius and depth extends this result across electron density, with fifty computed points and crystal-like states defined by satisfying both an orientational order threshold and a sixfold fraction threshold. The authors stress that these contours are finite-sample structural crossovers. They also note that the Bragg fraction and the first-shell coherence measure are orientation selective, so a rotated or reconstructed crystal can lose weight at the fixed reference vectors without having lost all order.

Stronger pinning raises the temperature for rearrangement

Heating at fixed impurity configurations separates two effects that are easy to conflate. The primary cage-escape bracket rises from k_B T/E_0 = 0.005-0.007 in the clean system to 0.020-0.030 at c = 4, where crystalline order has already disappeared at zero temperature. Random wells disrupt positional and orientational coherence, reducing the stability range of a recognizable crystal, but they also create barriers to particle rearrangement, allowing an already irregular configuration to remain caged to higher temperatures.

The paper is careful about what this mobility crossover is. It is a finite-observation-time crossover at lags 100 and 400 t_0, not a diffusion coefficient and not a thermodynamic glass transition; the authors state that there is no glass in the presence of explicit external disorder in the Hamiltonian, since the Hamiltonian lacks translational invariance. The reported bars span sampled crossing intervals and are not statistical confidence intervals. Sampling checks matter here: a separate clean N = 572 calculation gives different apparent thresholds than the primary N = 288 protocol, so the defensible clean scale is roughly theta = 0.005-0.010 rather than a high-precision number.

Larkin estimates track the measured coherence scale

The paper derives weak-pinning Larkin estimates for both models using clean elastic stiffness and the known disorder covariance, with no numerical disorder averaging. For the Hartree-Fock amplitude scan at fixed r_s, the weak-pinning scaling is R_c proportional to W^-1, with estimates of about 3.2a, 2.5a, 1.9a, and 2.3a at r_s = 3, 4, 8, and 16. At r_s = 4 and eta = 0.5 the estimate is about 5a, already larger than the simulated quantum cell, and at eta = 4 the formal value is about 0.6a, which the authors read as breakdown of a collective crystal description rather than a resolved sub-lattice elastic domain.

For the classical concentration scan, the pristine triangular Coulomb lattice has shear modulus mu_cl = 0.245065 (e^2/4 pi epsilon) n_e^{3/2}, and the compact wells at q = 1 satisfy 2R_p < a so independent Poisson impurities give uncorrelated forces at distinct pristine lattice sites. Substituting the actual well depth and radius gives R_c/a approximately 0.63/sqrt(c) and a linear Gaussian estimate of about 0.87/sqrt(c). At c = 1/8 the estimates are R_c about 1.8a and the Gaussian value about 2.5a, compared with the measured coherence of about 2.7a. The authors call this order-of-magnitude agreement without fitted parameters, not a validation of the prefactor, and note that a Larkin disk holds only about 1.5 wells at every c, so R_c simply tracks the mean well spacing. The two scans have different power laws because their random-force variances scale as W^2 and as pin density times well depth squared, respectively.

Large density modulation is not evidence of crystal order

One of the paper's sharpest points concerns a diagnostic that is easy to misuse. In the quantum study, the density contrast grows while crystalline organization is lost, so large modulation cannot serve as a proxy for crystal order. Both a crystal and an amorphous pinned assembly can have strongly modulated density, which is why the authors argue that spatial correlations and their Fourier spectra test crystalline organization while mobility measurements test rearrangement.

This distinction also drives the paper's closing argument against a specific claim in the literature. The authors conclude that saying extrinsic disorder enhances Wigner crystallization and lowers the critical r_s is misleading: disorder may enhance some short-range order, but the long-range order characteristic of a crystal is always suppressed, and the system becomes an amorphous solid. Their stated reason is that if increasing disorder enhanced crystallization, one could arbitrarily lower the critical r_s by adding impurities, which makes little sense because high impurity density must always lead to an amorphous structure. The paper's own evidence for this is the combination of a shrinking coherence length in both models and a rising rearrangement temperature in the classical model, which together describe a pinned amorphous solid rather than a better crystal.

What these results do and do not establish

The quantum calculation uses a mean-field Hartree-Fock solid, and the authors acknowledge that correlated calculations place clean crystallization at much lower densities than Hartree-Fock, with reported critical r_s values of 37 plus or minus 5, 31 with an uncertainty of 1, and 35 with an uncertainty of 1 from different Monte Carlo studies. They argue that the mean-field technique keeps the physics transparent even though the actual r_s regime of interest shifts, which they consider acceptable for a qualitative understanding of the interaction-disorder competition. The classical calculation is restricted to q = 1; at sufficiently high electron density, neighboring-site force correlations must be retained because 2R_p < a no longer holds.

Several reported quantities are explicitly operational rather than thermodynamic. The intra-cell coherence diagnostic is evaluated within half the supercell range, the classical structural contours are finite-sample crossovers, the mobility crossover depends on observation time, and the strong-disorder coherence scale is an order-of-scale result. The paper also notes that a reentrant ordered interval at weak disorder was predicted by Nelson and found in earlier work, but these trajectories do not resolve a reentrant ordered interval or a precise melting temperature. The most useful takeaway is therefore a comparison rather than a threshold: crystalline coherence and resistance to motion are distinct properties, and in this paper disorder suppresses the first while strengthening the second, so a pinned two-dimensional electron solid can look less like a crystal and still be harder to rearrange.

Sources used in this answer