Search NASA⌕ Search

SEARCH · Search NASA

Results for “crystallization kinetics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

Engineering and characterization of carbohydrate‐binding modules for imaging cellulose fibrils biosynthesis in plant protoplasts

Abstract Carbohydrate binding modules (CBMs) are noncatalytic domains that assist tethered catalytic domains in substrate targeting. CBMs have therefore been used to visualize distinct polysaccharides present in the cell wall of plant cells and tissues. However, most previous studies provide a qualitative analysis of CBM‐polysaccharide interactions, with limited characterization of engineered tandem CBM designs for recognizing polysaccharides like cellulose and limited application of CBM‐based probes to visualize cellulose fibrils synthesis in model plant protoplasts with regenerating cell walls. Here, we examine the dynamic interactions of engineered type‐A CBMs from families 3a and 64 with crystalline cellulose‐I and phosphoric acid swollen cellulose. We generated tandem CBM designs to determine various characteristic properties including binding reversibility toward cellulose‐I using equilibrium binding assays. To compute the adsorption ( nk on ) and desorption ( k off ) rate constants of single versus tandem CBM designs toward nanocrystalline cellulose, we employed dynamic kinetic binding assays using quartz crystal microbalance with dissipation. Our results indicate that tandem CBM3a exhibited the highest adsorption rate to cellulose and displayed reversible binding to both crystalline/amorphous cellulose, unlike other CBM designs, making tandem CBM3a better suited for live plant cell wall biosynthesis imaging applications. We used several engineered CBMs to visualize Arabidopsis thaliana protoplasts with regenerated cell walls using confocal laser scanning microscopy and wide‐field fluorescence microscopy. Lastly, we also demonstrated how CBMs as probe reagents can enable in situ visualization of cellulose fibrils during cell wall regeneration in Arabidopsis protoplasts.

09 BIOMASS FUELS↗

Structural, functional, and regulatory evaluation of a cysteine post-translationally modified Gcn5-related N -acetyltransferase

Polyamines within the cell are tightly regulated by spermidine/spermine N-acetyltransferase (SSAT) enzymes. While several SSATs have been investigated in different bacterial species, there is still a significant gap in knowledge about which proteins are functional SSATs in many organisms. For example, while it is known that Pseudomonas aeruginosa synthesizes the polyamine spermidine, the SSAT that acetylates this molecule and its importance in regulating intracellular polyamines remains unknown. We previously identified a candidate Gcn5-related N-acetyltransferase (GNAT) protein from P. aeruginosa (PA2271) that could fulfill this role since it acetylates spermidine, but no further studies were conducted. Here, we explored the structure/function relationship of the PA2271 protein by determining its X-ray crystal structure and performing enzyme kinetics assays. We also identified active site residues that are essential for catalysis and substrate binding. As the study progressed, we encountered results that led us to explore the importance of four cysteine residues on enzyme activity and disulfide bond formation or modification of cysteine residues. We found these cysteine residues in PA2271 are important for protein solubility and activity, and there is an interrelationship between cysteine residues that contribute to these effects. Furthermore, we also found disulfide bonds could form between C 121 and C 165 and speculate that these residues may contribute to redox regulation of PA2271 protein activity.

59 BASIC BIOLOGICAL SCIENCES↗

Size-Resolved Shape Evolution in Inorganic Nanocrystals Captured via High-Throughput Deep Learning-Driven Statistical Characterization

