Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regular hexagon”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Evaporative Assembly: Dual-Scale Nanostructures via Evaporative Assembly (Adv. Mater. Interfaces 7/2020)

In article number 1901954, Sunita Srivastava, Surita R. Bhatia, and co-workers demonstrate that evaporative assembly of functionalized nanoparticles can yield striking dual-scale hierarchical structures. Regular microscale stripes of nanoparticle monolayers with hexagonal nanoscale order are obtained on physically and chemically homogeneous substrates through evaporation of a suspension of DNAfunctionalized nanoparticles with a charged shell. The stripe width, spacing and nanoparticle ordering can be controlled by varying nanoparticle concentration and can be described by a simple analytical model. The results indicate that the interplay between “stick-slip” motion of the droplet contact line and Coulombic and steric nanoparticle interactions control the formation of the observed structures. Finally, this work demonstrates a simple cost-effective mask-free method for fabricating nanostructured 2D materials and metasurfaces for applications ranging from energy conversion/storage to optoelectronics and nanophotonics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1992

Artificial intelligence approach to planning the robotic assembly of large tetrahedral truss structures

An assembly planner for tetrahedral truss structures is presented. To overcome the difficulties due to the large number of parts, the planner exploits the simplicity and uniformity of the shapes of the parts and the regularity of their interconnection. The planning automation is based on the computational formalism known as production system. The global data base consists of a hexagonal grid representation of the truss structure. This representation captures the regularity of tetrahedral truss structures and their multiple hierarchies. It maps into quadratic grids and can be implemented in a computer by using a two-dimensional array data structure. By maintaining the multiple hierarchies explicitly in the model, the choice of a particular hierarchy is only made when needed, thus allowing a more informed decision. Furthermore, testing the preconditions of the production rules is simple because the patterned way in which the struts are interconnected is incorporated into the topology of the hexagonal grid. A directed graph representation of assembly sequences allows the use of both graph search and backtracking control strategies.

Homemdemello, Luiz S.↗
NASA NTRS1990

The ringed X-galaxy NGC 7020

The southern SO (sup +) galaxy NGC 7020 presents an unusual morphology: it includes a very regular outer ring which is completely detached and which envelops an inner ring/lens zone with a hexagon surrounding an X shape. The outer ring has a high contrast compared to those usually observed in barred galaxies, yet NGC 7020 is not obviously barred. The morphology of this galaxy poses an interesting puzzle in that the hexagonal/X zone is not a typical type of feature to find in the interior of such a regular ring. Instead, the zone bears a striking resemblance to the edge-on galaxy IC 4767, recently studied by Whitmore and Bell (1988 = WB88) and dubbed by them as the X-galaxy because its inner regions appear to be crossed by two distinct enhancements lined at plus or minus 22 deg with respect to the major axis. The observation of a similar phenomenon in NGC 7020 is interesting because of the suggestion by WB88 that X structures could be related to accretion of matter associated with a merger or tidal encounter between an SO and a small satellite galaxy. If this interpretation is correct for NGC 7020, then it has important implications for the nature of the outer ring. An alternative interpretation is that the inner hexagonal/X zone is a region where resonant periodic orbits in a weak bi-symmetric potential perturbation are influencing the morphology more strongly than might be expected. A brief summary of a more extensive paper (Buta 1990c = B90c) and a few other details concerning this interesting galaxy are given.

Buta, Ronald↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Application of a long short-term memory for deconvoluting conductance contributions at charged ferroelectric domain walls

Ferroelectric domain walls are promising quasi-2D structures that can be leveraged for miniaturization of electronics components and new mechanisms to control electronic signals at the nanoscale. Despite the significant progress in experiment and theory, however, most investigations on ferroelectric domain walls are still on a fundamental level, and reliable characterization of emergent transport phenomena remains a challenging task. Here, we apply a neural-network-based approach to regularize local I ( V )-spectroscopy measurements and improve the information extraction, using data recorded at charged domain walls in hexagonal (Er 0.99 ,Zr 0.01 )MnO 3 as an instructive example. Using a sparse long short-term memory autoencoder, we disentangle competing conductivity signals both spatially and as a function of voltage, facilitating a less biased, unconstrained and more accurate analysis compared to a standard evaluation of conductance maps. The neural-network-based analysis allows us to isolate extrinsic signals that relate to the tip-sample contact and separating them from the intrinsic transport behavior associated with the ferroelectric domain walls in (Er 0.99 ,Zr 0.01 )MnO 3 . Our work expands machine-learning-assisted scanning probe microscopy studies into the realm of local conductance measurements, improving the extraction of physical conduction mechanisms and separation of interfering current signals.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Analytical and Photogrammetric Characterization of a Planar Tetrahedral Truss

