Phase diagram of the low-density two-dimensional homogeneous electron gas

Evidence for electron-hole crystals in a Mott insulator

The coexistence of correlated electron and hole crystals enables the realization of quantum excitonic states, capable of hosting counterflow superfluidity and topological orders with long-range quantum entanglement. Here we report evidence for imbalanced electron-hole crystals in a doped Mott insulator, namely, α-RuCl3, through gate-tunable non-invasive van der Waals doping from graphene. Real-space imaging via scanning tunnelling microscopy reveals two distinct charge orderings at the lower and upper Hubbard band energies, whose origin is attributed to the correlation-driven honeycomb hole crystal composed of hole-rich Ru sites and rotational-symmetry-breaking paired electron crystal composed of electron-rich Ru-Ru bonds, respectively. Moreover, a gate-induced transition of electron-hole crystals is directly visualized, further corroborating their nature as correlation-driven charge crystals. The realization and atom-resolved visualization of imbalanced electron-hole crystals in a doped Mott insulator opens new doors in the search for correlated bosonic states within strongly correlated materials.

Charge-transfer contacts for the measurement of correlated states in high-mobility WSe2

Two-dimensional semiconductors, such as transition metal dichalcogenides, have demonstrated tremendous promise for the development of highly tunable quantum devices. Realizing this potential requires low-resistance electrical contacts that perform well at low temperatures and low densities where quantum properties are relevant. Here we present a new device architecture for two-dimensional semiconductors that utilizes a charge-transfer layer to achieve large hole doping in the contact region, and implement this technique to measure the magnetotransport properties of high-purity monolayer WSe2. We measure a record-high hole mobility of 80,000 cm2 V-1 s-1 and access channel carrier densities as low as 1.6 × 1011 cm-2, an order of magnitude lower than previously achievable. Our ability to realize transparent contact to high-mobility devices at low density enables transport measurements of correlation-driven quantum phases including the observation of a low-temperature metal-insulator transition in a density and temperature regime where Wigner crystal formation is expected and the observation of the fractional quantum Hall effect under large magnetic fields. The charge-transfer contact scheme enables the discovery and manipulation of new quantum phenomena in two-dimensional semiconductors and their heterostructures.

Imaging tunable Luttinger liquid systems in van der Waals heterostructures

One-dimensional (1D) interacting electrons are often described as a Luttinger liquid1-4 having properties that are intrinsically different from those of Fermi liquids in higher dimensions5,6. In materials systems, 1D electrons exhibit exotic quantum phenomena that can be tuned by both intra- and inter-1D-chain electronic interactions, but their experimental characterization can be challenging. Here we demonstrate that layer-stacking domain walls (DWs) in van der Waals heterostructures form a broadly tunable Luttinger liquid system, including both isolated and coupled arrays. We have imaged the evolution of DW Luttinger liquids under different interaction regimes tuned by electron density using scanning tunnelling microscopy. Single DWs at low carrier density are highly susceptible to Wigner crystallization consistent with a spin-incoherent Luttinger liquid, whereas at intermediate densities dimerized Wigner crystals form because of an enhanced magneto-elastic coupling. Periodic arrays of DWs exhibit an interplay between intra- and inter-chain interactions that gives rise to new quantum phases. At low electron densities, inter-chain interactions are dominant and induce a 2D electron crystal composed of phased-locked 1D Wigner crystal in a staggered configuration. Increased electron density causes intra-chain fluctuation potentials to dominate, leading to an electronic smectic liquid crystal phase in which electrons are ordered with algebraical correlation decay along the chain direction but disordered between chains. Our work shows that layer-stacking DWs in 2D heterostructures provides opportunities to explore Luttinger liquid physics.

Nematic Metal in a Multivalley Electron Gas: Variational Monte Carlo Analysis and Application to AlAs

