Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “theoretical limits”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Electron heating in 2D: Combining Fermi–Ulam acceleration and magnetic-moment non-adiabaticity in a mirror-configuration plasma

We analyze a new mechanism for the creation and confinement of energetic electrons in a mirror-configuration plasma. A Fermi–Ulam-type process, driven by end-localized coherent electrostatic oscillations, provides axial acceleration, while a natural non-adiabaticity of μ provides phase decorrelation and energy isotropization. This novel 2D combination causes the electron energy distribution function, calculated with a diffusive-loss model, to assume a Maxwellian shape with the μ non-adiabaticity, reducing loss-cone escape and annulling the absolute-barrier energy-limiting Chirikov criterion of lower dimensional models. Here, the theoretical predictions are compared with data from an experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Mind the gap: Bridging the divide between AI aspirations and the reality of autonomous microscopy

What does materials science look like in the “Age of Artificial Intelligence?” Each material’s domain—synthesis, characterization, and modeling—has a different answer to this question, motivated by unique challenges and constraints. This work focuses on the tremendous potential of autonomous characterization within electron microscopy. We present our recent advancements in developing domain-aware, multimodal models for microscopy analysis capable of describing complex atomic systems. We then address the critical gap between the theoretical promise of autonomous microscopy and its current practical limitations, showcasing recent successes while highlighting the necessary developments to achieve robust, real-world autonomy.

2D materials↗
DOE OSTI2020

Magnetic, thermal, and electronic-transport properties of EuMg 2 Bi 2 single crystals

The trigonal compound EuMg 2 Bi 2 has recently been discussed in terms of its topological band properties. These are intertwined with its magnetic properties. In this paper detailed studies of the magnetic, thermal, and electronic-transport properties of EuMg 2 Bi 2 single crystals are presented. The Eu + 2 spins 7/2 in EuMg 2 Bi 2 exhibit an antiferromagnetic (AFM) transition at a temperature T N = 6.7 K, as previously reported. By analyzing the anisotropic magnetic susceptibility χ data below T N in terms of molecular-field theory (MFT), the AFM structure is inferred to be a c -axis helix, where the ordered moments in the hexagonal a b -plane layers are aligned ferromagnetically in the a b plane with a turn angle between the moments in adjacent moment planes along the c axis of ≈ 120 ° . An alternate but less likely magnetic structure is a planar structure with nearest-neighbor Eu spins aligned at ≈ 120 ° with respect to each other, where these ordered-moment layers are stacked along the c axis. The magnetic heat capacity exhibits a λ anomaly at T N with evidence of dynamic short-range magnetic fluctuations both above and below T N . The high- T limit of the magnetic entropy is close to the theoretical value for spins 7/2. The in-plane electrical resistivity ρ ( T ) data indicate metallic character with a mild and disorder-sensitive upturn below T min = 23 K. An anomalous rapid drop in ρ ( T ) on cooling below T N as found in zero field is replaced by a two-step decrease in magnetic fields. The ρ ( T ) measurements also reveal an additional transition below T N in applied fields of unknown origin that is not observed in the other measurements and may be associated with an incommensurate to commensurate AFM transition. The dependence of T N on the c -axis magnetic field H ⊥ was derived from the field-dependent χ ( T ) , C p ( T ) , and ρ ( T ) measurements. This T N ( H ⊥ ) was found to be consistent with the prediction of MFT for a c -axis helix with S = 7 / 2 and was used to generate a phase diagram in the H ⊥ – T plane.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantum simulation of nonequilibrium dynamics and thermalization in the Schwinger model