Future space science missions are likely to require near-optical quality reflectors which are supported by a stiff truss structure. This support truss should conform closely with its intended shape to minimize its contribution to the overall surface error of the reflector. The current investigation was conducted to evaluate the planar surface accuracy of a regular tetrahedral truss structure by comparing the results of predicted and measured node locations. The truss is a 2-ring hexagonal structure composed of 102 equal-length truss members. Each truss member is nominally 2 meters in length between node centers and is comprised of a graphite/epoxy tube with aluminum nodes and joints. The axial stiffness and the length variation of the truss components were determined experimentally and incorporated into a static finite element analysis of the truss. From this analysis, the root mean square (RMS) surface error of the truss was predicted to be 0.11 mm (0004 in). Photogrammetry tests were performed on the assembled truss to measure the normal displacements of the upper surface nodes and to determine if the truss would maintain its intended shape when subjected to repeated assembly. Considering the variation in the truss component lengths, the measures rms error of 0.14 mm (0.006 in) in the assembled truss is relatively small. The test results also indicate that a repeatable truss surface is achievable. Several potential sources of error were identified and discussed.

Wu, K. Chauncey↗
DOE OSTI2025

Magnetic excitations from the hexagonal spin clusters in the 𝑆 = $\frac{1}{2}$ distorted honeycomb lattice antiferromagnet Cu 2 ⁢(pymca)⁢ 3 ⁢(ClO 4 )

Cu 2 ⁢(pymca) ⁢3 (ClO 4 ) (pymca: pyrimidine-2-carboxylate) consists of a slightly distorted honeycomb lattice of Cu 2+ spins, which shows no long-range magnetic order down to 0.6 K. A magnetization study revealed 1/3 and 2/3 plateau phases [A. Okutani et al., J. Phys. Soc. Jpn. 88, 013703 (2019)], which is not expected for regular honeycomb antiferromagnets. Inelastic neutron scattering experiments were performed using a powder sample to investigate the exchange interactions of this material. The spin excitations from the singlet ground state to the first three triplet states, predicted from the antiferromagnetic hexagonal spin cluster interacting with 3.9 meV, were observed. Using the exact diagonalization methods, the intercluster coupling was estimated from the excitation peak width to be about 20% of the intracluster interaction, which is consistent with the previously reported value. Finally, our exchange path model explains the anisotropic exchange interactions in the distorted honeycomb plane.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Amino Acids Essential for the Assembly of Cellulose Synthase Complexes

Over 40 years ago, Andrew Staehelin's group showed that plant cellulose microfibrils are synthesized by protein complexes arranged in hexagonal arrays called 'particle rosettes'. Plant cellulose synthases (CESAs) differ from their bacterial ancestral forms by insertion of three unique sequences responsible for the assembly into compact multimeric units that aggregate further into rosettes. We used site-directed mutagenesis to replace amino acids within these sequences predicted to be essential for assembly and developed an in vivo method to determine the ability of mutated CesA1 transgenes to complement an Arabidopsis temperature-sensitive root-swelling1 (rsw1) mutant. Replacement of a Cys residue in the Class-Specific Region (CSR) or Pro417 and Arg453 of the Plant-Conserved Region (P-CR) rendered an AtCesA1 transgene unable to complement the rsw1 mutation. Despite an expected role for Arg457 in the trimerization of CESA proteins, AtCesA1 transgenes with Arg457Ala mutations were able to fully restore the wild-type phenotype in rsw1. Staehelin observed a third order of assembly in Micrasterias dendiculata, where hundreds of rosettes form large hexagonal arrays during the synthesis of secondary wall cellulose macrofibrils. We found by SEM of cytosolic faces of Micrasterias cell fragments that these arrays are associated with fields of regularly spaced slime secretion pore complexes.