The two-dimensional electron gas is of fundamental importance in quantum many-body physics. We study a minimal extension of this model with C_{4} (as opposed to full rotational) symmetry and an electronic dispersion with two valleys with anisotropic effective masses. Electrons in our model interact via Coulomb repulsion, screened by distant metallic gates. Using variational Monte Carlo simulations, we find a broad intermediate range of densities with a metallic valley-polarized, spin-unpolarized ground state. Our results are of direct relevance to the recently discovered "nematic" state in AlAs quantum wells. For the effective mass anisotropy relevant to this system, m_{x}/m_{y}≈5.2, we obtain a transition from an anisotropic metal to a valley-polarized metal at r_{s}≈12 (where r_{s} is the dimensionless Wigner-Seitz radius). At still lower densities, we find a (possibly metastable) valley and spin-polarized state with a reduced electronic anisotropy.

Sublattice Structure and Topology in Spontaneously Crystallized Electronic States

The prediction and realization of the quantum anomalous Hall effect are often intimately connected to honeycomb lattices in which the sublattice degree of freedom plays a central role in the nontrivial topology. Two-dimensional Wigner crystals, on the other hand, form triangular lattices without sublattice degrees of freedom, resulting in a topologically trivial state. Here, we discuss the possibility of spontaneously formed honeycomb-lattice crystals that exhibit the quantum anomalous Hall effect. Starting from a single-band system with nontrivial quantum geometry, we derive the mean-field energy functional of a class of crystal states and express it as a model of sublattice pseudospins in momentum space. We find that nontrivial quantum geometry leads to extra terms in the pseudospin model that break an effective "time-reversal symmetry" and favor a topologically nontrivial pseudospin texture. When the effects of these extra terms dominate over the ferromagnetic exchange coupling between pseudospins, the anomalous Hall crystal state becomes energetically favorable over the trivial Wigner crystal state.

Non-equilibrium quantum domain reconfiguration dynamics in a two-dimensional electronic crystal and a quantum annealer

Relaxation dynamics of complex many-body quantum systems trapped into metastable states is a very active field of research from both the theoretical and experimental point of view with implications in a wide array of topics from macroscopic quantum tunnelling and nucleosynthesis to non-equilibrium superconductivity and energy-efficient memory devices. In this work, we investigate quantum domain reconfiguration dynamics in the electronic superlattice of a quantum material using time-resolved scanning tunneling microscopy and unveil a crossover from temperature to noisy quantum fluctuation dominated dynamics. The process is modeled using a programmable superconducting quantum annealer in which qubit interconnections correspond directly to the microscopic interactions between electrons in the quantum material. Crucially, the dynamics of both the experiment and quantum simulation is driven by spectrally similar pink noise. We find that the simulations reproduce the emergent time evolution and temperature dependence of the experimentally observed electronic domain dynamics.

The Emergence of the Hexagonal Lattice in Two-Dimensional Wigner Fragments

At very low density, the electrons in a uniform electron gas spontaneously break symmetry and form a crystalline lattice called a Wigner crystal. But which type of crystal will the electrons form? We report a numerical study of the density profiles of fragments of Wigner crystals from first principles. To simulate Wigner fragments, we use Clifford periodic boundary conditions and a renormalized distance in the Coulomb potential. Moreover, we show that high-spin restricted open-shell Hartree-Fock theory becomes exact in the low-density limit. We are thus able to accurately capture the localization in two-dimensional Wigner fragments with many electrons. No assumptions about the positions where the electrons will localize are made. The density profiles we obtain emerge naturally when we minimize the total energy of the system. We clearly observe the emergence of the hexagonal crystal structure, which has been predicted to be the ground-state structure of the two-dimensional Wigner crystal.

Signatures of Correlated Defects in an Ultraclean Wigner Crystal in the Extreme Quantum Limit