Precise size and shape control in nanocrystal synthesis is essential for utilizing nanocrystals in various industrial applications, such as catalysis, sensing, and energy conversion. However, traditional ensemble measurements often overlook the subtle size and shape distributions of individual nanocrystals, hindering the establishment of robust structure–property relationships. In this study, we uncover intricate shape evolutions and growth mechanisms in Co 3 O 4 nanocrystal synthesis at a subnanometer scale, enabled by deep-learning-assisted statistical characterization. By first controlling synthetic parameters such as cobalt precursor concentration and water amount then using high resolution electron microscopy imaging to identify the geometric features of individual nanocrystals, this study provides insights into the interplay between synthesis conditions and the sizedependent shape evolution in colloidal nanocrystals. Utilizing population-wide imaging data encompassing over 441,067 nanocrystals, we analyze their characteristics and elucidate previously unobserved size-resolved shape evolution. This high-throughput statistical analysis is essential for representing the entire population accurately and enables the study of the size dependency of growth regimes in shaping nanocrystals. Our findings provide experimental quantification of the growth regime transition based on the size of the crystals, specifically (i) for faceting and (ii) from thermodynamic to kinetic, as evidenced by transitions from convex to concave polyhedral crystals. Additionally, we introduce the concept of an “onset radius,” which describes the critical size thresholds at which these transitions occur. This discovery has implications beyond achieving nanocrystals with desired morphology; it enables finely tuned correlation between geometry and material properties, advancing the field of colloidal nanocrystal synthesis and its applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Machine-learning accelerated geometry optimization in molecular simulation

Geometry optimization is an important part of both computational materials and surface science because it is the path to finding ground state atomic structures and reaction pathways. These properties are used in the estimation of thermodynamic and kinetic properties of molecular and crystal structures. This process is slow at the quantum level of theory because it involves an iterative calculation of forces using quantum chemical codes such as density functional theory (DFT), which are computationally expensive and which limit the speed of the optimization algorithms. It would be highly advantageous to accelerate this process because then one could do either the same amount of work in less time or more work in the same time. Here, we provide a neural network (NN) ensemble based active learning method to accelerate the local geometry optimization for multiple configurations simultaneously. We illustrate the acceleration on several case studies including bare metal surfaces, surfaces with adsorbates, and nudged elastic band for two reactions. In all cases, the accelerated method requires fewer DFT calculations than the standard method. In addition, we provide an Atomic Simulation Environment (ASE)-optimizer Python package to make the usage of the NN ensemble active learning for geometry optimization easier.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational discovery of two-dimensional rare-earth iodides: promising ferrovalley materials for valleytronics

Two-dimensional Ferrovalley materials with intrinsic valley polarization are rare but highly promising for valley-based nonvolatile random access memory and valley filter devices. These ferromagnetic materials exhibit valleys at or near the Fermi level with intrinsic magnetism. The strong coupling between magnetism and spin–orbit coupling induces intrinsic valley polarization. Using Kinetically Limited Minimization, an unconstrained crystal structure prediction algorithm, and prototype sampling based on first-principles calculations, we have discovered new Ferrovalley materials, rare-earth iodides RI 2 , where R is a rare-earth element belonging to Sc, Y, or La-Lu, and I is Iodine. The rare-earth iodides are layered and demonstrate either 2H, 1T, or 1T d phase as the ground state in bulk, analogous to transition metal dichalcogenides (TMDCs). The calculated exfoliation energy of monolayers (MLs) is comparable to that of graphene and TMDCs, suggesting possible experimental synthesis. The MLs in the 2H phase exhibit ferromagnetism due to unpaired electrons in d and f orbitals. Throughout the rare-earth series, d bands have valley polarization at K and $\overline{K}$ points in the Brillouin zone in the vicinity of the Fermi level. Large intrinsic valley polarization in the range of 15–143 meV without external stimuli is observed in these Ferrovalley materials, which can be enhanced further by applying an in-plane bi-axial strain. These valleys can selectively be probed and manipulated for information storage and processing, potentially offering superior performance beyond conventional electronics and spintronics. Here we further show that the 2H ferromagnetic phase of RI 2 MLs possesses non-zero Berry curvature and exhibits anomalous valley Hall effect with considerable anomalous Hall conductivity. Our work will incite exploratory synthesis of the predicted Ferrovalley materials and their application in valleytronics and beyond.

2D materials↗

A sorghum ascorbate peroxidase with four binding sites has activity against ascorbate and phenylpropanoids