arabidopsis↗
NASA NTRS2019

Saturn Atmospheric Dynamics One Year After Cassini: Long-Lived Features and Time Variations in the Drift of the Hexagon

We examine Saturn's atmospheric dynamics with observations in the visible range from ground-based telescopes and Hubble Space Telescope (HST). We present a detailed analysis of observations acquired during 2018 obtaining drift rates of major meteorological systems from the equator to the North polar hexagon. A system of polar storms that appeared in the planet in March 2018 and remained active with a complex phenomenology at least until September is analyzed elsewhere [Sanchez-Lavega et al., 2019]. Many of the regular cloud features visible in 2018 are long-lived and can be identified in Saturn images in 2017, and in some cases, for up to a decade using also Cassini ISS images. Without considering the polar storms, the most interesting long-lived cloud systems are: i) A bright white spot in the Equatorial Zone that can be tracked continuously since 2014 with minimal changes in its zonal velocity, which was 444:3 +or- 3:1 ms(exp -1) in 2014 and 452:4 + or - 1:7 ms(exp -1) in 2018. This velocity is remarkably different from the zonal winds at the cloud level at its latitude during the Cassini mission, and is closer to zonal winds obtained at the time of the Voyagers flybys and to zonal winds from Cassini VIMS infrared images of the lower atmosphere. ii) A large long-lived Anticyclone Vortex, here AV, that formed after the Great White Spot of 2010-2011. This vortex has changed significantly in visual contrast, drift rate and latitude with minor changes in size over the last years. iii) A system of subpolar vortices at latitudes 60-65N present at least since 2011. These vortices and additional atmospheric features here studied follow drift rates consistent with zonal winds obtained by Cassini. We also present a study of the positions of the vertices of Saturn's North polar hexagon from 2015 to 2018. These measurements are compared with previous analyses during the Cassini mission (2007-2014), observations with HST in the 90s, and data from the Voyagers in 1980-1981 to explore the long term variability of the hexagon's drift rate. We find variations in the drift rate of the hexagon through these epochs that can not be fit by seasonal changes in the polar area. Instead, the different drift rates reinforce the role of the North Polar Spot that was present in the Voyager epoch and in the early 90s to cause a faster drift rate of the hexagon at that time compared with the current slower one.

R Hueso↗
DOE OSTI2024

Anisotropic physics-regularized interpretable machine learning of microstructure evolution

Anisotropic Physics-Regularized Interpretable Machine Learning Microstructure Evolution (APRIMME) is a general-purpose machine learning solution for grain growth simulations. In prior work, PRIMME employed a deep neural network to predict site-specific migration as a function of its neighboring sites to model normal, isotropic, grain growth behavior. This work aims to extend this method by incorporating grain boundary misorientation-based grain growth behavior. APRIMME is trained on anisotropic simulations created using the Monte Carlo-Potts (MCP) model. Furthermore, the results of this work are compared statistically using grain radius, number of sides per grain, mean neighborhood misorientations, and the standard deviation of triple junction dihedral angles, and are found to match in most cases. The exceptions are small and seem to be related to two causes: (1) the deterministic model of APRIMME is learning from the stochastic simulations of MCP, which seems to accentuate triple junction behaviors; and, (2) a bias against very small grains is made evident in a quicker decrease in grains than expected at the beginning of an APRIMME simulation. APRIMME is also evaluated for its general ability to capture anisotropic grain growth behavior by first investigating different test case initial conditions, including a circle grain, three grain, and hexagonal grain microstructures.

36 MATERIALS SCIENCE↗
NASA NTRS2011

Corrigendum to "Microstructural Characterization of Metal Foams: An Examination of the Applicability of the Theoretical Models for Modeling Foams"