Low-disorder two-dimensional electron systems in the presence of a strong, perpendicular magnetic field terminate at very small Landau level filling factors in a Wigner crystal (WC), where the electrons form an ordered array to minimize the Coulomb repulsion. The nature of this exotic, many-body, quantum phase is yet to be fully understood and experimentally revealed. Here we probe one of WC's most fundamental parameters, namely, the energy gap that determines its low-temperature conductivity, in record mobility, ultrahigh-purity, two-dimensional electrons confined to GaAs quantum wells. The WC domains in these samples contain ≃1000 electrons. The measured gaps are a factor of three larger than previously reported for lower quality samples, and agree remarkably well with values predicted for the lowest-energy, intrinsic, hypercorrelated bubble defects in a WC made of flux-electron composite fermions, rather than bare electrons. The agreement is particularly noteworthy, given that the calculations are done for disorder-free composite fermion WCs, and there are no adjustable parameters. The results reflect the exceptionally high quality of the samples, and suggest that composite fermion WCs are indeed more stable compared to their electron counterparts.

Stoner Ferromagnetism in a Momentum-Confined Interacting 2D Electron Gas

In this work we investigate the ground state of a momentum-confined interacting 2D electron gas, a momentum-space analog of an infinite quantum well. The study is performed by combining analytical results with a numerical exact diagonalization procedure. We find a ferromagnetic ground state near a particular electron density and for a range of effective electron (or hole) masses. We argue that this observation may be relevant to the generalized Stoner ferromagnetism recently observed in multilayer graphene systems. The collective magnon excitations exhibit a linear dispersion, which originates from a diverging spin stiffness.

Metal-Insulator Transition in a Semiconductor Heterobilayer Model

Transition metal dichalcogenide superlattices provide an exciting new platform for exploring and understanding a variety of phases of matter. The moiré continuum Hamiltonian, of two-dimensional jellium in a modulating potential, provides a fundamental model for such systems. Accurate computations with this model are essential for interpreting experimental observations and making predictions for future explorations. In this work, we combine two complementary quantum Monte Carlo (QMC) methods, phaseless auxiliary field quantum Monte Carlo and fixed-phase diffusion Monte Carlo, to study the ground state of this Hamiltonian. We observe a metal-insulator transition between a paramagnet and a 120° Néel ordered state as the moiré potential depth and the interaction strength are varied. We find significant differences from existing results by Hartree-Fock and exact diagonalization studies. In addition, we benchmark density-functional theory, and suggest an optimal hybrid functional which best approximates our QMC results.

Crystallization of polarons through charge and spin ordering transitions in 1T-TaS2

The interaction of electrons with the lattice in metals can lead to reduction of their kinetic energy to the point where they may form heavy, dressed quasiparticles-polarons. Unfortunately, polaronic lattice distortions are difficult to distinguish from more conventional charge- and spin-ordering phenomena at low temperatures. Here we present a study of local symmetry breaking of the lattice structure on the picosecond timescale in the prototype layered dichalcogenide Mott insulator 1T-TaS2 using X-ray pair-distribution function measurements. We clearly identify symmetry-breaking polaronic lattice distortions at temperatures well above the ordered phases, and record the evolution of broken symmetry states from 915 K to 15 K. The data imply that charge ordering is driven by polaron crystallization into a Wigner crystal-like state, rather than Fermi surface nesting or conventional electron-phonon coupling. At intermediate temperatures the local lattice distortions are found to be consistent with a quantum spin liquid state.

Synergistic correlated states and nontrivial topology in coupled graphene-insulator heterostructures

Graphene has aroused great attention due to the intriguing properties associated with its low-energy Dirac Hamiltonian. When graphene is coupled with a correlated insulating substrate, electronic states that cannot be revealed in either individual layer may emerge in a synergistic manner. Here, we theoretically study the correlated and topological states in Coulomb-coupled and gate-tunable graphene-insulator heterostructures. By electrostatically aligning the electronic bands, charge carriers transferred between graphene and the insulator can yield a long-wavelength electronic crystal at the interface, exerting a superlattice Coulomb potential on graphene and generating topologically nontrivial subbands. This coupling can further boost electron-electron interaction effects in graphene, leading to a spontaneous bandgap formation at the Dirac point and interaction-enhanced Fermi velocity. Reciprocally, the electronic crystal at the interface is substantially stabilized with the help of cooperative interlayer Coulomb coupling. We propose a number of substrate candidates for graphene to experimentally demonstrate these effects.

Interstitial-Induced Ferromagnetism in a Two-Dimensional Wigner Crystal