In this work, we present simulations of nonequilibrium dynamics of quantum field theories on digital quantum computers. As a representative example, we consider the Schwinger model, a (1+1)-dimensional U(1) gauge theory, coupled through a Yukawa-type interaction to a thermal environment described by a scalar field theory. We use the Hamiltonian formulation of the Schwinger model discretized on a spatial lattice. With the thermal scalar fields traced out, the Schwinger model can be treated as an open quantum system and its real-time dynamics are governed by a Lindblad equation in the Markovian limit. The interaction with the environment ultimately drives the system to thermal equilibrium. In the quantum Brownian motion limit, the Lindblad equation is related to a field theoretical Caldeira-Leggett equation. By using the Stinespring dilation theorem with ancillary qubits, we perform studies of both the nonequilibrium dynamics and the preparation of a thermal state in the Schwinger model using IBM's simulator and quantum devices. The real-time dynamics of field theories as open quantum systems and the thermal state preparation studied here are relevant for a variety of applications in nuclear and particle physics, quantum information and cosmology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A round-robin approach provides a detailed assessment of biomolecular small-angle scattering data reproducibility and yields consensus curves for benchmarking

Through an expansive international effort that involved data collection on 12 small-angle X-ray scattering (SAXS) and four small-angle neutron scattering (SANS) instruments, 171 SAXS and 76 SANS measurements for five proteins (ribonuclease A, lysozyme, xylanase, urate oxidase and xylose isomerase) were acquired. From these data, the solvent-subtracted protein scattering profiles were shown to be reproducible, with the caveat that an additive constant adjustment was required to account for small errors in solvent subtraction. Further, the major features of the obtained consensus SAXS data over the q measurement range 0–1 Å −1 are consistent with theoretical prediction. The inherently lower statistical precision for SANS limited the reliably measured q -range to <0.5 Å −1 , but within the limits of experimental uncertainties the major features of the consensus SANS data were also consistent with prediction for all five proteins measured in H 2 O and in D 2 O. Thus, a foundation set of consensus SAS profiles has been obtained for benchmarking scattering-profile prediction from atomic coordinates. Additionally, two sets of SAXS data measured at different facilities to q > 2.2 Å −1 showed good mutual agreement, affirming that this region has interpretable features for structural modelling. SAS measurements with inline size-exclusion chromatography (SEC) proved to be generally superior for eliminating sample heterogeneity, but with unavoidable sample dilution during column elution, while batch SAS data collected at higher concentrations and for longer times provided superior statistical precision. Careful merging of data measured using inline SEC and batch modes, or low- and high-concentration data from batch measurements, was successful in eliminating small amounts of aggregate or interparticle interference from the scattering while providing improved statistical precision overall for the benchmarking data set.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Large-signal Stability Analysis of Grid-forming Inverters with Equivalent-circuit Models

Here, this paper proposes an energy function-based direct method for large-signal stability assessment of grid-forming (GFM) inverters leveraging an equivalent-circuit representation of all involved control- and physical-layer dynamics. Three different primary controls, a standard inner-current outer-voltage cascaded-control architecture, output LCL filter, and reference-current saturation limiting are featured in the modeling and analysis framework. A composite energy function for the GFM inverter is obtained by summing up individual energy contributions gleaned from the circuit representation. The approach can readily be generalized to different primary controls, output-filter arrangements, and current limiters since it is based on a circuit-theoretic foundation. Numerical simulations validate the efficacy of the approach in estimating the critical clearing time following a large-signal disturbance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Simultaneous estimation of multiple eigenvalues with short-depth quantum circuit on early fault-tolerant quantum computers

We introduce a multi-modal, multi-level quantum complex exponential least squares (MM-QCELS) method to simultaneously estimate multiple eigenvalues of a quantum Hamiltonian on early fault-tolerant quantum computers. Our theoretical analysis demonstrates that the algorithm exhibits Heisenberg-limited scaling in terms of circuit depth and total cost. Notably, the proposed quantum circuit utilizes just one ancilla qubit, and with appropriate initial state conditions, it achieves significantly shorter circuit depths compared to circuits based on quantum phase estimation (QPE). Numerical results suggest that compared to QPE, the circuit depth can be reduced by around two orders of magnitude under several settings for estimating ground-state and excited-state energies of certain quantum systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Jipole: A Differentiable ipole-based Code for Radiative Transfer in Curved Spacetimes