Establishing the geometry of foam cells is useful in developing microstructure-based acoustic and structural models. Since experimental data on the geometry of the foam cells are limited, most modeling efforts use an idealized three-dimensional, space-filling Kelvin tetrakaidecahedron. The validity of this assumption is investigated in the present paper. Several FeCrAlY foams with relative densities varying between 3 and 15% and cells per mm (c.p.mm.) varying between 0.2 and 3.9 c.p.mm. were microstructurally evaluated. The number of edges per face for each foam specimen was counted by approximating the cell faces by regular polygons, where the number of cell faces measured varied between 207 and 745. The present observations revealed that 50-57% of the cell faces were pentagonal while 24-28% were quadrilateral and 15-22% were hexagonal. The present measurements are shown to be in excellent agreement with literature data. It is demonstrated that the Kelvin model, as well as other proposed theoretical models, cannot accurately describe the FeCrAlY foam cell structure. Instead, it is suggested that the ideal foam cell geometry consists of 11 faces with 3 quadrilateral, 6 pentagonal faces and 2 hexagonal faces consistent with the 3-6-2 Matzke cell

Raj. Sai V.↗
NASA NTRS2012

Wargaming in Both Rectilinear and Hexagonal Spaces

There are two main approaches to managing wargame entity interactions (movement, line of sight, area of effect, etc) freespace and gridded In the freespace approach, the units exist as entities in a continuous volume of (usually) Cartesian 3D space. They move in any direction (based on interaction with "terrain" that occupies the same space) and interact with each other based on references and displacements from their position in that space. In the gridded approach, space is broken up into (usually regular) shaped pieces. Units are considered to occupy the entire volume of one of these pieces, movement, line of sight, and other interactions are based on the relationships among the spaces rather than the absolute positions of the units themselves. Both approaches have advantages and drawbacks. The general issue that this discussion has addressed is that there is no "perfect" approach to implementing a wargaming battlespace. Each of them (and this extends to others not discussed) has different sets of advantages and disadvantages. Nothing will change that basic nature of the various approaches, nor would it be desirable to do so. Along with the advantages, the challenges define the feel of the game and focus the thinking of the players on certain aspects and away from others. The proposed approach to combining square and hexagonal approaches, which we will call the rhombus interface, leverages rhombuses constructed from equilateral triangles into which the hexagon can be decomposed to bridge the gap between the approaches, maintain relative consistency between the two as much as possible, and provide most of the feel of the hexagonal approach.

Hoover, Alex↗
NASA NTRS2010

Microstructural Characterization of Metal Foams: An Examination of the Applicability of the Theoretical Models for Modeling Foams

Establishing the geometry of foam cells is useful in developing microstructure-based acoustic and structural models. Since experimental data on the geometry of the foam cells are limited, most modeling efforts use the three-dimensional, space-filling Kelvin tetrakaidecahedron. The validity of this assumption is investigated in the present paper. Several FeCrAlY foams with relative densities varying between 3 and 15 percent and cells per mm (c.p.mm.) varying between 0.2 and 3.9 c.p.mm. were microstructurally evaluated. The number of edges per face for each foam specimen was counted by approximating the cell faces by regular polygons, where the number of cell faces measured varied between 207 and 745. The present observations revealed that 50 to 57 percent of the cell faces were pentagonal while 24 to 28 percent were quadrilateral and 15 to 22 percent were hexagonal. The present measurements are shown to be in excellent agreement with literature data. It is demonstrated that the Kelvin model, as well as other proposed theoretical models, cannot accurately describe the FeCrAlY foam cell structure. Instead, it is suggested that the ideal foam cell geometry consists of 11 faces with 3 quadrilateral, 6 pentagonal faces and 2 hexagonal faces consistent with the 3-6-2 cell.

Raj, S. V.↗
NASA NTRS2011

Microstructural Characterization of Metal Foams: An Examination of the Applicability of the Theoretical Models for Modeling Foams