The two-dimensional Wigner crystal (WC) occurs in the strongly interacting regime (r_{s}≫1) of the two-dimensional electron gas (2DEG). The magnetism of a pure WC is determined by tunneling processes that induce multispin ring-exchange interactions, resulting in fully polarized ferromagnetism for large enough r_{s}. Recently, Hossain et al. [Proc. Natl. Acad. Sci. U.S.A. 117, 32244 (2020)PNASA60027-842410.1073/pnas.2018248117] reported the occurrence of a fully polarized ferromagnetic insulator at r_{s}≳35 in an AlAs quantum well, but at temperatures orders of magnitude larger than the predicted exchange energies for the pure WC. Here, we analyze the large r_{s} dynamics of an interstitial defect in the WC, and show that it produces local ferromagnetism with much higher energy scales. Three hopping processes are dominant, which favor a large, fully polarized ferromagnetic polaron. Based on the above results, we speculate concerning the phenomenology of the magnetism near the metal-insulator transition of the 2DEG.

Tunable quantum criticalities in an isospin extended Hubbard model simulator

Studying strong electron correlations has been an essential driving force for pushing the frontiers of condensed matter physics. In particular, in the vicinity of correlation-driven quantum phase transitions (QPTs), quantum critical fluctuations of multiple degrees of freedom facilitate exotic many-body states and quantum critical behaviours beyond Landau's framework¹. Recently, moiré heterostructures of van der Waals materials have been demonstrated as highly tunable quantum platforms for exploring fascinating, strongly correlated quantum physics^2-22. Here we report the observation of tunable quantum criticalities in an experimental simulator of the extended Hubbard model with spin-valley isospins arising in chiral-stacked twisted double bilayer graphene (cTDBG). Scaling analysis shows a quantum two-stage criticality manifesting two distinct quantum critical points as the generalized Wigner crystal transits to a Fermi liquid by varying the displacement field, suggesting the emergence of a critical intermediate phase. The quantum two-stage criticality evolves into a quantum pseudo criticality as a high parallel magnetic field is applied. In such a pseudo criticality, we find that the quantum critical scaling is only valid above a critical temperature, indicating a weak first-order QPT therein. Our results demonstrate a highly tunable solid-state simulator with intricate interplay of multiple degrees of freedom for exploring exotic quantum critical states and behaviours.

Anisotropic Two-Dimensional Disordered Wigner Solid

The interplay between the Fermi sea anisotropy, electron-electron interaction, and localization phenomena can give rise to exotic many-body phases. An exciting example is an anisotropic two-dimensional (2D) Wigner solid (WS), where electrons form an ordered array with an anisotropic lattice structure. Such a state has eluded experiments up to now as its realization is extremely demanding: First, a WS entails very low densities where the Coulomb interaction dominates over the kinetic (Fermi) energy. Attaining such low densities while keeping the disorder low is very challenging. Second, the low-density requirement has to be fulfilled in a material that hosts an anisotropic Fermi sea. Here, we report transport measurements in a clean (low-disorder) 2D electron system with anisotropic effective mass and Fermi sea. The data reveal that at extremely low electron densities, when the r_{s} parameter, the ratio of the Coulomb to the Fermi energy, exceeds ≃38, the current-voltage characteristics become strongly nonlinear at small dc biases. Several key features of the nonlinear characteristics, including their anisotropic voltage thresholds, are consistent with the formation of a disordered, anisotropic WS pinned by the ubiquitous disorder potential.

The Magnetic Genome of Two-Dimensional van der Waals Materials

Magnetism in two-dimensional (2D) van der Waals (vdW) materials has recently emerged as one of the most promising areas in condensed matter research, with many exciting emerging properties and significant potential for applications ranging from topological magnonics to low-power spintronics, quantum computing, and optical communications. In the brief time after their discovery, 2D magnets have blossomed into a rich area for investigation, where fundamental concepts in magnetism are challenged by the behavior of spins that can develop at the single layer limit. However, much effort is still needed in multiple fronts before 2D magnets can be routinely used for practical implementations. In this comprehensive review, prominent authors with expertise in complementary fields of 2D magnetism (i.e., synthesis, device engineering, magneto-optics, imaging, transport, mechanics, spin excitations, and theory and simulations) have joined together to provide a genome of current knowledge and a guideline for future developments in 2D magnetic materials research.