Abstract In planta, H2O2 is produced as a by-product of enzymatic reactions and during defense responses. Ascorbate peroxidase (APX) is a key enzyme involved in scavenging cytotoxic H2O2. Here, we report the crystal structure of cytosolic APX from sorghum (Sorghum bicolor) (Sobic.001G410200). While the overall structure of SbAPX was similar to that of other APXs, SbAPX uniquely displayed four bound ascorbates rather than one. In addition to the ɣ-heme pocket identified in other APXs, ascorbates were bound at the δ-meso and two solvent-exposed pockets. Consistent with the presence of multiple binding sites, our results indicated that the H2O2-dependent oxidation of ascorbate displayed positive cooperativity. Bound ascorbate at two surface sites established an intricate proton network with ascorbate at the ɣ-heme edge and δ-meso sites. Based on crystal structures, steady-state kinetics, and site-directed mutagenesis results, both ascorbate molecules at the ɣ-heme edge and the one at the surface are expected to participate in the oxidation reaction. We provide evidence that the H2O2-dependent oxidation of ascorbate by APX produces a C2-hydrated bicyclic hemiketal form of dehydroascorbic acid at the ɣ-heme edge, indicating two successive electron transfers from a single-bound ascorbate. In addition, the δ-meso site was shared with several organic compounds, including p-coumaric acid and other phenylpropanoids, for the potential radicalization reaction. Site-directed mutagenesis of the critical residue at the ɣ-heme edge (R172A) only partially reduced polymerization activity. Thus, APX removes stress-generated H2O2 with ascorbates, and also uses this same H2O2 to potentially fortify cell walls via oxidative polymerization of phenylpropanoids in response to stress.

Plant Sciences↗

Gas Hydrate Film Growth in Microfluidic Channels for Carbon Dioxide Capture and Sequestration Applications

Gas or clathrate hydrates are a solid, crystalline compound composed of water and guest molecules that typically form at high pressure and low temperature conditions. Carbon dioxide (CO2) hydrates may be involved in several carbon dioxide capture and sequestration (CCS) applications, including CO2 pipeline transportation and CO2 offshore sequestration. Within these applications, the formation mechanism and kinetics must be well understood to manage the CCS processes, either by preventing or promoting hydrate formation. In this work, a high-pressure glass microfluidic reactor is used in tandem with visual microscopy and in-situ Raman spectroscopy to study both the morphological and kinetic behavior of gas hydrate crystals. Subcooling, pressure, and CO2 flow rate are investigated for their impact on the thickening behavior of pure CO2 hydrates, with flow rate being the only parameter to have a significant effect. Visual and Raman spectroscopy evidence show that both a dense hydrate layer and a porous hydrate layer form, and the latter may provide a path for mass transfer to continue hydrate crystallization. A first principles mass transfer model is developed to describe CO2 hydrate crystal thickening at the interface between gas and water. The impacts of gas impurities and channel wettability are also studied. This method is further applied to investigate the conversion of methane hydrate to CO2 hydrate for combined energy recovery and methane hydrate formation. The authors acknowledge the US Department of Energy Basic Energy Science award # DE-SC0022162.