Establishing the geometry of foam cells is useful in developing microstructure-based acoustic and structural models. Since experimental data on the geometry of the foam cells are limited, most modeling efforts use an idealized three-dimensional, space-filling Kelvin tetrakaidecahedron. The validity of this assumption is investigated in the present paper. Several FeCrAlY foams with relative densities varying between 3 and 15 percent and cells per mm (c.p.mm.) varying between 0.2 and 3.9 c.p.mm. were microstructurally evaluated. The number of edges per face for each foam specimen was counted by approximating the cell faces by regular polygons, where the number of cell faces measured varied between 207 and 745. The present observations revealed that 50 to 57 percent of the cell faces were pentagonal while 24 to 28 percent were quadrilateral and 15 to 22 percent were hexagonal. The present measurements are shown to be in excellent agreement with literature data. It is demonstrated that the Kelvin model, as well as other proposed theoretical models, cannot accurately describe the FeCrAlY foam cell structure. Instead, it is suggested that the ideal foam cell geometry consists of 11 faces with three quadrilateral, six pentagonal faces and two hexagonal faces consistent with the 3-6-2 Matzke cell. A compilation of 90 years of experimental data reveals that the average number of cell faces decreases linearly with the increasing ratio of quadrilateral to pentagonal faces. It is concluded that the Kelvin model is not supported by these experimental data.

Raj, S. V.↗
NASA NTRS1974

X-ray powder diffraction study of poly/carbon monofluoride/, CF/1.12/

Data from X-ray diffraction studies of the poly(carbon monofluoride) with empirical formula CF(1.09-1.15) are reported, and possible intercalation arrangements for the substance are discussed. The data do not conform to true hexagonal symmetry, indicating that the carbon atoms are not coplanar. Each bond angle of carbon is 118.8 deg, and the carbon-carbon distance is 1.47 A. The interlayer distance is 5.76 A. A total absence of (hkl) reflections in the X-ray pattern shows that the separate CF layers are not regularly arranged with respect to one another.

Mahajan, V. K.↗
DOE OSTI2020

Dual‐Scale Nanostructures via Evaporative Assembly

Abstract Dual‐scale hierarchical structures with regular microscale patterns and varying degree of nanoscale crystalline order are synthesized on physically and chemically homogeneous substrates by evaporative self‐assembly with a suspension of DNA‐functionalized nanoparticles (NPs) with a charged core shell. For a certain NP concentration range, periodic concentric rings in a stripe‐like micropattern are produced over macroscale surface areas by an NP monolayer with hexagonal lattice structure at the nanoscale. The stripe width, spacing, and nanoparticle ordering can be controlled by varying the NP concentration. The results indicate that the interplay between “stick‐slip” motion of the droplet contact line and coulombic and steric NP interactions control the formation of the observed structures. A simple analytical model is proposed to account for the experimental observations and guide the future design of different nanostructure morphologies. This work demonstrates a simple cost‐effective mask‐free method for fabricating large‐area nanostructured 2D materials and metasurfaces for applications ranging from energy conversion/storage to optoelectronics and nanophotonics.

Srivastava, Sunita↗
DOE OSTI2020

Influence of crystal structure on constitutive anisotropy of silica sand at particle-scale

In spite of recent advanced geotechnical research regarding the influence of particle-scale morphology and fabric on the constitutive behavior of silica sands, there is still no published literature that investigated the influence of crystal structure of sand particles on their constitutive behavior. The main constituent of silica sand particles is quartz mineral, which has a trigonal crystal system and hexagonal lattice. This paper investigates the influence of quartz crystal structure on the constitutive response of silica sand particles using 3D X-ray diffraction (3DXRD), synchrotron micro-computed tomography (SMT), and 3D finite element (FE) analysis. Here, unconfined uniaxial compression experiments on single particles of synthetic silica cubes exposed the notion of constitutive anisotropy caused by the crystal structure of quartz. Since the silica cubes have simple regular geometry (1 mm) 3 , they were loaded in different directions with respect to their crystal local planes. The crystal structure was also identified for a natural silica sand specimen composed of 2705 particles and loaded in confined uniaxial compression. Physical limitations of the experimental measurements made a strong case for the use of 3D FE analysis to study how the change in crystal local orientation of individual sand particles can fundamentally produce anisotropy in the constitutive response of natural sand specimens. The 3D FE analysis was validated to accurately model the crystal-based constitutive anisotropy in silica particles using the results of the silica cube experiments. The analysis was then executed on 3D FE meshes that closely matched the physical shape and fabric of the 2705 particles in the natural sand specimen. Both the experiments and 3D FE analysis showed that the crystal structure of quartz essentially causes a directional anisotropy in the constitutive behavior of silica particles.

58 GEOSCIENCES↗
DOE OSTI2021

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↗