Competing correlated states around the zero-field Wigner crystallization transition of electrons in two dimensions

The competition between kinetic energy and Coulomb interactions in electronic systems leads to complex many-body ground states with competing orders. Here we present zinc oxide-based two-dimensional electron systems as a high-mobility system to study the low-temperature phases of strongly interacting electrons. An analysis of the electronic transport provides evidence for competing correlated metallic and insulating states with varying degrees of spin polarization. Some features bear quantitative resemblance to quantum Monte Carlo simulation results, including the transition point from the paramagnetic Fermi liquid to Wigner crystal and the absence of a Stoner transition. At very low temperatures, we resolve a non-monotonic spin polarizability of electrons across the phase transition, pointing towards a low spin phase of electrons, and a two-order-of-magnitude positive magnetoresistance that is challenging to understand within traditional metallic transport paradigms. This work establishes zinc oxide as a platform for studying strongly correlated electrons in two dimensions.

Terahertz Fingerprint of Monolayer Wigner Crystals

The strong Coulomb interaction in monolayer semiconductors represents a unique opportunity for the realization of Wigner crystals without external magnetic fields. In this work, we predict that the formation of monolayer Wigner crystals can be detected by their terahertz response spectrum, which exhibits a characteristic sequence of internal optical transitions. We apply the density matrix formalism to derive the internal quantum structure and the optical conductivity of the Wigner crystal and to microscopically analyze the multipeak shape of the obtained terahertz spectrum. Moreover, we predict a characteristic shift of the peak position as a function of charge density for different atomically thin materials and show how our results can be generalized to an arbitrary two-dimensional system.

Wigner localization in two and three dimensions: An ab initio approach

In this work, we investigate the Wigner localization of two interacting electrons at very low density in two and three dimensions using the exact diagonalization of the many-body Hamiltonian. We use our recently developed method based on Clifford periodic boundary conditions with a renormalized distance in the Coulomb potential. To accurately represent the electronic wave function, we use a regular distribution in space of Gaussian-type orbitals and we take advantage of the translational symmetry of the system to efficiently calculate the electronic wave function. We are thus able to accurately describe the wave function up to very low density. We validate our approach by comparing our results to a semi-classical model that becomes exact in the low-density limit. With our approach, we are able to observe the Wigner localization without ambiguity.

Interaction-Driven Filling-Induced Metal-Insulator Transitions in 2D Moiré Lattices

Using a realistic band structure for twisted WSe_{2} materials, we develop a theory for the interaction-driven correlated insulators to conducting metals transitions through the tuning of the filling factor around commensurate fractional fillings of the moiré unit cell in the 2D honeycomb lattice, focusing on the dominant half-filled Mott insulating state, which exists for both long- and short-range interactions. We find metallic states slightly away from half-filling, as have recently been observed experimentally. We discuss the stabilities and the magnetic properties of the resulting insulating and metallic phases, and comment on their experimental signatures. We also discuss the nature of the correlated insulator states at the rational fractional fillings.

Wigner Crystallization of Electrons in a One-Dimensional Lattice: A Condensation in the Space of States

We study the ground state of a system of spinless electrons interacting through a screened Coulomb potential in a lattice ring. By using analytical arguments, we show that, when the effective interaction compares with the kinetic energy, the system forms a Wigner crystal undergoing a first-order quantum phase transition. This transition is a condensation in the space of the states and belongs to the class of quantum phase transitions discussed in [M. Ostilli and C. Presilla, J. Phys. A 54, 055005 (2021).JPAMB51751-811310.1088/1751-8121/aba144]. The transition takes place at a critical value r_{s}_{c} of the usual dimensionless parameter r_{s} (radius of the volume available to each electron divided by effective Bohr radius) for which we are able to provide rigorous lower and upper bounds. For large screening length these bounds can be expressed in a closed analytical form. Demanding Monte Carlo simulations allow to estimate r_{s}_{c}≃2.3±0.2 at lattice filling 3/10 and screening length 10 lattice constants. This value is well within the rigorous bounds 0.7≤r_{s}_{c}≤4.3. Finally, we show that if screening is removed after the thermodynamic limit has been taken, r_{s}_{c} tends to zero. In contrast, in a bare unscreened Coulomb potential, Wigner crystallization always takes place as a smooth crossover, not as a quantum phase transition.