Recent imaging of supermassive black holes by the Event Horizon Telescope has relied on exhaustive parameter-space searches, matching observations to large, precomputed libraries of theoretical models. As observational data become increasingly precise, the limitations of this computationally expensive approach grow more acute, creating a pressing need for more efficient methods. In this work, we present Jipole, an automatically differentiable (AD), ipole-based code for radiative transfer in curved spacetimes, designed to compute image gradients with respect to underlying model parameters. These gradients quantify how parameter changes—such as the black hole’s spin or the observer’s inclination—affect the image, enabling more efficient parameter estimation and reducing the number of required images. We validate Jipole against ipole in two analytical tests and then compare pixelwise intensity derivatives from AD with those from finite-difference methods. We then demonstrate the utility of these gradients by performing parameter recovery for an analytical model in three increasingly complex cases for the injected image: ideal, blurred, and blurred with added noise. In most cases, high-accuracy fits are obtained in only a few optimization steps, failing only in cases with extremely low signal-to-noise ratios. These results highlight the potential of AD-based methods to accelerate robust, high-fidelity model-data comparisons in current and future black hole imaging efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Requirements and Conceptual Design of Off-gas Systems for the Reprocessing of Metallic Fuels

An assessment has been conducted to determine how key regulations regarding volatile radionuclide emissions to the atmosphere may apply to the off-gas streams associated with electrochemical reprocessing. The scope of this assessment was based upon a generic electrochemical reprocessing scheme with a throughput rate of 200 MTIHM/y applied to metallic fuel discharged from a sodium fast reactor (SFR), but the findings are able to be translated to other advanced nuclear scenarios as merited. Air dispersion modeling was performed using the EPA CAP-88 model and evaluated the uncontrolled decontamination factors (DFs) that would be required to achieve regulatory compliance with the dose-based limits set forth by EPA regulation 40 CFR 190.10(a). These DFs were compared to those required by fuel cycle–based limits set forth by EPA regulation 40 CFR 190.10(b). Two theoretical sites with disparate climatological conditions were selected for air dispersion modeling (Idaho and Tennessee). The radionuclides modeled included 3 H, 85 Kr, 129 I, and selected alpha-emitting transuranic isotopes (referred to here as 239 Pu-TRU <1y ). It was found that the fuel cycle-based limits in 40 CFR 190.10(b) are most restrictive for 85 Kr and 239 Pu-TRU <1y , with DFs of 3 and 6.1E+09, respectively. The dose-based limit as derived from 40 CFR 190.10(a) could require mitigation of tritium in some scenarios, with an estimated DF of about 3 for the reference scenarios. The fuel cycle-based limit for 129 I resulted in a DF of about 240 for the reference scenario. The need for iodine mitigation based on dose to the public depended upon the physical form of iodine as either particulate or vapor-phase species. Emission of iodine from the facility as a vapor necessitated DFs of about 2 but emission as a particulate would require DFs >6,000 to meet thyroid dose-based limits. Effects of physical form on needed iodine mitigation are significant, but the understanding of speciation of iodine both during electrochemical reprocessing and after release to the atmosphere is limited. The electrochemical processing unit operations were evaluated to identify potential release points for the volatile radionuclides and to assess the potential for retention of the radionuclides within the process (thus decreasing the need for mitigation). Mitigation strategies for 3 H, 85 Kr, 129 I, and 239 Pu-TRU <1y were identified. In all cases, there are reasonably achievable pathways to regulatory compliance, although in some cases additional R&D is merited to verify the chemical speciation of these isotopes and to develop and demonstrate potential treatment technologies for this application. Whether or not additional off-gas controls (beyond common operations such as HEPA filtration and oxygen and moisture control) are needed for any of these regulated or volatile radionuclides depends on the (a) type of facility (NRC-regulated or DOE), (b) used fuel process rate, (c) used fuel burnup and composition, (d) speciation and retention of volatile radionuclides in the process and in the cell gas cleanup system, (e) site-specific parameters such as location, meteorology, stack height, and site boundaries, and (f) levels of conservatism and safety factors used in assessing compliance to air emissions regulations. Performance of this assessment revealed several areas where information is lacking or additional research is required in order to better determine if or what kinds of off-gas control might be needed. First, and most significantly, the understanding of the chemical speciation and physical form and partitioning of iodine during electrochemical processing operations is lacking and prevents the ability to accurately assess the potential iodine mitigation requirements. Future research in this area should be multifaceted and include thermodynamic modeling of iodine speciation in different process steps, experiments to quantify the kinetics of vapor-phase and melt-phase transitions, bench-scale experiments to determine the potential chemical and physical form of iodine emissions from the electrorefining process, and verification of iodine behavior with experiments utilizing operational facilities. Similarly, an improved understanding of iodine behavior in the environment after release from the facility stack will be required to refine dose estimations, as particulate and vapor-phase emissions result in significantly different doses to the MEI.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Tuning Catalytic Sites on Zr 6 O 8 Metal–Organic Framework Nodes via Ligand and Defect Chemistry Probed with tert -Butyl Alcohol Dehydration to Isobutylene