Wadsworth, Lindsey [Colorado School of Mines, Gold↗

Enzyme kinetics by GH7 cellobiohydrolases on chromogenic substrates is dictated by non‐productive binding: insights from crystal structures and MD simulation

Cellobiohydrolases (CBHs) in the glycoside hydrolase family 7 (GH7) ( EC3.2.1.176 ) are the major cellulose degrading enzymes both in industrial settings and in the context of carbon cycling in nature. Small carbohydrate conjugates such as p‐ nitrophenyl‐β‐ d ‐cellobioside (pNPC), p‐ nitrophenyl‐β‐ d ‐lactoside (pNPL) and methylumbelliferyl‐β‐ d ‐cellobioside have commonly been used in colorimetric and fluorometric assays for analysing activity of these enzymes. Despite the similar nature of these compounds the kinetics of their enzymatic hydrolysis vary greatly between the different compounds as well as among different enzymes within the GH7 family. Through enzyme kinetics, crystallographic structure determination, molecular dynamics simulations, and fluorometric binding studies using the closely related compound o‐nitrophenyl‐β‐ d ‐cellobioside (oNPC), in this work we examine the different hydrolysis characteristics of these compounds on two model enzymes of this class, TrCel7A from Trichoderma reesei and PcCel7D from Phanerochaete chrysosporium . Protein crystal structures of the E212Q mutant of TrCel7A with pNPC and pNPL, and the wildtype TrCel7A with oNPC, reveal that non‐productive binding at the product site is the dominating binding mode for these compounds. Enzyme kinetics results suggest the strength of non‐productive binding is a key determinant for the activity characteristics on these substrates, with PcCel7D consistently showing higher turnover rates ( k cat ) than TrCel7A, but higher Michaelis–Menten ( K M ) constants as well. Furthermore, oNPC turned out to be useful as an active‐site probe for fluorometric determination of the dissociation constant for cellobiose on TrCel7A but could not be utilized for the same purpose on PcCel7D, likely due to strong binding to an unknown site outside the active site.

59 BASIC BIOLOGICAL SCIENCES↗

In Situ Synchrotron X-ray Scattering Investigation of Cathodic ZIF-8 Deposition on Graphite Using 3D-Printed Cells

Electrochemical approaches provide unique means to fabricate thin films of metal organic frameworks (MOFs). However, the kinetics of electrochemical MOF deposition has not been quantitatively investigated so far. In this study, we report the first in situ measurements of electrochemical MOF growth using transmission synchrotron X-ray scattering. Electrochemical cells based on poly(lactic acid) with two windows were fabricated using fused-deposition modelling. Here, the resulting 3D-printed cells, the surface of which was coated with paraffin wax to prevent solvent percolation through the polymer material, were used to monitor the cathodic growth of zeolitic imidazolate framework-8 (ZIF-8) on graphite in a methanol solution containing ZnCl 2 and 2-methylimidazole (Hmim) at different cathodic potentials. The resulting time-resolved X-ray diffraction data showed a gradual increase in crystal size with negligible changes in crystal orientation during the cathodic ZIF-8 deposition. More importantly, the time-resolved data provided a means to quantitatively assess the kinetics of the cathodic ZIF-8 growth using the Gualtieri model, revealing that the cathodic potential and Hmim concentration affected crystal growth kinetics, but not nucleation kinetics. These ZIF-8 samples exhibited changes in X-ray diffraction pattern after being washed with methanol and dried in air, indicating that in situ measurements are essential to investigate mechanisms behind MOF electrodeposition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydration structure and water exchange kinetics at xenotime–water interfaces: implications for rare earth minerals separation

Hydration of surface ions gives rise to structural heterogeneity and variable exchange kinetics of water at complex mineral–water interfaces. In this work, we employ ab initio molecular dynamics (AIMD) simulations and water adsorption calorimetry to examine the aqueous interfaces of xenotime, a phosphate mineral that contains predominantly Y 3+ and heavy rare earth elements. Consistent with natural crystal morphology, xenotime is predicted to have a tetragonal prismatic shape, dominated by the {100} surface. Hydration of this surface induces multilayer interfacial water structures with distinct OH orientations, which agrees with recent crystal truncation rod measurements. The exchange kinetics between two adjacent water layers exhibits a wide range of underlying timescales (5–180 picoseconds), dictated by ion–water electrostatics. Adsorption of a bidentate hydroxamate ligand reveals that {100} xenotime surface can only accommodate monodentate coordination with water exchange kinetics strongly depending on specific ligand orientation, prompting us to reconsider traditional strategies for selective separation of rare-earth minerals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon

The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.

Science & Technology - Other Topics↗

Mapping the Complete Reaction Energy Landscape of a Metal–Organic Framework Phase Transformation

Crystalline materials undergo valuable phase transformations, and the energetic processes that underlie these transformations can be fully characterized through a combination of thermodynamic and kinetic studies. Here, we report the first complete reaction energy landscape of metal–organic framework (MOF) interpenetration, specifically in the phase transformation of NU-1200 to its doubly interpenetrated counterpart, STA-26. We characterized the thermodynamics of this phase transformation by pairing experiments with density functional theory (DFT) calculations. This analysis revealed that factors such as the increase in crystal density likely drive Zr- and Hf-NU-1200 to STA-26 interpenetration, while other chemical interactions such as steric repulsions prevent Th-NU-1200 from interpenetrating. Using time-resolved in situ X-ray diffraction, we monitored phase transformation reaction profiles and extracted quantitative kinetic information using the Avrami-Erofe’ev model. As a result, we obtained activation energies for the Zr- and Hf-NU-1200 transformations to Zr- and Hf-STA-26, respectively, revealing slower phase change kinetics for MOFs with stronger bonds. Finally, we paired the kinetic data with experimental observations to classify the mechanistic model of this phase transformation as partial dissolution. Here, we anticipate that this thermodynamic, kinetic, and mechanistic understanding will broadly inform further studies on the energetics of crystallization.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Pulling thin single crystal silicon wafers from a melt: The new leading-edge solar substrate

The Floating Silicon Method (FSM) has been established as a viable, stable method for growing single crystal ribbons directly from a silicon melt. With intense helium jet cooling to drive the linear progress of a [111] facet, pulled in the <110> direction, ribbons in the 0.6 – 3.0 mm thickness range can be grown at linear growth rates from 0.3 mm/s to >6 mm/s as reported in the literature. We report on recent progress towards growing (100) oriented ribbons with a net thickness of less than 200 microns and a ribbon width up to 18 cm using a stable, continuous process in the Leading Edge prototype furnace. Here, the 3D details of the single crystal growth are explained by extending the mechanics of the Limit Cycle Theory, with novel Internal Side Effect morphology described by a proposed Facet Flow Theory. Grown-in crystalline defect distributions are described as well as values of critical impurities like oxygen, carbon, dopants, and metals that are relevant for use as wafers for solar cells.

14 SOLAR ENERGY↗

Correlating wavelength dependence in LiMn2O4 cathode photo-accelerated fast charging with deformations in local structure

The growth in electrified transportation has benefited from the massive worldwide research efforts used to discover and improve electrode materials and electrolytes. Nevertheless, lithium-ion batteries still suffer from a slow-charging limitation. Recently, it has been demonstrated that white light illumination of LiMn2O4 provokes faster charging, improving the kinetics of delithiation without the use of nanostructured active materials. In this work, we probe the mechanism of photo-accelerated fast charging and show that Mn d-d electronic transitions occurring under red light illumination are largely responsible for the increased charging rate. It is further demonstrated through X-ray absorption spectroscopy methods that LiMn2O4 Mn-Mn bond distances shorten after d-electron excitation. The shrinkage in the crystal volume beneficially contributes to delithiation kinetics by lowering the resistance to lithium-ion conduction. Advanced materials that can absorb light to modulate their structure may provide us with a new mechanistic pathway to pursue for increasing charge transfer rates.

25 ENERGY STORAGE↗

Unravelling crystal growth of nanoparticles

Crystal growth from nanoscale constituents is a ubiquitous phenomenon in biology, geology and materials science. Numerous studies have focused on understanding the onset of nucleation and on producing high-quality crystals by empirically sampling constituents with different attributes and varying the growth conditions. However, the kinetics of post-nucleation growth processes, an important determinant of crystal morphology and properties, have remained underexplored due to experimental challenges associated with real-space imaging at the nanoscale. Here we report the imaging of the crystal growth of nanoparticles of different shapes using liquid-phase transmission electron microscopy, resolving both lateral and perpendicular growth of crystal layers by tracking individual nanoparticles. We observe that these nanoscale systems exhibit layer-by-layer growth, typical of atomic crystallization, as well as rough growth prevalent in colloidal systems. Surprisingly, the lateral and perpendicular growth modes can be independently controlled, resulting in two mixed crystallization modes that, until now, have received only scant attention. Combining analytical considerations with molecular dynamics and kinetic Monte Carlo simulations, we develop a comprehensive framework for our observations, which are fundamentally determined by the size and shape of the building blocks. Furthermore, these insights unify the understanding of crystal growth across four orders of magnitude in particle size and suggest novel pathways to crystal engineering.

Imaging techniques↗

Rotational Dynamics of an Amphidynamic Zirconium Metal–Organic Framework Determined by Dielectric Spectroscopy

A zirconium metal–organic framework with a difluorophenylene rotator bearing a permanent electric dipole of ~3.2 D was synthesized, and its rotational motion was analyzed by temperature- and frequency-dependent broadband dielectric spectroscopy. While solid-state NMR confirms fast rotation qualitatively, the dissipation factors measured between 113 and 153 K suggested an activation energy E a = 2.6 kcal/mol, but deviations from a single Debye relaxation suggested a dynamic process that cannot be accounted for by a well-defined potential with a single activation barrier. The dynamic heterogeneity of the dipolar rotor was confirmed by analysis in terms of a Cole–Cole relaxation, which suggested a mean barrier of ~1.9 kcal/mol, with a heterogeneity that decreases as temperature increases. Based on the single-crystal structure, we propose that the kinetic heterogeneity results from a temperature-dependent potential where rotation motion is mediated by the escape of the rotator from an energy well created by a double Ph–H···F–Ph hydrogen bond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Statistically averaged molecular dynamics simulations of hydrogen diffusion in magnesium and magnesium hydrides

Magnesium has a different crystal structure from its dihydride with hydrogenation leading to a phase transition from the hexagonal closely packed Mg into a tetragonal α-MgH 2 rutile type structure. Such materials exhibit complex hydrogen uptake and release kinetics because hydrogen diffusivities significantly change when the crystal structure changes. To provide a foundational understanding of (de)hydrogenation kinetics that is applicable to all stages of the reaction, we performed statistically averaged molecular dynamics simulations to derive hydrogen diffusivities as a function of temperature and hydrogen content for both magnesium and magnesium hydride. Furthermore, our studies confirm that hydrogen diffusivities in magnesium hydride are much lower than in magnesium, in agreement with experimental data. Additionally, we observe that in either magnesium or magnesium hydride, higher hydrogen compositions result in reduced diffusivities. The latter was not revealed by prior experiments, which were conducted at fixed hydrogen composition. Finally, we discover a non-Arrhenius behavior in magnesium hydride. The physical origin of this behavior is also discussed.

36 MATERIALS SCIENCE↗

Bilayer Wigner crystals in a transition metal dichalcogenide heterostructure

One of the first theoretically predicted manifestations of strong interactions in many-electron systems was the Wigner crystal, in which electrons crystallize into a regular lattice. The crystal can melt via either thermal or quantum fluctuations. Quantum melting of the Wigner crystal is predicted to produce exotic intermediate phases and quantum magnetism because of the intricate interplay of Coulomb interactions and kinetic energy. However, studying two-dimensional Wigner crystals in the quantum regime has often required a strong magnetic field or a moiré superlattice potential, thus limiting access to the full phase diagram of the interacting electron liquid. Here we report the observation of bilayer Wigner crystals without magnetic fields or moiré potentials in an atomically thin transition metal dichalcogenide heterostructure, which consists of two MoSe 2 monolayers separated by hexagonal boron nitride. We observe optical signatures of robust correlated insulating states at symmetric (1:1) and asymmetric (3:1, 4:1 and 7:1) electron doping of the two MoSe 2 layers at cryogenic temperatures. Here, we attribute these features to bilayer Wigner crystals composed of two interlocked commensurate triangular electron lattices, stabilized by inter-layer interaction. The Wigner crystal phases are remarkably stable, and undergo quantum and thermal melting transitions at electron densities of up to 6 × 10 12 per square centimetre and at temperatures of up to about 40 kelvin. Our results demonstrate that an atomically thin heterostructure is a highly tunable platform for realizing many-body electronic states and probing their liquid–solid and magnetic quantum phase transitions.

36 MATERIALS SCIENCE↗