Signatures of Wigner crystal of electrons in a monolayer semiconductor

When the Coulomb repulsion between electrons dominates over their kinetic energy, electrons in two-dimensional systems are predicted to spontaneously break continuous-translation symmetry and form a quantum crystal¹. Efforts to observe^2-12 this elusive state of matter, termed a Wigner crystal, in two-dimensional extended systems have primarily focused on conductivity measurements on electrons confined to a single Landau level at high magnetic fields. Here we use optical spectroscopy to demonstrate that electrons in a monolayer semiconductor with density lower than 3 × 10¹¹ per centimetre squared form a Wigner crystal. The combination of a high electron effective mass and reduced dielectric screening enables us to observe electronic charge order even in the absence of a moiré potential or an external magnetic field. The interactions between a resonantly injected exciton and electrons arranged in a periodic lattice modify the exciton bandstructure so that an umklapp resonance arises in the optical reflection spectrum, heralding the presence of charge order¹³. Our findings demonstrate that charge-tunable transition metal dichalcogenide monolayers¹⁴ enable the investigation of previously uncharted territory for many-body physics where interaction energy dominates over kinetic energy.

Observation of spontaneous ferromagnetism in a two-dimensional electron system

What are the ground states of an interacting, low-density electron system? In the absence of disorder, it has long been expected that as the electron density is lowered, the exchange energy gained by aligning the electron spins should exceed the enhancement in the kinetic (Fermi) energy, leading to a (Bloch) ferromagnetic transition. At even lower densities, another transition to a (Wigner) solid, an ordered array of electrons, should occur. Experimental access to these regimes, however, has been limited because of the absence of a material platform that supports an electron system with very high quality (low disorder) and low density simultaneously. Here we explore the ground states of interacting electrons in an exceptionally clean, two-dimensional electron system confined to a modulation-doped AlAs quantum well. The large electron effective mass in this system allows us to reach very large values of the interaction parameter [Formula: see text], defined as the ratio of the Coulomb to Fermi energies. As we lower the electron density via gate bias, we find a sequence of phases, qualitatively consistent with the above scenario: a paramagnetic phase at large densities, a spontaneous transition to a ferromagnetic state when [Formula: see text] surpasses 35, and then a phase with strongly nonlinear current-voltage characteristics, suggestive of a pinned Wigner solid, when [Formula: see text] exceeds [Formula: see text] However, our sample makes a transition to an insulating state at [Formula: see text], preceding the onset of the spontaneous ferromagnetism, implying that besides interaction, the role of disorder must also be taken into account in understanding the different phases of a realistic dilute electron system.

Itinerant-Electron Magnetism: The Importance of Many-Body Correlations

Do electrons become ferromagnetic just because of their repulsive Coulomb interaction? Our calculations on the three-dimensional electron gas imply that itinerant ferromagnetism of delocalized electrons without lattice and band structure, the most basic model considered by Stoner, is suppressed due to many-body correlations as speculated already by Wigner, and a possible ferromagnetic transition lowering the density is precluded by the formation of the Wigner crystal.

Landau-Quantized Excitonic Absorption and Luminescence in a Monolayer Valley Semiconductor