Metal–organic frameworks (MOFs) have drawn wide attention as candidate catalysts, but some essential questions about their nature and performance have barely been addressed. (1) How do OH groups on MOF nodes act as catalytic sites? (2) What are the relationships among these groups, node defects, and MOF stability, and how do reaction conditions influence them? (3) What are the interplays between catalytic properties and transport limitations? To address these questions, we report an experimental and theoretical investigation of the catalytic dehydration of tert -butyl alcohol (TBA) used to probe the activities of OH groups of Zr 6 O 8 nodes in the MOFs UiO-66 and MOF-808, which have different densities of vacancy sites and different pore sizes. The results show that (1) terminal node OH groups are formed as formate and/or acetate ligands present initially on the nodes react with TBA to form esters, (2) these OH groups act as catalytic sites for TBA dehydration to isobutylene, and (3) TBA also reacts to break node–linker bonds to form esters and thereby unzip the MOFs. The small pores of UiO-66 limit the access of TBA and the reaction with the formate/acetate ligands bound within the pores, whereas the larger pores of MOF-808 facilitate transport and favor reaction in the MOF interior. However, after removal of the formate and acetate ligands by reaction with methanol to form esters, interior active sites in UiO-66 become accessible for the reaction of TBA, with the activity depending on the density of defect sites with terminal OH groups. The number of vacancies on the nodes is important in determining a tradeoff between the catalytic activity of a MOF and its resistance to unzipping. Computations at the level of density functional theory show how the terminal OH groups on node vacancies act as Brønsted bases, facilitating TBA dehydration via a carbocation intermediate in an E1 mechanism; the calculations further illuminate the comparable chemistry of the unzipping.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

The transition from field emission to collisional space-charge limited current with nonzero initial velocity

Abstract Multiple electron emission mechanisms often contribute in electron devices, motivating theoretical studies characterizing the transitions between them. Previous studies unified thermionic and field emission, defined by the Richardson-Laue-Dushman (RLD) and Fowler–Nordheim (FN) equations, respectively, with the Child-Langmuir (CL) law for vacuum space-charge limited current (SCLC); another study unified FN and CL with the Mott-Gurney (MG) law for collisional SCLC. However, thermionic emission, which introduces a nonzero injection velocity, may also occur in gas, motivating this analysis to unify RLD, FN, CL, and MG. We exactly calculate the current density as a function of applied voltage over a range of injection velocity (i.e., temperature), mobility, and gap distance. This exact solution approaches RLD, FN, and generalized CL (GCL) and MG (GMG) for nonzero injection velocity under appropriate limits. For nonzero initial velocity, GMG approaches zero for sufficiently small applied voltage and mobility, making these gaps always space-charge limited by either GMG at low voltage or GCL at high voltage. The third-order nexus between FN, GMG, and GCL changes negligibly from the zero initial velocity calculation over ten orders of magnitude of applied voltage. These results provide a closed form solution for GMG and guidance on thermionic emission in a collisional gap.