We investigate Landau-quantized excitonic absorption and luminescence of monolayer WSe_{2} under magnetic field. We observe gate-dependent quantum oscillations in the bright exciton and trions (or exciton polarons) as well as the dark trions and their phonon replicas. Our results reveal spin- and valley-polarized Landau levels (LLs) with filling factors n=+0, +1 in the bottom conduction band and n=-0 to -6 in the top valence band, including the Berry-curvature-induced n=±0 LLs of massive Dirac fermions. The LL filling produces periodic plateaus in the exciton energy shift accompanied by sharp oscillations in the exciton absorption width and magnitude. This peculiar exciton behavior can be simulated by semiempirical calculations. The experimentally deduced g factors of the conduction band (g∼2.5) and valence band (g∼15) exceed those predicted in a single-particle model (g=1.5, 5.5, respectively). Such g-factor enhancement implies strong many-body interactions in gated monolayer WSe_{2}. The complex interplay between Landau quantization, excitonic effects, and many-body interactions makes monolayer WSe_{2} a promising platform to explore novel correlated quantum phenomena.

Temperature-driven narrowing of the insulating gap as a precursor of the insulator-to-metal transition: Implications for the electronic structure of solids

We present a microscopic picture rationalizing the surprisingly steep decrease in the bandgap with temperature in insulators, crystalline or otherwise. The gap narrowing largely results from fluctuations of long-wavelength optical phonons-when the latter are present-or their disordered analogs if the material is amorphous. We elaborate on this notion to show that possibly with the exception of weakly bound solids made of closed-shell electronic configurations, the existence of an insulating gap or pseudogap in a periodic solid implies that optical phonons must be present, too. This means that in an insulating solid, the primitive cell must have at least two atoms and/or that a charge density wave is present, with the possible exception of weakly bonded solids such as rare-gas or ferromagnetic Wigner crystals. As a corollary, a (periodic) elemental solid held together by nonclosed shell interactions and whose primitive unit contains only one atom will ordinarily be a metal, consistent with observation. Consequences of the present picture for Wigner solids are discussed. A simple field theory of the metal-insulator transition is constructed that directly ties long-wavelength optical vibrations with fluctuations of an order parameter for the metal-insulator transition. The order parameter is shown to have at least two components, yet no Goldstone mode arises as a result of the transition.

Ab Initio Lifetime and Concomitant Double-Excitation Character of Plasmons at Metallic Densities

The accurate calculation of excited state properties of interacting electrons in the condensed phase is an immense challenge in computational physics. Here, we use state-of-the-art equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD) to calculate the dynamic structure factor, which can be experimentally measured by inelastic x-ray and electron scattering. Our calculations are performed on the uniform electron gas at densities corresponding to Wigner-Seitz radii of r_{s}=5, 4, and 3 corresponding to the valence electron densities of common metals. We compare our results to those obtained using the random-phase approximation (RPA), which is known to provide a reasonable description of the collective plasmon excitation and which resums only a small subset of the polarizability diagrams included in EOM-CCSD. We find that EOM-CCSD, instead of providing a perturbative improvement on the RPA plasmon, predicts a many-state plasmon resonance, where each contributing state has a double-excitation character of 80% or more. This finding amounts to an ab initio treatment of the plasmon linewidth, which is in good quantitative agreement with previous diagrammatic calculations, and highlights the strongly correlated nature of lifetime effects in condensed-phase electronic structure theory.

New Reentrant Insulating Phases in Strongly Interacting 2D Systems with Low Disorder

The metal-insulator transition (MIT) in two-dimension (2D) was discovered by Kravchenko et al. more than two decades ago in strongly interacting 2D electrons residing in a Si-metal-oxide-semiconductor field-effect transistor (Si-MOSFET). Its origin remains unresolved. Recently, low magnetic field reentrant insulating phases (RIPs), which dwell between the zero-field (B = 0) metallic state and the integer quantum Hall (QH) states where the Landau-level filling factor υ > 1, have been observed in strongly correlated 2D GaAs hole systems with a large interaction parameter, rs, (~20–40) and a high purity. A new complex phase diagram was proposed, which includes zero-field MIT, low magnetic field RIPs, integer QH states, fractional QH states, high field RIPs and insulating phases (HFIPs) with υ < 1 in which the insulating phases are explained by the formation of a Wigner crystal. Furthermore, evidence of new intermediate phases was reported. This review article serves the purpose of summarizing those recent experimental findings and theoretical endeavors to foster future research efforts.

Phase diagram of a symmetric electron-hole bilayer system: a variational Monte Carlo study

We study the phase diagram of a symmetric electron-hole bilayer system at absolute zero temperature and in zero magnetic field within the quantum Monte Carlo approach. In particular, we conduct variational Monte Carlo simulations for various phases, i.e. the paramagnetic fluid phase, the ferromagnetic fluid phase, the anti-ferromagnetic Wigner crystal phase, the ferromagnetic Wigner crystal phase and the excitonic phase, to estimate the ground-state energy at different values of in-layer density and inter-layer spacing. Slater-Jastrow style trial wave functions, with single-particle orbitals appropriate for different phases, are used to construct the phase diagram in the (r _s , d) plane by finding the relative stability of trial wave functions. At very small layer separations, we find that the fluid phases are stable, with the paramagnetic fluid phase being particularly stable at [Formula: see text] and the ferromagnetic fluid phase being particularly stable at [Formula: see text]. As the layer spacing increases, we first find that there is a phase transition from the ferromagnetic fluid phase to the ferromagnetic Wigner crystal phase when d reaches 0.4 a.u. at r _s   =  20, and before there is a return to the ferromagnetic fluid phase when d approaches 1 a.u. However, for r _s   <  20 and [Formula: see text] a.u., the excitonic phase is found to be stable. We do not find that the anti-ferromagnetic Wigner crystal is stable over the considered range of r _s and d. We also find that as r _s increases, the critical layer separations for Wigner crystallization increase.

Framework for constructing generic Jastrow correlation factors

We have developed a flexible framework for constructing Jastrow factors which allows for the introduction of terms involving arbitrary numbers of particles. The use of various three- and four-body Jastrow terms in quantum Monte Carlo calculations is investigated, including a four-body van der Waals-like term, and anisotropic terms. We have tested these Jastrow factors on one- and two-dimensional homogeneous electron gases, the Be, B, and O atoms, and the BeH, H2O, N2, and H2 molecules. Our optimized Jastrow factors retrieve more than 90% of the fixed-node diffusion Monte Carlo correlation energy in variational Monte Carlo for each system studied.

Connecting the reentrant insulating phase and the zero-field metal-insulator transition in a 2D hole system

We present the transport and capacitance measurements of 10 nm wide GaAs quantum wells with hole densities around the critical point of the 2D metal-insulator transition (critical density p(c) down to 0.8 × 10(10)/cm2, r(s) ∼ 36). For metallic hole density p(c) < p < p(c) + 0.15 × 10(10)/cm2, a reentrant insulating phase (RIP) is observed between the ν = 1 quantum Hall state and the zero-field metallic state and it is attributed to the formation of pinned Wigner crystal. Through studying the evolution of the RIP versus 2D hole density, we show that the RIP is incompressible and continuously connected to the zero-field insulator, suggesting a similar origin for these two phases.

Methods for calculating forces within quantum Monte Carlo simulations

Atomic force calculations within the variational and diffusion quantum Monte Carlo methods are described. The advantages of calculating diffusion quantum Monte Carlo forces with the 'pure' rather than the 'mixed' probability distribution are discussed. An accurate and practical method for calculating forces using the pure distribution is presented and tested for the SiH molecule. The statistics of force estimators are explored and violations of the central limit theorem are found in some cases.

Continuum variational and diffusion quantum Monte Carlo calculations

This topical review describes the methodology of continuum variational and diffusion quantum Monte Carlo calculations. These stochastic methods are based on many-body wavefunctions and are capable of achieving very high accuracy. The algorithms are intrinsically parallel and well suited to implementation on petascale computers, and the computational cost scales as a polynomial in the number of particles. A guide to the systems and topics which have been investigated using these methods is given. The bulk of the article is devoted to an overview of the basic quantum Monte Carlo methods, the forms and optimization of wavefunctions, performing calculations under periodic boundary conditions, using pseudopotentials, excited-state calculations, sources of calculational inaccuracy, and calculating energy differences and forces.