42 ENGINEERING↗
DOE OSTI2026

Enhanced two-dimensional ferromagnetism in van der Waals β-UTe 3 monolayers

The discovery of local-moment magnetism in van der Waals (vdW) semiconductors down to the single-layer limit has led to a paradigm shift in the understanding of two-dimensional (2D) magnets. The incorporation of strong electronic and magnetic correlations in 2D vdW metals remains a sought-after platform to enable control of emergent quantum phases and to achieve more theoretically tractable microscopic models of complex materials. To date, however, there is limited success in the discovery of such metallic vdW platforms, and f-electron monolayers remain out of reach. Here, we demonstrate that strongly correlated β–uranium tritelluride (β-UTe 3 ) can be exfoliated to the monolayer limit. Unexpectedly, β-UTe 3 remains ferromagnetic in this limit with an enhanced ordering temperature of 35 kelvin, a factor of two larger than its bulk counterpart. Our work establishes β-UTe 3 as a materials platform for investigating and modeling correlated behavior in the monolayer limit and opens numerous avenues for quantum control with, e.g., strain engineering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Stress- and interface-compatible red phosphorus anode for high-energy and durable sodium-ion batteries

Sodium-ion battery is a promising candidate for energy storage application, but the absence of high-capacity and low-cost anode materials significantly limit their practical specific energy and cost. Red phosphorus (RP) possesses a high theoretical specific capacity but suffers from large volume change, low electronic conductivity and unstable solid-electrolyte interphase (SEI). Herein, a hierarchical micro/nanostructured antimony-doped RP/carbon anode was developed, which demonstrate extraordinary electrochemical performance with high initial Coulombic efficiency of ~90%, high areal (~1.7 mAh cm-2) capacity together with good cycle stability and rate capability. Combined experimental and computational studies consistently revealed that such a unique structural design can dramatically accommodate the mechanical stress and moreover effectively restrain the undesired decomposition of electrolyte solvents regardless of electrolytes formulation, resulting in superior structural integrity and thin/robust SEI formation during cycling. The present finding has offered an alternative strategy for stress management and interface engineering on high-capacity alloying-based anode materials.

Liu, Xiang↗
DOE OSTI2023

Estimates of the wavenumber wavelet power spectrum of magnetic fluctuations during magnetic reconnection

Fluctuation analyses of experimental observations generally lack high temporal resolution and are in frequency-space f, contrary to theoretical efforts in wavenumber-space k. This is due to the inherent limits of the Fourier transform, though it is prominent due to the ease of diagnostic implementation. Advances in wavelet-based analysis have provided relief due to its temporal resolution, but in its common use, is still hard to compare to theoretical models. By using the two-point correlation technique in conjunction with large data sets, a wavelet power spectrum in wavenumber-space can be created. Dubbed the wavenumber wavelet power spectrum, this spectrum relates wavenumber to power in time. Further, this analysis technique more closely connects characterizations of experimentally observed fluctuations with other system parameters and theoretical predictions. In this article, we develop the wavenumber wavelet power spectrum using magnetic fluctuations caused by tearing instability driven magnetic reconnection in reproducible, high temperature laboratory plasmas. These dynamic magnetic fluctuations generated in reversed field pinch plasmas are broadband, ranging from the low frequency, 10's of kHz, up to the ion gyroradii frequencies, 100's of kHz. The dominant fluctuations have poloidal and toroidal mode numbers (m,n)=(1,6−10) and can grow to 2%–3% of the mean magnetic field. During these reconnection events, ions, and electrons are energized, magnetic fluctuation amplitudes increase, plasma flow is halted, and the toroidal magnetic flux increases, all on a semi-periodic basis. The newly developed spectrum provides better temporal resolution of spectrum characteristics to correlate with these particle energization phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An Orbital Basis Set for Double Photoionization of Atoms and Molecules

The ab initio theoretical treatment of one-photon double photoionization processes has been limited to atoms and diatomic molecules by the challenges posed by large grid-based representations of the double ionized continuum wave function. To provide a path for extensions to polyatomics, an energy-adapted orbital basis approach is demonstrated that reduces the dimensions of such representations and simultaneously allows larger time steps in time-dependent computational descriptions of double ionization. Additionally, an algorithm that exploits the diagonal nature of the two-electron integrals in the grid basis and dramatically accelerates the transformation between grid and orbital representations is presented. Excellent agreement between the present results and benchmark theoretical calculations is found for H – and Be atoms, as well as the hydrogen molecule, including for the triply differential cross sections that relate the angular distribution and energy sharing of all of the particles in the molecular frame.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Machine learning molecular dynamics simulations toward exploration of high-temperature properties of nuclear fuel materials: case study of thorium dioxide

Predicting materials properties of nuclear fuel compounds is a challenging task in materials science. Their thermodynamical behaviors around and above the operational temperature are essential for the design of nuclear reactors. However, they are not easy to measure, because the target temperature range is too high to perform various standard experiments safely and accurately. Moreover, theoretical methods such as first-principles calculations also suffer from the computational limitations in calculating thermodynamical properties due to their high calculation-costs and complicated electronic structures stemming from f-orbital occupations of valence electrons in actinide elements. Here, we demonstrate, for the first time, machine-learning molecular-dynamics to theoretically explore high-temperature thermodynamical properties of a nuclear fuel material, thorium dioxide. The target compound satisfies first-principles calculation accuracy because f-electron occupation coincidentally diminishes and the scheme meets sampling sufficiency because it works at the computational cost of classical molecular-dynamics levels. We prepare a set of training data using first-principles molecular dynamics with small number of atoms, which cannot directly evaluate thermodynamical properties but captures essential atomistic dynamics at the high temperature range. Then, we construct a machine-learning molecular-dynamics potential and carry out large-scale molecular-dynamics calculations. Consequently, we successfully access two kinds of thermodynamic phase transitions, namely the melting and the anomalous λ transition induced by large diffusions of oxygen atoms. Furthermore, we quantitatively reproduce various experimental data in the best agreement manner by selecting a density functional scheme known as SCAN. Our results suggest that the present scale-up simulation-scheme using machine-learning techniques opens up a new pathway on theoretical studies of not only nuclear fuel compounds, but also a variety of similar materials that contain both heavy and light elements, like thorium dioxide.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Synthesis of Honeycomb-Structured Beryllium Oxide via Graphene Liquid Cells

Using high-resolution transmission electron microscopy and electron energy-loss spectroscopy, we demonstrate that beryllium oxide crystallizes in the planar hexagonal structure in a graphene liquid cell by a wet-chemistry approach. These liquid cells can feature van-der-Waals pressures up to 1 GPa, producing a miniaturized high-pressure container for the crystallization in solution. The thickness of as-received crystals is beyond the thermodynamic ultra-thin limit above which the wurtzite phase is energetically more favorable according to the theoretical prediction. The crystallization of the planar phase is ascribed to the near-free-standing condition afforded by the graphene surface. Our calculations show that the energy barrier of the phase transition is responsible for the observed thickness beyond the previously predicted limit. These findings open a new door for exploring aqueous-solution approaches of more metal-oxide semiconductors with exotic phase structures and properties in graphene-encapsulated confined cells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